| 查看: 1465 | 回复: 13 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
trend17新虫 (小有名气)
|
[求助]
考察算法抗噪性时在matlab里用awgn函数加入高斯白噪声的一些问题？？？？？ 已有2人参与
|
||
|
最近在写一篇算法的文章，在抗噪性分析这块我用的是awgn函数，Y=AWGN(X,10)，加入了信噪比是10的高斯白噪声，然后没有设定sigpower，所以应该默认就是0dBw了，然后我发现信噪比10db的噪声加入后序列幅值变化也不是特别大，而且我把所以幅值相加，发现加了噪声后总幅值反而变小了... 而且不知道是不是没有设定sigpower，所以就算是10db的噪声，抗噪性还是很好，不知道这边有没有大神也搞过这方面的，能给点加噪技巧的指点，多谢多谢！  jiazao.jpg |
回复此楼» 猜你喜欢
086502化学工程342求调剂
已经有4人回复
-
材料求调剂一志愿哈工大324
已经有7人回复
-
328求调剂
已经有7人回复
-
340求调剂
已经有5人回复
-
一志愿南昌大学324求调剂
已经有4人回复
-
求调剂推荐 材料 304
已经有15人回复
-
291求调剂
已经有19人回复
-
331环境科学与工程求调剂
已经有3人回复
-
289求调剂
已经有8人回复
-
315分求调剂
已经有7人回复
trend17
新虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 450.4
- 散金: 61
- 帖子: 71
- 在线: 17.1小时
- 虫号: 3816510
- 注册: 2015-04-17
- 性别: GG
- 专业: 多相流热物理学
|
恩恩，多谢，刚回去用了awgn(x,10,'measured')后发现10db的噪声确实是蛮强的... 话说本身awgn(x,10)里的x并没有normalized，用awgn(x,10,'measured')后matlab就帮忙估计了signal power，那么我不用'measured'，直接用awgn(x,10)后就没有估计signal power值了，然后x也没有normalized，这样子的加噪应该是不规范的吧？换句话说，您能不能从专业角度解释下用awgn(x,10,'measured')比awgn(x,10)的优点，我担心自己描述的不够好，只是因为它在加入噪声之前测定信号强度吗？ 真的真的麻烦了，万分感谢！ |
11楼2015-07-29 12:03:34
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
- 专业: 马、恩、列、斯思想研究
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
