数字信号处理--用双线性变换法设计IIR数字滤波器要点.pdf
《数字信号处理--用双线性变换法设计IIR数字滤波器要点.pdf》由会员分享,可在线阅读,更多相关《数字信号处理--用双线性变换法设计IIR数字滤波器要点.pdf(12页珍藏版)》请在三一文库上搜索。
1、1 数字信号处理实验内容 用双线性变换法设计 IIR 数字滤波器 一、实验目的: 1、熟悉用双线性变换法设计IIR 数字滤波器的原理与方法。 2、掌握数字滤波器的计算机仿真方法。 3、熟悉 Batterworth滤波器设计方法及特点 二、实验原理: (一) 、IIR 数字滤波器的设计步骤: 按照一定规则把给定的滤波器技术指标转换为模拟低通滤波器的技术指标; 根据模拟滤波器技术指标设计为响应的模拟低通滤波器; 跟据脉冲响应不变法和双线性不变法把模拟滤波器转换为数字滤波器; 如果要设计的滤波器是高通、带通或带阻滤波器,则首先把它们的技术指标转化为 模拟低通滤波器的技术指标,设计为数字低通滤波器,最
2、后通过频率转换的方法来得到所要 的滤波器。 在 MATLAB 中,经典法设计IIR数字滤波器主要采用以下步骤: IIR 数字滤波器设计步骤 (二)、用模拟滤波器设计数字滤波器的方法 1、冲激响应不变法: 冲激响应不变法是从时域出发,要求数字滤波器的冲激响应h(n) 对应于模拟滤波器h(t) 的 等间隔抽样。 优点:时域逼近良好;保持线性关系。 缺点:频域响应混叠。只适用于限带低通滤波器和带通滤波器 模拟滤波器原型 buttap,cheb1ap 频率变换 模拟离散化 bilinear,impin var IIR 数字滤波器 / T / T 3/T 3/T j jIm(z) Re(z) 1 S 平
3、面Z 平面 2 2、双线性变换法 优点:克服了频域混叠 缺点:高频时会引起畸变 1)冲激响应不变法impinvar 格式: BZ,AZ= impinvar(B,A,Fs ) 功能:把具有 B,A 模拟滤波器传递函数模型转换为采样频率为Fs 的数字滤波器的传递函数 模型 BZ,AZ ,Fs 默认值为1。 例:一个4 阶的 Butterworth模拟低通滤波器的系统函数如下: 1225 1 )( 234 ssss sH a 试用冲激响应不变法求出Butterworth模拟低通数字滤波器的系统函数。 num=1; den=1,sqrt(5),2,sqrt(2),1; num1,den1=impinv
4、ar(num,den) 2)双线性变换法bilinear 格式一: Zd,Pd,Kd= bilinear(Z,P,K,Fs) 功能:把模拟滤波器的零极点模型转换成数字滤波器的零极点模型,Fs 是采样频率 格式二: numd,dend= bilinear(num,den,Fs) 功能:把模拟滤波器的传递函数模型转换为数字滤波器的传递函数模型。 例:一个三阶的模拟Butterworth模拟低通滤波器的系统函数如下: 123 1 )( 23 sss sH, 试用双线性变换法求出数字Butterworth数字低通滤波器的 系统函数。 num=1; den=1,sqrt(3),sqrt(2),1; nu
5、m1,den1=bilinear(num,den,1) / T / T j jIm(z) Re(z) 1 S 平面 Z平面 j S 平面 1 1 1 1 SS TT 将整个平面压缩变换到平面一个的带状区域 3 3) IIR 数字滤波器的频率变换实现 步骤: 按一定的规则将数字滤波器的技术指标转换为模拟低通滤波器的技术指标 根据转换后的技术指标使用滤波器阶数函数,确定滤波器的最小阶数N和截止频率Wc 利用最小阶数N产生模拟低通原型 利用截止频率Wc把模拟低通滤波器原型转换为模拟低通、高通、带通、带阻滤波器 利用冲激响应不变法或双线性变换法把模拟滤波器转换为数字滤波器 表一IIR 滤波器阶次估计
6、函数名功能说明 buttord 计算 Butterworth 滤波器的阶次及截止频率 cheb1ord 计算 Chebyshev滤波器的阶次 cheb2ord 计算 Chebyshev滤波器的阶次 ellipord 计算椭圆滤波器的最小阶次 表二模拟低通滤波器原型设计 函数名功能说明 buttap Butterworth 模拟低通滤波器原型设计 cleb1ap Chebyshev模拟低通滤波器原型设计 cheb2ap Chebyshev模拟低通滤波器原型设计 ellipap 椭圆模拟低通滤波器原型设计 表三模拟滤波器变换函数 函数名功能说明 lp2bp 模拟低通转换为带通 lp2bs 模拟低通
7、转换为带阻 lp2hp 模拟低通转换为高通 lp2lp 改变模拟低通的截止频率 (三)、数字滤波器的设计 1 数字滤波器的设计参数 滤波器的4 个重要的通带、阻带参数为: p f:通带截止频率(Hz) s f:阻带起始频率(Hz) p R:通带内波动(dB) ,即通带内所允许的最大衰减; s R:阻带内最小衰减 设采样速率 (即奈奎斯特速率)为 N f,将上述参数中的频率参数转化为归一化角频率参数: p:通带截止角频率( rad/s) ,)2/( Npp ff; s:阻带起始角频率( rad/s) ,)2/( Nss ff 通过以上参数就可以进行离散滤波器的设计。 2、巴特沃斯滤波器设计 4
8、1) 巴特沃斯滤波器阶数的选择: 在已知设计参数 p,s,p R, s R之后,可利用“buttord”命令可求出所需要的滤 波器的阶数和3dB 截止频率,其格式为: n,Wn=buttordWp ,Ws, Rp,Rs,其中Wp,Ws,Rp, Rs 分别为通带截止频率、 阻带起始频率、通带内波动、阻带内最小衰减。返回值n 为滤波器的最低阶数,Wn 为 3dB 截止频率。 2) 巴特沃斯滤波器系数计算: 由巴特沃斯滤波器的阶数n 以及 3dB 截止频率Wn 可以计算出对应传递函数H(z)的分 子分母系数, MATLAB提供的命令如下: (a)巴特沃斯低通滤波器系数计算: b,a=butter(n
9、,Wn) ,其中 b 为 H(z)的分子多项式系数,a 为 H(z)的分母多项式系数 (b)巴特沃斯高通滤波器系数计算:b,a=butter(n ,Wn, High ) (c)巴特沃斯带通滤波器系数计算:b,a=butter(n ,W1 ,W2) ,其中 W1,W2 为截止频率, 是 2 元向量,需要注意的是该函数返回的是2*n 阶滤波器系数。 (d)巴特沃斯带阻滤波器系数计算:b,a=butter(ceil(n/2) , W1 ,W2 , stop ),其中 W1 , W2 为截止频率,是2 元向量,需要注意的是该函数返回的也是2*n 阶滤波器系数。 三、巴特沃斯滤波器设计实例: 例题 1:
10、 采样速率为8000Hz, 要求设计一个低通滤波器, p f=2100Hz , s f=2500Hz, p R=3dB, s R=25dB。 用直接设计法 程序如下: fn=8000;% 采样频率 fp=2100; %通带截止频率 fs=2500; % 阻带起始频率 Rp=3; %通带最大衰减 Rs=25;%阻带最小衰减 Wp=fp/(fn/2);% 计算归一化角频率 Ws=fs/(fn/2); n,Wn=buttord(Wp,Ws,Rp,Rs);%计算阶数和截止频率 b,a=butter(n,Wn);% 计算 H(z)分子、分母多项式系数 H,F=freqz(b,a,1000,8000);%
11、 计算 H(z)的幅频响应 ,freqz(b,a,计算点数 ,采样速率 ) subplot(2,1,1) plot(F,20*log10(abs(H) % 画出幅频特性图 xlabel(Frequency(Hz); ylabel(Magnitude(dB) title( 低通滤波器 ) axis(0 4000 -30 3);grid on 5 pha=angle(H)*180/pi; subplot(2,1,2) plot(F,pha);grid on % 画出相频特性图 xlabel(Frequency(Hz); ylabel(phase); 050010001500200025003000
12、35004000 -30 -20 -10 0 Frequency(Hz) M a g n i t u d e (d B ) 低通滤波器 05001000150020002500300035004000 -200 -100 0 100 200 Frequency(Hz) p h a s e 用双线性变换法 wp=2100*2*pi; %利用 ws=2500*2*pi; Rp=3; Rs=25; Fs=8000; Ts=1/Fs; % 选择滤波器的最小阶数 N,Wn=buttord(wp,ws,Rp,Rs,s);%创建 butterworth模拟滤波器 Z,P,K=buttap(N); % 把滤波
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 数字信号 处理 双线 变换 设计 IIR 数字滤波器 要点
链接地址:https://www.31doc.com/p-5208686.html