崔慧敏,郭宇龍,馬紀軍,關(guān)貴注,王宇寧
(北京遙測技術(shù)研究所 北京 100076)
作為一種以激光光束為載體直接進行信息傳輸?shù)男碌耐ㄐ欧绞?,空間激光通信憑借其傳輸速率高、通信容量大、保密性和抗干擾能力強等優(yōu)點,在國內(nèi)和國際上都得到了普遍重視[1,2],在軍用和民用領(lǐng)域也得到廣泛應(yīng)用[3,4]。但是激光的發(fā)射光束窄、發(fā)散角小[5],且通信距離遠,易受復(fù)雜大氣環(huán)境影響。因此,對激光通信終端進行高精度視軸指向控制,對建立和保持高質(zhì)量且穩(wěn)定的激光通信鏈路具有關(guān)鍵性作用。為保證激光通信終端始終跟瞄對方以保證通信鏈路的暢通,目前激光終端多采用粗精組合控制的方式實現(xiàn)視軸的高精度指向[6,7]。其中粗環(huán)的主要作用是實現(xiàn)目標的掃描捕獲和視軸初始指向[8]。受機械諧振,傳感器噪聲以及內(nèi)外部干擾因素的影響,粗跟蹤環(huán)對脫靶量的補償精度有限[9]。為進一步補償粗跟蹤的殘差,通常系統(tǒng)中包含由快速反射鏡構(gòu)成的精跟蹤環(huán)路,實現(xiàn)對目標激光終端信標光的高精度指向。
相比于轉(zhuǎn)動慣量較大的粗跟蹤環(huán)U型架負載,快速反射鏡慣量較小,帶寬較高,可以對輸入進行快速有效的響應(yīng)[10]。為了制定良好的粗精環(huán)控制切換策略,需要提高目標指向精度,最大化粗精環(huán)指向角度等效轉(zhuǎn)換精度。而目前針對解決該問題的相關(guān)文獻較少,且坐標轉(zhuǎn)換方法通常使用四元數(shù)法[11],通過某一點在已知坐標系下的坐標不斷乘以坐標旋轉(zhuǎn)矩陣得到該點在目標坐標系下的坐標。當經(jīng)過多次旋轉(zhuǎn)坐標轉(zhuǎn)換后,會導(dǎo)致最終的方程較為復(fù)雜,目標在已轉(zhuǎn)換坐標系下的位置較難求解。
為解決該問題,本文提出一種新的解決思路,以U型架為基準實現(xiàn)從中間向兩端進行坐標轉(zhuǎn)換。該方法相對于端到端的坐標轉(zhuǎn)換,計算簡單且精度較高。文中首先建立了激光終端各坐標系,并建立了坐標系之間的轉(zhuǎn)換關(guān)系矩陣,接著介紹了使用傳統(tǒng)的端到端坐標轉(zhuǎn)換方法,進行粗跟蹤環(huán)和精跟蹤環(huán)等效角度轉(zhuǎn)換過程,然后提出了一種新的解決思路,并對兩種方法的轉(zhuǎn)換結(jié)果精度和計算復(fù)雜度進行對比分析,最后給出了結(jié)論分析。
為了方便進行角度等效轉(zhuǎn)換,首先需要建立坐標系。跟蹤快速反射鏡坐標系和粗跟蹤環(huán)U型架坐標系的示意圖如圖1所示,其中M1反射鏡,庫德鏡1,庫德鏡2,庫德鏡3和第3反射鏡均處于U型架坐標系中??焖俜瓷溏R可以繞著x軸,y軸做旋轉(zhuǎn)運動。以粗跟蹤環(huán)U型架為基準建立空間坐標系x1y1z1,U型架可以繞著z1軸做方位運動,也可以繞著軸y1做俯仰運動。激光光束通過跟蹤快反鏡,M1反射鏡,3個庫德鏡和第三反射鏡反射后最終指向激光信標的位置。
圖1 U型架和跟蹤快速反射鏡坐標系Fig.1 The U-frame and tracking fast mirror coordinate systems
假設(shè)在快速反射鏡坐標系下有一束光入射到跟蹤快反鏡上,則快反運動引起的光線變化過程如下:
Step1:將跟蹤快反坐標系先繞x軸逆時針轉(zhuǎn)動α1(所對應(yīng)的轉(zhuǎn)換矩陣為R1,表示方法以下類同),再繞z軸逆時針轉(zhuǎn)動β1(R2),到達新的跟蹤快反位置。
Step2:快反坐標系,繞z軸順時針轉(zhuǎn)β1(R3),再繞x軸順時針轉(zhuǎn)α1(R4),接著繞z軸逆時針轉(zhuǎn)動135°(R5),再繞x軸順時針轉(zhuǎn)45°(R6),可得到M1反射鏡坐標系。
Step3:M1反射鏡坐標系,先繞x軸逆時針轉(zhuǎn)動90°(R7),再繞本體坐標系z軸逆時針轉(zhuǎn)動α3(R8_1),可得到庫德鏡3坐標系。
Step4:庫德鏡3坐標系,繞x軸順時針轉(zhuǎn)動180°(R9),可得到庫德鏡2坐標系。
Step5:庫德鏡2坐標系,繞x軸順時針轉(zhuǎn)動90°(R10),可得到庫德鏡1坐標系。
Step6:庫德鏡1坐標系,繞x軸逆時針轉(zhuǎn)動45°(R11),再繞z軸逆時針轉(zhuǎn)動135°(R12),再繞本體坐標系y軸順時針轉(zhuǎn)動β3(R13_1),得到終端本體坐標系。
假設(shè)一束光線在初始跟蹤快反坐標系下的坐標為[x0y0z0],將快反坐標系的輸入光線,經(jīng)過多次坐標變化和反射,可得到最終輸出光線在終端本體坐標系下坐標為[x1y1z1]:
其中R17_1為繞本體坐標系z軸順時針轉(zhuǎn)動α3,R18代表繞z軸逆時針轉(zhuǎn)動45°。
假設(shè)跟蹤快反鏡轉(zhuǎn)動角度α1=0,β1=0。U型架在初始位置基礎(chǔ)上繞z1軸逆時針轉(zhuǎn)動到新的角度α2,繞y1軸順時針轉(zhuǎn)動到新的角度β2。即按以下步驟進行旋轉(zhuǎn)變換:
Step1:在跟蹤快反鏡初始位置下,先繞z軸順時針轉(zhuǎn)β1(R3),再繞x軸順時針轉(zhuǎn)α1(R4),接著繞z軸逆時針轉(zhuǎn)動135°(R5),再繞x軸順時針轉(zhuǎn)45°(R6),到達M1反射鏡。
Step2:在當前坐標系下,先繞x軸逆時針轉(zhuǎn)動90°(R7),再繞主z軸逆時針轉(zhuǎn)動α2(R8_2),到達庫德鏡3。
Step3:在當前坐標系下,繞x軸順時針轉(zhuǎn)動180°(R9),到達庫德鏡2。
Step4:在當前坐標系下,繞x軸順時針轉(zhuǎn)動90°(R10),到達庫德鏡1。
Step5:在當前坐標系下,繞x軸逆時針轉(zhuǎn)動45°(R11),再繞z軸逆時針轉(zhuǎn)動135°(R12),再繞主y軸順時針轉(zhuǎn)動β2(R13_2),到達第三反射鏡。入射光線經(jīng)過多次反射,最終輸出光線在U型架坐標系下坐標為[x2y2z2]:
假設(shè)快速反射鏡在x,y方向的偏轉(zhuǎn)角度分別為Δx,Δy。在快速反射鏡坐標系下,入射光線坐標為[-1 0 0]。跟蹤快速反射鏡初始法線為。跟蹤快反運動后的法線和M1反射鏡法線分別為
則入射光線經(jīng)過快速反射鏡和M1反射鏡反射后,在XOY平面上產(chǎn)生的變化為
U型架上庫德鏡3、庫德鏡2、庫德鏡1和第三反射鏡初始法線分別為
第三反射鏡繞y軸轉(zhuǎn)動β角,旋轉(zhuǎn)矩陣為Ry,四個反射鏡均繞z軸旋轉(zhuǎn)α角,旋轉(zhuǎn)矩陣分別為Rz。旋轉(zhuǎn)后庫德鏡3、庫德鏡2、庫德鏡1和第三反射鏡的新的法線方向如下:
當有一束入射光線[0 0 1],經(jīng)過多個庫德鏡和反射鏡反射后,最終輸出矢量為
當輸入矢量為[ΔxΔy1]時,輸出矢量分別為
將矢量變換到U型架坐標系
假設(shè)RUout3=[x3y3z3]′,第三反射鏡矢量在U型架坐標系為
從整個推導(dǎo)過程可以很容易看出,方法二要比方法一計算量小,使用MATLAB對兩種方法的轉(zhuǎn)換結(jié)果進行計算,由于方程的解很長,這里不再一一列出,但方法二的解很明顯要比方法一的解簡單很多,且很多情況下,在相同的輸入條件下方法二可以很容易得到等效解,但方法一在該條件下卻無解。第一種從端到端的坐標轉(zhuǎn)換方法其本質(zhì)是通過坐標系的旋轉(zhuǎn)得到轉(zhuǎn)臺新角度,快反鏡角度和轉(zhuǎn)臺原來角度之間的關(guān)系(α2,β2)=f(α1,β1,α3,β3),從而進一步得到轉(zhuǎn)臺角度變化量和快反鏡方位、俯仰角度之間的直接轉(zhuǎn)換關(guān)系
從公式1和公式2可以看出,轉(zhuǎn)換過程需要多個矩陣的乘法運算,這對于對實時性要求較高的底層控制器來說,計算量太大。另外,由于激光波束寬度有限,對快反鏡的調(diào)節(jié)角度范圍通常在10個微弧度量級,過多的轉(zhuǎn)換過程會導(dǎo)致最終的計算誤差增大,嚴重情況下有可能造成目標跟蹤丟失。而本文提出的第二種從兩端到中間的坐標轉(zhuǎn)換方法,取U型架坐標系為中間參考坐標系,將不同條件下光線坐標變化同時映射到U型架坐標系中進行對比分析,并采用法線旋轉(zhuǎn),點積運算等手段,使得計算過程在每一階段都有精確的結(jié)果輸出,有效地減小了計算量,滿足工程應(yīng)用中對實時性的要求。
如在U型架方位角和俯仰角分別為30°、60°情況下,如果快反鏡在x方向和y方向的轉(zhuǎn)動量分別為0.01°和0.02°,入射光線坐標為(-1,0,0),經(jīng)過快反鏡和M1鏡反射后,光線在xoy平面產(chǎn)生的變化量為[0.0399,-0.0144,0.9991],經(jīng)過庫德鏡3反射后,該矢量在U型架坐標系的坐標為[-0.0187,-1,-0.0123],經(jīng)過庫德鏡2、庫德鏡1和第3反射鏡反射后,該矢量坐標為[0.0013,0.0223,1]。令快反鏡在x和y方向沒有轉(zhuǎn)動量,將該矢量轉(zhuǎn)換到U型架坐標系,從而得到方位軸和俯仰軸角度變化量為-0.0011547°和0.002°。而方法一由于轉(zhuǎn)換矩陣和逆矩陣較多,MATLAB中最終得到的計算公式很長且無法化簡,這對于處理能力有限的底層嵌入式控制器來說,無法實現(xiàn)工程應(yīng)用。
為實現(xiàn)激光終端高質(zhì)量穩(wěn)定通信,對激光終端進行粗精控制切換時,為解決一系列復(fù)雜的坐標轉(zhuǎn)換問題建立了合適的坐標系并采用兩種方法對角度轉(zhuǎn)換進行求解,并針對該激光終端研究對象,得出采用從兩端向中間參考坐標系進行轉(zhuǎn)換的方法比從端到端的坐標轉(zhuǎn)換方法的計算量小,且求解精度高的結(jié)論。這為激光通信終端的精確、平滑指向提供了理論基礎(chǔ),對促進我國空間激光通信事業(yè)的發(fā)展具有重要的意義。