发博文
个人资料
梦龙
梦龙
  • 博客等级:
  • 博客积分:1260
  • 博客访问:21,956
  • 关注人气:14
公告
谨以此怀念我的本科——
穷乡僻壤偏要惊天动地
偏安一隅也要心怀天下
 
不写稿很多年
再提笔已无力
图片播放器
留言
加载中…
友情链接

孙志岗的博客

计算机学院的老师好像都很喜欢写blog?哈哈~

刘挺

突然有个小学妹说刘老师是她偶像~

郧叶

很特殊的好友,至于有多特殊嘛,我也说不清道不明

评论
加载中…
访客
加载中…
好友
加载中…
博文
标签:

matlab

comsol

校园

分类: 学术交流
    如果说matlab在解偏微分方程时,性能不佳,那么comsol则很好地互补上了。当然,更好的消息就是这两个软件的连接比较简单,互相调用方便。
http://baike.baidu.com/view/656888.htm
    COMSOL公司是全球多物理场建模与仿真解决方案的提倡者和领导者,其旗舰产品COMSOL Multiphysics,使工程师和科学家们可以通过模拟,赋予设计理念以生命。它有无与伦比的能力,使所有的物理现象可以在计算机上完美重现。COMSOL的用户利用它提高了手机的接收性能,利用它改进医疗设备的性能并提供更准确的诊断,利用它使汽车和飞机变得更加安全和节能,利用它寻找新能源,利用它探索宇宙,甚至利用它去培养下一代的科学家。
    COMSOL Multiphysics起源于MATLAB的Toolbox,最初命名为Toolbox 1.0。后来改名为Femlab 1.0(FEM为有限元,LAB是取自于Matlab),这个名字也一直沿用到Femlab3.1。从2003年3.2a版本开始,正式命名为COMSOL Multiphysics。
    一看这两软件这么有渊源,就知道联合仿真,有戏。具体实现步骤如下
1.系统配置
阅读  ┆ 评论  ┆ 转载 ┆ 收藏 
    matlab能知道现在是几点钟吗?能知道现在的日期吗?
    可以通过调用几个特殊字符来实现哦~还能因此生成真随机数
1.基于系统时间的真随机数
1.1伪随机数
    先给个说明,matlab的随机数是伪随机数,比如你输入rand,输出一个数
    你关掉matlab,再开,再输入rand,出来的还是这个数。。。晕了吧。。。
    伪随机数,实际上是通过一个固定的、可以重复的计算方法产生的。计算机或计算器产生的随机数有很长的周期性。它们不真正地随机,因为它们实际上是可以计算出来的,但是它们具有类似于随机数的统计特征。这样的发生器叫做伪随机数发生器。
    其实matlab的随机数原理是有一组已经生成的伪随机数,可以用seed来说明:
    输入rand('seed',0),就是调用第零个seed的伪随机数!
1.2真随机数
    如何变成真随机数呢?只要把seed的序号,变成和系统时间有关就行了嘛!比如现在是凌晨一点,就嗲用第一个seed,凌晨两点就调用第二个seed
1.3clock
阅读  ┆ 评论  ┆ 转载 ┆ 收藏 
标签:

matlab

fft

傅里叶

校园

分类: 学术交流
注:数学本非我专业,以下说法欢迎纠错
1.傅里叶变换的定义
傅里叶变换从数学上的定义,为
F(w)=int(x*exp(-jwt),t=-inf...inf)
其中,int表示积分,t是时间,x是时域信号,inf表示无穷,exp表示指数运算。其含义说的是给一个无限长的时域信号和一个频点w,可以唯一确定一个复数F。于是,F和w就有了这种对应关系,考虑到F是个复数。F的绝对值和w关系叫幅频,F的幅角和w关系叫相频。
2.matlab的fft命令
help fft可以知道这个和数学上的傅里叶不一样,因为计算机是离散的!因为计算机的时域信号存储量是有限的!比如等时采样得到的信号,高频分量是不可能获得的,对于比较大的w将无法计算。于是,fft这样计算傅里叶变换:把时域信号进行周期延拓,取一组w,就是时域信号的周期及该周期的二分之一,三分之一,直到n分之一,其中n是一个周期内的采样点。这样做的结果,就是对一段有限长的时域信号,将其长度作为基频率,分析基频和高频含量。当然,能分析到的最高频为n次谐波,再高次谐波由于香农定理而无法体现。
3.写一个数学定义傅里叶变换的程序
将有限长时域信号不延拓,时域信号
阅读  ┆ 评论  ┆ 转载 ┆ 收藏 
标签:

