dsp應(yīng)用課程設(shè)計(jì)--自適應(yīng)濾波器的設(shè)計(jì)與實(shí)現(xiàn)

1、<p><b>　　DSP應(yīng)用技術(shù)課程</b></p><p><b>　　設(shè)計(jì)說(shuō)明書</b></p><p>　　起止日期： 2013年12月2日 至 2013年12月13日</p><p>　自適應(yīng)濾波器的設(shè)計(jì)與實(shí)現(xiàn)</p><p>　課題名稱自適應(yīng)濾波器的設(shè)計(jì)與實(shí)現(xiàn)人 數(shù)9<

2、;/p><p>　組 長(zhǎng)同組人員</p><p>　課題的主要內(nèi)容和要求1、主要內(nèi)容：1)、根據(jù)實(shí)驗(yàn)的經(jīng)驗(yàn)編寫程序，不懂的可查閱資料。 2)、程序完善、優(yōu)化，面向C6000開發(fā)調(diào)試。3)、能演示輸入輸出數(shù)據(jù)的完整圖形，顯示結(jié)果，并對(duì)程序或相關(guān)設(shè)計(jì)予以解釋。2、要求：1)、能夠熟練運(yùn)用CCS2.2仿真編譯環(huán)境進(jìn)行軟件仿真并檢查和改正程序錯(cuò)誤。2)、理解濾波器的原理及優(yōu)缺點(diǎn)，看有何可以改進(jìn)的

3、地方。3)、對(duì)程序進(jìn)行再次優(yōu)化，能夠避免一些錯(cuò)誤，得到更好的仿真結(jié)果。</p><p>　具體任務(wù)1、 翻閱書籍查找資料，熟悉開發(fā)軟件CCS2.2和濾波器的結(jié)構(gòu)和組成。2、 到網(wǎng)上收集資料，了解LSM算法及輸入和輸出信號(hào)的關(guān)系及濾波器的特性。3、 理解編寫的程序，并檢查錯(cuò)誤，操作程序，并且通過(guò)仿真軟件查看最后的輸出結(jié)果是否正確。</p><p>　時(shí)間安排與完成情況1、時(shí)間安排：1、組織動(dòng)

4、員與分組，分組討論，1天。2、資料查閱，功能分析，2天。3、程序設(shè)計(jì)，5天。4、驗(yàn)證與調(diào)試，2天。5、書寫課程設(shè)計(jì)說(shuō)明書，1天。2、完成情況;（1）程序已經(jīng)編好，并已經(jīng)排查掉錯(cuò)誤，能正常實(shí)現(xiàn)實(shí)驗(yàn)效果。</p><p><b>　　注：本頁(yè)可以多頁(yè)。</b></p><p><b>　　附錄：</b></p><p><

5、;b>　　程序清單：</b></p><p><b>　　//exp8a.c</b></p><p>　　#include "LP_coef.dat"</p><p>　　#define N0 48 /*自適應(yīng)濾波器的階數(shù)*/</p><p>　　#define N1 48 /*

6、未知濾波器的階數(shù)*/</p><p>　　#define Ns 128 /*輸入信號(hào)的數(shù)量*/</p><p>　　#define TENMU 0x1000 /*步長(zhǎng)*/</p><p>　　#pragma CODE_SECTION(main, "lms_code");</p><p>　　#pragma DATA_SEC

7、TION(fir_index, "lms_data");</p><p>　　#pragma DATA_SECTION(sys_index, "lms_data");</p><p>　　#pragma DATA_SECTION(w, "lms_coef");</p><p>　　#pragma DATA_

8、SECTION(s, "lms_data");</p><p>　　#pragma DATA_SECTION(f, "lms_data");</p><p>　　#pragma DATA_SECTION(in, "lms_in");</p><p>　　#pragma DATA_SECTION(d, &qu

9、ot;lms_out");</p><p>　　extern void init(int *, unsigned int);</p><p>　　extern void random(int *, unsigned int);</p><p>　　int w[N0], /*自適應(yīng)濾波器的系數(shù)*/</p><p>　　s[N0]

