第一篇：全國(guó)大學(xué)生數(shù)學(xué)建模經(jīng)驗(yàn)總結(jié)

全國(guó)大學(xué)生數(shù)學(xué)建模經(jīng)驗(yàn)總結(jié)

我堅(jiān)信:只有想不到的，沒有做不到的！

我是一名渝州學(xué)院大二學(xué)生，2010年加入學(xué)院建模隊(duì)，2011年9月參加了全國(guó)大學(xué)生數(shù)學(xué)建模競(jìng)賽，獲得江西省二等獎(jiǎng)，雖然接觸建模時(shí)間不是很長(zhǎng)，但建模給我?guī)淼膮s很多。

建模對(duì)很多人來說是很模糊的東西，但作為一名建模人就應(yīng)該擔(dān)任起對(duì)建模的責(zé)任。作為一名?？粕?，我知道自己相對(duì)于別人起跑線低，所以進(jìn)入建模隊(duì)后以嚴(yán)格的紀(jì)律來要求自己，別人懂的我得懂，別人不懂的我也要弄懂，不為別的，只因?yàn)橐鼍鸵鲎詈茫?/p>

我雖在有些方面較強(qiáng)于別人，但我上課時(shí)還是認(rèn)真聽取老師的講解，世上沒有相同的人，每個(gè)人的想法思路也不可能完全相同，何不把別人的借鑒過來為己所用，多一個(gè)思路就證明多一條出路，多一條出路解題時(shí)分析題目的能力自然會(huì)比別人想的全面！所以在有些時(shí)候我還是把自己的思路講給同學(xué)們參考，一起討論解題的最好辦法！

記得高三班主任在送我上大學(xué)的時(shí)候說過“不論什么時(shí)候多從別人的角度出發(fā)，凡事不要只為自己，你有足夠的能力去做好任何事，多從事物的本身出發(fā)去考慮，要不做就不要做，但做就要做最好，以后的一切你自己把握！”擔(dān)任過班長(zhǎng)、學(xué)習(xí)委員、各科課代表的我曾對(duì)大學(xué)的錄取不屑一顧，是人才總有發(fā)揮的地方，不論我的選擇怎樣家人、老師總會(huì)支持我。人生是用來闖的，不做錯(cuò)事是不可能的，但做錯(cuò)事后我會(huì)勇敢承認(rèn)，人無圣賢，孰能無過！

加入建模隊(duì)后我認(rèn)識(shí)到，不是任何錯(cuò)誤都可以犯的，也許某步棋沒走對(duì)，那滿盤就盡輸了！所以不論每次訓(xùn)練我都警告自己失誤不是每個(gè)人都犯得起的，建模建的是我們的思維，我們分析問題的能力，處理問題的方法，文字的表達(dá)能力。沒有人會(huì)去當(dāng)面問你解題的思路，我們唯一能做的就是把自己的思路想法用精煉的文字表達(dá)出來，這也鍛煉了我對(duì)人處事的方式，以前總會(huì)想差不多就可以了，不去追求理解事物的本質(zhì)，但建模不允許我們這么做，任何問題都要求我們刨根問底，對(duì)題目不理解何來的思路可循？

單純從快遞的運(yùn)輸方面舉例來說，如果讓我們?cè)O(shè)計(jì)一種最快、最低廉的運(yùn)送方案，大多數(shù)人會(huì)考慮到運(yùn)輸設(shè)備的調(diào)用、運(yùn)輸路線的選擇、接受地點(diǎn)的設(shè)置等方面，但對(duì)于從建模走過的人來說，我們會(huì)考慮的更深一層，例如該選用什么運(yùn)輸方式可以讓運(yùn)費(fèi)最低，可以結(jié)合當(dāng)?shù)亟煌ㄇ闆r，利用單文件多種運(yùn)輸方式結(jié)合進(jìn)行遞送，通過對(duì)各種運(yùn)輸方式了解后，在保證運(yùn)送時(shí)間相當(dāng)?shù)那闆r下，選用最低費(fèi)用的運(yùn)輸方式，當(dāng)然得結(jié)合實(shí)際情況。

總之，建模讓我從另一個(gè)方面看世界，能讓我更深入的分析理解問題，以致讓我用最好的方法進(jìn)行處理！

要做就做最好，不論什么時(shí)候只要確立了目標(biāo)就要堅(jiān)持到底！

第二篇：全國(guó)大學(xué)生數(shù)學(xué)建模大賽經(jīng)驗(yàn)總結(jié)

比賽須知：

1、比賽期間參賽隊(duì)員可以使用各種圖書資料、計(jì)算機(jī)和軟件，在國(guó)際互聯(lián)網(wǎng)上瀏覽。

2、但是不得與隊(duì)外任何人員（包括在網(wǎng)上）討論，否則按舞弊處理。比賽期間不能大聲喧嘩，不能影響別的隊(duì)。

6、數(shù)學(xué)建模大賽時(shí)間安排：

發(fā)題時(shí)間：2014年9月12日至9月15日上午8點(diǎn)

參加數(shù)學(xué)建模競(jìng)賽的十大秘訣誠(chéng)信是最重要的數(shù)學(xué)建模競(jìng)賽是考查學(xué)生研究能力和實(shí)踐能力的一場(chǎng)綜合性比賽，有很多方面的知識(shí)和能力可以考查，但其中我覺得最重要的是誠(chéng)信。我感到中國(guó)在這方面的教育還遠(yuǎn)遠(yuǎn)不夠，我知道有很多同學(xué)寫論文并不是實(shí)事求是地去做，而是編造數(shù)據(jù)、修改結(jié)論，明明自己沒法編程實(shí)現(xiàn)卻硬說自己做出來了，還編了一些數(shù)據(jù)。這些行為也許能夠騙過評(píng)委，也許可以因“此”而獲獎(jiǎng)，但是這對(duì)他們將來是很不利的,希望能夠引起足夠的注意。團(tuán)隊(duì)合作是能否獲獎(jiǎng)的關(guān)鍵

在三天的比賽中，團(tuán)隊(duì)交流所占用的時(shí)間可能會(huì)超過一半。在一個(gè)小組中，出現(xiàn)意見不一是非常正常的，如果一個(gè)隊(duì)意見完全一致，我想他們肯定不會(huì)拿獎(jiǎng)。出現(xiàn)分歧的時(shí)候應(yīng)當(dāng)如何解決是很關(guān)鍵的，甚至直接決定你是否可以獲獎(jiǎng)，我的建議是“妥協(xié)”，這似乎是個(gè)貶義詞，但我的意思是說不要總認(rèn)為自己的觀點(diǎn)是正確的，多聽聽別人的觀

點(diǎn)，在兩者之間謀求共同點(diǎn)。如果三個(gè)人都是自傲類型的人，也許每個(gè)人都非常強(qiáng)，但一旦合作，分歧就無法解決，做出來的就是一團(tuán)糟，也就是說“三個(gè)諸葛亮頂不上一個(gè)臭皮匠”。我奉勸這樣的話最好別組成一隊(duì)了。合作在競(jìng)賽前就應(yīng)當(dāng)培養(yǎng)，比如一塊兒做模擬題什么的，充分利用每個(gè)人的優(yōu)點(diǎn)，也可以張三準(zhǔn)備圖論，李四準(zhǔn)備最優(yōu)化方法，然后幾天后大家一塊交流，這些都是可以磨合團(tuán)隊(duì)之間的關(guān)系的。通常在比賽時(shí)，三個(gè)人的分工是明確的，一個(gè)是領(lǐng)軍人物，主要是構(gòu)建整個(gè)問題的框架并提出有創(chuàng)意的idea，自然其他部分比如論文寫比如程序設(shè)計(jì)比如計(jì)算他也能參加，應(yīng)該算是一名全能型的人物;第二個(gè)是算手，顧名思義，主司計(jì)算方面的問題，比如編程計(jì)算一個(gè)微積分或者手工計(jì)算一條最優(yōu)路徑等。優(yōu)秀的團(tuán)隊(duì)算手一般會(huì)精通(是精通不是入門)一個(gè)軟件的應(yīng)用，比如C比如MATLAB比如LINGO;最后一個(gè)是寫手，主要工作在于論文的寫作和潤(rùn)色上。好的論文要讓人一眼就明了其中的意思，所以寫手的工作還是需要一定的技巧的。當(dāng)然，最重要的還是三個(gè)隊(duì)員之間的討論和交流，同心協(xié)力，在整個(gè)比賽過程中形成一種良好的交流氛圍。時(shí)間和體力的問題

競(jìng)賽中時(shí)間分配也很重要，分配不好可能完不成論文，所以開始時(shí)要大致做一下安排。不必分的太細(xì)，比如第一天做第一小題，第二天做第二小題，這樣反而會(huì)有壓力，一切順其自然。開始階段不忙寫作，可以將一些小組討論的要點(diǎn)記錄下來，不要太工整，隨便寫一下，到第三天再開始寫論文也不遲的。也不要到第三天晚上才開始。另外要

說的就是體力要跟上，三天一般睡眠只有不到10 個(gè)小時(shí)，所以沒有體力是不行的，建議是賽前熬夜編程幾次，既訓(xùn)練了自己的建模能力，也達(dá)到了訓(xùn)練體力的目的，賽前鍛煉身體我覺得沒什么用處，多熬夜就行了，但比賽前一天可不許熬。重視摘要

摘要是論文的門面，摘要寫的不好評(píng)委后面就不會(huì)去看了，自然只能給個(gè)成功參賽獎(jiǎng)。摘要首先不要寫廢話，也不要照抄題目的一些話，直奔主題，要寫明自己怎樣分析問題，用什么方法解決問題，最重要的是結(jié)論是什么要說清楚，在中國(guó)的競(jìng)賽中結(jié)論如果正確一般得獎(jiǎng)是必然的，如果不正確的話評(píng)委可能會(huì)繼續(xù)往下看，也可能會(huì)扔在一邊，但不寫結(jié)論的話就一定不會(huì)得獎(jiǎng)了，這一點(diǎn)不比美國(guó)競(jìng)賽，所以要認(rèn)真寫。摘要至少需要琢磨兩個(gè)小時(shí)，不要輕視了它的重要性。很有必要多看看優(yōu)秀論文的摘要是如何寫的，并要作為賽前準(zhǔn)備的內(nèi)容之一。論文寫作要正規(guī)

論文一定要大致按照摘要、問題重述、模型假設(shè)、符號(hào)說明、問題分析、(建立、分析、求解模型)、模型檢驗(yàn)、參考文獻(xiàn)、附錄等等的方式來寫。一篇論文結(jié)構(gòu)上如果失敗的話，比賽也一定不會(huì)成功，一般初評(píng)會(huì)先淘汰一些結(jié)構(gòu)失敗的文章，如果論文沒有好的結(jié)構(gòu)，內(nèi)容再好也沒有用。論文前面的結(jié)構(gòu)一般都不會(huì)變，后面可以按照實(shí)際情況來安排，省略的部分可以有結(jié)果說明、靈敏度分析、其他模型、模型擴(kuò)展、優(yōu)缺點(diǎn)分析等等，多看些優(yōu)秀論文就知道還有哪些形式了。附

錄可以貼一些算法流程圖或比較大的結(jié)果或圖表等等。分析問題要認(rèn)真

