疑问中，不知道怎么加噪声

lehehea

用duffing振子检测微弱正弦信号的程序
主程序为：
global w;
w=1；tspan=[0:0.01:1000]；x0=[1,1]；[t,z]=ode45('f',tspan,x0);plot(z(:,1),z(:,2))%相图
子程序f为：
function dx=f(t,x)
global w;
k=0.5;f=0.717282；
y=randn(1,length(t)); y=0.00005*y;%加入的白噪声？？？是这样加入吗？
dx(1)=w*x(2);
dx(2)=w*(-k*x(2)+x(1)^3-x(1)^5+f*cos(w*t)+y);%
dx=[dx(1);dx(2)];
但是这样并不能得到很低的信噪比，文献中提到的-100多的信噪比更不同提了，这是怎么回事？
到底该怎么加噪声啊？在线等牛人给解答，多谢

octopussheng

用wgn试试

lehehea

我用了，但是不能得到那么低的噪声，谢谢