s306080135.m

对于小波去噪

%装载语音信号
N=1024;
s=wavread('srcpu1.wav',N);
figure(1);
subplot(311);
plot(1:N,s,'LineWidth',2);
title('原始信号');
xlabel('时间n');
ylabel('幅值A');
s=s+0.001*randn(1,N)';
%用小波db3对s进行5层分解
level=5;
[c,l]=wavedec(s,level,'db3');
%选用全局阈值进行信号去噪处理
thr=5;
[sd,csd,lsd,perf0,perf12]=wdencmp('gbl',c,l,'db3',level,thr,'h',l);
%figure(1);
subplot(312);
plot(s,'LineWidth',2);
title('加噪声后的信号');
%xlabel(‘时间n’)
ylabel('幅值A');
subplot(313);
plot(sd,'LineWidth',2);
title('去噪后的信号');
xlabel('时间n');
ylabel('幅值A');