【正文】 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。 fSum=0。 for ( i=0。iFIRNUMBER。i++ ) { fSum+=(fXn[i]*fHn[i])。 } return(fSum)。}程序2include 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 *ebsr = 0xa01。 *egcr = 0x220。 *egcr = 0X220。 *ce01 = 0X3000。 *ce11 = 0X1fff。 *ce21 = 0x1fff。 *ce31 = 0x1fff。 *emirst = 0。 *sdc1 = 0X5958。 *sdc2 = 0X38F。 *init = 0。} 程序3includevoid 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。 *TMCR_MM =0x000。}程序4wst=*pi。wp=*pi。deltaw=wstwp。 %過渡帶寬△w的計算 N0=ceil(11*pi/deltaw)。 %按布萊克曼窗計算所需的濾波器長度 N=N0+mod(N0+1,2)。 %為了實現(xiàn)第一類偶對稱濾波器，應(yī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è)計的脈沖響應(yīng)應(yīng)為理想脈沖響應(yīng)與窗函數(shù)乘積 [db,mag,pha,grd,w]=myfreqz(h,[1])。 %對設(shè)計結(jié)果進行檢驗 dw=2*pi/2000。 %頻率分辨率 Rp=(min(db(1:wp/dw+1))) %檢驗通帶波動 As=round(max(db(wst/dw+1:501))) %檢驗最小阻帶衰減 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。實際脈沖響應(yīng)39。) %繪制實際脈沖響應(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。)。ylabel(39。相位(Φ)39。)程序5function hd=ideal(wc,M)。al=(M1)/2。n=[0:(M1)]。m=nal+eps。hd=sin(wc*m)./(pi*m)。程序5function [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)。25