IIR和FIR设置低通滤波器.docx

、题目 设原始信号为 x(n) cos( n ) 0.01* sin(n)，n 1,2, ,1000 100 采用 IIR 滤波器和 FIR滤波器设计低通滤波器，比较两类滤波器 的滤波结果。 要求：采用 MATLAB语言实现，并分析结果。 、设计过程 首先对原始信号进行频谱分析， 确定滤波器参数。 通过分析加噪 信号的频谱，噪声信号为cos( n ) cos( n ) 0.01*sin(n) 信号的频谱，噪声信号为 100 ， 原始信号为 100 。 所以确定参数 ωp=0.1π，ωs=0.25 π；αp=0.08 ，αs=8 ；采样间隔 T=8s 1、IIR滤波器设计 按照数字滤波器技术指标(通带边界频率 Wp 通带最大衰减 阻 带截止频率 阻带最小衰减)要求设计一个过渡模拟低通滤波器 Ha(s)，再使用脉冲响应不变法或双线性变换法将 Ha(s)转换成数字低 通滤波器的系统函数 H(z)。本题采用双线性变换法， 具体转换关系为： 2 wp Ωp tan p p T 2 Ωs 2tan ws s T 2 由参数设计巴特沃斯滤波器，根据 [N,wc]=buttord(wp,ws,Rp,As, s)计算模拟滤波器的阶数及截止频率， 由 [B,A]=butter(N,wc,Rp,s) 计算模拟滤波器系统函数的分子和分母多项式的系数向量 B 和 A。 2、FIR滤波器设计 （ 1）根据对阻带衰减及过渡带的指标要求， 选择窗函数的类型， 并估计窗口的长度 N。 （ 2）构造希望逼近的频率响应函数 Hd（ e^jw ） （ 3）计算 hd（n） （ 4）加窗得到设计结果： h（n） =hd（n）*w （n） 采用 [n,Wn,beta,type]=kaiserord（f,a,dev）函数来估计滤波器的阶 数，根据阻带衰减及过渡带的指标要求，选择窗函数的类型（本题使 用方法二及凯塞窗滤波器） 。这里 f 对应的频率， a=[1,0]为 f 指定的 各个频带上的幅值向量，一般只有 0 和 1 表示。 Devs用于指定各个 频带输出滤波器的频率响应与其期望幅值之间的最大输出误差或者 偏差 三、结果图 图1 ：无噪信号与原始信号 图 2 ：原始信号频谱与加噪信号频谱 图 3：原始信号与 IIR 滤波器的频谱对比图图 4：原始信号与 FIR 滤波器的频谱对比图图 5 图 3：原始信号与 IIR 滤波器的频谱对比图 图 4：原始信号与 FIR 滤波器的频谱对比图 图 5：原始信号与 IIR、 FIR滤波后信号对比图 结果分析： 图 6：加噪信号与 IIR、FIR 滤波后频谱对比图 A.相同技术指标下滤波器阶数 IIR 滤波器阶数： 6； FIR滤波器阶数： 14。 由此可以得到， 为了达到相同的技术指标，即通带截止频率、阻 带起始频率、通带内衰减、阻带内衰减， FIR滤波器所需的阶数更高 因为 FIR滤波器传输函数的极点固定在原点，所以只能用较高的阶数 达到高的选择性；对于同样的滤波器设计指标， FIR 滤波器所要求的 阶数可以比 IIR滤波器高 5~10 倍，结果，成本较高，信号延时也较大； 如果按相同的选择性和相同的线性要求来说，则 IIR 滤波器就必须加 全通网络进行相位较正，同样要大增加滤波器的节数和复杂性。 B.滤波后信号频谱图的分析 由图 6，FIR 滤波器主要采用非递归结构，不论在理论上还是在 实际的有限精度运算中都不存在稳定性问题， 运算误差也较小。 此外， FIR滤波器可以采用快速付里叶变换算法，在相同阶数的条件下，运 算速度可以快得多。可以看出，经两种滤波器滤波后输出的信号频谱， FIR滤波后信号的频谱更接近实际无噪声信号的频谱。 clear; close all ; clc n=1:8:8000; %以 8 为间隔取 1000 个点 s(n)=cos(n/100); %噪声信号 x(n)=cos(n/100)+0.01*sin(n); %原始信号 n1=1:512; %对采样后序列进行 512 点傅里叶变换 fk=1:8000/512:8000; %横坐标转化为频率 %无噪信号和原始信号 figure(1) subplot(1,2,1) plot(n,s(n)) xlabel( n ); ylabel( s(n) ); title( 无噪信号 ); hold on subplot(1,2,2) plot(n,x(n)) xlabel( n ); ylabel( x(n) ); title( 原始信号 ); fw=fft(s(n),512); fw1=fft(x(n),512); mag=abs(fw); mag1=abs(fw1); %原始信号频谱分析 figure(2) subplot(1,2,1) plot(fk/8