【正文】 OK_Callback(hObject, eventdata, handles)hObject OK控件的句柄(詳情請參見 GCBO)eventdata 保留在MATLAB里面的參數(shù)，它的作用是方便后面的程序使用handles GUI界面的句柄，可以通過它獲取整個界面的信息（詳情請參見GUIDATA）清除坐標系原有的內(nèi)容cla獲得程序使用者輸入的小球質量、剛開始小球的位置與速度，設置萬有引力常量默認值以及程序的運動時間g = *10^11。white39。defaultUicontrolBackgroundColor39。amp。)。white39。defaultUicontrolBackgroundColor39。amp。)。white39。defaultUicontrolBackgroundColor39。amp。)。white39。defaultUicontrolBackgroundColor39。amp。)。white39。defaultUicontrolBackgroundColor39。amp。)。white39。defaultUicontrolBackgroundColor39。amp。)。white39。defaultUicontrolBackgroundColor39。amp。)。white39。defaultUicontrolBackgroundColor39。amp。)。white39。defaultUicontrolBackgroundColor39。amp。)。white39。defaultUicontrolBackgroundColor39。amp。)。white39。defaultUicontrolBackgroundColor39。amp。)。white39。defaultUicontrolBackgroundColor39。amp。)。white39。defaultUicontrolBackgroundColor39。amp。)。white39。defaultUicontrolBackgroundColor39。amp。)。white39。defaultUicontrolBackgroundColor39。amp。)。white39。defaultUicontrolBackgroundColor39。amp。)。white39。defaultUicontrolBackgroundColor39。amp。)。white39。defaultUicontrolBackgroundColor39。amp。)。white39。defaultUicontrolBackgroundColor39。amp。)。white39。defaultUicontrolBackgroundColor39。amp。)。white39。defaultUicontrolBackgroundColor39。amp。)。此前指令更改了handles，必須靠guidata指令才能把更新了的handles加以保存，以供后續(xù)使用guidata(hObject,handles)。 = 0。,0)。String39。set(,39。,0)。String39。set(,39。 = 0。將地球的參數(shù)保存為“GUI數(shù)據(jù)”形式，以便共享 = *10^23。,)。String39。set(,39。,)。String39。 = 0。 = 0。set(,39。,0)。String39。set(,39。,)。此前指令更改了handles，必須靠guidata指令才能把更新了的handles加以保存，以供后續(xù)使用guidata(hObject, handles)。 = 0。初始化第三個天體的質量，位置和速度 = 0。 = 0。 = 0。 = 0。（2）三體問題的編碼1o三個天體各個參數(shù)的初始化部分。圖31 三體問題受力分析 運動規(guī)律及迭代方程一、 萬有引力定律f12=gm1m2x1x22+y1y22+z1z22f23=gm2m3x2x32+y2y32+z2z32f31=gm3m1x3x12+y3y12+z3z12二、 牛頓第二定律a1=f12+f31m1a2=f12f23m2a3=f23f31m3三、 速度迭代p1=p1+a1tp2=p2+a2tp3=p3+a3t四、 位置迭代d1=d1+p1t+12a1t2d2=d2+p2t+12a2t2d3=d3+p3t+12a3t2 求解時還須將加速度a，速度迭代方程，位置迭代方程按直角坐標系分解如附錄1。 受力分析如圖31所示，每個天體都受到其他兩個天體的萬有引力作用。2 模型假設1) 忽略相對論效應的影響；2) 不考慮其它天體對系統(tǒng)的作用，所處宇宙空間是各向同性的；3) 忽略宇宙空間粒子、潮汐力的影響；4) 假設運動過程中沒有能量損失；5) 將天體視為質點，不考慮它們在運動的過程中彼此之間的相互碰撞；3 模型建立 具體時刻的狀態(tài)因為三個天體的位置和速度是隨時間變化的，現(xiàn)在取很短的時間間隔t，假定這段時間間隔內(nèi)系統(tǒng)參數(shù)沒有變化，則可以考慮從系統(tǒng)每一時刻的狀態(tài)計算后下一時刻的狀態(tài)。同時，本文利用MATLAB的GUI功能，設計出簡潔友好的界面，并且通過動畫的方法將天體的位置和速度在空間直角坐標系上實時的表示出來，作出了天體物理實驗用計算機模擬來進行的嘗試。因此，除了解決一個重要的難題，更重要的是在解決問題的過程當中產(chǎn)生了好幾種嶄新的數(shù)學思想。截止到目前為止，三體問題特解的族數(shù)被擴充到了16族。這三類方法各有利弊。天體物理中的三體問題可以簡單的被描述為：三個天體在他們之間的相互之間的萬有引力作用下，給出各個天體的初始條件，包括三個天體的質量、初始位置和初始速度，確定它們的后續(xù)運動，包括位置和速度隨時間的變化。在1744年，這個問題被Euler完全的解決了。在1710年，Johann這部分的天體力學以牛頓的萬有引力理論開始。因此，可以用現(xiàn)代動力天文學的基礎理論來研究這個問題。 牛頓曾經(jīng)用幾何學原理解決了只存在相互引力作用的兩個球的運動這一兩體問題。Bernoulli因為他關于分析處理兩體問題的研究獲得了法國科學院獎。在經(jīng)典力學，量子力學和天文學中，三體問題是一個出名的，還沒有解決的問題，它被Whittaker描述為在所有動力學問題中最出名的問題，Hilbert認為它是一個完美數(shù)學問題的標準。探究三體問題的方式大抵能夠分為三類：第一類是分析方法，其基本原理是把天體的坐標和速度展開為時間或其它小參數(shù)的級數(shù)形式的近似解析表達式，討論天體的坐標和軌道要素隨時間的變化規(guī)律；第二類是定性方法，選用微分方程的定性理論來研究很長一段時間內(nèi)三體運動的宏觀運動規(guī)律和全局性質；第三類是數(shù)值方法，這是直接依據(jù)微分方程的計算方式得出天體在特定時刻的詳細位置和速度。在十八和十九世紀，幾乎所有著名的數(shù)學家們都曾試圖解決這個問題，但是問題的進展不大。多年來，人們對這個問題進行了大量的研究，從1750年到20世紀初，關于這個問題有800多篇援引了許多杰出的數(shù)學家和天文學家成果的論文發(fā)表。 本文的研究內(nèi)容本文決定通過計算機的數(shù)值運算，利用經(jīng)典力學中的牛頓三大定律和萬有引力定律迭代求解天體運動中的三體問題。同時，還可以在不同的觀察角度了解它們的運動過程。因此，我們可以將時間分成許多小部分，用計算機迭代求解，得到隨時間變化的系統(tǒng)的各個參量和系統(tǒng)的末狀態(tài)。在萬有引力f1f2f31作用下，天體天體天體3的加速度分別為f12+f31m1，f12f23m2，f23f31m3。本程序的編碼和界面設計部分用MATLAB進行開發(fā)。 = 0。 = 0。 = 0。 = 0。 = 0。 = hObject。String39。,0)。set(,39。String39。,0)。將太陽的參數(shù)保存為“GUI數(shù)據(jù)”形式，以便共享 = *10^29。 = 0。set(,39。String39。,0)。set(,39。String39。,0)。 = 0。在圖形界面的第三個天體的文本輸入框內(nèi)顯示默認的輸入值，此輸入值為月球的參數(shù)。set(,39。String39。,0)。set(,39。String39。 = 0。 = 0。String39。設置完所有屬性后，在對象的創(chuàng)建過程中執(zhí)行以下代碼function FirstBodyX_CreateF(hObject, eventdata, handles)hObject FirstBodyX 控件的句柄(詳情請參見 GCBO)eventdata 保留在MATLAB里面的參數(shù)，它的作用是方便后面的程序使用handles GUI界面的句柄，它會在調(diào)用完所有的CreateF函數(shù)后創(chuàng)建提示: 可編輯文本框控件通常在窗口上有一個白色的背景，更多詳情請參見ISPC 和COMPUTER.if ispc amp。), get(0,39。,39。String39。設置完所有屬性后，在對象的創(chuàng)建過程中執(zhí)行以下代碼function FirstBodyY_CreateF(hObject, eventdata, handles)hObject FirstBodyY 控件的句柄(詳情請參見 GCBO)eventdata 保留在MATLAB里面的參數(shù)，它的作用是方便后面的程序使用handles GUI界面的句柄，它會在調(diào)用完所有的CreateF函數(shù)后創(chuàng)建提示: 可編輯文本框控件通常在窗口上有一個白色的背景，更多詳情請參見ISPC和COMPUTER.if ispc amp。), get(0,39。,39。String39。設置完所有屬性后，在對象的創(chuàng)建過程中執(zhí)行以下代碼function FirstBodyZ_CreateF(hObject, eventdata, handles)hObject FirstBodyZ 控件的句柄(詳情請參見 GCBO)eventdata 保留在MATLAB里面的參數(shù)，它的作用是方便后面的程序使用handles GUI界面的句柄，它會在調(diào)用完所有的CreateF函數(shù)后創(chuàng)建提示: 可編輯文本框控件通常在窗口上有一個白色的背景，更多詳情請參見ISPC 和COMPUTER.if ispc amp。), get(0,39。,39。String39。設置完所有屬性后，在對象的創(chuàng)建過程中執(zhí)行以下代碼function FirstBodyU_CreateF(hObject, eventdata, handles)設置完所有屬性后，在對象的創(chuàng)建過程中執(zhí)行以下代碼function FirstBodyX_CreateF(hObject, eventdata, handles)hObject FirstBodyU 控件的句柄(詳情請參見 GCBO)eventdata 保留在MATLAB里面的參數(shù)，它的作用是方便后面的程序使用handles GUI界面的句柄，它會在調(diào)用完所有的CreateF函數(shù)后創(chuàng)建提示: 可編輯文本框控件通常在窗口上有一個白色的背景，更多詳情請參見ISPC 和COMPUTER.if ispc amp。), get(0,39。,39。String39。設置完所有屬性后，在對象的創(chuàng)建過程中執(zhí)行以下代碼function FirstBodyV_CreateF(hObject, eventdata, handles)hObject FirstBodyV 控件的句柄(詳情請參見 GCBO)eventdata 保留在MATLAB里面的參數(shù)，它的作用是方便后面的程序使用handles GUI界面的句柄，它會在調(diào)用完所有的CreateF函數(shù)后創(chuàng)建提示: 可編輯文本框控件通常在窗口上有一個白色的背景，更多詳情請參見ISPC 和COMPUTER.if ispc amp。), get(0,39。,39。String39。設置完所有屬性后，在對象的創(chuàng)建過程中執(zhí)行以下代碼function FirstBodyW_CreateF(hObject, eventdata, handles)hObject FirstBodyW 控件的句柄(詳情請參見 GCBO)eventdata 保留在MATLAB里面的參數(shù)，它的作用是方便后面的程序使用handles GUI界面的句柄，它會在調(diào)用完所有的CreateF函數(shù)后創(chuàng)建提示: 可編輯文本框控件通常在窗口上有一個白色的背景，更多詳情請參見ISPC 和COMPUTER.if ispc amp。), get(0,39。,39。String39。設置完所有屬性后，在對象的創(chuàng)建過程中執(zhí)行以下代碼function FirstBodyM_CreateF(hObject, eventdata, handles)hObject FirstBodyM 控件的句柄(詳情請參見 GCBO)eventdata 保留在MATLAB里面的參數(shù)，它的作用是方便后面的程序使用handles GUI界面的句柄，它會在調(diào)用完所有的CreateF函數(shù)后創(chuàng)建提示: 可編輯文本框控件通常在窗口上有一個白色的背景，更多詳情請參見ISPC和 COMPUTER.if ispc amp。), get(0,39。,39。String39。設置完所有屬性后，在對象的創(chuàng)建過程中執(zhí)行以下代碼function SecondBodyX_CreateF(hObject, eventdata, handles)hObject SecondBodyX控件的句柄(詳情請參見 GCBO)eventdata 保留在MATLAB里面的參數(shù)，它的作用是方便后面的程序使用handles GUI界面的句柄，它會在調(diào)用完所有的CreateF函數(shù)后創(chuàng)建提示: 可編輯文本框控件通常在窗口上有一個白色的背景，更多詳情請參見ISPC和COMPUTER.if ispc amp。), get(0,39。,39。String39。設置完所有屬性后，在對象的創(chuàng)建過程中執(zhí)行以下代碼function SecondBodyY_CreateF(hObjec