我晕头了，还望请教

harejavahill

经常我们提到信噪比，都知道是越大越好，可是有谁能告诉我，污染语音去噪后的信噪比怎么计算啊？是用去噪后的语音功率/噪音功率吗？但去噪前的信噪比是去噪前的语音功率/噪音功率，这两个信噪比比较起来，发现是后者大，这岂不是说明，去噪后的语音还没有去噪后好啊？？？？
我有点晕头了，还望请教！！！

eight

这句话本身是正确的，可能你没有理解好。信噪比的定义在不同领域可能不一样，但都是指信号的有用成份与噪声的强弱对比。另外，去噪后，计算其信噪比时公式中分母对应的噪声决不是去噪前的噪声，不要搞错了。

[ 本帖最后由 eight 于 2007-3-3 02:58 编辑 ]

harejavahill

谢谢，其实你说的这个我也知道，但就是不知道怎么求，概念有些模糊，还望指教，例如：
原始语音信号：x1
加入白噪音后的语音信号s=x1+0.04*randn(size(x1));
去噪前的信噪比是snr=10log10(power1/power2)
其中power1表示s的功率，power2表示噪音0.04*randn(size(x1))的功率，去噪前的信噪比这样求是正确的吗？？？？？

然后用阀值法对s进行去噪，得到了去噪后的语音信号c
那么这个时候的信噪比怎么求呢？？？
snr=snr=10log10(power1/power2)
其中power1表示c的功率，那么 power2表示谁的功率呢，怎么求啊，难道是用c的功率减去x1的功率吗？？？
还请高手点拨！！！！！！！谢谢

eight

原帖由 harejavahill 于 2007-3-3 10:15 发表
谢谢，其实你说的这个我也知道，但就是不知道怎么求，概念有些模糊，还望指教，例如：
原始语音信号：x1
加入白噪音后的语音信号s=x1+0.04*randn(size(x1));
去噪前的信噪比是snr=10log10(power1/power2)
其 ...

你还没有看懂我写的“信号的有效成分”的意思，去噪前，power1是x1而不是s。至于去噪后，我认为应该是：power1依然是x1，power2是(c-x1)的功率，建议看看关于信噪比的精华贴

harejavahill

十分感谢eight，但是我按照你所说的那样，在matlab中编程，比较去噪前与去噪后的信噪比，结果令人十分失望。程序如下：
去噪前：
fs=22050;                  %语音信号采样频率为22050
x1=wavread('d:\huzhoushifanxueyuan.wav'); %读取语音信号的数据，赋给变量x1   
noise=0.04*randn(size(x1));
s=x1+noise;
sigpower=sum(abs(x1(:)).^2)/length(x1(:));
nospower=sum(abs(noise(:)).^2)/length(noise(:));
snr=10*log10(sigpower/nospower);
disp(snr)
得到的结果为10.9
去噪后：
fs=22050;                  %语音信号采样频率为22050
x1=wavread('d:\huzhoushifanxueyuan.wav'); %读取语音信号的数据，赋给变量x1   
noise=0.04*randn(size(x1));
s=x1+noise;
xd=wden(s,'sqtwolog','s','one',3,'sym8');
%对上述信号进行图示
sigpower=sum(abs(x1(:)).^2)/length(x1(:));
noise=xd-x1;
nospower=sum(abs(noise(:)).^2)/length(noise(:));
snr=10*log10(sigpower/nospower);
disp(snr)
得到的结果为：5.6
这是在matlab中运行的结果，发现去噪后的信噪比反而变小了，很奇怪，希望指教，谢谢！！！！！！！

eight

原帖由 harejavahill 于 2007-3-3 11:25 发表
十分感谢eight，但是我按照你所说的那样，在matlab中编程，比较去噪前与去噪后的信噪比，结果令人十分失望。程序如下：
去噪前：
fs=22050;                  %语音信号采样频率为22050
x1=wavread('d:\huzh ...

1. 你前后两次读取文件得到x1，两者是否一样的？
2. 你前后两次用randn生成噪声，显然不对，随机的结果令两次噪声不一样，snr结果不可信

eight

1. 
   
2. x1=sin(1:0.05:10);
   
3. noise=0.04*randn(size(x1));
   
4. s=x1+noise;
   
5. sigpower=sum(abs(x1(:)).^2)/length(x1(:));
   
6. nospower=sum(abs(noise(:)).^2)/length(noise(:));
   
7. snr=10*log10(sigpower/nospower);
   
8. disp(snr)
   
9. xd=wden(s,'sqtwolog','s','one',3,'sym8');
   
10. %对上述信号进行图示
    
11. noise=xd-x1;
    
12. nospower=sum(abs(noise(:)).^2)/length(noise(:));
    
13. snr=10*log10(sigpower/nospower);
    
14. disp(snr)

复制代码

以上代码我运行了6次，每次显示的信噪比都是后者比前者大，没有问题

[ 本帖最后由 eight 于 2007-3-3 11:51 编辑 ]

harejavahill

这个倒是奇怪了，两次读取的x1是一样的，其中x1是我自己用话筒在pc上录制的一段语音。我也按照你的程序运行了几次，的确是对的，可运行我的程序就不一样了呢。
我就把你的x1改成了我自己录制的语音信号，可结果还是去噪后的信噪比比去噪前的小，如下程序：
fs=22050;                  %语音信号采样频率为22050
x1=wavread('d:\huzhoushifanxueyuan.wav'); %读取语音信号的数据，赋给变量x1   
noise=0.04*randn(size(x1));
s=x1+noise;
sigpower=sum(abs(x1(:)).^2)/length(x1(:));
nospower=sum(abs(noise(:)).^2)/length(noise(:));
snr=10*log10(sigpower/nospower);
disp(snr)
xd=wden(s,'sqtwolog','s','one',3,'sym8');
%对上述信号进行图示
noise=xd-x1;
nospower=sum(abs(noise(:)).^2)/length(noise(:));
snr=10*log10(sigpower/nospower);
disp(snr)
还望多多指教，谢谢！~！！
因为我不是信号专业的，所以有很多不是很懂，还望eight耐心指教！！

eight

原帖由 harejavahill 于 2007-3-3 11:52 发表
这个倒是奇怪了，两次读取的x1是一样的，其中x1是我自己用话筒在pc上录制的一段语音。我也按照你的程序运行了几次，的确是对的，可运行我的程序就不一样了呢。
我就把你的x1改成了我自己录制的语音信号，可结果 ...

1. 试试修改噪声的强度，即调整0.04这个值
2. 用其他去噪方法试试，或者调整wden函数的参数。有可能这个去噪方法对你的信号不适用
3. plot x1 和 xd，一般去噪效果有否改善可以通过肉眼分辨出来

[ 本帖最后由 eight 于 2007-3-3 11:58 编辑 ]

harejavahill

谢谢eight,当我调整了噪音的强度后，发现只有在一定的区域里，去噪后的信噪比比去噪前的信噪比要大，这是为什么啊？？希望指教。
而且信噪比都比较小，一般在2.0左右，这是正确的吗？？