基于Matlab的数字滤波.pptVIP

Magnitudespecifications选项，可以定义幅值衰减的情况。例如设计带通滤波器时，可以定义wstop1（频率fstop1处的幅值衰减）、wpass（通带范围内的幅值衰减）、wstop2（频率fstop2处的幅值衰减）。当采用窗函数设计时，通带截止频率处的幅值衰减固定为6db，所以不必定义。为了更好的说明各种设计方法之间的差别，采用控制变量的方法，在其他条件不变的前提下，分别采用窗函数法（凯泽窗、汉宁窗、汉明窗、布莱克窗窗）、波纹逼近法、最小乘方法对FIR数字低通滤波器进行设计并比较，最终得出结论，说明各种设计方法的适用情况。97阶凯泽窗基于Matlab的数字滤波器设计数字滤波器概述1数字滤波器的分类1.2数字滤波器的表示方法数字滤波器的优点Matlab中的数字滤波器设计1软件实现方式2硬件实现方式基于Matlab的数字滤波器设计1数字滤波器概述数字滤波器不用物理元器件构成，而是按照某种算法编写一段程序或模块，对数字信号进行处理加工，从而达到滤波的要求。数字滤波器是实现滤波过程的一种数字信号处理系统，具有离散时间系统的基本特征。1.1数字滤波器的分类数字滤波器按冲击响应分按结构特点分按频率特性分按频率特性分数字滤波器高通滤波器带阻滤波器低通滤波器带通滤波器数字滤波器无限长冲击响应滤波器(IIR)：对单位冲激的输入信号的响应为无限长序列有限长冲击响应滤波器(FIR)：对单位冲激的输入信号的响应为有限长序列按冲击响应分按结构特点分数字滤波器数字滤波器递归型滤波器非递归型滤波器无限长冲击响应滤波器有限长冲击响应滤波器1.2数字滤波器的表示方法logo1.差分方式表示其中输入信号为x(t),输出信号为y(t),ai,bi为常数传递函数表示在差分方程的基础上，方程取Z变换，化简可得0102单位响应冲击响应输入单位冲击时间序列的滤波器输出序列01任意输入信号时的数字滤波器输出02其中，H(jwTs)和频率f构成幅频特性，ψ(w)和频率f构成相频特性。02频率特性是数字滤波器对正弦输入序列的响应01频率特性表示2数字滤波器的优点方便的做到分时复用高度的灵活性性能稳定无阻抗匹配问题有两种方法：1基于滤波程序和算法的软件实现方式。2基于数字部件模块的硬件实现方式。适用于知道需要滤除的谐波的次数适用于硬件搭建模型0103023Matlab中的数字滤波器设计01以简单的减法滤波器作为介绍：减法滤波器差分方程为02其幅频特性方程为软件实现方式根据想消除的谐波次数，来确定参数的值。可以根据下式确定式中，为1,2,3……。这种方法虽然可以达到滤波效果，但是需要人为计算，在编写滤波程序之前，需要知道滤波器的差分方程，较为麻烦。例假设一个系统的采样频率为1200Hz，要滤除直流分量和4,8,12次谐波。分析：由于奈奎斯特采样定理的限制，最多只能消除12次谐波，当采用减法滤波器时，经过计算所需差分方程为壹贰在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);幅频特性叠加直流分量和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);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在Matlab中，有现成的滤波器模块，使用时只需设置相应参数，使用方