simulink

pid控制

校园

分类: 学术交流
    学过控制理论的大概都有体会,数学推导太复杂,都怨恨拉普拉斯吧~其实有了matlab后,很多过程得到了简化,比如即将要介绍的simulink自动设计反馈控制器。先普及一下基本知识——
    任何一个系统,不管是线性还是非线性,都可以写成如下形式:
dx=f(x,u)
y=g(x,u)
其中x是列向量,dx是列向量对时间t的导数。u是系统的输入,y是系统的输出,都是列向量,在此假设系统比较简单,是个单输入单输出系统(SISO),也就是说,u和y都是一维的,也就是标量。
    所谓设计控制器,就是设计一个函数h
u=h(y,yref,x)
其中yref是y的目标值。将这个方程带入原来的系统,形成闭环系统
dx=f(x,h(y,yref,x))
y=g(x,h(y,yref,x))
考虑到x一般无法直接测量,实际系统中只能得到y,因此最常用的一种h是
e=yref-y
h=kp*e+ki*int(e)+kd*de
其中kp、ki、kd是固定不变的常数,int表示积分,de表示微分。这样的h也被称为(串联负反馈)pid控制。而其设计难点,则是这三个
阅读  ┆ 评论  ┆ 转载 ┆ 收藏 
    今天上微博,看见有个数学题,1、5、8、11找规律,填出后三个数。然后还说这是小学生一两年级的寒假作业,其本意是说明现在的作业偏离教育本质,已经妖魔化。但撇开题目的出处,这其实是数学上的重要问题:拟合、预测。
    拟合和预测的方法很多,matlab的拟合工具箱cftool就非常强大,用cftool求解这个数学题,大致就能体会着工具箱的用法了。步骤如下:
1.打开matlab,输入x=[1 2 3 4];y=[1 5 8 11];这其实就是输入题目中给出的样本数据。此时如果输入plot(x,y),则可以看见x和y的关系,直观上感觉是个近似线性,所以可以给出答案为1、5、8、11、14、17、20。当然,这是直观感受,答案是猜的,猜得对不对不清楚。
2.输入cftool,弹出curve fitting tool的gui界面。
3.如图所示做些设置。这些设置归结为两条:其一,告诉matlab,自变量和应变量分布是多少,也就是此问题中,x和y的样本是多少。由于左边看见了xyz三个data框,因此推测cftool可以支持一维和二维的拟合。其二是告诉matlab需要使用的拟合方法,比如我选择的三阶多项式拟合,打开框框,可以有其他选择。
阅读  ┆ 评论  ┆ 转载 ┆ 收藏 
(2012-01-21 19:52)
标签:

新年

许愿

校园

分类: 心声
元旦新年错过了许愿,可不能再错过春节许愿啦!写下龙年的愿望(排名不分先后),等明年看能实现多少吧~
1.学术交流几次。作为科研男,最重要的就是做点成绩,然后同行认可,顺便看看别人在做什么,掺和一下就变成学科交叉了。读研快两年了,做了点东西,是时候整理一下,登陆学术圈吧
2.不要恋爱。感谢我的女朋友,快三十年了还坚持没出现,让我专心科研~原来只是一句玩笑话,现在突然变得现实了~
3.在水源、论坛等等社交net上好好干吧~现在是水源mathtools和hit的板主,希望能为板面做点什么
4.老爸老妈亲朋好友们,身体健康!顺顺利利!
5.如果有闲暇的话,一个人
阅读  ┆ 评论  ┆ 转载 ┆ 收藏 
标签:

matlab

示波器

tmtool

校园

分类: 学术交流
    最近新到一个泰克的示波器,用matlab读取之,在对面实验室大牛帮助之下,一天就把程序写得差不多了,然后花了三天时间找出其中一个非常低级的错误,总计四天搞定。如果你也有类似的任务,可以参考如下步骤:
    1.启动matlab,输入tmtool,弹出新窗口。
    2.左边,test&measurement窗口下,打开hardware树,选择示波器与电脑的连接方式。比如我选择的是visa下的usb连接。之前可能需要安装usb驱动程序,买示波器的时候会有盘。据说有些情况下也可以不装驱动,因为matlab有自带的!
    3.中间主窗口,scan,进行搜索,搜索到示波器后,就可以按下“connect”进行连接了,如图所示


 &nbs
