dsp浮點(diǎn)實(shí)驗(yàn)二-免費(fèi)閱讀
【正文】 ? 數(shù)組 inputx 保持不變 ，將數(shù)組 inputy 分別改變?yōu)閿?shù)據(jù)文件 、（三角波）重復(fù) 3),4)，來(lái)考察序列的互相關(guān)功率譜。 .VAR inputy[M]= 。 F4=F4*F13,F2=dM(I0,M1)。 F6=F2。 RTS。 end: /******************* 計(jì)算功率譜的子程序 *************************/ psd: B0=output。 I0=inputx。 .SECTION/PM pm_rsti。 /* DM區(qū)變量定義 */ .VAR inputx[M]= 。 ? 數(shù)組 inputx 保持不變 （方波），將數(shù)組 inputy 分別改變?yōu)閿?shù)據(jù)文件、 square6 （三角波）重復(fù) 3),4)，來(lái)考察序列的互相關(guān)結(jié)果。 .VAR inputy[LENGTH_Y]= 。 /* 內(nèi)循環(huán)次數(shù)遞減 */ F9=PASS F15, F0=DM(I1,M1), F5=PM(I9,M9)。 F13=F0*F5。 M1=1。根據(jù)相關(guān)原理，兩個(gè)輸入序列長(zhǎng)度為 M=N=64，輸出序列長(zhǎng)度為 M+N1=127。 3. 典型實(shí)驗(yàn)結(jié)果 0 25 500 in p u t 0 25 500250 圖 三角波及其 DFT變換后的幅度譜 0 25 500 in p u t 0 25 500 圖 高斯白噪聲及其 DFT變換后的的幅度譜 4. 實(shí)驗(yàn)要求 1) 利用 plot功能，觀察輸入數(shù)據(jù) x、輸出數(shù)據(jù)的實(shí)部 real 、輸出數(shù)據(jù)的虛部 imag和求模后的 modul的波形； 2) 觀察各種波形的 DFT運(yùn)算和求模運(yùn)算的結(jié)果，比較并且說(shuō)明不同輸入數(shù)據(jù)對(duì)應(yīng)的DFT運(yùn)算輸出和求模運(yùn)算輸出的結(jié)果； 3) 在實(shí)驗(yàn)報(bào)告中記錄 square64的 DFT輸入和輸出數(shù)據(jù)及其波形，并 利用數(shù)字信號(hào)處理課程所學(xué)知識(shí)分析結(jié)果。 2）選擇輸入數(shù)據(jù)文件 編譯程序前，在下述程序行中將輸入數(shù)據(jù)文件改為 ： .VAR input[N]= 。 F4=RSQRTS F13, F7=F3。 F2=DM(I0,M1)。 /*調(diào)用子程序 ，計(jì)算 DFT的模值 *／ end： IDLE； /********* 計(jì)算 DFT子程序 ****************************/ dft： /*dft子程序 *／ B8＝ sine； ／ *將 sine放到 I8基址，并設(shè)定循環(huán)基址 *／ L8＝＠ sine； B9=sine； ／ *Cos表的循環(huán)基址 *／ I9＝ sine＋ N/4； ／ *通過(guò) sine＋ N／ 4得到 cos放到 I9基址 *／ L9=＠ sine； ／ *循環(huán)長(zhǎng)度． *／ I10=0； ／ *用來(lái)調(diào)整 M8的步長(zhǎng) *／ L10＝ 0； F15＝ 0； ／ *在每一外循環(huán)時(shí)，對(duì) F F9清 0*／ LCNTR＝ N， DO outer UNTIL LCE； ／ *外循環(huán) *／ F8＝ PASS F15， M8＝ I10； F9=PASS F15， F0＝ DM（ I0， M1）， F5=PM（ I9， M8）； ／ *上二句是清 F F9寄存器，設(shè)定 M8步長(zhǎng) ，并讀取 X（ n）和 COS值 *／ F12=F0*F5， F4＝ PM（ I8， M8）； LCNTR=N－ l， DO inner UNTIL LCE。在最后一個(gè)數(shù)據(jù)（數(shù)組中）被讀取后，地址寄存器中的地址會(huì)自動(dòng)變?yōu)閿?shù)組的起始地址。分別用來(lái)存放計(jì)算結(jié)果的實(shí)部和虛部。 通過(guò)本練習(xí)使學(xué)生加深對(duì)各種基本波形的 DFT 的理解， 加深對(duì)數(shù)字信號(hào)處理理論的理解。 執(zhí)行菜單命令 \Debug\Step Into 單步運(yùn)行程序，打開(kāi) Register File 和對(duì)應(yīng)的存儲(chǔ)區(qū)，觀察寄存器和存儲(chǔ)器的變化。 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ù)文件 、 、 、 、 以及初始化文件 加入工程中。 I1=inputy+LENGTH_Y1。 /* 指向 inputy 的結(jié)尾 */ F9=PASS F15, F0=DM(I1,M2), F5=PM(I9,M9)。 /* 等待 */ /********* 計(jì)算時(shí)域相關(guān)性的子程序 ****************************/ conv: I2=output。 L0=0。 .SECTION/PM pm_rsti。 卷積計(jì)算子程序中，將計(jì)算按照輸出長(zhǎng)度分成兩段完成，在程序中利用 left 和 right 兩個(gè)二重循環(huán)實(shí)現(xiàn)。實(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)處理理論的理解。根據(jù)卷積原理，兩個(gè)輸入序列長(zhǎng)度為 M=N=64，輸出序列長(zhǎng)度為 M+N1=127。 /* 定義數(shù)據(jù)存儲(chǔ)區(qū)中的變量 */ .VAR inputx[LENGTH_X]= 。 .SECTION/PM pm_code。 copydata:DM(I3,M1)=F0。 /* 用于給累加器清零 */ /* 計(jì)算前 LENGTH_X 個(gè)數(shù)據(jù) */ LCNTR=LENGTH_X, DO left UNTIL LCE。 /* 存儲(chǔ)計(jì)算結(jié)果 ,I2 遞增 */ /* 計(jì)算后面 LENGTH_Y1 個(gè)數(shù)據(jù) */ R1=LENGTH_Y。 inner2: F9=F9+F13, F0=DM(I1,M2),F5=PM(I9,M9)?？赏ㄟ^(guò) 點(diǎn)擊按鈕命令或選擇菜單命令 \Project\Build Project 來(lái)完成編譯。 3） 在實(shí)驗(yàn)報(bào)告中記錄 square1 和 square2 卷積的輸入和輸出數(shù)據(jù)及其波形 , 并 利用數(shù)字信號(hào)處理課程所學(xué)知識(shí)分析結(jié)果。由于 cos（ 2nkπ/N） = sin(2nk
畢業(yè)設(shè)計(jì)相關(guān)推薦
鄂ICP備17016276號(hào)-1