matlab声音信号的采集与滤波处理.pdfVIP

. 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 威尼斯的泪 .wav); %读取语音信号的数据，赋给变量 x1 sound(x1,22050); %播放语音信号 y1=fft(x1,1024); %对信号做 1024 点 FFT变换 word 资料 . 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); 实验效果 word 资料 . word 资料 . （2）实验程序 clear; fs=22050; 威尼斯的泪 .wav); f=fs*(0:511)/1024; t=0:1/22050:(length(x1)-1)/22050; % 将所加噪声信号的点数调整到与原始信号相同 %Au=1 d=[0.5*cos(2*pi*1000*t)]; % 噪声为 1kHz 的余弦信号 x2=x1+d; %sound(x1,8000); %pause(50); sound(x2,22050); % 播放加噪声后的语音信号 y2=fft(x2,1024); f