dsp浮點(diǎn)實(shí)驗(yàn)二(編輯修改稿)
2024-09-26 15:10 本頁面
【文章內(nèi)容簡介】 j*S]； k＝ 0 to N－ l n＝ 0 這里： C＝ COS（ 2*PI*k*n/N）， S=sin（ 2*pi*k*n/N）， j=sqrt（ 1） */ ＃ include “ ” ＃ define N 64 .SECTION／ DM dm_data； ／ *變量定義在 DM數(shù)據(jù)段 *／ .VAR Input［ N］ =“ ”； ／ * 輸入數(shù)據(jù) *／ .VAR real[N] ＝“ ”； ／ *輸出數(shù)據(jù)的實(shí)部 *／ .VAR imagp[N]＝“ ”； ／ *輸出數(shù)據(jù)的虛部 *／ .SECTION／ PM pm_data； .VAR sine［ N］ =“ sin64． dat”； ／ *定義正弦旋轉(zhuǎn)因子表，放在 PM數(shù)據(jù)段 *／ ／ *余弦旋轉(zhuǎn)因子表由正弦表移動(dòng) N／ 4而得 *／ .VAR modul[N]= “” .SECTION／ PM pm_rsti； NOP； USTAT2＝ 0x108421； ／ *USTAT2為用戶自定義寄存器 *／ DM（ WAIT） =USTAT2； ／ *WAIT寄存器用來設(shè)定與外設(shè)數(shù)據(jù)傳 輸時(shí)的等待周期及狀態(tài)數(shù) *／ JUMP Start； .SECTION／ PM pm_code； /*程序代碼段 *／ START： M1=l； ／ *地址修改寄存器 */ M9=1； B0=input 。 ／ *設(shè)定循環(huán)基址 *／ L0=＠ input； ／ *輸入數(shù)據(jù)循環(huán)長度 *／ I1=imag； ／ *虛部輸出地址 *／ L1=0； CALL dft（ DB）； ／ *調(diào)用子程序，計(jì)算 DFT*／ I2=real 。 ／ *實(shí)部輸出地址 *／ L2=0； nop。 CALL modl。 /*調(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； ／ *通過 sine＋ N／ 4得到 cos放到 I9基址 *／ L9=＠ sine； ／ *循環(huán)長度． *／ I10=0； ／ *用來調(diào)整 M8的步長 *／ L10＝ 0； F15＝ 0； ／ *在每一外循環(huán)時(shí)，對 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步長 ，并讀取 X（ n）和 COS值 *／ F12=F0*F5， F4＝ PM（ I8， M8）； LCNTR=N－ l， DO inner UNTIL LCE。 ／ *內(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） 。 F13＝ F0*F4 , F9＝ F9＋ F12； ／ * F8中放虛部值， F9中放實(shí)部值 *／ F8=F8－ F13， DM（ I2， M1）＝ F9； ／ *存實(shí)部值 *／ MODIFY(I10, M9)； ／ *使旋轉(zhuǎn)因子表的尋址步長加一 *／ outer： DM（ I1， M1） =F8 ／ *存虛部值 *／ 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++，建立新工程，將位于將 DSP_exp\float\unit_2\DFT_MOD 目錄下的源文件（ ）添加到工程中，同時(shí)也將備用的數(shù)據(jù)文件 、 、 以及初始化文件 加入工程中。 2）選擇輸入數(shù)據(jù)文件 編譯程序前，在下述程序行中將輸入數(shù)據(jù)文件改為 ： .VAR input[N]= 。 其中 為方波信號。 3）編譯程序 分析、理解源程序，在此基礎(chǔ)上再編譯運(yùn)行程序?？赏ㄟ^點(diǎn)擊按鈕命令或選擇菜單命令 \Project\Build Project 來完成編譯。本程序的 sessions 為 ADSP21062 ADSP2106x Simulator。 4) 觀察結(jié)果 ? 程序調(diào)試時(shí)，可在程序中設(shè)置多個(gè)斷點(diǎn)，來觀察輸入數(shù)據(jù) x、輸出數(shù)據(jù)的實(shí)部 real 、輸出數(shù)據(jù)的虛部 imag和求模后的數(shù)據(jù) mudul的值； ? 打開 Mem
畢業(yè)設(shè)計(jì)相關(guān)推薦
鄂ICP備17016276號-1