基于定點(diǎn)DSP的軟件鎖相環(huán)的設計和實(shí)現

1 軟件鎖相環(huán)

1.1 軟件鎖相環(huán)的結構

圖1表示軟件數字接收機中的解調器。它包括一個(gè)由改進(jìn)的costas環(huán)路^[1]構成的載波跟蹤環(huán)路。

采樣后的中頻信號經(jīng)過(guò)數字混頻,濾掉高頻分量,通過(guò)改進(jìn)的costas環(huán)路產(chǎn)生控制信號,控制數控振蕩器(NCO)得到新的本振參考信號。其中的相位檢測器和環(huán)路濾波器結構如圖2所示。

1.2 軟件鎖相環(huán)的設計

由圖2可見(jiàn),SPLL的計算由計算相位誤差和更新環(huán)路中間變量、輸出控制信號兩部分組成。算法描述可用偽碼表示:

算法描述中的中斷周期就是環(huán)路采樣時(shí)間間隔。

中斷發(fā)生后,第一步讀取基帶同相項數據和正交項數據

A是基帶信號幅度,θ_e是相位誤差。第二步計算硬判決的同相數據乘以相位誤差。

由式(3)、(5)和(6)得:

即把同相數據硬判決后結果乘以正交項數據后再乘以歸一化因子K_norm。K_norm初始值由I_baseband和Q_baseband的初始值決定,

由于定時(shí)恢復環(huán)路和AGC(自動(dòng)增益控制)環(huán)路的作用,K_norm在解調過(guò)程中近似保持恒定。第三步更新環(huán)路中間變量s_pll(n),

f₀是固定的NCO中心頻率。至此,一次完整的SPLL計算完成。

在第三步和第四步計算中,環(huán)路濾波器系數C₁,C₂可以通過(guò)環(huán)路采樣時(shí)間間隔T(或者環(huán)路更新時(shí)間間隔)、環(huán)路自由頻率ω_n及環(huán)路阻尼系數ζ確定,如下兩式^[2]:

K_d為相位檢測器的增益,由于在實(shí)現軟件鎖相環(huán)時(shí),基帶信號的同相和正交分量都經(jīng)過(guò)歸一化處理,故K_d=1;K₀為數控振蕩器的增益,K₀=2πT。T為調制數據速率的倒數(1/76800),ζ一般取0.707。在啟動(dòng)載波恢復之前有一個(gè)頻率捕獲過(guò)程,通過(guò)1024點(diǎn)的FFT,可保證接收的中頻信號與本振信號之間的頻率差Δf₀最大不超過(guò)75Hz。另外系統設計要求SPLL的捕獲時(shí)間(pull in time)小于50ms,由此可以確定環(huán)路自由頻率ω_n。ω_n>150(radHz),取ω_n=300(radHz)。將K₀、K_d、ζ、ω_n、T代入式(3)、(4)得,C₁=67.3273309,C₂=0.1859953。

2 48位定點(diǎn)擴展精度算法

圖3描述的算法在TMS320VC5510上用C語(yǔ)言直接利用浮點(diǎn)運算實(shí)現時(shí),只能通過(guò)C編譯器產(chǎn)生模擬浮點(diǎn)運算的定點(diǎn)指令。這種方法效率很低,每次環(huán)路計算需耗費875個(gè)指令周期。在調制數據速率為76.8kbps的數字接收機中,需要67.2MIPS的運算量。為了降低環(huán)路計算的運算量,同時(shí)保持浮點(diǎn)運算具有動(dòng)態(tài)范圍大、精度高的優(yōu)點(diǎn),筆者提出了一種48位定點(diǎn)擴展精度計算的方法。參加運算的每個(gè)操作數由三個(gè)16位定點(diǎn)數W₂、W₁、W₀級聯(lián)表示,其中高16位為二進(jìn)制補碼的整數部分,低32位為二進(jìn)制補碼的小數部分,符號位在最高位,也可稱(chēng)為Q15.32格式,如圖3所示。

一個(gè)Q15.32數的表示范圍是(-32768,32768),小數分辨率是遠遠超過(guò)16位定點(diǎn)表示的精度,以下用加(ADD_48)、減(SUB_48)、乘(MULT_48)三種基本運算來(lái)說(shuō)明定點(diǎn)擴展精度算法。操作數X由X₂、X₁、X₀構成,操作數Y由Y₂、Y₁、Y₀構成,結果W由W₂、W₁、W₀構成。執行48位加法運算時(shí),W₂W₁W₀=X₂X₁X₀+Y₂Y₁Y₀,首先把小數部分X₁X₀和Y₁Y₀相加,結果保存到W₁W₀中,產(chǎn)生的進(jìn)位位CARRY與X₂、Y₂相加,結果保存到W₂。執行48位減法運算時(shí),W₂W₁W₀=X₂X₁X₀-Y₂Y₁Y₀,首先X₁X₀減去Y₁Y₀,結果保存到W₁W₀,產(chǎn)生借位位BORROW,再由X₂減去Y₂和借位位BORROW,結果保存到W₂。兩個(gè)Q15.32數相乘時(shí),乘積是一個(gè)Q30.64數,出于前面實(shí)現SPLL時(shí)對動(dòng)態(tài)范圍和計算精度的要求,該Q30.64數可以雙向截位為Q15.32的48位定點(diǎn)數。具體做法是保留符號位和整數部分的低15位以及小數部分的高32位。48位定點(diǎn)數的乘法由圖4所示。

除了以上加、減、乘三種基本運算外,48位定點(diǎn)擴展精度算法還包括取負(NEG_48)、數據拷貝(MOVE_48)兩種操作。取負操作即將X₁X₀取負,結果保存到W₁W₀,產(chǎn)生借位位BORROW,再用0減去X₂和借位位BORROW,結果保存到W₂;數據拷貝,即把X₁X₀拷貝到W₁W₀,X₂拷貝到W₂。

在TMS320VC5510可編程DSP的基礎上,利用48位定點(diǎn)擴展精度算法實(shí)現SPLL。在實(shí)現過(guò)程中,采取了模塊化的思路。首先,把SPLL整個(gè)環(huán)路計算封裝成一個(gè)可調用的C語(yǔ)言函數。函數參數包括C₁、C₂、K_norm、基帶信號的I及Q分量、環(huán)路中間變量、調整頻率。DSP的中斷例程(ISR)可以直接調用環(huán)路計算函數,而且通過(guò)輸入不同的C₁、C₂,適用于不同的載波恢復環(huán)路中。另外,在函數內部用匯編語(yǔ)言進(jìn)行編程,以充分利用DSP的計算能力,把48位定點(diǎn)擴展精度算法的五個(gè)基本操作封裝成用匯編指令寫(xiě)的宏(macro),對照計算流程,調用這些宏,完成SPLL的核心計算部分。經(jīng)統計,每次環(huán)路計算需132個(gè)指令周期,總的運算量10.1MIPS,是浮點(diǎn)算法運算量(67.2MIPS)的14%。

用48位擴展精度算法實(shí)現軟件接收機中的SPLL,解決了浮點(diǎn)算法運算量大的問(wèn)題,同時(shí)還具備浮點(diǎn)算法動(dòng)態(tài)范圍大、精度高的優(yōu)點(diǎn),已經(jīng)成功應用于“創(chuàng )新一號”小衛星地面手持低功耗通信終端中。另外,本文提出的SPLL實(shí)現算法,通過(guò)修改環(huán)路濾波器系數,也可以應用在其他軟件接收機中,具有很好的擴展性。