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

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

一源程序：

1.计算理想带通滤波器单位抽样响应ha(n)的M函数文件dbpfr.m

functionhd=dbpfr(omega,N)

dw=omega(2)-omega(1);

alfas=40;

M=ceil((alfas-7.95)/(2.286*dw));

ifmod(M,2)==0;

N=M-1;

elseN=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));

ifmod(M,2)==0

N=M-1;

elseN=M;

end

ifalfas=50

beta=0.1102*(alfas-8.7)

elseifalfas=21

beta=0

elsebeta=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=[wk1wk2wk3];

3.计算贝塞尔函数的参考程序besl.m

functions=besl(x)

eps=10^(-9);

n=1;s=1;

D=1;

whileD(eps*s)

T=x/(2*n);

n=n+1;

D=D.*T.^2;

s=s+D;

end

4.TTLBQ.m程序文件.

通带允许起伏:Rp=-1dBwp1=0.3*piwp2=0.5*pi

阻带衰减:As=40dBws1=0.15*piws2=0.65*pi

clc;

clear;

Rp=-1;

alfas=40;

omega=[0.15*pi0.3*pi0.5*pi0.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(理想幅频响应曲线);

二截图

三思考题：

1．窗函数的指标包括：旁瓣峰值幅度、主瓣宽度、过渡带宽度、阻带最小衰减

2．根据阻带衰减可以选择窗函数的类型，根据过渡带宽计算窗函数的长度。

3．N取奇数和偶数时，有区别，实验中主要使用N为奇数。

4．一般选取4~9范围内的值，越大，凯泽窗窗宽越窄，频谱的旁瓣就越小，但是主瓣

宽度也相应增加。