基于dsp芯片的自適應濾波器實現(xiàn)課程設計(文件)
【正文】 y(n),這樣就可以從z(n)中減去y(n),得到所要提取的信號s(n)的估計值e(n),表示為：e(n)=z(n)y(n)=s(n)+d(n)y(n)。e(n)有兩種作用：一是得到信號s(n)的最佳估計；二是用于調整濾波器系數(shù)的誤差信號。通常有兩種判據(jù)條件：最小均方誤差判據(jù)和最小二乘法判據(jù)。理想信號 與濾波器輸出之差的期望值最小，并且根據(jù)這個判據(jù)來修改權系數(shù)。由式(1)可見，F(xiàn)IR濾波算法實際上是一種乘法累加運算。由表達式可以看出 輸出是兩矢量的內積。其完整代碼是int FIRLMS(int *nx,float *nh,int nError,int nCoeffNumber){ int i,r。inCoeffNumber。 } r/=128。這里實現(xiàn)自適應算法中的公式。189。這里用到了輸入信號數(shù)組xx[],輸出信號數(shù)組rr[],誤差數(shù)組wc[]，以及濾波器權系數(shù)h[]。iINPUTNUMBER。 }用FIRLMS函數(shù)得到了一個濾波后的樣點信號值nLastOutput 存放在輸出信號數(shù)組rr[]中。但是同時還有很多的問題還有待于進一步深入研究，我將在今后的學習工作中，要加強學習，不斷進取。float h[COEFFNUMBER],fU。 fU=。 for ( i=0。 } for ( i=COEFFNUMBER+1。 wc[i]=rr[i]xx[i]。 r=0。 nh[i]+=fWork。}/****************************************************************************** func .h 頭文件******************************************************************************/void interrupt tint( void )。void TIME_init(void)。 // DSP operating clock int clkin。 // DSP clocks per usec int nullloopclk。 ~(mask))define SetMask(addr,mask) addr = (addr | (mask)) // Read and write data at addr, only bits in mask are affecteddefine ReadMask(addr,mask) (addr amp。 ~(mask)) | (data firstbit(mask))extern int firstbit(unsigned short mask)。 ~(mask)) | (data)// Read and write data at addr, shift data so that bit 0 of data is aligned// with lowest set bit in mask. Only bits in mask are modified. Example// with address 0x1000 initial contents = 0x4007:// WriteField(0x1000, 0x56, 0x0ff0) will set address 0x1000 to 0x4567.define ReadField(addr,mask) ((addr amp。extern DSPCLK dspclk。 // PLL multiplier int plldiv。void SetDSPPLL(unsigned int uPLL)。void SDRAM_init( void )。 } r/=128。inCoeffNumber。}int FIRLMS(int *nx,float *nh,int nError,int nCoeffNumber){ int i,r。i++ ) { nLastOutput=FIRLMS(xx+i,h,nLastOutputxx[i1],
外語相關推薦
鄂ICP備17016276號-1