阅读  ┆ 评论  ┆ 转载 ┆ 收藏 
标签:

哈工大

校园

分类: 日记
    离开母校一年半了,总嚷嚷着要回去看看,哪怕面对读研生活的身不由己,我也从未停止过想象母校重游。借着项目间隙,一周时间,我请假的极限,我回来了。2011年的最后几天,我愿陪在母校身旁。
    24日,凌晨多,上海,躺着,兴奋。几个小时后,飞机,哈尔滨,降落,寒风显得意外熟悉,反而少了几分刺骨。对着网络大吼几声:“我在哈哈哈哈哈哈哈哈哈哈尔滨啦啦啦啦啦啦啦”,感觉就像大一离家入学时一样兴奋。好友们一阵羡慕嫉妒恨。趁着圣诞节平安夜,我履行自己的承诺,给记者团一位女生送苹果,于是,我的圣诞节苹果记录改写为收到两个,送出一个,目标收支平衡。晚上和寝室兄弟们小聚,大米说签了工作,在上海,天南地北而来,天南地北而散,却总有相遇的时候。
    25日,和春哥去的食堂,母校食堂依然给力,物价仿佛停在两年前。和记者团同学去了极乐寺,烧香拜佛,心诚则灵。这样,我也算去过哈尔滨名胜古迹了。晚上去了医大,见一位记者团女生,居然去看了场电影,我可是很少看电影的工科男呀,忆往昔,青涩记忆。
    26日,回校报编辑部。下飞机的时候,
阅读  ┆ 评论  ┆ 转载 ┆ 收藏 
标签:

2011总结

校园

分类: 日记
元旦的时候,许下好多心愿,实现了多少呢?先回顾下愿望吧~
http://blog.sina.com.cn/s/blog_4eff3a0e0100nblk.html
    1.结婚
没实现。而且好像还越来越遥远了,心中的她原来只是个错觉。。。
    2.新电脑
部分实现。实验室的新电脑的确是到位了,不过是大家一起用
    3.英语
实现啦。哈哈,终于克服英语障碍了,居然会用英语写长篇大论,还在交流会上和老外秀英语~
    4.减肥
没实现。这个愿望。。。我已经放弃了。。。
    5.自行车
实现啦。哈哈,学会自行车了!不过天太冷,把车给别人,自己坐校车吧
    6.回哈尔滨
实现啦。这个愿望实现啦!虽然是年底实现的~
    7.天下太平
部分实现。据说好多人都上船了
    8.和谐小家
实现啦。家里好和谐的
    9.中大奖
没实现。因为没去买彩票。。
阅读  ┆ 评论  ┆ 转载 ┆ 收藏 
标签:

2011总结

校园

分类: 日记
    今天在群里突然有人提到说,大家2010年度总结都很勤快,2011都懒了。。。转念一想,是啊,2011就这么浑浑噩噩地结束了——照例应该写个年度十大新闻,排名不分先后:
1.因为打篮球过于拼命,双脚骨折,床上躺了3个月,幸好研究生没啥课,可以请长假
2.失去了挂科的权利,因为研究生课程都修完了。。。至于以后的读博之类,都没有考试
3.步莎姐后尘,相亲几次。。。是不是我表现得太饥渴了?结果班里同学帮我各种介绍,不过都不太靠谱哎。。。其实是我不靠谱吧,习惯了孤单,就继续孤单吧~曾一度让我错认为是幸福的人啊~
4.投了篇论文,算是正式进入学术圈吧。。。虽然我一直对这种未成熟的论文持保留意见,不过导师说练练手,咱就练练吧
阅读  ┆ 评论  ┆ 转载 ┆ 收藏 
  

新浪BLOG意见反馈留言板 不良信息反馈 电话:4006900000 提示音后按1键(按当地市话标准计费) 欢迎批评指正

新浪简介 | About Sina | 广告服务 | 联系我们 | 招聘信息 | 网站律师 | SINA English | 会员注册 | 产品答疑

新浪公司 版权所有