湖南大学数字信号处理实验报告四 - 图文(2)

[Bi,Ai]=residuez(r(N-1),p(N-1),[]);%处理单根 B(K,:)=[real(Bi) 0];A(K,:)=[real(Ai) 0]; else

for i=1:2:N-2

pi=p(i:i+1,:);%取出一对极点 ri=p(i:i+1,:);%取出一对对应留数 [Bi,Ai]=residuez(ri,pi,[]); B(fix((i+1)/2),:)=real(Bi); A(fix((i+1)/2),:)=real(Ai); end end

cplxcomp函数：

function I=cplxcomp(p1,p2) I=[];

for j=1:length(p2) for i=1:length(p1)

if (abs(p1(i)-p2(j))<0.0001) I=[I,i]; end; end; end; I=I';

实验结果：

c. FIR 滤波器的方程和系统函数分别为

y(n)?b0x(n)?b1x(n?1)?b2x(n?2)?b3x(n?3)?b4x(n?4)

画出该滤波器的直接II型结构图。

(1) 当b0?1,b1??2,b2?6,b3??1,b4?3时，求系统的单位冲激响应，并绘制响应的波形。

(2) 当b0?2,b1??1,b2?4,b3??3,b4??2时，求系统的单位冲激响应，并绘制响应的波形。

（1） 系统1：

系统二