10、, /*自適應(yīng)濾波延遲樣品緩沖*/</p><p>　　f[N1], /*未知系統(tǒng)延遲樣品緩沖*/</p><p>　　in[Ns], /*輸入樣本*/</p><p>　　d[Ns]; /*未知系統(tǒng)輸出緩沖區(qū)*/</p><p>　　unsigned int fir_index,sys_index;</p>

11、<p>　　void main()</p><p><b>　　{</b></p><p>　　int d,y,e,ue,j,i;</p><p>　　init(w,N0); /*初始化自適應(yīng)濾波器系數(shù)*/</p><p>　　init(s,N0); /*初始化自適應(yīng)濾波器的delay-line*/

12、</p><p>　　init(f,N1); /*未知的過(guò)濾delay-line初始化*/</p><p>　　fir_index=0; /*未知的過(guò)濾delay-line初始指標(biāo)*/</p><p>　　sys_index=0; /*初始化自適應(yīng)濾波技術(shù)delay-line指數(shù)*/</p><p>　　for(;;) /*產(chǎn)生兩

13、過(guò)濾器和樣品*/</p><p><b>　　{</b></p><p>　　random(in,Ns);</p><p>　　for(i=0;i<Ns;i++)</p><p><b>　　{</b></p><p>　　f[0]=in[i]; /*得到信號(hào)數(shù)據(jù)*/

14、</p><p><b>　　d=0;</b></p><p>　　for(j=0;j<N1;j++)</p><p>　　d +=_smpy(LP_coef[j],f[j]);/*計(jì)算d（n）*/</p><p>　　for(j=N1-1;j>0;j--)</p><p>　　f[j

15、] =f[j-1]; /*信號(hào)未知區(qū)中緩沖更新*/</p><p>　　s[0]=in[i]; /*將生成的偽隨機(jī)噪聲信號(hào)存入d</p><p>　　_sys[]數(shù)組中*/</p><p><b>　　y=0;</b></p><p>　　for(j=0;j<N0;j++)</p><p>

16、;　　y += _smpy(w[j],s[j]); /*計(jì)算LMS的輸出*/</p><p>　　e = d-y; /*計(jì)算誤差信號(hào)*/</p><p>　　ue = _smpy(TENMU,e);</p><p>　　for(j=0;j<N0;j++)</p><p>　　w[j] +=_smpy(u

17、e,s[j]); /*更新自適應(yīng)濾波器的系數(shù)*/</p><p>　　for(j=N0-1;j>0;j--)</p><p>　　s[j]=s[j-1]; /*更新信號(hào)緩沖*/</p><p><b>　　}</b></p><p><b>　　}</b></p>

18、<p><b>　　}</b></p><p>　　// init.c :初始化一個(gè)為零的數(shù)組</p><p>　　#pragma CODE_SECTION(init,"Ims_code");</p><p>　　void init(int *ptr,unsigned int N){</p>&

19、lt;p>　　unsigned int i;</p><p>　　for(i=N;i>0;i--)</p><p>　　*ptr++ =0;</p><p><b>　　}</b></p><p>　　// random.c:產(chǎn)生零均值隨機(jī)噪聲</p><p>　　#define U

20、SE_DATA 0</p><p>　　/**使用隨機(jī)數(shù)據(jù)文件作為輸入*/</p><p>　　#pragma CODE_SECTION(random,"Ims_code");</p><p>　　#pragma DATA_SECTION(randdata,"Ims_data");</p><p>　　

21、#pragma DATA_SECTION(i,"Ims_cinit");</p><p>　　#include "randdata.dat"</p><p>　　static unsigned int i=0;</p><p>　　void random(int* x,unsigned int N){</p>&

22、lt;p>　　unsigned int t;</p><p>　　for(t=N;t>0;t--){</p><p>　　*x++=randdata[i];</p><p><b>　　i++;</b></p><p>　　if(i==10000)</p><p><b>　

23、　i=0;</b></p><p><b>　　}</b></p><p><b>　　}</b></p><p>　　//LP_ coef.dat:</p><p>　　#pragma DATA_SECTION(LP_coef,"Ims_data");</p&

24、gt;<p>　　int LP_coef[48]={</p><p>　　-6.28,46,27,-35,-100,-93,26,91,240,52,-291,-497,-278,337,888,773,-210,-1486,-1895,</p><p>　　-442,2870,6793,9445,9445,6793,2870,-442,-1895,-1486,-210,77