【正文】 LAB里面的參數(shù)，它的作用是方便后面的程序使用handles GUI界面的句柄，可以通過它獲取整個界面的信息（詳情請參見GUIDATA）獲得當(dāng)前輸入框的輸入字符tempData=get(hObject,39。BackgroundColor39。endfunction ThirdBodyM_Callback(hObject, eventdata, handles)hObject ThirdBodyM控件的句柄(詳情請參見 GCBO)eventdata 保留在MATLAB里面的參數(shù)，它的作用是方便后面的程序使用handles GUI界面的句柄，可以通過它獲取整個界面的信息（詳情請參見GUIDATA）獲得當(dāng)前輸入框的輸入字符tempData=get(hObject,39。BackgroundColor39。end3o開始計算三個天體的實時位置和速度使用三個天體的各個參數(shù)進(jìn)行迭代運算，同時畫出三個天體的運行動畫。v1 = 。w2 = 。控制迭代時間time = 100000。,[1 0 0],39。,39。.39。 ，用k來控制迭代進(jìn)行的次數(shù)for k=0::time設(shè)置步長t=100000。 以下為對y軸方向的計算ay1=f12(2)/m1。 z1=z1+w1*t+1/2*az1*t^2。,y1/10^10,39。zdata39。,[0 0 1],39。,39。.39。，用k來控制迭代進(jìn)行的次數(shù)for k=0::time設(shè)置步長t=100000。 以下為對y軸方向的計算ay2=f23(2)/m2。 z2=z2+w2*t+1/2*az2*t^2。,y2/10^10,39。zdata39。,[1 0 0],39。,39。.39。，用k來控制迭代進(jìn)行的次數(shù)for k=0::time設(shè)置步長t=100000。 以下為對y軸方向的計算ay1=+f31(2)/m1。 z1=z1+w1*t+1/2*az1*t^2。,y1/10^10,39。zdata39。,39。Color39。erasemode39。,39。)。坐標(biāo)值迭代x1=x1+u1*t+1/2*ax1*t^2。ay3=(f23(2)f31(2))/m3。az2=(f12(3)f23(3))/m2。 重置三個小球的位置set(h,39。xdata39。,x3/10^10,39。Value39。改變觀察視角view([az,el])。end用于控制方位角的文本輸入框的回調(diào)函數(shù)，用于同步方位角滑鍵的值，同時改變觀察者的方位角：function azEditText_Callback(hObject, eventdata, handles)hObject azEditText控件的句柄(詳情請參見 GCBO)eventdata 保留在MATLAB里面的參數(shù)，它的作用是方便后面的程序使用handles GUI界面的句柄，可以通過它獲取整個界面的信息（詳情請參見GUIDATA）獲得當(dāng)前輸入框的輸入字符。String39。defaultUicontrolBackgroundColor39。)。設(shè)置完所有屬性后，在對象的創(chuàng)建過程中執(zhí)行以下代碼 function elSlider_CreateF(hObject, eventdata, handles)hObject elSlider控件的句柄(詳情請參見 GCBO)eventdata 保留在MATLAB里面的參數(shù)，它的作用是方便后面的程序使用handles GUI界面的句柄，它會在調(diào)用完所有的CreateF函數(shù)后創(chuàng)建提示: 滑塊控件通常有一個淺灰色的背景if isequal(get(hObject,39。同時，將字符轉(zhuǎn)換為雙精度數(shù)el=str2double(get(hObject,39。))。)) set(hObject,39。,39。設(shè)置完所有屬性后，在對象的創(chuàng)建過程中執(zhí)行以下代碼function elEditText_CreateF(hObject, eventdata, handles)hObject elEditText控件的句柄(詳情請參見 GCBO)eventdata 保留在MATLAB里面的參數(shù)，它的作用是方便后面的程序使用handles GUI界面的句柄，它會在調(diào)用完所有的CreateF函數(shù)后創(chuàng)建提示: 可編輯文本框控件通常在窗口上有一個白色的背景，更多詳情請參見ISPC 和 COMPUTER.if ispc amp。))。), get(0,39。String39。BackgroundColor39。改變觀察視角view([az,el])。String39。BackgroundColor39。將字符轉(zhuǎn)換為雙精度數(shù)el=str2double(get(,39。,y3/10^10,39。ydata39。,x1/10^10,39。 z1=z1+w1*t+1/2*az1*t^2。y2=y2+v2*t+1/2*ay2*t^2。x3=x3+u3*t+1/2*ax3*t^2。萬有引力定律 f12=g*m1*m2/(sqrt((x1x2)^2+(y1y2)^2+(z1z2)^2))^3*[x1x2,y1y2,z1z2]。,39。normal39。Marker39。,50,39。 刷新畫面drawnow 控制速度pause()endelse 設(shè)置天體顏色、點型、大小等參數(shù)，每次循環(huán)都重新設(shè)置整個畫面h=line(39。,z1/10^10)。 w1=w1+az1*t。y1=y1+v1*t+1/2*ay1*t^2。 以下為對x軸方向的計算牛頓第二定律ax1=+f31(1)/m1。MarkerSize39。)。,39。 刷新畫面drawnow 控制速度pause()endelseif m2 == 0 amp。,z2/10^10)。w2=w2+az2*t。y2=y2+v2*t+1/2*ay2*t^2。 以下為對x軸方向的計算牛頓第二定律ax2=f23(1)/m2。MarkerSize39。)。,39。 刷新畫面drawnow 控制速度pause()endelseif m1 == 0 amp。,z1/10^10)。 w1=w1+az1*t。y1=y1+v1*t+1/2*ay1*t^2。 以下為對x軸方向的計算牛頓第二定律ax1=f12(1)/m1。MarkerSize39。)。,39。否則，程序做三體運動if m3 == 0 amp。x3 = 。獲得輸入框中第二個天體的各個參數(shù) m2 = 。如果其中一個天體的質(zhì)量為0，那么程序?qū)⒆鰞审w問題的運算：function OK_Callback(hObject, eventdata, handles)hObject OK控件的句柄(詳情請參見 GCBO)eventdata 保留在MATLAB里面的參數(shù)，它的作用是方便后面的程序使用handles GUI界面的句柄，可以通過它獲取整個界面的信息（詳情請參見GUIDATA）清除坐標(biāo)系原有的內(nèi)容cla獲得程序使用者輸入的小球質(zhì)量、剛開始小球的位置與速度，設(shè)置萬有引力常量默認(rèn)值以及程序的運動時間g = *10^11。defaultUicontrolBackgroundColor39。)。defaultUicontrolBackgroundColor39。)。defaultUicontrolBackgroundColor39。)。defaultUicontrolBackgroundColor39。)。defaultUicontrolBackgroundColor39。)。defaultUicontrolBackgroundColor39。)。defaultUicontrolBackgroundColor39。)。defaultUicontrolBackgroundColor39。)。defaultUicontrolBackgroundColor39。)。defaultUicontrolBackgroundColor39。)。defaultUicontrolBackgroundColor39。)。defaultUicontrolBackgroundColor39。)。defaultUicontrolBackgroundColor39。)。defaultUicontrolBackgroundColor39。)。defaultUicontrolBackgroundColor39。)。defaultUicontrolBackgroundColor39。)。defaultUicontrolBackgroundColor39。)。defaultUicontrolBackgroundColor39。)。defaultUicontrolBackgroundColor39。)。defaultUicontrolBackgroundColor39。)。defaultUicontrolBackgroundColor39。)。 = 0。String39。,0)。set(,39。將地球的參數(shù)保存為“GUI數(shù)據(jù)”形式，以便共享 = *10^23。String39。,)。 = 0。set(,39。String39。,)。 = 0。 = 0。 = 0。圖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，速度迭代方程，位置迭代方程按直角坐標(biāo)系分解如附錄1。2 模型假設(shè)1) 忽略相對論效應(yīng)的影響；2) 不考慮其它天體對系統(tǒng)的作用，所處宇宙空間是各向同性的；3) 忽略宇宙空間粒子、潮汐力的影響；4) 假設(shè)運動過程中沒有能量損失；5) 將天體視為質(zhì)點，不考慮它們在運動的過程中彼此之間的相互碰撞；3 模型建立 具體時刻的狀態(tài)因為三個天體的位置和速度是隨時間變化的，現(xiàn)在取很短的時間間隔t，假定這段時間間隔內(nèi)系統(tǒng)參數(shù)沒有變化，則可以考慮從系統(tǒng)每一時刻的狀態(tài)計算后下一時刻的狀態(tài)。因此，除了解決一個重要的難題，更重要的是在解決問題的過程當(dāng)中產(chǎn)生了好幾種嶄新的數(shù)學(xué)思想。這三類方法各有利弊。在1744年，這個問題被Euler完全的解決了。這部分的天體力學(xué)以牛頓的萬有引力理論開始。 牛頓曾經(jīng)用幾何學(xué)原理解決了只存在相互引力作用的兩個球的運動這一兩體問題。在經(jīng)典力學(xué)，量子力學(xué)和天文學(xué)中，三體問題是一個出名的，還沒有解決的問題，它被Whittaker描述為在所有動力學(xué)問題中最出名的問題，Hilbert認(rèn)為它是一個完美數(shù)學(xué)問題的標(biāo)準(zhǔn)。在十八和十九世紀(jì)，幾乎所有著名的數(shù)學(xué)家們都曾試圖解決這個問題，但是問題的進(jìn)展不大。 本文的研究內(nèi)容本文決定通過計算機的數(shù)值運算，利用經(jīng)典力學(xué)中的牛頓三大定律和萬有引力定律迭代求解天體運動中的三體問題。因此，我們可以將時間分成許多小部分，用計算機迭代求解，得到隨時間變化的系統(tǒng)的各個參量和系統(tǒng)的末狀態(tài)。本程序的編碼和界面設(shè)計部分用MATLAB進(jìn)行開發(fā)。 = 0。 = 0。 = hObject。,0)。String39。將太陽的參數(shù)保存為“GUI數(shù)據(jù)”形式，以便共享 = *10^29。set(,39。,0)。String39。 = 0。set(,39。,0)。String39。 = 0。設(shè)置完所有屬性后，在對象的創(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。,39。設(shè)置完所有屬性后，在對象的創(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。,39。設(shè)置完所有屬性后，在對象的創(chuàng)建過程中執(zhí)行以下代碼function FirstBodyZ_CreateF(hObject, eventdata, handles)hObject FirstBodyZ 控件的句柄(詳情請參見 GCBO)eventdata 保留在MATLAB里面的參數(shù)，它的作用是方便后面的程序使用handles GUI界面的句柄，它會在調(diào)用完所有的CreateF函數(shù)后創(chuàng)建提示: 可編輯文本框控件通常在窗口上有一個白色的背景，更多詳情