数字信号处理DSP 课程设计报告

1,4);stem(x3);title('fft相角')Р所得语音频谱分析图见图3.3.1Р图3.3.1 声音信号的频谱分析Р(第一幅是时域波形;采样1是FFT后波形;采样2是频域幅度谱;采样3是频域相位谱)Р3.4滤波器设计Р设计指标:通带截止频率为1100Hz,阻带截止频率为1200Hz,通带波纹为1dB,阻带波纹为20dB,用脉冲(冲激)响应不变法设计的一个满足上述指标的巴特沃思IIR滤波器Р%巴特沃思低通滤波器РFs=8000; %采样频率Рwp=1100*2/Fs;Рws=1200*2/Fs;%根据采样频率将滤波器边界进行转化РRp=1;РRs=20;%通带波纹和阻带波纹РNn=128;Р[N,Wn]=buttord(wp,ws,Rp,Rs);%求滤波器的最小阶数和3dB固有频率Р[b,a]=butter(N,Wn)%设计BUTTERWORTH低通滤波器,返回滤波器的系数矩阵Рfigure(2);Р[H,f]=freqz(b,a,Nn,Fs);%用Nn点求得频率特性,H为N个点处的频率响应复值输出向量;F为与第N点处对应的频率值f(Hz)Рsubplot(2,1,1);plot(f,20*log10(abs(H)));title('滤波器特性')%20log10(abs(H))DB为频响幅值曲线幅值Рxlabel('频率/Hz');ylabel('振幅/dB');grid on;%滤波器幅频特性Рsubplot(2,1,2);plot(f,180/pi*unwrap(angle(H)));%angle(H)为频响相位曲线相位值Рxlabel('频率/Hz');ylabel('相位/^o');grid on;%滤波器相频特性Р巴特沃思低通滤波器的相位谱和振幅谱见图3.4.1Р图3.4.1 利用脉冲响应不变法设计的数字巴特沃斯滤波器(w(单位:)Р3.5 信号滤波信号