【導(dǎo)讀】通過本次課程設(shè)計(jì)，要求學(xué)生對(duì)DSP這門課程所學(xué)的基礎(chǔ)知識(shí)、理論加以理解，初步培養(yǎng)學(xué)生對(duì)工程設(shè)計(jì)的獨(dú)立工作能力，掌握電子系統(tǒng)設(shè)計(jì)的一般方法。培養(yǎng)學(xué)生的動(dòng)手能力，以及團(tuán)隊(duì)合作能力。通過DSP處理器通過A/D采集音頻信號(hào)，在CCS軟件中分析音頻信號(hào)的頻譜圖，DSP芯片主要是將轉(zhuǎn)化成數(shù)字信號(hào)的語音信號(hào)用DSP算法。對(duì)其進(jìn)行處理，并將處理后的信號(hào)傳送到輸出端。沖激響應(yīng)濾波器,在結(jié)構(gòu)上是非遞歸型的。它可以在幅度特性隨意設(shè)計(jì)的同時(shí),保證精確。嚴(yán)格的線性相位。結(jié)構(gòu)中極性震蕩等不穩(wěn)定現(xiàn)象，誤差較小。波器，它在通帶上具有單位增益和線性相位，在阻帶上具有零響應(yīng)。聲音信號(hào)經(jīng)ADC轉(zhuǎn)換后得到的數(shù)字信號(hào)，輸入到數(shù)據(jù)緩沖區(qū)中，

　　

【正文】 Delay(int numbers) { int i,j。 for(i=0。i4000。i++) for(j=0。jnumbers。j++)。 } void kfft(pr,pi,n,k,fr,fi,l,il) Uint16 n,k,l,il。 double pr[],pi[],fr[],fi[]。 { 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]=。 基于 FIR的語音信號(hào)濾波 第 19 頁 共 22 頁 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])*。 } return。 基于 FIR的語音信號(hào)濾波 第 20 頁 共 22 頁 } void main(void) { int i,j,n,k,p,s=0。 double x[128],pr[128],pi[128],fr[128],fi[128],mo[128]。 double fs,fstop,r,rm。 int xm,zm,ym。 n=128。 // PLL, WatchDog, enable Peripheral Clocks // This example function is found in the file. InitSysCtrl()。 // Setup only the GP I/O only for McBSP functionality EALLOW。 // Allow access to EALLOW protected registers =0x3F0F。 //as spi and mcbsp pin EDIS。 // Disable access to EALLOW protected registers // Disable CPU interrupts DINT。 // Initialize PIE control registers to their default state. // The default state is all PIE interrupts disabled and flags // are cleared. InitPieCtrl()。 // Disable CPU interrupts and clear all CPU interrupt flags: IER = 0x0000。 IFR = 0x0000。 // Initialize the PIE vector table with pointers to the shell Interrupt // Service Routines (ISR). InitPieVectTable()。 spi_init()。 aic23_init()。 //初始化 TLV320AIC23，設(shè)置內(nèi)部寄存器 mcbsp_init()。 // For this example, only initialize the Mcbsp fs = 8000。 fstop = 3000。 npass = fstop/fs。 for (i=0。 i=m。 i++) { xmid[i]=0。 } 基于 FIR的語音信號(hào)濾波 第 21 頁 共 22 頁 for(。) { //模擬音頻輸入 for(i=0。i=n1。i++) { mcbsp_read_rdy()。 //MCBSP 接收一個(gè)數(shù)據(jù) 32 位 read_data1 = 。 //右聲道 PA[i]=read_data1。 } //原信號(hào) FFT 轉(zhuǎn)換 for (j=0。 j=4。 j++) { for (i=0。 i=n1。 i++) { xm=PA[i]。 x[i]=xm/。 pr[i]=x[i]。 pi[i]=0。 } kfft(pr,pi,128,7,fr,fi,0,1)。 for (i=0。i=n1。i++) { mo[i] = sqrt(fr[i]*fr[i]+fi[i]*fi[i])。 zm = (int)(mo[i]*)。 PC[i] = zm。 } PC[0]=0。 k++。 } //原信號(hào)濾波 firdes(m, npass)。 for (s=0。 s=4。 s++) { for (i=0。 i=n1。 i++) { xm = PA[i]。 u = xm/。 for (p=0。 p=m。 p++) { xmid[mp] = xmid[mp1]。 } xmid[0] = u。 r = 0。 基于 FIR的語音信號(hào)濾波 第 22 頁 共 22 頁 rm= 0。 for (j=0。 j=m。 j++) { r = xmid[j] * h[j]。 rm = rm + r。 } y = rm。 ym = (int)( * y)。 PB[i] = ym。 } } //濾波后的 FFT 轉(zhuǎn)換 for (j=0。 j=4。 j++) { for (i=0。 i=n1。 i++) { xm=PB[i]。 x[i]=xm/。 pr[i]=x[i]。 pi[i]=0。 } kfft(pr,pi,128,7,fr,fi,0,1)。 for (i=0。i=n1。i++) { mo[i] = sqrt(fr[i]*fr[i]+fi[i]*fi[i])。 zm = (int)(mo[i]*)。 PD[i] = zm。 } PD[0]=0。 k++。 } //原信號(hào)右聲道輸出，濾波后的左聲道輸出 For(i=0。 i=n1。 i++) { mcbsp_write_rdy(PA[i], PB[i])。 } } }