语音信号去噪.doc

文档介绍：语音信号去噪 matlab实验问题
 
题目要求:
   ,并对录制的信号进行采样;
   ;
   ,采用窗函数法和双线性变换法设计滤波器,
   并划出滤波器的频域响应;
   ,画出滤波后信号的时域波形和频谱,
   并对滤波前后的信号进行对比,分析信号的变化;
   ;
   。
   报告要求:
  根据自己录制的语音信号选择合适的滤波器,并按照题目要求完成报告;
参考程序:
程序1:用MATLAB对原始语音信号进行分析,画出它的时域波形和频谱
fs=22050;                  %语音信号采样频率为22050
x1=wavread('h:\课程设计2\'); %读取语音信号的数据,赋给变量x1
sound(x1,22050);           %播放语音信号
y1=fft(x1,1024);           %对信号做1024点FFT变换
f=fs*(0:511)/1024;
figure(1)
plot(x1)                   %做原始语音信号的时域图形
title('原始语音信号');
xlabel('time n');
ylabel('fuzhi n');
figure(2)
freqz(x1)                  %绘制原始语音信号的频率响应图
title('频率响应图')
figure(3)
subplot(2,1,1);
plot(abs(y1(1:512)))       %做原始语音信号的FFT频谱图
title('原始语音信号FFT频谱')
subplot(2,1,2);
plot(f,abs(y1(1:512)));
title('原始语音信号频谱')
xlabel('Hz');
ylabel('fuzhi');
 
程序2:给原始的语音信号加上一个高频余弦噪声,频率为5kHz。画出加噪后的语音信号时域和频谱图,与原始信号对比,可以很明显的看出区别。
fs=22050;
x1=wavread('h:\课程设计2\');
f=fs*(0:511)/1024;
t=0:1/22050:(size(x1)-1)/22050; %将所加噪声信号的点数调整到与原始信号相同
Au=;
d=[Au*cos(2*pi*5000*t)]'; %噪声为5kHz的余弦信号
x2=x1+d;
sound(x2,22050); %播放加噪声后的语音信号
y2=fft(x2,1024);
figure(1)
plot(t,x2)
title('加噪后的信号');
xlabel('time n');
ylabel('fuzhi n');
figure(2)
subplot(2,1,1);
plot(f,abs(y1(1:512)));
title('原始语音信号频谱');
xlabel('Hz');
ylabel('fuzhi');
subplot(2,1,2