| 查看: 1494 | 回复: 13 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
trend17新虫 (小有名气)
|
[求助]
考察算法抗噪性时在matlab里用awgn函数加入高斯白噪声的一些问题？？？？？ 已有2人参与
|
||
|
最近在写一篇算法的文章，在抗噪性分析这块我用的是awgn函数，Y=AWGN(X,10)，加入了信噪比是10的高斯白噪声，然后没有设定sigpower，所以应该默认就是0dBw了，然后我发现信噪比10db的噪声加入后序列幅值变化也不是特别大，而且我把所以幅值相加，发现加了噪声后总幅值反而变小了... 而且不知道是不是没有设定sigpower，所以就算是10db的噪声，抗噪性还是很好，不知道这边有没有大神也搞过这方面的，能给点加噪技巧的指点，多谢多谢！  jiazao.jpg |
回复此楼» 猜你喜欢
291求调剂
已经有7人回复
-
300求调剂
已经有8人回复
-
求调剂推荐
已经有7人回复
-
289 分105500药学专硕求调剂(找B区学校)
已经有4人回复
-
0854求调剂
已经有13人回复
-
初试324 中药学 一志愿天中医 求调剂
已经有3人回复
-
药学求调剂
已经有14人回复
-
327求调剂
已经有27人回复
-
急需调剂
已经有5人回复
-
271求调剂
已经有33人回复
ll550
木虫 (职业作家)
- 应助: 259 (大学生)
- 金币: 4586.3
- 散金: 230
- 红花: 34
- 帖子: 3331
- 在线: 767小时
- 虫号: 3180771
- 注册: 2014-05-05
- 专业: 马、恩、列、斯思想研究
trend17: 回帖置顶 2015-08-03 21:10:43
|
NOISE=NOISE-mean(NOISE); randn产生的是服从高斯分布[-1,1]的mean =0，var =1的随机数。 mean(NOISE)本身就是0，这句是多余的。 NOISE=sqrt(noise_variance)/std(NOISE)*NOISE; 其次，std(NOISE)本身是1，这句也是多余的。 说道awgn这个函数，本身awgn(x,10)，你的x并没有normalized，所以，你应该用awgn(x,10,'measured')才对。 如果非要用awgn(x,10)，x= x./norm(x)才对。 或者 awgn(x,10) |
5楼2015-07-29 10:41:39
superjintao
木虫 (知名作家)
- 应助: 20 (小学生)
- 金币: 18211.7
- 红花: 12
- 帖子: 5146
- 在线: 233.6小时
- 虫号: 457038
- 注册: 2007-11-12
- 性别: GG
- 专业: 信号理论与信号处理
3楼2015-07-29 07:00:10
trend17
新虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 450.4
- 散金: 61
- 帖子: 71
- 在线: 17.1小时
- 虫号: 3816510
- 注册: 2015-04-17
- 性别: GG
- 专业: 多相流热物理学
|
求幅值和我倒也就是看看的，主要是我觉得考察抗噪性的话加了10db的噪声应该算比较大了，但是序列幅值变化也不是很大，所以担心这样子就算最后检验抗噪性好也没什么说服力... 另外反而用以下程序加入10db噪声后序列幅值就变化的比较大 X=data(1:20000,1); %载入时间序列 SNR=10; %信噪比 NOISE=randn(size(X)); NOISE=NOISE-mean(NOISE); signal_power = 1/length(X)*sum(X.*X); noise_variance = signal_power / ( 10^(SNR/10) ); NOISE=sqrt(noise_variance)/std(NOISE)*NOISE; Y=X+NOISE; 因为实验室这边没人搞这个的，都在自己摸索，还望能帮忙解释下两者区别 |
4楼2015-07-29 09:05:32
ll550
木虫 (职业作家)
- 应助: 259 (大学生)
- 金币: 4586.3
- 散金: 230
- 红花: 34
- 帖子: 3331
- 在线: 767小时
- 虫号: 3180771
- 注册: 2014-05-05
- 专业: 马、恩、列、斯思想研究
6楼2015-07-29 10:42:43
