MATLAB在信号处理中的应用13618.pptVIP

4.8.3 IIR数字滤波器的频率变换设计法 1. IIR数字滤波器的频率变换设计法的基本原理根据滤波器设计要求，设计模拟原型低通滤波器，然后进行频率变换，将其转换为相应的模拟滤波器高通、带通等，最后利用冲激响应不变法或双线性变换法，将模拟滤波器数字化成相应的数字滤波器。 图4.39 IIR数字滤波器MATLAB设计步骤流程图 1．MATLAB的典型设计 利用在MATLAB设计IIR数字滤波器可分以下几步来实现 1按一定规则将数字滤波器的技术指标转换为模拟低通滤波器的技术指标； 2根据转换后的技术指标使用滤波器阶数函数，确定滤波器的最小阶数N和截止频率Wc； 3利用最小阶数N产生模拟低通滤波原型； 4利用截止频率Wc把模拟低通滤波器原型转换成模拟低通、高通、带通或带阻滤波器； 5利用冲激响应不变法或双线性不变法把模拟滤波器转换成数字滤波器。 [例4-39] 设计一个数字信号处理系统，它的采样率为Fs＝100Hz，希望在该系统中设计一个Butterworth型高通数字滤波器，使其通带中允许的最小衰减为0.5dB，阻带内的最小衰减为40dB，通带上限临界频率为30Hz，阻带下限临界频率为40Hz。 MATLAB源程序设计如下：%把数字滤波器的频率特征转换成模拟滤波器的频率特征wp30*2*pi;ws40*2*pi;rp0.5;rs40;Fs100;[N,Wc]buttordwp,ws,rp,rs,s;%选择滤波器的最小阶数[Z,P,K]buttapN;%创建Butterworth低通滤波器原型[A,B,C,D]zp2ssZ,P,K;%零－极点增益模型转换为状态空间模型[AT,BT,CT,DT]lp2hpA,B,C,D,Wc;　 %实现低通向高通的转变[num1,den1]ss2tfAT,BT,CT,DT;%状态空间模型转换为传递函数模型%运用双线性变换法把模拟滤波器转换成数字滤波器[num2,den2]bilinearnum1,den1,100;[H,W]freqznum2,den2;%求频率响应plotW*Fs/2*pi,absH;grid;%绘出频率响应曲线xlabel频率/Hz; ylabel幅值 程序运行结果如图4.40所示。 2．MATLAB的直接设计 图4.39 IIR数字滤波器MATLAB设计步骤流程图 [例4-41] 试设计一个带阻IIR数字滤波器，其具体的要求是：通带的截止频率：wp1＝650Hz、wp2＝850Hz；阻带的截止频率：ws1＝700Hz、ws2＝800Hz；通带内的最大衰减为rp＝0.1dB；阻带内的最小衰减为rs＝50dB；采样频率为Fs＝2000Hz。MATLAB源程序设计如下： wp1650;wp2850;ws1700;ws2800;rp0.1;rs50;Fs2000; wp[wp1,wp2]/Fs/2;ws[ws1,ws2]/Fs/2;　 %利用Nyquist频率频率归一化 [N,wc]ellipordwp,ws,rp,rs,z;%求滤波器阶数 [num,den]ellipN,rp,rs,wc,stop;%求滤波器传递函数 [H,W]freqznum,den;%绘出频率响应曲线 plotW*Fs/2*pi,absH;grid; xlabel频率/Hz;ylabel幅值 该程序运行后的幅频响应曲线如图4.42所示。 * 第 4 章 MATLAB在信号处理中的应用 第4章 MATLAB在信号处理中的应用 4.1 信号及其表示 4.2 信号的基本运算 4.3 信号的能量和功率 4.4 线性时不变系统 4.5 线性时不变系统的响应 4.6 线性时不变系统的频率响应 4.7 傅里叶Fourier变换 4.8 IIR数字滤波器的设计方法 4.9 FIR数字滤波器设计 4.1 信号及其表示 4.1.1连续时间信号的表示连续时间信号：时间变化连续。如yxt 离散时间信号序列：时间离散，如xnTxt|tnT. 4.1.2工具箱中的信号产生函数 电压控制振荡器 vco 产生高斯单脉冲信号 gmonopuls 产生高斯正弦脉冲信号 gauspuls 产生Dirichlet或周期sinc函数 diric 产生调频余弦信号 chirp 产生非周期的三角波信号 tripuls 产生sinc函数波形 sinc 产生非周期的方波信号 rectpule 产生方波信号 square 产生冲激串 pulstran 产生锯齿波或三角波信号 sawtooth 功能 函数名 功能 函数名 4.1.3离散时间信号的表示在MATLAB中，离散时间信号xn的表示：需用一个向量x表示序列幅值，用另一个等长的定位时间变量n，才能完整地表示一