【文章內(nèi)容簡(jiǎn)介】 0xf800。// asm( ssbx intm )。 //Disable all mask interrupts //// *(unsigned int*)IMR=0x0。///**/ *(unsigned int*)IFR=0xffff。// asm( nop )。 asm( nop )。 asm( nop )。 } /************************************************************ 函數(shù)名稱 : void Delay(int numbers) 函數(shù)說明 : 延時(shí) 輸入?yún)?shù) : numbers 輸出參數(shù) : 無************************************************************/void Delay(UINT16 numbers){ UINT16 i,j。 for(i=0。i4000。i++) for(j=0。jnumbers。j++)。 }//// 函數(shù)名稱 : void xint2_init(void)// 函數(shù)說明 : 初始化XINT2// 輸入?yún)?shù) : 無// 輸出參數(shù) : 無//void xint2_init() //外部中斷2初始化子程序{ *(unsigned int*)IMR=0x0004。//使能int2中斷 /* bit 15 1: XINT2 flag write 1 to clear bit 143 0: reserved bit 2 1: XINT2 Polarity 1 rising eage bit 1 0: XINT2 Priority 0 High priority bit 0 1: XINT2 Enable 1 Enable interrupt*/ asm( rsbx INTM)。//開總中斷}//// 函數(shù)名稱 : void kfft(double pr[Length],double pi[Length],int n,int k,double fr[Length],double fi[Length],int l,int il)// 函數(shù)說明 : 基2快速傅立葉變換子程序// 輸入?yún)?shù) : // 輸出參數(shù) : 無//void kfft(double pr[Length],double pi[Length],int n,int k,double fr[Length],double fi[Length],int l,int il){ int it,m,is,i,j,nv,l0。 double p,q,s,vr,vi,poddr,poddi。 for (it=0。 it=n1。 it++) { m=it。 is=0。 for (i=0。 i=k1。 i++) { j=m/2。 is=2*is+(m2*j)。 m=j。} fr[it]=pr[is]。 fi[it]=pi[is]。 } pr[0]=。 pi[0]=。 p=(*n)。 pr[1]=cos(p)。 pi[1]=sin(p)。 if (l!=0) pi[1]=pi[1]。 for (i=2。 i=n1。 i++) { p=pr[i1]*pr[1]。 q=pi[i1]*pi[1]。 s=(pr[i1]+pi[i1])*(pr[1]+pi[1])。 pr[i]=pq。 pi[i]=spq。 } for (it=0。 it=n2。 it=it+2) { vr=fr[it]。 vi=fi[it]。 fr[it]=vr+fr[it+1]。 fi[it]=vi+fi[it+1]。 fr[it+1]=vrfr[it+1]。 fi[it+1]=vifi[it+1]。 } m=n/2。 nv=2。 for (l0=k2。 l0=0。 l0) { m=m/2。 nv=2*nv。 for (it=0。 it=(m1)*nv。 it=it+nv) for (j=0。 j=(nv/2)1。 j++) { p=pr[m*j]*fr[it+j+nv/2]。 q=pi[m*j]*fi[it+j+nv/2]。 s=pr[m*j]+pi[m*j]。 s=s*(fr[it+j+nv/2]+fi[it+j+nv/2])。 poddr=pq。 poddi=spq。 fr[it+j+nv/2]=fr[it+j]poddr。 fi[it+j+nv/2]=fi[it+j]poddi。 fr[it+j]=fr[it+j]+poddr。 fi[it+j]=fi[it+j]+poddi。 } } if (l!=0) for (i=0。 i=n1。 i++) { fr[i]=fr[i]/(*n)。 fi[i]=fi[i]/(*n)。 } if (il!=0) for (i=0。 i=n1。 i++) { pr[i]=sqrt(fr[i]*fr[i]+fi[i]*fi[i])。 if (fabs(fr[i])*fabs(fi[i])) { if ((fi[i]*fr[i])0) pi[i]=。 else pi[i]=。 } else pi[i]=atan(fi[i]/fr[i])*。 } }void firdes(double npass)//fir濾波器，求出單位脈沖響應(yīng){ int pi=。 int t。 for (t=0。 t51。 t++) { if (t == ((511)/2)) { h[t]=2*npass。 } else { h[t] = sin((t(511)/)*npass*pi*2)/(pi*(t(511)/))。//fir濾波器的函數(shù) } } } void Convolveok( //卷積函數(shù) double *x, //原始輸入數(shù)據(jù) double *h, //沖擊響應(yīng) double *y, //卷積輸出結(jié)果 UINT16 length, // 序列長(zhǎng)度 UINT16 fLen) { UINT16 m,p,j。 double r,rm。 double xmean = 。 double xmid[100]。 h[0] = 。 h[50] = 。 for(m=0。mfLen。m++) { xmid[m] = 。 } for (m=0。 mlength。 m++) { xmean = x[m] + xmean。 } xmean = *xmean/length。 for (m=0。 mlength。 m++) { x[m] = x[m] xmean。 } for (m=0。 mlength。 m++) { for (p=0。 pfLen。 p++) { xmid[fLenp1] = xmid[fLenp2]。 } xmid[0] = x[m]。 r = 。 rm= 。 for (j=0。 jfLen。 j++) { r = xmid[j] * h[j]。 rm = rm + r。 } y[m] = rm。 }}//// 函數(shù)名稱 : void delay_50ms(void) // 函數(shù)說明 : 50ms延遲// 輸入?yún)?shù) : 無// 輸出參數(shù) : 無//void delay_50ms(){ UINT16 i,j。 for(i=0。i=1000。i++)。 for(j=0。j=2000。j++)。 //延遲250*1000*CLKOUT=500000*CLKOUT //1/CLKOUT= } //// 函數(shù)名稱 : void delay_20ms(void) // 函數(shù)說明