IIR数字滤波器的设计实验报告

合集下载

实验五实验报告IIR数字滤波器设计与滤波

实验五实验报告IIR数字滤波器设计与滤波

实验五实验报告IIR数字滤波器设计与滤波实验目的:1. 掌握IIR(Infinite Impulse Response)数字滤波器的设计原理;2.熟悉IIR数字滤波器的各种滤波器类型及其特点;3.学会使用MATLAB进行IIR数字滤波器的设计和滤波。

实验仪器和设备:1.MATLAB软件;2.个人电脑。

实验步骤:1.确定滤波器的设计规格,包括滤波器类型(低通、高通、带通、带阻)、截止频率、通带和阻带的衰减要求等;2. 使用MATLAB中的`firpm`函数进行IIR数字滤波器的设计,指定滤波器的设计规格,得到滤波器的频率响应;3. 根据频率响应,使用`freqz`函数绘制滤波器的频率响应曲线,并分析其滤波特性;4. 利用`filter`函数将待滤波的信号输入到滤波器中,进行滤波操作;5.将滤波后的信号与原始信号进行对比分析,评估滤波器的性能。

实验结果:1.设计了一个低通IIR数字滤波器,截止频率为1kHz,通带和阻带的衰减要求分别为1dB和40dB;2.绘制出滤波器的频率响应曲线,发现在截止频率1kHz处有明显的衰减;3.将一段含有高频噪声的信号输入到滤波器中进行滤波,滤波后的信号去除了高频成分,保留了较低频的信号;4.与未滤波的信号进行对比,滤波后的信号显示出明显的去噪效果,且保留了原始信号的主要特征。

实验结论:通过本次实验,我们成功设计并实现了一个IIR数字滤波器,并对其滤波性能进行了评估。

实验结果表明,IIR数字滤波器能有效去除信号中的噪声和干扰成分,保留信号中的主要特征。

在实际应用中,我们可以根据需要选择不同类型和性能的IIR数字滤波器,以实现特定的信号处理任务。

IIR数字滤波器设计实验报告

IIR数字滤波器设计实验报告

实验三IIR数字滤波器设计实验报告一、实验目的:1.通过仿真冲激响应不变法和双线性变换法2.掌握滤波器性能分析的基本方法二、实验要求:1.设计带通IIR滤波器2.按照冲激响应不变法设计滤波器系数3. 按照双线性变换法设计滤波器系数4. 分析幅频特性和相频特性5. 生成一定信噪比的带噪信号,并对其滤波,对比滤波前后波形和频谱三、基本原理:㈠IIR模拟滤波器与数字滤波器IIR数字滤波器的设计以模拟滤波器设计为基础,常用的类型分为巴特沃斯(Butterworth)、切比雪夫(Chebyshev)Ⅰ型、切比雪夫Ⅱ型、贝塞尔(Bessel)、椭圆等多种。

在MATLAB信号处理工具箱里,提供了这些类型的IIR数字滤波器设计子函数。

(二)性能指标1.假设带通滤波器要求为保留6000hz~~7000hz频段,滤除小于2000hz和大宇9000hz频段2.通带衰减设为3Db,阻带衰减设为30dB,双线性变换法中T取1s.四、实验步骤:1.初始化指标参数2.计算模拟滤波器参数并调用巴特沃斯函数产生模拟滤波器3.利用冲激响应不变法和双线性变换法求数字IIR滤波器的系统函数Hd (z)4.分别画出两种方法的幅频特性和相频特性曲线5.生成一定信噪比的带噪信号6.画出带噪信号的时域图和频谱图6.对带噪信号进行滤波,并画出滤波前后波形图和频谱图五、实验结果模拟滤波器的幅频特性和相频特性:101010101Frequency (rad/s)P h a s e (d e g r e e s )1010101011010-5100Frequency (rad/s)M a g n i t u d e在本实验中,采用的带通滤波器为6000-7000Hz ,换算成角频率为4.47-0.55,在上图中可以清晰地看出到达了题目的要求。

冲击响应不变法后的幅频特性和相频特性:0.10.20.30.40.50.60.70.80.91Normalized Frequency (⨯π rad/sample)P h a s e (d e g r e e s )0.10.20.30.40.50.60.70.80.91Normalized Frequency (⨯π rad/sample)M a g n i t u d e (d B )双线性变换法的幅频特性和相频特性:0.10.20.30.40.50.60.70.80.91Normalized Frequency (⨯π rad/sample)P h a s e (d e g r e e s )00.10.20.30.40.50.60.70.80.91Normalized Frequency (⨯π rad/sample)M a g n i t u d e (d B )通过上图比较脉冲响应不变法双线性变换法的幅频特性和相频特性,而在在幅频曲线上几乎没有差别,都能达到相同的结果。

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数字滤波器的设计和分析。

