论文中仿真程序，我怎么编都不一样。求各位大神指点一下。

woshiqiao

图片中是论文中的仿真实例。。
下面是我的程序：
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);

我找不到哪里错了。可是就是和文章中的不一样（相差很大，我做出的是一条折线）。还有这个是一个周期的程序。怎么弄成文章中的那样，显示多个周期的？
本人新手。很多不了解。诚信求教。望各位大神指点。

yugang2010

你直接将数列加起来就可以了。[s(t) s(t) s(t).....]，matlab里[]可以将数列加起来。

yummyking

学习了。我想问这个

1325729927

猫头鹰先生

公式编辑不对，可分开S1+S2，时间轴为0-0.16，你是0.008！
祝好！

woshiqiao

猫头鹰先生 发表于 2015-6-3 13:15
公式编辑不对，可分开S1+S2，时间轴为0-0.16，你是0.008！
祝好！

能给我具体代码吗？公式编辑不对在哪？我是按照论文里面的一个周期编的啊。。只是不知怎么扩展到N个周期里面去。

yummyking

这种样式的信号。怎么求故障频率，我已经有了编程思路，楼主加我球球号462427140。你给我说怎么求故障频率，以及其他各种频率，我给你画图。

woshiqiao

其实，类似的仿真程序，我有。可是不是自己编的。诊断方法还在研究中。我只想知道这道题我的程序出错在哪。还有顺便练练手，学习一下下次碰到这样的该怎么弄。就是应为你上面那个仿真例子。我怀疑我找到论文中的那个仿真例子出错了。exp（）这个函数里面是不是也应该是负数？

yummyking

上午画的，可以这个怎么求故障频率啊

ChaChing

不清楚所附的式子含意, 但直觉怪怪的
凭猜测给参考

1. T=1/25000; f1=3000; f2=8000; t=0:T:0.008; t0=mod(t*T,0.01);
   
2. s=exp(-800*(t-0.002)).*(sin(2*pi*f1*t)+1.2*sin(2*pi*f2*t)).*(t>=0.002&t<=0.008);
   
3. plot(t,s);

复制代码

woshiqiao

你这个有噪声没？要是没噪声直接用FFT就可以把故障频率测出来吧。。有噪声就先消噪再用其他方法求就行。你随便搜一篇论文应该就可以找到测故障频率的方法。毕竟这是仿真信号。比较有规律。比较容易测出。

woshiqiao

ChaChing 发表于 2015-6-4 15:33
不清楚所附的式子含意, 但直觉怪怪的
凭猜测给参考

谢谢你，感觉是这样的。那个怎么把这一个周期的弄成N个周期的？

woshiqiao

woshiqiao 发表于 2015-6-4 16:18
谢谢你，感觉是这样的。那个怎么把这一个周期的弄成N个周期的？

大神，二楼是一个周期的（题目中给的条件是0-0.008s的）。怎么弄成多个周期的？比如我想看0-0.16s的图怎么编程（就像图中显示的那样）？？

唯芥末

starryman

无法下载或查看附件