一般競(jìng)賽題目自己肯定沒有見過，而且我發(fā)現(xiàn)近些年來的賽題都不是書上哪個(gè)模型可以直接套成功的，很多根本就沒有固定的模型可以參考，所以分析問題不是一個(gè)去找書本的過程，依賴書本就意味著自己的思想被束縛起來。可以完全按照自己的分析去完成，平時(shí)練習(xí)的時(shí)候?qū)W習(xí)的是一種方法，通過以前學(xué)到的方法來解決，不是套用書本來解決，沒有模型套怎么辦，只有靠自己去實(shí)際分析。我估計(jì)在前面說的五點(diǎn)也許會(huì)有三分之一的隊(duì)可以做到，而且可以做的很好，但是這一點(diǎn)上就需要真本事了，平時(shí)多努力，比賽發(fā)揮正常，這一點(diǎn)做好是沒有問題的。編程求解是重要手段

美國(guó)競(jìng)賽時(shí)，美國(guó)學(xué)生中的論文很多是編程數(shù)據(jù)的說明，比如99 年

A 題行星撞地球那題，他們也能夠模擬出撞擊后果，這對(duì)我們來說簡(jiǎn)直是不可思議的。美國(guó)學(xué)生實(shí)踐能力較強(qiáng)，而中國(guó)學(xué)生擅長(zhǎng)理論分析，所以我把編程放在了分析的后面是有中國(guó)特色的。數(shù)學(xué)建模競(jìng)賽特別強(qiáng)調(diào)計(jì)算機(jī)編程解決實(shí)際問題的能力，最近幾年尤其強(qiáng)調(diào)，編程方面的能力不是一朝一夕可以練成的，需要長(zhǎng)期刻苦的訓(xùn)練，常用的工具有MATLAB、Mathematica、C/C++ 等等，一個(gè)人只需要會(huì)一門語言就行了，但需要精通它。比如要畫柱狀圖該怎么做，要用Floyd 算法怎么辦，賽前不準(zhǔn)備是沒有辦法在比賽中很好運(yùn)用的，因此每個(gè)常用的算法都自己去編程實(shí)現(xiàn)一下。模型的假設(shè)與模型的建立

評(píng)委看完摘要后緊接著就是看模型假設(shè)了，有一個(gè)萬能的方法就是可以抄題目中可以作為假設(shè)的幾句話，這樣會(huì)給人留下好的印象，畢竟說明你審題了。但不能全抄，要加上自己的一些假設(shè)。一般假設(shè)用文字描述就行了，最好不要太具體了，一些重要參數(shù)不要被定死只能取某些值，否則會(huì)讓人感覺論文的局限性較強(qiáng)。模型的建立是根據(jù)你對(duì)問題分析而來的，提出的數(shù)學(xué)符號(hào)和建立模型最好要比較接近，在同一頁(yè)最好，以便評(píng)委可以對(duì)照符號(hào)來看，數(shù)學(xué)公式要嚴(yán)謹(jǐn)，推導(dǎo)要嚴(yán)密，這些都反映了參賽者的數(shù)學(xué)素質(zhì)和能力，即使你推導(dǎo)不對(duì)，別人看到你的陣勢(shì)也首先會(huì)誤以為你是對(duì)的。那么多的試卷，評(píng)委不可能順著你的公式一直推下去，但你要寫得顯得有數(shù)學(xué)修養(yǎng)才行。9 圖文表并貌可以增色

我聽說一個(gè)不確切的信息是評(píng)委老師喜歡用MATLAB 編程的論文，不知道有沒有這回事，但這說明了老師需要看一個(gè)具有圖或表在其中的論文，一篇如果像政治書那樣寫的論文估計(jì)沒有人會(huì)對(duì)它感興趣的，尤其是科技論文。MATLAB 編程之所以受到青睞是因?yàn)镸ATLAB 提供的圖形處理能力很強(qiáng)大。圖表的說明性特別強(qiáng)，如果結(jié)論有很多數(shù)據(jù)的話，最好做成圖表的形式加以說明，會(huì)令你的論文更有說服力，也更容易受到評(píng)委的好評(píng)。其他

其他內(nèi)容還是有很多的，說也說不完，挑幾個(gè)重要的講。比如不要上網(wǎng)討論，網(wǎng)上的人水平參差不齊，你不知道誰是對(duì)的，而且很多人想

得獎(jiǎng)，不會(huì)告訴你正確的，反而騙你說相反的，有時(shí)真理往往掌握在少數(shù)人手里。還有就是論文寫作中靈敏度分析不要寫太多，大致說明一下就可以了，不要喧賓奪主。最后想到的就是要使用數(shù)學(xué)公式編輯器來寫論文，不要用什么上下標(biāo)來表示，論文字體用小四，分標(biāo)題用四號(hào)黑體等等。

第三篇：2014全國(guó)大學(xué)生數(shù)學(xué)建模競(jìng)賽

嫦娥三號(hào)軟著陸軌道設(shè)計(jì)與控制策略

摘要

隨著月球探測(cè)任務(wù)的發(fā)展,未來月球探測(cè)考察目標(biāo)將主要是 復(fù)雜地形特性的高科學(xué)價(jià)值區(qū)域。為了能夠安全地在這些遍布巖石、的區(qū)域內(nèi)完成高精度軟著陸，這就要求導(dǎo)航和控制系統(tǒng)具有較強(qiáng)的自主性和實(shí)時(shí)性。本文針對(duì)最終著陸段安全、精確的需求，對(duì)月球軟著陸導(dǎo)航與控制方法進(jìn)行較深入研究，主要內(nèi)容包括：

首先，提出一種基于單幀圖像信息的障礙檢測(cè)方法。該方法根據(jù)著陸區(qū)內(nèi)障礙成像的特點(diǎn)，通過匹配相應(yīng)的陰影區(qū)與光照區(qū)完成對(duì)巖石、彈坑的檢測(cè)，利用圖像灰度方差對(duì)粗糙區(qū)域進(jìn)行提取：在檢測(cè)出故障信息的基礎(chǔ)上，選取安全著陸點(diǎn)以保證軟著陸任務(wù)的成功。

其次，給出一種基于矢量觀測(cè)信息的自主光學(xué)導(dǎo)航方法。該方法利用光學(xué)相機(jī)和激光測(cè)距儀測(cè)量值構(gòu)建著陸點(diǎn)相對(duì)著陸器的矢量信息，結(jié)合著陸器的姿態(tài)信息確定著陸器的位置。為了消除測(cè)量噪聲帶來的干擾，利用擴(kuò)展Kalman濾波理論設(shè)計(jì)了導(dǎo)航濾波器。

再次，提出一種李雅普諾夫函數(shù)障礙規(guī)避制導(dǎo)方法。該方法通過對(duì)狀態(tài)函數(shù)、危險(xiǎn)地形勢(shì)函數(shù)的設(shè)計(jì)，以滿足平移過程中減低障礙威脅與精確定點(diǎn)著陸器，設(shè)計(jì)PWPF（調(diào)頻調(diào)寬）調(diào)節(jié)器實(shí)現(xiàn)定推理等效變推力控制效果。

最后，針對(duì)采用變推力主發(fā)動(dòng)機(jī)的月球著陸器，提出一種垂直軟著陸控制方法。該方法采用標(biāo)稱控制與閉環(huán)控制相結(jié)合的方式，規(guī)劃標(biāo)稱軌跡以保證著陸器到達(dá)著陸點(diǎn)時(shí)其下降速度、加速度亦為零，設(shè)計(jì)閉環(huán)控制器產(chǎn)生附加控制量消除初始偏差、著陸器質(zhì)量變化的干擾，以保證著陸器沿標(biāo)稱軌跡到達(dá)著陸點(diǎn)。

本文分別對(duì)所提出的最終著陸段導(dǎo)航與控制方法進(jìn)行數(shù)學(xué)仿真以驗(yàn)證個(gè)方法的可行性。仿真結(jié)果表明，本文多給出導(dǎo)航方法能夠達(dá)到較高的性能指標(biāo)，滿足在危險(xiǎn)區(qū)域?qū)崿F(xiàn)高精度軟著陸的需要。

關(guān)鍵詞： 月球軟著陸；自主導(dǎo)航與控制；障礙檢測(cè)；規(guī)避制導(dǎo)；適量測(cè)量

一、問題重述

嫦娥三號(hào)于2013年12月2日1時(shí)30分成功發(fā)射，12月6日抵達(dá)月球軌道。根據(jù)計(jì)劃，嫦娥三號(hào)將在北京時(shí)間12月14號(hào)在月球表面實(shí)施軟著陸。嫦娥三號(hào)如何實(shí)現(xiàn)軟著陸以及能否成功成為外界關(guān)注焦點(diǎn)。嫦娥三號(hào)在著陸準(zhǔn)備軌道上的運(yùn)行質(zhì)量為2.4t，其安裝在下部的主減速發(fā)動(dòng)機(jī)是目前中國(guó)航天器上最大推力的發(fā)動(dòng)機(jī)，能夠產(chǎn)生1500N到7500N的可調(diào)節(jié)推力，進(jìn)而對(duì)嫦娥三號(hào)實(shí)現(xiàn)精準(zhǔn)控制。其比沖（即單位質(zhì)量的推進(jìn)劑產(chǎn)生的推力）為2940m/s，可以滿足調(diào)整速度的控制要求。在四周安裝有姿態(tài)調(diào)整發(fā)動(dòng)機(jī)，在給定主減速發(fā)動(dòng)機(jī)的推力方向后，能夠自動(dòng)通過多個(gè)發(fā)動(dòng)機(jī)的脈沖組合實(shí)現(xiàn)各種姿態(tài)的調(diào)整控制。嫦娥三號(hào)的預(yù)定著陸點(diǎn)為19.51W，44.12N，海拔為-2641m。嫦娥三號(hào)將在近月點(diǎn)15公里處以拋物線下降，相對(duì)速度從每秒1.7公里逐漸降為零。整個(gè)過程大概需要十幾分鐘的時(shí)間。在距月面100米處時(shí)，嫦娥三號(hào)要進(jìn)行短暫的懸停，掃描月面地形，避開障礙物，尋找著陸點(diǎn)。之后，嫦娥三號(hào)在反推火箭的作用下繼續(xù)慢慢下降，直到離月面4米高時(shí)再度懸停。此時(shí)，關(guān)掉反沖發(fā)動(dòng)機(jī)，探測(cè)器自由下落。

嫦娥三號(hào)在高速飛行的情況下，要保證準(zhǔn)確地在月球預(yù)定區(qū)域內(nèi)實(shí)現(xiàn)軟著陸，關(guān)鍵問題是著陸軌道與控制策略的設(shè)計(jì)。其著陸軌道設(shè)計(jì)的基本要求：著陸準(zhǔn)備軌道為近月點(diǎn)15km，遠(yuǎn)月點(diǎn)100km的橢圓形軌道；著陸軌道為從近月點(diǎn)至著陸點(diǎn)，其軟著陸過程共分為6個(gè)階段，分別為著陸準(zhǔn)備軌道、主減速段、快速調(diào)整段、粗避障段、精避障段、緩速下降階段，要求滿足每個(gè)階段在關(guān)鍵點(diǎn)所處的狀態(tài)；盡量減少軟著陸過程的燃料消耗。

根據(jù)上述的基本要求，請(qǐng)你們建立數(shù)學(xué)模型解決下面的問題：

（1）確定著陸準(zhǔn)備軌道近月點(diǎn)和遠(yuǎn)月點(diǎn)的位置，以及嫦娥三號(hào)相應(yīng)速度的大小與方向。（2）確定嫦娥三號(hào)的著陸軌道和在6個(gè)階段的最優(yōu)控制策略。

