实验报告模板-实验四.doc

实验四 数字滤波器的结构与MATLAB实现 课程名称：数字信号处理 实验日期： 2017.6.2 班级： 数媒1403 姓名： 李欣 学号： 1030514327 一、实验目的掌握：(1) IIR滤波器结构及其MATLAB实现; (2) FIR滤波器结构及其MATLAB实现二、实验要求按照实验材料内容，利用MATLAB编程实现相关任务。三、实验步骤 （以实验材料第二部分“FIR滤波器结构及其MATLAB实现”为例，给出实验步骤描述、相关代码等。）解：FIR直接型系统函数表示FIR结构的直接形式由包含系数 n b 的行向量b 描述。在MATLAB中，把分母矢量a置为1，这种结构用filter函数实现。直接型的代码实现：impseq.mfunction n,x = impseq(np,ns,nf)%数字信号处理教程MATLAB释义与实现% 单个脉冲序列生成子程序% 电子工业出版社出版 陈怀琛编著 2004年9月% 产生 x(n) = delta(n-np); ns = n,np Mb=b zeros(1,N-M);elseif M Na= a zeros(1,M-N);N=M;elseNM=0end%K=floor(N/2);B=zeros(K,3);A=zeros(K,3);if K*2=N;b=b 0;a=a 0;end%broots=cplxpair(roots(b);aroots=cplxpair(roots(a);for i=1:2:2 *KBrow=broots(i:1:i+1,:);Brow=real(poly(Brow);B(fix(i+1)/2),:)=Brow;Arow=aroots(i:1:i+1,:);Arow=real(poly(Arow);A(fix(i+1)/2),:)=Arow;endcastfiltr.mfunction y=casfiltr(b0,B,A,x);% IIR各FIR滤波器的级联型的实现% % y=casfiltr(b0,B,A,x);% y=输出序列% b0=级联型的增益系数% B=包含各bk和K乘3维实系数矩阵%A=包含各ak的K乘3维实系数矩阵% x=输入序列%K,L=size(B);N=length(x);w=zeros(K+1,N);w(1,:)=x;for i=1:1:Kw(i+1,:)=filter(B(i,:),A(i,:),w(i,:);endy=b0 * w(K+1,:);b0,B,A=dir2cas(b,1);h=casfiltr(b0,B,A,delta);y=casfiltr(b0,B,A,x);subplot(2,2,3);stem(h);title(级联型h(n);subplot(2,2,4);stem(y);title(级联型y(n);解：若给定直接型的h(n)或H(k )，可以由dir2fs函数转换为频率采样结构。dir2fs.mfunctionC,B,A=dir2fs(h)%C=并联部分增益的行向量%B=按行排列的分子系数矩阵%A=按行排列的分母系数矩阵%h=FIR滤波器的单位脉冲响应向量M=length(h);H=fft(h,M);magH=abs(H);phaH=angle(H);if(M=2*floor(M/2)L=M/2-1;A1=1,-1,0;1,1,0;C1=real(H(1),real(H(L+2);elseL=(M-1)/2;A1=1,-1,0;C1=real(H(1);endk=1:L;B=zeros(L,2);A=ones(L,3);A(1:L,2)=-2*cos(2*pi*k/M);A=A;A1;B(1:L,1)=cos(phaH(2:L+1);B(1:L,2)=cos(phaH(2:L+1)-(2*pi*k/M);C=2*magH(2:L+1),C1;h=1,2,3,2,1/9;