实验四iir数字滤波器的设计实验报告
iir滤波器实验报告
iir滤波器实验报告IIR滤波器实验报告引言：数字信号处理是现代通信、音频处理和图像处理等领域中不可或缺的技术。
滤波器作为数字信号处理的重要组成部分，被广泛应用于信号去噪、频率分析和信号重建等方面。
本实验旨在通过设计和实现一个IIR滤波器，探究其在信号处理中的应用和性能。
一、背景知识1.1 数字滤波器数字滤波器是一种能够改变信号频率特性的系统，可以通过去除或增强特定频率的成分来实现信号处理的目的。
根据其传递函数的特点，数字滤波器可以分为FIR（有限脉冲响应）和IIR（无限脉冲响应）滤波器。
1.2 IIR滤波器IIR滤波器是一种具有无限脉冲响应的滤波器，其传递函数中包含了反馈回路。
相比于FIR滤波器，IIR滤波器具有更窄的转换带宽和更陡峭的滚降特性，能够更好地逼近理想滤波器的频率响应。
二、实验目的本实验旨在通过设计和实现一个IIR滤波器，探究其在信号处理中的应用和性能。
具体实验目标如下：1. 理解IIR滤波器的原理和设计方法；2. 掌握IIR滤波器的设计过程和参数选择；3. 分析IIR滤波器在不同输入信号下的性能表现。
三、实验设计与实施3.1 IIR滤波器的设计在本实验中，我们选择了巴特沃斯滤波器作为IIR滤波器的设计模型。
巴特沃斯滤波器具有最平坦的幅频响应特性和最小的群延迟，适用于许多实际应用场景。
首先，我们需要确定滤波器的阶数和截止频率。
阶数决定了滤波器的复杂度和性能，而截止频率则决定了滤波器的频率响应特性。
根据实际需求和信号特性，我们选择了一个二阶巴特沃斯低通滤波器，截止频率为1kHz。
其次，根据巴特沃斯滤波器的设计公式，我们可以计算出滤波器的传递函数和巴特沃斯极点的位置。
通过极点的选择和配置，我们可以调整滤波器的频率响应和滚降特性。
3.2 IIR滤波器的实施根据设计得到的传递函数和极点位置，我们可以使用MATLAB等工具进行IIR滤波器的实施和验证。
具体步骤如下：1. 根据巴特沃斯滤波器的传递函数公式，计算出滤波器的系数；2. 使用MATLAB的filter函数，将待处理的信号输入滤波器，得到滤波后的输出信号；3. 对比输入和输出信号的频谱特性，分析滤波器的性能。
实验四IIR数字滤波器设计实验报告
实验四IIR数字滤波器设计实验报告
为了实现信号的滤波处理，IIR（或称为滤波器）数字滤波器是一种常用的信号处理
技术。
本次实验就是探究IIR数字滤波器的设计和分析。
在实验开始前，对于IIR数字滤波器有所了解，它是一种无限级别功能的数字滤波器，其功能强大，可以实现任意自定义系数的滤波器。
在预处理实验中，便首先采用Matlab
工具搭建了IIR数字滤波器的框架，考虑到本次滤波处理内容，本次采用的是Chebyshev
类型的等离子体，其滤波效果要求超过50dB，进一步完善了对于设计工作的要求。
经过Chebyshev Type I等离子体的设计，确定了系统的结构，并设定了15个滤波器，接着从设定的各项参量入手，从而确定系统各项参量，运用梯形图确定根位置，并使用MATLAB中的filter函数进行系统模拟，得到经历处理后系统输出信号与未经处理时对比，结果显示滤波效果达到了相应预期要求。
在实验中，IIR数字滤波器的设计让我深刻体会到了系统滤波的重要性以及十分强大
的功能。
而它的实现，又显示了精确的数字处理技术在信号处理中的重要作用，使得研究
信号处理时，得以有效和准确地对信号进行分辨和滤波处理。
IIR数字滤波器设计实验报告
IIR数字滤波器的设计一、实验目的1、了解IIR数字滤波器的工作原理和作用2、掌握IIR数字滤波器的两种设计方法3、掌握使用MATLAB形成IIR数字滤波器二、实验内容有三首音乐，第一首为正常音质的音乐。
第二首为被加了紧邻原音乐的干扰的音乐。
第三首为被加了远离原音乐干扰的音乐。
要求设计IIR数字滤波器将被干扰的音乐恢复成不受干扰的音乐。
三、实验步骤步骤1: 将实际模拟低通滤波器指标转化为归一化模拟低通滤波器指标λs, αs, αp步骤2: 确定归一化模拟低通滤波器的系统函数Ha(p)步骤3: 由Ha(p)确定实际模拟低通滤波器的系统函数Ha(s)步骤4: 由Ha(s)确定的参数利用MATLAB形成IIR数字滤波器四、实验方法1、脉冲不变相应法：Matlab提供了脉冲不变响应法的库函数:[bz,az]=impinvar(b,a,Fs);表示将分子向量为b，分母向量为a的模拟滤波器通过脉冲响应不变法转换为分子向量为bz，分母向量为az的数字滤波器，采样频率为Fs，单位Hz。
2、双线性变换法：Matlab提供了双线性变换法的库函数:[bz,az]=bilinear(b,a,Fs);表示将分子向量为b，分母向量为a的模拟滤波器通过双线性变换法转换为分子向量为bz，分母向量为az的数字滤波器，采样频率为Fs，单位Hz。
五、实验程序与结果MATLAB代码：clear all;[s1,Fs,bits]=wavread('F:\music2-1.wav');s2=wavread('F:\music2-2.wav');s3=wavread('F:\music2-3.wav');t=(0:length(s1)-1)/Fs; % 计算数据时刻N=length(s1);if mod(N,2)==0;N=N;else s1(N)=[];N=N-1;end;fx=(0:N/2)*Fs/N;%%%%%%%%信号1%%%%%%%%figure(1);subplot(2,1,1);plot(t,s1); %绘制原音乐波形图subplot(2,1,2);s1f=fft(s1);plot(fx,abs(s1f(1:N/2+1)));%%%%%%%%信号2%%%%%%%%figure(2);subplot(2,1,1);plot(t,s2); % 绘制受紧邻原音乐的干扰的音乐的波形图subplot(2,1,2);s2f=fft(s2);plot(fx,abs(s2f(1:N/2+1)));%%%%%%%%信号3%%%%%%%%figure(3);subplot(2,1,1);plot(t,s3); %绘制受远离原音乐的干扰的音乐的波形图subplot(2,1,2);s3f=fft(s3);plot(fx,abs(s3f(1:N/2+1)));%%%%%%%%%滤波器设计1%%%%%%%Wp1=[2*8000*pi/Fs,2*10000*pi/Fs];Ws1=[2*8500*pi/Fs,2*9500*pi/Fs]; Rp=3;Rs=30;Wp11=2*Fs*tan(Wp1/2);Ws11=2*Fs*tan(Ws1/2);[N1,Wn1]=buttord(Wp11,Ws11,Rp,Rs,'s');[b,a]=butter(N1,Wn1,'stop','s');[bz,az]=bilinear(b,a,Fs);Y1=filter(bz,az,s2);figure(4);subplot(2,1,1);plot(t,Y1); %绘制滤波后波形图subplot(2,1,2);Yf1=fft(Y1);plot(abs(Yf1));wavwrite(Y1,Fs,bits,'F:\ music2-2lvbo.wav');%%%%%%%%%滤波器设计2%%%%%%%Wp2=2*7000*pi;Ws2=2*9000*pi;Rp2=3;Rs2=30;[N2,Wn2]=buttord(Wp2,Ws2,Rp2,Rs2,'s');[b2,a2]=butter(N2,Wn2,'s');[bz2,az2]=impinvar(b2,a2,Fs);Y2=filter(bz2,az2,s3);figure(5);subplot(2,1,1);plot(t,Y2); % 绘制滤波后波形图subplot(2,1,2);Yf2=fft(Y2);plot(abs(Yf2)); % 绘制滤波后波形图%%%%%%%%%%%%%%%%%%%%%%W=linspace(0,pi,pi*16000);Hz1=freqz(bz2,az2,W);Hz2=freqz(bz,az,W);figure(6);subplot(2,1,1);plot(abs(Hz1));subplot(2,1,2);plot(abs(Hz2));wavwrite(Y2,Fs,bits,'F:\music2-3lvbo.wav');。
数字信号处理实验报告四IIR数字滤波器设计及软件实现
数字信号处理实验报告四IIR数字滤波器设计及软件实现实验目的：本实验的目的是了解IIR数字滤波器的设计原理和实现方法，通过MATLAB软件进行数字滤波器设计和信号处理实验。
一、实验原理IIR数字滤波器是一种使用有限数量的输入样本和前一次输出值的滤波器。
它通常由差分方程和差分方程的系数表示。
IIR滤波器的特点是递归结构，故其频率响应是无限长的，也就是说它的频率响应在整个频率范围内都是存在的，而不像FIR滤波器那样只有在截止频率处才有响应。
根据设计要求选择合适的滤波器类型和滤波器结构，然后通过对滤波器的模型进行参数化，设计出满足滤波要求的IIR滤波器。
常见的IIR滤波器设计方法有模拟滤波器设计方法和数字滤波器设计方法。
在本实验中，我们主要使用数字滤波器设计方法，即离散时间滤波器设计方法。
二、实验内容（一）设计IIR数字滤波器的步骤：1.确定滤波器类型：根据滤波要求选择合适的滤波器类型，如低通滤波器、高通滤波器、带通滤波器、带阻滤波器等。
2.确定滤波器的阶数：根据滤波要求确定滤波器的阶数。
阶数越高，滤波器的频率响应越陡峭，但计算复杂度也越高。
3. 设计滤波器原型：根据滤波要求，设计滤波器的原型。
可以选择Butterworth滤波器、Chebyshev滤波器、Elliptic滤波器等作为原型。
4.选择滤波器结构：根据计算机实现条件和算法复杂度，选择合适的滤波器结构。
常见的滤波器结构有直接形式I、直接形式II、级联形式等。
5.参数化滤波器模型：根据原型滤波器的差分方程，选择合适的参数化方法。
常见的参数化方法有差分方程法、极点/零点法、增益法等。
6.根据参数化的滤波器模型，计算出所有的滤波器系数。
（二）用MATLAB软件实现IIR数字滤波器设计：1.打开MATLAB软件，并创建新的脚本文件。
2. 在脚本文件中，使用MATLAB提供的滤波器设计函数，如butter、cheby1、ellip等，选择合适的滤波器类型进行设计。
实验四IIR数字滤波器设计及软件实现实验报告(word文档)
实验四 IIR 数字滤波器设计及软件实现实验报告实验四 IIR 数字滤波器设计及软件实现实验报告一、实验目的（1）熟悉用双线性变换法设计IIR 数字滤波器的原理与方法；（2）学会调用 MATLAB信号办理工具箱中滤波器设计函数（或滤波器设计解析工具fdatool）设计各种IIR 数字滤波器，学会依照滤波需求确定滤波器指标参数。
（3）掌握 IIR 数字滤波器的 MATLAB实现方法。
（3）经过观察滤波器输入输出信号的时域波形及其频谱，建立数字滤波的看法。
二、实验原理设计 IIR 数字滤波器一般采用间接法（脉冲响应不变法和双线性变换法），应用最广泛的是双线性变换法。
基本设计过程是：①先将给定的数字滤波器的指标变换成过渡模拟滤波器的指标；②设计过渡模拟滤波器；③将过渡模拟滤波器系统函数变换成数字滤波器的系统函数。
MATLAB信号办理工具箱中的各种IIR 数字滤波器设计函数都是采用双线性变换法。
第六章介绍的滤波器设计函数butter、cheby1、cheby2和ellip能够分别被调用来直接设计巴特沃斯、切比雪夫1、切比雪夫 2 和椭圆模拟和数字滤波器。
本实验要求读者调用如上函数直接设计IIR 数字滤波器。
本实验的数字滤波器的MATLAB实现是指调用MATLAB信号办理工具箱函数filter对给定的输入信号x(n) 进行滤波，获取滤波后的输出信号y(n ）。
三、实验内容及步骤（1）调用信号产生函数mstg 产生由三路控制载波调幅信号相加构成的复合信号st ，该函数还会自动绘图显示st 的时域波形和幅频特点曲线，如图所示。
由图可见，三路信号时域混叠无法在时域分别。
但频域是分其他，因此能够经过滤波的方法在频域分别，这就是本实验的目的。
图三路调幅信号st 的时域波形和幅频特点曲线（ 2）要求将st 中三路调幅信号分别，经过观察st 的幅频特点曲线，分别确定能够分实验四 IIR 数字滤波器设计及软件实现实验报告离 st 中三路控制载波单频调幅信号的三个滤波器 （低通滤波器、 带通滤波器、 高通滤波器）的通带截止频率和阻带截止频率。
iir滤波器设计实验报告
iir滤波器设计实验报告IIR滤波器设计实验报告一、实验目的本实验旨在通过实际操作，掌握IIR滤波器的基本设计方法，了解滤波器性能参数对滤波效果的影响，加深对滤波器理论的理解。
二、实验原理IIR滤波器（Infinite Impulse Response）是一种离散时间滤波器，其系统函数具有无限长的时间响应。
IIR滤波器设计方法主要包括冲激响应不变法和双线性变换法。
本实验采用冲激响应不变法进行设计。
三、实验步骤1. 确定滤波器性能参数：根据实际需求，确定滤波器的类型（低通、高通、带通、带阻）、通带边缘频率、阻带边缘频率、通带波动和阻带衰减等性能参数。
2. 计算滤波器系数：根据冲激响应不变法，利用Matlab编程计算滤波器系数。
具体过程包括定义系统函数、计算冲激响应、计算频率响应等步骤。
3. 编写滤波器程序：根据计算出的滤波器系数，编写IIR滤波器程序。
程序应实现输入信号的滤波处理，并输出滤波后的信号。
4. 测试滤波器性能：对编写的滤波器程序进行测试，观察其滤波效果，分析性能参数对滤波效果的影响。
5. 优化滤波器性能：根据测试结果，对滤波器性能参数进行调整，优化滤波效果。
四、实验结果及分析通过本次实验，我们成功地设计并实现了IIR滤波器。
在测试过程中，我们观察到了滤波器对不同频率信号的过滤效果，并分析了性能参数对滤波效果的影响。
具体来说，通带边缘频率决定了滤波器对低频信号的过滤程度，阻带边缘频率则影响对高频信号的过滤程度。
通带波动和阻带衰减则分别反映了滤波器在通带和阻带的波动程度和衰减程度。
通过对这些性能参数的调整，我们可以实现对不同类型信号的有效过滤。
五、实验总结通过本次实验，我们深入理解了IIR滤波器的工作原理和设计方法，掌握了Matlab编程在滤波器设计中的应用。
实验过程中，我们不仅学会了如何根据实际需求选择合适的性能参数，还学会了如何调整这些参数以优化滤波效果。
此外，我们还观察到了不同性能参数对滤波效果的影响，加深了对滤波器理论的理解。
实验四 IIR数字滤波器设计
实验四 IIR数字滤波器设计一、实验目的1、熟悉用双线性变换法设计IIR数字滤波器的原理与方法；2、熟悉用脉冲响应不变法设计IIR数字滤波器的原理与方法。
二、实验内容1、已知低通滤波器的指标为：通带边缘频率：0.4π，Ap=0.5dB；阻带边缘频率：0.6π，As=50dB；a、采用脉冲响应不变法设计巴特沃斯，T=1.画出幅度响应和数字滤波器的脉冲响应h（n）；b、b、采用脉冲响应不变法设计巴特沃斯，T=1.画出幅度响应和数字滤波器的脉冲响应h（n）。
2、用双线性变换法设计低通滤波器，满足技术指标wp=0.2π，Ap=0.25dB；ws=0.4π，As=50dB，并对方波信号进行滤波，画出滤波前后的波形图并进行简要分析3、设计一个数字高通滤波器H（z），它用在结构xa（t）A/D H（z) D/A ya（t）中，满足下列要求：a、采样速率为10kHZ；b、阻带边缘频率为1.5kHZ，衰减为40dB；c、通带边缘频率为2kHZ，衰减为3dB；d、单调的通带和阻带。
4、设计一个带阻滤波器，要求通带上下截止频率为0.8π，0.2π，通带内衰减不大于1dB，阻带起始频率为0.7π，0.4π，阻带内衰减不小于30dB。
设计巴特沃斯带阻滤波器并画出该数字高通滤波器的幅度响应和脉冲响应。
三、实验程序及解释和实验分析及图形1a、clear;close all; %清屏wp=0.4*pi;ws=0.6*pi;Ap=0.5As=50;T=1;Fs=1/T;OmegaP=wp/T;OmegaS=ws/T;[cs,ds]=afd_butt(OmegaP,OmegaS,Ap,As);[b,a]=impinvar(cs,ds,Fs);[h,w]=freqz(b,a);subplot(2,2,1);plot(w/pi,abs(h));title('幅度响应');grid;subplot(2,2,2);plot(w/pi,angle(h));title('相位响应');grid;subplot(2,2,2);plot(w/pi,angle(h));title('相位响应');grid;subplot(2,2,3);plot(w/pi,20*log(abs(h)));title('幅度响应dB');grid;n=[0:1:59];imp=[1;zeros(59,1)];y=filter(b,a,imp);subplot(2,2,4);plot(n,y);title('脉冲响应');grid;分析：因为w=Ωt,w与Ω呈线性关系，所以其相位响应图是线性的，如相位响应图所示；因为其设计的是低通滤波器，所以会把其高频部分滤掉，留下低频，从而不是特别陡峭，如脉冲响应图所示。
matlab实验报告 IIR数字滤波器设计
实验报告姓名：李鹏博 实验名称： IIR 数字滤波器设计 学号：2011300704 课程名称： 数字信号处理 班级：03041102 实验室名称： 航海西楼303 组号： 1 实验日期： 2014.06.20一、实验目的、要求掌握IIR 数字滤波器设计的冲激响应不变法和双线性变换法。
掌握IIR 数字滤波器的计算机编程实现方法，即软件实现。
二、实验原理为了从模拟滤波器设计IIR 数字滤波器，必须先设计一个满足技术指标的模拟滤波器，然后将其数字化，即从s 平面映射到z 平面，得到所需的数字滤波器。
虽然IIR 数字滤波器的设计本质上并不取决于连续时间滤波器的设计，但是因为在许多应用中，数字滤波器就是用来模仿模拟滤波器功能的，所以由模拟滤波器转化为数字滤波器是很自然的。
因此，由模拟滤波器设计数字滤波器的方法准确、简便，是目前最普遍采用的方法。
三、实验环境PC 机，Windows XP ，office 2003，Matlab 软件。
四、实验过程、数据记录、分析及结论实验过程1．编程设计滤波器，用冲激响应不变法设计IIR 数字滤波器。
2．编程设计滤波器，用双线性变换法设计IIR 数字滤波器。
3．求脉冲响应、频率响应以及零极点。
4．编程滤波，求滤波器输出，完成对不同频率的多个正弦信号的滤波。
实验步骤根据所给定的技术指标进行指标转换。
112c c f πΩ=，222c c f πΩ=，112s s f πΩ=，222s s f πΩ=，21p c c B Ω==Ω-Ω，221222s s s s s B Ω-ΩΩΩ=Ω，3,18p s αα=-=-。
根据指标设计Butterworth 模拟低通滤波器。
调用函数[n,wn]=buttord(wp,ws,rp,rs,’s ’)确定阶次。
调用函数[zl,pl,kl]=buttap(n)，求低通原型的模型。
调用函数[bl,al]=zp2tf(zl,pl,kl)实现模型转换。
实验四 IIR数字滤波器的设计实验报告
数字信号处理实验报告实验四 IIR数字滤波器的设计学生姓名张志翔班级电子信息工程1203班学号***********指导教师2015.4.29实验四 IIR 数字滤波器的设计一、实验目的：1. 掌握双线性变换法及脉冲响应不变法设计IIR 数字滤波器的具体设计方法及其原理，熟悉用双线性变换法及脉冲响应不变法设计低通、高通和带通IIR 数字滤波器的MATLAB 编程。
2. 观察双线性变换及脉冲响应不变法设计的滤波器的频域特性，了解双线性变换法及脉冲响应不变法的特点。
3. 熟悉Butterworth 滤波器、切比雪夫滤波器和椭圆滤波器的频率特性。
二、实验原理: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） =0.3KHz, δ=0.8Db, =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--------++幅频响应图：分析：由图可知，切比雪夫滤波器幅频响应是通带波纹，阻带单调衰减的。
实验四IIR数字滤波器设计及软件实现实验报告
实验四IIR数字滤波器设计及软件实现实验报告一、实验目的（1）熟悉用双线性变换法设计IIR数字滤波器的原理与方法；（2）学会调用MATLAB信号处理工具箱中滤波器设计函数（或滤波器设计分析工具fdatool）设计各种IIR数字滤波器，学会根据滤波需求确定滤波器指标参数。
（3）掌握IIR数字滤波器的MATLAB实现方法。
（3）通过观察滤波器输入输出信号的时域波形及其频谱，建立数字滤波的概念。
二、实验原理设计IIR数字滤波器一般采用间接法（脉冲响应不变法和双线性变换法），应用最广泛的是双线性变换法。
基本设计过程是：①先将给定的数字滤波器的指标转换成过渡模拟滤波器的指标；②设计过渡模拟滤波器；③将过渡模拟滤波器系统函数转换成数字滤波器的系统函数。
MATLAB信号处理工具箱中的各种IIR数字滤波器设计函数都是采用双线性变换法。
第六章介绍的滤波器设计函数butter、cheby1 、cheby2 和ellip可以分别被调用来直接设计巴特沃斯、切比雪夫1、切比雪夫2和椭圆模拟和数字滤波器。
本实验要求读者调用如上函数直接设计IIR数字滤波器。
本实验的数字滤波器的MATLAB实现是指调用MATLAB信号处理工具箱函数filter对给定的输入信号x(n)进行滤波，得到滤波后的输出信号y(n）。
三、实验内容及步骤（1）调用信号产生函数mstg产生由三路抑制载波调幅信号相加构成的复合信号st，该函数还会自动绘图显示st的时域波形和幅频特性曲线，如图10.4.1所示。
由图可见，三路信号时域混叠无法在时域分离。
但频域是分离的，所以可以通过滤波的方法在频域分离，这就是本实验的目的。
图10.4.1三路调幅信号st的时域波形和幅频特性曲线（2）要求将st中三路调幅信号分离，通过观察st的幅频特性曲线，分别确定可以分离st中三路抑制载波单频调幅信号的三个滤波器（低通滤波器、带通滤波器、高通滤波器）的通带截止频率和阻带截止频率。
要求滤波器的通带最大衰减为0.1dB,阻带最小衰减为60dB。
实验四 IIR数字滤波器的设计与滤波
实验四 IIR 数字滤波器的设计与滤波一、巴特沃斯模拟滤波器的设计1． 模拟滤波器的设计参数模拟滤波器的4个重要的通带、阻带参数为：p f 或Omegap ：通带截止频率 s f 或Omegas ：阻带截至频率p R ：通带内波动(dB)，即通带内所允许的最大衰减；s R ：阻带内最小衰减通过以上参数就可以进行模拟滤波器的设计。
2． 巴特沃斯模拟滤波器设计1） 巴特沃斯滤波器阶数的选择：在已知设计参数p f ，s f ，p R ，s R 之后，可利用“buttord ”命令可求出所需要的滤波器的阶数和3dB 截止频率，其格式为：[N ，Omegac]=buttord[fp ，fs ，Rp ，Rs ，’s ’]，其中fp ，fs ，Rp ，Rs 分别为通带截止频率、阻带起始频率、通带内波动、阻带内最小衰减。
返回值N 为滤波器的最低阶数，Wc 为3dB 截止频率。
2） 巴特沃斯滤波器系数计算：由巴特沃斯滤波器的阶数N 以及3dB 截止频率Omegac 可以计算出对应传递函数H(z)的分子分母系数，MATLAB 提供的命令如下：● 巴特沃斯低通滤波器系数计算：[b ，a]=butter(N,Omegac)，其中b 为H(z)的分子多项式系数，a 为H(z)的分母多项式系数● 巴特沃斯高通滤波器系数计算：[b ，a]=butter(N,Omegac,’High ’)● 巴特沃斯带通滤波器系数计算：[b ，a]=butter(N ，[Omega1，Omega2])，其中[Omega1，Omega2]为通带截止频率，是2元向量，需要注意的是该函数返回的是2N 阶滤波器系数。
● 巴特沃斯带阻滤波器系数计算：[b ，a]=butter(N ，[Omega1，Omega2]，’stop ’)，其中[Omega1，Omega2]为通带截止频率，是2元向量，需要注意的是该函数返回的也是2N 阶滤波器系数。
二、巴特沃斯数字滤波器的设计1. 数字滤波器的设计参数滤波器的4个重要的通带、阻带参数为：p f ：通带截止频率（Hz ） s f ：阻带起始频率（Hz ）p R ：通带内波动（dB ），即通带内所允许的最大衰减； s R ：阻带内最小衰减设采样速率（即奈奎斯特速率）为N f ，将上述参数中的频率参数转化为归一化频率参数：p ω：归一化通带截止频率，)2//(N p p f f =ω；s ω：归一化阻带截至频率，)2//(N s s f f =ω通过以上参数就可以进行数字滤波器的设计。
iir数字滤波器的设计实验报告
iir数字滤波器的设计实验报告IIR数字滤波器的设计实验报告引言数字滤波器是数字信号处理中的重要组成部分，用于去除信号中的噪声、滤波、频率分析等。
在数字滤波器中，IIR（Infinite Impulse Response）滤波器是一种常见且广泛应用的滤波器类型。
本实验旨在设计一个IIR数字滤波器，并通过实验验证其性能。
一、实验目的本实验的目标是设计一个IIR数字滤波器，实现对输入信号的滤波功能。
具体而言，我们将通过以下步骤完成实验：1. 确定滤波器的滤波类型（低通、高通、带通或带阻）和截止频率。
2. 设计滤波器的传递函数。
3. 使用Matlab或其他数学软件进行滤波器的频率响应和时域响应分析。
4. 利用实验数据对滤波器进行性能评估。
二、实验原理IIR数字滤波器的设计基于差分方程，其传递函数可以表示为：H(z) = (b0 + b1*z^(-1) + b2*z^(-2) + ... + bn*z^(-n)) / (1 + a1*z^(-1) +a2*z^(-2) + ... + am*z^(-m))其中，b0、b1、...、bn和a1、a2、...、am是滤波器的系数。
滤波器的阶数为max(m, n)。
根据滤波器的滤波类型和截止频率，可以确定这些系数的具体值。
三、实验步骤1. 确定滤波器的类型和截止频率。
例如，我们选择设计一个低通滤波器，截止频率为1kHz。
2. 根据所选滤波器类型和截止频率，计算滤波器的传递函数。
3. 使用Matlab或其他数学软件进行滤波器的频率响应和时域响应分析。
可以绘制滤波器的幅频响应曲线和相频响应曲线，以及滤波后的信号波形。
4. 利用实验数据对滤波器进行性能评估。
可以通过输入不同频率的信号，观察滤波器的效果，并计算滤波器的截止频率、增益和相位特性等参数。
四、实验结果与分析通过实验，我们得到了设计的低通滤波器的频率响应和时域响应曲线。
在频率响应曲线中，我们可以观察到滤波器在截止频率附近的衰减特性，以及在截止频率以下的通过特性。
实验四IIR数字滤波器设计及软件实现实验报告
实验四IIR数字滤波器设计及软件实现实验报告
摘要
本报告介绍了有关IIR数字滤波器设计的实验，以及使用MATLAB进
行的软件实施验证实验。
实验结果表明，IIR滤波器的设计和实施过程中，模糊C不做任何处理，也能实现意料之外的良好滤波效果。
1.介绍
本文介绍了实验四的IIR数字滤波器设计与软件实现实验。
在完成本
实验之前，学生完成了实验一，实验二和实验三，分别设计了低通滤波器、带通滤波器和高通滤波器。
在本实验中，学生将总结前三个实验的知识，
设计和实施一个二阶高通滤波器，以及一个四阶带阻滤波器。
2.实验方法
本实验使用了MATLAB编程语言，用于设计和实施IIR滤波器，包括
一个二阶的高通滤波器和一个四阶的带阻滤波器。
首先，选择预定义的滤
波器系统函数，并调整其参数，以实现特定的滤波器性能。
然后，针对调
整好的滤波器，编写MATLAB代码，实施设计的滤波器。
3.实验结果
（1）二阶高通滤波器
二阶高通滤波器的设计参数如下：
参数，值
-----------------，----------
截止频率，0.25Hz
最小插入损耗，0dB 最大衰减率，40dB。
IIR数字滤波器的设计实验报告
IIR数字滤波器的设计一、实验目的：掌握冲激相应不变法和双线性变换法设计IIR数字滤波器的原理和方法；观察冲激相应不变法和双线性变换法设计IIR数字滤波器的频率特性；了解冲激相应不变法和双线性变换法的特点和区别。
二、实验原理：无限长单位冲激响应（IIR）数字滤波器的设计思想：a)设计一个合适的模拟滤波器b)利用一定的变换方法将模拟滤波器转换成满足预定指标的数字滤波器切贝雪夫I型：通带中是等波纹的，阻带是单调的切贝雪夫II型：通带中是单调的，阻带是等波纹的1．用冲击响应不变法设计一个低通切贝雪夫I型数字滤波器通带上限截止频率为400Hz阻带截止频率为600Hz通带最大衰减为0.3分贝阻带最小衰减为60分贝抽样频率1000Hz2．用双线性变换法设计切贝雪夫II型高通滤波器通带截止频率2000Hz阻带截止频率1500Hz通带最大衰减0.3分贝阻带最小衰减50分贝抽样频率20000Hz四、实验程序：1）Wp=2*pi*400;Ws=2*pi*600;Rp=0.3;Rs=60;Fs=1000;[N,Wn]=cheb1ord(Wp,Ws,Rp,Rs,'s'); [Z,P,K]=cheb1ap(N,Rp);[A,B,C,D]=zp2ss(Z,P,K);[At,Bt,Ct,Dt]=lp2lp(A,B,C,D,Wn); [num1,den1]=ss2tf(At,Bt,Ct,Dt);[num2,den2]=impinvar(num1,den1,Fs); [H,W1]=freqs(num1,den1);figure(1)subplot(2,1,1);semilogx(W1/pi/2,20*log10(abs(H)));grid; xlabel(' 频率/ Hz');ylabel(' 模拟滤波器幅值(db)');[H,W2]=freqz(num2,den2,512,'whole',Fs); subplot(2,1,2);plot(W2,20*log10(abs(H)));grid;xlabel(' 频率/ Hz');ylabel(' 数字滤波器幅值(db)'); [H,W]=freqz(num2,den2,512,'whole'); figure(2)subplot(2,1,1);8plot(W/2/pi,20*log10(abs(H)));grid; xlabel(' 数字角频率/ pi');ylabel(' 数字滤波器幅值(db)'); %xlabel(' 频率/ Hz'); %ylabel(' 幅值'); subplot(2,1,2);%stem(W/pi,abs(H));grid;%figure(3)%semilogx(W1,W2);%grid;plot(W/2/pi,angle(H)/pi);grid;xlabel(' 数字角频率/ pi');ylabel(' 相角/ pi');2）Wp=2*pi*2000;Ws=2*pi*1500;Rp=0.3;Rs=50;Fs=20000;[N,Wn]=cheb2ord(Wp,Ws,Rp,Rs,'s'); [Z,P,K]=cheb2ap(N,Rs);[A,B,C,D]=zp2ss(Z,P,K);[At,Bt,Ct,Dt]=lp2hp(A,B,C,D,Wn); [num1,den1]=ss2tf(At,Bt,Ct,Dt); [num2,den2]=bilinear(num1,den1,Fs); [H,W]=freqz(num2,den2);subplot(2,1,1);plot(W*Fs/2/pi,20*log10(abs(H)));grid; xlabel(' 数字角频率/ Hz');ylabel(' 幅值(dby)');subplot(2,1,2);plot(W*Fs/2/pi,angle(H)/pi);grid; xlabel(' 数字角频率/ Hz');ylabel(' 相角/ pi');五、仿真图形：六、实验分析：1.在第一个实验中模拟滤波器的频率响应曲线延伸无穷，而数字滤波器只给出500Hz一下的频谱是什么原因？答：模拟滤波器的频率响应从负无穷至正无穷，经过采样后，变到（0,2π）。
iir数字滤波器实验报告
iir数字滤波器实验报告IIR数字滤波器实验报告引言：数字滤波器是数字信号处理中重要的组成部分，它可以对信号进行滤波和去噪，提取出我们所需要的信息。
在本次实验中，我们将重点研究和实验IIR数字滤波器的性能和应用。
一、IIR数字滤波器的原理IIR（Infinite Impulse Response）数字滤波器是一种递归滤波器，它的输出不仅与当前输入有关，还与之前的输入和输出有关。
IIR滤波器的结构可以由巴特沃斯、切比雪夫等滤波器设计方法得到。
与FIR（Finite Impulse Response）数字滤波器相比，IIR滤波器具有更低的计算复杂度和更好的频率响应特性。
二、IIR数字滤波器的设计在本次实验中，我们选择了巴特沃斯滤波器作为IIR滤波器的设计方法。
巴特沃斯滤波器是一种理想的低通滤波器，具有平坦的通带和陡峭的阻带。
通过选择不同的阶数和截止频率，我们可以得到不同性能的滤波器。
三、IIR数字滤波器的性能评估为了评估IIR数字滤波器的性能，我们进行了一系列实验。
首先，我们使用MATLAB软件进行了滤波器的设计和模拟。
通过绘制滤波器的频率响应曲线和幅度响应曲线，我们可以直观地了解滤波器的性能。
其次，我们使用真实的信号进行了滤波实验。
通过对比滤波前后信号的波形和频谱图，我们可以评估滤波器的去噪和频率特性。
四、IIR数字滤波器的应用IIR数字滤波器在实际应用中具有广泛的用途。
例如，语音信号处理中常用的降噪算法就是基于IIR滤波器的。
此外，IIR滤波器还可以用于信号增强、图像处理等领域。
通过调整滤波器的参数，我们可以实现不同的滤波效果，满足不同应用场景的需求。
五、实验结果与讨论在本次实验中，我们设计了一个二阶巴特沃斯低通滤波器，截止频率为1kHz。
通过MATLAB软件模拟和实际信号滤波实验，我们得到了滤波前后信号的波形和频谱图。
实验结果表明，滤波器能够有效地去除高频噪声，并保留低频信号的主要成分。
同时，滤波器的频率响应也符合设计要求，具有良好的通带和阻带特性。
实验四 IIR数字滤波器设计
图I 5阶Butterworth 数字高通滤波器试验四IIR 数字滤波器的设计与MATLAB 实现一、试验目的：1、要求把握∏R 数字滤波器的设计原理、方法、步骤。
2、能够依据滤波器设计指标进行滤波器设计。
3、把握数字巴特沃斯滤波器和数字切比雪夫滤波器的设计原理和步骤。
二、试验原理：∏R 数字滤波器的设计方法：频率变换法、数字域直接设计以及计算机帮助等。
这里只介绍频率变换法。
由模拟低通滤波器到数字低通滤波器的转换，基本设计 过程：1、将数字滤波器的设计指标转换为模拟滤波器指标2、设计模拟滤波器G （S ）3、将G （S ）转换为数字滤波器H （Z ）在低通滤波器设计基础上，可以得到数字高通、带通、带阻滤波器的设计流程如 下：1、给定数字滤波器的设计要求（高通、带通、带阻）2、转换为模拟（高通、带通、带阻）滤波器的技术指标3、转换为模拟低通滤波器的指标4、设计得到满意3步骤中要求的低通滤波器传递函数5、通过频率转换得到模拟（高通、带通、带阻）滤波器6、变换为数字（高通、带通、带阻）滤波器三、标准数字滤波器设计函数MATLAB 供应了一组标准的数字滤波器设计函数，大大简化了滤波器设计过程。
1 > butter例题1设计一个5阶Butterworth 数字高通滤波器，阻带截止频率为250Hz ,设 采样频率为IKHz.I k H J-∣H ∏ t er (5. 250/500.' high')L z, ∣>, kJ but i er(5t 250 500, , ∣∣ i glιt)f r eqz (b 1 5 I 2, I 000)50 100 150 200 250 300 350 400 450 500 Frequency (Hz) o o o o opo 1 3 in 3 3w=⅛e2 50 100 150 200 250 300 350 400 450 500 Fιequetιcy (Hz) - A ・ > A ・o o o o o o o o o 力 o o 1 -23 < 京⅛cy.⅛)φseud2、chebyl 和cheby2例题2设，十一个7阶chebyshevll型数字低通滤波器，截止频率为3000Hz,Rs=30dB,采样频率为IKHz。
实验四 IIR数字滤波器的设计
实验四 IIR数字滤波器的设计----b1a2bc00-6eaf-11ec-9f14-7cb59b590d7d实验四-iir数字滤波器的设计实验四 IIR数字滤波器的设计1。
实验目的（1）掌握脉冲响应不变法和双线性变换法设计iir数字滤波器的具体方法和原理，熟悉双线性变换法和脉冲响应不变法设计低通、带通iir数字滤波器的计算机编程；（2）观察双线性变换法和脉冲响应不变法设计的数字滤波器的频域特性，了解双线性变换法和脉冲响应不变法的特点和区别；（3）熟悉butterworth滤波器、chebyshev滤波器和椭圆滤波器的频率特性。
2.实验原理和方法iir数字滤波器的设计方法可以概括为如图所示，本实验主要掌握iir滤波器的第一种方法，即利用模拟滤波器设计iir数字滤波器，这是iir数字滤波器设计最常用的方法。
利用模拟滤波器设计，需要将模拟域的ha(s)转换为数字域h(z)，最常用的转换方法为脉冲响应不变法和双线性变换法。
IIR数字滤波器设计方法由模拟滤波器设计IIR数字滤波器线性设计IIR数字滤波器冲激响应不变量法双线性变换法零极点累积测试法频域近似法时域近似法（1）冲激响应不变量法用数字滤波器的单位脉冲响应序列h(n)模仿模拟滤波器的冲激响应ha(t)，让h(n)正好等于ha(t)的采样值，即h（n）？医管局（新界）其中t为采样间隔。
如果以ha(s)及h(z)分别表示ha(t)的拉氏变换及h(n)的z变换，则h（z）|z？美国东部时间？1吨？？Kha（s？j2？tk）在matlab中，可用函数impinvar实现从模拟滤波器到数字滤波器的脉冲响应不变映射。
（2）双线性变换法s平面与z平面之间满足下列映射关系2秒？21? zt1？Z1.1还是Z？t2t？s?ss平面的虚轴单值地映射于z平面的单位圆上，s平面的左半平面完全映射到z平面的单位圆内。
双线性变换不存在频率混叠问题。
在MATLAB中，通过函数双线性可以实现模拟滤波器到数字滤波器的双线性变换映射。
无限冲激响应(IIR)数字滤波器的设计实验报告
数字信号处理实验报告姓名：寇新颖 学号：20100304026 专业：电子信息科学与技术实验四 无限冲激响应(IIR)数字滤波器的设计一、实验目的1．掌握双线性变换法及冲激响应不变法设计IIR 数字滤波器的具体设计方法及其原理，熟悉用双线性变换法及冲激响应不变法设计低通IIR 数字滤波器的计算机编程。
2．观察双线性变换及冲激响应不变法设计的滤波器的频域特性，了解双线性变换法及冲激响应不变法的特点。
3．熟悉Butterworth 滤波器的频率特性。
二、实验原理1．利用模拟滤波器设计IIR 数字滤波器方法（1）根据所给出的数字滤波器性能指标计算出相应的模拟滤波器的设计指标。
（2）根据得出的滤波器性能指标设计出相应的模拟滤波器的系统函数H(S)。
（3）根据得出的模拟滤波器的系统函数H(S)，经某种变换得到对该模拟滤波器相应的数字仿真系统——数字滤波器。
将模拟滤波器转换成数字滤波器的实质是，用一种从s 平面到z 平面的映射函数将Ha(s)转换成H(z)。
对这种映射函数的要求是：(1) 因果稳定的模拟滤波器转换成数字滤波器，仍是因果稳定的。
(2)数字滤波器的频率响应模仿模拟滤波器的频响，s 平面的虚轴映射z 平面的单位圆，相应的频率之间成线性关系。
冲激响应不变法和双线性变换法都满足如上要求。
2．冲激响应不变法用数字滤波器的单位脉冲响应序列h(n)模仿模拟滤波器的冲激响应h a (t),让h(n)正好等于h a (t)的采样值，即h(n)=h a (nT)，其中T 为采样间隔。
3．双线性变换法s 平面与z 平面之间满足以下映射关系：1111--+-=z z s s 平面的虚轴单值地映射于z 平面的单位圆上，s 平面的左半平面完全映射到z 平面的单位圆内。
双线性变换不存在混叠问题。
双线性变换时一种非线性变换)2/(ωtg =Ω，这种非线性引起的幅频特性畸变可通过预畸而得到校正。
以低通数字滤波器为例，将设计步骤归纳如下：(1)确定数字滤波器的性能指标：通带临界频率f p 、阻带临界频率f s ；通带内的最大衰减A p ；阻带内的最小衰减A s ；(2)确定相应的数字角频率，ωp =2πf p ；ωs =2πf s ；(3)计算经过预畸的相应模拟低通原型的频率，)2/(ωtg =Ω；(4)根据Ωp 和Ωs 计算模拟低通原型滤波器的阶数N ，并求得低通原型的传递函数H a (s)；(5)用上面的双线性变换公式代入H a (s)，求出所设计的传递函数H(z)；(6)分析滤波器特性，检查其指标是否满足要求。
实验四IIR数字滤波器的设计实验报告
实验四IIR数字滤波器的设计实验报告实验四：IIR数字滤波器的设计实验目的：1.了解IIR数字滤波器的基本原理和设计流程；2.学习使用MATLAB进行IIR数字滤波器的设计；3.实际设计一个IIR数字滤波器，并对输入信号进行滤波处理。
实验设备：1.计算机2.MATLAB软件实验原理：IIR数字滤波器是一种非线性滤波器，可以通过差分方程的形式表示。
其特点是具有无穷长的单位脉冲响应，即滤波器对输入信号的响应是无限长的。
IIR数字滤波器的设计一般包括两个方面：滤波器的结构和滤波器的参数。
其中，滤波器的结构包括滤波器的拓扑结构和级联结构，滤波器的参数包括滤波器的截止频率、通带增益、阻带衰减等。
实验步骤：1.确定滤波器的类型（低通滤波器、高通滤波器、带通滤波器等）；2.根据滤波器的要求，设计滤波器的截止频率、通带增益、阻带衰减等参数；3.根据滤波器的类型和参数，选择合适的滤波器结构和滤波器参数；4.使用MATLAB软件进行IIR数字滤波器的设计，编写相应的代码；5.载入输入信号，并对输入信号进行滤波处理；6.分析输出信号的频谱特性和时域波形。
实验结果：通过实验，我们成功设计了一个IIR数字滤波器，并对输入信号进行了滤波处理。
实验结果显示，滤波器能够有效地去除输入信号中的高频噪声，得到了更清晰的输出信号。
输出信号的频谱特性和时域波形符合设计要求。
实验结论：IIR数字滤波器是一种常用的数字滤波器，具有较好的滤波效果和较低的计算复杂度。
通过实验，我们深入了解了IIR数字滤波器的设计原理和流程，并成功应用于实际信号处理中。
实验结果表明，IIR数字滤波器能够有效地去除输入信号中的噪声，提取出所需的信号信息。
这对于信号处理和通信系统设计具有重要意义。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
数字信号处理实验报告实验四 IIR数字滤波器的设计学生姓名张志翔班级电子信息工程1203班学号指导教师实验四 IIR数字滤波器的设计一、实验目的：1. 掌握双线性变换法及脉冲响应不变法设计IIR 数字滤波器的具体设计方法及其原理，熟悉用双线性变换法及脉冲响应不变法设计低通、高通和带通IIR 数字滤波器的MATLAB 编程。
2. 观察双线性变换及脉冲响应不变法设计的滤波器的频域特性，了解双线性变换法及脉冲响应不变法的特点。
3. 熟悉Butterworth 滤波器、切比雪夫滤波器和椭圆滤波器的频率特性。
二、实验原理: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） =, δ=, =, 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,,20,'s'); %给定通带（wp）和阻带(ws)边界角频率，通带波动波动，阻带最小衰减20dB，求出最低阶数和通带滤波器的通带边界频率Wn[B,A]=cheby1(N,,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 =den = 1系统函数：1234 12340.0304 -0.1218z 0.1827z-0.1218z0.0304zH(z)=1.0000+1.3834z+1.4721z+ 0.8012z+0.2286z--------++幅频响应图：分析：由图可知，切比雪夫滤波器幅频响应是通带波纹，阻带单调衰减的。
δ=，fr=,At=30Db,满足设计要求（2）fc=,δ=1dB,fr=,At=25dB,T=1ms;分别用脉冲响应不变法及双线性变换法设计一Butterworth数字低通滤波器，观察所设计数字滤波器的幅频特性曲线，记录带宽和衰减量，检查是否满足要求。
比较这两种方法的优缺点。
MATLAB源程序：T = ;fs = 1000;fc = 200;fr = 300;wp1 = 2*pi*fc;wr1 = 2*pi*fr;[N1,wn1] = buttord(wp1,wr1,1,25,'s')[B1,A1] = butter(N1,wn1,'s');[num1,den1] = impinvar(B1,A1,fs);%脉冲响应不变法[h1,w] = freqz(num1,den1);wp2 = 2*fs*tan(2*pi*fc/(2*fs))wr2 = 2*fs*tan(2*pi*fr/(2*fs))[N2,wn2] = buttord(wp2,wr2,1,25,'s')[B2,A2] = butter(N2,wn2,'s');[num2,den2] = bilinear(B2,A2,fs);%双线性变换法[h2,w] = freqz(num2,den2);f = w/(2*pi)*fs;plot(f,20*log10(abs(h1)),'-.',f,20*log10(abs(h2)),'-');axis([0,500,-100,10]);grid;xlabel('频率/Hz ');ylabel('幅度/dB')title('巴特沃思数字低通滤波器');legend('脉冲相应不变法','双线性变换法',1);结果分析：脉冲响应不变法的低通滤波器系统函数:num1 0den1 1987654321876543210.0004z 0.0060z 0.0450z 0.2045z 0.6309z 1.3869z 2.2053z 2.5324z 1.9199z 1 3.6569z 0.0002z 0.0075z 0.0611z 0.1444z 0.0995z 0.0153z 0.0002z 2.3647------------------+-+-+-+-++++++++-=)(z H双线性变换法设计的低通滤波器系统函数:num2den2 16543216543210025.00208.01501.02989.09130.06019.010176.01072.02681.03575.02681.00.10720179.0)(------------+-+-+-++++++=z z z z z z z z z z z z z H分析：脉冲响应不变法的频率变化是线性的，数字滤波器频谱响应出现了混叠，影响了过渡带的衰减特性，并且无传输零点；双线性变化法的频率响应是非线性的，因而消除了频谱混叠，在f=500Hz出有一个传输零点。
脉冲响应不变法的一个重要特点是频率坐标的变换是线性的，ω＝ΩΤ，ω与Ω是线性关系：在某些场合，要求数字滤波器在时域上能模仿模拟滤波器的功能时，如要实现时域冲激响应的模仿，一般使用脉冲响应不变法。
脉冲响应不变法的最大缺点：有频谱周期延拓效应，因此只能用于带限的频响特性，如衰减特性很好的低通或带通,而高频衰减越大，频响的混淆效应越小，至于高通和带阻滤波器,由于它们在高频部分不衰减，因此将完全混淆在低频响应中，此时可增加一保护滤波器，滤掉高于的频带，再用脉冲响应不变法转换为数字滤波器，这会增加设计的复杂性和滤波器阶数，只有在一定要满足频率线性关系或保持网络瞬态响应时才采用。
双线性变换法的主要优点是S平面与Z平面一一单值对应，s平面的虚轴(整个jΩ)对应于Z平面单位圆的一周，S平面的Ω=0处对应于Z平面的ω=0处，Ω= ∞处对应于Z平面的ω= π处,即数字滤波器的频率响应终止于折叠频率处，所以双线性变换不存在混迭效应。
双线性变换缺点: Ω与ω成非线性关系，导致：a. 数字滤波器的幅频响应相对于模拟滤波器的幅频响应有畸变，(使数字滤波器与模拟滤波器在响应与频率的对应关系上发生畸变)。
b. 线性相位模拟滤波器经双线性变换后，得到的数字滤波器为非线性相位。
c.要求模拟滤波器的幅频响应必须是分段恒定的，故双线性变换只能用于设计低通、高通、带通、带阻等选频滤波器。
(3)利用双线性变换法分别设计满足下列指标的Butterworth型、Chebyshev型和椭圆型数字低通滤波器，并作图验证设计结果：fc= ，δ≤ ,fr=2kHz ， At≥40dB, fs=8kHz，比较这种滤波器的阶数。
MATLAB源程序：clear all;wc=2*pi*1200;wr=2*pi*2000;rp=;rs=40;fs=8000;w1=2*fs*tan(wc/(2*fs));w2=2*fs*tan(wr/(2*fs));[Nb,wn]=buttord(w1,w2,rp,rs,'s') %巴特沃思[B,A]=butter(Nb,wn,'s');[num1,den1]=bilinear(B,A,fs);[h1,w]=freqz(num1,den1);[Nc,wn]=cheb1ord(w1,w2,rp,rs,'s') %切比雪夫[B,A]=cheby1(Nc,rp,wn,'s');[num2,den2]=bilinear(B,A,fs);[h2,w]=freqz(num2,den2);[Ne,wn]=ellipord(w1,w2,rp,rs,'s') %椭圆型[B,A]=ellip(Ne,rp,rs,wn,'low','s');[num3,den3]=bilinear(B,A,fs);[h3,w]=freqz(num3,den3);f=w/(2*pi)*fs;plot(f,20*log10(abs(h1)),'-',f,20*log10(abs(h2)),'--',f,20* log10(abs(h3)),':');axis([0,4000,-100,10]);grid;xlabel('Frequency in Hz'); ylabel('Gain in dB');title('三种数字低通滤波器');legend('巴特沃思数字低通滤波器','切比雪夫数字低通滤波器','椭圆数字低通滤波器',3);巴特沃思数字低通滤波器的系统函数系数：num1=den1=切比雪夫数字低通滤波器的系统函数系数：num2=den2= 1椭圆数字低通滤波器的系统函数系数：num3=den3= 1程序结果图：分析：设计结果表明，巴特沃思数字低通滤波器、切比雪夫数字低通滤波器、椭圆数字低通滤波器的阶数分别是9、5、4阶。
可见，对于给定的阶数，椭圆数字低通滤波器的阶数最少(换言之，对于给定的阶数，过渡带最窄)，就这一点来说，他是最优滤波器。
由图表明，巴特沃思数字低通滤波器过渡带最宽，幅频响应单调下降；椭圆数字低通滤波器过渡带最窄，并具有等波纹的通带和阻带响应；切比雪夫数字低通滤波器的过渡带介于两者之间。