数字信号处理综合设计.doc

数字信号处理综合设计 1.语音信号的频谱分析 要求首先画出语音信号的时域波形;然后对语音信号进行频谱分析,在MA TLAB中,可以利用函数fft对信号进行快速付立叶变换,得到信号的频谱特性;从而加深对频谱特性的理解。 fs=20000; %语音信号采样频率为20000 x1=wavread(d:\lianxi.wav,20000; %读取语音信号的数据,赋给变量x1 sound(x1,20000; %播放语音信号 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; 设计数字滤波器和画出频率响应 根据语音信号的特点给出有关滤波器的性能指标:1低通滤波器性能指标,fp=1000Hz, fc=1200 Hz,As=100dB,Ap=1dB;2高通滤波器性能指标,fc=2800 Hz,fp=3000 Hz As=100dB,Ap=1dB;3带通滤波器性能指标,fp1=1200 Hz,fp2=3000 Hz,fc1=1000 Hz, fc2=3200 Hz,As=100dB,Ap=1dB。要求学生首先用窗函数法设计上面要求的三种滤波器,在MA TLAB中,可以利用函数fir1设计FIR滤波器,然后在用双线性变换法设计上面要求的三种滤波器;之后再利用函数butter和cheby1设计上面要求的三种IIR滤波器。最后,利用MATLAB中的函数freqz画出各滤波器的频率响应。 巴特沃思低通1 wp=0.1*pi; ws=0.5673*pi; Rp=1; Rs=100; Fs=20000; Ts=1/Fs; wp1=2/Ts*tan(wp/2; %将模拟指标转换成数字指标 ws1=2/Ts*tan(ws/2; [N,Wn]=buttord(wp1,ws1,Rp,Rs,s; %选择滤波器的最小阶数 [Z,P,K]=buttap(N; %创建butterworth模拟滤波器 [Bap,Aap]=zp2tf(Z,P,K; [b,a]=lp2lp(Bap,Aap,Wn; [bz,az]=bilinear(b,a,Fs; %用双线性变换法实现模拟滤波器到数字滤波器的转换[H,W]=freqz(bz,az; %绘制频率响应曲线 figure(1 plot(W*Fs/(2*pi,abs(H grid xlabel(频率/Hz ylabel(频率响应幅度 title(Butterworth f1=filter(bz,az,x1; 巴特沃思高通2 x1=wavread(d:\lianxi.wav,20000; t=0:1/20000:(size(x1-1/20000; wp=2*pi*3000; ws=2*pi*1555; Rp=1; Rs=100; Fs=20000; Ts=1/Fs; [N,Wn]=buttord(wp,ws,Rp,Rs,s; %选择滤波器的最小阶数 [Z,P,K]=buttap(N; %创建butterworth模拟滤波器 [Bap,Aap]=zp2tf(Z,P,K; [b,a]=lp2hp(Bap,Aap,Wn; [bz,az]=bilinear(b,a,Fs; %用双线性变换法实现模拟滤波器到数字滤波器的转换[H,W]=freqz(bz,az; %绘制频率响应曲线 figure(1 plot(W*Fs/(2*pi,abs(H grid xlabel(频率/Hz ylabel(频率响应幅度 title(Butterworth f1=filter(bz,az,x1; 巴特沃思带通3 wp=2*pi*[1200,3000]; ws=2*pi*[212,5546]; Rp=1; Rs=100; Fs=20000; Ts=1/Fs; [N,Wn]=buttord(wp,ws,Rp,Rs,s; %选择滤波器的最小阶数 [B,A]=butter(N,Wn,s %创建buttord模拟滤波器 [Bz,Az]=bilinear(B,A,Fs; %用双线性变换法实现模拟滤波器到数字滤波器的转换 [H,W]=freqz(Bz,Az; %绘制频率响应曲线