第19章模域频率变换法.pptVIP

一、实验目的　　(1)加深对模拟域频率变换法的了解。　　(2)掌握使用模拟低通滤波器原型进行频率变换及设计低通、高通、带通、带阻滤波器的方法。　　(3)了解MATLAB有关模拟域频率变换的子函数及其使用方法。 二、实验涉及的MATLAB子函数　　1.lp2lp　　功能：低通到低通模拟滤波器变换。这种变换是使用butter、cheby1、cheby2、ellip函数设计数字低通滤波器的一个步骤。　　调用格式：　　［bt，at］＝lp2lp(b，a，W0)；将传递函数表示的截止频率为1rad/s的模拟低通滤波器原型变换成截止频率为W0的低通滤波器。 　　［At，Bt，Ct，Dt］＝lp2lp(A，B，C，D，W0)；将连续状态方程表示的低通滤波器原型变换成截止频率为W0的低通滤波器。 　　2.lp2hp　　功能：低通到高通模拟滤波器变换。这种变换是使用butter、cheby1、cheby2、ellip函数设计数字高通滤波器的一个步骤。　　调用格式：　　［bt，at］＝lp2hp(b，a，W0)；将传递函数表示的截止频率为1rad/s的模拟低通滤波器原型变换成截止频率为W0的高通滤波器。　　［At，Bt，Ct，Dt］＝lp2hp(A，B，C，D，W0)；将连续状态方程表示的低通滤波器原型变换成截止频率为W0的高通滤波器。 　　3.lp2bp　　功能：低通到带通模拟滤波器变换。这种变换是使用butter、cheby1、cheby2、ellip函数设计数字带通滤波器的一个步骤。　　调用格式：　　［bt，at］＝lp2bp(b，a，W0，BW)；将传递函数表示的截止频率为1rad/s的模拟低通滤波器原型变换成中心频率为W0、带宽为BW的带通滤波器。　　［At，Bt，Ct，Dt］＝lp2bp(A，B，C，D，W0，BW)；将连续状态方程表示的低通滤波器原型变换成中心频率为W0、带宽为BW的带通滤波器。 　　如果已知被设计的滤波器低端截止频率为W1，高端截止频率为W2，则可以计算出W0和BW：　　W0＝sqrt(W1*W2)；BW＝W2－W1； 　　4.lp2bs　　功能：低通到带阻模拟滤波器变换。这种变换是使用butter、cheby1、cheby2、ellip函数设计数字带阻滤波器的一个步骤。　　调用格式：　　［bt，at］＝lp2bs(b，a，W0，BW)；将传递函数表示的截止频率为1rad/s的模拟低通滤波器原型变换成中心频率为W0、带宽为BW的带阻滤波器。 　　［At，Bt，Ct，Dt］＝lp2bs(A，B，C，D，W0，BW)；将连续状态方程表示的低通滤波器原型变换成中心频率为W0、带宽为BW的带阻滤波器。　　如果已知被设计的滤波器低端截止频率为W1，高端截止频率为W2，则可以计算出W0和BW：　　W0＝sqrt(W1*W2)；BW＝W2－W1； 　　5.set　　功能：设置图形对象属性。(本实验仅介绍设置坐标刻度的功能。)　　调用格式：　　set(gca，￠Xtick￠，xs，￠Ytick￠，ys)；二维坐标刻度设置。　　set(gca，￠Xtick￠，xs，￠Ytick￠，ys，￠Ztick￠，zs)；三维坐标刻度设置。　　xs，ys，zs可以是任何合法的实数向量，分别决定了x，y，z轴的刻度。 三、实验原理　　1.模拟域频率变换法　　IIR数字滤波器的设计，通常采用模拟域频率变换法和数字域频率变换法来实现。本实验以模拟域频率变换法为主，介绍MATLAB语言辅助设计的方法。　　模拟域频率变换法的基本设计思想是：先进行频率变换，后进行数字化变换。即在模拟低通滤波器原型设计好以后，通过频率变换，将模拟低通滤波器原型变换成实际的模拟低通、高通、带通、带阻滤波器。再通过相应的变换方法，将模拟滤波器变换成数字滤波器。 　　依靠MATLAB提供的相应子函数，采用模拟域频率变换法进行IIR数字滤波器设计的具体步骤如图19-1所示。从图中我们可以明确本实验要学习的模拟域频率变换在IIR数字滤波器的设计中所处的位置。 　　 图19-1 模拟域频率变换法设计IIR数字滤波器 　　2.由模拟滤波器原型设计模拟低通滤波器　　MATLAB提供的lp2lp子函数可用于模拟滤波器原型到实际的模拟低通滤波器的转换。　　例19-1 用频率变换法设计一个巴特沃斯模拟低通滤波器，要求(同例18-2)：通带截止频率fp＝2 kHz，通带最大衰减Rp≤1 dB，阻带截止频率fs＝5 kHz，阻带最小衰减As≥20 dB。　　解 编写程序如下：　　fp＝2000；Omgp＝2*pi*fp； %输入实际滤波器的通带截止频率 　　fs＝5000；Omgs＝2*pi*fs；%输入实际滤波器的阻带截