Matlab模拟三个天体运动

原创

matlab科研助手 2022-03-30 10:51:14 博主文章分类：物理应用 ©著作权

©著作权归作者所有：来自51CTO博客作者matlab科研助手的原创作品，请联系作者获取转载授权，否则将追究法律责任

1 简介

这是一个三体运动模拟软件,描述自由空间中存在的两至三颗可视作质点的星体在相互之间的万有引力作用下所做的运动。基于牛顿经典力学原理，采用迭代法进行数值模拟，实时显示星体的位置以及轨迹。

2 部分代码

%
%%
function varargout = Threebody(varargin)
% THREEBODY MATLAB code for Threebody.fig
%      THREEBODY, by itself, creates a new THREEBODY or raises the existing
%      singleton*.
%
%      H = THREEBODY returns the handle to a new THREEBODY or the handle to
%      the existing singleton*.
%
%      THREEBODY('CALLBACK',hObject,eventData,handles,...) calls the local
%      function named CALLBACK in THREEBODY.M with the given input arguments.
%
%      THREEBODY('Property','Value',...) creates a new THREEBODY or raises the
    %将选中的运动类型的数据传递到三行显示栏中去
    %将代码中预的设数据写入第一行显示栏(第一颗星)
    set(handles.Edit_m1,'String',num2str(star(1).M*10^(-19)));
    set(handles.Edit_x1,'String',num2str(star(1).X));
    set(handles.Edit_y1,'String',num2str(star(1).Y));
    set(handles.Edit_z1,'String',num2str(star(1).Z));
    set(handles.Edit_u1,'String',num2str(star(1).U*10^(-3)));
    set(handles.Edit_v1,'String',num2str(star(1).V*10^(-3)));
    set(handles.Edit_w1,'String',num2str(star(1).W*10^(-3)));
    %将代码中预的设数据写入第二行显示栏(第二颗星)
    set(handles.Edit_m2,'String',num2str(star(2).M*10^(-19)));
    set(handles.Edit_x2,'String',num2str(star(2).X));
    set(handles.Edit_y2,'String',num2str(star(2).Y));
    set(handles.Edit_z2,'String',num2str(star(2).Z));
    set(handles.Edit_u2,'String',num2str(star(2).U*10^(-3)));
    set(handles.Edit_v2,'String',num2str(star(2).V*10^(-3)));
    set(handles.Edit_w2,'String',num2str(star(2).W*10^(-3)));
    %将代码中预的设数据写入第三行显示栏(第三颗星)
    set(handles.Edit_m3,'String',num2str(star(3).M*10^(-19))); 
    set(handles.Edit_x3,'String',num2str(star(3).X));     
    set(handles.Edit_y3,'String',num2str(star(3).Y));   
    set(handles.Edit_z3,'String',num2str(star(3).Z));
    set(handles.Edit_u3,'String',num2str(star(3).U*10^(-3))); 
    set(handles.Edit_v3,'String',num2str(star(3).V*10^(-3))); 
    set(handles.Edit_w3,'String',num2str(star(3).W*10^(-3)));
    %保存数据
    guidata(gcf,handles);
end
%%
% --------------------------------------------------------------------
function normal_Callback(hObject, eventdata, handles)
%%
% hObject    handle to normal (see GCBO)
% eventdata  reserved - to be defined in a future version of MATLAB
% handles    structure with handles and user data (see GUIDATA)
view([-37.5,30]); 
%%
% --------------------------------------------------------------------
function XOY_view_Callback(hObject, eventdata, handles)
%%
% hObject    handle to XOY_view (see GCBO)
% eventdata  reserved - to be defined in a future version of MATLAB
% handles    structure with handles and user data (see GUIDATA)
view([-90,90]); 
%%
% --------------------------------------------------------------------
function YOZ_view_Callback(hObject, eventdata, handles)
%%
% hObject    handle to YOZ_view (see GCBO)
% eventdata  reserved - to be defined in a future version of MATLAB
% handles    structure with handles and user data (see GUIDATA)
view([-90,0]);
%%
% --------------------------------------------------------------------
function ZOX_view_Callback(hObject, eventdata, handles)
%%
% hObject    handle to ZOX_view (see GCBO)
% eventdata  reserved - to be defined in a future version of MATLAB
% handles    structure with handles and user data (see GUIDATA)
view([180,0]);
%%
% --------------------------------------------------------------------
function Version_Callback(hObject, eventdata, handles)
%%
% hObject    handle to Version (see GCBO)
% eventdata  reserved - to be defined in a future version of MATLAB
% handles    structure with handles and user data (see GUIDATA)
hs = msgbox({'软件操作说明:'
   ''
   '    Version: 1.0 '
   ''
   '    Author:'
   ''
   '    Chen Tianyang, Li Hao'
   ''
   '    Data:2017-10-25'
   ''},'Version Information');
%改变字体大小
ht = findobj(hs, 'Type', 'text');
set(ht,'FontSize',10);
%改变对话框大小
set(hs, 'Resize', 'on'); 
%%
% --------------------------------------------------------------------
function UserGuide_Callback(hObject, eventdata, handles)
% hObject    handle to UserGuide (see GCBO)
% eventdata  reserved - to be defined in a future version of MATLAB
% handles    structure with handles and user data (see GUIDATA)
hs = msgbox({'软件操作说明:'
   ''
   ''},'UserGuide');
%改变字体大小
ht = findobj(hs, 'Type', 'text');
set(ht,'FontSize',10);
%改变对话框大小
set(hs, 'Resize', 'on');