【正文】 輸入或者設(shè)置一個非標(biāo)準(zhǔn)初始狀態(tài)，然后進(jìn)行一系列運行。 運用上述方法得到的敏感系數(shù)矩陣給出如下： （只給出 S1， S S3略） 3 4 3 1 2 24 5 5 2 3 33 5 4 2 2 21 45 .5 0 2 1 0 2 .0 8 0 1 0 1 .0 5 0 1 0 1 .4 1 8 1 0 4 .0 2 4 1 0 8 .9 3 9 103 .8 5 0 1 0 1 .4 0 1 1 0 7 .3 0 1 1 0 1 .0 0 1 1 0 3 .2 1 0 1 0 4 .0 3 0 101 .6 9 2 1 0 6 .4 1 1 1 0 3 .2 4 0 1 0 4 .4 0 7 1 0 1 .2 3 9 1 0 1 .8 3 3 108 .3 6 2 1 0S? ? ? ? ? ?? ? ? ? ? ?? ? ? ? ? ??? ? ? ? ? ? ? ? ? ?? ? ? ? ? ? ? ?? ? ? ? ? ? ??? 4 3 1 2 14 5 5 2 2 33 5 4 2 2 22 .5 7 0 1 0 1 .8 6 2 1 0 5 .5 8 0 1 0 8 .7 4 2 1 0 1 .4 1 4 1 05 .8 6 0 1 0 1 .4 1 0 1 0 1 .3 1 2 1 0 3 .9 3 6 1 0 1 .1 9 6 1 0 9 .9 0 1 102 .5 7 5 1 0 7 .9 0 2 1 0 5 .7 1 0 1 0 1 .7 3 2 1 0 2 .6 9 0 1 0 4 .5 7 7 10? ? ? ? ?? ? ? ? ? ?? ? ? ? ? ?????? ? ? ? ? ? ? ?? ? ? ? ? ? ? ? ?? ? ? ? ? ? ? ? ??? A A3: 1 ? ?? 、 2, 3 0 .1 2 9 7 2 .1 3 2 9j? ?? A2： 1 ? ? 、 2,3 0 .6 7 6 1 1 .8 9 7 8j? ? ? ? 18 由于數(shù)值仿真的起始點選為（ 1,0， 1），靠近平衡點（ ,0， ），仿真實驗中混沌系統(tǒng)的基頻 0? =，基周期為 002 / 2. 94 44 3TS????。在十進(jìn)制蟻群算法中增加局部搜索策略提高了搜索精度，且能較好的找出推力方向角與時間的關(guān)系。 七、參考文獻(xiàn) [1] 董文永 劉進(jìn) 丁健利 朱福喜，《最優(yōu)化技術(shù)與數(shù)學(xué)建?！罚迦A大學(xué)出版社。 [X,Y]=meshgrid(1:y,1:x)。)。在本題求嫦娥三號軟著陸軌道設(shè)計和控制策 略上，應(yīng)學(xué)會與地球衛(wèi)星情形類比，找出相同和不同之處，合理假設(shè)，進(jìn)行建模。 (2) 采用了依靠月表圖像特征點匹配跟蹤方法后 , 導(dǎo)航系統(tǒng)的水平位置誤差得到修正 ,最終相對目標(biāo)著陸點的導(dǎo)航誤差不大于 2 m。由參考文獻(xiàn)可知，其相互關(guān)系可表示為： 1 2 3f i bc bsp S p S q S q? ? ? （ 27） 其中， S S2和 S3分別表示相對于 ,i bc bsp q q 的誤差敏感系數(shù)矩陣。第 j個觀測量的總估計誤差 jq 由以下四部分組成 ( t )( t ) ( t ) ( t ) ( t )1 0 0 1 0 0jb s jn sjjj jb c jn cqqq q Q q Q? ? ? ? （ 26） 針對主制動這一特定操作階段，上述四部分誤差具有如下特性： jbcq — 第 j 個觀測量的測量誤差，恒為常值，其分布服從零均值高斯分布； jbsq — 第 j 個觀測量的刻度因素誤差系數(shù)，恒為常值，其分布服從零均值高斯分布； jncq — 第 j 個觀測量的隨機(jī)誤差，其為一高斯白噪聲； jnsq — 第 j 個觀測量的刻度因素隨機(jī)誤差系數(shù)，其為一高斯白噪聲。 , 即當(dāng)前 星下點附近特征點可較少 , 但目標(biāo)著陸點附近特征點應(yīng)較多 , 以適應(yīng)著陸器高度不斷下降 , 導(dǎo)航相機(jī)所覆蓋的月面越來越小的變化趨勢。我們需要對附錄 3 和附錄4 的高程圖在 matlab 軟件中精確分析每個區(qū)域塊是否平坦，同樣需要調(diào)整姿態(tài)發(fā)動機(jī)，達(dá)到著陸器水平運動的目的。 用函數(shù)逼近法進(jìn)行參數(shù)化的相關(guān)參數(shù)設(shè)置為：設(shè)將軌跡離散化劃為 6 段，那么待優(yōu)化參數(shù)共 8 個參數(shù)，即 7 個推力方向角和 1 個終端時刻 ft 。 用于多元連續(xù)函數(shù)的 DACA 算法與一元函數(shù)類似。采用函數(shù)逼近法進(jìn)行參數(shù)化。 （ 11） 由圖 521 在 3000 米處經(jīng)緯度坐標(biāo)為 ,由運動方程和逆推法，得近月點的經(jīng)緯度坐標(biāo)為 ， ,距月球表面上 15km。,( ??tf 由方程（ 2）中的攝動加速度 ?F 形成。 五、 模型 的建立與求解 準(zhǔn)備工作 由于涉及到專業(yè)術(shù)語，為此我們首先給出必要的說明： ：又稱軌道參數(shù)，是用來描述天體在其軌道運行狀態(tài)的一組參數(shù)，通常情況下指的是用萬有引力定律描述天體按圓錐曲線運動時 所必需的 6 個參數(shù)。 二、問題 假設(shè) 1．假設(shè)月球為均勻引力場； 2．假設(shè)著陸器軟著陸過程中的月球自轉(zhuǎn)忽略； 3．假設(shè)著陸器在一個固定的鉛垂面內(nèi)運動； 4．在模型一中，假設(shè)僅考慮地球引力攝動和月球非球形引力對著陸器的影響； 5．在模型三中，假設(shè)不考慮攝動項對著陸器的影響； 二、 主要符號說明一覽表 4 符號 符 符號的含義及單位 M 月球的質(zhì)量（ kg） m 著陸器的質(zhì)量 (kg) R 月球平均半徑 (km) eR 月球赤道平均半徑 (km) F 主發(fā)動機(jī)的推力 （ N） r 著陸器與月心距離（ km） v 著陸器的徑向速度 (km/s) ? 著陸器極角（ rad） ? 月球引力常數(shù)（ 22 /kgmN? ） ? 主發(fā)動機(jī)推力方向角，其定義為 F與當(dāng)?shù)厮椒较驃A角（ rad） ISP 主發(fā)動機(jī)比沖（ m/s） AR 橢圓軌道近月點半徑（ km） BR 橢圓軌道遠(yuǎn)月點半徑（ km） J 優(yōu)化性能指標(biāo) 四 、問題 分析 5 題目要求根據(jù)探測器繞月橢圓軌道的半徑及著陸過程的各個階段的狀態(tài)，求出遠(yuǎn)月點和近月點的位置，我們知道在從距離月球表面