fir滤波器线性相位 用凯泽窗设计线性相位带通FIR滤波器.doc
文本预览下载声明
fir滤波器线性相位 用凯泽窗设计线性相位带通FIR滤波器
用凯泽窗设计线性相位带通FIR滤波器 一 源程序: 1.计算理想带通滤波器单位抽样响应ha(n)的M函数文件dbpfr.m function hd=dbpfr(omega,N) dw=omega(2)-omega(1); alfas=40; M=ce…
毛泽东思想邓小平理论和三个代表重要思想概论考试习题 一、单项选择题: 1、在中国共产党的历史上,第一次使用“思想路线”这一概念的是( ) A陈独秀 B李大钊 C毛泽东 D刘少奇 2、明确提出“关于主观主义”的问题的文献是( ) A《古田会议决议》 B…
人们常说,教师是太阳底下最光辉的职业。从事这个行业已经15年了,担任班主任工作也已经11年了。在这十多年中有苦也有乐,尤其是从事班主任工作,看着一批一批的学生经过自己的精心培育,升入高一年级,心里比蜜还甜。在自己的教育生涯中,我用自己的爱心去精心浇灌…
用凯泽窗设计线性相位带通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 alfasbeta=0
else beta=0.5842*(alfas-21). .4+0.07886*(alfas-21)
end
end
n=[1:N-2];
x1=beta*sqrt(1-(1-2.*n./(N-1)).)
x2=beta;
s1=besl(x1);
s2=besl(x2);
wk1=1./s2;
wsclc;
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(‘理想幅频响应曲线’);
二 截图
三 思考题:
1. 窗函数的指标包括:旁瓣峰值幅度、主瓣宽度、过渡带宽度、阻带最小衰减
2. 根据阻带衰减可以选择窗函数的类型,根据过渡带宽计算窗函数的长度。
3. N取奇数和偶数时,有区别,实验中主要使用N为奇数。
4. ?一般选取4~9范围内的值,?越大,凯泽窗窗宽越窄,频谱的旁瓣就越小,但是主瓣宽度也相应增加。
用凯泽窗设计线性相位带通FIR滤波器 一 源程序: 1.计算理想带通滤波器单位抽样响应ha(n)的M函数文件dbpfr.m function hd=dbpfr(omega,N) dw=omega(2)-omega(1); alfas=40
显示全部