基于Matlab的数字滤波课件

基于Matlab的数字滤波器设计

基于Matlab的数字滤波器设计1基于Matlab的数字滤波器设计1.数字滤波器概述1.1数字滤波器的分类1.2数字滤波器的表示方法2.数字滤波器的优点3Matlab中的数字滤波器设计3.1软件实现方式3.2硬件实现方式基于Matlab的数字滤波器设计1.数字滤波器概述21数字滤波器概述数字滤波器不用物理元器件构成，而是按照某种算法编写一段程序或模块，对数字信号进行处理加工，从而达到滤波的要求。数字滤波器是实现滤波过程的一种数字信号处理系统，具有离散时间系统的基本特征。1数字滤波器概述数字滤波器不用物理元器件构成，而是按照某31.1数字滤波器的分类按结构特点分

按冲击响应分

按频率特性分

数字滤波器

1.1数字滤波器的分类按结构特点分按冲击响应分按频率特性4按频率特性分数字滤波器

低通滤波器高通滤波器带通滤波器带阻滤波器

按频率特性分数字滤波器低通滤波器高通滤波器带通滤波器带阻滤5按冲击响应分数字滤波器无限长冲击响应滤波器(IIR)：对单位冲激的输入信号的响应为无限长序列

有限长冲击响应滤波器(FIR)：对单位冲激的输入信号的响应为有限长序列按冲击响应分数字滤波器无限长冲击响应滤波器(IIR)：有限长6按结构特点分数字滤波器数字滤波器递归型滤波器

非递归型滤波器无限长冲击响应滤波器有限长冲击响应滤波器按结构特点分数字滤波器数字滤波器递归型滤波器非递归型滤波71.2数字滤波器的表示方法1.差分方式表示

其中输入信号为x(t),输出信号为y(t),ai,bi为常数1.2数字滤波器的表示方法1.差分方式表示其中输入信号为82.传递函数表示在差分方程的基础上，方程取Z变换，化简可得2.传递函数表示在差分方程的基础上，方程取Z变换，化简可得93.单位响应冲击响应输入单位冲击时间序列的滤波器输出序列任意输入信号时的数字滤波器输出3.单位响应冲击响应任意输入信号时的数字滤波器输出104.频率特性表示频率特性是数字滤波器对正弦输入序列的响应

其中，H(jwTs)和频率f构成幅频特性，ψ(w)和频率f构成相频特性。4.频率特性表示其中，H(jwTs)和频率f构成幅频特性，ψ112数字滤波器的优点性能稳定

高度的灵活性

无阻抗匹配问题

方便的做到分时复用

2数字滤波器的优点性能稳定高度的灵活性无阻抗匹配问题123Matlab中的数字滤波器设计

有两种方法：1基于滤波程序和算法的软件实现方式。

2基于数字部件模块的硬件实现方式。适用于知道需要滤除的谐波的次数适用于硬件搭建模型3Matlab中的数字滤波器设计有两种方法：适用于13软件实现方式以简单的减法滤波器作为介绍：减法滤波器差分方程为

其幅频特性方程为软件实现方式以简单的减法滤波器作为介绍：其幅频特性方程为14根据想消除的谐波次数，来确定参数的值。可以根据下式确定式中，为1,2,3……。这种方法虽然可以达到滤波效果，但是需要人为计算，在编写滤波程序之前，需要知道滤波器的差分方程，较为麻烦。根据想消除的谐波次数，来确定参数的值。可以根据下式确定15例假设一个系统的采样频率为1200Hz，要滤除直流分量和4,8,12次谐波。

