数字信号Mtlab作业.docVIP

生物医学信号处理大作业 题目：布莱克曼窗设计FIR滤波器处理语音信号 学生姓名 学院名称 精密仪器与光电子工程 专 业 生物医学工程?班 学 号 3011202312 题目要求： 录制自己的一段语音：“天津大学精密仪器与光电子工程学院生物医学工程X班XXX, College of precision instrument and opto-electronics engineering, biomedical engineering”，时间控制在15秒到30秒左右；利用wavread函数对自己的语音进行采样，记住采样频率。 （1）求原始语音信号的特征频带（比如谱峰位置）：可以分别对一定时间间隔内，求功率谱（傅里叶变换结果取模的平方）并画出功率谱。 （2）根据语音信号频谱特点，设计FIR或IIR滤波器，分别画出滤波器幅频和相频特性曲线。说明滤波器特性参数。用设计的滤波器对信号滤波，画出滤波前后信号的频谱图。用sound函数回放语音信号,说明利用高通/低通/带通滤波后的效果，不同特征频带被滤除后分别有什么效果。 （3）求出特征频段语音信号随时间变化的曲线（每隔一定时间求一次功率谱，连接成曲线，即短时 FFT）。 （4）选做：语谱图：横轴为时间，纵轴为频率，灰度值大小表示功率谱值的大小。(提示，可以采用spectrogram函数) (5) 选做：分析自己的语音频谱特点，比如中英文发音的区别。 求原始语音信号的特征频带 代码： clc clear all File=(I:\Matlab大作业\work.wav) [Y1,Fs,nbits]=wavread(File); Y1 = Y1(:,1);%取双声道中的一个声道 l=length(Y1); Y=Y1(1:l); Fs=44100; figure(1), plot(Y); title(原声音信号波形) for i = 1:18; x = Y(Fs*(i-1)+1 : Fs*i); z = fft(x , Fs); p = z .*conj(z)/Fs; %功率谱 figure(i); plot(p); xlabel(频率Hz) ylabel(功率) axis([0 2500 0 1]); end 1.1 原始语音信号 1.2 以1秒为时间间隔，求出17秒的语音信号的功率谱如下： 根据语音信号频谱特点，设计FIR或IIR滤波器，分别画出滤波器幅频和相频特性曲线。说明滤波器特性参数。用设计的滤波器对信号滤波，画出滤波前后信号的频谱图。用sound函数回放语音信号,说明利用高通/低通/带通滤波后的效果，不同特征频带被滤除后分别有什么效果。 代码： %==========低通滤波器==========% fp=1300; fs=1500; % 输入给定指标 % 用窗函数法设计滤波器 wc=(fp+fs)/Fs; %理想低通滤波器截止频率(关于pi归一化） B=2*pi*(fs-fp)/Fs; %过渡带宽度指标 N=ceil(12*pi/B); %blackman窗的长度N h1=fir1(N-1,wc,low,blackman(N)); %用窗函数法设计法的绘图部分 [h11,w1]=freqz(h1,1); figure(2) subplot(2,1,1); plot(w1*Fs/pi/2,20*log10(abs(h11))); title(低通滤波器); grid on; axis([0,2500,-150,10]); xlabel(频率/Hz);ylabel(幅度/dB) subplot(2,1,2) plot(w1*Fs/pi/2,angle(h11)); grid on; axis([0,1000,-4,4]); xlabel(频率/Hz); ylabel(角度/rad); %sound(ywt,Fs); %==========高通滤波器==========% fp=300; fs=100; % 输入给定指标 %窗函数法设计滤波器 wc=(fp+fs)/Fs; %理想高通滤波器截止频率(关于pi归一化） B=2*pi*(fp-fs)/Fs; %过渡带宽度指标 N=ceil(12*pi/B); %blackman窗的长度N h2=fir1(N-1,wc,high,blackman(N)); %=====以下为绘图部分===== [h22,w2]=freqz(h2,1); figure(3)