凯泽窗带通滤波器MATLAB,用凯泽窗设计线性相位带通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<=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('理想幅频响应曲线'); 二 截图

三 思考题：

1． 窗函数的指标包括：旁瓣峰值幅度、主瓣宽度、过渡带宽度、阻带最小衰减 2． 根据阻带衰减可以选择窗函数的类型，根据过渡带宽计算窗函数的长度。 3． N取奇数和偶数时，有区别，实验中主要使用N为奇数。

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

电子信息工程

专业课程设计任务书

学生姓名

xx

专业班级

电信

学号

xxxx

题

目

用双线性变换法设计原型低通为巴特沃兹型的数字IIR高通滤波器

课题性质

其它

课题来源

自拟课题

指导教师

xxx

同组姓名

xxxx

主要内容

用双线性变换法设计原型低通为巴特沃兹型的数字IIR高通滤波器，要求通带边界频率为500Hz，阻带边界频率分别为400Hz，通带最大衰减1dB，阻带最小衰减40dB，抽样频率为2000Hz，用MATLAB画出幅频特性，画出并分析滤波器传输函数的零极点；

信号经过该滤波器，其中300Hz，600Hz，滤波器的输出是什么？用Matlab验证你的结论并给出的图形。

任务要求

1、掌握用双线性变换法设计原型低通为巴特沃兹型的数字IIR高通滤波器的原理和设计方法。

2、求出所设计滤波器的Z变换。

3、用MATLAB画出幅频特性图。

4、验证所设计的滤波器。

参考文献

1、程佩青著，《数字信号处理教程》，清华大学出版社，2001

2、Sanjit

K.

Mitra著，孙洪，余翔宇译，《数字信号处理实验指导书(MATLAB版)》，电子工业出版社，2005年1月

3、郭仕剑等，《MATLAB

7.x数字信号处理》，人民邮电出版社，2006年

4、胡广书，《数字信号处理

理论算法与实现》，清华大学出版社，2003年

审查意见

指导教师签字：

教研室主任签字：

年

月

日

说明：本表由指导教师填写，由教研室主任审核后下达给选题学生，装订在设计(论文)首页

1

需求分析

当今，数字信号处理(DSP：Digtal

Signal

Processing)技术正飞速发展，它不但自成一门学科，更是以不同形式影响和渗透到其他学科：它与国民经济息息相关，与国防建设紧密相连；它影响或改变着我们的生产、生活方式，因此受到人们普遍的关注。

数字化、智能化和网络化是当代信息技术发展的大