【正文】 理數(shù)據(jù)的數(shù)量級 不同 ，最后 在整個降落區(qū)域的范圍內(nèi)搜索最優(yōu)著陸點 ；由于在緩速下降和自由落體階段中，發(fā)動機已經(jīng)關閉，故僅對其做簡單物理分析。 本文所研究的問題二主要為過程的最優(yōu)控制和建立嫦娥三號軟著陸軌道。 本文所研究的問題一主要為基礎計算和物理知識 ,首先我們需要根據(jù) 預定的著陸點的經(jīng)緯度確定軌道 ，然后通過拋物線的運動計算出在月球著陸時的水平路程，然后計算出偏移角度，據(jù)此確定近月點的經(jīng)緯度 ， 而嫦娥三號 的著陸軌道為過月球中心點的橢圓軌道，所以 遠月點的經(jīng)緯度和近月點對稱， 則可以由近月點計算出遠月點 的經(jīng)緯度 。 （ 2） 確定嫦娥三號的著陸軌道和在 6個階段的最 優(yōu)控制策略。其著陸軌道設計的基本要求：著陸準備軌道為近月點 15km，遠月點 100km 的橢圓形軌道；著陸軌道為從近月點至著陸點，其軟著陸過程共分為 6 個階段，要求滿足每個階段在關鍵點所處的狀態(tài)；盡量 減少軟著陸過程的燃料消耗 。嫦娥三 號的 預定著陸點 為 ， ， 海拔 為 2641m。嫦娥三號在著陸準備軌道上的運行質(zhì)量為 ，其安裝在下部的主減速發(fā)動機能夠產(chǎn)生 1500N到 7500N的可調(diào)節(jié)推力，其比沖（即 單位質(zhì)量的推進劑產(chǎn)生的推力 ）為 2940m/s， 可以 滿足調(diào)整 速度 的 控制 要求。 然后進行敏感性分析 ，主要利用蒙特卡洛分析著陸軌道的粗避障階段和精避障階段月面不同地形高度，對嫦娥三號降落時所需調(diào)整概率大小的影響，接著分析嫦娥三號著陸占地面積大小對著陸調(diào)整概率 P的影響。首先以模型二為基礎進行誤差分析，當 主減速階段的推力 ??????????????、初始質(zhì)量 ??2變化時，計算嫦娥三號質(zhì)量 ??2和燃料消耗速率的變化趨勢。 其余 階段 僅對其做簡單物理分析后繪制出六個階段的著陸軌道。 模型二主要 對主減速階段和快速調(diào)整階段進行初步分析。最后在軟著陸的橢圓軌道上，由動力勢能和重力勢能的變化，計算出嫦娥三號在遠月點和近月點的速度分別為 1700??/??和 1615??/??， 沿軌道切線方向 。 E，176。 W， 176。 。首先以月球中心為坐標原點建立空間坐標系，根據(jù)計算 的作用力 可知地球影響較小，故忽略不計。 20xx 高教社杯全國大學生數(shù)學建模競賽 編 號 專 用 頁 賽區(qū)評閱編號（由賽區(qū)組委會評閱前進行編號）： 賽區(qū)評閱記錄（可供賽區(qū)評閱時使用）： 評 閱 人 評 分 備 注 全國統(tǒng)一編號（由賽區(qū)組委會送交全國前編號）： 全國評閱編號（由全國組委會評閱前進行編號）： 嫦娥三號軟著路軌道設計與控制策略 摘要 本文主要為分階段研究嫦娥三號的軟著陸軌道設計與最優(yōu)控制策略。 建立模型一 確定 近月點 和遠月點的 位置 ，以及嫦娥三號速度 大小與方向 。然后將嫦娥三號軟著陸看作拋物線的運動過程，計算在最大推力下的減速運動，求得月面偏移距離為 km，由此計算出偏移角度為 176。從而得出近月點和遠月點的經(jīng)緯度分別為 （ 176。 N） 和 （ 176。 S） 。 建立模型二和模型 三 確定著陸軌道和在 6 個階段的最 優(yōu)控制策略 。 模型三分六個階段確定軌道和最優(yōu)控制策略，主減速階段建立目標函數(shù) 燃料 ???， 假設推力 ??????????????最大 ,將 最優(yōu)燃耗軟著陸問題轉(zhuǎn)化為最短時間控制問題，然后 采用擬牛頓法和四階 Admas 預測 校正 得到 ??? = ； 快速調(diào)整階段采用重力轉(zhuǎn)彎制導， 在假設條件下對嫦娥三號進行受力分析，得到嫦娥三號的動力學模型，然后通過開關控制得到燃耗最優(yōu)控制，并畫出仿真圖；粗避障階段采用多項式制導， 通過初始狀態(tài)和末端狀態(tài)反解多項式系數(shù)進而求取標稱軌跡 ， 然后 將避障區(qū)域網(wǎng)格化，比較網(wǎng)格內(nèi)的方差大小確定最優(yōu)區(qū)域范圍 ； 精避障階段 需 在滿足本文提出的避障原則式 (18)下 搜索全局最優(yōu)解，以網(wǎng)格區(qū)域總體得分 ζ作為目標函數(shù) ， 得到最優(yōu)區(qū)域為坐標 (544,278)附近 ，并以螺旋搜索法搜索安全半徑的個數(shù) 。 建立模型四 做相應的誤差分析和敏感性分析 。再以模型三為基礎進行分析， 對初始高度變化前后主減速階段的 ??????????????的 偏角 ??和 ??1和著陸軌道進行對比分析并計算誤差。 關鍵字： 拋物線、燃料 ?、擬牛頓法、 Admas、網(wǎng)格化、蒙特卡洛模擬 1 嫦娥三號于 20xx年 12月 2日 1時 30分成功發(fā)射， 12月 6日 抵達 月球軌道 。在四周安裝有姿態(tài)調(diào)整發(fā)動機，在給定主減速發(fā)動機的推力方向后，能夠自動通過多 個發(fā)動機的脈沖組合實現(xiàn) 各種 姿態(tài) 的調(diào)整控制。 嫦娥三號在高速飛行的情況下，要保證準確地在月球預定區(qū)域內(nèi)實現(xiàn)軟著陸，關鍵問題是著陸軌道與控制策略的設計。 根據(jù)上述的基本要求，請你們建立數(shù)學模型解決下面的問題： （ 1） 確定著陸 準備軌道近月點 和遠月點的 位置 ，以及嫦娥三號相應速度的大小與方向。 （ 3）對于你們設計的著陸軌道和控制策略做相應的誤差分析和敏感性分析。最后因為在著陸軌道上 衛(wèi)星的能量守恒， 則可以通過勢能和動能的轉(zhuǎn)換來計算嫦娥三號的速度和方向 。 因為嫦娥三號的軟著陸主要分為六個階段，所以此問應分為六個階段來求解。最后通過整個分析得出總的著陸軌道。 然后進行敏感度分析，即 對 著陸軌道的粗避障階段和精避障階段月面不同地形高度對嫦娥三號 降落時所需調(diào)整概率大小的影響， 最后分析嫦娥三號著陸占地面積大小對著陸調(diào)整概率 P的影響。11為引力常量， m、 M 分別為兩物體質(zhì)量， R為兩物體距離 ， 為兩物體間的作用力 = m ： ??為物體質(zhì)量， 為物體在 ??作用下產(chǎn)生的加速度 ????： 軟著陸起始速度 ??： 加速度 X： 平拋產(chǎn)生的距離 ??????： 物體的動能（ i = 1,2?） ??????： 物體的重力勢能（ i = 1,2?） ??????????????： 嫦娥三號的推力 λ ： 偏好系數(shù) ζ： 降落地點總體得分 ??： 第 ??段離散段的平均加速度 由于本文使用參數(shù)和公式較多，其他公式和符號在 具體 模型中再做說明。 表 1 地球及月球在近月點和遠月點的重力加速度 （單位 ： / ） 近月點 遠月點 地球 月球 由上表可知，地球在近月點和遠月點的 重力加速度 數(shù)值很小 ， 即地球?qū)︽隙? 三號與月球影響很小， 故可忽略不計。 模型的分析 根據(jù)附件 2給出的軟著陸過程 示意圖 ，即嫦娥三號將在近月點 15 公里處以拋物線下降，相對速度從每秒 公里逐漸降為零。 利用 matlab 繪制嫦娥三號繞月飛行的三維動態(tài)圖，更直觀的反應嫦娥三號的環(huán)月飛行 ，如 圖 3（ 源程序見附錄）： 圖 2 嫦娥三號繞月軌道坐標 圖 圖 3 嫦娥三號環(huán)月飛行 同時由附件二所給的 嫦娥三號著陸區(qū)域和著陸點示意圖 可知，只要保證嫦娥三號的著陸區(qū)域在虹灣著陸區(qū)，則認為著陸成功 。 作為軟著陸 的繞月軌道 。一個是繞月 的向心力，一個是 與繞行面相切的力，則 選擇 最終狀態(tài)為繞赤道運行 更為準確 。 圖 4 嫦娥三號 繞月飛行軌道分析 模型的 建立與 計算 據(jù)了解，嫦娥三號主發(fā)動機是目前中國航天器上最大推力的發(fā)動機，能夠產(chǎn)生從 1500 牛到 7500 牛的可調(diào)節(jié)推力，故可根據(jù)推力范圍求取嫦娥三號的加速度范圍。 主減速段 看作 平拋運動 ： 起始速度 ???? = 103 ?? ??? 加速度的取值范圍 ≤ ?? ≤ ?? ??2? 4 平拋產(chǎn)生的距離 x = ????22????= ????（ ?? = ?? ??2? ） 圖 5 嫦娥三號拋物示意圖 由上圖，并結(jié)合計算所得的拋物距離，得到準備著陸的點與軟著陸點相差176。 由附件所給條件可知距離月球 表面 15km時，速度的大小為 1700??/??，則此速度看作近月點速度， 在穩(wěn)定的軌道下， 從近月 點到遠月點可看作 重力勢能和動能 相互轉(zhuǎn)換的過程 ， 而遠月點距離地球表面為 100km，可以計算重力勢能的變化，即可算出遠月點的速度： ????2 = ????1 ?(?