matlab声音的滤波处理,MATLAB声音信号的采集与滤波处理

Matlab 实验报告(题目二)

(题目二)声音信号的采集与滤波处理(采用IIR滤波器或FIR滤波器)

参考资料：信号的采集、数字信号处理及滤波实例

要求：(1)采集声音信号或打开已录好的声音文件，并显示其信号图与频域图。

(2)根据信号的特点，选用合适的滤波器，给定滤波器的规一化性能指标(参考指标,实际中依据每个同学所叠加噪声情况而定)例如：通带截止频率wp=0.25*pi, 阻通带截止频率ws=0.3*pi; 通带最大衰减Rp=1 dB; 阻带最小衰减Rs=15 dB，对信号进行滤波。

在Matlab中，可以利用函数fir1设计FIR滤波器，可以利用函数butte，cheby1和ellip设计IIR滤波器；利用Matlab中的函数freqz画出各滤波器的频率响应，滤波器设计完后，用filter函数用这些数字滤波器对含噪语音信号分别进行滤波处理。

(3)还原音乐信号，并画出其时域图与频域图，并与原始信号比较，且回放音乐信号。

(1)打开一个自己录制的音乐文件进行实验，这是实验程序：

fs=22050; %语音信号采样频率为22050

x1=wavread('e:\威尼斯的泪.wav'); %读取语音信号的数据，赋给变量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');

实验效果