【正文】 BackgroundColor39。改變觀察視角view([az,el])。String39。BackgroundColor39。將字符轉(zhuǎn)換為雙精度數(shù)az=str2double(get(,39。)) set(hObject,39。))。同時，將字符轉(zhuǎn)換為雙精度數(shù)az=str2double(get(hObject,39。設置完所有屬性后，在對象的創(chuàng)建過程中執(zhí)行以下代碼function azSlider_CreateF(hObject, eventdata, handles)hObject azSlider控件的句柄(詳情請參見 GCBO)eventdata 保留在MATLAB里面的參數(shù)，它的作用是方便后面的程序使用handles GUI界面的句柄，它會在調(diào)用完所有的CreateF函數(shù)后創(chuàng)建提示: 滑塊控件通常有一個淺灰色的背景if isequal(get(hObject,39。)。ydata39。,x2/10^10,39。xdata39。az3=(f23(3)f31(3))/m3。y1=y1+v1*t+1/2*ay1*t^2。x2=x2+u2*t+1/2*ax2*t^2。，用k來控制迭代進行的次數(shù)for k=0::time設置步長t=100000。.39。,39。,[0 0 1],39。MarkerSize39。,z3/10^10)。zdata39。z3=z3+w3*t+1/2*az3*t^2。ay3=f31(2)/m3。萬有引力定律 f31=g*m3*m1/(sqrt((x3x1)^2+(y3y1)^2+(z3z1)^2))^3*[x3x1,y3y1,z3z1]。,39。normal39。Marker39。,z3/10^10)。zdata39。z3=z3+w3*t+1/2*az3*t^2。ay3=f23(2)/m3。萬有引力定律 f23=g*m2*m3/(sqrt((x2x3)^2+(y2y3)^2+(z2z3)^2))^3*[x2x3,y2y3,z2z3]。,39。normal39。Marker39。,z2/10^10)。zdata39。z2=z2+w2*t+1/2*az2*t^2。ay2=f12(2)/m2。萬有引力定律 f12=g*m1*m2/(sqrt((x1x2)^2+(y1y2)^2+(z1z2)^2))^3*[x1x2,y1y2,z1z2]。,39。normal39。Marker39。 加網(wǎng)格線grid on坐標軸顯示的空間范圍axis([100 100 100 100 100 100])當m1,m2,m3其中有一個為0時，程序做兩體運動。獲得輸入框中第三個天體的各個參數(shù) m3 = 。w1 = 。這里由于坐標值的數(shù)量級是1010，x，y，z坐標軸的范圍為100~100，所以在畫圖的時候，計算出來的各個坐標值都要除以1010。), get(0,39。String39。), get(0,39。String39。), get(0,39。String39。), get(0,39。String39。), get(0,39。String39。), get(0,39。String39。), get(0,39。String39。), get(0,39。String39。), get(0,39。String39。), get(0,39。String39。), get(0,39。String39。), get(0,39。String39。), get(0,39。String39。), get(0,39。String39。), get(0,39。String39。), get(0,39。String39。), get(0,39。String39。), get(0,39。String39。), get(0,39。String39。), get(0,39。String39。), get(0,39。String39。 = 0。set(,39。String39。在圖形界面的第三個天體的文本輸入框內(nèi)顯示默認的輸入值，此輸入值為月球的參數(shù)。,0)。set(,39。String39。 = 0。,0)。set(,39。String39。 = 0。 = 0。 = 0。在萬有引力f1f2f31作用下，天體天體天體3的加速度分別為f12+f31m1，f12f23m2，f23f31m3。同時，還可以在不同的觀察角度了解它們的運動過程。多年來，人們對這個問題進行了大量的研究，從1750年到20世紀初，關于這個問題有800多篇援引了許多杰出的數(shù)學家和天文學家成果的論文發(fā)表。探究三體問題的方式大抵能夠分為三類：第一類是分析方法，其基本原理是把天體的坐標和速度展開為時間或其它小參數(shù)的級數(shù)形式的近似解析表達式，討論天體的坐標和軌道要素隨時間的變化規(guī)律；第二類是定性方法，選用微分方程的定性理論來研究很長一段時間內(nèi)三體運動的宏觀運動規(guī)律和全局性質(zhì)；第三類是數(shù)值方法，這是直接依據(jù)微分方程的計算方式得出天體在特定時刻的詳細位置和速度。Bernoulli因為他關于分析處理兩體問題的研究獲得了法國科學院獎。因此，可以用現(xiàn)代動力天文學的基礎理論來研究這個問題。在1710年，Johann天體物理中的三體問題可以簡單的被描述為：三個天體在他們之間的相互之間的萬有引力作用下，給出各個天體的初始條件，包括三個天體的質(zhì)量、初始位置和初始速度，確定它們的后續(xù)運動，包括位置和速度隨時間的變化。截止到目前為止，三體問題特解的族數(shù)被擴充到了16族。同時，本文利用MATLAB的GUI功能，設計出簡潔友好的界面，并且通過動畫的方法將天體的位置和速度在空間直角坐標系上實時的表示出來，作出了天體物理實驗用計算機模擬來進行的嘗試。 受力分析如圖31所示，每個天體都受到其他兩個天體的萬有引力作用。（2）三體問題的編碼1o三個天體各個參數(shù)的初始化部分。 = 0。初始化第三個天體的質(zhì)量，位置和速度 = 0。此前指令更改了handles，必須靠guidata指令才能把更新了的handles加以保存，以供后續(xù)使用guidata(hObject, handles)。set(,39。,0)。 = 0。String39。set(,39。,)。 = 0。String39。set(,39。,0)。此前指令更改了handles，必須靠guidata指令才能把更新了的handles加以保存，以供后續(xù)使用guidata(hObject,handles)。amp。white39。amp。white39。amp。white39。amp。white39。amp。white39。amp。white39。amp。white39。amp。white39。amp。white39。amp。white39。amp。white39。amp。white39。amp。white39。amp。white39。amp。white39。amp。white39。amp。white39。amp。white39。amp。white39。amp。white39。amp。white39。z1 = 。u2 = 。v3 = 。 m2 ~=0設置天體顏色、點型、大小等參數(shù)，每次循環(huán)都重新設置整個畫面h=line(39。,50,39。Marker39。normal39。速度值迭代u1=u1+ax1*t。 以下為對z軸方向的計算az1=f12(3)/m1。,x1/10^10,39。ydata39。 m2 ~=0設置天體顏色、點型、大小等參數(shù)，每次循環(huán)都重新設置整個畫面i=line(39。,40,39。Marker39。normal39。速度值迭代u2=u2+ax2*t。 以下為對z軸方向的計算az2=f23(3)/m2。,x2/10^10,39。ydata39。 m1 ~=0設置天體顏色、點型、大小等參數(shù)，每次循環(huán)都重新設置整個畫面h=line(39。,50,39。Marker39。normal39。速度值迭代u1=u1+ax1*t。 以下為對z軸方向的計算az1=+f31(3)/m1。,x1/10^10,39。ydata39。,39。)。MarkerSize39。,[0 1 0],39。,39。ax2=(f12(1)f23(1))/m2。 以下為對y軸方向的計算ay1=(f12(2)+f31(2))/m1。v3=v3+ay3*t。w2=w2+az2*t。,z1/10^10)。set(l,39。關閉界面close。String39。BackgroundColor39。,az)。BackgroundColor39。end用于控制仰角的滑鍵的回調(diào)函數(shù)，用于同步仰角文本輸入框中顯示的數(shù)值，同時改變觀察者的仰角：function elSlider_Callback(hObject, eventdata, handles)hObject elSlider控件的句柄(詳情請參見 GCBO)eventdata 保留在MATLAB里面的參數(shù)，它的作用是方便后面的程序使用handles GUI界面的句柄，可以通過它獲取整個界面的信息（詳情請參見GUIDATA）獲得當前滑鍵的值el=get(hObject,39。,num2str(el))。,[.9 .9 .9])。獲得當前的方位角az=str2double(get(,39。), get(0,39。)。 isequal(get(hObject,39。Value39。)) set(hObject,39。改變仰角輸入文本框的顯示數(shù)值，此數(shù)值為當前滑鍵所代表的值set(,39。white39。amp。改變方位角滑鍵的值set(,39。defaultUicontrolBackgroundColor39。))。,z3/10^10)。zdata39。,y1/10^10,39。z3=z3+w3*t+1/2*az3*t^2。 v1=v1+ay1*t。u2=u2+ax2*t。 f31=g*m3*m1/(sqrt((x3x1)^2+(y3y1)^2+(z3z1)^2))^3*[x3x1,y3y1,z3z1]。,40,39。l=line(39。.39。,39。,[1 0 0],39。xdata39。 重置三個小球的位置set(h,39。 v1=v1+ay1*t。坐標值迭代x1=x1+u1*t+1/2*ax1*t^2。erasemode39。Color39。,39。 m3 ~= 0 amp。xdata39。 重置三個小球的位置set(i,39。 v2=v2+ay2*t。坐標值迭代x2=x2+u2*t+1/2*ax2*t^2。erasemode39。Color39。,39。 m3 ~= 0 amp。xdata39。 重置兩個小球的位置set(h,39。 v1=v1+ay1*t。坐標值迭代x1=x1+u1*t+1/2*ax1*t^2。erasemode39。Color39。,39。 m1 ~= 0 amp。z3 = 。y2 = 。x1 = 。BackgroundColor39。此前指令更改了handles，必須靠guidata指令才能把更新了的handles加以保存，以供后續(xù)使用guidata(hObject,handles)。BackgroundColor39。此前指令更改了handles，必須靠guidata指令才能把更新了的handles加以保存，以供后續(xù)使用guidata(hObject,handles)。BackgroundColor39。此前指令更改了handles，必須靠guidata指令才能把更新了的handles加以保存，以供后續(xù)使用guidata(hObject,handles)。BackgroundColor39。此前指令更改了handles，必須靠guidata指令才能把更新了的handles加以保存，以供后續(xù)使用guidata(hObject,handles)。BackgroundColor39。此前指令更改了handles，必須靠guidata指令才能把更新了的handles加以保存，以供后續(xù)使用guidata(hObject,handles)。BackgroundColor39。此前指令更改了handles，必須靠guidata指令才能把更新了的handles加以保存，以供后續(xù)使用guidata(hObject,handles)。BackgroundColor39。此前指令更改了handles，必須靠guidata指令才能把更新了的handles加以保存，以供后續(xù)使用guidata(hObject,handles)。BackgroundColor39。此前指令更改了handles，必須靠guidata指令才能把更新了的handles加以保存，以供后續(xù)使用guidata(hObject,handles)。BackgroundColor39。此前指令更改了handles，必須靠guidata指令才能把更新了的handles加以保存，以供后續(xù)使用guidata(hObject,handles)。BackgroundColor39。此前指令更改了handles，必須靠guidata指令才能把更新了的handles加以保存，以供后續(xù)使用guidata(hObject,handles)。BackgroundColor39。此前指令更改了handles，必須靠guidata指令才能把更新了的handles加以保存，以供后續(xù)使用guidata(hObject,handles)。BackgroundColor39。此前指令更改了handles，必須靠guidata指令才能把更新了的handles加以保存，以供后續(xù)使用guidata(hObject,handles)。BackgroundColor39。此前指令更改了handles，必須靠guidata指令才能把更新了的hand