关于一个振动实验的matlab仿真问题

screw_520

最近,编写了一个随机共振的简单程序,参数都是选得别人做过的实验参数,但是仿真只成功了一次....让人很费解,不知道有没有哪位高手可以帮小弟看下这个问题!~

u=1;
A=0.3;
f=0.01;
D=0.5;
fs=5;
x1=[1:0.2:1000];
y1=A*cos(2*pi*f*x1);
y2=awgn(y1,D);
x=[1:1000];
y0=0;
[x,y]=ode45(@lzw_func,[1,1000],y0,0.2);
subplot(3,1,1),plot(x1,y1);
subplot(3,1,2),plot(x1,y2);
subplot(3,1,3),plot(x,y);

其中lzw_func函数的定义如下:
function dy=lzw_func(x,y);
u=1;
A=0.3;
f=0.01;
D=0.5;
fs=5;
x=[1:0.2:1000];
y1=A*cos(2*pi*f*x);
y2=awgn(y1,D);
dy=u*y-y^3+y2;

实验内容就是在一个余弦信号中加入强度为D的高斯白噪声,然后使用龙阶-库塔算法来进行计算.....但是试着修改程序了N次,在成功了一次后再没成功了....

请高手帮忙看下,如果有疑问请发下帖,本人QQ:178855161,经常在线,有兴趣爱好方向相同的朋友可以加好友讨论下!~谢谢各位!

octopussheng

看看这个帖子，希望有帮助

http://forum.vibunion.com/forum/ ... 3076&highlight=