基于最小二乘法模型的FAST 反射面形狀調(diào)節(jié)問題研究

羅忠旭，秦小華，陳 旭

（新疆農(nóng)業(yè)大學，新疆 烏魯木齊 830000）

1 問題重述

射電天文學對天體物理研究的多個領(lǐng)域做出了巨大貢獻，引起越來越多天文學家的重視，世界各國都在努力突破毫米波射電望遠鏡的技術(shù)困難，以建造更大的高精度的天線和實現(xiàn)對較短波段的觀測[1]?！爸袊煅邸笔且患?00 m 口徑球面射電望遠鏡（Five-hundred-meter Aperture Spherical radio Telescope，縮寫FAST），是目前世界上最大、最靈敏的單口徑射電望遠鏡[2]。主動變形反射面是FAST 主要創(chuàng)新點之一，中國天文界公認整體索網(wǎng)支撐是作為大射電望遠鏡（FAST）主動反射面的最佳支撐結(jié)構(gòu)[3]。其工作原理是通過下拉索與促動器的配合來完成基準球面到工作拋物面的轉(zhuǎn)換，使被觀測體發(fā)出的電磁波經(jīng)反射后聚集在饋源艙處。工作拋物面隨著物體的移動而不斷變化，形成一個動態(tài)的“照明區(qū)域”，實現(xiàn)動態(tài)接收。天線對于射電望遠鏡來說，主要是把微弱的宇宙射電輻射信號收集起來，傳送到接收機中去。收集到的信號越強烈，得到的訊息也就越豐富[4]。因此，調(diào)節(jié)成拋物面是主動反射面技術(shù)的關(guān)鍵。

根據(jù)主動反射面技術(shù)建立數(shù)學模型和算法解決如下問題。

問題：關(guān)于反映觀測天體S 的方位角度α，仰角β 分別等于0°和90°，即被觀測天體位于基準球面正上方時，確定理想拋物面。如圖1 所示。

圖1 天體S 方位角與仰角示意圖

2 問題分析

對于問題，當方位角α、仰角β 分別等于0°和90°時，確定理想的拋物面，建立了以下3 種模型。模型1：理想拋物面的頂點與基準態(tài)球面的最低點重合；模型2：理想拋物面的頂點與基準態(tài)球面的最低點不重合時，但是理想拋物面的口徑端點在基準態(tài)球面上；模型3，理想拋物面的一般情況。

3 模型假設(shè)

（1）基準態(tài)下，所有主索節(jié)點均位于基準球面上。

（2）認為反射面板是無孔的。

（3）電磁波信號及反射信號均視為直線傳播。

（4）三角形反射面板面積相等且是均質(zhì)的。

4 符號說明

表1 符號說明

5 模型的建立與求解

5.1 問題模型的建立與求解

模型1：

由圖2 可知，F(xiàn)=0.466R=139.986 4 m，p=2F=279.6 m，R 為基準態(tài)球面的半徑。

圖2 模型1 示意圖

可以確定一個理想拋物面，其母線方程，即zOx 坐標面的拋物線方程：

將數(shù)據(jù)代入方程（2），可得理想拋物面方程：

模型2：當理想拋物面的頂點與基準態(tài)球面的最低點不重合，但口徑端點落在基準態(tài)球面上時，如圖3 所示。

圖3 模型2 示意圖

設(shè)理想拋物面的母線，即zOx 坐標面的拋物線的方程：

由幾何關(guān)系可求得該拋物線與基準態(tài)球面的一個交點，將這點坐標，和已知條件R、F 代入拋物線方程中，可求解p，從而求出此拋物線方程，將該拋物線繞z 軸旋轉(zhuǎn)得理想拋物面。流程圖如圖4 所示。

圖4 模型2 計算流程圖

如圖3 所示，在△CMN 中，

下面將利用“FAST 基準態(tài)的2 226 個主索節(jié)點空間位置坐標”，基于最小二乘法原理，通過建立和求解優(yōu)化模型，計算得出方程（8）中的參數(shù)d1，從而確定最終的理想拋物面方程。

通過分析照明區(qū)域的口徑長度范圍，在照明范圍內(nèi)的主索節(jié)點滿足x2+y2≤1502條件。通過Excel 設(shè)置條件語句，篩選出滿足的706 個節(jié)點，即照明范圍內(nèi)共有706個主索節(jié)點。但是由于python 軟件的技術(shù)限制，無法同時對706 個節(jié)點進行數(shù)據(jù)分析，通過系統(tǒng)抽樣法，按照5∶1 比例，簡化程序的運算量，篩選得到118 個最有代表性的節(jié)點，進行分析計算參數(shù)。

首先做如下假設(shè)：

設(shè)（xi，yi，zi）來表示初始狀態(tài)第i 個主索節(jié)點坐標，其中i 表示“附件1”的從上往下第i 個節(jié)點。用（xi′，yi′，zi′）來表示（xi，yi，zi）對應(yīng)的FAST 在工作狀態(tài)的各個主索節(jié)點坐標。di為由點（xi，yi，zi）變化到點（xi′，yi′，zi′）時，對應(yīng)的徑向伸縮量，即基準態(tài)主索節(jié)點到最終的工作態(tài)主索節(jié)點之間的距離，且伸縮范圍為-0.6～0.6 m。

現(xiàn)在，對基準態(tài)第i 個主索節(jié)點，其坐標為（xi，yi，zi），將其單位化為

所以（xi′，yi′，zi′）與（xi，yi，zi）的坐標關(guān)系：

再根據(jù)所篩選出來的118 個主索節(jié)點，定義118 個徑向伸縮變量di（i=1，2，...118），基于最小二乘法原理，設(shè)目標方程：

當|di|≤0.6 時，在這個約束條件下，聯(lián)立所有方程組，從而就得到了目標函數(shù)Drt，其自變量為di（經(jīng)篩選存于“照明區(qū)域主索節(jié)點坐標.xlsx”，i 為1 到118）。然后，用python 求這個多元函數(shù)最小值，從而找到最優(yōu)的工作拋物面。

5.2 模型的優(yōu)點

模型的優(yōu)點在于，問題中的模型3 里考慮了706 個變量的影響，把拋物面全部的點都考慮進去了，讓求出的理想拋物面更加精確?；谶@些節(jié)點坐標數(shù)據(jù)，建立了基于最小二乘法的優(yōu)化模型，得到了相對于模型1 和模型2，更加精確的理想拋物面方程。