dsp浮點(diǎn)實(shí)驗(yàn)二-資料下載頁(yè)
2025-08-12 15:10本頁(yè)面
【導(dǎo)讀】本實(shí)驗(yàn)主要通過(guò)四個(gè)練習(xí)來(lái)了解如何利用SIMULATOR實(shí)現(xiàn)基本的信號(hào)處理方法。兩個(gè)輸入序列的長(zhǎng)度相同，分別寫(xiě)入。程序編譯運(yùn)行時(shí)將數(shù)據(jù)讀入存儲(chǔ)器中，并對(duì)其進(jìn)行卷積運(yùn)算，計(jì)算結(jié)果存放。于另一個(gè)緩沖區(qū)中。卷積的輸出指定到數(shù)組output中。程序主要由主程序和卷積子程序二個(gè)模塊組成。個(gè)二重循環(huán)實(shí)現(xiàn)。序列X采用正常序?qū)ぶ?，利用指令F5=DM讀取X，其中地址增量寄存器M9＝1，而序列Y是逆序?qū)ぶ?，利用指令。列Y的最后一個(gè)元素。啟動(dòng)運(yùn)行VisualDSP++，運(yùn)行菜單命令\Project\New來(lái)新建一個(gè)工程，工程名稱(chēng)自定義，、、以及初始化文件加入工程中。其中為單個(gè)沖擊脈沖信號(hào)，為正弦信號(hào)。分析、理解源程序，在此基礎(chǔ)上再編譯運(yùn)行程序?？赏ㄟ^(guò)點(diǎn)擊按鈕命令或選擇菜單命令。本程序的sessions為ADSP-21062ADSP-2106xSimulator。執(zhí)行菜單命令\Debug\StepInto單步運(yùn)行程序，打開(kāi)RegisterFile和對(duì)應(yīng)的存儲(chǔ)區(qū)，觀(guān)。察寄存器和存儲(chǔ)器的變化。數(shù)組inputx保持為，將數(shù)組inputy分別改變?yōu)閿?shù)據(jù)文件(方波)、
【正文】 uare6 （三角波）重復(fù) 3),4)，來(lái)考察序列的互相關(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 圖 正弦信號(hào)及其自相關(guān) 3. 實(shí)驗(yàn)要求 1） 利用 plot 功能，觀(guān)察各種輸入數(shù)據(jù)、輸出數(shù)據(jù)的的波形； 2） 按照實(shí)驗(yàn)步驟完成實(shí)驗(yàn)，觀(guān)察和比較自相關(guān)和互相關(guān)的計(jì)算結(jié)果，說(shuō)明它們的差別。說(shuō)明輸入序列和輸出序列的長(zhǎng)度之間的關(guān)系； 3） 在實(shí)驗(yàn)報(bào)告中記錄 trig64的 自相關(guān)、它與 sin64的互相關(guān)輸入和輸出數(shù)據(jù)及其波形、并 利用數(shù)字信號(hào)處理課程所學(xué)知識(shí)分析結(jié)果。 練習(xí)四 —— 利用 相關(guān)函數(shù)計(jì)算功率譜 1. 程序說(shuō)明 本練習(xí)綜合了練習(xí)二和練習(xí)三的程序，參考程序?yàn)?。程序先對(duì)輸入信號(hào)進(jìn)行相關(guān)運(yùn)算，然后對(duì)相關(guān)運(yùn)算的結(jié)果進(jìn)行 DFT 運(yùn)算。 DFT 運(yùn)算后進(jìn)行求模，所得的結(jié)果為該輸入信號(hào)的功率譜估計(jì)。 下面是練習(xí)四的主程序段和功率譜計(jì)算子程序，有關(guān)相關(guān)運(yùn)算和 DFT子程序及其說(shuō)明請(qǐng)參閱練習(xí)二和練習(xí)三的有關(guān)程序段。 include /* 存儲(chǔ)器映射的 IOP 寄存器定義 */ define N 64 /* 符號(hào)常數(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。 /* 復(fù)位向量段 */ NOP。 USTAT2= 0x108421。 /* 復(fù)位后程序從這 里開(kāi)始執(zhí)行 */ DM(WAIT)=USTAT2。 /* 設(shè)置外部存儲(chǔ)器訪(fǎng)問(wèn)為 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: /******************* 計(jì)算功率譜的子程序 *************************/ 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。 /********* 求模計(jì)算子程序 ****************************/* 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. 實(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\psd目錄下的源文件加入到工程中，同時(shí)加入數(shù)據(jù)文件 、 、 ，以及初始化文件 。 2）選擇輸入數(shù)據(jù)文件 編譯程序前，在下述程序行中將輸入數(shù)據(jù)文件改為 ： .VAR inputx[M]= 。 .VAR inputy[M]= 。 其中 為三角波信號(hào)。 3) 編譯程序 分析、理解源程序，在此基礎(chǔ)上再編譯運(yùn)行程序?？赏ㄟ^(guò)點(diǎn)擊按鈕命令或選擇菜單命令 \Project\Build Project 來(lái)完成編譯。本程序的 sessions 為 ADSP21062 ADSP2106x Simulator。 4) 觀(guān)察結(jié)果 按 F5運(yùn)行程序，然后執(zhí)行菜單命令 \View\Debug Windows\Plot 命令，分別作出輸入數(shù)據(jù)與輸出數(shù)據(jù) real[]、 imag[]和 modul[]的圖形。 執(zhí)行菜單命令 \Debug\Step Into 單步運(yùn)行程序，打開(kāi) Register File 和對(duì)應(yīng)的存儲(chǔ)區(qū)，觀(guān)察寄存器和存儲(chǔ)器的變化。 5) 改變數(shù)據(jù)再運(yùn)行 ? 將數(shù)組 inputx和數(shù)組 inputy同時(shí)改變?yōu)閿?shù)據(jù)文件 (方波 )、 （高斯白噪聲），重復(fù) 3),4)，來(lái)考察序列的功率譜。 ? 數(shù)組 inputx 保持不變 ，將數(shù)組 inputy 分別改變?yōu)閿?shù)據(jù)文件 、（三角波）重復(fù) 3),4)，來(lái)考察序列的互相關(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. 實(shí)驗(yàn)要求 1) 利用 plot 功能，觀(guān)察各種輸入數(shù)據(jù)、輸出數(shù)據(jù)的的波形； 2) 按照實(shí)驗(yàn)步驟完成實(shí)驗(yàn)，觀(guān)察和比較自相關(guān)譜和互相關(guān)譜的計(jì)算結(jié)果，說(shuō)明它們的差別。 3) 在實(shí)驗(yàn)報(bào)告中記錄 square32 的自相關(guān)譜、它與 trig32 的互相關(guān)譜的輸入和輸出數(shù)據(jù)及其波形、并 利用數(shù)字信號(hào)處理課程所學(xué)知識(shí)分析結(jié)果。
畢業(yè)設(shè)計(jì)相關(guān)推薦
鄂ICP備17016276號(hào)-1