课程设计
课程设计名称：数字信号处理课程设计
专业班级：
学生姓名：
学号：
指导教师：
课程设计时间：2014-6-16至2014-6-20
电子信息工程专业课程设计任务书
说明：本表由指导教师填写，由教研室主任审核后下达给选题学生，装订在设计（论文）首页
一、需求分析
用双线性变换法设计原型低通为切比雪夫I 型的数字IIR 高通滤波器，要求通带边界频率为500Hz ，阻带边界频率分别为400Hz ，通带最大衰减1dB ，阻带最小衰减40dB ，抽样频率为2000Hz ，用MATLAB 画出幅频特性，画出并分析滤波器传输函数的零极点；
信号)2s i n ()2s i n ()()()(2121t f t f t x t x t x ππ+=+=经过该滤波器，其中
=1f 300Hz ，=2f 600Hz ，滤波器的输出)(t y 是什么？用Matlab 验证你的结论并给出)(),(),(),(21t y t x t x t x 的图形。
二、设计内容
IIR 数字滤波器的设计通常采用脉冲响应不变法和双线性变换法进行设计。而我选择的课程设计是采用双线性变换法设计切比雪夫I 型的数字IIR 高通滤波器。
三、设计原理
1.双线性变换法
脉冲响应不变法的主要缺点是频谱交叠产生的混淆，这是从S 平面到Z 平面的标准变换z=e 的多值对应关系导致的。为了克服冲激响应法可能产生的频率响应的混叠失真，这是因为从S 平面到Ｚ平面是多值的映射关系所造成的。为了克服这一缺点，可以采用非线性频率压缩方法，将整个频率轴上的频率范围压缩到-π/T ～π/T 之间，再用z =e sT 转换到Z 平面上。也就是说，第一步先将整个S 平面压缩映射到S 1平面的-π/T ～π/T 一条横带里；第二步再通过标准变换关系
z =e s 1T 将此横带变换到整个Z 平面上去。这样就使S 平面与Z 平面建立了一一对应的单值关系，消除了多值变换性，也就消除了频谱混叠现象，映射关系如图1所示：
Z 平面
S 1平面
S 平面
图1 双线性变换的映射关系
稳定性分析：
（1.1） 令s =σ+j ω，则有
（1.2） 为了将S 平面的整个虚轴j Ω压缩到S1平面j Ω1轴上的-π/T 到π/T 段上，
可以通过以下的正切变换实现 ： 式中,T 仍是采样间隔。 优点
与脉冲响应不变法相比，双线性变换的主要优点：靠频率的严重非线性关系得到S 平面与Z 平面的单值一一对应关系，整个j Ω轴单值对应于单位圆一周，这个关系就是式 所表示的，其中ω和Ω为非线性关系。如图图中看到，在零频率附近，Ω～ω接近于线性关系，Ω进一步增加时，ω增长变得缓慢，(ω终止于折叠频率处)，所以双线性变换不会出现由于高频部分超过折叠频率而混淆到低频部分去的现象。 缺点
双线性变换法的缺点：Ω与ω的非线性关系，导致数字滤波器的幅频响应相对于模拟滤波器的幅频响应有畸变，(使数字滤波器与模拟滤波器在响应与频率的对应关系上发生畸变)。例如，一个模拟微分器，它的幅度与频率是线性关系，但通过双线性变换后，就不可能得到数字微分器。
1
1112--+-=
z z T s s T s T z -+=?222
22
2)/2()/2(ωσωσ+-++=T T z ??? ??Ω=Ω2tan 21T T
为有效通带截止频率，表示与通带波纹有关的参量，值越大通带不动切比雪夫滤波器的振幅平方特性如图所示：
图2 切比雪夫滤波器的振幅平方特性
N为偶数，cos2()=1，得到min，
，（1.5）
N为奇数，cos2(，得到max，
（1.6）
有关参数的确定:
a、通带截止频率：
预先给定;
b、与通带波纹有关的参数,通带波纹表示成
（1.7）
所以,
,
给定通带波纹值
分贝数后，可求得
。
c、阶数N:由阻带的边界条件确定。
、A2为事先给定的边界条件,即在阻带中的频率点处
,要求滤波器频响衰
减到1/A2以上。
(1-8)
(1-9)
(1-10)
（1-11）3.数字滤波器
数字滤波器由数字乘法器、加法器和延时单元组成的一种算法或装置。数字滤波器的功能是对输入离散信号的数字代码进行运算处理，以达到改变信号频谱
的目的。
图3 数字滤波器的方框图
数字滤波器是一个离散时间系统（按预定的算法，将输入离散时间信号(对应数字频率)转换为所要求的输出离散时间信号的特定功能装置）。应用数字滤波器处理模拟信号（对应模拟频率）时，首先须对输入模拟信号进行限带、抽样和模数转换。数字滤波器输入信号的数字频率（2π*f/fs,f 为模拟信号的频率，fs 为采样频率，注意区别于模拟频率),按照奈奎斯特抽样定理，要使抽样信号的频谱不产生重叠，应小于折叠频率（ws/2=π)，其频率响应具有以2π为间隔的周期重复特性，且以折叠频率即ω=π点对称。为得到模拟信号，数字滤波器处理的输出数字信号须经数模转换、平滑。数字滤波器具有高精度、高可靠性、可程控改变特性或复用、便于集成等优点。数字滤波器在语言信号处理、图像信号处理、医学生物信号处理以及其他应用领域都得到了广泛应用。
四、设计思路
根据以上IIR 数字滤波器设计方法，下面运用双线性变换法基于MATLAB 设计一个IIR 高通滤波器。
1.在设计高通滤波器之前，首先要根据工程实际的需要确定滤波器的技术指标ωp 和ωs ；
2.将数字高通滤波器的技术指标转换成模拟高通滤波器的技术指标Ωp ，Ωs ，转
换公式为 ；
3.利用频率变换λ=1/η将模拟高通滤波器技术指标转换成归一化模拟低通滤波器G （p ）的技术指标；
??
? ??Ω=Ω2tan 21T T
4.设计模拟低通滤波器G （P ）,并去归一化得：H （s ）=G （p ）
5.采用双线性变换将模拟低通滤波器H （s ）转换成数字低通滤波器H （z ）：
（1.12） 6.采用频带变换，将数字低通滤波器转换成所需类型的数字高通滤波器。
五、程序流程图
↓
↓
↓
↓
↓
↓
1
1
112)
()(--+-==z z T s s H z H
六、程序源代码
clc ； clear all;
Fp=500;%通带截止频率 Fs=400;%组带截止频率 Ap=1;%通带最大衰减 As=40;%
阻带最小衰减 Ft=2000;%抽样频率 Wp=2*pi*Fp/Ft; Ws=2*pi*Fs/Ft; wp=tan(Wp/2); ws=tan(Ws/2); ws1=1 %归一化
wp1=wp/ws %估计滤波器的阶数
[N,Wn]=cheb1ord(ws1,wp1,Ap,As,'s'); %设计滤波器[B,A]=cheby1(N,1,Wn,'s');
[BT,AT]=lp2hp(B,A,wp);
[num,den]=bilinear(BT,AT,0.5)
[z,p,k]=tf2zp(num,den) %显示传输函数
disp('分子系数:');disp(num);
disp('分母系数:');disp(den); %计算增益响应
w = 0:pi/255:pi;
h = freqz(num,den,w);
g = 20*log10(abs(h)); %绘制增益响应
figure;
plot(w/pi,g);grid;%绘制切比雪夫高通滤波器幅频特性axis([0 1 -100 1]);
xlabel('omega/\pi');
ylabel('增益/dB');
title(‘切比雪夫高通滤波器');
figure;
zplane(z,p); %绘制传输函数零极点
title(‘传输函数的零极点');
f1=300;f2=6000;
t=0:0.0001:1
x1=sin(2*pi*f1*t);
x2=sin(2*pi*f2*t);
x=x1+x2;
figure;
subplot(2,2,1)%绘制x1的波形plot(x1);grid on;
axis([0,50*pi,-3,3]);
xlabel('t');ylabel('x1(t)'); title('x1的波形');
subplot(2,2,2) %绘制x1的波形plot(x2);grid on;
axis([0,50*pi,-3,3]);
xlabel('t');ylabel('x2(t)'); title('x2的波形');
subplot(2,2,3)%绘制输入x的波形plot(x);grid on;
axis([0,50*pi,-3,3]);
xlabel('t');ylabel('x(t)');
title(‘输入信号x的波形')
%X=fft(x);
y=filter(num,den,x);%数字滤波器输出
subplot(2,2,4);%绘制输出y的波形
plot(real(y));
grid on;
axis([0,50*pi,-3,3]);
xlabel(‘t');ylabel(‘y');
title(‘滤波器输出y的波形');
七、仿真结果图
八、参考资料
[1] 程佩青.数字信号处理教程[M]，北京：清华大学出版社，2001
[2] Sanjit K. Mitra.孙洪，余翔宇译，数字信号处理实验指导书（MATLAB版）[M]，北京：电子工业出版社，2005年1月
[3] 郭仕剑等.MATLAB 7.x数字信号处理[M]，北京：人民邮电出版社，2006年
[4] 胡广书.数字信号处理理论算法与实现[M]，北京：清华大学出版社，2003年
[5] 杨毅明.数字信号处理[M]，北京：机械工业出版社，2012年
九、设计心得
经过上学期数字信号处理课程的学习，我对DSP有了初步的了解。数字信号处理是研究用数字方法对信号进行分析、变换、滤波、检测、调制、解调以及快速算法的一门技术学科。数字滤波器的实用型式很多，大略可分为有限冲激响应型和无限冲激响应型两类，可用硬件和软件两种方式实现。在硬件实现方式中，它由加法器、乘法器等单元所组成，这与电阻器、电感器和电容器所构成的模拟滤波器完全不同。数字信号处理系统很容易用数字集成电路制成，显示出体积小、稳定性高、可程控等优点。数字滤波器也可以用软件实现。软件实现方法是借助于通用数字计算机按滤波器的设计算法编出程序进行数字滤波计算。
由于上学期的DSP课程是双语教学，加之课程难度也不小，因此，学起来十分的吃力。这次的课程设计让我有一次地接触到了数字信号处理和MATLAB 软件。刚开始的时候还是感觉无从下手，因为之前的基本功不扎实，但是随着老师的讲解和自己私下查找的资料，渐渐掌握到了门路。这次的课设加深了我对课堂抽象概念的理解,巩固了课堂上所学的理论知识,并能很好地理解与掌握数字信号处理中的基本概念、基本原理、基本分析方法。同时掌握编程方法和解决实际问题的技巧。
目录 1. 数字滤波器的设计任务及要求 (2) 2. 数字滤波器的设计及仿真 (2) 2.1数字滤波器的设计 (3) 2.2数字滤波器的性能分析 (3) 3. 数字滤波器的实现结构对其性能影响的分析 (8) 3.1数字滤波器的实现结构一及其幅频响应 (10) 3.2数字滤波器的实现结构二及其幅频响应 (12) 3.3 数字滤波器的实现结构对其性能影响的小结 (12) 4. 数字滤波器的参数字长对其性能影响的分析 (13) 4.1数字滤波器的实现结构一参数字长及幅频响应特性变化 4.2数字滤波器的实现结构二参数字长及幅频响应特性变化 4.3 数字滤波器的参数字长对其性能影响的小结 (16) 5. 结论及体会 (16) 5.1 滤波器设计、分析结论 (16) 5.2 我的体会 (16) 5.3 展望 (16)
1.数字滤波器的设计任务及要求 1. 设计说明 每位同学抽签得到一个四位数，由该四位数索引下表确定待设计数字滤波器的类型及其设计方法，然后用指定的设计方法完成滤波器设计。 要求：滤波器的设计指标： 低通： (1)通带截止频率πrad (id) pc 32 ln = ω (2)过渡带宽度πrad ) (i d 160 10log tz ≤?ω (3)滚降dB αroll 60= 其中，i d — 抽签得到那个四位数（学号的最末四位数），本设计中i d =0201。 2. 滤波器的初始设计通过手工计算完成； 3. 在计算机辅助计算基础上分析滤波器结构对其性能指标的影响（至少选择两种以上合适的滤波器 结构进行分析）； 4. 在计算机辅助计算基础上分析滤波器参数的字长对其性能指标的影响； 5. 以上各项要有理论分析和推导、原程序以及表示计算结果的图表； 6. 课程设计结束时提交设计说明书。 2.数字滤波器的设计及仿真 2.1数字滤波器（编号0201）的设计 数字滤波器是数字信号处理的重要工具之一，它通过数值运算处理改变输入信号所含频率成分的相对比例或者滤出某些频率成分的数字器件或程序，而数字滤波器处理精度高、体积小、稳定、重量轻、灵活、不存在阻抗匹配问题，可以实现模拟滤波器无法实现的特殊功能。 本次课程设计使用MATLAB 信号处理箱和运用切比雪夫法设计数字滤波器，将手工计算一个切比雪夫I 型的IIR 的低通模拟滤波器的系统函数，并在MATLAB 的FDATool 设计工具分析其性能指标。
目录 1.题目.......................................................................................... .2 2.要求 (2) 3.设计原理 (2) 3.1 数字滤波器基本概念 (2) 3.2 数字滤波器工作原理 (2) 3.3 巴特沃斯滤波器设计原理 (2) 3.4脉冲响应不法 (4) 3.5实验所用MA TLAB函数说明 (5) 4.设计思路 (6) 5、实验内容 (6) 5.1实验程序 (6) 5.2实验结果分析 (10) 6.心得体会 (10) 7.参考文献 (10)
一、题目：巴特沃斯数字低通滤波器 二、要求：利用脉冲响应不变法设计巴特沃斯数字低通滤波器，通带截止频率100HZ,采样频率1000HZ ，通带最大衰减为0.5HZ ，阻带最小衰减为10HZ ，画出幅频、相频相应相应曲线。并假设一个信号x(t)=sin(2*pi*f1*t)+sin(2*pi*f2*t),其中f1=50HZ,f2=200HZ 。用此信号验证滤波器设计的正确性。 三、设计原理 1、数字滤波器的基本概念 所谓数字滤波器，是指输入、输出均为数字信号，通过数值运算处理改变输入信号所含频率成分的相对比例，或者滤波器除某些频率成分的数字器件或程序，因此，数字滤波的概念和模拟滤波相同，只是的形式和实现滤波方法不同。正因为数字滤波通过数值运算实现滤波，所以数字滤波处理精度高、稳定、体积小、质量轻、灵活、不存在阻抗匹配问题，可以实验模拟滤波器无法实现的特殊滤波功能。如果要处理的是模拟信号，可通过A\DC 和D\AC,在信号形式上进行匹配转换，同样可以使用数字滤波器对模拟信号进行滤波。 2、数字滤波器的工作原理 数字滤波器是一个离散时间系统，输入x(n)是一个时间序列，输出y(n)也是一个时间序列。如数字滤波器的系统函数为H(Z),其脉冲响应为h(n),则在时间域内存在下列关系 y(n)=x(n) h(n) 在Z 域内，输入输出存在下列关系 Y(Z)=H(Z)X(Z) 式中，X(Z),Y(Z)分别为输入x(n)和输出y(n)的Z 变换。 同样在频率域内，输入和输出存在下列关系 Y(jw)=X(jw)H(jw) 式中，H(jw)为数字滤波器的频率特性，X(jw)和Y(jw)分别为x(n)和y(n)的频谱。w 为数字角频率，单位rad 。通常设计H(jw)在某些频段的响应值为1，在某些频段的响应为0.X(jw)和H(jw)的乘积在频率响应为1的那些频段的值仍为X(jw)，即在这些频段的振幅可以无阻碍地通过滤波器，这些频带为通带。X(jw)和H(jw)的乘积在频段响应为0的那些频段的值不管X(jw)大小如何均为零，即在这些频段里的振幅不能通过滤波器，这些频带称为阻带。 一个合适的数字滤波器系统函数H(Z)可以根据需要输入x(n)的频率特性，经数字滤波器处理后的信号y(n)保留信号x(n)中的有用频率成分，去除无用频率成分。 3、巴特沃斯滤波器设计原理 （1）基本性质 巴特沃斯滤波器以巴特沃斯函数来近似滤波器的系统函数。巴特沃斯滤波器是根据幅频特性在通频带内具有最平坦特性定义的滤波器。 巴特沃思滤波器的低通模平方函数表示1 () ΩΩ+ =Ωc N /22 a 11 ) (j H
题目滤波器设计与实现 学生姓名王朝全学号201005020228 专业班级电信1002 设计内容与要求一、设计内容： 设计一个模拟高通切比雪夫Ⅰ滤波器，技术指标：通带截止频率2000Hz，通带最大衰减1dB；过渡带700Hz，阻带最小衰减60dB，画出其幅度谱和相位谱。 二、设计要求： 1 设计报告一律按照规定的格式，使用A4纸，格式、封面统一给出模版。 2 报告内容 （1）设计题目及要求 （2）设计原理(包括滤波器工作原理、涉及到的matlab函数的说明) （3）设计内容（设计思路，设计流程、仿真结果） （4）设计总结（收获和体会） （5）参考文献 （6）程序清单 起止时间2012年12 月3日至2012年12月11 日指导教师签名2012年12月2日 系（教研室）主任签 名 年月日 学生签名年月日
一、设计题目及要求 设计一个模拟高通切比雪夫Ⅰ滤波器，技术指标：通带截止频率 2000Hz ，通带最大衰减1dB ；过渡带700Hz ，阻带最小衰减60dB ，画出其幅度谱和相位谱。 二、设计原理 (包括滤波器工作原理、涉及到的matlab 函数的说明) 2.1切比雪夫滤波器的内容： 2.1.1、切比雪夫滤波器的幅度平方函数： 从公式中我们可以发现它需要计算的 ε和N ， 110 1.0-=p a ε 边带频率p Ω，通带最大衰减p α=0.1dB 、阻带最大衰减s α=50dB 、 阻带截止频率s Ω，它们满足： ) /cosh() /110cosh(11.0p s a ar ar N s ΩΩ-≥ -ε 根据公式可以求出最小的切比雪夫I 型滤波器阶数的N 最小值。 2.1.2、切比雪夫滤波器的三个参量： Ωc 通带截止频率，ε表征通带内波纹大小由： 2 222 1 ()() 1( )a N p A H j C εΩ=Ω= Ω +Ω
******************* 实践教学 ******************* 兰州理工大学 计算机与通信学院 2013年春季学期 信号处理课程设计 题目：切比雪夫I型低通滤波器设计 专业班级：通信工程三班 姓名： 学号： 指导教师：蔺莹 成绩：
摘要 本次课程设计将完成一个数字切比雪夫低通IIR 滤波器的设计，利用双线性变换和冲激响应不变法完成设计，并利用MATLAB 进行仿真。 已知数字滤波器的性能指标为：通带截止频率为: 0.4,1,0.45,15p p s P R dB R dB ω=π=ω=π=通带波动为阻带波动为，要求设计满足以上技术指标的切比雪夫I 型低通滤波器。绘制出理想冲激响应和实际冲激响应结果图。并且给出幅度响应结果图。 关键字：数字滤波器 切比雪夫 双线性变换 冲激响应不变
前言 (1) 一．数字滤波器 (2) 1．1 数字滤波器的概念 (2) 1．2数字滤波器的分类 (2) 1．3 IIR数字滤波器设计原理 (3) 二．切比雪夫滤波器 (5) 三．双线性变换法 (8) 四．脉冲响应不变法 (12) 五．切比雪夫低通滤波器的设计 (15) 5．1 程序流程图 (15) 5．2 设计步骤 (15) 六．总结 (18) 七．参考文献 (19) 致谢 (20) 附录 (21)
随着信息时代和数字世界的到来，数字信号处理已成为当今一门极其重要的学科和技术领域。目前数字信号处理在通信、语音、图像、自动控制、雷达、军事、航空航天、医疗和家用电器等众多领域得到了广泛的应用。在数字信号处理中起着重要的作用并已获得广泛应用的是数字滤波器（DF，Digital Filter），根据其单位冲激响应函数的时域特性可分为两类：无限冲激响应IIR（Infinite Impulse Response）滤波器和有限冲激响应FIR（Finite Impulse Response）滤波器。与FIR 滤波器相比，IIR的实现采用的是递归结构，极点须在单位圆内，在相同设计指标下，实现IIR滤波器的阶次较低，即所用的存储单元少，从而经济效率高。MATLAB是英文MATrix LABoratory(矩阵实验室)的缩写。它是美国的MathWorks公司推出的一套用于科学计算和图形处理可视化、高性能语言与软件环境。MATLAB的信号处理工具箱是专门应用于信号处理领域的专用工具箱，它的两个基本组成就是滤波器的设计与实现部分以及谱分析部分。工具箱提供了丰富而简便的设计，使原来繁琐的程序设计简化成函数的调用。只要以正确的指标参数调用相应的滤波器设计程序或工具箱函数，便可以得到正确的设计结果，使用非常方便。
一、设计一个切比雪夫Ⅱ型低通滤波器 wp=0.2*pi; %通带边界频率； ws=0.4*pi; %阻带截止频率； rp=1; %通带最大衰减； rs=80; %阻带最小衰减； Fs=1000 %假设抽样脉冲1000hz [N,Wn]=cheb2ord(wp,ws,rp,rs,'s'); %Chebyshev II型滤波器参数计算（模拟域）； [Z,P,K]=cheby2(N,rs,Wn,'s'); %构造Chebyshev II型滤波器（零极点模型）；[H,W]=zp2tf(Z,P,K); %将零极点模型转化成传递函数的模型； figure(1); freqs(H,W); %在Figure1上显示滤波器的幅频响应及相频响应；[P,Q]=freqs(H,W); %返回滤波器的冲击响应的复数形式； figure(2); plot(Q*Fs/(2*pi),abs(P));grid; %在Figure2上显示幅频特性曲线； xlabel('频率/Hz'); ylabel('幅值');
二、设计一个高通Chebyshow型数字滤波器 wp=100;ws=80;Fs=300;rp=1;rs=45; %数字滤波器的各项指标； WP=100*2*pi; %把数字滤波器的频率特征转换成模拟滤波器的频率特征； WS=300*2*pi; [N,Wn]=cheb2ord(WP,WS,rp,rs,'s'); %Chebyshev II型滤波器参数计算（模拟域）； [Z,P,K]=cheb2ap(N,rs); %创建Chebyshev滤波器原型； [A,B,C,D]=zp2ss(Z,P,K); %表达式从零极点增益形式转换成状态方程形式； [AA,BB,CC,DD]=lp2hp(A,B,C,D,Wn); %实现低通到高通滤波器类型的转换；[a,b,c,d]=bilinear(AA,BB,CC,DD,Fs); %采用双线性变换法，从模拟高通到数字高通； [P,Q]=ss2tf(a,b,c,d); %表达式从状态方程形形式转换成传输函数形式；figure(1); freqz(P,Q); %绘出频率响应； [H,W]=freqz(P,Q); figure(2); plot(W*Fs/(2*pi),abs(H));grid; xlabel('频率/Hz'); ylabel('幅值');
设计题目基于切比雪夫I型的数字高通滤波器的设计 设计要求 设计一个9阶切比雪夫I型高通滤波器,通带纹波为10dB,下边界频率为400 / rad s ,并绘出其幅频响应曲线 设计过程1．系统设计方案 1.1 Matlab的简介和主要功能： 简介：MATLAB 是一种用于算法开发、数据可视化、数据分析以及数值计算的高级技术计算语言和交互式环境。使用 MATLAB，您可以较使用传统的编程语言（如 C、C++ 和 Fortran）更快地解决技术计算问题。 MATLAB 的应用范围非常广，包括信号和图像处理、通讯、控制系统设计、测试和测量、财务建模和分析以及计算生物学等众多应用领域。附加的工具箱（单独提供的专用 MATLAB 函数集）扩展了 MATLAB 环境，以解决这些应用领域内特定类型的问题。 MATLAB 提供了很多用于记录和分享工作成果的功能。可以将您的 MATLAB 代码与其他语言和应用程序集成，来分发您的 MATLAB 算法和应用。 主要功能:1.此高级语言可用于技术计算 2.此开发环境可对代码、文件和数据进行管理 3.交互式工具可以按迭代的方式探查、设计及求解问题 4.数学函数可用于线性代数、统计、傅立叶分析、筛选、优化以及数值 积分等 5.二维和三维图形函数可用于可视化数据 6.各种工具可用于构建自定义的图形用户界面 7.各种函数可将基于 MATLAB 的算法与外部应用程序和语言（如 C、 C++、Fortran、Java、COM 以及 Microsoft Excel）集成 1.2 开发算法和应用程序 开发算法和应用程序 MATLAB 提供了一种高级语言和开发工具，使您可以迅速地开发并分析算法和应用程序。
DSP 设计滤波器报告 姓名：张胜男 班级：07级电信（1）班 学号：078319120 一·低通滤波器的设计 （一）实验目的：掌握IIR 数字低通滤波器的设计方法。 （二）实验原理： 1、滤波器的分类 滤波器分两大类：经典滤波器和现代滤波器。 经典滤波器是假定输入信号)(n x 中的有用成分和希望取出的成分各自占有不同的频带。这样，当)(n x 通过一个线性系统（即滤波器）后可讲欲去除的成分有效的去除。 现代滤波器理论研究的主要内容是从含有噪声的数据记录（又称时间序列）中估计出信号的某些特征或信号本身。 经典滤波器分为低通、高通、带通、带阻滤波器。每一种又有模拟滤波器（AF ）和数字滤波器（DF ）。对数字滤波器，又有IIR 滤波器和FIR 滤波器。 IIR DF 的转移函数是： ∑∑=-=-+==N k k k M r r r z a z b z X z Y z H 10 1)()()( FIR DF 的转移函数是： ∑-=-=10)()(N n n z n h z H FIR 滤波器可以对给定的频率特性直接进行设计，而IIR 滤波器目前最通用的方法是利用已经很成熟的模拟滤波器的设计方法进行设计。 2、滤波器的技术要求 低通滤波器： p ω：通带截止频率（又称通带上限频率） s ω：阻带下限截止频率 p α：通带允许的最大衰减 s α：阻带允许的最小衰减 （p α，s α的单位dB ） p Ω：通带上限角频率 s Ω：阻带下限角频率 （s p p T ω=Ω，s s s T ω=Ω）即 C p p F ωπ2=Ω C s s F ωπ2=Ω 3、IIR 数字滤波器的设计步骤：
华北科技学院课程设计任务书 2013 — 2014 学年第二学期 电子信息工程学院（系、部）通信工程专业 B111 班级课程名称：移动通信 设计题目：巴特沃斯、切比雪夫滤波器的仿真与实现完成期限：自16 周至 18 周共 3 周
目录 1．前言 (3) 1.1 MATLAB (3) 1.2 滤波器的概念 (5) 1.2.1滤波器的原理 (6) 1.2.2理想滤波器与实际滤波器 (6) 1.2.3 滤波器的分类 (7) 2.设计目的 (9) 3.设计原理 (9) 3.1.模拟低通滤波器的设计指标及逼近方法 (9) 3.2.巴特沃斯低通滤波器的设计方法 (10) 3.3.切比雪夫滤波器的设计方法 (14) 4.详细设计与系统分析 (21) 4.1程序设计 (21) 4.1.1巴特沃斯滤波器 (21) 4.1.2切比雪肤滤波器 (23) 4.2同一滤波器不同参数的比较 (25) 4.2.1巴特沃斯滤波器 (25) 4.2.2切比雪夫滤波器 (27) 4.3不同滤波器同一阶数的比较 (30) 4.3.1低通滤波器 (30) 4.3.2高通滤波器 (30) 4.3.3带通滤波器 (31) 4.3.4带阻滤波器 (31) 5.心得体会 (32) 6.参考文献 (32)
摘要：利用MATLAB设计滤波器,可以按照设计要求非常方便地调整设计参数,极大地减轻了设计的工作量,有利于滤波器设计的最优化。MATLAB因其强大的数据处理功能被广泛应用于工程计算,其丰富的工具箱为工程计算提供了便利,利用MATLAB信号处理工具箱可以快速有效地设计各种数字滤波器,设计简单方便。本文介绍了在MATLAB R2009a 环境下滤波器设计的方法和步骤。关键词：滤波器，MATLAB 1．前言 1.1 MATLAB MATLAB是美国MathWorks公司开发的一种功能极其强大的高技术计算语言和内容极其丰富的软件库，集数值计算、矩阵运算和信号处理与显示于一身。该软件最初是由美国教授Cleve Moler 创立的。1980年前后，他在教线性代数课程时，发现用其他高级语言编程时极不方便，便构思开发了MATLAB，即矩阵实验室（Matrix Laboratory）。该软件利用了当时代表数值线性代数领域最高水平的EISPACK和LINPACK两大软件包，并且利用Fortran 语言编写了最初的一套交互式软件系统，MATLAB的最初版本便由此产生了。 最初的MATLAB由于语言单一，只能进行矩阵的运算，绘图也只能用原始的描点法，内部函数只有几十个，因此功能十分简单。1984年该公司推出了第一个MATLAB的商业版，并用C语言作出了全部改写。现在的MATLAB程序是MathWorks公司用C语言开发的，第一版由steve Bangert主持开发编译解释程序，Steve Kleiman完
摘要 随着信息与数字时代的到来,数字信号处理已成为当今一门极其重要的学科与技术领域。在现代通信系统中,由于信号中经常混有各种复杂成分,因此很多信号的处理都就是基于滤波器而进行的。所以,数字滤波器在数字信号处理中起着举足轻重的作用。而数字滤波器的设计都要以模拟滤波器为基础的,这就是因为模拟滤波器的理论与设计方方法都已发展的相当成熟, 且有典型的模拟滤波器供我们选择。,如巴特沃思滤波器、切比雪夫滤波器等。 本次课程设计将运用MATLAB设计一个基于切比雪夫低通滤波器,并出所设计滤波器的幅度及幅度衰减特性。 关键词:模拟低通滤波切比雪夫
1课题描述 数字滤波器就是数字信号处理的重要工具之一,它通过数值运算处理改变输入信号所含频率成分的相对比例或者滤出某些频率成分的数字器件或程序,而数字滤波器处理精度高、体积小、稳定、重量轻、灵活、不存在阻抗匹配问题,可以实现模拟滤波器无法实现的特殊功能。故本课题使用MATLAB 信号处理箱与运用切比雪夫法设计数字低通滤波器。 2设计原理 2、1切比雪夫滤波器介绍 在巴特沃兹滤波器中,幅度响应在通带与阻带内都就是单调的。因此,若滤波器的技术要求就是用最大通带与阻带的逼近误差来给出的话,那么,在靠近通带低频端与阻带截止频率以上的部分都会超出技术指标。一种比较有效的途径就是使逼近误差均匀地分布于通带或阻带内,或同时在通带与阻带内都均匀分布,这样往往可以降低所要求的滤波器阶次。通过选择一种具有等波纹特性而不就是单调特性的逼近方法可以实现这一点。切比雪夫型滤波器就具有这种性质:其频率响应的幅度既可以在通带中就是等波纹的,而在阻带中就是单调的(称为I 型切比雪夫滤波器),也可以在通带中就是单调的,而在阻带中就是等波纹的(称为II 型切比雪夫滤波器)。I 型切比雪夫滤波器的幅度平方函数就是 2|)(|Ωj H C =)/(11 22c N C ΩΩ+ε (2、1) 式中为N 阶切比雪夫多项式,定义为 )cos cos()(1x N x C N -= (2、2) 从定义切比雪夫多项式可以直接得出由 )(x C N 与)(1x C N -求)(1x C N +的递推公式。将三角恒等式代入 (2、2)式,得 )(1x C N +=2x -)(x C N ) (1x C N - (2、3) 从 (2、2)式我们注意到,当0
实验2_FIR滤波器设计_切比雪夫一致逼近法实时数字信号处理 FIR数字滤波器设计 姓名: 专业: 学号: 日期: 目录 第1章 FIR数字滤波器设计原 理 ...................................................... 3 1.1 FIR数字滤波器设计方 法 ..................................................................... .............. 3 1.2 FIR数字滤波器的切比雪夫逼近法设计原 理 (3) 第2章 FIR数字滤波器设计 (4) 2.1 语音文件滤波处理的需求分 析 ..................................................................... . (4) 2.2 FIR陷波滤波器参数设 计 ..................................................................... . (4)
第3章 FIR数字滤波器滤波实验与结果分 析 (5) 3.1 MATLAB数字滤波实验...................................................................... (5) 3.1.1 滤波实验流程设 计 ..................................................................... .. (5) 3.1.2 语音文件噪声添 加 ..................................................................... .. (5) 3.1.3 语音文件滤波处 理 ..................................................................... ....................... 6 3.2 语音文件滤波实验结 果 ..................................................................... . (6) 3.2.1 滤波实验结 果 ..................................................................... (6) 3.2.2 滤波实验结果分 析 ..................................................................... .. (8) 第1章 FIR数字滤波器设计原理 1.1 FIR数字滤波器设计方法
目录1.题目...................................................................... (2) 2.要求...................................................................... . (2) 3.设计原理...................................................................... .. (2) 3.1数字滤波器基本概念 (2) 3.2数字滤波器工作原理 (2) 3.3巴特沃斯滤波器设计原理 (2) 3.4脉冲响应不法...................................................................... . (4) 3.5实验所用MATLAB函数说明 (5)
4.设计思路...................................................................... (6) 5、实验内容...................................................................... .. (6) 5.1实验程序...................................................................... (6) 5.2实验结果分析...................................................................... (10) 6.心得体会...................................................................... .. (10) 7.参考文献...................................................................... .. (10) 一、题目：巴特沃斯数字低通滤波器 二、要求：利用脉冲响应不变法设计巴特沃斯数字低通滤波器，通带截止频率100HZ,采样频率1000HZ，通带最大衰减为0.5HZ，阻带最小衰减为10HZ，画出幅频、相频相应相应曲线。并假设一个信号x(t)=sin(2*pi*f1*t)+sin(2*pi*f2*t),其中f1=50HZ,f2=200HZ。用此信号验证滤波器设计的正确性。 三、设计原理 1、数字滤波器的基本概念 所谓数字滤波器，是指输入、输出均为数字信号，通过数值运算处理改变输入信号所含频率成分的相对比例，或者滤波器除某些频率成分的数字器件或程序，因此，数字滤波的概念和模拟滤波相同，只是的形式和实现滤波方法不同。正因为数字滤波通过数值运算实现滤波，所以数字滤波处理精度高、稳定、体积小、质量轻、灵活、不存在阻抗匹配问题，可以实验模拟滤波器无法实现的特殊滤
切比雪夫1型数字滤波器的设计及滤波过程 切比雪夫1型低通模拟滤波器的幅度平方函数为： ) ( 11 )(22 2 2 | )(|Ω ΩΩ + == Ωp N C j H A a ε 其中ε表示通带内幅度波动的程度，ε越大，波动幅度也越大。 110 1.0-=Ap ε ) (x C N 称为N 阶切比雪夫多项式。 1、滤波器设计及结果如下 IIR-DF 滤波器设计（切比雪夫1型） （1） 切比雪夫1型低通数字滤波器的损耗函数曲线、滤波分离出的DSB 信号的时 域波形及其频谱如下： 0.1 0.2 0.3 0.4 0.50.6 0.7 0.8 0.9 1 -100 -90-80-70-60-50-40-30-20-10 0w/π 幅度/d B 切比雪夫1型模拟低通滤波器的幅频响应曲线
00.010.020.03 0.040.050.060.070.08 -1 -0.500.5 1t/s y 1(t ) y1(t)的时域波形 f/Hz 幅度 y1(t)的频谱 其中阶数N=7 （2） 切比雪夫1型带通数字滤波器的损耗函数曲线、滤波分离出的DSB 信号的时 域波形及其频谱如下：
0.1 0.2 0.3 0.4 0.50.6 0.7 0.8 0.9 1 -100 -90-80-70-60-50-40-30-20-10 0w/ 幅度/d B 切比雪夫1型带通数字滤波器幅频响应曲线 0.01 0.02 0.03 0.040.050.06 0.07 0.08 -1-0.500.5 1t/s y 2(t ) y2(t)的时域波形 200 400 600 800 10001200 1400 1600 1800 2000 f/Hz 幅度 y2(t)的频谱 其中阶数N=8
课程设计 课程名称：数字信号处理 题目编号： 0202 题目名称：切比雪夫Ⅱ型IIR低通滤波器 专业名称：电子信息工程 班级：电子1204班 学号： 20124470411 学生姓名：刘春阳 任课教师：黄国玉 2015年09月30日
课程设计任务书
目录 1. 数字滤波器的设计任务及要求（编号202） (2) 2. 数字滤波器的设计及仿真 (3) 2.1数字滤波器（编号202）的设计 (3) 2.2数字滤波器（编号202）的性能分析 (6) 3. 数字滤波器的实现结构对其性能影响的分析 (7) 3.1数字滤波器的实现结构一（直接型）及其幅频响应 (8) 3.2数字滤波器的实现结构二（级联型）及其幅频响应 (10) 3.3 数字滤波器的实现结构对其性能影响的小结 (10) 4. 数字滤波器的参数字长对其性能影响的分析 (11) 4.1数字滤波器的实现结构一（直接型）参数字长及幅频响应特性变化 (12) 4.2数字滤波器的实现结构二（级联型）参数字长及幅频响应特性变化 (14) 4.3 数字滤波器的参数字长对其性能影响的小结 (16) 5. 结论及体会 (16) 5.1 滤波器设计、分析结论 (16) 5.2 我的体会 (16) 5.3 展望 (17)
1.数字滤波器的设计任务及要求（0202） 每位同学抽签得到一个四位数，由该四位数索引下表一确定待设计数字滤波器的类型及其设计方法， 然后用指定的设计方法完成滤波器设计。 要求： （1）滤波器设计指标：通带截止频 pc ln ()32 d rad i πω= , 过渡带宽度 10 tz () 160 log d rad i πω?≤ ，滚降roll 60dB α=; 其中， id —抽签得到那个四位数（题目编号) （2）滤波器的初始设计通过手工计算完成； （3）在计算机辅助计算基础上分析滤波器结构对其性能指标的影响（至少选择两种 以上合适的滤波器结构进行分析）； （4）在计算机辅助计算基础上分析滤波器参数的字长对其性能指标的影响； （5） 以上各项要有理论分析和推导、原程序以及表示计算结果的图表； （6）课程设计结束时提交设计说明书。 2.数字滤波器的设计及仿真 2.1 数字滤波器（编号202）的设计 随着信息和数字时代的到来，数字信号处理已成为当今一门极其重要的学科和技术领域。在现代通信系统中，由于信号中经常混有各种复杂的成分，因此很多信号的处理都是基于滤波器而进行的。所以数字滤波器在数字信号处理中起着举足轻重的作用。而数字滤波器的设计都要以模拟滤波器为基础的，这是因为模拟滤波器的设计方法都已发展的相当成熟，且有典型的模拟滤波器供我们选择。如巴特沃斯滤波器、切比雪夫滤波器等。本次课程设计将手工计算一个切比雪夫II 型的IIR 的低通模拟滤波器的系统函数，并在MATLAB 的FDATool 设计工具分析其性能指标。
切比雪夫滤波器设计和仿真 摘 要：滤波器是一种常见的电路形式，在电子线路中有广泛的应用。滤波器的设计在这些领域中是必不可缺的。滤波器是一种能使有用频率信号通过而同时抑制（或衰减）无用频率信号的电子电路或装置。按照不同的频域或时域特性要求，可分巴特沃斯（Butterworth ）、切比雪夫型（Chebyshev ）、贝塞尔型 （Bessel)椭圆型（Elliptic ），这些都是属于模拟低通滤波器。切比雪夫型滤 波器的特点是通带内增益有起伏。本文介绍的是借助Multisim 10仿真实现二阶切比雪夫低通有源滤波器的设计。 关键词：滤波器；频域或时域特性；切比雪夫；Multisim10 Chebyshev LPF's design and emulation Abstrac ：Filter would completely eliminate signals above the cutoff frequency, and perfectly pass signals below the cutoff frequency . In real filters, various rade-offs are made to get optimum performance for a given application. There are manybooks that provide information on popular filter types like the Butterworth, Bessel, and Chebyshev filters, just to name few. Chebyshev filters are designed to have ripple in the pass-band, but steeper roll off after the cutoff frequency. Cutoff frequency is defined as the frequency at which the response falls below the ripple band. For a given filter order, a steeper cutoff can be achieved by allowing more pass-band ripple. The transient response of a Chebyshev filter to a pulse input shows more overshoot and ringing than a Butterworth filter. Key Words:Filter ；popular filter types ；Chebyshev function ；Multisim 10 引言 随着现代科学技术的发展，滤波技术在通信、测试、信号处理、数据采集和实时控制等领域都得到了广泛的应用。滤波器的设计在这些领域中是必不可缺的，有时甚至是至关重要的环节。比如说，在通信领域，常常利用种滤波器来抑制噪声，去除干扰，以提高信噪比。随着电子计算机的普及和材料科学的进步，各种各样的滤波器的辅助设计软件也得以不断推出，设计人员可以选择高功效的滤波 、管路敷设技术通过管线不仅可以解决吊顶层配置不规范高中资料试卷问题，而且可保障各类管路习题到位。在管路敷设过程中，要加强看护关于管路高中资料试卷连接管口处理高中资料试卷弯扁度固定盒位置保护层防腐跨接地线弯曲半径标高等，要求技术交底。管线敷设技术包含线槽、管架等多项方式，为解决高中语文电气课件中管壁薄、接口不严等问题，合理利用管线敷设技术。线缆敷设原则：在分线盒处，当不同电压回路交叉时，应采用金属隔板进行隔开处理；同一线槽内，强电回路须同时切断习题电源，线缆敷设完毕，要进行检查和检测处理。、电气课件中调试对全部高中资料试卷电气设备，在安装过程中以及安装结束后进行 高中资料试卷调整试验；通电检查所有设备高中资料试卷相互作用与相互关系，根据生产工艺高中资料试卷要求，对电气设备进行空载与带负荷下高中资料试卷调控试验；对设备进行调整使其在正常工况下与过度工作下都可以正常工作；对于继电保护进行整核对定值，审核与校对图纸，编写复杂设备与装置高中资料试卷调试方案，编写重要设备高中资料试卷试验方案以及系统启动方案；对整套启动过程中高中资料试卷电气设备进行调试工作并且进行过关运行高中资料试卷技术指导。对于调试过程中高中资料试卷技术问题，作为调试人员，需要在事前掌握图纸资料、设备制造厂家出具高中资料试卷试验报告与相关技术资料，并且了解现场设备高中资料试卷布置情况与有关高中资料试卷电气系统接线等情况，然后根据规范与规程规定，制定设备调试高中资料试卷方案。 、电气设备调试高中资料试卷技术电力保护装置调试技术，电力保护高中资料试卷配置技术是指机组在进行继电保护高中资料试卷总体配置时，需要在最大限度内来确保机组高中资料试卷安全，并且尽可能地缩小故障高中资料试卷破坏范围，或者对某些异常高中资料试卷工况进行自动处理，尤其要避免错误高中资料试卷保护装置动作，并且拒绝动作，来避免不必要高中资料试卷突然停机。因此，电力高中资料试卷保护装置调试技术，要求电力保护装置做到准确灵活。对于差动保护装置高中资料试卷调试技术是指发电机一变压器组在发生内部故障时，需要进行外部电源高中资料试卷切除从而采用高中资料试卷主要保护装置。
设计一个高通滤波器，并检验它的性能 采样率为10kHZ 阻带边缘为1.5Khz,衰减为40bB 通带边缘为2kHz，波纹为3Db 采用切比雪夫2型滤波器 Fs=1e4 ; fs=1.5e3; fp=2e3; As=40; Rp=3 ;wp=2*fp/Fs; ws=2*fs/Fs; 用MA TLAB工具： 设计的滤波器性能 >>Fs=1e4; >>fs=1.5e3; >>fp=2e3; >>As=40; >>Rp=3; >>wp=2*fp/Fs; >>ws=2*fs/Fs; >>[N,wn]=cheb2ord(wp,ws,Rp,As); >>[b,a]=cheby2(N,As,wn,'high'); >>[db,mag,pha,grd,w]=freqz_m(b,a); >> subplot(2,2,1);plot(w/pi,mag); >> axis([0,1,0,1]); >> setX([0 0.3 0.4 1]); >>setY([0.01 0.7279 1]) >> title('Magnitude Response'); >> subplot(2,2,2);plot(w/pi,db); >> axis([0 1 -70 0]) >> setX([0 0.3 0.4 1]) >> setY([-40 -2.7589]) >> title('Magnitude Response in dB'); 输入信号x=cos(0.6*pi*n); 取200个抽样值，为了便于观看，在画图是裁剪掉100个值
不直观，直接做频谱图 为了说明该滤波器的性能，做一组对比。换上另一种频率的信号x=cos(0.2*pi*n); 由上图可以看到该信号完全被滤去了
目录 1摘要 (3) 2设计原理 (4) 2.1 切比雪夫滤波器介绍 (4) 2.2滤波器的分类 (5) 2.3 模拟滤波器的设计指标 (6) 3切比雪夫I型滤波器 (7) 3.1 切比雪夫滤波器的设计原理 (7) 3.2切比雪夫滤波器的设计步骤 (10) 3.3 用matlab编程设计切比雪夫低通滤波器 (11) 4 总结 (18) 5 参考文献 (18)
摘要 随着信息和数字时代的到来，数字信号处理已成为当今一门极其重要的学科和技术领域。在现代通信系统中，由于信号中经常混有各种复杂成分，因此很多信号的处理都是基于滤波器而进行的。所以，数字滤波器在数字信号处理中起着举足轻重的作用。 而数字滤波器的设计都要以模拟滤波器为基础的，这是因为模拟滤波器的理论和设计方方法都已发展的相当成熟，且有典型的模拟滤波器供我们选择。，如巴特沃思滤波器、切比雪夫滤波器等。 本次课程设计将运用MATLAB设计一个基于切比雪夫低通滤波器，并出所设计滤波器的幅度及幅度衰减特性。 关键词：滤波切比雪夫模拟低通 1切比雪夫滤波器介绍
在巴特沃兹滤波器中，幅度响应在通带和阻带内都是单调的。因此，若滤波器的技术要求是用最大通带和阻带的逼近误差来给出的话，那么，在靠近通带低频端和阻带截止频率以上的部分都会超出技术指标。一种比较有效的途径是使逼近误差均匀地分布于通带或阻带内，或同时在通带和阻带内都均匀分布，这样往往可以降低所要求的滤波器阶次。通过选择一种具有等波纹特性而不是单调特性的逼近方法可以实现这一点。切比雪夫型滤波器就具有这种性质：其频率响应的幅度既可以在通带中是等波纹的，而在阻带中是单调的（称为I 型切比雪夫滤波器），也可以在通带中是单调的，而在阻带中是等波纹的（称为II 型切比雪夫滤波器），其中切比雪夫II 型滤波器又称为逆切比雪夫滤波器。 I 型切比雪夫滤波器的幅度平方函数是 2|)(|Ωj H C = )/(11 2 2c N C ΩΩ+ε (2.1) 其中ε是一个小于1的正数，它与通带波纹有关，ε越大，波纹也越大，式中为N 阶切比雪夫多项式，定义为 )cos cos()(1x N x C N -= (2.2) 当N 大于或等于1时，从定义切比雪夫多项式可以直接得出由 ) (x C N 和 )(1x C N -求)(1x C N +的递推公式。将三角恒等式代入 (2.2)式，得 )(1x C N +=2x -)(x C N )(1x C N - (2.3) 从 (2.2)式我们注意到，当0
IIR数字低通滤波器 一、设计目的 课程设计是理论教学的延伸，是掌握所学知识的一种重要手段，对于贯彻理论联系实际、提高教育质量，培养合格人才等具有特殊作用本次课程设计一方面通过MA TLAB仿真设计内容，使学生加深对理论知识理解的同时增强其逻辑斯维的能力，另一方面对课堂教学中的理论知识做一个总结和补充。 二、设计要求 2.1 简述设计目的 2.2 阐述设计原理 2.3 按步骤设计滤波器，给出系统函数 2.4 用MATLAB语言编程、绘制幅频特性曲线 2.5 定性分析两种设计方法的滤波器的性能。比较优缺点，并判定设计是否能满足要求 三、设计原理 3.1 巴特沃斯滤波器原理 由于已知指标，故可求出滤波器的阶数N，由式知，求出归一化极点，将代入，得到归一化传输函数。也可以根据N查表得到归一化传输函数。然后再将去归一化。将代入，得到实际的滤波器传输函数Ha(S)。这里3dB截止频率可以按照或。这样即可设计出低通巴特沃斯滤波器。巴特沃斯滤波器的幅度响应在通带内具有最平坦的特性，且在通带和阻带内幅度的特性，是单调变化的。模拟巴特沃斯滤波器的幅度平方函数为=2N，式中N称为滤波器的阶数，为角频率，在处幅度响应的平方为。 3.2 双线性变换法工作原理
使数字滤波器的频率响应与模拟滤波器的频率响应相似。冲激响应不变不得法、阶跃响应不变法：时域模仿逼近缺点是产生频率响应的混叠失真双线性变换法也是一种由S平面到z平面的映射过程，双线性变换法与脉冲响应不变法不同，它是一种从S 平面到z平面简单映射。双线性变换中数字域与频率和模拟频率之间的非线性关系限制了它的应用范围，只有当非线性失真是允许的或能被裣时，才能采用双线性变换法，通常低通、高通、带通和带阻等滤波器等具有分段恒定的频率特性，可以采用预畸变的方法来补偿频率畸变，因此可以采用双线性变换设计方法。 3.3 脉冲响应不变法工作原理 冲激响应不变法遵循的准则是使数字滤波器的单位取样响应与参照的模拟滤波器的脉冲响应的取样值完全一样，即h(n)=ha(nT),其中T为取样周期。实际是由模拟滤波器转换成为数字滤波器，就是要建立模拟系统函数Ha（S）与数字系统函数H（z）之间的关系。脉冲响应不变法是从S平面映射到z平面，这种映射不是简单的代数映射，而是S平面的每一条宽为的横带重复地映射到整个z平面。 四、按步骤设计滤波器 4.1用脉冲响应不变法设计低通滤波器 4.1.1数字低通的技术指标为 4.1.2模拟低通的技术指标为 4.1.3设计巴特沃斯低通滤波器。先计算阶数N及3dB截止频率。 取N=9.。将和代入公式，得到3dB截止频率，此值满足通带技术要求，同时给阻带衰减留一定余量，这对防止频率混叠有一定好处。根据阶数N=9,查表得到归一化传输函数为 为去归一化，将代入中，得到实际的传输函数为 4.1.4 用脉冲响应不变法将转换成如下：