（3）對(duì)于你們?cè)O(shè)計(jì)的著陸軌道和控制策略做相應(yīng)的誤差分析和敏感性分析。

二、問題分析

對(duì)于問題一：

嫦娥三號(hào)從15公里左右的高度下降到月球表面，在這一過程中不考慮月球表面太陽(yáng)風(fēng)的影響，忽略月球的自轉(zhuǎn)速度引起的科氏力的影響，由于下降時(shí)間比較短也不考慮太陽(yáng)、地球?qū)︽隙鹑?hào)的攝動(dòng)影響，嫦娥三號(hào)水平速度要從1.692km/s降為0m/s由于3000m處時(shí)嫦娥三號(hào)已經(jīng)基本位于著陸點(diǎn)上方，所以此時(shí)假設(shè)在3000m處的速度只存在豎直向下的速度而不存在水平分速度，因?yàn)榻德錅p速時(shí)間比較短只有垂直于月面的方向運(yùn)動(dòng)才能實(shí)現(xiàn)，所以在確定著陸點(diǎn)位置和著陸軌跡時(shí)應(yīng)當(dāng)考慮燃料最優(yōu)情況下推力最大，方向自由的方法即取F?7500N建立主減速段動(dòng)力學(xué)模型。

三、符號(hào)說明

四、模型假設(shè)

對(duì)于問題一：

忽略月球的自傳和太陽(yáng)、地球?qū)︽隙鹑?hào)衛(wèi)星的引力攝動(dòng) 月球近似為一個(gè)質(zhì)量均勻的標(biāo)準(zhǔn)球體 將嫦娥三號(hào)是為一個(gè)質(zhì)點(diǎn)

主減速忽略動(dòng)作調(diào)整所產(chǎn)生的燃料消耗段不考慮太陽(yáng)風(fēng)的影響

五、模型建立與求解

5.1問題一的建模與求解 解法一： 假設(shè)嫦娥三號(hào)在t時(shí)刻在遠(yuǎn)月點(diǎn)開始緩慢下降，在n時(shí)刻到達(dá)近月點(diǎn)，整個(gè)過程遵循開普勒第三定律，即

v0?0

在t時(shí)刻有：v1?2??R1????? ??R0?R0?R1?r0 R0?r1?r2 其中v1：遠(yuǎn)月點(diǎn)速度

v2：近月點(diǎn)速度

R0：遠(yuǎn)月點(diǎn)月心距

R1：近月點(diǎn)月心距（已知月球的半徑為1738千米）

R0?1738?100?1838km

R1?1738?15?1753km 在t1時(shí)刻處v2? k?2??R1??? ?R0?R0?R1??R0?0.512k?0.488 R0?R1利用能量平衡式求得近地點(diǎn)速度為

2?0.512?49012（）?1.692km/s（沿切線方向）v2?，比當(dāng)?shù)氐沫h(huán)境速度17531.672km/s大?vk?0.0196km/s，徑向速度vk?0。

1同理解得v1?1.6139km/s（沿切線方向）

vri?0

解得主減速段動(dòng)力學(xué)模型的建立：

根據(jù)題意，在橫向飛行的水平距離遠(yuǎn)遠(yuǎn)小于月球半徑的平均值，所以可以將整個(gè)減速段過程簡(jiǎn)化為水平和豎直方向運(yùn)動(dòng)方程，根據(jù)牛頓第二定律、速度計(jì)算公式有：

ax?Tx may?tTymTxt?a

?1.692km/s ?m?0?Qdt0??Ty???a?dt?57m/s t?0??m??Qdt?0??t?T22x?Ty2??7500N

v2?2at?S

運(yùn)用matlab編程解得S?451810.4m； 其中 ax：水平方向加速度

ay：豎直方面加速度

a：月球表面重力加速度a? Tx：推力的水平方向分力

Ty：推力的豎直方向分力

t:主減速段時(shí)間

S:嫦娥三號(hào)主減速段水平位移

Q:嫦娥三號(hào)發(fā)動(dòng)機(jī)燃料秒消耗率

根據(jù)已知資料得到嫦娥三號(hào)著陸過程中緯度改變，經(jīng)度基本不變，月球赤緯和地球緯度一樣也分為南北各90個(gè)分度,又因?yàn)樵虑驑O區(qū)半徑為1735.843km，所以每一個(gè)緯度的豎直高度差為19.2871

4g 6千米。即近月點(diǎn)位置坐標(biāo)為?19.0464W,28.9989N?海拔15km,遠(yuǎn)月點(diǎn)位置坐標(biāo)為?160.9536E,28.9989S?海拔100km。

解法2:軌跡方程法。

眾所周知,太陽(yáng)系中的八大行星都在按照各自的橢圓軌道繞太陽(yáng)進(jìn)行公轉(zhuǎn),太陽(yáng)位于橢圓的一個(gè)焦點(diǎn)上,行星的運(yùn)動(dòng)遵循開普勒三定律,筆者發(fā)現(xiàn),在各類物理競(jìng)賽中,常會(huì)涉及到天體運(yùn)動(dòng)速度的計(jì)算,本文擬從能量和行星運(yùn)動(dòng)的軌跡方程兩個(gè)不同的角度來探索行星在近日點(diǎn)和遠(yuǎn)日點(diǎn)的速度。

該解法的指導(dǎo)思想是對(duì)橢圓的軌跡方程求導(dǎo),并結(jié)合一般曲線的曲率半徑通式求出近日點(diǎn)和遠(yuǎn)日點(diǎn)的曲率半徑表達(dá)式,然后利用萬有引力提供向心力列方程求解。如圖1所示,橢圓的軌跡方程為

x2y2?2?1 ?5? 2ba將?5?式變形為

a2x2?b2y2?a2b2 ?6?

根據(jù)隱函數(shù)的求導(dǎo)法則將?6?式對(duì)x求導(dǎo)有

2a2x?2b2yy??0 ?7? 即

a2xy???2 ?8?

by將?7?式再次對(duì)x求導(dǎo)得

2a2?2b2(y?y??yy??)?0 ?9? 將?8?、?9?兩式聯(lián)立得

a2b2y2?a4x2 ?10? y???-43by根據(jù)曲率半徑公式有 r?(1?y?)?11? ??y122 將?8?、?10?、?11?式聯(lián)立并將A點(diǎn)坐標(biāo)A(0，a)代入可得A點(diǎn)的曲率半徑為

b2RA? ?12?

a根據(jù)橢圓的對(duì)稱性,遠(yuǎn)日點(diǎn)B的曲率半徑為

b2RB?RA? ?13?

a 由于在A、B兩點(diǎn)行星運(yùn)行速度方向與萬有引力方向垂直,萬有引力只改變速度方向,并不改變速度大小,故分別根據(jù)萬有引力提供向心力得

GMmmvA ?14? ?(a?c)2RAGMmmvB ?15? ?2(a?c)RB將?13?至?15?式聯(lián)立可得 22vA?bGMbGM，vB? ??a?caa?ca

5.2問題二的建模與求解 模型一：動(dòng)力學(xué)模型

典型的月球軟著陸任務(wù)中,探測(cè)器一般首先發(fā)射到100km的環(huán)月停泊軌道,然后根據(jù)所選定的著陸位置,在合適的時(shí)間給著陸器一個(gè)有限脈沖,使得著陸器轉(zhuǎn)入近月點(diǎn)(在著落位置附近)為15km,遠(yuǎn)月點(diǎn)為100km的月球橢圓軌道,這一階段稱為霍曼轉(zhuǎn)移段。當(dāng)著陸器運(yùn)行到近月點(diǎn)時(shí),制動(dòng)發(fā)動(dòng)機(jī)開始工作,其主要任務(wù)是抵消著陸器的初始動(dòng)能和勢(shì)能,使著陸器接觸地面時(shí),相對(duì)月面速度為零,即實(shí)現(xiàn)所謂的軟著陸,這一階段稱為動(dòng)力下降段。著陸器的大部分燃料都是消耗在此階段,所以月球軟著陸軌跡優(yōu)化主要是針對(duì)動(dòng)力下降段這一階段。由于月球表面附近沒有大氣,所以在飛行器的動(dòng)力學(xué)模型中沒有大氣阻力項(xiàng)。而且從15km左右的軌道高度軟著陸到月球表面的時(shí)間比較短,一般在幾百秒的范圍內(nèi),所以諸如月球引力非球項(xiàng)、日月引力攝動(dòng)等影響因素均可忽略不計(jì),所以這一過程可以在二體模型下描述。其示意圖如圖1所示,其中o為月球質(zhì)心,x軸方向?yàn)橛稍滦闹赶蛑懫鞯某跏嘉恢?y軸方向?yàn)槌跏嘉恢弥懫魉俣确较颉?/p>

圖 1 月球軟著陸極坐標(biāo)系

其動(dòng)力學(xué)方程如下: r??v ????

v??(F/m)sin???/r?r

2?2 ????((F/m)cos??2v?)/r

m???F/ISP

在上式中r為著陸器與月心距離,v為著陸器徑向速度,?為著陸器極角,?為著陸器極角角速度,?為月球引力常數(shù),F著陸器制動(dòng)發(fā)動(dòng)機(jī)推力,m為著陸器質(zhì)量,?為制動(dòng)發(fā)動(dòng)機(jī)推力方向角,其定義為F與當(dāng)?shù)厮椒较驃A角,ISP為制動(dòng)發(fā)動(dòng)機(jī)比沖。根據(jù)動(dòng)力下降段的起點(diǎn)位置可以確定動(dòng)力學(xué)方程初始條件,由于起點(diǎn)處于霍曼轉(zhuǎn)移軌道的近地點(diǎn),故其初始條件為: r0?rp

?0?0

v0?0 ?0?1rp?rp(2ra)ra?rp其中rp和ra分別為霍曼轉(zhuǎn)移段的近地點(diǎn)半徑和遠(yuǎn)地點(diǎn)半徑。

終端條件為實(shí)現(xiàn)軟著陸, 即

rf?R

vf?0

?f?0

其中R為月球半徑,終端條件中對(duì)終端極角?f及終端時(shí)間tf無約束。

優(yōu)化變量為制動(dòng)發(fā)動(dòng)機(jī)推力方向角?(t)。

優(yōu)化的性能指標(biāo)為在滿足上述初始條件和終端條件的前提下, 使著陸過程中燃料消耗最少,即

J??m(t)dt

t0f設(shè)計(jì)主減速段制導(dǎo)控制律 2動(dòng)力下降段燃料最優(yōu)精確著陸問題描述 2.1 燃料最優(yōu)精確著陸問題

著陸器運(yùn)動(dòng)方程：考慮采用變推力發(fā)動(dòng)機(jī)情況，有

r?v

.v?g?a

(1)

a?Tmm??aT..其中r?[rhrxry]T，v?[vhvxvy]T分別表示著陸器相對(duì)期望著陸點(diǎn)的位置和速度矢量；T為推力器提供的推力矢量，幅值為 T，對(duì)應(yīng)控制加速度矢量 a；g為火星的重力加速度矢量，此處認(rèn)為是常值；m為著陸器質(zhì)量，對(duì)應(yīng)推力器質(zhì)量排除系數(shù)?。指標(biāo)函數(shù)：考慮燃料消耗

min(m0?mf)???min?0fTdt

