【導(dǎo)讀】，先彈出數(shù)據(jù)后，再SP+1。，必須先設(shè)置，后使用。例5-1設(shè)計(jì)一存儲(chǔ)空間為100個(gè)單元的堆棧。例5-3寫求解直線方程的程序，計(jì)算y=mx+b。MPY@x,A；將x地址的內(nèi)容與T中的m相乘，傳送速度比加載和存儲(chǔ)指令要快；傳送數(shù)據(jù)不需要通過(guò)累加器；可以尋址程序存儲(chǔ)器；pmad為16位立即數(shù)程序存儲(chǔ)器地址；Ymem從CB數(shù)據(jù)總線上讀出。mmrx,mmry為AR0～AR7或SP；MMR為任何一個(gè)存儲(chǔ)器映象寄存器；例5-6將數(shù)組x[5]初始化為{1,2,3,4,5}。利用雙操作數(shù)指令可以節(jié)省機(jī)器周期。越多，節(jié)省的機(jī)器周期數(shù)也越多。周期，節(jié)省的總機(jī)器周期數(shù)=1T*N=NT?？偩€，得到32位操作數(shù)。這樣，有兩種數(shù)據(jù)排列方法：。低地址存放低16位操作數(shù)。

　　

【正文】 A右移 16位，即存高字節(jié)） ‖LD *AR4,B ； B=1（ 7FFFh左移 16位放在 B的高字節(jié)） MASR *AR2+,*AR3+,B,A ； A=1x2/72,T=x2， AR2指向 d_temp， ； AR3指向 c2 （湊整運(yùn)算為結(jié)果加 215再對(duì) 150位清 0 ） MPYA A ； A=T*A=x2（ 1x2/72） STH A,*AR2 ；（ d_temp） =x2（ 1x2/72） MASR *AR2,*AR3+,B,A ； A=1x2/42（ 1x2/72）； T=x2（ 1x2/72） 63 MPYA *AR2+ ； B=x2（ 1x2/42（ 1x2/72）） ST B,*AR2 ；（ d_temp） =x2（ 1x2/42（ 1x2/72）） ‖LD *AR4,B ； B=1 MASR *AR2,*AR3+,B,A ； A=1x2/20（ 1x2/42（ 1x2/72）） MPYA *AR2+ ； B=x2（ 1x2/20（ 1x2/42（ 1x2/72））） ST B,*AR2 ； （ d_temp） =B=x2（ 1x2/20（ 1x2/42（ 1x2/72））） ‖LD *AR4,B ； B=1 MASR *AR2,*AR3+,B,A； A=1x2/6（ 1x2/20（ 1x2/42（ 1x2/72））） MPYA d_x ； B=x（ 1x2/6（ 1x2/20（ 1x2/42（ 1x2/72）））） STH B,d_sinx ； sin（ theta） RET .end 64 2. 用 z變換計(jì)算余弦值 下面采用 z變換和反 z變換來(lái)求 cos（ x ） ，其的方法為： 根據(jù) z變換定義 ， 序列 x(n) 的變換公式為： 反 z變換為 x(n)=163。1[X(z)]可用留數(shù)法 、 部分分式展開法和長(zhǎng)除法求得 。 則 cos(ω 0n)u(n)的 z變換為： H（ z） =163。[cos(n)u(n)]= 163。[ejω 0n+ejω 0nu(n)] 2165 其中 C＝ cosω 0， A=2cosω 0T， B=1。 ω 0為余弦輸出信號(hào)的頻率 ， T為離散余弦序列的采樣頻率 。 如果以該函數(shù)設(shè)計(jì)一離散時(shí)間系統(tǒng) ， 為則其單位沖擊響應(yīng)就是余弦輸出信號(hào) 。 此時(shí)的輸出序列Y（ k） 為 H（ z） 的反 z變換 。 66 Y（ k） =163。1[H（ z） ]=AY[k1]+BY[k2]+X[k]+CX[k1] 當(dāng) k=1時(shí) Y（ k） =Y（ 1） =AY[2]+BY[3]+X[1]+CX[2]=0 當(dāng) k=0時(shí) Y（ k） =Y（ 0） =AY[1]+BY[2]+X[0]+CX[1]=0+0+1+0=1 當(dāng) k=1時(shí) Y（ k） =Y（ 1） =AY[0]+BY[1]+X[1]+CX[0]=1+0+0+C=A+C 當(dāng) k=2時(shí) Y（ k） =Y（ 2） =AY[1]+BY[0]+X[2]+CX[1]=AY[1]+BY[0] 當(dāng) k=3時(shí) Y（ k） =Y（ 3） =AY[2]+BY[1]+X[3]+CX[2]=AY[2]+BY[1] 當(dāng) k=n時(shí) Y（ k） =Y（ n） =AY[n1]+BY[n2] 在 k＞ 2以后 ， Y（ k） 能用 Y[k1]和 Y[k2]算出 ， 這是一個(gè)遞歸的差分方程 。 如果按第七章實(shí)驗(yàn)二的方式產(chǎn)生余弦信號(hào) ，對(duì)應(yīng)的初始化程序?yàn)椋? 初始化 y[1]和 y[2]： SSBX FRCT ；置 FRCT=1， 準(zhǔn)備進(jìn)行小數(shù)乘法運(yùn)算 ST INIT_A,AA ；將常數(shù) A裝入變量 AA ST INIT_B,BB ；將常數(shù) B裝入變量 BB ST INIT_C,CC ；將常數(shù) C裝入變量 CC 67 PSHD CC ；將變量 CC壓入堆棧 POPD Y2 ；初始化 Y2=CC LD AA, A ；裝 AA到 A累加器 ADD CC, A ； A累加器 =AA+CC STH A,Y2 ； Y2=Y[1]=AA+CC LD AA,T ；裝 AA到 T寄存器 MPY Y2,A ； Y2乘系數(shù) A， 結(jié)果 Y[1]*AA放入 A累加器 ADD BB, A ； A累加器 =Y[2]=Y[1]*AA+BB*Y[0] STH A,Y1 ；將 A累加器中 Y[2]的高 16位存入變量 Y1=Y[2] 68 以后的遞推過(guò)程由中斷服務(wù)程序完成 Y[3]到 Y[n]運(yùn)算 ， 相應(yīng)的程序片段為： LD BB,T ；將系數(shù) B裝入 T寄存器 MPY Y2,A ； Y2乘系數(shù) B， 結(jié)果 BB*Y[1]放入 A累加器 LTD Y1 ； 將 Y1=Y[2]裝入 T， 同時(shí)復(fù)制到 Y2， Y[2]退化為 Y[1] MAC AA,A ；完成新余弦數(shù)據(jù)的計(jì)算 ， A累加器中為 ； Y1*AA+Y2*BB或 Y[3]＝ AA*Y[2]+BB*Y[1] STH A,1,Y1 ；將新數(shù)據(jù)存入 Y1， 因所有系數(shù)都除過(guò) 2， ；所以在保存結(jié)果時(shí)左移一位 ， 恢復(fù)數(shù)據(jù)正常大小 。 STH A,1,Y0 ；將新正弦數(shù)據(jù)存入 Y0 69 第四節(jié) 用 DSP實(shí)現(xiàn) FIR濾波器 1． FIR濾波器基本概念 1. FIR濾波器沒(méi)有反饋回路，因此它是無(wú)條件穩(wěn)定系統(tǒng)，其單位沖激響應(yīng) h（ n）是一個(gè)有限長(zhǎng)序列。 2. FIR濾波算法實(shí)際上是一種乘法累加運(yùn)算。 要 點(diǎn) 70 2. FIR濾波器中 z1的實(shí)現(xiàn) 對(duì)于 N級(jí)的 FIR濾波器 ， 在數(shù)據(jù)存儲(chǔ)器中開辟一個(gè)稱之為滑窗的 N個(gè)單元的緩沖區(qū) ， 存放最新的 N個(gè)輸入樣本；從最老的樣本開始 ， 每讀一個(gè)樣本后 ， 將此樣本向下移位 ， 讀完最后一個(gè)樣本后 ， 輸入最新樣本至緩沖區(qū)的頂部 。 用線性緩沖區(qū)實(shí)現(xiàn) z1的優(yōu)點(diǎn)是 ， 新老數(shù)據(jù)在存儲(chǔ)器中存放的位置直接明了 。 特點(diǎn) (1)用線性緩沖區(qū)法實(shí)現(xiàn) z1 71 C54x片內(nèi)沒(méi)有 I/O資源 ， CPU通過(guò)外部譯碼可以尋址 64K的 I/O單元 。 有兩條實(shí)現(xiàn)輸入和輸出的指令： PORTR PA,Smem ；將為 PA的端口內(nèi)容送 ；數(shù)據(jù)存儲(chǔ)器 Smem PORTW Smem,PA ；將地址為 Smem的數(shù)據(jù) ；存儲(chǔ)器內(nèi)容送端口 PA 數(shù)據(jù)的輸入 /輸出 72 N=6的線性緩沖區(qū)存儲(chǔ)器圖 ????50ii )in(xa)n(y73 使用存儲(chǔ)器延時(shí)指令 DELAY， 可以將數(shù)據(jù)存儲(chǔ)單元中的內(nèi)容向較高地址的下一單元傳送 。 實(shí)現(xiàn) z1的運(yùn)算指令為： DELAY Smem ； （ Smem） → Seme+1， 即數(shù)據(jù)存儲(chǔ) ；器單元的內(nèi)容送下一高地址單元 DELAY *AR2 ； AR2指向源地址 ， 即將 AR2所指單 ；元內(nèi)容復(fù)制到下一高地址單元中 存儲(chǔ)器的延時(shí)操作 延時(shí)指令與其它指令 的 結(jié)合 LT+DELAY→LTD 指令 ；單數(shù)據(jù)存儲(chǔ)器的值裝入 ； T寄存器并送下一單元延時(shí) MAC+DELAY→MACD 指令 ；操作數(shù)與程序存儲(chǔ)器值相乘 ；后累加并送下一單元延時(shí) 74 在數(shù)據(jù)存儲(chǔ)器中開辟一個(gè)稱之為滑窗的 N個(gè)單元的緩沖區(qū) ， 滑窗中存放最新的 N個(gè)輸入樣本；每次輸入新樣本時(shí) ， 以新樣本改寫滑窗中的最老的數(shù)據(jù) ， 而滑窗中的其它數(shù)據(jù)不作移動(dòng)；利用片內(nèi)BK（ 循環(huán)緩沖區(qū)長(zhǎng)度 ） 寄存器對(duì)滑窗進(jìn)行間接尋址 ， 循環(huán)緩沖區(qū)地址首尾相鄰 。 利用循環(huán)緩沖區(qū)實(shí)現(xiàn) Z1的優(yōu)點(diǎn)是不需要移動(dòng)數(shù)據(jù) ，不存在一個(gè)機(jī)器周期中要求能一次讀和一次寫的數(shù)據(jù)存儲(chǔ)器 ， 因而可以將循環(huán)緩沖區(qū)定位在數(shù)據(jù)存儲(chǔ)器的任何位置 （ 線性緩沖區(qū)要求定位在DARAM） 。 特點(diǎn) （ 2）用循環(huán)緩沖區(qū)法實(shí)現(xiàn) z1 75 N=6的循環(huán)緩沖區(qū)存儲(chǔ)器圖 ????50ii )in(xa)n(y76 3． FIR濾波器的實(shí)現(xiàn)方法 例 525 編寫 N=5， y（ n）=a0*x（ n） +a1*x（ n1）+a2*x（ n2） +a3*x（ n3）+a4*x（ n4） 的計(jì)算程序 。 先將系數(shù) a0～ a4存放在數(shù)據(jù)存儲(chǔ)器中 ， 然后設(shè)置線性緩沖區(qū) ， 用以存放輸入和輸出數(shù)據(jù) 。 (1)用線性緩沖區(qū)和直接尋址方法實(shí)現(xiàn) FIR 線性緩沖區(qū)安排 77 .title “” ；定義源程序名 .mmregs ；定義存儲(chǔ)器映象寄存器 .def start ；定義語(yǔ)句標(biāo)號(hào) start .bss y,1 ；為結(jié)果 y預(yù)留 1個(gè)單元的空間 XN .usect “XN”,1 ；在自定義的未初始化段 “ XN” XNM1 .usect “XN”,1 ；中保留 5個(gè)單元的空間 XNM2 .usect “XN”,1 XNM3 .usect “XN”,1 XNM4 .usect “XN”,1 A0 .usect “A0”,1 ；在自定義的未初始化段 “ A0” A1 .usect “A0”,1 ；中保留 5個(gè)單元的空間 A2 .usect “A0”,1 A3 .usect “A0”,1 A4 .usect “A0”,1 PA0 .set 0 ；定義 PA0為輸出端口 PA1 .set 1 ；定義 PA1為輸入端口 78 .data table: .word 1*32768/10 ；假定程序空間有五個(gè)參數(shù) .word 3*32768/10 .word 5*32768/10 .word 3*32768/10 .word 1*32768/10 .text start: .SSBX FRCT ；設(shè)置進(jìn)行小數(shù)相乘 STM A0,AR1 ； 將數(shù)據(jù)空間用于放參數(shù)的首地址送 AR1 RPT 4 ；重復(fù)下條指令 5次傳送 MVPD table,*AR1+ ；傳送程序空間的參數(shù)到數(shù)據(jù)空間 LD XN,DP ；設(shè)置數(shù)據(jù)存儲(chǔ)器頁(yè)指針的起始位置 PORTR PA1,@XN ； 從數(shù)據(jù)輸入端口 I/O輸入最新數(shù)據(jù) x（ n） FIR1: LD @XNM4,T ； x（ n4） → T 79 MPY @A4,A ； a4*x（ n4） → A LTD @XNM3 ； x（ n3） → T,x（ n3） → x（ n4） MAC @A3,A ； A+a3*x（ n3） → A LTD @XNM2 ； x（ n2） → T,x（ n2） → x（ n3） MAC @A2,A ； A+a2*x（ n2） → A LTD @XNM1 ； x（ n1） → T,x（ n1） → x（ n2） MAC @A1,A ； A+a1*x（ n1） → A LTD @XN ； x（ n） → T,x（ n） → x（ n1） MAC @A0,A ； A+a0*x（ n） → A STH A,@y ；保存 y（ n） 的高字節(jié) PORTW @y,PA0 ；輸出 y（ n） BD FIR1 ；執(zhí)行完下條指令后循環(huán) PORTR PA1,@XN ；輸入 x（ n） .en