数字信号处理_上机作业3.doc

目 录 数字信号处理第三次上机实验报告 1 实验目的 1 实验一 1 一、实验题目 1 二、 程序 1 1、脉冲响应不变法 1 2、双线性变换法 2 3、 波形验证程序 3 三、 运行结果 3 1、 脉冲响应不变法 3 2、 双线性变换法 5 四、 结论 6 实验二 6 一、实验题目 6 二、程序 6 三、运行结果 7 1、滤波器的单位脉冲响应 7 2、 频率响应特性曲线 8 四、 结论 8 数字信号处理第三次上机实验报告 实验目的 掌握基于MATLAB通过脉冲响应不变法、双线性变换法设计数字滤波器（IIR）的方法及步骤。 用含有低频与高频成分的混合信号验证所设计的滤波器。 掌握频率采样法设计严格线性相位数字滤波器（FIR）的方法及步骤。 实验一 一、实验题目 用脉冲响应不变法或双线性变换法设计一个数字低通巴特沃斯滤波器。满足一下参数：通带截止频率：，通带最大衰减：，阻带截止频率：，阻带最大衰减：。 绘制该滤波器的幅频响应与相频响应曲线图； 通过一个含有低频与高频成份的混合信号，进行滤波验证。 程序 1、脉冲响应不变法 wp=0.2*pi;%通带截止频率 ws=0.3*pi;%阻带截止频率 rp=1;%通带最大衰减 rs=15;%阻带最小衰减 [N,wc]=buttord(wp,ws,rp,rs,s);%计算对应模拟滤波器的阶数N及3dB衰减频率 [Bs,As]=butter(N,wc,s);%模拟滤波器传递函数分子分母多项式的系数 [Bz,Az]=impinvar(Bs,As);%用脉冲响应不变法将模拟滤波器转换为数字滤波器 w=linspace(0,pi,500); [H,w1]=freqz(Bz,Az,w); P=20*log10(angle(H)); H=20*log10(abs(H)); subplot(211) plot(w/pi,H),grid on title(幅频响应曲线) xlabel(\omega/\pi) ylabel(|H(e^j^\omega)|/dB) subplot(212) plot(w/pi,P),grid on title(相频响应曲线) xlabel(\omega/\pi) ylabel(相位\Phi/dB) 2、双线性变换法 wp=0.2;%通带截止频率 ws=0.3;%阻带截止频率 rp=1;%通带最大衰减 rs=15;%阻带最小衰减 [N,wc]=buttord(wp,ws,rp,rs); [Bz,Az]=butter(N,wc); w=linspace(0,pi,500); [H,w1]=freqz(Bz,Az,w); P=20*log10(angle(H)); H=20*log10(abs(H)); figure(1) subplot(211) plot(w/pi,H),grid on title(幅频响应曲线) xlabel(\omega/\pi) ylabel(|H(e^j^\omega)|/dB) subplot(212) plot(w/pi,P),grid on title(相频响应曲线) xlabel(\omega/\pi) ylabel(相位\Phi/dB) 波形验证程序 t=0:0.001:2-0.001;%连续自变量 f1=25;f2=250; x=sin(2*pi*f1*t)+sin(2*pi*f2*t); y=filter(Bz,Az,x); figure(2) subplot(211)%原信号 plot(t,x) xlabel(t/s),ylabel(幅度) title(原信号波形) axis([0 100/f2 min(x)-0.2 max(x)+0.2]) subplot(212)%滤波后信号波形 plot(t,y) xlabel(t/s),ylabel(幅度) title(滤波后信号波形) axis([0 100/f2 min(y)-0.2 max(y)+0.2]) 运行结果 脉冲响应不变法 ①滤波器特性曲线 ②滤波验证 设,，让该信号通过低通滤波器，则高频成份被滤除。波形如下图所示： 双线性变换法 ①滤波器特性曲线 ②滤波验证 结论 根据脉冲响应不变法设计的数字低通巴特沃斯滤波器，在通带截止频率处的衰减为0.92dB，在阻带截止频率处的衰减为15dB，均满足设计指标。 根据双线性变换法设计的数字低通巴特沃斯滤波器，在通带截止频率处的衰减为0.564dB，在阻带截止频率处的衰减为15dB，均满足设计指标。 从混合信号的验证中也可以看出，所设计的滤波器能够完成低频滤波的功能。 实验二 一、实验题目 采用频率采样法设计一个严格线性相位FIR数字低通滤波器，要求通带截止频率，过渡带宽度，阻带最小衰减，允许在过渡带上设置采样点。画