【導讀】本實驗主要通過四個練習來了解如何利用SIMULATOR實現(xiàn)基本的信號處理方法。兩個輸入序列的長度相同，分別寫入。程序編譯運行時將數(shù)據(jù)讀入存儲器中，并對其進行卷積運算，計算結(jié)果存放。于另一個緩沖區(qū)中。卷積的輸出指定到數(shù)組output中。程序主要由主程序和卷積子程序二個模塊組成。個二重循環(huán)實現(xiàn)。序列X采用正常序?qū)ぶ?，利用指令F5=DM讀取X，其中地址增量寄存器M9＝1，而序列Y是逆序?qū)ぶ?，利用指令。列Y的最后一個元素。啟動運行VisualDSP++，運行菜單命令\Project\New來新建一個工程，工程名稱自定義，、、以及初始化文件加入工程中。其中為單個沖擊脈沖信號，為正弦信號。分析、理解源程序，在此基礎(chǔ)上再編譯運行程序?？赏ㄟ^點擊按鈕命令或選擇菜單命令。本程序的sessions為ADSP-21062ADSP-2106xSimulator。執(zhí)行菜單命令\Debug\StepInto單步運行程序，打開RegisterFile和對應(yīng)的存儲區(qū)，觀。察寄存器和存儲器的變化。數(shù)組inputx保持為，將數(shù)組inputy分別改變?yōu)閿?shù)據(jù)文件(方波)、

　　

【正文】 uare6 （三角波）重復 3),4)，來考察序列的互相關(guān)結(jié)果。 3. 典型結(jié)果 0 25 500 in p u t 0 100 2000 圖 高斯白噪聲及其的自相關(guān) 0 25 500 in p u t 0 100 2000 圖 正弦信號及其自相關(guān) 3. 實驗要求 1） 利用 plot 功能，觀察各種輸入數(shù)據(jù)、輸出數(shù)據(jù)的的波形； 2） 按照實驗步驟完成實驗，觀察和比較自相關(guān)和互相關(guān)的計算結(jié)果，說明它們的差別。說明輸入序列和輸出序列的長度之間的關(guān)系； 3） 在實驗報告中記錄 trig64的 自相關(guān)、它與 sin64的互相關(guān)輸入和輸出數(shù)據(jù)及其波形、并 利用數(shù)字信號處理課程所學知識分析結(jié)果。 練習四 —— 利用 相關(guān)函數(shù)計算功率譜 1. 程序說明 本練習綜合了練習二和練習三的程序，參考程序為 。程序先對輸入信號進行相關(guān)運算，然后對相關(guān)運算的結(jié)果進行 DFT 運算。 DFT 運算后進行求模，所得的結(jié)果為該輸入信號的功率譜估計。 下面是練習四的主程序段和功率譜計算子程序，有關(guān)相關(guān)運算和 DFT子程序及其說明請參閱練習二和練習三的有關(guān)程序段。 include /* 存儲器映射的 IOP 寄存器定義 */ define N 64 /* 符號常數(shù) */ define M 32 .SECTION/DM dm_data。 /* DM區(qū)變量定義 */ .VAR inputx[M]= 。 .VAR inputy[M]= 。 .VAR output[N]=。 .VAR real[N]= 。 .VAR imag[N]= 。 .SECTION/PM pm_data。 /* PM 區(qū)變量定義 */ .VAR sine[N]= 。 /* 旋轉(zhuǎn)因子表 */ .VAR modul[N]= 。 .SECTION/PM pm_rsti。 /* 復位向量段 */ NOP。 USTAT2= 0x108421。 /* 復位后程序從這 里開始執(zhí)行 */ DM(WAIT)=USTAT2。 /* 設(shè)置外部存儲器訪問為 0等待 */ JUMP start。 .SECTION/PM pm_code。 start: M1=1。 M9=1。 I0=inputx。 L0=0。 I1=inputy。 L8=0。 CALL relate。 CALL psd。 CALL modl。 IDLE。 end: /******************* 計算功率譜的子程序 *************************/ psd: B0=output。 L0=@output。 /* 建立循環(huán)緩沖區(qū) */ I1=imag。 L1=0。 CALL dft (DB)。 /* 延遲調(diào)用 DFT子程序 */ I2=real。 /* 在被調(diào)用函數(shù)執(zhí)行前首先執(zhí)行 */ L2=0。 nop。 RTS。 /********* 求模計算子程序 ****************************/* modl: I0=real。 L0=0。 I1=imag。 L1=0。 I8=modul。 L8=0。 F2=DM(I0,M1)。 F6=F2。 F3=DM(I1,M1)。 F7=F3。 F8=F2*F6。 F13=F3*F7。 LCNTR=N, DO mod_cal UNTIL LCE。 F8=F2*F6,F13=F8+F13,F3=DM(I1,M1)。 F4=RSQRTS F13, F7=F3。 F4=F4*F13,F2=dM(I0,M1)。 F13=F3*F7,F6=F2。 mod_cal:PM(I8,M9)=F4。 nop。 nop。 RTS。 2. 實驗步驟 1) 調(diào)入程序 啟動運行 VisualDSP++，運行菜單命令 \Project\New來新建一個工程，工程名稱自定義，正確設(shè)置工程的各個選項，將 DSP_exp\float\unit_2\psd目錄下的源文件加入到工程中，同時加入數(shù)據(jù)文件 、 、 ，以及初始化文件 。 2）選擇輸入數(shù)據(jù)文件 編譯程序前，在下述程序行中將輸入數(shù)據(jù)文件改為 ： .VAR inputx[M]= 。 .VAR inputy[M]= 。 其中 為三角波信號。 3) 編譯程序 分析、理解源程序，在此基礎(chǔ)上再編譯運行程序?？赏ㄟ^點擊按鈕命令或選擇菜單命令 \Project\Build Project 來完成編譯。本程序的 sessions 為 ADSP21062 ADSP2106x Simulator。 4) 觀察結(jié)果 按 F5運行程序，然后執(zhí)行菜單命令 \View\Debug Windows\Plot 命令，分別作出輸入數(shù)據(jù)與輸出數(shù)據(jù) real[]、 imag[]和 modul[]的圖形。 執(zhí)行菜單命令 \Debug\Step Into 單步運行程序，打開 Register File 和對應(yīng)的存儲區(qū)，觀察寄存器和存儲器的變化。 5) 改變數(shù)據(jù)再運行 ? 將數(shù)組 inputx和數(shù)組 inputy同時改變?yōu)閿?shù)據(jù)文件 (方波 )、 （高斯白噪聲），重復 3),4)，來考察序列的功率譜。 ? 數(shù)組 inputx 保持不變 ，將數(shù)組 inputy 分別改變?yōu)閿?shù)據(jù)文件 、（三角波）重復 3),4)，來考察序列的互相關(guān)功率譜。 3. 典型結(jié)果 0 10 20 300 in p u t 0 25 500 . 0 * 1 002 . 5 * 1 0 4 圖 三角波的自相關(guān)功率譜 0 10 20 300 in p u t 0 25 500100 圖 高斯白噪聲的自相關(guān)功率譜 4. 實驗要求 1) 利用 plot 功能，觀察各種輸入數(shù)據(jù)、輸出數(shù)據(jù)的的波形； 2) 按照實驗步驟完成實驗，觀察和比較自相關(guān)譜和互相關(guān)譜的計算結(jié)果，說明它們的差別。 3) 在實驗報告中記錄 square32 的自相關(guān)譜、它與 trig32 的互相關(guān)譜的輸入和輸出數(shù)據(jù)及其波形、并 利用數(shù)字信號處理課程所學知識分析結(jié)果。