【正文】 ? 程序調(diào)試時，可在程序中設(shè)置多個斷點，來觀察輸入數(shù)據(jù) x、輸出數(shù)據(jù)的實部 real 、輸出數(shù)據(jù)的虛部 imag和求模后的數(shù)據(jù) mudul的值； ? 打開 Memory\Two Column 窗口，在窗口中用鼠標右鍵點擊，在彈出的菜單中選擇用浮點數(shù)（ Floating Point 32bit）格式顯示數(shù)據(jù)?？赏ㄟ^點擊按鈕命令或選擇菜單命令 \Project\Build Project 來完成編譯。 其中 為方波信號。 2. 實驗步驟 1) 調(diào)入程序 啟動運行 VisualDSP++，建立新工程，將位于將 DSP_exp\float\unit_2\DFT_MOD 目錄下的源文件（ ）添加到工程中，同時也將備用的數(shù)據(jù)文件 、 、 以及初始化文件 加入工程中。 NOP。 mod_cal: PM(I8,M9)=F4。 F4=F4*F13,F2=DM(I0,M1)。 F8=F2*F6,F13=F8+F13,F3=DM(I1,M1)。 F13=F3*F7。 F7=F3。 F6=F2。 L8=0。 L1=0。 L0=0。 ／ *內(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） 。 CALL modl。 ／ *設(shè)定循環(huán)基址 *／ L0=＠ input； ／ *輸入數(shù)據(jù)循環(huán)長度 *／ I1=imag； ／ *虛部輸出地址 *／ L1=0； CALL dft（ DB）； ／ *調(diào)用子程序，計算 DFT*／ I2=real 。例子中在每一內(nèi)循環(huán)后，將 I10＋ l，即使 M8值加 1。 數(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變化，即對每一個 k值均有 N次相乘，將其相加，得到一個數(shù)的 DFT變換值。故在內(nèi)循環(huán)中，無需將數(shù)組的起始地址復位。 DFT 數(shù)據(jù)處理中采用雙重循環(huán)來實現(xiàn)。 即： Call dft（ DB）； 采用遲延跳轉(zhuǎn)方式，即該指令后的兩條指令在子程序正式執(zhí)行前被調(diào)用執(zhí)行，來提高指令的執(zhí)行效率。 I I2中分別存放輸出數(shù)據(jù) imag和 real的地址。 3）在 DM數(shù)據(jù)段定義了 兩數(shù)組： real[N]和 imag[N]。由于 cos（ 2nkπ/N） = sin(2nkπ/N+π/2)，所以 cos［ N］＝ sine［ N＋ N/4］即可由正弦值數(shù)組中跳過 N／ 4個數(shù)來獲得其余弦值。 1. 程序說明 1） 程序中，輸入的 64點數(shù)據(jù)放在數(shù)據(jù)文件 、 、 ，并賦給DM段中的 input[N]數(shù)組。輸入數(shù)據(jù)放在幾個文件中，程序編譯運行時將數(shù)據(jù)讀入存儲器中，并對其進行 DFT 運算，計算結(jié)果存放于輸出緩沖區(qū)中。它僅存在于離散點，因此是離散譜。 3） 在實驗報告中記錄 square1 和 square2 卷積的輸入和輸出數(shù)據(jù)及其波形 , 并 利用數(shù)字信號處理課程所學知識分析結(jié)果。 ? 數(shù)組 inputx改變?yōu)?，將數(shù)組 inputy 改變?yōu)閿?shù)據(jù)文件 (方波 )， 重復 3),4)，來考察程序的執(zhí)行結(jié)果。 6） 改變數(shù)據(jù)再運行 ? 數(shù)組 inputx保持為 ，將數(shù)組 inputy 分別改變?yōu)閿?shù)據(jù)文件 (方波 )、（三角波）重復 3),4)，來考察程序的執(zhí)行結(jié)果。 4）觀察結(jié)果 按 F5運行程序，然后執(zhí)行菜單命令 \View\Debug Windows\Plot 命令，分別作出輸入數(shù)據(jù) inputx[]和 inputy[] 與輸出數(shù)據(jù) output[]的圖形?？赏ㄟ^ 點擊按鈕命令或選擇菜單命令 \Project\Build Project 來完成編譯。 其中 為單個沖擊脈沖信號， 為正弦信號。 2）選擇輸入數(shù)據(jù)文件 編譯程序前，在下述程序行中將輸入數(shù)據(jù)文件改為 和 ： .VAR inputx[LENGTH_X]= 。 /* 存儲計算結(jié)果 ,I2 遞增 */ RTS。 inner2: F9=F9+F13, F0=DM(I1,M2),F5=PM(I9,M9)。 LCNTR=R1, DO inner2 UNTIL LCE。 /* 指向 inputy 的結(jié)尾 */ R1=R11。 I9=I2。 /* 存儲計算結(jié)果 ,I2 遞增 */ /* 計算后面 LENGTH_Y1 個數(shù)據(jù) */ R1=LENGTH_Y。 inner1: F9=F9+F13, F0=DM(I1,M2), F5=PM(I9,M9)。 LCNTR=LENGTH_Y, DO inner1 UNTIL LCE。 I1=inputy+LENGTH_Y1。 /* 用于給累加器清零 */ /* 計算前 LENGTH_X 個數(shù)據(jù) */ LCNTR=LENGTH_X, DO left UNTIL LCE。 M9=1。 M1=1。 /* 調(diào)用子程序 ,進行卷積計算 */ end: IDLE。 copydata:DM(I3,M1)=F0。 /* 把數(shù)組 inputx復制到 output 的后半段（位于 LENGTH_Y－ 1 開始） */ LCNTR=LENGTH_X, do copydata until lce。 I3=output+LENGTH_Y1。 I0=inputx。 .SECTION/PM pm_code。 /* 復位后程序從這里開始執(zhí)行 */ DM(WAIT)=USTAT2。 /* 復位向量段 */ NOP。 .VAR output[LENGTH_X+LENGTH_Y1]= 。 /* 定義數(shù)據(jù)存儲區(qū)中