声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

查看: 2676|回复: 11

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

[复制链接]
发表于 2009-9-14 20:50 | 显示全部楼层 |阅读模式

马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。

您需要 登录 才可以下载或查看,没有账号?我要加入

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

使用道具 举报

发表于 2009-9-14 22:50 | 显示全部楼层
help rand,randn,imnoise

评分

1

查看全部评分

 楼主| 发表于 2009-9-15 21:22 | 显示全部楼层

回复 沙发 friendchj 的帖子

本帖最后由 牛小贱 于 2014-6-14 21:53 编辑

谢谢你,加噪声的这些函数,但是不知道怎么用在分析噪声对混沌振子的影响上。我把我的程序列出来:
主程序为:
  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)吗?
发表于 2009-9-16 00:27 | 显示全部楼层
不懂LZ的专业东东!
但y=y/std(y)错误:diveded by zero是一定的, 因y=randn(1,length(t))仅是一个值, std(y)就一定为零!

评分

1

查看全部评分

 楼主| 发表于 2009-9-20 19:35 | 显示全部楼层

回复 地板 ChaChing 的帖子

y=randn(1,length(t))不是产生一行,length(t)列的高斯随机数吗?它的标准差怎么是0呢?
发表于 2009-9-21 08:45 | 显示全部楼层
ode45呼叫函数f(t,x)是一次一个t和x!
length(t)=1
 楼主| 发表于 2009-9-23 09:35 | 显示全部楼层

回复 6楼 ChaChing 的帖子

哦,谢谢啊,这个我给弄混了
那如果我想加入0均值的高斯噪声,应该怎么办呢?
发表于 2009-9-23 10:58 | 显示全部楼层
真的不太懂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);

评分

1

查看全部评分

 楼主| 发表于 2009-9-24 18:53 | 显示全部楼层

回复 8楼 ChaChing 的帖子

谢谢啊,我先看看
发表于 2012-4-16 11:02 | 显示全部楼层
你好,我也是在做这方面,能不能给我这个加噪的程序做个参考,我想加非高斯的,我的邮箱64471059@qq.com,谢谢
发表于 2012-7-12 19:27 | 显示全部楼层
x0错了吧
发表于 2014-6-14 16:42 | 显示全部楼层
ChaChing 发表于 2009-9-23 10:58
真的不太懂LZ的专业东东! 但改成这样是有输出的
global w;
w=1; tspan=[0:0.01:1000]; x0=[1,1];

噪声好像没有加上去啊。
您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

QQ|小黑屋|Archiver|手机版|联系我们|声振论坛

GMT+8, 2024-11-11 17:02 , Processed in 0.070468 second(s), 24 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

快速回复 返回顶部 返回列表