用双线性变换法设计低通数字滤波器.docVIP

用双线性变换法设计低通数字滤波器 设计目的：利用双线性变换法实现模拟和数字滤波器幅度特性 设计内容： 设计低通数字滤波器，要求在通带内频率低于0.2πrad时，容许幅度误差在1dB以内；在频率0.3π到π之间的阻带衰减大于15dB。指定模拟滤波器采用巴特沃斯低通滤波器。试分别用双线性变换法设计数字低通滤波器。 设计原理： 一、通过模拟滤波器设计IIR数字滤波器 IIR数字滤波器在Z域中的传递函数为一有理分式，即 而模拟滤波器在S域中的传递函数为 可见H(z)与H(s)具有相同的形式，利用线性映照的方法，可以把S平面上的模拟滤波器映照成Z平面上的IIR数字滤波器。 按技术要求设计一个模拟滤波器，然后按一定的映照关系将中成数字滤波器的要求，必须对由复变量S到复变量Z直接的映照提出如下要求： ⑴因果温度的模拟滤波器转换成数字滤波器后仍是因果温度的，因此，映射应使S平面的左半平面映射为Z平面的单位圆内部，。 ⑵数字滤波器的幅频特性应与模拟滤波器的幅频特性一致，故S平面的虚轴线性映射到Z平面的单位圆上，即频率轴要对应。下图表明了上述映照关系 模拟滤波器到数字滤波器的转换可以在时域实现，也可以在频域实现。时域转换法是使数字滤波器的时域响应与模拟滤波器的时域采样值相等，具体方法有：冲激不变法、阶跃不变法和匹配Z变换法。频域变换法是使数字滤波器在范围内的幅度特性与模拟滤波器在范围内的幅度特性一致，具体方法有：双线性变换法，微分映照法。但上述方法中得到广泛的只有冲激响应不变法和双线性变换法。 二、用双线性变换法设计IIR低通滤波器 冲激响应不变法的主要缺点是会产生频率混叠现象，使数字滤波器的频响偏移模拟滤波器的频响。产生的原因是模拟低通的最高截止频率超过了折叠频率,在数字化后产生了频谱混叠，再通过标准映射关系，结果在附近形成了频谱混叠。为了克服一缺点，可以采用非线性频率压缩方法，将整个频率轴上的范围压缩到直接，再用转换到z平面上。设，经过非线性频率压缩到，表示，这里用正切变换实现频率压缩： 式中T仍是采样间隔，当从经过0变化到时，则由经过0变化到，实现了s平面上正虚轴完全压缩到平面上虚轴的直接的转换，这样便有 在通过转换到z平面上，得到 上式称为双性变换。从s平面到平面具有非线性频率压缩的功能，因此不可能产生频率混叠现象，这 是双线性变换法比脉冲响应不变法最大的优点。 优点：避免了频率响应的混迭现象s 平面与 z 平面为单值变换 缺点：除了零频率附近，与之间严重非线性线性相位模拟滤波器与非线性相位数字滤波器要求模拟滤波器的幅频响应为分段常数型，不然会产生畸变，预畸变。 设计步骤： 一、基本设置： 确定数字低通滤波器的技术指标：通带截止频率ωp、通 带衰减(p、阻带截止频率ωs、阻带衰减(s。 将数字低通滤波器的技术指标转换成模拟低通滤波器的 技术指标。 如果采用脉冲响应不变法，边界频率的转换关系为 如果采用双线性变换法，边界频率的转换关系为 (3) 按照模拟低通滤波器的技术指标设计模拟低通滤波器。 (4) 将模拟滤波器Ha(s)，从s平面转换到z平面，得到数字低通 滤波器系统函数H(z)。 二、用双线性变换法设计数字低通滤波器。 ① 数字低通技术指标仍为 ωp=0.2πrad, ( p=1dB; ωs=0.3πrad, ( s=15dB ② 模拟低通的技术指标为 ③ 设计巴特沃斯低通滤波器。阶数N计算如下： 取N=6。并求得：Ωc=0.7662rad/s 根据N=6，查表6.2.1得到的归一化传输函数Ha(p)与脉冲响应不变法得到的相同。为去归一化，将p=s/Ωc代入Ha(p)，得实际的Ha(s) ④ 用双线性变换法将Ha(s)转换成数字滤波器H(z) 本实验设计的模拟和数字滤波器幅度特性分别如图6.4.6 结果分析：由MATLAB仿真软件可以得到该实验的频谱特性如图所示： (一)、 （二）、 设计程序： % 用双线性变换法设计DF T=1; Fs=1/T; wpz=0.2;wsz=0.3; wp=2*tan(wpz*pi*/2); ws=2*tan(wsz*pi/2);rp=1;rs=1; %预畸变校正转换指标 ［N, wc］=buttord(wp, ws, rp, rs, ′s′); 　　　 %设计过渡模拟滤波器 ［B, A］=butter(N, wc, ′s′); ［Bz, Az］=bili