分析：由于奈奎斯特采样定理的限制，最多只能消除12次谐波，当采用减法滤波器时，经过计算所需差分方程为例假设一个系统的采样频率为1200Hz，要滤除直流分量16在Matlab的M脚本下编写滤波函数，绘出幅频特性。%设置减法滤波器的传递函数系数a1=1,b1=[100000-1];f=0:1:600;h1=abs(freqz(b1,a1,f,1200));%由传递函数系数确定传递函数的幅频特性H1=h1/max(h1);%绘出幅频特性plot(f,H1);xlabel('f/Hz');ylabel('H1');在Matlab的M脚本下编写滤波函数，绘出幅频特性。17幅频特性幅频特性18叠加直流分量和4、8次谐波，然后滤波N=24;t1=(0:0.02/N:0.04);m=size(t1);%基波电压va=100*sin(2*pi*50*t1);%叠加直流分量和4,8次谐波分量va1=35+100*sin(2*pi*50*t1)+30*sin(4*pi*100*t1)+10*sin(8*pi*100*t1);叠加直流分量和4、8次谐波，然后滤波19%采用减法滤波器滤掉va1中的直流分量和4,8次谐波分量Y=zeros(1,6);forjj=7:m(2)Y(jj)=(va1(jj)-va1(jj-6))/1.414;end%输出波形,每周期采样24个点，一个点相当于15度，6个点为90度，移项90度的波形和原始波形叠加后波形幅值会增大倍，为了便于比较，叠加后的波形幅值除以。plot(t1,va,'-ro',t1,va1,'-bs',t1,Y,'-g*');xlabel('t/s');ylabel('v/v');gridon%采用减法滤波器滤掉va1中的直流分量和4,8次谐波分量20基于Matlab的数字滤波课件21硬件实现方式在Matlab中，有现成的滤波器模块，使用时只需设置相应参数，使用方便。在Matlab中Simulink命令下SingleProcessingBlockset\Filtering\FilterDesigns目录下有DigitalFilterDesign模块可以设计用户所需要的数字滤波器。硬件实现方式在Matlab中，有现成的滤波器模块，使用时只22参数设计区特性区滤波器转换（TransFormFilter）

实现模型（RealizeModel）

零极点编辑器（Pole-zeroEditor）

导入滤波器（ImportFilter）

设计滤波器（DesignFilter）

参数设计区特性区滤波器转换（TransFormFilter23Responsetype（滤波器类型）选项，包括lowpass（低通）、highpass（高通）、bandpass（带通）、bandstop（带阻）和特殊的FIR滤波器Designmethod（设计方法）选项，包括IIR滤波器的butterworth（巴特沃思）法、chebyshevtypei（切比雪夫i型）法、chebyshevtypeii（切比雪夫ii型）法、elliptic（椭圆滤波器）法和FIR滤波器的equiripple法、least-squares（最小乘方）法、window（窗函数）法Responsetype（滤波器类型）选项，包括lowpa24Filterorder（滤波器阶数）选项，定义滤波器的阶数，包括specifyorder（指定阶数）和minimumorder（最小阶数）。在specifyorder中填入所要设计的滤波器的阶数（n阶滤波器，specifyorder＝n-1），如果选择minimumorder则matlab根据所选择的滤波器类型自动使用最小阶数。Filterorder（滤波器阶数）选项，定义滤波器的阶数25Frenquencyspecifications选项，可以详细定义频带的各参数，包括采样频率fs和频带的截止频率。它的具体选项由filtertype选项和designmethod选项决定，bandpass（带通）滤波器需要定义fstop1（下阻带截止频率）、fpass1（通带下限截止频率）、fpass2（通带上限截止频率）、fstop2（上阻带截止频率），而lowpass（低通）滤波器只需要定义fstop1、fpass1。采用窗函数设计滤波器时，由于过渡带是由窗函数的类型和阶数所决定的，所以只需要定义通带截止频率，而不必定义阻带参数。Frenquencyspecifications选项，可以26Magnitudespecifications选项，可以定义幅值衰减的情况。例如设计带通滤波器时，可以定义wstop1（频率fstop1处的幅值衰减）、wpass（通带范围内的幅值衰减）、wstop2（频率fstop2处的幅值衰减）。当采用窗函数设计时，通带截止频率处的幅值衰减固定为6db，所以不必定义。Magnitudespecifications选项，可以定27为了更好的说明各种设计方法之间的差别，采用控制变量的方法，在其他条件不变的前提下，分别采用窗函数法（凯泽窗、汉宁窗、汉明窗、布莱克窗窗）、波纹逼近法、最小乘方法对FIR数字低通滤波器进行设计并比较，最终得出结论，说明各种设计方法的适用情况。为了更好的说明各种设计方法之间的差别，采用控制变量的方2897阶凯泽窗97阶凯泽窗2997阶汉明窗97阶汉明窗3097阶汉宁窗97阶汉宁窗3197阶布莱克窗97阶布莱克窗32波纹逼近法波纹逼近法33最小乘方法最小乘方法34结论在同样滤波器阶数情况下，窗函数设计法设计方法相对简便，通带内稳定性高，过渡带宽度随窗函数种类变化，其中汉宁窗、汉明窗、布莱克曼窗函数法的过渡带相比凯泽窗较窄，但凯泽窗最小阻带衰减最高；波纹逼近法设计出的滤波器相对其他几种方法通带内不够平稳；最小乘方法相比于其他设计方法误差较低，通带内平稳度较高并且最小阻带衰减较高。通过分析推广到实际设计中去：对通带内平稳度要求较高时优先选用窗函数法；对最小阻带衰减要求较高时可优先选择最小乘方法；要求阶数或计算量较低时可优先选择波纹逼近法。结论在同样滤波器阶数情况下，窗函数设计法设计方法相对简35示例在小电流接地系统中注入83.3Hz的正弦信号，将其滤除。分析：采用波纹逼近法，采样频率1000Hz。假设输入电流为示例36基于Matlab的数字滤波课件37基于Matlab的数字滤波课件38基于Matlab的数字滤波课件39谢谢！谢谢！40基于Matlab的数字滤波器设计

