用凯泽窗设计线性相位带通FIR滤波器由刀豆文库小编整理,希望给你工作、学习、生活带来方便,猜你可能喜欢“fir线性相位滤波器”。
用凯泽窗设计线性相位带通FIR滤波器
一 源程序:
1.计算理想带通滤波器单位抽样响应ha(n)的M函数文件dbpfr.m function hd=dbpfr(omega,N)dw=omega(2)-omega(1);alfas=40;M=ceil((alfas-7.95)/(2.286*dw));if mod(M,2)==0;N=M-1;else N=M;end
n=[0:N-1];wc1=1/2*(omega(1)+omega(2));wc2=1/2*(omega(3)+omega(4));alpha=(N-1)/2;hd=sin(wc2*(n-alpha))./(pi.*(n-alpha))-sin(wc1*(n-alpha))./(pi.*(n-alpha));jj=find(n==alpha);hd(jj)=(wc2-wc1)/pi;2.计算凯泽窗函数wk(n)的M函数文件kwnd.m function[wk,N]=kwnd(omega,alfas)dw=omega(2)-omega(1);alfas=40;M=ceil((alfas-7.95)./(2.286*dw));if mod(M,2)==0 N=M-1;else N=M;end if alfas>=50 beta=0.1102*(alfas-8.7)else if alfas(eps*s)T=x/(2*n);n=n+1;D=D.*T.^2;s=s+D;end 4.TTLBQ.m程序文件.通带允许起伏:Rp=-1dB wp1=0.3*pi wp2=0.5*pi 阻带衰减: As
三 思考题:
1. 窗函数的指标包括:旁瓣峰值幅度、主瓣宽度、过渡带宽度、阻带最小衰减 2. 根据阻带衰减可以选择窗函数的类型,根据过渡带宽计算窗函数的长度。3. N取奇数和偶数时,有区别,实验中主要使用N为奇数。
4. 一般选取4~9范围内的值,越大,凯泽窗窗宽越窄,频谱的旁瓣就越小,但是主瓣宽度也相应增加。