dsp浮點(diǎn)實(shí)驗(yàn)二-展示頁(yè)
2024-09-02 15:10本頁(yè)面
【正文】 的變量 */ .VAR inputx[LENGTH_X]= 。而序列 Y 是逆序?qū)ぶ罚弥噶頕0=DM(I1,M2)讀 取 Y(mk)，其中地址增量寄存器 M2＝－ 1，開(kāi)始計(jì)算時(shí)寄存器 I1 指向序列 Y 的最后一個(gè)元素。 Left 循環(huán)計(jì)算 m=0~N1， right 循環(huán)計(jì)算 m=N~N+M1。在主程序中，先將輸入數(shù)組 X 復(fù)制到輸出數(shù)組的后半部分，然后調(diào)用卷積子程序 conv 計(jì)算卷積。根據(jù)卷積原理，兩個(gè)輸入序列長(zhǎng)度為 M=N=64，輸出序列長(zhǎng)度為 M+N1=127。程序編 譯運(yùn)行時(shí)將數(shù)據(jù)讀入存儲(chǔ)器中，并對(duì)其進(jìn)行卷積運(yùn)算，計(jì)算結(jié)果存放于另一個(gè)緩沖區(qū)中。這是一個(gè)較簡(jiǎn)單的計(jì)算兩個(gè)輸入序列的卷積程序。 實(shí)驗(yàn)二 —— 用 SIMULATOR 模擬實(shí)現(xiàn)數(shù)字信號(hào)處理 本實(shí)驗(yàn)主要通過(guò)四個(gè)練習(xí)來(lái)了解如何利用 SIMULATOR實(shí)現(xiàn)基本的信號(hào)處理方法。實(shí)驗(yàn)包括以下幾個(gè)練習(xí)： 練習(xí)一：時(shí)域卷積運(yùn)算 練習(xí)二： DFT運(yùn)算 練習(xí)三：時(shí)域相關(guān)運(yùn)算 練習(xí)四：利用相關(guān)函數(shù)計(jì)算信號(hào)的功率譜 練習(xí)一 —— 時(shí)域卷積運(yùn)算 練習(xí)一的目的是使學(xué)生掌握 VisualDSP++的使用，同時(shí)加深對(duì)數(shù)字信號(hào)處理理論的理解。兩個(gè)輸入序列的長(zhǎng)度相同，分別寫(xiě)入到幾個(gè)文件中。 1． 程序說(shuō)明 設(shè)輸入序列 X的長(zhǎng)度為 M， Y的長(zhǎng)度為 N，則程序按照因果系統(tǒng)計(jì)算卷積： ??? ?? 10 )()()( Mk kmYkXmZ , 10 ???? MNm 程序中一個(gè)輸入序列 X 指定到數(shù)組 inputx中，另一個(gè)輸入序列 Y指定到數(shù)組 inputy 中，卷積的輸出指定到數(shù)組 output 中。 程序主要由主程序和卷積子 程序（位于 pm_code 段中）二個(gè)模塊組成。 卷積計(jì)算子程序中，將計(jì)算按照輸出長(zhǎng)度分成兩段完成，在程序中利用 left 和 right 兩個(gè)二重循環(huán)實(shí)現(xiàn)。在計(jì)算卷積時(shí)，序列 X 采用正常序?qū)ぶ罚弥噶?F5=DM(I9,M9)讀取 X(k)，其中地址增量寄存器 M9＝ 1，開(kāi)始計(jì)算時(shí)寄存器 I9 指向序列 X 的第一個(gè)元素。 實(shí)驗(yàn)測(cè)試程序主要模塊的源程序如下： /******************************************************************** ADSP2106x 程序功能：計(jì)算 X[]與 Y[]的卷積 ******************************************************************/ include /* 存儲(chǔ)器映射的 IOP 寄存器定義 */ define LENGTH_X 64 /* 符號(hào)常數(shù) */ define LENGTH_Y 64 .SECTION/DM dm_data。 .VAR inputy[LENGTH_Y]= 。 .SECTION/PM pm_rsti。 USTAT2= 0x108421。 /* 設(shè)置外部存儲(chǔ)器訪問(wèn)為 0 等待 */ JUMP start。 start: M1=1。 L0=0。 L3=0。 F0=DM(I0,M1)。 CALL conv。 /* 等待 */ /********* 計(jì)算時(shí)域相關(guān)性的子程序 ****************************/ conv: I2=output。 M2=1。 F15=0。 I9=I2。 /* 指向 inputy 的結(jié)尾 */ F9=PASS F15, F0=DM(I1,M2), F5=PM(I9,M9)。 F13=F0*F5。 /* F0=數(shù)組 inputy , F5=數(shù)組 inputx */ left: DM(I2,M1)=F9。 /* 內(nèi)循次數(shù)計(jì)數(shù)器 */ LCNTR=LENGTH_Y1, DO right UNTIL LCE。 I1=inputy+LENGTH_Y1。 /* 循環(huán)次數(shù)遞減 */ F9=PASS F15, F0=DM(I1,M2), F5=PM(I9,M9)。 F13=F0*F5。 /* F0=數(shù)組 inputy , F5=數(shù)組 inputx */ right: DM(I2,M1)=F9。 2． 實(shí)驗(yàn)步驟 1）調(diào)入程序 啟動(dòng)運(yùn)行 VisualDSP++，運(yùn)行菜單命令 \Project\New來(lái)新建一個(gè)工程，工程名稱(chēng)自定義，正確設(shè)置工程的各個(gè)選項(xiàng)，將 DSP_exp\float\unit_2\Conv 目錄下的源文件（ 和）加入到工程中，同時(shí)也將備用的數(shù)據(jù)文件 、 、 、 、 以及初始化文件 加入工程中。 .VAR inputy[LENGTH_Y]= 。 3）編譯程序 分析、理解源程序，在此基礎(chǔ)上再編譯運(yùn)行程序。本程序的 sessions 為 ADSP21062 ADSP2106x Simulator。 執(zhí)行菜單命令 \Debug\Step Into 單步運(yùn)行程序，打開(kāi) Register File 和對(duì)應(yīng)的存儲(chǔ)區(qū)，觀察寄存器和存儲(chǔ)器的變化。 ? 數(shù)組 inputx 改變?yōu)?（ 4 個(gè)沖擊脈沖），將數(shù)組 inputy 分別改變?yōu)閿?shù)據(jù)文件， (方波 )、 （三角波）重復(fù) 3),4)，來(lái)考察程序的執(zhí)行結(jié)果。 3. 典型結(jié)果 0 25 500 in p u t 0 100 2000 圖 正弦波與沖擊脈沖的卷積 0 25 50010 in p u t 0 100 20002500 圖 兩個(gè)不同寬度的矩形脈沖的卷積 0 25 50010 in p u t 0 100 2000250 圖 沖擊脈沖序列與矩形脈沖的卷積 4．實(shí)驗(yàn)要求 1） 利用 plot 功能，觀察各種輸入數(shù)據(jù)、輸出數(shù)據(jù)的的波形； 2） 按照實(shí)驗(yàn)步驟完成實(shí)驗(yàn)，比較給定的幾種數(shù)據(jù)的卷積結(jié)果，說(shuō)明幾種信號(hào)卷積的差別。 練習(xí)二 — — DFT 計(jì)算 離散傅里葉變換 (DFT)的計(jì)算公式為： ?????? ??? ?? ????? NnkjNnknxenxkX NnNnkjNn ??? 2s i n2c os)()()( 10210 其中 x(n)是輸入的數(shù)字序列， n=0?N1； X(k)是輸入數(shù)字序列的 DFT 系數(shù)， k=0?N1。 通過(guò)本練習(xí)使學(xué)生加深對(duì)各種基本波形的 DFT 的理解， 加深對(duì)數(shù)字信號(hào)處理理論的理解。 本練習(xí)采用 64 點(diǎn)離散傅里葉變換， DFT 程序流圖如圖 所示。 2） 正弦旋轉(zhuǎn)因子 sin(2nkπ/N)的值預(yù)先計(jì)算，并放在 PM數(shù)據(jù)段的 sin[N]數(shù)組中
畢業(yè)設(shè)計(jì)相關(guān)推薦
鄂ICP備17016276號(hào)-1