用凯泽窗设计线性相位带通FIR滤波器.docVIP

用凯泽窗设计线性相位带通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=21 beta=0 else beta=0.5842*(alfas-21).^0.4+0.07886*(alfas-21) end end n=[1:N-2]; x1=beta*sqrt(1-(1-2.*n./(N-1)).^2) x2=beta; s1=besl(x1); s2=besl(x2); wk1=1./s2; wk2=s1./s2; wk3=1./s2; wk=[wk1 wk2 wk3]; 3.计算贝塞尔函数的参考程序besl.m function s=besl(x) eps=10^(-9); n=1;s=1; D=1; while D(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=40dB ws1=0.15*pi ws2=0.65*pi clc; clear; Rp =-1; alfas=40; omega=[0.15*pi 0.3*pi 0.5*pi 0.65*pi]; [wk,N]=kwnd(omega,alfas) hd=dbpfr(omega,N) figure(1),clf; stem(hd);grid; title( hd(n)); figure(2),clf; stem(wk);grid; title(w(n)); hn=hd.*wk figure(3),clf stem(hn);grid; title(hd(n)); Hn=fft(hn,128); figure(4),clf subplot(1,2,1); plot(abs(Hn(1:64)));grid; title(幅频响应曲线); subplot(1,2,2); plot(20.*log(abs(Hn(1:32))));grid; title(对数幅频响应曲线); Hd=fft(hd,128); figure(5),clf plot(abs(Hd(1:65)));grid; title(理想幅频响应曲线); 二 截图 三 思考题： 窗函数的指标包括：旁瓣峰值幅度、主瓣宽度、过渡带宽度、阻带最小衰减 根据阻带衰减可以选择窗函数的类型，根据过渡带宽计算窗函数的长度。 N取奇数和偶数时，有区别，实验中主要使用N为奇数。 ?一般选取4~9范围内的值，?越大，凯泽窗窗宽越窄，频谱的旁瓣就越小，但是主瓣宽度也相应增加。