基于Matlab的数字滤波器算法的设计论文.docVIP

基于MATLAB的数字滤波器设计及滤波仿真 1.实验原理: ①正弦信号： ②巴特沃斯模拟滤波器的设计：[N,Wc]=buttord(Wp,Ws,ap,as,s) [N,Wc]=buttord(wp，ws，αp，αs) 用于计算巴特沃斯数字滤波器的阶数N和3dB截止频率wc。调用参数wp，ws分别为数字滤波器的通带、阻带截止频率的归一化值，要求：0≤wp≤1，0≤ws≤1。1表示数字频率pi。 αp，αs分别为通带最大衰减和组带最小衰减(dB)。 当ws≤wp时，为高通滤波器; 当wp和ws为二元矢量时，为带通或带阻滤波器，这时wc也是二元向量。N，wc作为butter函数的调用参数。 坛[N,Ωc]=buttord(Ωp，Ωs，αp，αs，‘s’) 用于计算巴特沃斯模拟滤波器的阶数N和3dB截止频率Ωc。 Ωp，Ωs，Ωc均为实际模拟角频率 [BS,AS]=butter(N,Wc,s) [B,A] = BUTTER(N,Wn)，设计一个阶数为n，频率为Wn的低通滤波器; [B,A] = BUTTER(N,Wn,ftype)可以设计高通，带阻滤波器，其中ftype参数的形式可以指定何种滤波器， ftype为‘high’时，设计一个阶数为n，频率为Wn的高通滤波器;ftype为‘stop’时，得到滤波器阶数 为2*n，频率范围为Wn = [W1,W2]的带阻滤波器; ③双线性变换设计IIR滤波器： [BZ,AZ]=bilinear(BS,AS,1/T) ④信号滤波 Y =filter(B,A,X) 输入X为滤波前序列，Y为滤波结果序列，B/A 提供滤波器系数，B为分子， A为分 母整个滤波过程是通过下面差分方程实现的： a(1)*y(n) = b(1)*x(n) + b(2)*x(n-1) + ... + b(nb+1)*x(n-nb) - a(2)*y(n-1) - ... - a(na+1)*y(n-na) [Y,Zf] = FILTER(B,A,X,Zi)，输入X为滤波前序列，Y为滤波结果序列，B/A 提供滤波器系数，B为分子， A为分母，并输入Zi指定X的初始状态，Zf为最终状态矢量2.Simulink的使用方法： 模块库中的模块可以直接用鼠标进行拖曳（选中模块，按住鼠标左键不放）而放到模型窗口中进行处理。 在模型窗口中，选中模块，则其4个角会出现黑色标记。此时可以对模块进行以下的基本操作。 移动：选中模块，按住鼠标左键将其拖曳到所需的位置即可。若要脱离线而移动，可按住shift键，再进行拖曳。 复制：选中模块，然后按住鼠标右键进行拖曳即可复制同样的一个功能模块。 删除：选中模块，按Delete键即可。若要删除多个模块，可以同时按住Shift键，再用鼠标选中多个模块，按Delete键即可。也可以用鼠标选取某区域，再按Delete键就可以把该区域中的所有模块和线等全部删除。 转向：为了能够顺序连接功能模块的输入和输出端，功能模块有时需要转向。在菜单Format中选择Flip Block旋转180度，选择Rotate Block顺时针旋转90度。或者直接按Ctrl+F键执行Flip Block，按Ctrl+R键执行Rotate Block。 改变大小：选中模块，对模块出现的4个黑色标记进行拖曳即可。 模块命名：先用鼠标在需要更改的名称上单击一下，然后直接更改即可。名称在功能模块上的位置也可以变换180度，可以用Format菜单中的Flip Name来实现，也可以直接通过鼠标进行拖曳。Hide Name可以隐藏模块名称。 颜色设定： Format菜单中的Foreground Color可以改变模块的前景颜色，Background Color可以改变模块的背景颜色；而模型窗口的颜色可以通过Screen Color来改变。 参数设定：用鼠标双击模块，就可以进入模块的参数设定窗口，从而对模块进行参数设定。参数设定窗口包含了该模块的基本功能帮助，为获得更详尽的帮助，可以点击其上的help按钮。通过对模块的参数设定， 就可以获得需要的功能模块。 属性设定：选中模块，打开Edit菜单的Block Properties可以对模块进行属性设定。包括Description属性、 Priority优先级属性、Tag属性、Open function属性、Attributes format string属性。其中Open function属性是一个很有用的属性，通过它指定一个函数名，则当该模块被双击之后，Simulink就会调用该函数执行，这种函数在MATLAB中称为回调函数。 模块的输入输出信号：模块处理的信号包括标量信号和向量信号；标量信号是一种单