基于Matlab的数字滤波器设计41基于Matlab的数字滤波器设计1.数字滤波器概述1.1数字滤波器的分类1.2数字滤波器的表示方法2.数字滤波器的优点3Matlab中的数字滤波器设计3.1软件实现方式3.2硬件实现方式基于Matlab的数字滤波器设计1.数字滤波器概述421数字滤波器概述数字滤波器不用物理元器件构成，而是按照某种算法编写一段程序或模块，对数字信号进行处理加工，从而达到滤波的要求。数字滤波器是实现滤波过程的一种数字信号处理系统，具有离散时间系统的基本特征。1数字滤波器概述数字滤波器不用物理元器件构成，而是按照某431.1数字滤波器的分类按结构特点分

按冲击响应分

按频率特性分

数字滤波器

1.1数字滤波器的分类按结构特点分按冲击响应分按频率特性44按频率特性分数字滤波器

低通滤波器高通滤波器带通滤波器带阻滤波器

按频率特性分数字滤波器低通滤波器高通滤波器带通滤波器带阻滤45按冲击响应分数字滤波器无限长冲击响应滤波器(IIR)：对单位冲激的输入信号的响应为无限长序列

有限长冲击响应滤波器(FIR)：对单位冲激的输入信号的响应为有限长序列按冲击响应分数字滤波器无限长冲击响应滤波器(IIR)：有限长46按结构特点分数字滤波器数字滤波器递归型滤波器

非递归型滤波器无限长冲击响应滤波器有限长冲击响应滤波器按结构特点分数字滤波器数字滤波器递归型滤波器非递归型滤波471.2数字滤波器的表示方法1.差分方式表示

其中输入信号为x(t),输出信号为y(t),ai,bi为常数1.2数字滤波器的表示方法1.差分方式表示其中输入信号为482.传递函数表示在差分方程的基础上，方程取Z变换，化简可得2.传递函数表示在差分方程的基础上，方程取Z变换，化简可得493.单位响应冲击响应输入单位冲击时间序列的滤波器输出序列任意输入信号时的数字滤波器输出3.单位响应冲击响应任意输入信号时的数字滤波器输出504.频率特性表示频率特性是数字滤波器对正弦输入序列的响应

其中，H(jwTs)和频率f构成幅频特性，ψ(w)和频率f构成相频特性。4.频率特性表示其中，H(jwTs)和频率f构成幅频特性，ψ512数字滤波器的优点性能稳定

高度的灵活性

无阻抗匹配问题

方便的做到分时复用

2数字滤波器的优点性能稳定高度的灵活性无阻抗匹配问题523Matlab中的数字滤波器设计

有两种方法：1基于滤波程序和算法的软件实现方式。

2基于数字部件模块的硬件实现方式。适用于知道需要滤除的谐波的次数适用于硬件搭建模型3Matlab中的数字滤波器设计有两种方法：适用于53软件实现方式以简单的减法滤波器作为介绍：减法滤波器差分方程为