(2)邊界條件：即初始條件和終端條件

r(0)?r0,v(0)?v0,m(0)?m0,r(tf)?v(tf)?[000]

(3)控制約束：考慮發(fā)動(dòng)機(jī)一旦啟動(dòng)不能關(guān)閉，存在最大和最小推力約束

0?T1?T?T

2(4)狀態(tài)約束：為避免在著陸前撞擊到火星地表，需確保整個(gè)下降段位于火星地平面以上，即

rh?0

（5）進(jìn)一步地，若著陸區(qū)域附近表面崎嶇不平，僅僅確保地表約束不能滿足需求時(shí)，可以考慮下降傾角約束，即將著陸器下降軌線約束到以著陸點(diǎn)為頂點(diǎn)的圓錐體內(nèi)

2.2 等效后燃料最優(yōu)精確著陸問題 定義等效變換變量

Ttrx2?ry2rh?tan?alt

（6）

u?a?T

?m

（7）

??Tmz?lnm??等效著陸器運(yùn)動(dòng)方程： ?.??r??0I3.?.??

y??v??00??.??00?z????其中p?[u?T0??r??0?v???I0?????30??z????0?7*7?0??u?g??0????Acy?Bc(p?g4)

（8）????????],g4?[gTT?0]T

t指標(biāo)函數(shù)：

min?0f?(t)dt

（9）

邊界條件：同式(3)。

控制約束：由文獻(xiàn)[10]可知，控制約束(4)可等效表示為

u??1T1e?z0[1?(z?z0)?(z?z0)2]???T2e?z0[1?(z?z0)]

（10）（11）

2狀態(tài)約束：地表約束同式(5)，傾角約束(6)可等效表示為

T

Sy?cy?0

（12）

其中

?0100000?S???

0010000??c???tan?alt

T000000?

3.燃料最優(yōu)精確著陸問題的離散化及變換 3.1 等效燃料最優(yōu)精確著陸問題的離散化

首先將整個(gè)飛行時(shí)間均分成 n 段（對(duì)應(yīng) n +1 個(gè)點(diǎn)），每段步長(zhǎng)為?t，離散化后的著陸器運(yùn)動(dòng)方程為：yk?1?Ayk?B(pk?g4)

其中A?R7?7,B?R7?4分別為離散系統(tǒng)的系統(tǒng)矩陣和輸入矩陣

12A?e?tAc?I3??tAc??tAc??

2?t?t112B??e??t?s?AcBcds??esAcds?Bc??tBc??tBc??t2Bc??

0026其中I3為三階單位陣。

有系統(tǒng)性質(zhì)可知，整個(gè)控制時(shí)域內(nèi)系統(tǒng)狀態(tài)滿足 y3?Ay2?B?p2?g4??A3y0?A2B?p0?g4??AB?p1?g4??B?p2?g4??yn?Ayn?1?B?pn?1?g4??Any0?An?1B?p0?g4????AB?pn?2?g4??B?pn?1?g4?y1?Ay0?B?p0?g4?y2?Ay1?B?p1?g4??A2y0?AB?p0?g4??B?pn?2?g4??B?p1?g4?

為表達(dá)方便，令

?y0??p0???0??A0??y??p?????1??1??1??1??A? ,p??p2?，????2???A2? Y??y2?????????????????????n?????yn??7?n?1??1?pn??4?n?1??1??n????A?7?n?1??7??0??0????B?1??????AB???2???2??3??AB????????n?1?A??n????則(15)可等價(jià)于

0???0??0?????B?0?1???????2??AB?B?B000???????2? ?ABB00???3??A?AB?B??????????0????n?1???A???AB?B?A2BABB????n????7?n?1??4?n?1???000000Y??y0??p??g4

分別定義如下常值矩陣：

最終可得離散化后的燃料最優(yōu)化問題如下： 指標(biāo)函數(shù)：式(9)可表示為

邊界條件：式(3)可表示為

控制約束：式(10)和式(11)分別可表示為

狀態(tài)約束：式(5)和式(12)分別可表示為

含有 p個(gè)線性約束和 q個(gè)二階錐約束的最優(yōu)化問題的標(biāo)準(zhǔn)形式為 指標(biāo)函數(shù)

min(?Tx)滿足約束

DTx?f?0Ax?ci?b?dinTiTi

（k=1,?，n）

n*pp其中x?R為待優(yōu)化向量，??R，線性約束參數(shù)D?R,f?R，二階錐約束參數(shù)維數(shù)n(Ai,bi,ci,di)由相應(yīng)約束確定

則式(17)～式(23)可最終轉(zhuǎn)換為如下最優(yōu)化問題： 指標(biāo)函數(shù)：min(vpp)滿足：

初值約束:MxΨ0p?Mx(Ψ0y0)?A0g4?r0末值約束:MxΨ0p?Mx(Ψ0y0)?A0g4控制約束:Murkp?v?rkp 控制上限:?(vzΨk?TT?TTv0?T?0

?0

T1vr)p?1?vTz(Φky0?Akg4)?z0,z?0 ?z0?kT2e 控制下限：

4數(shù)值仿真結(jié)果與分析本節(jié)以某火星著陸器為例，計(jì)算了典型初始條件下滿足各種約束的燃料最優(yōu)精確著陸軌跡。其中探測(cè)器各參數(shù)分別取為：m0?2000kg,g?[?3.711400]ms2,c?2kms,T1?1.3kN,T2?13kN.。著陸器初始位置矢量r0= [1500,-600, 800] m，初始速度矢量v0= [-30, 10, 40]m/s，傾角?alt=86°。二階錐優(yōu)化問題可以通過大量免費(fèi)的優(yōu)化工具求解，如 CSDP、DSDP、OpenOpt、SeDuMi、SDPA、SDPLR等。本文選用 SDPT3 進(jìn)行計(jì)算，通過執(zhí)行線性搜索確定燃料最優(yōu)下降時(shí)間tf為 43s，圖 1 給出了相應(yīng)的最優(yōu)著陸軌跡、下降速度、加速度、控制推力、推力幅值以及探測(cè)器質(zhì)量變化曲線。

由優(yōu)化結(jié)果可以看出，探測(cè)器在給定時(shí)間飛行并軟著陸到指定位置，且在整個(gè)下降過程始終與火星地表保持一定的安全距離，驗(yàn)證了下降傾角約束的有效性。其推力幅值曲線呈現(xiàn)“最大-最小-最大”的最優(yōu)控制形式，不過為了保持發(fā)動(dòng)機(jī)始終處于點(diǎn)火狀態(tài)，在中間段對(duì)應(yīng)最小推力約束，這與文獻(xiàn)中的分析結(jié)論一致。此外，通過利用如 TOMLAB 等商業(yè)最優(yōu)控制軟件進(jìn)行復(fù)核計(jì)算，也驗(yàn)證了此計(jì)算結(jié)果的燃料最優(yōu)性能。

*

圖 1 給定初始條件下火星著陸器動(dòng)力下降段燃料最優(yōu)計(jì)算結(jié)果

需要注意到，此燃料最優(yōu)軌跡的獲取對(duì)著陸器的實(shí)時(shí)在線計(jì)算性能提出了較高的要求，經(jīng)測(cè)試，無論使用何種優(yōu)化工具，計(jì)算給定飛行任務(wù)時(shí)間的最優(yōu)軌跡均需數(shù)秒，而全局最優(yōu)則需要數(shù)十秒甚至更長(zhǎng)，這在實(shí)際任務(wù)中是不允許的。因此，可行的方案是通過在地面計(jì)算大量的燃料最優(yōu)軌跡，并尋找規(guī)律，選取關(guān)鍵路徑點(diǎn)狀態(tài)存儲(chǔ)到著陸器計(jì)算機(jī)中，通過在線查表或者在利用對(duì)計(jì)算量要求較小的反饋制導(dǎo)律完成安全著陸任務(wù)。

因此，為了研究探測(cè)器燃料最優(yōu)軌跡特性，選取相同的探測(cè)器參數(shù)，暫不考慮推力器最小幅值約束和傾斜角約束（但考慮地表約束），固定初始高度為 1500m，初始位置水平方向從-8000m 到 8000m 內(nèi)取值，分別選取各種不同的初始速度，可得燃料最優(yōu)精確著陸軌跡簇如圖 2 所示。

圖 2 各種不同初始速度對(duì)應(yīng)的火星著陸器動(dòng)力下降段燃料最優(yōu)軌跡簇

1)對(duì)任意探測(cè)器初始位置，特定初始速度對(duì)應(yīng)的燃料最優(yōu)著陸軌跡在末端必然收斂到一個(gè)固定的近似圓錐體內(nèi)。

2)取決于探測(cè)器初始位置和速度的關(guān)系，燃料最優(yōu)軌跡有兩種形式：S 型和 C 型，其中 S 型主要對(duì)應(yīng)于期望著陸點(diǎn)位置水平距離較大情況。3)當(dāng)探測(cè)器初始水平速度為零時(shí)，圓錐體軸線垂直于火星地表，所有最優(yōu)軌線關(guān)于該軸線中心對(duì)稱。4)初始速度的大小也直接影響到任務(wù)的可靠性，因此需要在超聲速進(jìn)入段和降落傘減速段將著陸器速度下降到合理范圍內(nèi)。

上述結(jié)論對(duì)上注探測(cè)器關(guān)鍵點(diǎn)的選取有著較強(qiáng)的指導(dǎo)意義，比如基于最優(yōu)軌線的斜率對(duì)路徑點(diǎn)合并、基于最優(yōu)軌線簇的對(duì)稱性對(duì)上注軌線進(jìn)行等效延伸、或者嘗試僅將 S 型和 C 型的轉(zhuǎn)折點(diǎn)作為路徑點(diǎn)等，這樣可以大大降低探測(cè)器自主存儲(chǔ)與計(jì)算需求，進(jìn)而有效提升任務(wù)的可靠性。重力轉(zhuǎn)彎軟著陸過程

對(duì)于最終著陸點(diǎn)，假設(shè)探測(cè)器的下降軌跡在一平面內(nèi)，且月球引力場(chǎng)為垂直于月面XY的均勻引力場(chǎng)，引力加速度g沿-Z，如圖1所示，制動(dòng)推力方向沿探測(cè)器的本體軸z。重力轉(zhuǎn)彎軟著陸過程中探測(cè)器質(zhì)心動(dòng)力學(xué)方程可表示為

上式中各變量的物理意義如圖1中所示，其中m>0為探測(cè)器質(zhì)量；k>0為制動(dòng)發(fā)動(dòng)機(jī)比沖；u表示制動(dòng)發(fā)動(dòng)機(jī)的秒耗量

可通過一定的機(jī)構(gòu)加以調(diào)節(jié)，故作為軟著陸問題的控制變量。假定制動(dòng)發(fā)動(dòng)機(jī)的最大推力與初始質(zhì)量比大于月面引力加速度，并且制動(dòng)推進(jìn)系統(tǒng)能夠在一定的初始條件下將探測(cè)器停止月面上。

重力轉(zhuǎn)彎過程中，探測(cè)器的高度、速度和姿態(tài)角度可由雷達(dá)高度表、多普勒雷達(dá)及慣性儀表測(cè)得。令軟著陸初始條件探測(cè)器到達(dá)月面時(shí)速度減小到給定的值，故終端條件自由。軟著陸燃耗最優(yōu)問題的描述 對(duì)于最終著陸段，可假設(shè)

