数字信号处理课程设计-用双线性变换法设计高通和带通数字滤波器

实验三 用双线性变换法设计IIR 滤波器一、实验目的1、了解两种工程上最常用的变换方法：脉冲响应不变法和双线性变换法。

2、掌握双线性变换法设计IIR 滤波器的原理及具体设计方法，熟悉双线性法设计低通、带通和高通IIR 滤波器的计算机程序。

3、观察用双线性变换法设计的滤波器的品与特性，并与脉冲响应不变法相比较，了解双线性变换法的特点。

4、熟悉用双线性变换法涉及数字Butterworth 和Chebyshev 滤波器的全过程。

5、了解多项式呈几何多项式乘方运算的计算机变成方法。

二、实验原理与方法1、确定数字滤波器的性能指标。

这些指标包括：通带、阻带临界频率pf 、s f ；通带内的最大衰减pα；阻带内的最小衰减s α；采样周期T 。

2、确定相应的数字频率，Tf T f s s p pπωπω2,2==。

3、计算经过频率预畸的相应参考模拟低通原型的频率)2(),2(s s p p tg tg ωω=Ω=Ω。

4、计算低通原型阶数N ，计算3db 归一化频率c Ω，从而求得低通原型的传递函数)(s H a 。

5、用上表中所列变换公式1111--+-=zz s ，代入)(s H a ，求得数字滤波器传递函数)(z H =1111|)(--+-=zz s a s H 。

6、分析滤波器频域特性，检查其指标是否满足要求。

三、实验内容及步骤1、采样频率为1HZ，设计一个Chebyshev高通数字滤波器其中通带临界频率fp=0.3HZ,通带内衰减小于0.8db(α=0.8db),阻带临界频率fs=0.2HZ,p阻带内衰减大于20db(α=20db)。

求这个数字滤波器的传递函数H(z)，输s出它的幅频特性，观察其通带衰减和阻带衰减是否满足要求。

高通数字滤波器的设计f=1; fp=3/10; fs=2/10; Rp=0.8; Rs=20;[n,Wn]=cheb1ord(2*fp/f,2*fs/f,Rp,Rs);[b,a]=cheby1(n,Rp,Wn,'high');freqz(b,a,512,1)b =0.0262 -0.1047 0.1570 -0.1047 0.0262a =1.0000 1.5289 1.6537 0.9452 0.2796n =4 Wn =0.60002、采样频率为1HZ ，设计一个数字低通滤波器，要求其通带临界频率fp=0.2HZ ，通带内衰减小于1db(dbp1=α),阻带临界频率fs=0.3HZ,阻带内衰减25db(dbs 25=α)。

目录一、摘要 (3)二、设计思想 (3)2.1 IIR数字滤波器设计思路 (3)2.2设计IIR数字滤波器的两种方法 (4)2.3双线性变换法的基本原理 (5)2.4用双线性变换法设计IIR数字滤波器的步骤 (6)三、程序源代码和运行结果 (6)3.1低通滤波器 (6)3.2高通滤波器 (9)3.3带通滤波器 (12)3.4带阻滤波器 (14)四、网络结构图 (17)五、心得体会 (19)参考文献 (19)一、摘要数字滤波器是具有一定传输选择特性的数字信号处理装置,其输入、输出均为数字信号,实质上是一个由有限精度算法实现的线性时不变离散系统。

它的基本工作原理是利用离散系统特性对系统输入信号进行加工和变换,改变输入序列的频谱或信号波形,让有用频率的信号分量通过,抑制无用的信号分量输出。

数字滤波器和模拟滤波器有着相同的滤波概念,根据其频率响应特性可分为低通、高通、带通、带阻等类型,与模拟滤波器相比,数字滤波器除了具有数字信号处理的固有优点外,还有滤波精度高(与系统字长有关)、稳定性好(仅运行在0与l两个电平状态)、灵活性强等优点。

数字滤波器按单位脉冲响应的性质可分为无限长单位脉冲响应滤波器IIR和有限长单位脉冲响应滤波器(FIR)两种。

本文介绍IIR数字滤波器的设计。

二、设计思想2.1 IIR数字滤波器设计思路IIR 数字滤波器可用一个n阶差分方程表示y(n)=Σb r x(n-r)+Σa k y(n-k)或用它的Z域系统函数:对照模拟滤波器的传递函数:不难看出,数字滤波器与模拟滤波器的设计思路相仿,其设计实质也是寻找一组系数{b,a},去逼近所要求的频率响应,使其在性能上满足预定的技术要求;不同的是模拟滤波器的设计是在S平面上用数学逼近法去寻找近似的所需特性H(S),而数字滤波器则是在Z平面寻找合适的H(z)。

IIR数字滤波器的单位响应是无限长的,而模拟滤波器一般都具有无限长的单位脉冲响应,因此与模拟滤波器相匹配。

双线性变换法在数字滤波器上应用冲击响应不变法是使数字滤波器在时域上模仿模拟滤波器的冲击响应，但是它的缺点是产生频响的混叠失真，这是因为从s平面到z平面不是一一映射关系。

为了克服这个缺点，可以采用双线性变换法。

既然冲击响应不变法是将一条横带变换到整个z平面上去，可以将s平面压缩变换到某一中介s1平面的一条横带里，再通过标准变换关系z=exp(s1*T)将此带变换到整个z平面上去，这样就使s平面与z平面之间建立一一对应的单值关系，消除了多值变换性。

为了将s平面的jΩ轴压缩到s1平面的jΩ1轴上的-pi/T到pi/T一段上，可以通过以下的正切变换来实现：Ω=c*tan(Ω1*T/2)这里c是任意常数。

这样当Ω1由-pi/T经0变化到pi/T时，Ω由-∞经过0变化到+∞，也映射到了整个jΩ轴。

将这个关系延拓到整个s平面和s1平面，则可以得到s=c*tanh(s1*T/2)=c*[1-exp(-s1*T)] /[1+exp(-s1*T)]再将s1平面通过标准变换关系映射到z平面，即令z=exp(s1*T)，并且，通常取c=2/T，得到s=(2/T) * [1-z^(-1)] / [1+z^(-1)]同样对z求解，得到z=[1+((2/T)*s] / [1-((2/T)*s]这样的变换叫做双线性变换。

为了验证这种映射具有s平面的虚轴映射到z平面单位圆上的特性，考虑z=exp(j*ω)，得s=(2/T) * [1-exp(-j*ω)] / [1+exp(-j*ω)]=(2/T) * j*[sin(ω/2)]/ cos(ω/2)]=(2/T) * j * tan(ω/2)=jΩ除了使s平面的虚轴映射到单位圆上之外，s平面的左半部分映射到单位圆的内部，s平面的右半部分映射到单位圆的外部。

观察式子z=[1+((2/T)*s] / [1-((2/T)*s]，发现s的实部为负时，因子[1+((2/T)*s] / [1-((2/T)*s]的幅度小于1，相当于单位圆的内部。

数字信号处理课程设计--双线性变换法设计数字高通滤器《数字信号处理》课程设计报告题目双线性变换法设计数字高通滤器学院信息工程学院专业通信工程班级1004学号2010013598 2010013643 学生姓名陈涛周亮指导教师刘振二0一二年十二月目录一前言 (1)二课程设计的目的和意义 (2)三课程设计题目描述及要求 (2)四详细设计过程 (2)1.设计思想 (2)2. 设计原理 (3)2.1 巴特沃斯型滤波器原理 (3)2.2 双线性变换法原理 (6)3. 设计过程 (7)3.1巴特沃斯模拟滤波器的设计 (7)3.2巴特沃斯模拟滤波器转变为数字高通滤波器 (7)五总结 (10)1.设计过程中遇到的问题及解决办法 (10)2.总结与体会 (10)参考文献 (11)一前言数字信号处理的目的是对真实世界的连续模拟信号进行测量或滤波。

因此在进行数字信号处理之前需要将信号从模拟域转换到数字域，这通常通过模数转换器实现。

而数字信号处理的输出经常也要变换到模拟域，这是通过数模转换器实现的。

数字信号处理的算法需要利用计算机或专用处理设备如数字信号处理器（DSP）和专用集成电路（ASIC）等。

数字信号处理技术及设备具有灵活、精确、抗干扰强、设备尺寸小、造价低、速度快等突出优点，这些都是模拟信号处理技术与设备所无法比拟的。

巴特沃斯滤波器是电子滤波器的一种。

巴特沃斯滤波器的特点是通频带的频率响应曲线最平滑。

这种滤波器最先由英国工程师斯替芬·巴特沃斯（Stephen Butterworth）在1930年发表在英国《无线电工程》期刊的一篇论文中提出的。

数字滤波器是具有一定传输选择特性的数字信号处理装置,其输入、输出均为数字信号,实质上是一个由有限精度算法实现的线性时不变离散系统。

此课程设计主要介绍了用双线性不变法设计数字高通滤波器，要求采用巴特沃斯型滤波器，介绍了设计步骤，然后在MATLAB 环境下进行了仿真与调试，实现了设计目标。

电气与信息工程学院数字信号处理实验报告学生姓名班级电子信息工程学号指导教师2019.12实验四 IIR 数字滤波器的设计一、实验目的：1. 掌握双线性变换法及脉冲响应不变法设计IIR 数字滤波器的具体设计方法及其原理，熟悉用双线性变换法及脉冲响应不变法设计低通、高通和带通IIR 数字滤波器的MATLAB 编程。

2. 观察双线性变换及脉冲响应不变法设计的滤波器的频域特性，了解双线性变换法及脉冲响应不变法的特点。

3. 熟悉巴特沃思滤波器、切比雪夫滤波器和椭圆滤波器的频率特性。

二、实验原理:1． 脉冲响应不变法用数字滤波器的单位脉冲响应序列 模仿模拟滤波器的冲激响应 ,让 正好等于 的采样值，即 ，其中 为采样间隔，如果以 及 分别表示 的拉式变换及 的Z 变换，则)2(1)(m T j s H T z H m a e z sT ∑∞-∞==+=π2．双线性变换法S 平面与z 平面之间满足以下映射关系：);(,2121,11211ωωσj re z j s s T s T z z z T s =+=-+=+-⋅=-- s 平面的虚轴单值地映射于z 平面的单位圆上，s 平面的左半平面完全映射到z 平面的单位圆内。

双线性变换不存在混叠问题。

双线性变换是一种非线性变换，这种非线性引起的幅频特性畸变可通过预畸而得到校正。

三、实验内容及步骤：实验中有关变量的定义:fc 通带边界频率；fr阻带边界频率；δ通带波动；At 最小阻带衰减；fs采样频率；T采样周期上机实验内容：（1）fc=0.3KHz,δ=0.8dB,fr=0.2KHz, At=20dB,T=1ms;设计一个切比雪夫高通滤波器，观察其通带损耗和阻带衰减是否满足要求。

MATLAB源程序：wp=2*1000*tan(2*pi*300/(2*1000));ws=2*1000*tan(2*pi*200/(2*1000));[N,wn]=cheb1ord(wp,ws,0.8,20,'s'); %给定通带（wp）和阻带(ws)边界角频率，通带波动波动0.8，阻带最小衰减20dB，求出最低阶数和通带滤波器的通带边界频率Wn[B,A]=cheby1(N,0.5,wn,'high','s');%给定通带（wp）和阻带(ws)边界角频率，通带波动[num,den]=bilinear(B,A,1000);[h,w]=freqz(num,den);f=w/(2*pi)*1000;plot(f,20*log10(abs(h)));axis([0,500,-80,10]);grid;xlabel('频率');ylabel('幅度/dB')程序结果num = 0.0304 -0.1218 0.1827 -0.1218 0.0304 den = 1 1.3834 1.4721 0.8012 0.2286系统函数：123412340.0304 -0.1218z 0.1827z-0.1218z0.0304z H(z)=1.0000+1.3834z+1.4721z+ 0.8012z+0.2286z--------++幅频响应图：分析：由图可知，切比雪夫滤波器幅频响应是通带波纹，阻带单调衰减的。

西安邮电大学数字信号处理课内上机报告专业班级： 学生姓名： 学号(班内序号):年 月 日——————————————————————————装订线————————————————————————————————报告份数：题目：设计一个数字带通滤波器，通带范围为0.3pi rad到0.4pi，通带最大衰减为3db，阻带最小衰减为18db，0.2pi 以下和0.5pi以上范围为阻带。

采用切比雪夫型，利用双线性变换法设计之源程序％所设计的数字滤波器的指标Ts = 0.1，Fs=1/Ts，Rp = 3，Rs = 18;wp1=0.3*pi，wp2=0.4*pi;ws1=0.2*pi，ws2=0.5*pi;％频率的预畸变Wp1=(2/T)*tan(wp1/2);Wp2=(2/T)*tan(wp2/2);Wp=[Wp1,Wp2]; ％模拟滤波器的通带截止频率Ws1=(2/T)*tan(ws1/2);Ws2=(2/T)*tan(ws2/2);Ws=[Ws1,Ws2]; ％模拟滤波器的阻带截止频率BW=Ws2-Ws1; ％模拟滤波器的带宽% BW=Wp2-Wp1;Omegaw0=sqrt(Ws1*Ws2); ％模拟滤波器的中心频率% Omegaw0=sqrt(Wp1*Wp2);％求模拟低通滤波器的阶数与边缘频率[N,OmegaC]=cheb2ord(Wp,Ws,Rp,Rs,'s');%[N,OmegaC]=cheb1ord(Wp,Ws,Rp,Rs,'s')％求切比雪夫模拟低通滤波器的零、极点与增益[z0,p0,k0]=cheb2ap(N,Rs);%[z0,p0,k0]=cheb1ap(N,Rp);%设计归一化的模拟原型带通滤波器%求原型滤波器的分子系数AnalogB=k0*real(poly(z0));%求原型滤波器的分母系数AnalogA=real(poly(p0));%模拟低通到模拟带通的分子、分母系数的变换[BandB,BandA]=lp2bp(AnalogB,AnalogA,Omegaw0,BW);%双线性变换：模拟带通与数字带通的分子分母系数的变换[DigitalB,DigitalA] = bilinear(BandB,BandA,Fs);%变为二阶节级联结构[sos,G] = tf2sos(DigitalB,DigitalA);%求数字带通滤波器的幅频、相频特性、及其群延迟％求数字带通滤波器的幅频特性[Hz,Wz]=freqz(DigitalB,DigitalA,1024,'whole');％将数字带通滤波器的幅频特性转化为分贝表示dbHz=20*log10((abs(Hz)+eps)/max(abs(Hz)));％求数字带通滤波器的相频特性φ＝angle(Hz)％求数字带通滤波器的群延迟特性grd = grpdelay(DigitalB,DigitalA,Wz);％作图subplot(2,3,1);plot(Wz/pi,abs(Hz));title('幅频响应');xlabel(''),ylabel('幅度:|Hz|');axis([0,1,0,1.1]);set(gca,'XTickMode','manual','XTick',[0,ws1/pi,wp1/pi,wp2/pi,ws2/pi,1]);grid; subplot(2,3,4);plot(Wz/pi,dbHz);title('模值(dB)');xlabel('频率（单位：\pi）');ylabel('分贝（dB）');axis([0,1,-40,5]);set(gca,'XTickMode','manual','XTick',[0,ws1/pi,wp1/pi,wp2/pi,ws2/pi,1]);set(gca,'YTickmode','manual','YTick',[-50,-30,-2,0]);set(gca,'YTickLabelMode','manual','YTickLabels',['50';'30';' 2';' 0']);grid; subplot(2,3,2);plot(Wz/pi,angle(Hz)/pi);title('相频响应');xlabel('');ylabel('单位：\pi ');axis([0,1,-1,1]);set(gca,'XTickMode','manual','XTick',[0,ws1/pi,wp1/pi,wp2/pi,ws2/pi,1]);grid; subplot(2,3,5);title(零极点图');ylabel('单位：\dB');xlabel('单位：\pi ');zplane(DigitalB,DigitalA);axis([-1.1,1.1,-1.1,1.1]);subplot(2,3,3); plot(Wz/pi,grd); title('群延迟')xlabel('频率（单位：\pi）'); ylabel('样本'); axis([0,1,0,8])set(gca,'XTickMode','manual','XTick',[0,ws1/pi,wp1/pi,wp2/pi,ws2/pi,1]); %画高刻度线set(gca,'YTickmode','manual','YTick',[0:0.5:10]); gridset(gcf,'color','w') % 置图形背景色为白色幅频响应幅度:|H z|模值(dB)频率：（单位：π）分贝（d B）相频响应单位:πReal PartI m a g i n a r y P a r t群延迟频率（单位：pi ）样本Chebyshev II 型IIR 数字带通滤波器幅频响应幅度:|H z|模值(dB)频率：（单位：π）分贝（d B）相频响应单位:πReal PartI m a g i n a r y P a r t群延迟频率（单位：pi ）样本Chebyshev I 型 IIR 带通滤波器指导教师评语：实 验 成 绩: 指导(辅导)教师 :——————————————————————————装订线————————————————————————————————。

《数字信号处理》教学大纲课程类型：专业课总学时：通信工程专业70；信息工程专业64讲课学时：通信工程专业60；信息工程专业54实践学时：通信工程专业10；信息工程专业10一、课程的目的与任务本课程讲授数字信号处理的基本理论和基本分析方法，并且进行理论与算法的实践。

要求学生掌握离散时间信号与系统的基本理论，掌握离散时间系统的时域分析与 Z变换及离散傅立叶变换和快速傅里叶变换的理论计算法；掌握IIR和FIR数字滤波器的结构、理论和设计方法，为学生毕业后从事数字技术及其工程应用提供必要的训练。

二、课程有关说明《数字信号处理》是通信工程专业和信息工程专业的专业课，课程的内容包括：线性时不变离散时间系统的基础知识、数学模型（差分方程）及其求解，Z变换，离散傅立叶变换（DFT）理论及应用，快速傅立叶变换（FFT），无限长单位脉冲响应（IIR）数字滤波器设计，有限长单位脉冲响应（FIR）数字滤波器设计等内容。

除了理论教学外，还配有一定数量的上机实验。

数字信号处理在理论上所涉及的范围及其广泛。

高等数学、随机过程、复变函数等都是其数学基本工具。

电路理论、信号与系统等是其理论基础。

其算法及实现（硬件和软件）与计算机学科和微电子技术密不可分。

学生应该认真学习以上的知识，更好地掌握数字信号处理的基本理论、算法和实现技能。

主要教学方式：教师主讲，答疑、课堂讨论为辅，并结合实验教学。

考核评分方式：闭卷考试三、教学内容绪论（2学时）本章应掌握：数字信号处理的基本概念。

熟悉：数字信号处理系统的基本组成。

了解：数字信号处理的学科概貌、学科特点、实际应用、发展方向和实现方法。

第一章时域离散信号和时域离散系统（4学时）第一节时域离散信号本节应掌握：序列的运算，即移位、翻褶、和、积、累加、差分、时间尺度变换、卷积和等；序列的周期性。

熟悉：几种常用序列，即单位抽样序列、单位阶跃序列、矩形序列、实指数序列、复指数序列、正弦序列。

了解：用单位抽样序列来表示任意序列。

实验五、双线性变换法设计IIR 数字滤波器一、实验目的：1、熟悉用双线性变换法设计IIR 数字滤波器的原理与方法。

2、掌握数字滤波器的计算机仿真方法。

3、熟悉Batterworth 滤波器设计方法及特点 二、实验原理（一）、IIR 数字滤波器的设计步骤：① 按照一定规则把给定的滤波器技术指标转换为模拟低通滤波器的技术指标； ② 根据模拟滤波器技术指标设计为响应的模拟低通滤波器；③ 跟据脉冲响应不变法和双线性不变法把模拟滤波器转换为数字滤波器；④ 如果要设计的滤波器是高通、带通或带阻滤波器，则首先把它们的技术指标转化为模拟低通滤波器的技术指标，设计为数字低通滤波器，最后通过频率转换的方法来得到所要的滤波器。

在MATLAB 中，经典法设计IIR 数字滤波器主要采用以下步骤：IIR 数字滤波器设计步骤（二）、用模拟滤波器设计数字滤波器的方法 1、冲激响应不变法：冲激响应不变法是从时域出发，要求数字滤波器的冲激响应h (n ) 对应于模拟滤波器h (t ) 的等间隔抽样。

优点：时域逼近良好；保持线性关系。

缺点：频域响应混叠。

只适用于限带低通滤波器和带通滤波器2、双线性变换法优点：克服了频域混叠模拟滤波器原型 buttap,cheb1ap频率变换 模拟离散化 bilinear,impin varIIR 数字滤波器/Tπ/T π-3/Tπ3/Tπ-j ΩσjIm (z)Re(z)1S 平面Z 平面1S ~S T Tππ-将整个平面压缩变换到平面一个的带状区域缺点：高频时会引起畸变1）冲激响应不变法impinvar格式：[BZ,AZ]= impinvar （B,A,Fs ）功能：把具有[B,A]模拟滤波器传递函数模型转换为采样频率为Fs 的数字滤波器的传递函数模型[BZ,AZ]，Fs 默认值为1。

例：一个4阶的Butterworth 模拟低通滤波器的系统函数如下：12251)(234++++=s s s s s H a试用冲激响应不变法求出Butterworth 模拟低通数字滤波器的系统函数。

郑州轻工业学院课程设计说明书题目:姓名：院< 系）：专业班级:学号：指导教师:成绩:时间：年月日至年月曰郑州轻工业学院课程设计任务书题目基于双线性变换法的IIR数字高通滤波器设计专业、班级电子信息工程08级2班学号姓名主要内容、基本要求、主要参考资料等：主要内容：首先依据给定的性能指标，采用双线性变换法设计IIR数字高通滤波器；然后利用MATLAB软件的wavread函数读取.wav格式的语音信号，并利用所设计的滤波器对音频信号进行滤波处理，画出滤波前后信号的时域波形及频谱；最后回放语音信号，分析滤波前后的语音变化。

基本要求：1滤波器技术指标为：f p=3200Hz。

A p=1dB。

f s=3000Hz。

A s=100dB2、采用双线性变换法设计IIR数字低通滤波器；3、掌握利用wavread函数读取.wav格式语音信号的方法；4、对语音信号进行滤波，并画出滤波前后信号的时域波形及频谱；5、回放语音信号，分析滤波前后的语音变化。

主要参考资料：1、从玉良•数字信号处理原理及其MATLAB实现[M].北京：电子工业出版社.2009.72、胡广书.数字信号处理理论、算法与实现[M].北京：清华大学出版社.2003,8完成期限：2018624— 2018628指导教师签名：课程负责人签名：2018年6月24日摘要随着信息时代和数字世界的到来，数字信号处理已成为当今一门极其重要的学科和技术领域。

在数字信号处理中起着重要的作用并已获得广泛应用的是数字滤波器＜DF, DigitalFilter ）。

数字滤波器是一种用来过滤时间离散信号的数字系统，通过对抽样数据进行数学处理来达到频域滤波的目的。

本次主要利用matlab的信号处理工具箱设计一个数字高通IIR滤波器，并用此滤波器处理一段音频信号。

文中主要介绍了用双线性变换法设计切比雪夫I型高通数字滤波器的实现方法。

关键字：数字信号处理数字滤波器切比雪夫I双线性变换MATLAB1数字滤波器1.1数字滤波器介绍数字滤波器是具有一定传输选择特性的数字信号处理装置，其输入、输出均为数字信号，实质上是一个由有限精度算法实现的线性时不变离散系统。

一、实验目的1．了解工程上两种最常用的变换方法：脉冲响应不变法和双线性变换法。

2．掌握双线性变换法设计IIR 滤波器的原理及具体设计方法，熟悉用双线性设计法设计低通、带通和高通IIR 数字滤波器的计算机程序.3．观察用双线性变换法设计的滤波器的频域特性，并与脉冲响应不变法相比较,了解双线性变换法的特点。

4．熟悉用双线性变换法设计数字Butterworth 和Chebyshev 滤波器的全过程。

5．了解多项式乘积和多项式乘方运算的计算机编程方法。

二、实验原理与方法从模拟滤波器设计IIR 数字滤波器具有四种方法：微分－差分变换法、脉冲响应不变法、双线性变换法、z 平面变换法。

工程上常用的是其中的两种:脉冲响应不变法、双线性变换法。

脉冲响应不变法需要经历如下基本步骤：由已知系统传输函数H(S)计算系统冲激响应h(t)；对h(t)等间隔采样得到h （n ）=h （n T)；由h （n ）获得数字滤波器的系统响应H （Z)。

这种方法非常直观，其算法宗旨是保证所设计的IIR 滤波器的脉冲响应和模拟滤波器的脉冲响应在采样点上完全一致。

而双线性变换法的设计准则是使数字滤波器的频率响应与参考模拟滤波器的频率响应相似。

脉冲响应不变法一个重要的特点是频率坐标的变换是线性的（），其确定是有频谱的周期延拓效应，存在频谱混叠的现象。

为了克服脉冲响应不变法可能产生的频谱混叠，提出了双线性变换法，它依靠双线性变换式：， ， 其中 ,建立其S 平面和Z 平面的单值映射关系，数字域频率和模拟域频率的关系是： , （3－1） 由上面的关系式可知,当时，终止在折叠频率处,整个轴单值的对应于单位圆的一周。

因此双线性变换法不同于脉冲响应不变法，不存在频谱混叠的问题。

从式（3－1）还可以看出，两者的频率不是线性关系。

这种非线性关系使得通带截至频率、过渡带的边缘频率的相对位置都发生了非线性畸变。

这种频率的畸变可以通过预畸变来校正。

用双线性变换法设计数字滤波器时，一般总是先将数字滤波器的个临界频率经过式（3－1）的频率预畸变,求得相应参考模拟滤波器的个临界频率，然后设计参考模拟滤波器的传递函数，最后通过双T Ω=ω1111--+-=z z s s s z -+=11Ω+=j s σωj re z =)2/(ωtg =Ω)(2Ω=arctg ω∞→Ωωπω=Ωj线性变换式求得数字滤波器的传递函数。