用MATLAB编程分析噪声对混沌振子的影响的时候怎么加入噪声?

lehehea

用MATLAB编程分析噪声对混沌振子的影响的时候怎么加入噪声?
研究duffing振子检测微弱信号，但是不知道噪声怎么加入，求高人指点啊

friendchj

help rand,randn,imnoise

lehehea

谢谢你，加噪声的这些函数，但是不知道怎么用在分析噪声对混沌振子的影响上。我把我的程序列出来：
主程序为：

1. global w;
   
2. w=1；
   
3. tspan=[0:0.01:1000]；
   
4. x0=[1,1；
   
5. [t,z]=ode45('f',tspan,x0);
   
6. plot(z(:,1),z(:,2))%相图
   
7. 子程序f为：
   
8. function dx=f(t,x)
   
9. global w;
   
10. k=0.5;
    
11. f=0.717282；
    
12. y=randn(1,length(t));
    
13. y=y/std(y);
    
14. y=0.00005*y;%噪声均值为0，标准差为0.000005
    
15. dx(1)=w*x(2);
    
16. dx(2)=w*(-k*x(2)+x(1)^3-x(1)^5+f*cos(w*t)+y);%
    
17. dx=[dx(1);dx(2)];

复制代码

这样加噪声对吗？为什么我的程序运行时总是出现：y=y/std(y)错误：diveded by zero.
还有，李月、杨宝俊著《混沌振子检测引论》里第55页说到可检测的信噪比为-111.46dB，这时的周期信号功率和噪声方差分别是什么？周期信号功率是取的0.5*10^(-16)吗？

ChaChing

不懂LZ的专业东东!
但y=y/std(y)错误：diveded by zero是一定的, 因y=randn(1,length(t))仅是一个值, std(y)就一定为零!

lehehea

y=randn(1,length(t))不是产生一行，length(t)列的高斯随机数吗？它的标准差怎么是0呢？

ChaChing

ode45呼叫函数f(t,x)是一次一个t和x!
length(t)=1

lehehea

哦，谢谢啊，这个我给弄混了
那如果我想加入0均值的高斯噪声，应该怎么办呢？

ChaChing

真的不太懂LZ的专业东东! 但改成这样是有输出的
global w;
w=1; tspan=[0:0.01:1000]; x0=[1,1];
[t,z]=ode45('f',tspan,x0); plot(z(:,1),z(:,2))%相图

function dx=f(t,x)
global w;
k=0.5; f=0.717282; y=0.00005*randn; dx = zeros(2,1);
dx(1)=w*x(2);
dx(2)=w*(-k*x(2)+x(1)^3-x(1)^5+f*cos(w*t)+y);

lehehea

谢谢啊，我先看看

64471059

你好,我也是在做这方面，能不能给我这个加噪的程序做个参考，我想加非高斯的，我的邮箱64471059@qq.com，谢谢

cumtchenyy

x0错了吧

yzsldj

ChaChing 发表于 2009-9-23 10:58
真的不太懂LZ的专业东东! 但改成这样是有输出的
global w;
w=1; tspan=[0:0.01:1000]; x0=[1,1];

噪声好像没有加上去啊。