為一小角度。由此可將系統(tǒng)方程（1）化簡(jiǎn)為

要設(shè)計(jì)制導(dǎo)律實(shí)現(xiàn)軟著陸，就是使

著陸時(shí)間

對(duì)于月球軟著陸的燃耗最優(yōu)控制問題，其性能指標(biāo)可表示為

對(duì)于系統(tǒng)（2）的軟著陸過程，燃耗最優(yōu)問題等價(jià)于著陸時(shí)間最優(yōu)問題，性能指標(biāo)為

在月球重力轉(zhuǎn)彎軟著陸過程中，如果存在一個(gè)推力控制程序?qū)⑻綔y(cè)器從初始條件轉(zhuǎn)移到終端條件，并使性能指標(biāo)（3）或（4）式最大，則稱這個(gè)推力程序?yàn)檐浿懭己淖顑?yōu)或時(shí)間最優(yōu)制導(dǎo)律。根據(jù)pontryagin極大值原理，系統(tǒng)的哈密頓函數(shù)及其對(duì)u的偏導(dǎo)數(shù)為

使哈密頓函數(shù)（5）式達(dá)到極大地控制輸入u就是最優(yōu)控制，科表示為。

如果存在一個(gè)有限區(qū)間

則最優(yōu)控制u（t）取值不能由哈密頓函數(shù)確定。此時(shí)如果最優(yōu)解存在，則稱為奇異解，（8）式稱為奇異條件。

最優(yōu)制導(dǎo)問題的性質(zhì)：1）對(duì)于自治系統(tǒng)（2）的時(shí)間最優(yōu)控制問題，沿最優(yōu)軌跡其哈密頓函數(shù)滿足

將其對(duì)時(shí)間求導(dǎo)并將（2c）和（6c）式代入，得

另外，由于自由，根據(jù)橫截條件有3）根據(jù)（6a）式。又由（9）式可得T（t）=0，4）根據(jù)極大值原理，系統(tǒng)的狀態(tài)變量和共軛變量都是時(shí)間的連續(xù)可微函數(shù)，將切換函數(shù)對(duì)時(shí)間求導(dǎo)，利用（2），（6）式和性質(zhì)2）得 軟著陸最優(yōu)控制中奇異條件的分析

對(duì)于月球重力轉(zhuǎn)彎軟著陸問題，最優(yōu)制導(dǎo)律具有兩個(gè)很好的性質(zhì)。

定理一。月球重力轉(zhuǎn)彎軟著陸系統(tǒng)（2）的燃耗最優(yōu)制導(dǎo)或時(shí)間最優(yōu)制導(dǎo)問題不存在奇異條件。證明。用反證法，假設(shè)存在奇異條件，則在某個(gè)閉區(qū)間設(shè)，并由（5）式得

。根據(jù)反正假將（10）式兩邊對(duì)時(shí)間求導(dǎo)，并將（2）和（6）式代入化簡(jiǎn)得性質(zhì)2），并考慮到或者情形1.得

下面證明這兩種情形均與反證假設(shè)矛盾。根據(jù)式

及性質(zhì)2）可知，由性質(zhì)3）必有

根據(jù)

是時(shí)間t的斜率非零的線性函數(shù)，m和情形2.1）若定，根據(jù)橫截條件有在區(qū)間內(nèi)為常數(shù)。這與反證假設(shè)矛盾。

。下面再分三種情況進(jìn)行分析。

又因?yàn)?/p>

不與此時(shí)由（6b）式有反證假設(shè)矛盾。2）若盾。3），與反證假設(shè)矛又

因

為

因此有成立，這與

此時(shí)（10）式在上根據(jù)定理一，重力轉(zhuǎn)彎軟著陸的最優(yōu)制導(dǎo)律是一種開關(guān)（Bang-Bang）控制，只須控制發(fā)動(dòng)機(jī)開關(guān)，不需要調(diào)節(jié)推力的大小。

定理2.對(duì)于月球重力轉(zhuǎn)彎軟著陸過程，其開關(guān)控制器的最優(yōu)推力程序（7）最多進(jìn)行一次切換。

證明。只要證明最多只在一個(gè)時(shí)間點(diǎn)成立即可。軟著陸系統(tǒng)（2）在最優(yōu)推力控制程序（7）的作用下，按最后軌跡降落。由性質(zhì)3）知，為常數(shù)。根據(jù)性質(zhì)4），若嚴(yán)格單調(diào)，因而在上至多有一個(gè)零點(diǎn)，即至多進(jìn)行一次切換；若，則上為常數(shù)。由定理1，5 軟著陸最優(yōu)開關(guān)制導(dǎo)律

不可能在任何區(qū)間上成立，故必有既沒有切換點(diǎn)。

對(duì)于最優(yōu)推力控制程序（7），其切換函數(shù)中含有共軛變量，它是一個(gè)關(guān)于狀態(tài)變量的穩(wěn)式表達(dá)式。為實(shí)現(xiàn)實(shí)時(shí)制導(dǎo)，需求出關(guān)于狀態(tài)變量的切換函數(shù)來。

根據(jù)定理一和定理二，重力轉(zhuǎn)彎軟著陸最優(yōu)控制程序沒有奇異值狀態(tài)，并且在著陸過程中最多切換一次，其工作方式有4種：1）全開；2）全關(guān)；3）先開有關(guān)；4）先關(guān)后開。對(duì)于方式1）軟著陸起始點(diǎn)即是開機(jī)點(diǎn)；方式2），3）不能實(shí)現(xiàn)軟著陸；最后一種是通常情況下的最優(yōu)著陸方式，即探測(cè)器先做無制動(dòng)下降，然后打開發(fā)動(dòng)機(jī)軟著陸到月面。設(shè)開機(jī)時(shí)刻為到發(fā)動(dòng)機(jī)工作時(shí)間為

式，在區(qū)間

內(nèi)積分，并考慮

將（11）式中的對(duì)數(shù)按泰勒展開，忽略

并令

消掉T得到切換函數(shù)為

由切換函數(shù)（12）式可以看出，速度、位置的誤差和制動(dòng)發(fā)動(dòng)機(jī)推動(dòng)的將直接影響著陸的效果。一種方法是將終端高度從到達(dá)月面時(shí)實(shí)現(xiàn)軟著陸設(shè)置為離月面還有幾米時(shí)實(shí)現(xiàn)軟著陸。另一種方法是考慮制動(dòng)過程由一個(gè)主發(fā)動(dòng)機(jī)和一組小推力發(fā)動(dòng)機(jī)共同完成，通過調(diào)整開啟的小發(fā)動(dòng)機(jī)的數(shù)量，來實(shí)現(xiàn)變推力降落。具體地，令切換函數(shù)為

式中各符號(hào)的含義如圖2所示

關(guān)機(jī)點(diǎn)可取為2m,可取為20m，可取為1m/s。為實(shí)現(xiàn)著陸的最優(yōu)性，減速度

取為

其中T如（12）式中所示，m0為探測(cè)器的初始質(zhì)量。

圖三為最優(yōu)著陸過程與其改進(jìn)方法按圖2降落的次優(yōu)著陸過程的對(duì)比圖。由此圖中可看出，改進(jìn)方法提高了著陸的安全性，當(dāng)探測(cè)器的初始質(zhì)量mo=350kg，發(fā)動(dòng)機(jī)著陸過程多消耗燃料2.2kg。

時(shí)，改進(jìn)方法比最優(yōu)

（a）

（b）

問題三 協(xié)方差分析方法的基本原理 對(duì)于如下非線性函數(shù)關(guān)系

y?f?x1,x2??xn?（1）

可以使用一階泰勒級(jí)數(shù)展開對(duì)其進(jìn)行線性化，有

y??y?f??f?f?x1????xn???x1?xn?（2）?x1?xn其中，??x1??xn?為x1??xn的高階項(xiàng)。從而得到線性化方程

?y???f?xi（3）i?1?xin或表示為

?Y?P?X(4)

這里 P 是偏導(dǎo)數(shù)矩陣: Pi??f（5）?xi若自變量?x1???xn是隨機(jī)變量，則線性化方程的函數(shù)?y的協(xié)方差矩陣為：

E?Y??YT?EP?X?XTPT?PE?X?XTPT(6)即 ??????Cy?PCXPT(7)式中Cx是自變量的協(xié)方差矩陣；Cy是函數(shù)?Y的協(xié)方差矩陣。

協(xié)方差矩陣中對(duì)角線元素是方差，非對(duì)角線元素為協(xié)方差。顯然，只要求出傳遞矩陣 P ,便可確定源誤差與欲求量誤差之間的關(guān)系。若給定各種源誤差，如發(fā)動(dòng)機(jī)安裝誤差、敏感器測(cè)量誤差或發(fā)動(dòng)機(jī)推力和點(diǎn)火時(shí)間等誤差時(shí)，便可以分析其對(duì)目標(biāo)軌道誤差的影響以及對(duì)控制系統(tǒng)精度的影響，進(jìn)一步對(duì)各系統(tǒng)及元部件提出適當(dāng)?shù)木纫?。?jì)算向月飛行軌道誤差的協(xié)方差迭代方程

考慮到軌道參數(shù)的誤差之相對(duì)于軌道參數(shù)的標(biāo)稱值是小量，因此可以將軌道運(yùn)動(dòng)方程進(jìn)行線性化，從而得到能夠反映軌道參數(shù)偏差量的傳播關(guān)系的誤差方程。在應(yīng)用雙二體模型且在地球影響球范圍內(nèi)時(shí)，對(duì)軌道運(yùn)動(dòng)產(chǎn)生攝動(dòng)影響的各項(xiàng)，如月球引力攝動(dòng)、太陽(yáng)引力攝動(dòng)、大氣阻力攝動(dòng)和太陽(yáng)光壓攝動(dòng)等對(duì)誤差方程的影響很小，因此在誤差方程中將它們忽略掉。反映軌道位置和速度誤差的線性化方程如下：

?????v??r???g??（8）??v????r??rT?u???r，其中u?為地球引力常數(shù)。式中 g?r????3rr?rx2?ry2?rz2（9）

寫成狀態(tài)方程形式：

?????0I???r??r???????????(10)??v??G0???v??????????g式中 G??T

?r??0I???r?令F?????G0??,X????v?（11）

????則式（9）變?yōu)?/p>

??F?X(12)X下面推導(dǎo)矩陣 F 的表達(dá)式：

??g??u??G??T??T??3r??r?r?r?????u???u???r?r?T?3????3??T?r?r??r??r????u????u????u????u???r???3???3???3????3I3??rr?rr?r?y??z?r?????x???r（13）

式中 r x，r y 和 r z 是探測(cè)器在地心慣性坐標(biāo)系里的軌道位置坐標(biāo)。則G??u?3??T(I?rr)(14)332rr?rx2rxry?rx????T??2rr??ry??rxryrz???ryrxry??r??rzrxrzry?z??rxrz??ryrz?（15）2?rz??

將式（15）、（14）代入（10），得： ?0?0??02?-u?rx(1?32)F??r3r??3u?rxry?r5v??3u?rxrz?r5?

積分式（11），得到： 0003u?rxryr520003u?rxrzr53u?rzryr5210000ry-u?(1?3)32rr3u?rzryr5-u?rz(1?3)0r3r200?10??01??00?（16）

??00??00???

X??t??eF?tX?0?

（17）式中

(F?t)2(F?t)3(F?t)4(F?t)ne?I?F?t??????2!3!4!n!

