FIR滤波器解析.pptVIP

滤波器部分 电子信息工程一班 数字信号讨论课 二组 目录 一：线性相位滤波器特点 二：利用窗函数设计FIR滤波 器 三：吉布斯效应及改进措施 四：理想低通滤波器和实际低通滤波器特点 张三 五：模拟信号采集 六：信号频域分析 七：实现信号滤波 张三 线性相位滤 波器特点 第一类线性相位条件： h(n)=h(N-1-n) 0≤n≤N-1 即h(n)关于(N-1)/2成偶对称 第二类线性相位条件： h(n)=-h(N-1-N) 0≤n≤N-1 即h(n)关于(N-1)/2成奇对称 窗函数设计滤波器 IDFT DFT 图形解释 张三 时 域 相 乘 频 域 卷 积 C语言实现窗函数设 计FIR滤波器 #includegraphics.h #includestdio.h #includeconio.h #includemath.h #define PI 3.1415926 #define MAX 2048 #define N 51 float Hr[MAX],Hi[MAX],Hk[MAX], Hkxw[MAX],Hd[MAX],wn[MAX]; main() { int i,m,k,n; void GRAPH(); void DFT_FFT(); void HDN(),JXC(),HNC(),HMC(),BLKMC(); m=8; n=(int)pow(2,m);//FFT的长度 HDN(n); printf(Please input xh:(1-Rectangle; 2-Hanning;3-Haming;4-Blackman):); scanf(%d,k); switch(k) { case 1:JXC(n);break; case 2:HNC(n);break; case 3:HMC(n);break; case 4:BLKMC(n);break; } for(i=0;in;i++) { Hr[i]=Hd[i]*wn[i]; Hi[i]=0.0; } 滤波器类型 张三 void HDN(int n) { int i; float a; a=(N-1)/2； for(i=0;iN;i++) { if(i==a)Hd[i]=0.25; else Hd[i]=sin((PI/4.0)*(i-a))/(float)(PI*(i-a)); } for(i=N;in;i++) {Hd[i]=0.0;} } {Hd[i]=0.0;} } 偶函数---第一类线性相位 N=奇数， 可实现 低通 带通 高通 带阻 矩形窗对不同滤波器的实现 低通 第一类第一种情况 void HDN(int n) { int i; float a; a=(N-1)/2； for(i=0;iN;i++) { if(i==a)Hd[i]=0.25; else Hd[i]=sin((PI/4.0)*(i-a)) /(float)(PI*(i-a))； } for(i=N;in;i++) {Hd[i]=0.0;} } 矩形窗对不同滤波器的实现 第一类第一种情况 高通 void HDN(int n) { int i; float a; a=(N-1)/2； for(i=0;iN;i++) { if(i==a)Hd[i]=0.75; else Hd[i]=-1.21*sin((3*PI/7.5)* (i-a))/(float)(PI*(i-a)); } for(i=N;in;i++) {Hd[i]=0.0;} } 矩形窗对不同滤波器的实现 第一类第一种情况 带通 void HDN(int n) { int i; float a; a=(N-1)/2； for(i=0;iN;i++) { if(i==a)Hd[i]=0.25; else Hd[i]=(sin((PI/2.0)*(i-a)) -sin((PI/4.0)*(i-a)))/(float)(PI*(i-a)); } for(i=N;in;i++) {Hd[i]=0.0;} } w(n)=1/2[1-cos(2πn/(N-1))]RN(n) 低通， 阻带衰减约 为50dB Hanning窗 矩 形 窗 低通，截止 频率为π/4， 阻带衰减约为20dB w(n)=RN(n) w(n)=0.54-0.46cos[2πn/(N-1)]RN(n) 低通，