【導(dǎo)讀】信息、無(wú)線電、自動(dòng)控制等得到廣泛應(yīng)用。本文基于DSP完成了FIR低通濾波器。介紹了FIR數(shù)字濾波的特點(diǎn)和設(shè)計(jì)方法以及FIR低通濾波器的設(shè)計(jì)過程。最終基于DSP在TMS320VC5009平臺(tái)上用C語(yǔ)言編程在CCS軟件上。50KHz的FIR低通濾波器。

　　

【正文】 3] 王獻(xiàn)峰，石東 .基于 DSP 的 FIR 數(shù)字濾波器的設(shè)計(jì)與實(shí)現(xiàn) . 天津科技大學(xué)電子信息與自動(dòng)化學(xué)院，天津 300222. [4] 冉茂華 .基于 DSP的 FIR數(shù)字濾波器設(shè)計(jì) ,武漢理工大學(xué) .20xx. [5] 劉萬(wàn)松 .基于 DSP的 FIR數(shù)字濾波器設(shè)計(jì)，貴州大學(xué) .20xx. [6] 孫宗瀛 . DSP 原理設(shè)計(jì)與應(yīng)用 [M]. 北京：清華大學(xué)出版社， 20xx. [7] 王獻(xiàn)峰，石東 . 基于 DSP 的 FIR 數(shù)字濾波器的設(shè)計(jì)與實(shí)現(xiàn) . 天津科技大學(xué)電子信息與自動(dòng)化學(xué)院，天津 300222. [8]張雄偉 ,曹鐵勇 . DSP 芯片的原理與開發(fā)應(yīng)用 . 電子工業(yè)出版社 . 20xx. 21 附錄 程序 1 include include include include define FIRNUMBER 25 define SIGNAL1F 1000 define SIGNAL2F 4500 define SAMPLEF 10000 define PI float InputWave()。 float FIR()。 float fHn[FIRNUMBER]={ ,, ,, ,, , }。 float fXn[FIRNUMBER]={ }。 float fInput,fOutput。 float fSignal1,fSignal2。 float fStepSignal1,fStepSignal2。 float f2PI。 int i。 float fIn[256],fOut[256]。 int nIn,nOut。 main() { nIn=0。 nOut=0。 f2PI=2*PI。 22 fSignal1=。 fSignal2=PI*。 fStepSignal1=2*PI/30。 fStepSignal2=2*PI*。 while ( 1 ) { fInput=InputWave()。 fIn[nIn]=fInput。 nIn++。 nIn%=256。 fOutput=FIR()。 fOut[nOut]=fOutput。 nOut++。 /* break point */ if ( nOut=256 ) { nOut=0。 } } } float InputWave() { for ( i=FIRNUMBER1。i0。i ) fXn[i]=fXn[i1]。 fXn[0]=sin((double)fSignal1)+cos((double)fSignal2)/。 fSignal1+=fStepSignal1。 if ( fSignal1=f2PI ) fSignal1=f2PI。 fSignal2+=fStepSignal2。 if ( fSignal2=f2PI ) fSignal2=f2PI。 return(fXn[0])。 } float FIR() { float fSum。 23 fSum=0。 for ( i=0。iFIRNUMBER。i++ ) { fSum+=(fXn[i]*fHn[i])。 } return(fSum)。 } 程序 2 include void SDRAM_init( void ) { ioport unsigned int *ebsr =(unsigned int *)0x6c00。 ioport unsigned int *egcr =(unsigned int *)0x800。 ioport unsigned int *emirst=(unsigned int *)0x801。 //ioport unsigned int *emibe =(unsigned int *)0x802。 ioport unsigned int *ce01 =(unsigned int *)0x803。 //ioport unsigned int *ce02 =(unsigned int *)0x804。 //ioport unsigned int *ce03 =(unsigned int *)0x805。 ioport unsigned int *ce11 =(unsigned int *)0x806。 //ioport unsigned int *ce12 =(unsigned int *)0x807。 //ioport unsigned int *ce13 =(unsigned int *)0x808。 ioport unsigned int *ce21 =(unsigned int *)0x809。 //ioport unsigned int *ce22 =(unsigned int *)0x80A。 //ioport unsigned int *ce23 =(unsigned int *)0x80B。 ioport unsigned int *ce31 =(unsigned int *)0x80C。 //ioport unsigned int *ce32 =(unsigned int *)0x80D。 //ioport unsigned int *ce33 =(unsigned int *)0x80E。 ioport unsigned int *sdc1 =(unsigned int *)0x80F。 //ioport unsigned int *sdper =(unsigned int *)0x810。 //ioport unsigned int *sdt =(unsigned int *)0x811。 ioport unsigned int *init =(unsigned int *)0x812。 ioport unsigned int *sdc2 =(unsigned int *)0x813。 //*ebsr = 0x221。//0xa01 24 *ebsr = 0xa01。 *egcr = 0x220。 *egcr = 0X220。 *ce01 = 0X3000。 *ce11 = 0X1fff。 *ce21 = 0x1fff。 *ce31 = 0x1fff。 *emirst = 0。 *sdc1 = 0X5958。 *sdc2 = 0X38F。 *init = 0。 } 程序 3 include void CLK_init() { ioport unsigned int *clkmd。 clkmd=(unsigned int *)0x1c00。 *clkmd =0x2033。 // 0x2033。//0x2413。// 200MHz=0x2513 } void SetDSPPLL(unsigned int uPLL) { ioport unsigned int *clkmd。 clkmd=(unsigned int *)0x1c00。 *clkmd =uPLL。 } void TMCR_reset( void ) { ioport unsigned int *TMCR_MGS3=(unsigned int *)0x07FE。 ioport unsigned int *TMCR_MM =(unsigned int *)0x07FF。 *TMCR_MGS3 =0x510。 25 *TMCR_MM =0x000。 } 程序 4 wst=*pi。 wp=*pi。 deltaw=wstwp。 %過渡帶寬△ w的計(jì)算 N0=ceil(11*pi/deltaw)。 %按布萊克曼窗計(jì)算所需的濾波器長(zhǎng)度 N=N0+mod(N0+1,2)。 %為了實(shí)現(xiàn)第一類偶對(duì)稱濾波器，應(yīng)確保其長(zhǎng)度 N為奇數(shù) n=[0:1:N1]。 wc=(wst+wp)/2。 %截止頻率取為兩邊緣頻率的平均值 hd=ideal(wc,N)。 %求理想脈沖響應(yīng) wdbla=(blackman(N))39。 %求窗函數(shù) h=hd.*wdbla。 %設(shè)計(jì)的脈沖響應(yīng)應(yīng)為理想脈沖響 應(yīng)與窗函數(shù)乘積 [db,mag,pha,grd,w]=myfreqz(h,[1])。 %對(duì)設(shè)計(jì)結(jié)果進(jìn)行檢驗(yàn) dw=2*pi/20xx。 %頻率分辨率 Rp=(min(db(1:wp/dw+1))) %檢驗(yàn)通帶波動(dòng) As=round(max(db(wst/dw+1:501))) %檢驗(yàn)最小阻帶衰減 figure(1) subplot(2,2,1)。 stem(n,wdbla)。 title(39。布萊克曼窗 39。) %繪制布萊克曼窗 axis([0 N1 0 ])。ylabel(39。w(n)39。)。 text(N+1,0,39。n39。) subplot(2,2,2)。 stem(n,h)。 title(39。實(shí)際脈沖響應(yīng) 39。) %繪制實(shí)際脈沖響應(yīng) axis([0 N1 ])。 xlabel(39。n39。)。 ylabel(39。h(n)39。) subplot(2,2,3)。 plot(w/pi,db)。 %繪制幅度響應(yīng) 。 title(39。幅度響應(yīng)（單位： dB)39。)。 grid。 axis([0 1 150 10])。 xlabel(39。頻率（單位： pi)39。)。 ylabel(39。分貝數(shù) 39。) set(gca,39。XTickMode39。,39。manual39。,39。XTick39。,[0,,1]) set(gca,39。YTickMode39。,39。manual39。,39。YTick39。,[60,0]) subplot(2,2,4)。 plot(w,pha)。 title(39。相位響應(yīng) 39。)。 %繪制相頻特性 axis([0 2 4 4])。 xlabel(39。頻率（單位 :pi)39。)。 26 ylabel(39。相位 (Φ )39。) 程序 5 function hd=ideal(wc,M)。 al=(M1)/2。 n=[0:(M1)]。 m=nal+eps。 hd=sin(wc*m)./(pi*m)。 程序 5 function [db,mag,pha,grd,w]=myfreqz(b,a)。 N=1000。 [H,w]=freqz(b,a,N,39。whole39。)。 H=(H(1:1:501))39。w=(w(1:1:501))39。 mag=abs(H)。 db=20*log10((mag+eps)/max(mag))。 pha=angle(H)。 grd=grpdelay(b,a,w)。