MATLAB在数字信号处理中的应用.pdf

第4章 MATLAB在数字信号处理中的应用 MATLAB 中的信号处理工具箱提供了很多 数字信号处理中需要用到的函数及解决方法。 在数字信号处理中所学习到的问题，诸如滤波 器的设计、自适应滤波，维纳滤波、卡尔曼滤 波等理论都可以通过MATLAB仿真得到实现和验 证。本章提供了12个MATLAB在数字信号处理课 程学习中的应用实例，可作为学习数字信号处 理课程的参考。 第4章 MATLAB在数字信号处理中的应用 知 识 架 构 第4章 MATLAB在数字信号处理中的应用 4.1 IIR 带通滤波器设计 1．实验目的 (1)掌握脉冲响应不变法设计IIR数字滤波器的具体 设计方法。 (2)熟悉脉冲响应不变法设计低通滤波器的仿真。 2．实验原理 脉冲响应不变法是从滤波器的脉冲响应出发， 使数字滤波器的单位脉冲响应序列h(n)模仿模拟滤 波器的冲击响应ha(t)，使h(n)正好等于ha(t)的采 样值，即： 第4章 MATLAB在数字信号处理中的应用 4.1 IIR 带通滤波器设计 T为采样周期。若以Ha(s)及H(z)分别表示ha(t) 的拉氏变换及h(n)的Z变换，即： 根据采样序列Z变换与模拟信号拉氏变换的关系得: 上式表明，采用脉冲响应不变法将模拟滤波器 变换为数字滤波器时，它所完成的S平面到Z平面的 变换，正是以前讨论的拉氏变换到Z变换的标准变 换关系。 第4章 MATLAB在数字信号处理中的应用 4.1 IIR 带通滤波器设计 脉冲响应不变法特别适用于用部分分式表达的 传递函数，模拟滤波器的传递函数若只有单阶极 点，且分母的阶数高于分子阶数，则可表达为部分 分式形式： 其拉氏反变换为： 对ha(t)采样就得到数字滤波器的单位脉冲响应序 列： 再对h(n)取Z变换，得到数字滤波器的传递函数： 第4章 MATLAB在数字信号处理中的应用 4.1 IIR 带通滤波器设计 第二个求和为等比级数之和： 要收敛的话，必有： 所以有： 比较部分分式形式的Ha(s)和上式H(z)可以看 到，把S平面上的极点变换到Z平面上对应的极点， 而Ha(s)与H(z)中部分分式所对应的系数不变。如 果模拟滤波器是稳定的，则所有极点都在S左半平 面那么变换后H(z)的极点也都在单位圆内，因此数 字滤波器保持稳定。 第4章 MATLAB在数字信号处理中的应用 4.1 IIR 带通滤波器设计 3．仿真思路 在MATLAB中，可以用下列函数辅助设计IIR数字 滤波器。 (1)利用cheb1ord可以确定低通原型和切比雪夫滤波 器的阶数和截止频率。 (2) [num，den]=cheby1(N，Wn) ，[num，den]=cheby2(N，Wn) (切比雪夫1 型和2 型)可以进行滤波器的设计。 (3)利用impinvar可以完成脉冲响应不变法的模拟滤波 器到数字滤波器的转换。 4．程序代码 第4章 MATLAB在数字信号处理中的应用 4.1 IIR 带通滤波器设计 clear wp=6*pi*10^3;ws=9*pi*10^3;ap=1,as=15; Fs=30*10^3; wp1=wp/Fs;ws1=ws/Fs; [N,WC]=cheb1ord(wp,ws,ap,as,s); [b,a]=cheby1(N,ap,WC,‘s’) ； [bz,az]=impinvar(b,a,Fs)； w0=[wp1,ws1] ； Hx=freqz(bz,az,w0) ； [H,W]=freqz(bz,az)； dbHx=-20*log10(abs(Hx)/max(abs(H))) ；plot(W,abs(H)); xlabel(相对频率);ylabel(幅频); grid 第4章 MATLAB在数字信号处理中的应用 4.1 IIR 带通滤波器设计 5．运行结果与分析 结论：由图4.2可知，切比雪夫Ⅰ型滤波器的振 幅特性在通带内是等波纹的，在阻带内是单调的。 图 4.2： 6．思考题 (1) 带通滤波器的幅度有什么特点？图4.2中滤波器的下降沿 怎么样才能更陡峭？ (2