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