DSP实验报告6.docVIP

DSP第六、七次实验报告 1．实验目的： （1）进一步熟悉Matlab实验环境和语言。 （2）熟悉各种滤波器的结构及Matlab实现语言。 （3）掌握用冲击响应不变法和双线性变换法设计IIR滤波器的方法。 （4）掌握用窗函数法和频率抽样法设计FIR滤波器的方法。 2．实验内容及总结： 1.滤波器结构： IIR滤波器各种结构 直接型结构 例如直接型滤波器系统函数，则有系数向量a=[1,a1,a2,a3],b=[b0,b1,b2]，利用： Y=filter[b,a,x] 求信号x(n)通过此滤波器的输出。 由系统函数或差分方程求系统的二阶分式（含一阶分式）的级联结构 将例如的系统函数重写为二阶分式节的级联型，利用： [sos,G]=tf2sos(b,a) 由二阶分式的级联结构转换成系统函数的直接结构 是第二步的逆运算，调用函数： [b,a] = sos2tf(sos) 可以求得系数向量a,b，从而得到H(z) 由系统函数求部分分式展开（留数及其极点计算） 即求z反变换的部分分式展开法，利用： [r,p,c]=residuez(b,a) 其中极点为p，留数为r，直接项系数为c。 由r,p,c求系统函数 即第4步的逆运算，利用： [b,a]=residuez(r,p,c) 由直接型结构转换为并联型结构 需开发函数： [C,B,A]=tf2par(b,a) 其中，b,a为直接型的系数向量，C,B,A为并联型实系数向量，基本思想是：1.反复调用[r,p,c]=residuez(b,a)求出极点及留数；2.利用cplxpair函数把极点、留数对按复共轭极点-留数对，实极点-留数对的顺序排列；3.开发cplxcomp函数，保证极点和留数相互对应；4.调用[b,a]=residuez(r,p,c)计算并联二阶节的分子分母。 由并联型结构转换成直接型结构 开发函数： [b,a]=par2tf(C,B,A) 为[C,B,A]=tf2par(b,a)的逆函数。 FIR滤波器各种结构 从FIR滤波器的直接结构转换为级联型结构 与IIR相似，调用函数[sos,G]=tf2sos(b,a)，并令a=1，即： [sos,G]=tf2sos(b,1) 从FIR滤波器的直接结构转换为频率抽样结构 开发函数[C,B,A]=tf2fs(h)： function [C,B,A]=tf2fs(h) N=length(h); H=fft(h,N); R=input(‘r= ’); %输入r,r=1 magH=abs(H);phaH=angle(H)’; %求H(k)的幅度和相位 if (N==2*floor(N/2)) %N为偶数 L=N/2-1; A1=[1,-r,0;1,r,0]; %设置z=r及z=-r两个实极点，A1为其系数 C1=[real(H(1)),real(H(L+2))]; %两个实极点处H(k)的样值 Else %N为奇数 L=(N-1)/2; A1=[1,-r,0]; %设置z=r处的实极点，A1为其系数 C1=[real(H(1))]; %z=1处H(k)的值 end k=[1:L]’; B=zeros(L,2); A=zeros(L,3); A(1:L,1) =1; %计算分母系数 A(1:L,3)=-2*r*cos(2*pi*k/N); A(1:L,3) =r.^2 A=[A;A1]; B(1:L,1)=cos(phaH(2:L+1)); %计算分子系数 B(1:L,2)=cos(phaH(2:L+1)-(2*pi*k/N)); C=[2*magH(2:L+1),C1]’; %计算增益系数 格型滤波器结构 全零点格型滤波器 利用函数tf2latc，将FIR直接型系数向量b转换为FIR格型滤波器的反射系数向量，将b的首项归一化： k=tf2latc(b/b(1)) 逆函数，即由反射系数向量求直接型系数向量： b=latc2tf(k) 全极点格型滤波器 利用函数： k=tf2latc(a) 其中k的次序与全零点的完全颠倒，k1在最后一级，km在第一级。 零极点格型梯形滤波器 利用函数： [k,c]=tf2latc(b,a) 滤波器设计： 冲击响应不变法设计IIR滤波器 [例7.20] 采用冲击响应不变法设计一个数字巴特沃斯低通滤波器，要求：f1kHz，幅度特性下降小于1dB，频率大于f=1.5kHz的阻带内衰减大于15dB，抽样频率为fs=10kHz。 代码： clc;clear all OmegaP=2*pi*1000;OmegaS=2*pi*1500; %巴特沃斯低通滤波器技术指标 Rp=1;As=15;Fs=10*10^3; %抽样频率