实验七 数字滤波器结构.docVIP

实验七 滤波器设计 实验目的： 掌握滤波器结构类型转换的Matlab实现方法；掌握用滤波器直接形式、级联形式、并联形式、格型形式以及频率采样形式的Matlab实现。了解和学习组合滤波的结构类型特性；进一步理解和巩固理论知识，提高分析和解决实际问题的能力。 实验原理： 1．滤波器数学模型之间的转换 2．滤波器结构类型之间的转换 三．实验内容： 程序用到的自定义子函数程序： 1.IIR系统级联转直接型子函数： function [b,a]=cas2dir(b0,B,A) [K,L]=size(B); b=[1]; a=[1]; for i=1:1:K b=conv(b,B(i,:)); a=conv(a,A(i,:)); end b=b*b0; 2.IIR系统直接转级联型子函数： function [b0,B,A]=dir2cas(b,a); b0=b(1);b=b/b0; a0=a(1);a=a/a0; b0=b0/a0; M=length(b);N=length(a); if NM b=[b zeros(1,N-M)]; elseif MN a=[a zeros(1,M-N)];N=M; else NM=0; end K=floor(N/2);B=zeros(K,3);A=zeros(K,3); if K*2==N; b=[b 0]