（18）iN?t??Fi.()i!i?0F?t取前 6 階截?cái)?，即?/p>

eF?t??ti???F??i!??

（19）i?0??6i

得到計(jì)算誤差方程的迭代方程：

X?ti??t??eF?tX?ti?

（20）

eF?t相當(dāng)于式（4）中的 P 陣，由于誤差方程是時(shí)變方程，因此每一步迭代都需要重新計(jì)算 P 陣，計(jì)算 P 陣需要利用標(biāo)稱軌道參數(shù)數(shù)據(jù)。

進(jìn)一步根據(jù)式（7），得到協(xié)方差矩陣的迭代方程：

T

Ci?1?PCPiii

（21）向月飛行軌道誤差的協(xié)方差分析

引起軌道誤差的誤差源主要是導(dǎo)航誤差，包括位 置 誤 差 和 速 度 誤 差。其 中 ： 位 置 誤 差 ：?r??rx,?ry,?rz,?rx,?ry,?rz分別為在地心慣性坐標(biāo)系中 X 軸、Y 軸、Z 軸的分量。速度誤差：?v??vx,?vy,?vz,?vx,?vy,?vz分別是在地心慣性坐標(biāo)系 X 軸、Y 軸、Z 軸的分量。向月飛行軌道的初始軌道位置和速度誤差由運(yùn)載火箭的發(fā)射入軌精度決定，若探測(cè)器在飛行途中進(jìn)行軌道修正，則經(jīng)過軌道修正以后的軌道位置誤差將由導(dǎo)航誤差決定，速度誤差將由姿態(tài)誤差和制導(dǎo)誤差決定。

上述誤差決定了軌道誤差協(xié)方差分析的計(jì)算初始條件，表 1 給出了在不進(jìn)行中途軌道修正情況????下，在地心慣性坐標(biāo)系里，初始軌道位置誤差和初始速度誤差對(duì)軌道終點(diǎn)的位置和速度誤差的影響。圖 1 和圖 2 給出了在算例三中探測(cè)器從近地軌道入軌點(diǎn)開始至進(jìn)入月球軌道為止軌道位置的相應(yīng)的軌道位置和速度總誤差（3σ）的時(shí)間歷程。

表 1 初始軌道位置和速度誤差

對(duì)軌道終點(diǎn)誤差的影響

圖 1 軌道位置總誤差時(shí)間歷程（3σ）

圖 2 速度總誤差時(shí)間歷程（3σ）基于敏感系數(shù)矩陣的制導(dǎo)誤差分析

在月球軟著陸主制動(dòng)段，影響制導(dǎo)精度的誤差源主要有偏離標(biāo)準(zhǔn)飛行軌跡的初始條件誤差和導(dǎo)航與控制傳感器誤差。初始條件誤差由主制動(dòng)段以前的任務(wù)決定，傳感器誤差則由導(dǎo)航系統(tǒng)和傳感器本身決定。此外，影響制導(dǎo)精度的因素還包括月球自轉(zhuǎn)、月球不規(guī)則攝動(dòng)等誤差，對(duì)它們的研究可單獨(dú)進(jìn)行，這里暫不做介紹。2.1 誤差模型建立

2.1.1 初始狀態(tài)誤差模型

記著陸器的實(shí)際初始狀態(tài)為Xi，標(biāo)準(zhǔn)初始狀態(tài)為Xn,則定義初始狀態(tài)偏差xi為

xi?Xi?Xn

(7)對(duì)于主制動(dòng)段這一特定的飛行過程，這些偏差都是確定的；而針對(duì)整個(gè)月球探測(cè)任務(wù)，這些偏差就變得具有隨機(jī)性。在本文中，假定xi 的所有元素均服從零均值高斯分布，相互不獨(dú)立，其相關(guān)性取決于前一階段任務(wù)的特性。2.1.2 傳感器誤差模型

由于只研究誤差對(duì)制導(dǎo)律的影響，所以這里假設(shè)需要測(cè)量的量均可由導(dǎo)航系統(tǒng)直接測(cè)得，誤差大小

???????均考慮為典型誤差值。由上一目設(shè)計(jì)的制導(dǎo)律可以看出，需要由導(dǎo)航與控制傳感器測(cè)量的量主要為著陸器相對(duì)于著陸場(chǎng)坐標(biāo)系的位置、速度和加速度。定義待測(cè)量量Q為

?Q??X其估計(jì)值記為Q，則傳感器誤差定義為 ???YZUVWA?

T

q?Q?Q

(8)那么，單個(gè)測(cè)量量的估計(jì)誤差模型可用誤差向量 q的第j(j =1，2?7)個(gè)元素qj 來表示。由參考文獻(xiàn)[5]可知，第 j個(gè)觀測(cè)量的總估計(jì)誤差qj 由以下四部分組成

~?~???-?~qjbsqjn?st???????qt?q?Qt?qt?Qj?t?

(9)jjbcjnc

j100100~~~~~針對(duì)主制動(dòng)這一特定操作階段，上述四部分誤差具有如下特性：

qjbc—第 j 個(gè)觀測(cè)量的測(cè)量誤差，恒為常值，其分布服從零均值高斯分布； qjbs—第 j 個(gè)觀測(cè)量的刻度因素誤差系數(shù)，恒為常值，其分布服從零均值高斯分布； qjnc—第 j 個(gè)觀測(cè)量的隨機(jī)誤差，其為一高斯白噪聲；

qjns

—第 j 個(gè)觀測(cè)量的刻度因素隨機(jī)誤差系數(shù)，其為一高斯白噪聲。

2.2 制導(dǎo)誤差分析

由于采用閉環(huán)制導(dǎo)，制導(dǎo)控制系統(tǒng)對(duì)隨機(jī)誤差具有一定魯棒性，所以本文將著重對(duì)初始偏差和類似于qjbc和qjbs這樣的傳感器常值誤差進(jìn)行仿真研究，分析它們對(duì)制導(dǎo)精度的影響。2.2.1 誤差分析系統(tǒng)建立

誤差分析系統(tǒng)框圖如圖 1 所示，下面將對(duì)其結(jié)構(gòu)進(jìn)行分析。~~~~~~

圖 1 誤差分析系統(tǒng)結(jié)構(gòu)圖

圖中所示初始狀態(tài)偏差實(shí)際上是加在相應(yīng)積分器中。

由前面的分析可知，觀測(cè)量的實(shí)際輸出值受到初始狀態(tài)偏差、傳感器測(cè)量誤差以及傳感器刻度因素誤差的影響，故誤差分析系統(tǒng)模擬程序的實(shí)際輸入應(yīng)包含以下幾部分(以 X通道為例)：

X?X?xi?xbc???~xbsX

(10)100~~

其中，X為觀測(cè)量的實(shí)際輸出值，X 為標(biāo)準(zhǔn)值，xi 為初始狀態(tài)偏差(只在初始時(shí)刻存在)，xbc 為傳感器測(cè)量偏差，xbs為傳感器刻度因素誤差系數(shù)。由圖 1 可以看出，為了更準(zhǔn)確地表示傳感器誤差模型，這里考慮了傳感器的動(dòng)態(tài)性能，其傳遞函數(shù)設(shè)為一階慣性環(huán)節(jié)1?1?Ts?，其中，T 為傳感器時(shí)間常數(shù)，因傳感器的不同而取不同值。

由誤差分析系統(tǒng)結(jié)構(gòu)框圖可以看出，其輸入量主要包括：標(biāo)準(zhǔn)初始狀態(tài)向量、初始狀態(tài)偏差、傳感器測(cè)量誤差、傳感器刻度因素誤差系數(shù)、傳感器時(shí)間常數(shù)、期望終端狀態(tài)；輸出量為加入誤差前后的仿真終端狀態(tài)向量。2.2.2 誤差敏感系數(shù)矩陣求取

在有形如(7)式誤差輸入的情況下，首先根據(jù)圖 1 生成一個(gè)模擬整個(gè)閉環(huán)制導(dǎo)控制系統(tǒng)的數(shù)字仿真程序，然后運(yùn)行該程序，對(duì)比程序輸出即可得到誤差敏感系數(shù)矩陣。具體運(yùn)行過程如下：

第一步：將傳感器誤差設(shè)置為零，初始狀態(tài)設(shè)置為標(biāo)準(zhǔn)值，運(yùn)行模擬程序。這一步稱為標(biāo)準(zhǔn)運(yùn)行。第二步： 將其中一個(gè)傳感器誤差設(shè)置為非零輸入或者設(shè)置一個(gè)非標(biāo)準(zhǔn)初始狀態(tài)，然后進(jìn)行一系列運(yùn)行。

第三步： 將第二步運(yùn)行的系統(tǒng)輸出和標(biāo)準(zhǔn)運(yùn)行的系統(tǒng)輸出進(jìn)行比較即可確定各誤差源的影響。如X 通道標(biāo)準(zhǔn)初始偏差為xi，輸入該誤差前后，X 通道終端狀態(tài)分別為X0 和X1，則 X 通道對(duì)標(biāo)準(zhǔn)初始偏差xi的敏感性可用(X1?X0)/xi來反映。

通過這種方法，可得到一組反映月球軟著陸主制動(dòng)段終端總誤差向量pf和兩個(gè)傳感器誤差向量~??~~qbc、qbs以及初始狀態(tài)偏差向量pi之間關(guān)系的誤差敏感系數(shù)矩陣。由參考文獻(xiàn)[6]可知，其相互關(guān)系可表示為

??~~pf?S1pi?S2qbc?S3qbs（11）

其中，S1、S2和S3分別表示相對(duì)于pi、qbc和qbs的誤差敏感系數(shù)矩陣。

終端誤差向量能用這種形式表示的假設(shè)條件是動(dòng)力學(xué)的線性化必須在標(biāo)準(zhǔn)軌跡區(qū)域內(nèi)。驗(yàn)證該假設(shè)條件的方法有兩種： 擴(kuò)大輸入誤差仿真法和復(fù)合仿真法，這里略去其驗(yàn)證過程。2.2.3 誤差分析

假設(shè)導(dǎo)航系統(tǒng)采用常規(guī)慣性測(cè)量單元，表 1 列出了其典型誤差值，其中，位置誤差能保持在10數(shù)量級(jí)，速度在10數(shù)量級(jí)，加速度為 10g 數(shù)量級(jí)。1-52?~~

運(yùn)用上述方法得到的敏感系數(shù)矩陣給出如下：

