小波分析.m

来自「音频处理Matlab程序」· M 代码 · 共 30 行

clear all;%在噪声环境下语音信号的增强
sound=wavread('Blip.wav');%语音信号的读入
cound=length(sound);
noise=0.05*randn(1,cound);
y=sound'+noise;%随机噪声的加载
[thr,sorh,keepapp]=ddencmp('den','wv',y);%获取噪声的阈值
yd=wdencmp('gbl',y,'db4',2,thr,sorh,keepapp);%对信号进行消噪
subplot(2,2,1);plot(sound);
xlabel('采样点');
ylabel('模拟输入 (V)');
title('原始语音信号');
subplot(2,2,3);plot(yd);
xlabel('采样点');
ylabel('模拟输入 (V)');
title('去噪后的语音信号');
[C,L]=wavedec(sound,3,'haar');%用小波函数haar对信号进行3层分解
alpha=1.5;
[thr,nkeep]=wdcbm(C,L,alpha);%获取信号压缩的阈值
[cp,cxd,lxd,per1,per2]=wdencmp('lvd',C,L,'haar',3,thr,'s');%对信号进行压缩
subplot(2,2,2);plot(sound);
xlabel('采样点');
ylabel('模拟输入 (V)');
title('原始语音信号');
subplot(2,2,4);plot(cp);
xlabel('采样点');
ylabel('模拟输入 (V)');
title('压缩后的语音信号');