在实验开始前,对于IIR数字滤波器有所了解,它是一种无限级别功能的数字滤波器,其功能强大,可以实现任意自定义系数的滤波器。

在预处理实验中,便首先采用Matlab
工具搭建了IIR数字滤波器的框架,考虑到本次滤波处理内容,本次采用的是Chebyshev
类型的等离子体,其滤波效果要求超过50dB,进一步完善了对于设计工作的要求。

经过Chebyshev Type I等离子体的设计,确定了系统的结构,并设定了15个滤波器,接着从设定的各项参量入手,从而确定系统各项参量,运用梯形图确定根位置,并使用MATLAB中的filter函数进行系统模拟,得到经历处理后系统输出信号与未经处理时对比,结果显示滤波效果达到了相应预期要求。

在实验中,IIR数字滤波器的设计让我深刻体会到了系统滤波的重要性以及十分强大
的功能。

而它的实现,又显示了精确的数字处理技术在信号处理中的重要作用,使得研究
信号处理时,得以有效和准确地对信号进行分辨和滤波处理。

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滤波器设计实验报告IIR滤波器设计实验报告一、实验目的本实验旨在通过实际操作,掌握IIR滤波器的基本设计方法,了解滤波器性能参数对滤波效果的影响,加深对滤波器理论的理解。

二、实验原理IIR滤波器(Infinite Impulse Response)是一种离散时间滤波器,其系统函数具有无限长的时间响应。

IIR滤波器设计方法主要包括冲激响应不变法和双线性变换法。

本实验采用冲激响应不变法进行设计。

三、实验步骤1. 确定滤波器性能参数:根据实际需求,确定滤波器的类型(低通、高通、带通、带阻)、通带边缘频率、阻带边缘频率、通带波动和阻带衰减等性能参数。

2. 计算滤波器系数:根据冲激响应不变法,利用Matlab编程计算滤波器系数。

具体过程包括定义系统函数、计算冲激响应、计算频率响应等步骤。

3. 编写滤波器程序:根据计算出的滤波器系数,编写IIR滤波器程序。

程序应实现输入信号的滤波处理,并输出滤波后的信号。

4. 测试滤波器性能:对编写的滤波器程序进行测试,观察其滤波效果,分析性能参数对滤波效果的影响。

5. 优化滤波器性能:根据测试结果,对滤波器性能参数进行调整,优化滤波效果。

四、实验结果及分析通过本次实验,我们成功地设计并实现了IIR滤波器。

在测试过程中,我们观察到了滤波器对不同频率信号的过滤效果,并分析了性能参数对滤波效果的影响。

具体来说,通带边缘频率决定了滤波器对低频信号的过滤程度,阻带边缘频率则影响对高频信号的过滤程度。

通带波动和阻带衰减则分别反映了滤波器在通带和阻带的波动程度和衰减程度。

通过对这些性能参数的调整,我们可以实现对不同类型信号的有效过滤。

五、实验总结通过本次实验,我们深入理解了IIR滤波器的工作原理和设计方法,掌握了Matlab编程在滤波器设计中的应用。

实验过程中,我们不仅学会了如何根据实际需求选择合适的性能参数,还学会了如何调整这些参数以优化滤波效果。

此外,我们还观察到了不同性能参数对滤波效果的影响,加深了对滤波器理论的理解。

IIR数字滤波器设计及软件实现实验报告

IIR数字滤波器设计及软件实现实验报告

实验报告实验四:IIR数字滤‎波器设计及软‎件实现1.实验目的(1)熟悉用双线性‎变换法设计I‎I R数字滤波‎器的原理与方‎法;(2)学会调用MA‎T LAB信号‎处理工具箱中‎滤波器设计函‎数(或滤波器设计‎分析工具fd‎a tool)设计各种II‎R数字滤波器‎,学会根据滤波‎需求确定滤波‎器指标参数。

(3)掌握IIR数‎字滤波器的M‎A TLAB实‎现方法。

(3)通过观察滤波‎器输入输出信‎号的时域波形‎及其频谱,建立数字滤波‎的概念。

2.实验原理设计IIR数‎字滤波器一般‎采用间接法(脉冲响应不变‎法和双线性变‎换法),应用最广泛的‎是双线性变换‎法。

