小波变换语音消噪.docxVIP

改进阈值函数进行语音信号消噪，但是在程序运行过程中频频报错。本人经验不足调试不出，希望求得各位指改进函数表达式附图clearall;clc;closeall;fs=8000;%语音信号采样频率为8000xx=wavread(lw1.wav);x1=xx(:,1);%取单声道t=(0:length(x1)-1)/8000;y1=fft(x1,2048);%对信号做2048点FFT变换f=fs*(0:1023)/2048;figure(1)plot(t,x1)%做原始语音信号的时域图形y=awgn(x1,10,measured);%加10db的高斯白噪声[snr,mse]=snrmse(x1,y)%求得信噪比均方误差figure(2)plot(t,y)%做加噪语音信号的时域图形[c,l]=wavedec(y,3,db1);%多尺度一维分解%用db1小波对信号进行3层分解并提取系数a3=appcoef(c,l,db1,3);%a2=appcoef(c,l,db1,2);%a1=appcoef(c,l,db1,1);d3=detcoef(c,l,3);d2=detcoef(c,l,2);d1=detcoef(c,l,1);thr1=thselect(d1,rigrsure);%阈值获取，使用Stein的无偏风险估计原理thr2=thselect(d2,rigrsure);thr3=thse