dsp浮點實驗二-wenkub
2022-09-01 15:10:32 本頁面
【正文】 I9 指向序列 X 的第一個元素。 程序主要由主程序和卷積子 程序（位于 pm_code 段中）二個模塊組成。兩個輸入序列的長度相同，分別寫入到幾個文件中。 實驗二 —— 用 SIMULATOR 模擬實現(xiàn)數(shù)字信號處理 本實驗主要通過四個練習來了解如何利用 SIMULATOR實現(xiàn)基本的信號處理方法。程序編 譯運行時將數(shù)據(jù)讀入存儲器中，并對其進行卷積運算，計算結(jié)果存放于另一個緩沖區(qū)中。在主程序中，先將輸入數(shù)組 X 復(fù)制到輸出數(shù)組的后半部分，然后調(diào)用卷積子程序 conv 計算卷積。而序列 Y 是逆序?qū)ぶ罚弥噶頕0=DM(I1,M2)讀 取 Y(mk)，其中地址增量寄存器 M2＝－ 1，開始計算時寄存器 I1 指向序列 Y 的最后一個元素。 .VAR output[LENGTH_X+LENGTH_Y1]= 。 /* 復(fù)位后程序從這里開始執(zhí)行 */ DM(WAIT)=USTAT2。 I0=inputx。 /* 把數(shù)組 inputx復(fù)制到 output 的后半段（位于 LENGTH_Y－ 1 開始） */ LCNTR=LENGTH_X, do copydata until lce。 /* 調(diào)用子程序 ,進行卷積計算 */ end: IDLE。 M9=1。 I1=inputy+LENGTH_Y1。 inner1: F9=F9+F13, F0=DM(I1,M2), F5=PM(I9,M9)。 I9=I2。 LCNTR=R1, DO inner2 UNTIL LCE。 /* 存儲計算結(jié)果 ,I2 遞增 */ RTS。 其中 為單個沖擊脈沖信號， 為正弦信號。 4）觀察結(jié)果 按 F5運行程序，然后執(zhí)行菜單命令 \View\Debug Windows\Plot 命令，分別作出輸入數(shù)據(jù) inputx[]和 inputy[] 與輸出數(shù)據(jù) output[]的圖形。 ? 數(shù)組 inputx改變?yōu)?，將數(shù)組 inputy 改變?yōu)閿?shù)據(jù)文件 (方波 )， 重復(fù) 3),4)，來考察程序的執(zhí)行結(jié)果。它僅存在于離散點，因此是離散譜。 1. 程序說明 1） 程序中，輸入的 64點數(shù)據(jù)放在數(shù)據(jù)文件 、 、 ，并賦給DM段中的 input[N]數(shù)組。 3）在 DM數(shù)據(jù)段定義了 兩數(shù)組： real[N]和 imag[N]。 即： Call dft（ DB）； 采用遲延跳轉(zhuǎn)方式，即該指令后的兩條指令在子程序正式執(zhí)行前被調(diào)用執(zhí)行，來提高指令的執(zhí)行效率。故在內(nèi)循環(huán)中，無需將數(shù)組的起始地址復(fù)位。例子中在每一內(nèi)循環(huán)后，將 I10＋ l，即使 M8值加 1。 CALL modl。 L0=0。 L8=0。 F7=F3。 F8=F2*F6,F13=F8+F13,F3=DM(I1,M1)。 mod_cal: PM(I8,M9)=F4。 2. 實驗步驟 1) 調(diào)入程序 啟動運行 VisualDSP++，建立新工程，將位于將 DSP_exp\float\unit_2\DFT_MOD 目錄下的源文件（ ）添加到工程中，同時也將備用的數(shù)據(jù)文件 、 、 以及初始化文件 加入工程中?？赏ㄟ^點擊按鈕命令或選擇菜單命令 \Project\Build Project 來完成編譯。 5）改變數(shù)據(jù)再運行 ? 數(shù)組 input 分別改變?yōu)閿?shù)據(jù)文件 （三角波）、 (高斯白噪聲 )，重復(fù)3),4)，來考察程序的執(zhí)行結(jié)果。這是一個較簡單的計算輸入序列 的相關(guān)程序。 程序中一個輸入序列 X指定到數(shù)組 inputx中，另一個輸入序列 Y指定到數(shù)組 inputy中，相關(guān)的輸出指定到數(shù)組 output 中。 卷積計算子程序中，將計算按照輸出長度分成兩段完成，在程序中利用 left 和 right 兩個二重循環(huán)實現(xiàn)。 下面是練習三的計算子程序： /******************** 計算時域相關(guān)性的子程序 ********************/ relate: I2=output。 /* 用于給累加器清零 */ /* 下面的 left循環(huán)計算前面 LENGTH_X 個相關(guān)數(shù)據(jù) */ LCNTR=LENGTH_X, DO left UNTIL LCE。 LCNTR=LENGTH_Y, DO inner1 UNTIL LCE。 /* 存儲計算結(jié)果 ,I2遞增 */ /* 下面的 right循環(huán)計算后面 LENGTH_Y1 個相關(guān)數(shù)據(jù) */ R1=LENGTH_Y。 /* 指向 inputy */ R1=R11。 inner2: F9=F9+F13, F0=DM(I1,M1),F5=PM(I9,M9)。 2）選擇輸入數(shù)據(jù)文件 編譯程序前，在下述程序行中將輸入數(shù)據(jù)文件改為 ： .VAR inputx[LENGTH_X]= ?？赏ㄟ^點擊按鈕命令或選擇菜單命令 \Project\Build Project 來完成編譯。 5) 改變數(shù)據(jù)再運行 ? 將數(shù)組 inputx和數(shù)組 inputy同時改變?yōu)閿?shù)據(jù)文件 (正弦波 )、 square64(方波 )、（三角波） （高斯白噪聲），重復(fù) 3),4)，來考察序列的自相關(guān)結(jié)果。 練習四 —— 利用 相關(guān)函數(shù)計算功率譜 1. 程序說明 本練習綜合了練習二和練習三的程序，參考程序為 。 include /* 存儲器映射的 IOP 寄存器定義 */ define N 64 /* 符號常數(shù) */ define M 32 .SECTION/DM dm_data。 .VAR real[N]= 。 /* 旋轉(zhuǎn)因子表 */ .VAR modul[N]= 。 /* 復(fù)位后程序從這 里開始執(zhí)行 */ DM(WAIT)=USTAT2。 M9=1。 L8=0。 IDLE。 L1=0。 nop。 I1=imag。 F2=DM(I0,M1)。 F8=F2*F6。 F4=RSQRTS F13, F7=F3。 nop。 2）選擇輸入數(shù)據(jù)文件 編譯程序前，在下述程序行中將輸入數(shù)據(jù)文件改為 ： .VAR inputx[M]= ?？赏ㄟ^點擊按鈕命令或選擇菜單命令 \Project\Build Project 來完成編譯。 5) 改變數(shù)據(jù)再運行 ? 將數(shù)組 inputx和數(shù)組 inputy同時改變?yōu)閿?shù)據(jù)文件 (方波 )、 （高斯白噪聲），重復(fù) 3),4)，來考察序列的功率譜。 。 3. 典型結(jié)果 0 10 20 300 in p
畢業(yè)設(shè)計相關(guān)推薦
鄂ICP備17016276號-1