?5.502?10-3?-4-3.850?10??1.692?10-3S1??-3?8.362?10?-5.860?10-4?-3?-2.575?10?-2.080?10-4-1.050?10-31.418?10-11.401?10-57.301?10-5-1.001?10-26.411?10-53.240?10-4-4.407?10-2-2.570?10-4-1.862?10-3-5.580?10-11.410?10-57.902?10-51.312?10-55.710?10-4-1.157?10-38.100?10-53.936?10-21.732?10-2-2.743?1017.746?10-1-4.024?10-2-8.939?10-2??3.210?10-34.030?10-3?1.239?10-21.833?10-2?-2-1?8.742?101.414?10?-1.196?10-2-9.901?10-3??-2-2-2.690?10-4.577?10??-6.812?10-1-8.695?10-2-5.203?1002.110?10-14.235?10-16.170?10-3-3.281?1008.202?10-2-5.760?10-35.633?10-1-3.489?102??2.443?101?4.401?102??-9.833?102?6.864?101??23.020?10???-9.859?10-1-1.154?10-3?-40?-3.130?10-1.000?10?-1.379?10-33.560?10-4S2??-2-3?-5.402?101.540?10?1.045?10-31.864?10-3?-34.770?10-4??4.598?109.999?10-13.408?100-7.210?10-43.504?1005.000?10-55.643?10-3-1.527?10-19.368?10-1-6.721?10-1-1.306?10-1?-5.6314?100?-28.479?10??3.730?10-1S3??0?-8.924?10?4.619?10-1?0??2.033?10-5.494?10-1-3.533?10-1-2.810?1001.600?10-31.692?10-16.755?10-18.996?10-1-2095?10-12.473?10-21.664?10-1-1.027?1007.165?10-23.344?100-1.112?1008.613?10-17.852?1003.246?100-1.618?1003.540?10-14.982?10-17.670?10-1-1.122?100-2.397?100-2.380?10-1-3.650?100-2.563?100??2.556?10-1-4.291?10-2?3.401?100-1.888?10-1??-5.103?100-3.230?10-1?3.566?10-12.256?10-1??0-1-7.005?109.930?10??A1、A3:?1??2.759?2,3?0.1297?j2.1329 A2:?1?1.552?2,3??0.6761?j1.8978

由于數(shù)值仿真的起始點(diǎn)選為(1，0，-1)，靠近平衡點(diǎn)(1.5，0，-1.05)，仿真實(shí)驗(yàn)中混沌系統(tǒng)的基頻w0=2.1329，基周期為為T0?2??0?2.9443S。由前面的數(shù)值仿真實(shí)驗(yàn)知要使 Chua’s混沌系統(tǒng)保持其類隨機(jī)性，仿真步長(zhǎng)選在(0.0001，0.7)較為合適，用基周期來表達(dá)即為?129940T015T0? ,15T0?內(nèi)，綜觀三個(gè)連續(xù)混沌系統(tǒng)仿真步長(zhǎng)的理論計(jì)算，我們可以統(tǒng)一選取?15000T0這樣即可以提高仿真運(yùn)算速度，又可以使混沌吸引子的形狀和類隨機(jī)性不發(fā)生變化，這個(gè)選擇范圍也與通常連續(xù)混沌系統(tǒng)數(shù)值仿真步長(zhǎng)的經(jīng)驗(yàn)取值相吻合六、模型結(jié)果及分析

七、結(jié)果分析

八、模型評(píng)價(jià)與改進(jìn)方向

九、參考文獻(xiàn)

第四篇：全國(guó)大學(xué)生數(shù)學(xué)建模感想

數(shù)學(xué)建模有感

這一暑假，參加學(xué)校組織的學(xué)習(xí)了 “數(shù)學(xué)建?！?。這門課程，數(shù)學(xué)建模這門課程與其它的數(shù)學(xué)課程有一個(gè)很大的區(qū)別，數(shù)學(xué)建模是把現(xiàn)實(shí)世界中的實(shí)際問題加以提煉，抽象為數(shù)學(xué)模型，求出模型的解，驗(yàn)證模型的合理性．并用該數(shù)學(xué)模型所提供的解答來解釋現(xiàn)實(shí)問題，它是一門將數(shù)學(xué)綜合應(yīng)用到實(shí)際中解決實(shí)際問題的學(xué)科。建立模型就是對(duì)于現(xiàn)實(shí)中的原型．為了某個(gè)特定目的，作出一些必要的簡(jiǎn)化和假設(shè)。運(yùn)后適當(dāng)?shù)臄?shù)學(xué)工具得到一個(gè)數(shù)學(xué)結(jié)構(gòu)、它或者能解釋特定現(xiàn)象的現(xiàn)實(shí)狀態(tài)，或者能預(yù)測(cè)到對(duì)象的未來狀況，或者能提供處理對(duì)象的最優(yōu)決策或控制，而現(xiàn)實(shí)生活中很多的問題都可以通過建立數(shù)學(xué)模型來解決，所以學(xué)習(xí)數(shù)學(xué)建模，徹底讓我體會(huì)到了數(shù)學(xué)在生活中的重要之處。

一次建模，終生受益，這話一點(diǎn)也不假。在沒有接觸過數(shù)學(xué)建模這門學(xué)科之前，我一直認(rèn)為數(shù)學(xué)是一門純理論的學(xué)科，但是數(shù)學(xué)建模卻能把它應(yīng)用到實(shí)際中去，并用它去解決很多來自日常生活及經(jīng)濟(jì)、工程、理、化、生、醫(yī)等學(xué)科中的問題。

我覺得數(shù)學(xué)建模能為學(xué)生提供自主學(xué)習(xí)、自主探索、自主提出問題、自主解決問題的機(jī)會(huì)，培養(yǎng)學(xué)生的數(shù)學(xué)觀念、科學(xué)態(tài)度和合作精神，激發(fā)學(xué)生的學(xué)習(xí)興趣，培養(yǎng)學(xué)生認(rèn)真求實(shí)、崇尚真理、追求完美、講究效益、聯(lián)系實(shí)際的學(xué)習(xí)態(tài)度和學(xué)習(xí)習(xí)慣。它能提高學(xué)生應(yīng)用所學(xué)的數(shù)學(xué)知識(shí)解決實(shí)際問題的能力，從過去強(qiáng)調(diào)數(shù)學(xué)知識(shí)的“有用、可用”，到使學(xué)生所學(xué)知識(shí)的“想用、能用和會(huì)用”，讓學(xué)生更多自主的實(shí)踐，把學(xué)習(xí)知識(shí)、應(yīng)用知識(shí)、探索發(fā)現(xiàn)、使用計(jì)算機(jī)、培養(yǎng)良好的科學(xué)態(tài)度與思維品質(zhì)更好地結(jié)合起來，使學(xué)生在問題解決的過程中得到學(xué)數(shù)學(xué)、用數(shù)學(xué)的實(shí)際體驗(yàn)，加深對(duì)數(shù)學(xué)的理解。綜合上述可見，開展數(shù)學(xué)建?；顒?dòng)是非常有必要的。應(yīng)該在學(xué)校大力推廣，讓更多同學(xué)在參與中受益。

通過這次數(shù)學(xué)建模學(xué)到了很多東西，首先，是從現(xiàn)實(shí)生活中發(fā)現(xiàn)問題，這就需要我們用心觀察；然后就是解決問題的方法，由于我們不可能在課堂上學(xué)到所有知識(shí)，很多東西是要我們自學(xué)的，這就培養(yǎng)了我們自學(xué)能力，還有自己解決問題的能力：最后，因?yàn)槟Ｐ鸵⒃谡鎸?shí)數(shù)據(jù)上，就要求我們要有實(shí)事求是的態(tài)度了。數(shù)學(xué)建模雖然就告一段落了，學(xué)到的方法知識(shí)卻伴隨著我們以后的學(xué)習(xí)工作?！皵?shù)學(xué)建?！闭n已使“數(shù)學(xué)建?！钡暮诵乃枷肷钌钤谖乙庾R(shí)中扎根，使我在今后學(xué)習(xí)中越來越善于發(fā)現(xiàn)問題并用數(shù)學(xué)知識(shí)創(chuàng)造性的去分析解決問題。數(shù)學(xué)建模課對(duì)學(xué)生思維能力的訓(xùn)練和思維方法上的引導(dǎo)，這就是它的主要魅力所在。一年一度的全國(guó)數(shù)學(xué)建模大賽在在9 月13 日上午8 點(diǎn)拉開戰(zhàn)幕，各隊(duì) 將在3 天72 小時(shí)內(nèi)對(duì)一個(gè)現(xiàn)實(shí)中的實(shí)際問題進(jìn)行模型建立，求解和分析，確定

題目后，我們隊(duì)三人分頭行動(dòng)，一人去圖書館查閱資料，一人在網(wǎng)上搜索相關(guān)信 息，一人建立模型，通過三人的努力，在前兩天中建立出兩個(gè)模型并編程求解，經(jīng)過艱苦的奮斗，終于在第三天完成了論文的寫作，在這三天里我感觸很深，現(xiàn) 將心得體會(huì)寫出，希望與大家交流。1.團(tuán)隊(duì)精神：

團(tuán)隊(duì)精神是數(shù)學(xué)建模是否取得好成績(jī)的最重要的因素，一隊(duì)三個(gè)人要相互支 持，相互鼓勵(lì)。切勿自己只管自己的一部分（數(shù)學(xué)好的只管建模，計(jì)算機(jī)好的只 管編程，寫作好的只管論文寫作），很多時(shí)候，一個(gè)人的思考是不全面的，只有

大家一起討論才有可能把問題搞清楚，因此無論做任何板塊，三個(gè)人要一起齊心 才行，只靠一個(gè)人的力量，要在三天之內(nèi)寫出一篇高水平的文章幾乎是不可能的。2.有影響力的leader：

在比賽中，leader 是很重要的，他的作用就相當(dāng)與計(jì)算機(jī)中的CPU，是全隊(duì) 的核心，如果一個(gè)隊(duì)的leader 不得力，往往影響一個(gè)隊(duì)的正常發(fā)揮，就拿選題來說，有人想做A 題，有人想做B 題，如果爭(zhēng)論一天都未確定方案的話，可能就沒有足夠時(shí)間完成一篇論文了，又比如，當(dāng)隊(duì)中有人信心動(dòng)搖時(shí)（特別是第三天，人可能已經(jīng)心力交瘁了），leader 應(yīng)發(fā)揮其作用，讓整個(gè)隊(duì)伍重整信心，否則可能導(dǎo)致隊(duì)伍的前功盡棄。3.合理的時(shí)間安排：

做任何事情，合理的時(shí)間安排非常重要，建模也是一樣，事先要做好一個(gè)規(guī) 劃，建模一共分十個(gè)板塊（摘要，問題提出，模型假設(shè)，問題分析，模型假設(shè)，模型建立，模型求解，結(jié)果分析，模型的評(píng)價(jià)與推廣，參考文獻(xiàn)，附錄）。你每 天要做完哪幾個(gè)板塊事先要確定好，這樣做才會(huì)使自己游刃有余，保證在規(guī)定時(shí) 間內(nèi)完成論文，以避免由于時(shí)間上的不妥，以致于最后無法完成論文。4.正確的論文格式：

論文屬于科學(xué)性的文章，它有嚴(yán)格的書寫格式規(guī)范，因此一篇好的論文一定 要有正確的格式，就拿摘要來說吧，它要包括6 要素（問題,方法,模型,算法,結(jié)論,特色），它是一篇論文的概括，摘要的好壞將決定你的論文是否吸引評(píng)委的目光，但聽閱卷老師 說，這次有些論文的摘要里出現(xiàn)了大量的圖表和程序，這都是不符合論文格式的，這種論

文也不會(huì)取得好成績(jī)，因此我們寫論文時(shí)要端正態(tài)度，注意書寫格式。5.論文的寫作：

我個(gè)人認(rèn)為論文的寫作是至關(guān)重要的，其實(shí)大家最后的模型和結(jié)果都差不

多，為什么有些隊(duì)可以送全國(guó)，有些隊(duì)可以拿省獎(jiǎng)，而有些隊(duì)卻什么都拿不到，這關(guān)鍵在于論文的寫作上面。一篇好的論文首先讀上去便使人感到邏輯清晰，有 條例性，能打動(dòng)評(píng)委；其次，論文在語言上的表述也很重要，要注意用詞的準(zhǔn)確 性；另外，一篇好的論文應(yīng)有閃光點(diǎn)，有自己的特色，有自己的想法和思考在里 面，總之，論文寫作的好壞將直接影響到成績(jī)的優(yōu)劣。

