基于DSP的带阻滤波器设计毕业设计论文

基于DSP的带阻滤波器设计

参考文献

[1]潘松，黄继业，王国栋编著. 现代DSP技术. 西安：西安电子科技大学出版社, 2003. [2]赵刚，黄建明，隋燕等编著. 基于数字滤波器设计的讨论[J]. 天津：南开大学学报(自然科学版)，2003.26(3):15218.

[3]张雄伟，徐亮编著. DSP集成开发与应用实例. 北京：电子工业出版社，2002，28-33. [4]程佩青编著. 数字信号处理教程[M]. 北京：清华大学出版社，2000，3~15.

[5]余成波，杨菁，杨如民编著. 数字信号处理及MATLAB实现. 北京：清华大学出版社， 2005，90~127.

[6]李正周编著. MATLAB数字信号处理与应用. 北京：清华大学出版社，2008，93~119. [7]孙宗瀛编著. DSP原理设计与应用[M]. 北京：清华大学出版社，2002. [8]邹彦编著. DSP原理及应用. 北京：电子工业出版社，2004，40~43，279~330.

[9]张雄伟，陈亮，徐光辉编著. DSP集成开发与应用实例[M]. 北京：电子工业出版社，2002. [10]清源科技编著.DSP应用开发教程系列TMS32054xDSP硬件开发教程.北京：机械工业 出版社，2003，132-155.

[11]TMS320C54X DSP REFERRENCE SET VOL.3:Algebraic Instruction Set. Texas

Instruments Incorporated. 2001.

[12]程佩青编著. 数字滤波与快速傅里叶变换. 北京：清华大学出版社，1988. [13]胡广书编著. 数字信号处理理论、算法与实现. 北京：清华大学出版社，1997. [14]TMS320C54x Code Composer Studio Tutorial. Texas Instruments Incorporated. 2000. [15]张雄伟，曹铁勇编著. DSP芯片的原理与开发应用. 北京：电子工业出版社，2000. [16]周霖编著. DSP算法设计与系统方案. 北京：国防工业出版社，2004，7.

32

附录1：原理图

附 录

基于DSP的带阻滤波器设计

33

基于DSP的带阻滤波器设计

附录2：MATLAB程序 clear all; fpl=20; fph=80; fsl=40; fsh=60; Fs=10000; Wpl=2*pi*fpl/Fs Wph=2*pi*fph/Fs Wsl=2*pi*fsl/Fs Wsh=2*pi*fsh/Fs

tr_width=min((Wsl-Wpl),(Wph-Wsh)); %过渡带宽度 N=ceil(6.2*pi/tr_width) %滤波器长度 n=0:1:N-1;

Wcl=(Wsl+Wpl)/2; %理想低通滤波器的截止长度 Wch=(Wsh+Wph)/2; alpha=(N-1)/2; m=n-alpha+eps;

hd=[sin(Wcl*m)+sin(pi*m)-sin(Wch*m)]./(pi*m);

w_han=(hanning(N))'; %汉宁窗

h=hd.*w_han; %截取得到实际的单位脉冲响应 [H,w]=freqz(h,[1],1000,'whole'); H=(H(1:1:501))'; w=(w(1:1:501))';

mag=abs(H); %绝对幅值响应 db=20*log10((mag+eps)/max(mag)); %相对幅值响应 pha=angle(H); %相位响应 delta_w=2*pi/1000;

Ap=-(min(db(1:1:Wpl/delta_w+1))) %实际通带波纹 As=-round(max(db(Wsl/delta_w+1:1:Wsh/delta_w+1))) %实际阻带波纹

34

基于DSP的带阻滤波器设计

subplot(211) plot(w/pi,db)

title('汉宁窗幅度响应(dB)') axis([0,0.5,-50,10]) subplot(212) plot(w/pi,pha/pi) title('汉宁窗相位响应') axis([0,1,-5,5])

35