小波阈值语音消噪clc;clear;fs=wavread('G:\c12345.wav');y=fs(1000:1200

小波阈值语音消噪
clc;
clear;
fs=wavread('G:c12345.wav');
y=fs(1000:12000);
N=length(y);
figure(1);
subplot(111);
plot(y);
ylabel('幅值 A');
title('原始信号');
s=awgn(y,20,'measured');%加入高斯白噪声
figure(2);
plot(s);
ylabel('幅值 A');
title('加噪信号');
wname='db3';%选db3小波基
lev=5;%5层分解
[c,l]=wavedec(s,lev,wname);
a5=appcoef(c,l,wname,lev);
d5=detcoef(c,l,5);
d4=detcoef(c,l,4);
d3=detcoef(c,l,3);
d2=detcoef(c,l,2);
d1=detcoef(c,l,1);
cD=[d1 d2 d3 d4 d5];%运行到这里出错.提示d1,d2,d3,d4,d5维数不同不兼容.