其幅频特性方程为软件实现方式以简单的减法滤波器作为介绍：其幅频特性方程为54根据想消除的谐波次数，来确定参数的值。可以根据下式确定式中，为1,2,3……。这种方法虽然可以达到滤波效果，但是需要人为计算，在编写滤波程序之前，需要知道滤波器的差分方程，较为麻烦。根据想消除的谐波次数，来确定参数的值。可以根据下式确定55例假设一个系统的采样频率为1200Hz，要滤除直流分量和4,8,12次谐波。

分析：由于奈奎斯特采样定理的限制，最多只能消除12次谐波，当采用减法滤波器时，经过计算所需差分方程为例假设一个系统的采样频率为1200Hz，要滤除直流分量56在Matlab的M脚本下编写滤波函数，绘出幅频特性。%设置减法滤波器的传递函数系数a1=1,b1=[100000-1];f=0:1:600;h1=abs(freqz(b1,a1,f,1200));%由传递函数系数确定传递函数的幅频特性H1=h1/max(h1);%绘出幅频特性plot(f,H1);xlabel('f/Hz');ylabel('H1');在Matlab的M脚本下编写滤波函数，绘出幅频特性。57幅频特性幅频特性58叠加直流分量和4、8次谐波，然后滤波N=24;t1=(0:0.02/N:0.04);m=size(t1);%基波电压va=100*sin(2*pi*50*t1);%叠加直流分量和4,8次谐波分量va1=35+100*sin(2*pi*50*t1)+30*sin(4*pi*100*t1)+10*sin(8*pi*100*t1);叠加直流分量和4、8次谐波，然后滤波59%采用减法滤波器滤掉va1中的直流分量和4,8次谐波分量Y=zeros(1,6);forjj=7:m(2)Y(jj)=(va1(jj)-va1(jj-6))/1.414;end%输出波形,每周期采样24个点，一个点相当于15度，6个点为90度，移项90度的波形和原始波形叠加后波形幅值会增大倍，为了便于比较，叠加后的波形幅值除以。plot(t1,va,'-ro',t1,va1,'-bs',t1,Y,'-g*');xlabel('t/s');ylabel('v/v');gridon%采用减法滤波器滤掉va1中的直流分量和4,8次谐波分量60基于Matlab的数字滤波课件61硬件实现方式在Matlab中，有现成的滤波器模块，使用时只需设置相应参数，使用方便。在Matlab中Simulink命令下SingleProcessingBlockset\Filtering\FilterDesigns目录下有DigitalFilterDesign模块可以设计用户所需要的数字滤波器。硬件实现方式在Matlab中，有现成的滤波器模块，使用时只62参数设计区特性区滤波器转换（TransFormFilter）

实现模型（RealizeModel）

零极点编辑器（Pole-zeroEditor）

导入滤波器（ImportFilter）

设计滤波器（DesignFilter）

参数设计区特性区滤波器转换（TransFormFilter63Responsetype（滤波器类型）选项，包括lowpass（低通）、highpass（高通）、bandpass（带通）、bandstop（带阻）和特殊的FIR滤波器Designmethod（设计方法）选项，包括IIR滤波器的butterworth（巴特沃思）法、chebyshevtypei（切比雪夫i型）法、chebyshevtypeii（切比雪夫ii型）法、elliptic（椭圆滤波器）法和FIR滤波器的equiripple法、least-squares（最小乘方）法、window（窗函数）法Responsetype（滤波器类型）选项，包括lowpa64Filterorder（滤波器阶数）选项，定义滤波器的阶数，包括specifyorder（指定阶数）和minimumorder（最小阶数）。在specifyorder中填入所要设计的滤波器的阶数（n阶滤波器，specifyorder＝n-1），如果选择minimumorder则matlab根据所选择的滤波器类型自动使用最小阶数。Filterorder（滤波器阶数）选项，定义滤波器的阶数65Frenquencyspecifications选项，可以详细定义频带的各参数，包括采样频率fs和频带的截止频率。它的具体选项由filtertype选项和designmethod选项决定，bandpass（带通）滤波器需要定义fstop1（下阻带截止频率）、fpass1（通带下限截止频率）、fpass2（通带上限截止频率）、fstop2（上阻带截止频率），而lowpass（低通）滤波器只需要定义fstop1、fpass1。采用窗函数设计滤波器时，由于过渡带是由窗函数的类型和阶数所决定的，所以只需要定义通带截止频率，而不必定义阻带参数。Frenquencyspec