论文中仿真程序,我怎么编都不一样。求各位大神指点一下。
图片中是论文中的仿真实例。。下面是我的程序:
T=1/25000;
f1=3000;f2=8000;
t=0:T:0.008;
t0=mod(t*T,0.01);
s=exp(800*t0).*(sin(2*pi*f1*T*t.*(t>=0.002&t<=0.008))+1.2*sin(2*pi*f2*T*t.*(t>=0.002&t<=0.008)));
plot(t,s);
我找不到哪里错了。可是就是和文章中的不一样(相差很大,我做出的是一条折线)。还有这个是一个周期的程序。怎么弄成文章中的那样,显示多个周期的?
本人新手。很多不了解。诚信求教。望各位大神指点。
你直接将数列加起来就可以了。,matlab里[]可以将数列加起来。 学习了。我想问这个 {:{03}:} 公式编辑不对,可分开S1+S2,时间轴为0-0.16,你是0.008!
祝好! 猫头鹰先生 发表于 2015-6-3 13:15
公式编辑不对,可分开S1+S2,时间轴为0-0.16,你是0.008!
祝好!
能给我具体代码吗?公式编辑不对在哪?我是按照论文里面的一个周期编的啊。。只是不知怎么扩展到N个周期里面去。 这种样式的信号。怎么求故障频率,我已经有了编程思路,楼主加我球球号462427140。你给我说怎么求故障频率,以及其他各种频率,我给你画图。 其实,类似的仿真程序,我有。可是不是自己编的。诊断方法还在研究中。我只想知道这道题我的程序出错在哪。还有顺便练练手,学习一下下次碰到这样的该怎么弄。就是应为你上面那个仿真例子。我怀疑我找到论文中的那个仿真例子出错了。exp()这个函数里面是不是也应该是负数? 上午画的,可以这个怎么求故障频率啊 本帖最后由 ChaChing 于 2015-6-4 15:34 编辑
不清楚所附的式子含意, 但直觉怪怪的
凭猜测给参考
T=1/25000; f1=3000; f2=8000; t=0:T:0.008; t0=mod(t*T,0.01);
s=exp(-800*(t-0.002)).*(sin(2*pi*f1*t)+1.2*sin(2*pi*f2*t)).*(t>=0.002&t<=0.008);
plot(t,s); 你这个有噪声没?要是没噪声直接用FFT就可以把故障频率测出来吧。。有噪声就先消噪再用其他方法求就行。你随便搜一篇论文应该就可以找到测故障频率的方法。毕竟这是仿真信号。比较有规律。比较容易测出。 ChaChing 发表于 2015-6-4 15:33
不清楚所附的式子含意, 但直觉怪怪的
凭猜测给参考
谢谢你,感觉是这样的。那个怎么把这一个周期的弄成N个周期的? woshiqiao 发表于 2015-6-4 16:18
谢谢你,感觉是这样的。那个怎么把这一个周期的弄成N个周期的?
大神,二楼是一个周期的(题目中给的条件是0-0.008s的)。怎么弄成多个周期的?比如我想看0-0.16s的图怎么编程(就像图中显示的那样)?? {:{10}:} 无法下载或查看附件
页:
[1]
2