【正文】 環(huán)都重新設(shè)置整個(gè)畫面h=line(39。Color39。,[1 0 0],39。Marker39。,39。.39。,39。MarkerSize39。,50,39。erasemode39。,39。normal39。)。l=line(39。Color39。,[0 1 0],39。Marker39。,39。.39。,39。MarkerSize39。,40,39。erasemode39。,39。normal39。)。，用k來控制迭代進(jìn)行的次數(shù)for k=0::time設(shè)置步長t=100000。萬有引力定律 f31=g*m3*m1/(sqrt((x3x1)^2+(y3y1)^2+(z3z1)^2))^3*[x3x1,y3y1,z3z1]。 以下為對x軸方向的計(jì)算牛頓第二定律ax1=+f31(1)/m1。ax3=f31(1)/m3。坐標(biāo)值迭代x1=x1+u1*t+1/2*ax1*t^2。x3=x3+u3*t+1/2*ax3*t^2。速度值迭代u1=u1+ax1*t。u3=u3+ax3*t。 以下為對y軸方向的計(jì)算ay1=+f31(2)/m1。ay3=f31(2)/m3。y1=y1+v1*t+1/2*ay1*t^2。y3=y3+v3*t+1/2*ay3*t^2。 v1=v1+ay1*t。v3=v3+ay3*t。 以下為對z軸方向的計(jì)算az1=+f31(3)/m1。az3=f31(3)/m3。 z1=z1+w1*t+1/2*az1*t^2。z3=z3+w3*t+1/2*az3*t^2。 w1=w1+az1*t。w3=w3+az3*t。 重置三個(gè)小球的位置set(h,39。xdata39。,x1/10^10,39。ydata39。,y1/10^10,39。zdata39。,z1/10^10)。set(l,39。xdata39。,x3/10^10,39。ydata39。,y3/10^10,39。zdata39。,z3/10^10)。 刷新畫面drawnow 控制速度pause()endelse 設(shè)置天體顏色、點(diǎn)型、大小等參數(shù)，每次循環(huán)都重新設(shè)置整個(gè)畫面h=line(39。Color39。,[1 0 0],39。Marker39。,39。.39。,39。MarkerSize39。,50,39。erasemode39。,39。normal39。)。i=line(39。Color39。,[0 0 1],39。Marker39。,39。.39。,39。MarkerSize39。,40,39。erasemode39。,39。normal39。)。l=line(39。Color39。,[0 1 0],39。Marker39。,39。.39。,39。MarkerSize39。,40,39。erasemode39。,39。normal39。)。，用k來控制迭代進(jìn)行的次數(shù)for k=0::time設(shè)置步長t=100000。萬有引力定律 f12=g*m1*m2/(sqrt((x1x2)^2+(y1y2)^2+(z1z2)^2))^3*[x1x2,y1y2,z1z2]。 f23=g*m2*m3/(sqrt((x2x3)^2+(y2y3)^2+(z2z3)^2))^3*[x2x3,y2y3,z2z3]。 f31=g*m3*m1/(sqrt((x3x1)^2+(y3y1)^2+(z3z1)^2))^3*[x3x1,y3y1,z3z1]。 以下為對x軸方向的計(jì)算牛頓第二定律ax1=(f12(1)+f31(1))/m1。ax2=(f12(1)f23(1))/m2。ax3=(f23(1)f31(1))/m3。坐標(biāo)值迭代x1=x1+u1*t+1/2*ax1*t^2。x2=x2+u2*t+1/2*ax2*t^2。x3=x3+u3*t+1/2*ax3*t^2。速度值迭代u1=u1+ax1*t。u2=u2+ax2*t。u3=u3+ax3*t。 以下為對y軸方向的計(jì)算ay1=(f12(2)+f31(2))/m1。ay2=(f12(2)f23(2))/m2。ay3=(f23(2)f31(2))/m3。y1=y1+v1*t+1/2*ay1*t^2。y2=y2+v2*t+1/2*ay2*t^2。y3=y3+v3*t+1/2*ay3*t^2。 v1=v1+ay1*t。v2=v2+ay2*t。v3=v3+ay3*t。 以下為對z軸方向的計(jì)算az1=(f12(3)+f31(3))/m1。az2=(f12(3)f23(3))/m2。az3=(f23(3)f31(3))/m3。 z1=z1+w1*t+1/2*az1*t^2。z2=z2+w2*t+1/2*az2*t^2。z3=z3+w3*t+1/2*az3*t^2。 w1=w1+az1*t。w2=w2+az2*t。w3=w3+az3*t。 重置三個(gè)小球的位置set(h,39。xdata39。,x1/10^10,39。ydata39。,y1/10^10,39。zdata39。,z1/10^10)。set(i,39。xdata39。,x2/10^10,39。ydata39。,y2/10^10,39。zdata39。,z2/10^10)。set(l,39。xdata39。,x3/10^10,39。ydata39。,y3/10^10,39。zdata39。,z3/10^10)。 刷新畫面drawnow 控制速度pause()endend4o關(guān)閉按鈕的回調(diào)函數(shù)，用來結(jié)束程序，關(guān)閉圖形界面function Close_Callback(hObject, eventdata, handles)hObject Close控件的句柄(詳情請參見 GCBO)eventdata 保留在MATLAB里面的參數(shù)，它的作用是方便后面的程序使用handles GUI界面的句柄，可以通過它獲取整個(gè)界面的信息（詳情請參見GUIDATA）清除清空當(dāng)前畫圖窗口，以方便畫新的圖形clf。關(guān)閉界面close。5o改變觀察者的視角，可以從各個(gè)角度觀察三體的運(yùn)行狀態(tài)用于控制方向角的滑鍵的回調(diào)函數(shù)，用于同步方位角文本輸入框中顯示的數(shù)值，同時(shí)改變觀察者的觀察方位角：function azSlider_Callback(hObject, eventdata, handles)hObject azSlider控件的句柄(詳情請參見 GCBO)eventdata 保留在MATLAB里面的參數(shù)，它的作用是方便后面的程序使用handles GUI界面的句柄，可以通過它獲取整個(gè)界面的信息（詳情請參見GUIDATA）獲得當(dāng)前滑鍵的值az=get(hObject,39。Value39。)。將字符轉(zhuǎn)換為雙精度數(shù)el=str2double(get(,39。String39。))。改變方位角輸入文本框的顯示數(shù)值，此數(shù)值為當(dāng)前滑鍵所代表的值set(,39。String39。,num2str(az))。改變觀察視角view([az,el])。設(shè)置完所有屬性后，在對象的創(chuàng)建過程中執(zhí)行以下代碼function azSlider_CreateF(hObject, eventdata, handles)hObject azSlider控件的句柄(詳情請參見 GCBO)eventdata 保留在MATLAB里面的參數(shù)，它的作用是方便后面的程序使用handles GUI界面的句柄，它會(huì)在調(diào)用完所有的CreateF函數(shù)后創(chuàng)建提示: 滑塊控件通常有一個(gè)淺灰色的背景if isequal(get(hObject,39。BackgroundColor39。), get(0,39。defaultUicontrolBackgroundColor39。)) set(hObject,39。BackgroundColor39。,[.9 .9 .9])。end用于控制方位角的文本輸入框的回調(diào)函數(shù)，用于同步方位角滑鍵的值，同時(shí)改變觀察者的方位角：function azEditText_Callback(hObject, eventdata, handles)hObject azEditText控件的句柄(詳情請參見 GCBO)eventdata 保留在MATLAB里面的參數(shù)，它的作用是方便后面的程序使用handles GUI界面的句柄，可以通過它獲取整個(gè)界面的信息（詳情請參見GUIDATA）獲得當(dāng)前輸入框的輸入字符。同時(shí)，將字符轉(zhuǎn)換為雙精度數(shù)az=str2double(get(hObject,39。String39。))。改變方位角滑鍵的值set(,39。Value39。,az)。獲得當(dāng)前的仰角el=str2double(get(,39。String39。))。改變觀察視角view([az,el])。設(shè)置完所有屬性后，在對象的創(chuàng)建過程中執(zhí)行以下代碼function azEditText_CreateF(hObject, eventdata, handles)hObject azEditText控件的句柄(詳情請參見 GCBO)eventdata 保留在MATLAB里面的參數(shù)，它的作用是方便后面的程序使用handles GUI界面的句柄，它會(huì)在調(diào)用完所有的CreateF函數(shù)后創(chuàng)建提示: 可編輯文本框控件通常在窗口上有一個(gè)白色的背景，更多詳情請參見ISPC 和 COMPUTER.if ispc amp。amp。 isequal(get(hObject,39。BackgroundColor39。), get(0,39。defaultUicontrolBackgroundColor39。)) set(hObject,39。BackgroundColor39。,39。white39。)。end用于控制仰角的滑鍵的回調(diào)函數(shù)，用于同步仰角文本輸入框中顯示的數(shù)值，同時(shí)改變觀察者的仰角：function elSlider_Callback(hObject, eventdata, handles)hObject elSlider控件的句柄(詳情請參見 GCBO)eventdata 保留在MATLAB里面的參數(shù)，它的作用是方便后面的程序使用handles GUI界面的句柄，可以通過它獲取整個(gè)界面的信息（詳情請參見GUIDATA）獲得當(dāng)前滑鍵的值el=get(hObject,39。Value39。)。將字符轉(zhuǎn)換為雙精度數(shù)az=str2double(get(,39。String39。))。改變仰角輸入文本框的顯示數(shù)值，此數(shù)值為當(dāng)前滑鍵所代表的值set(,39。String39。,num2str(el))。改變觀察視角view([az,el])。設(shè)置完所有屬性后，在對象的創(chuàng)建過程中執(zhí)行以下代碼 function elSlider_CreateF(hObject, eventdata, handles)hObject elSlider控件的句柄(詳情請參見 GCBO)eventdata 保留在MATLAB里面的參數(shù)，它的作用是方便后面的程序使用handles GUI界面的句柄，它會(huì)在調(diào)用完所有的CreateF函數(shù)后創(chuàng)建提示: 滑塊控件通常有一個(gè)淺灰色的背景if isequal(get(hObject,39。BackgroundColor39。), get(0,39。defaultUicontrolBackgroundColor39。)) set(hObject,39。BackgroundColor39。,[.9 .9 .9])。end用于控制仰角的文本輸入框的回調(diào)函數(shù)，用于同步仰角滑鍵的值，同時(shí)改變觀察者仰角：function elEditText_Callback(hObject, eventdata, handles)hObject elEditText控件的句柄(詳情請參見 GCBO)eventdata 保留在MATLAB里面的參數(shù)，它的作用是方便后面的程序使用handles GUI界面的句柄，可以通過它獲取整個(gè)界面的信息（詳情請參見GUIDATA）獲得當(dāng)前輸入框的輸入字符。同時(shí)，將字符轉(zhuǎn)換為雙精度數(shù)el=str2double(get(hObject,39。String39。))。改變仰角滑鍵的值set(,39。Value39。,el)。獲得當(dāng)前的方位角az=str2double(get(,39。String39。))。改變觀察視角view([az,el])。設(shè)置完所有屬性后，在對象的創(chuàng)建過程中執(zhí)行以下代碼function elEditText_CreateF(hObject, eventdata, handles)hObject elEditText控件的句柄(詳情請參見 GCBO)eventdata 保留在MATLAB里面的參數(shù)，它的作用是方便后面的程序使用handles GUI界面的句柄，它會(huì)在調(diào)用完所有的CreateF函數(shù)后創(chuàng)建提示: 可編輯文本框控件通常在窗口上有一個(gè)白色的背景，更多詳情請參見ISPC 和 COMPUTER.if ispc amp。amp。 isequal(get(hObject,39。BackgroundColor39。), get(0,39。defaultUicontrolBackgroundColor39。)) set(hObject,39。BackgroundColor39。,39。white39。)。end用于恢復(fù)默認(rèn)方位角和仰角的默認(rèn)值按鈕的回調(diào)函數(shù)，仰角恢復(fù)為默認(rèn)值300，同時(shí)更新方位角滑鍵的值，方位