6.算法的設(shè)計(jì)：算法的設(shè)計(jì)的好壞將直接影響運(yùn)算速度的快慢，建議大家多用數(shù)學(xué)軟件（Mathematice,Matlab,Maple, Mathcad,Lindo,Lingo,SAS 等），這里提供十種數(shù)學(xué)以上便是我這次參加這次數(shù)學(xué)建模競(jìng)賽的一點(diǎn)心得體會(huì)，只當(dāng)貽笑大方，不過就數(shù)學(xué)建模本身而言，它是魅力無窮的，它能夠鍛煉和考查一個(gè)人的綜合素質(zhì)，也希望廣大同學(xué)能夠積極參與到這項(xiàng)活動(dòng)當(dāng)中來。

第五篇：全國(guó)大學(xué)生數(shù)學(xué)建模競(jìng)賽

全國(guó)大學(xué)生數(shù)學(xué)建模競(jìng)賽

1、數(shù)模競(jìng)賽的起源與歷史

數(shù)模競(jìng)賽是由美國(guó)工業(yè)與應(yīng)用數(shù)學(xué)學(xué)會(huì)在1985年發(fā)起的一項(xiàng)大學(xué)生競(jìng)賽活動(dòng)，目的在于激勵(lì)學(xué)生學(xué)習(xí)數(shù)學(xué)的積極性，提高學(xué)生建立數(shù)學(xué)模型和運(yùn)用計(jì)算機(jī)技術(shù)解決實(shí)際問題的綜合能力，鼓勵(lì)廣大學(xué)生踴躍參加課外科技活動(dòng)，開拓知識(shí)面，培養(yǎng)創(chuàng)精神及合作意識(shí)，推動(dòng)大學(xué)數(shù)學(xué)教學(xué)體系、教學(xué)內(nèi)容和方法的改革。我國(guó)大學(xué)生數(shù)學(xué)建模競(jìng)賽是由教育部高教司和中國(guó)工業(yè)與數(shù)學(xué)學(xué)會(huì)主辦、面向全國(guó)高等院校的、每年一屆的通訊競(jìng)賽。其宗旨是：創(chuàng)新意 識(shí)、團(tuán)隊(duì)精神、重在參與、公平競(jìng)爭(zhēng)。1992載在中國(guó)創(chuàng)辦，自從創(chuàng)辦以來，得到了教育部高教司和中國(guó)工業(yè)與應(yīng)用數(shù)學(xué)協(xié)會(huì)的得力支持和關(guān)心，呈現(xiàn)出迅速的發(fā)展發(fā)展勢(shì)頭，就2003年來說，報(bào)名階段須然受到“非典”影響，但是全國(guó)30個(gè)?。ㄊ小⒆灾螀^(qū)）及香港的637所院校就有5406隊(duì)參賽，在職業(yè)技術(shù)學(xué)院增加更快，參賽高校由2002年的1067所上升到了2003年的1410所?？梢哉f：數(shù)學(xué)建模已經(jīng)成為全國(guó)高校規(guī)模最大課外科技活動(dòng)。

2、什么是數(shù)學(xué)建模

數(shù)學(xué)建模（Mathematical Modelling）是一種數(shù)學(xué)的思考方法，是“對(duì)現(xiàn)實(shí)的現(xiàn)象通過心智活動(dòng)構(gòu)造出能抓住其重要且有

用的特征的表示，常常是形象化的或符號(hào)的表示?！睆目茖W(xué)，工程，經(jīng)濟(jì)，管理等角度看數(shù)學(xué)建模就是用數(shù)學(xué)的語言和方法，通過抽象，簡(jiǎn)化建立能近似刻畫并“解決”實(shí)際問題的一種強(qiáng)有力的數(shù)學(xué)工具。顧名思義，modelling一詞在英文中有“塑造藝術(shù)”的意思，從而可以理解從不同的側(cè)面，角度去考察問題就會(huì)有不盡的數(shù)學(xué)模型，從而數(shù)學(xué)建模 的創(chuàng)造又帶有一定的藝術(shù)的特點(diǎn)。而數(shù)學(xué)建模最重要的特點(diǎn)是要接受實(shí)踐的檢驗(yàn)，多次修改模型漸趨完善的過程。

3、競(jìng)賽的內(nèi)容：

競(jìng)賽題目一般來源于工程技術(shù)和管理科學(xué)等方面經(jīng)過適當(dāng)簡(jiǎn)化加工的實(shí)際問題，不要求參賽者預(yù)先掌握深入的專門知識(shí)，只需要學(xué)過普通高校的數(shù)學(xué)課程。題目有較大的靈活性供參賽者發(fā)揮其創(chuàng)造能力。參賽者應(yīng)根據(jù)題目要求，完成一篇包括模型假設(shè)、建立和求解、計(jì)算方法的設(shè)計(jì)和計(jì)算機(jī)實(shí)現(xiàn)、結(jié)果的分析和檢驗(yàn)、模型的改進(jìn)等方面的論文（即答卷）。競(jìng)賽評(píng)獎(jiǎng)以假設(shè)的合理性、建模的創(chuàng)造性、結(jié)果的正確性和文字表述的清晰程度為主要標(biāo)準(zhǔn)。

4、競(jìng)賽的步驟

建模是一種十分復(fù)雜的創(chuàng)造性勞動(dòng)，現(xiàn)實(shí)世界中的事物形形色色，五花八門，不可能用一些條條框 框規(guī)定出各種模型如何具體建立，這里只是大致歸納一下建模的一般步驟和原則：

1）模型準(zhǔn)備：首先要了解問題的實(shí)際背景，明確題目的要求，收集各種必要的信息．

2）模型假設(shè)：為了利用數(shù)學(xué)方法，通常要對(duì)問題做必要的、合理的假設(shè)，使問題的主要特征凸現(xiàn)出來，忽略問題的次要方面。

3）模型構(gòu)成：根據(jù)所做的假設(shè)以及事物之間的聯(lián)系，構(gòu)造各種量之間的關(guān)系,把它問題化

4）模型求解：利用已知的數(shù)學(xué)方法來求解上一步所得到的數(shù)學(xué)問題，此時(shí)往往還要作出進(jìn)一步的簡(jiǎn)化或假設(shè)。注意要盡量采用簡(jiǎn)單的數(shù)學(xué)工具。

5）模型分析：對(duì)所得到的解答進(jìn)行分析，特別要注意當(dāng)數(shù)據(jù)變化時(shí)所得結(jié)果是否穩(wěn)定。

6）模型檢驗(yàn)：分析所得結(jié)果的實(shí)際意義，與實(shí)際情況進(jìn)行比較，看是否符合實(shí)際，如果不夠理想，應(yīng)該修改、補(bǔ)充假設(shè)，或重新建模，不斷完善。

7）模型應(yīng)用：所建立的模型必須在實(shí)際應(yīng)用中才能產(chǎn)生效益，在應(yīng)用中不斷改進(jìn)和完善。

5、模型的分類

按模型的應(yīng)用領(lǐng)域分類： 生物數(shù)學(xué)模型、醫(yī)學(xué)數(shù)學(xué)模型、地質(zhì)數(shù)學(xué)模型、數(shù)量經(jīng)濟(jì)學(xué)模型、數(shù)學(xué)社會(huì)學(xué)模型

按是否考慮隨機(jī)因素分類 ：確定性模型、隨機(jī)性模型按是否考慮模型的變化分類 ：靜態(tài)模型、動(dòng)態(tài)模型按應(yīng)用離散方法或連續(xù)方法 ：離散模型、連續(xù)模型

按建立模型的數(shù)學(xué)方法分類 ：幾何模型、微分方程模型、圖

論模型、規(guī)劃論模型、馬氏鏈模型 按人們對(duì)事物發(fā)展過程的了解程度分類 ：

白箱模型： 指那些內(nèi)部規(guī)律比較清楚的模型。如力學(xué)、熱學(xué)、電學(xué)以及相關(guān)的工程技術(shù)問題。

灰箱模型： 指那些內(nèi)部規(guī)律尚不十分清楚，在建立和改善模型方面都還不同程度地有許多工作要做的問題。如氣象學(xué)、生態(tài)學(xué)經(jīng)濟(jì)學(xué)等領(lǐng)域的模型。

黑箱模型：

指一些其內(nèi)部規(guī)律還很少為人們所知的現(xiàn)象。如生命科學(xué)、社會(huì)科學(xué)等方面的問題。但由于因素眾多、關(guān)系復(fù)雜，也可簡(jiǎn)化為灰箱模型來研究。

6、數(shù)學(xué)建模應(yīng)用

今天，在國(guó)民經(jīng)濟(jì)和社會(huì)活動(dòng)的以下諸多方面，數(shù)學(xué)建模都有著非常具體的應(yīng)用。

1分析與設(shè)計(jì)： 例如描述藥物濃度在人體內(nèi)的變化規(guī)律以分析藥物的療效；建立跨音速空氣流和激波的數(shù)學(xué)模型，用數(shù)值模擬設(shè)計(jì)新的飛機(jī)翼型。預(yù)報(bào)與決策： 生產(chǎn)過程中產(chǎn)品質(zhì)量指標(biāo)的預(yù)報(bào)、氣象預(yù)報(bào)、人口預(yù)報(bào)、經(jīng)濟(jì)增長(zhǎng)預(yù)報(bào)等等，都要有預(yù)報(bào)模型。使經(jīng)濟(jì)效益最大的價(jià)格策略、使費(fèi)用最少的設(shè)備維修方案，是決策模型的例子。3 控制與優(yōu)化： 電力、化工生產(chǎn)過程的最優(yōu)控制、零件設(shè)計(jì)中的參數(shù)優(yōu)化，要以數(shù)學(xué)模型為前提。建立大系統(tǒng)控制與優(yōu)化的數(shù)

學(xué)模型，是迫切需要和十分棘手的課題。規(guī)劃與管理 生產(chǎn)計(jì)劃、資源配置、運(yùn)輸網(wǎng)絡(luò)規(guī)劃、水庫(kù)優(yōu)化調(diào)度，以及排隊(duì)策略、物資管理等，都可以用運(yùn)籌學(xué)模型解決 報(bào)名時(shí)間：從大賽的通知文稿發(fā)出后，就可以報(bào)名了，報(bào)名截止時(shí)間一般在開始比賽的前7到10天。

競(jìng)賽時(shí)間：每年的9月的第三個(gè)星期五上午8時(shí)至下一個(gè)星期一上午8時(shí)。

報(bào)名方式：如果有分賽區(qū)（每個(gè)賽區(qū)應(yīng)至少有6所院校的20個(gè)隊(duì)參加），就聯(lián)系分賽區(qū)報(bào)名，沒有分賽區(qū)，則直接向主委會(huì)報(bào)名。

大學(xué)生以隊(duì)為單位參賽，每隊(duì)3人（須屬于同一所學(xué)校），專業(yè)不限。競(jìng)賽分本科、?？苾山M進(jìn)行，本科生參加本科組競(jìng)賽，專科生參加?？平M競(jìng)賽（也可參加本科組競(jìng)賽），研究生不得參加。每隊(duì)可設(shè)一名指導(dǎo)教師（或教師組）。