音乐信号的滤波处理MATLAB课程设计.docVIP

河北北方学院信息科学与工程学院 《MATLAB与仿真》 期末作业 课程名称 音乐信号的滤波处理 实验学期 2013 至 2014 学年 第 2 学期 学生所在系部 信息科学与工程学院 年级 专业班级 学生姓名 学号 任课教师 成绩 年级与专业 11级电子信息工程 姓名 日期 2014.3.26 实验室 B601 课程题目 音乐信号的滤波处理 一、课题 1、选取一段MP3文件，格式为.WAV格式； 2、做波形显示以及fft变换 3、在语音信号上分别叠加均匀白噪声和高斯白噪声,使信噪比为（学号）dB； 4、对叠加噪声前后的信号进行频谱分析，确定降噪的滤波器指标； 5、根据滤波器指标利用窗函数法设计滤波器，在Matlab平台下编写程序，观察分析滤波器的幅频特性、相频特性和群延时，以及滤波前后信号时域特性和频域特性。 二、设备与环境 软件环境： MATLAB7.0 三、课题内容 原语音信号分析： 1.1实验程序 [x1,fs,bits]=wavread('蔡琴 - 渡口 - ape 无损音乐 高质量.wav'); %读取语音信号的数据，赋给变量x1 sound(x1,fs); %播放语音信号 y1=fft(x1,4096); %对信号做4096点FFT变换 f=fs*(0:511)/4096; 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、叠加高斯白噪声,使信噪比为75dB 2.1实验程序 [x1,fs,nbits]=wavread('蔡琴 - 渡口 - ape 无损音乐 高质量.wav'); y=awgn(x1,75); sound(y); figure; subplot(2,1,1); plot(y); title('加噪语音信号的时域波形','fontweight','bold'); S=fft(y,4096); %傅里叶变换 subplot(2,1,2); plot(abs(y)); title('加噪语音信号的频域波形','fontweight','bold'); 2.2仿真结果 3、窗函数法设计滤波器 实验程序 [y,fs,bits]=wavread('蔡琴 - 渡口 - ape 无损音乐 高质量.wav'); n = length (y) ; %求出语音信号的长度 s=awgn(y,75); S=fft(s); %傅里叶变换 Ft=4000; Fp=600; Fs=1800; wp=2*Fp/Ft; ws=2*Fs/Ft; rp=1; rs=50; p=1-10.^(-rp/20); %通带阻带波纹 q=10.^(-rs/20); fpts=[wp ws]; mag=[1 0]; dev=[p q]; [n21,wn21,beta,ftype]=kaiserord(fpts,mag,dev);%由kaiserord求滤波器的阶数和截止频率 b21=fir1(n21,wn21,kaiser(n21+1,beta)); %由fir1设计滤波器 z21=fftfilt(b21,s); sound(z21); m21=fft(z21); %求滤波后的信号 figure(4); subplot(2,2,1); plot(abs(S),'g'); title('滤波前信号的频谱','fontweight','bold'); grid; subplot(2,2,2); plot(abs(m21),'r'); title('滤波后信号的频谱','fontweight'