数字信号处理作业二（副本）.docVIP

南昌大学实验报告 学生姓名： 王庆 学 号： 55023010031 专业班级： 中兴101班 实验类型：□验证 □ 综合 (设计 □ 创新 实验日期： 实验成绩： 作业二 模拟滤波器转换为数字滤波器 设计目的 掌握数字滤波器的设计方法； 2、熟悉使用MATLAB语言进行编程实现。 二、设计内容 1、试用脉冲响应不变法，将实验一中的H(s)转换成H(z)； 2、画出H(z)的幅频特性和相频特性； 3、验证和的数值关系是否满足变换方法的频率映射关系； 4、用M文件实现上述H(z)：输入信号为9KHz+(+10)KHz,求输出的时域信号y(n)以及画出。 设计条件 计算机、MATLAB软件环境。 设计原理及过程 1、模拟滤波器到数字滤波器的转换可以在时域内实现，也可在频域内实现。时域转换法是使数字滤波器的时域响应与模拟滤波器的时域采样值相等，具体方法有：脉冲响应不变法、阶跃不变法和匹配z变换法。这里我们采用的是脉冲不变响应法。脉冲响应不变法是依据数字滤波器的单位冲激响应h(n)与模拟滤波器的单位冲激响应h(t)在采样点上的值相等，即得到的变换关系。 2、如果给定模拟滤波器的传输函数H (s)，则所求的数字滤波器的系统函数H(z)为： 。 3、暂设模拟滤波器的传输函数只有单阶极点，且分母多项式的阶次高于分子多项式的阶次，将Ha(s)用部分分式表示：（式中为Ha(s)的单阶极点） （2·1） 将Ha(s)进行逆拉氏变换得到ha(t)： ，式中u(t)是单位阶跃函数。 对ha(t)进行等间隔采样，采样间隔为T，得到： 对上式进行z变换，得到数字滤波器的系统函数: （2·2） 对比(2．1)式和(2．2)式，Ha(s)的极点映射到z平面，即极点变成，系数Ai不变化。 4、下面分析从模拟滤波器转换到数字滤波器，s平面和z平面之间的映射关系，从而找到这种转换方法的优点。这里我们以采样信号作为桥梁，推导其映射关系。设的采样信号用表示，即 对进行拉氏变换，得到： 式中ha(nT)是在采样点t=nT式的幅度值，它与序列h(n)幅度值相等，即h(n)=ha(nT)。因此得到： 上式表示采样信号的拉氏变换与相应序列的z变换之间的映射关系用下式表示： （2·3） 我们知道模拟信号ho(t)的傅里叶变换和其采样信号的傅里叶变换之间的关系满足： 将代入上式，得： （2·4） 由上两式得到： （2·5） 此式表明将模拟信号ha(t)的拉氏变换在s平面上沿虚轴按照周期延拓后，按照(2．3)式映射关系映射到z平面上，即(2．3)式可称为标准映射关系。 设计步骤 根据式子和由作业一求得的二阶VCVS网络的系统函数为，经计算得 式中T是采样间隔，若T选取过大，则会产生频率混叠现象。这里选取T=s，经过整理可得： 在matlab上画出H（z）的幅频特性和相频特性如下： （1）仿真程序如下： exp(-6500*10^(-7)) a1=-1.998,a2=0.987,b1=571.35; a=[1,a1,a2]; b=[0,b1]; figure,freqz(b,a) （2）仿真波形如下： 验证映射关系 设，，按照(2．3)式，得到： 因此得到：， 那么就有 上式关系说明，s平面的虚轴映射z平面的单位圆，s平面左半平面映射z平面单位圆内，s平面右半平面映射z平面单位圆外。这说明如果H。(s)因果稳定，转换后得到的H(z)也是因果稳定