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