matlab三体运动（仿真动画）

学学matlab动画，三体运动还是有点问题

clc;clear;close;
figure('name','三体运动');
axis equal
grid on
hold on
colordef black
view(3)
xlabel('x')
ylabel('y')
zlabel('z')
hold on
%初始化参数设置(质量、初始位置、初始速度、万有引力常数、时间步长)
m1=rand()*10^23;
m2=rand()*10^23;
m3=rand()*10^23;pos1=[rand()*10^7,rand()*10^7,rand()*10^7];
pos2=[rand()*10^7,rand()*10^7,rand()*10^7];
pos3=[rand()*10^7,rand()*10^7,rand()*10^7];v1=[rand()*1000,rand()*1000,rand()*1000];
v2=[rand()*1000,rand()*1000,rand()*1000];
v3=[rand()*1000,rand()*1000,rand()*1000];
G=6.67*10^(-11);
dt=0.005;%行星和轨迹（定义行星位置和轨迹变量）
planet1=plot3(pos1(1),pos1(2),pos1(3),'b:.','markersize',20);
planet2=plot3(pos2(1),pos2(2),pos2(3),'r:.','markersize',20);
planet3=plot3(pos3(1),pos3(2),pos3(3),'y:.','markersize',20);h1=animatedline('color','b');
h2=animatedline('color','r');
h3=animatedline('color','y');
frame=0;
while true%不断循环计算距离、相互作用力、速度、加速度、位置r12=normest(pos2-pos1);r23=normest(pos3-pos2);r13=normest(pos3-pos1);F12_val=G*m1*m2/(r12^2);F23_val=G*m2*m3/(r23^2);F13_val=G*m1*m3/(r13^2);F12_dir=(pos2-pos1)/normest(pos2-pos1);F23_dir=(pos3-pos2)/normest(pos3-pos2);F13_dir=(pos3-pos1)/normest(pos3-pos1);a1=((F12_val)*(F12_dir)+(F13_val)*(F13_dir))/m1;a2=((-F12_val)*(F12_dir)+(-F23_val)*(F23_dir))/m2;a3=((-F13_val)*(F13_dir)+(-F23_val)*(F23_dir))/m3;v1=v1+a1*dt;v2=v2+a2*dt;v3=v3+a3*dt;pos1=pos1+v1*dt+1/2*a1*dt^2;pos2=pos2+v2*dt+1/2*a2*dt^2;pos3=pos3+v3*dt+1/2*a3*dt^2;frame=frame+1;if frame==10000frame=0;% set(planet1,'position',[pos1(1),pos1(2),pos1(3)]);%set(planet3,'position',pos3(1),'YData',pos3(2),'ZData',pos3(3));addpoints(h1,pos1(1),pos1(2),pos1(3));addpoints(h2,pos2(1),pos2(2),pos2(3));addpoints(h3,pos3(1),pos3(2),pos3(3));drawnowendif r12==0||r23==0||r13==0       %判断是否相撞breakend
end

matlab三体运动（仿真动画）相关推荐

matlab仿真动画,用matlab制作简单仿真动画
用matlab制作简单仿真动画,并生成.avi格式的电影文件 MATALB知识点2008-05-01 13:17:54 阅读152 评论0 字号:大中小订阅第一种形式:利用for循环,在一定时间内控 ...
matlab行星运动轨迹仿真动画,利用Matlab可视化功能实现微分方程求解行星运动轨迹...
利用Matlab可视化功能实现微分方程求解行星运动轨迹 1.背景在物理学璀璨的发展史上,物理学家花了很长时间研究我们头顶浩瀚的星空,试图探究星星的运行模式,以及地球自身的运动模式.其中不乏像亚里士多 ...
MATLAB | 经典力学框架下的三体运动数值模拟软件
经典力学框架下的三体运动数值模拟软件 1 背景说明 2 实验目的 3 基本原理及思路 4 成果展示 4.1 用户界面展示 4.2 所需文件 4.3 软件运行效果 4.3.1 不等质双星 4.3.2 近 ...
matlab 图片制作动画制作,MATLAB作图之制作动画：单摆运动仿真
MA TLAB作图之制作动画:单摆运动仿真%制作动画 %挂摆横梁 plot([-0.2;0.2],[0;0],'color','y','linestyle','-',... 'linewidth',1 ...
MATLAB 模拟三体运动教程
MATLAB 模型三体运动Three_body_movement(无行星) clear all close all clc%--------data--------% m1 = rand() * 10 ...
matlab模拟三体运动_matlab三体运动
闲的蛋疼,学学matlab动画,三体运动还是有点问题(提示line没xdate属性,脑瓜也疼了),不多废话上代码 clc;clear;close; figure('name','三体运动'); axi ...
matlab与vrep联合仿真,基于MATLAB与V-REP的机器人加工轨迹生成与运动仿真.PDF
基于MATLAB与V-REP的机器人加工轨迹生成与运动仿真第卷第期 ( ) 57 5 厦门大学学报自然科学版 Vol.57 No.5 年月 ( ) 2018 9 Se.2018 Journ ...
matlab平行因子_基于MATLAB某客车悬置系统优化与运动仿真
本文对一CNG燃气客车悬置系统建立6阶振动数学模型,基于MATLAB分析振动频率与解耦率并对系统固有频率的配置以及解耦率设计优化程序,运用NX CAE运动仿真校核MATLAB分析结果的准确性,为悬置系 ...
【机械仿真】基于matlab GUI曲柄摇杆机构运动仿真【含Matlab源码 1608期】
一.获取代码方式获取代码方式1: 完整代码已上传我的资源:[机械仿真]基于matlab GUI曲柄摇杆机构运动仿真[含Matlab源码 1608期] 点击上面蓝色字体,直接付费下载,即可. 获取代码 ...

matlab三体运动（仿真动画）

matlab三体运动（仿真动画）相关推荐

最新文章

热门文章