工程数字信号处理62.pptVIP

数字滤波技术 FIR滤波器 FIR滤波器又称为有限脉冲响应滤波器，FIR就来源于名词“有限脉冲响应”的英文单词首字符缩写。 FIR传函： 脉冲响应函数： h(k)是一个有限长度的序列，故FIR滤波器称为有限脉冲响应滤波器。 数字滤波技术 FIR滤波器的滤波公式： y(k)=h(0)x(k)+h(1)x(k+1) +...+h(N) x(k+N) y(i)=0.09x(k-3)+0.15x(k-2)+0.20x(k-1)+0.22x(k) + 0.20x(k+1)+0.15x(k+2)+0.09x(k+3) 零相移处理 y(k)=h(0)x(k-N/2)+h(1)x(k-N/2+1) +...+h(N) x(k+N/2) 数字滤波技术 Dim x(2048),y(2048) For K = 0 To 512 x(k)=data(k) y(k)=0.00001 Next For K = 4 To 509 y(k)=0.09*x(k-3)+0.15*x(k-2) +0.2*x(k-1)+0.22*x(k) +0.2*x(k+1)+0.15*x(k+2) +0.09*x(k+3) Next SetInterval 104,dt Document.setarrayline 104,512,y 数字滤波技术 FIR滤波器系数脉冲响应函数设计法 对离散化系统，用归一化频率表示v=F/Fs： 频率特性富里叶级数表达： ？参考波形合成与分解 数字滤波技术 为简化计算，设Hd为偶函数(也可为奇函数)，则： n 取前Q项近似为： 富里叶变换的奇偶特性 令 ，取Q项近似： 数字滤波技术 令： 例如，当Q=5时有： 数字滤波技术 低通滤波器 : 高通滤波器： 数字滤波技术 带通滤波器： 带阻滤波器 ： 数字滤波技术 例如，低通滤波器设FN＝5kHz，f1=1kHz ,有： 取Q=5，可以得到： 数字滤波技术 N=？ fc=? N2=Int(N/2) h(0)=fc For i=1 To N2 Step 1 f=i*pi h(i)=sin(fc*f)/f Next 低通滤波器设计算法： 奇数处理 数字滤波技术 数字滤波技术 理想滤波器Cn是无限长序列，但FIR系数h(n)是有限的，最简单办法是直接截取。这样的截断在时域有跳变，导致滤波器特性变差。可以用一个向两边衰减的窗函数与Cn相乘：h(n)=w(n)C(n)，减小系数截断带来的时域有跳变造成的误差。 FIR滤波器窗函数设计法 数字滤波技术 设计中常用的几种窗函数： 矩形窗： w(n)=RN(n)； Hanning窗： Hamming窗： Blackmen窗： 数字滤波技术 N=？ fc=? N2=Int(N/2) h(0)=fc For i=1 To N2 Step 1 f=i*pi h(i)=sin(fc*f)/f Next f1= 2.0*pi/(N - 1.0) For i=0 To N2 Step 1 h(i)=h(i)*(0.54-0.46*cos(f1*(i+N2))) Next sum=h(0) For i=1 To N2 Step 1 sum=sum+2*h(i) Next For i=0 To N2 Step 1 h(i)=h(i)/sum Next 数字滤波技术 数字滤波技术 FIR滤波器频率取样/优化设计法： Fmin=F-F(C1,C2,C3 ,.. ) 数字滤波技术 MatLab: F = [0 0.2 0.4 0.6 0.8 0.9]; A = [1 1 0.6 0.4 0.2 0.0]; b = firls(24,F,A); [H,f] = freqz(b,1,512,2); plot(f,abs(H)) grid on legend(Ideal,firls Design) Disp(b) 数字滤波技术 实验：用MatLab或DRVI设计一个方波信号发生器，然后用窗函数法设计一个FIR低通/高通/带通滤波器，只保留需要的谐波成分。 数字滤波技术 IIR滤波器又称无限脉冲响应滤波器，IIR来源于名词“无限脉冲响应”的英文单词首字符缩写。 h(k)是一个无限长度的序列，故IIR滤波器称为无限脉冲响应滤波器。 IIR传函： 脉冲响应函数(长除法)： IIR滤波器 数字滤波技术 Butterworth Chebyshev Elliptic kaiser 种类 数字滤波技术 MatLab中的滤波功能 Fs = 44100; dt=1.0/Fs; T =2;