基于Matla的数字滤波.pptVIP

基于Matlab的数字滤波器设计 基于Matlab的数字滤波器设计 1.数字滤波器概述 1.1数字滤波器的分类 1.2 数字滤波器的表示方法 2.数字滤波器的优点 3 Matlab中的数字滤波器设计 3.1软件实现方式 3.2硬件实现方式 1数字滤波器概述 数字滤波器不用物理元器件构成，而是按照某种算法编写一段程序或模块，对数字信号进行处理加工，从而达到滤波的要求。数字滤波器是实现滤波过程的一种数字信号处理系统，具有离散时间系统的基本特征。 1.1数字滤波器的分类 按频率特性分 按冲击响应分 按结构特点分 1.2 数字滤波器的表示方法 1.差分方式表示 2.传递函数表示 3.单位响应冲击响应 输入单位冲击时间序列的滤波器输出序列 4.频率特性表示 频率特性是数字滤波器对正弦输入序列的响应 2 数字滤波器的优点 3 Matlab中的数字滤波器设计 有两种方法： 1基于滤波程序和算法的软件实现方式。 2基于数字部件模块的硬件实现方式。 软件实现方式 以简单的减法滤波器作为介绍： 减法滤波器差分方程为 根据想消除的谐波次数，来确定参数的值。可以根据下式确定 例 假设一个系统的采样频率为1200Hz，要滤除直流分量和4,8,12次谐波 。 分析：由于奈奎斯特采样定理的限制，最多只能消除12次谐波 ，当采用减法滤波器时，经过计算所需差分方程为 在Matlab的M脚本下编写滤波函数，绘出幅频特性。 %设置减法滤波器的传递函数系数 a1=1,b1=[1 0 0 0 0 0 -1]; f=0:1:600; h1=abs(freqz(b1,a1,f,1200)); %由传递函数系数确定传递函数的幅频特性 H1=h1/max(h1); %绘出幅频特性 plot(f,H1); xlabel(f/Hz);ylabel(H1); 叠加直流分量和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); %采用减法滤波器滤掉va1中的直流分量和4,8次谐波分量 Y=zeros(1,6); for jj=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); grid on 硬件实现方式 在Matlab中，有现成的滤波器模块，使用时只需设置相应参数，使用方便。 在Matlab中Simulink命令下Single Processing Blockset\Filtering\Filter Designs目录下有Digital Filter Design模块可以设计用户所需要的数字滤波器。 为了更好的说明各种设计方法之间的差别，采用控制变量的方法，在其他条件不变的前提下，分别采用窗函数法（凯泽窗、汉宁窗、汉明窗、布莱克窗窗）、波纹逼近法、最小乘方法对FIR数字低通滤波器进行设计并比较，最终得出结论，说明各种设计方法的适用情况。 97阶凯泽窗 97阶汉明窗 97阶汉宁窗 97阶布莱克窗 波纹逼近法 最小乘方法 结论 在同样滤波器阶数情况下，窗函数设计法设计方法相对简便，通带内稳定性高，过渡带宽度随窗函数种类变化，其中汉宁窗、汉明窗、布莱克曼窗函数法的过渡带相比凯泽窗较窄，但凯泽窗最小阻带衰减最高；波纹逼近法设计出的滤波器相对其他几种方法通带内不够平稳；最小乘方法相比于其他设计方法误差较低，通带内平稳度较高并且最小阻带衰减较高。通过分析推广到实际设计中去：对通带内平稳度要求较高时优先选用窗函数法；对最小阻带衰减要求较高时可优先选择最小乘方法；要求阶数或计算量较低时可优先选择波纹逼近法。 示例 在小电流接地系统中注入83.3Hz的正弦信号，将其滤除。 分析：采用波纹逼近法，采样频率1000Hz。 假设输入电流为 谢谢！ * * * * 按结构特点分 按冲击响应分 按频率特性分 数字滤波器 数字滤波器 低通滤波器 高通滤波器 带通滤波器 带阻滤波器 数字滤波器 无限长冲击响应滤波器(IIR)： 对