dsp浮點(diǎn)實(shí)驗(yàn)二(完整版)
【正文】 3=DM(I1,M1)。 2. 實(shí)驗(yàn)步驟 1) 調(diào)入程序 啟動(dòng)運(yùn)行 VisualDSP++，建立新工程，將位于將 DSP_exp\float\unit_2\DFT_MOD 目錄下的源文件（ ）添加到工程中，同時(shí)也將備用的數(shù)據(jù)文件 、 、 以及初始化文件 加入工程中。 5）改變數(shù)據(jù)再運(yùn)行 ? 數(shù)組 input 分別改變?yōu)閿?shù)據(jù)文件 （三角波）、 (高斯白噪聲 )，重復(fù)3),4)，來(lái)考察程序的執(zhí)行結(jié)果。 程序中一個(gè)輸入序列 X指定到數(shù)組 inputx中，另一個(gè)輸入序列 Y指定到數(shù)組 inputy中，相關(guān)的輸出指定到數(shù)組 output 中。 下面是練習(xí)三的計(jì)算子程序： /******************** 計(jì)算時(shí)域相關(guān)性的子程序 ********************/ relate: I2=output。 LCNTR=LENGTH_Y, DO inner1 UNTIL LCE。 /* 指向 inputy */ R1=R11。 2）選擇輸入數(shù)據(jù)文件 編譯程序前，在下述程序行中將輸入數(shù)據(jù)文件改為 ： .VAR inputx[LENGTH_X]= 。 5) 改變數(shù)據(jù)再運(yùn)行 ? 將數(shù)組 inputx和數(shù)組 inputy同時(shí)改變?yōu)閿?shù)據(jù)文件 (正弦波 )、 square64(方波 )、（三角波） （高斯白噪聲），重復(fù) 3),4)，來(lái)考察序列的自相關(guān)結(jié)果。 include /* 存儲(chǔ)器映射的 IOP 寄存器定義 */ define N 64 /* 符號(hào)常數(shù) */ define M 32 .SECTION/DM dm_data。 /* 旋轉(zhuǎn)因子表 */ .VAR modul[N]= 。 M9=1。 IDLE。 nop。 F2=DM(I0,M1)。 F4=RSQRTS F13, F7=F3。 2）選擇輸入數(shù)據(jù)文件 編譯程序前，在下述程序行中將輸入數(shù)據(jù)文件改為 ： .VAR inputx[M]= 。 5) 改變數(shù)據(jù)再運(yùn)行 ? 將數(shù)組 inputx和數(shù)組 inputy同時(shí)改變?yōu)閿?shù)據(jù)文件 (方波 )、 （高斯白噪聲），重復(fù) 3),4)，來(lái)考察序列的功率譜。 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 功能，觀察各種輸入數(shù)據(jù)、輸出數(shù)據(jù)的的波形； 2) 按照實(shí)驗(yàn)步驟完成實(shí)驗(yàn)，觀察和比較自相關(guān)譜和互相關(guān)譜的計(jì)算結(jié)果，說(shuō)明它們的差別。 其中 為三角波信號(hào)。 F13=F3*F7,F6=F2。 F3=DM(I1,M1)。 /********* 求模計(jì)算子程序 ****************************/* modl: I0=real。 L0=output。 L0=0。 /* 復(fù)位向量段 */ NOP。 .VAR inputy[M]= 。 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 功能，觀察各種輸入數(shù)據(jù)、輸出數(shù)據(jù)的的波形； 2） 按照實(shí)驗(yàn)步驟完成實(shí)驗(yàn)，觀察和比較自相關(guān)和互相關(guān)的計(jì)算結(jié)果，說(shuō)明它們的差別。 其中 為方波信號(hào)。 LCNTR=R1, DO inner2 UNTIL LCE。 inner1: F9=F9+F13, F0=DM(I1,M1), F5=PM(I9,M9)。 M9=1。 程序主要由主程序和相關(guān)子程序（位于 pm_code 段中）二個(gè)模塊組成。 練習(xí)三 —— 時(shí)域 相關(guān)計(jì)算 1. 程序說(shuō)明 本練習(xí)三完成序列的自相關(guān)計(jì)算。 其中 為方波信號(hào)。 F4=F4*F13,F2=DM(I0,M1)。 F6=F2。 ／ *內(nèi)循環(huán) *／ F13＝ F0*F4， F9＝ F9+F12， F0＝ DM（ I0， M1）， F5=PM（ I9， M8）； inner： F12=F0*F5， F8=F8－ F13， F4＝ PM（ I8， M8） 。 數(shù)據(jù)初始化執(zhí)行 N 次？結(jié)束F 13= X ( N ) S I N E( N )F 8= F 8 F 13F 9= F 9+ F 12旋轉(zhuǎn)因子的 M8=M8+ 1存 F9 、 F8 值F 8= 0 ； F 9= 0F 12= X ( 0) C O S ( 0)執(zhí)行 N 1 次 ?F 13= X ( n ) S I N E( n )F 9= F 9+ F 12F 12= X ( n + 1) C O S ( n + 1)F 8= F 8 F 13n = n + 1yesnoyesno 圖 DFT程序流圖 6）程序中因 sin（ 2nkπ/N）中的 k、 n都是變量，所以在內(nèi)循環(huán)中 n從 0到 N變化，即對(duì)每一個(gè) k值均有 N次相乘，將其相加，得到一個(gè)數(shù)的 DFT變換值。 I I2中分別存放輸出數(shù)據(jù) imag和 real的地址。輸入數(shù)據(jù)放在幾個(gè)文件中，程序編譯運(yùn)行時(shí)將數(shù)據(jù)讀入存儲(chǔ)器中，并對(duì)其進(jìn)行 DFT 運(yùn)算，計(jì)算結(jié)果存放于輸出緩沖區(qū)中。 6） 改變數(shù)據(jù)再運(yùn)行 ? 數(shù)組 inputx保持為 ，將數(shù)組 inputy 分別改變?yōu)閿?shù)據(jù)文件 (方波 )、（三角波）重復(fù) 3),4)，來(lái)考察程序的執(zhí)行結(jié)果。 2）選擇輸入數(shù)據(jù)文件 編譯程序前，在下述程序行中將輸入數(shù)據(jù)文件改為 和 ： .VAR inputx[LENGTH_X]= 。 /* 指向 inputy 的結(jié)尾 */ R1=R11。 LCNTR=LENGTH_Y, DO inner1 UNTIL LCE。 M1=1。 I3=output+LENGTH_Y1。 /* 復(fù)位向量段 */ NOP。 Left 循環(huán)計(jì)算 m=0~N1， right 循環(huán)計(jì)算 m=N~N+M1。這是一個(gè)較簡(jiǎn)單的計(jì)算兩個(gè)輸入序列的卷積程序。 1． 程序說(shuō)明 設(shè)輸入序列 X的長(zhǎng)度為 M， Y的長(zhǎng)度為 N，則程序按照因果系統(tǒng)計(jì)
畢業(yè)設(shè)計(jì)相關(guān)推薦
鄂ICP備17016276號(hào)-1