数字信号处理实验报告(00001).docVIP

数字信号处理实验报告 郑州航空工业管理学院 《数字信号处理》实验报告 专 业 电子信息工程 学 号 姓 名 将其从横截型转换为级联型结构，并画出各种结构的流程图。 实验程序 num=[0,10,0]; den=[1,-3,2]; [z,p,k]=tf2zp(num,den) [b,a]=zp2tf(z,p,k) num=[1.9,2.5,2.5,1.9]; den=[1,-6,5,-0.4]; [sos,g]=tf2sos(num,den) [b,a]=sos2tf(sos,g) z=[1,-3]; p=[2,-4]; k=5; [sos,g]=zp2sos(z,p,k) [z1,p1,k1]= sos2zp(sos,g) %直接型到并联型的转换 function [C,B,A]=dir2par(num,den) M=length(num);N=length(den); [r1,p1,C]=residuez(num,den);%先求系统的单根p1，对应的留数r1及直接项C p=cplxpair(p1eps);%用配对函数cplxpair由p1找共轭复根p I=cplxcomp(p1,p); %找p1变为p时的排序变化 r=r1(I); %让r1的排序变化为r，保持与极点对应 %变换为二阶子系统 K=floor(N/2);B=zeros(K,2);A=zeros(K,3);%二阶子系统变量的初始化 if K*2==N; %N为偶数，A(z)的次数为奇,有一个因式是一阶的 for i=1:2:N-2 Brow=r(i:1:i+1,:); %取出一对留数 Arow=p(i:1:i+1,:); %取出一对对应的极点 %二个留数极点转为二阶子系统分子分母系数 [Brow,Arow]=residuez(Brow,Arow,[]); B(fix((i+1)/2),:)=real(Brow);%取Brow的实部,放入系数矩阵B的相应行 A(fix((i+1)/2),:)=real(Arow);%取Arow的实部,放入系数矩阵A的相应行 end; [Brow,Arow]=residuez(r(N-1),p(N-1),[]);%处理实单根 B(K,:)=[real(Brow),0];A(K,:)=[real(Arow),0]; else %N为奇数，A(z)的次数为偶,所有因式都是二阶的 for i=1:2:N-1 Brow=r(i:1:i+1,:); %取出一对留数 Arow=p(i:1:i+1,:); %取出一对对应的极点 %二个留数极点转为二阶子系统分子分母系数 [Brow,Arow]=residuez(Brow,Arow,[]); B(fix((i+1)/2),:)=real(Brow);%取Brow的实部,放入系数矩阵B的相应行 A(fix((i+1)/2),:)=real(Arow);%取Arow的实部,放入系数矩阵A的相应行 end end num =[8 -4 11 -2]; den =[1 -1.25 0.75 -0.125]; [C,B,A]=dir2par(num,den) C = 16 B = -16.0000 20.0000 8.0000 0 A = 1.0000 -1.0000 0.5000 1.0000 -0.2500 0 试验结果分析 实验二 用冲激响应不变法设计IIR数字滤波器 实验目的 加深对冲激响应不变法设计IIR数字滤波器的基本原理的理解； 掌握用冲激响应不变法设计数字低通、带通滤波器的设计； 了解MATLAB有关冲激响应不变法的常用子函数。 实验原理 1. 设计思想 是将系统从s平面到z平面的一种映射方法，使数字滤波器的单位冲激响应序列模仿模拟滤波器的单位冲激响应，其变换关系式为。 2. 设计步骤 （(1)确定数字滤波器的性能指标：通带临界频率、阻带临界频率；通带内的最大衰减；阻带内的最小衰减； ?(2)确定相应的数字角频率； ? (3)将数字滤波器的数字指标转换为模拟滤波器的指标； ? (4)根据计算模拟低通原型滤波器的阶数N，并求得低通原型的传递函数； ? (5