数字信号处理实验指导书(三).docVIP

“数字信号处理”实验指导书（三） 一、实验课程编码：105003 二、实验课程名称：数字信号处理 三、实验项目名称：数字滤波器结构的MATLAB实现 四、实验目的 掌握应用MATLAB实现IIR、FIR数字滤波器结构的方法，即熟悉应用MATLAB实现IIR、FIR数字滤波器结构的相关函数。 五、主要设备 安装有MATLAB软件的电脑 六、实验内容 编写MATLAB 程序，实现以下题目： 1．已知某四阶IIR滤波器的系统函数如下 分别实现系统的级联型结构和并联型结构，并求系统的单位脉冲响应和单位阶跃响应。 2. 用直接型结构、级联型结构实现系统函数 并求系统的单位脉冲响应和单位阶跃响应。 3. 用直接型结构、并联型结构实现系统函数 并求系统的单位脉冲响应和单位阶跃响应。 4. 实现下列线性相位FIR滤波器的级联型结构，求系统的单位脉冲响应和单位阶跃响应。 （1） （2） 七、实验步骤 1、熟悉与离散信号频谱分析相关的MATLAB函数（参考附录1）； 2、通过运行附录2中提供的例题，熟悉用MATLAB实现数字滤波器结构的基本方法； 3、根据“六、实验内容”中各个题目的要求，编写MATLAB程序代码，调试程序，分析并保存结果。 八、实验结果 对实验练习题编写MATLAB程序并运行，在计算机上输出仿真结果。 附录1 主要的相关MATLAB函数 1. y＝filter（b,a,x） 2. 〔z,p,k〕＝tf2zp（num，den） 3. 〔r,p,k〕＝residuez（num，den） 附录2 例题 例1、用直接型实现系统函数为 的IIR数字滤波器，求单位脉冲响应和单位阶跃响应。 [MATLAB程序]： %用直接型结构实现IIR数字滤波器 b=[1,-3,11,27,18]; a=[16,12,2,-4,-1]; N=25; h=impz(b,a,N); %直接型单位脉冲响应 x=[ones(1,5),zeros(1,N-5)]; %单位阶跃信号 y=filter(b,a,x); % 直接型输出信号 subplot(1,2,1); stem(h,.); line([0,30],[0,0]); title(直接型单位脉冲响应); subplot(1,2,2); stem(y,.); line([0,30],[0,0]); title(直接型单位阶跃响应); [运行结果]： 例2、用级联型结构实现系统函数为 的IIR数字滤波器，求单位脉冲响应和单位阶跃响应。 [MATLAB程序]： %用级联型结构实现直接形式的系统函数 N=25; %delta=impseq(0,0,N); delta=[1,zeros(1,N-1)]; b=[1,-3,11,27,18]; a=[16,12,2,-4,-1]; [b0,B,A]=dir2cas(b,a); h=casfilter(b0,B,A,delta); x=[ones(1,5),zeros(1,N-5)]; y=casfilter(b0,B,A,x); subplot(1,2,1); stem(h,.); line([0,25],[0,0]); title(级联型单位脉冲响应); subplot(1,2,2); stem(y,.); line([0,25],[0,0]); title(级联型单位阶跃响应); [运行结果]： 例3、用并联型结构实现系统函数 [MATLAB程序]： %用并联型结构实现直接形式的系统函数 b=[1,-3,11,27,18]; a=[16,12,2,-4,-1]; N=25; %delta=impseq(0,0,N); delta=[1,zeros(1,N-1)]; [C,B,A]=dir2par(b,a); x=[ones(1,5),zeros(1,N-5)]; h=parfilter(C,B,A,delta); %并联型单位脉冲响应 y=parfilter(C,B,A,x); %并联型单位阶跃响应 subplot(1,2,1); stem(h,.); line([0,25],[0,0]); title(并联型单位脉冲响应); subplot(1,2,2); stem(y,.); line([0,25],[0,0]); title(并联型单位阶跃响应); [运行结果]： 例题中的函数: function [b0,B,A]=dir2cas(b,a) %直接型结构到级联型结构的转换子程序 %Direct-form to Cascade-form conversion in z-plane %(complex pair version) % [b0,B,A]=dir2