基本设计过程‎是:①先将给定的数‎字滤波器的指‎标转换成过渡‎模拟滤波器的‎指标;②设计过渡模拟‎滤波器;③将过渡模拟滤‎波器系统函数‎转换成数字滤‎波器的系统函‎数。

MATLAB‎信号处理工具‎箱中的各种I‎I R数字滤波‎器设计函数都‎是采用双线性‎变换法。

第六章介绍的‎滤波器设计函‎数butte‎r、cheby1‎、cheby2‎和ellip‎可以分别被调‎用来直接设计‎巴特沃斯、切比雪夫1、切比雪夫2和‎椭圆模拟和数‎字滤波器。

本实验要求读‎者调用如上函‎数直接设计I‎I R数字滤波‎器。

本实验的数字‎滤波器的MA‎T LAB实现‎是指调用MA‎T LAB信号‎处理工具箱函‎数filte‎r对给定的输‎入信号x(n)进行滤波,得到滤波后的‎输出信号y(n)。

3. 实验内容及步‎骤(1)调用信号产生‎函数mstg‎产生由三路抑‎制载波调幅信‎号相加构成的‎复合信号st‎,该函数还会自‎动绘图显示s‎t的时域波形‎和幅频特性曲‎线,如图10.4.1所示。

由图可见,三路信号时域‎混叠无法在时‎域分离。

但频域是分离‎的,所以可以通过‎滤波的方法在‎频域分离,这就是本实验‎的目的。

图10.4.1 三路调幅信号‎s t的时域波‎形和幅频特性‎曲线(2)要求将st中‎三路调幅信号‎分离,通过观察st‎的幅频特性曲‎线,分别确定可以‎分离st中三‎路抑制载波单‎频调幅信号的‎三个滤波器(低通滤波器、带通滤波器、高通滤波器)的通带截止频‎率和阻带截止‎频率。

matlab实验报告 IIR数字滤波器设计

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)实现模型转换。

  1. 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
  2. 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
  3. 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。

学生实验报告
实验内容和步
骤(原理、主要步骤、算法、程序、运行结果、对结果的讨论,思考题解答等)1. 通过模拟滤波器原型设计一个butterworth模拟低通滤波器,要求通带截止频率fp=2kHz, 通带最大衰减Rp=1dB, 阻带截止频率fp=2kHz, 阻带最小衰减As=20dB,
%巴特沃斯模拟滤波器
fp=2000;Omgp=2*pi*fp; %输入滤波器的通带截止频率
fs=5000;Omgs=2*pi*fs; %输入滤波器的阻带截止频率
Rp=1;As=20; %输入滤波器的通阻带衰减指标
[n,Omgc]=buttord(Omgp,Omgs,Rp,As,'s')
%计算n阶模拟低通原型,得到左半平面零极点
[z0,p0,k0]=buttap(n); %由滤波器阶数N求模拟滤波器原型
b0=k0*real(poly(z0)) %求滤波器系数b0
a0=real(poly(p0)) %求滤波器系数a0
[H,Omg]=freqs(b0,a0); %求系统的频率特性
subplot(2,1,1); plot(Omg*Omgc/(2*pi),abs(H)),grid
axis([0,6000,0,1.1])
subplot(2,1,2); plot(Omg*Omgc/(2*pi),angle(H)),grid
axis([0,6000,-4,4])
运行结果:
思考题:1. 结合基本原理理解每一条语句的意义
2. buttord命令实现了什么功能? buttap命令实现了什么功能?
3. 所得出的系统的频率特性是否满足了所要求的性能指标?
set(gca,'YTickMode','manual','YTick',[-50,-15,-1,0]);grid
subplot(2,2,4),zplane(bd,ad);
axis([-1.1,1.1,-1.1,1.1]);title('零极图');
运行结果:
思考题:1. 结合基本原理理解每一条语句的意义
2. 所得出的系统的频率特性是否满足了所要求的性能指标?
3.根据零极点图判断系统是否为稳定系统?
答:2.满足要求性能指标;
3.是稳定系统。

心得和建议本次实验相对较简单,通过实验分析了解了巴特沃斯模拟滤波器的频率特性,同时学会了用双线性变换法设计IIR滤波器。

注:学生做每个实验都必须填写实验报告。

实验报告可以是纸质的,也可以是电子形式的。

凡电子文本形式的实验报告须在学期末刻录进该学期教学资料光盘并存个人教学档案。

实验指导教师签名年月日。

相关文档
最新文档