matlab模拟三体运动_Matlab模拟三体运动

m1=1.5e30;%1号星的质量-红

m2=3.01e30;%2号星的质量-蓝

m3=3.5e30;%3号星的质量-绿

h1=animatedline('MaximumNumPoints',2000,'color','r');%动画

h2=animatedline('MaximumNumPoints',2000,'color','b');

h3=animatedline('MaximumNumPoints',2000,'color','g');

G=6.67e-11;%引力常量

v1=[20,40];%三颗恒星的初始速度

v2=[-30,20];

v3=[27,-30];

pos1=[0,1.5e18];%三颗恒星的坐标

pos2=[1.5e18,3e18];

pos3=[3e18,0];

delt=0.1;%观测时间间隔

for i=1:10000

addpoints(h1,pos1(1),pos1(2));

addpoints(h2,pos2(1),pos2(2));

addpoints(h3,pos3(1),pos3(2));

drawnow;

pos1=pos1+v1*delt;%更新位置

pos2=pos2+v2*delt;

pos3=pos3+v3*delt;

r12=normest(pos2-pos1);

r23=normest(pos3-pos2);

r31=normest(pos1-pos3);

f12_len=G*m1*m2/(r12^2);

f23_len=G*m2*m3/(r23^2);

f31_len=G*m3*m1/(r31^2);

f12_dir=(pos2-pos1)/r12;%计算万有引力

f23_dir=(pos3-pos2)/r23;

f31_dir=(pos1-pos3)/r31;

f12=f12_len*f12_dir;

f23=f23_len*f23_dir;

f31=f31_len*f31_dir;

v1=v1+(f12-f31)*delt;

v2=v2+(f23-f12)*delt;

v3=v3+(f31-f23)*delt;

end

clear;

matlab模拟三体运动_Matlab模拟三体运动相关推荐

三体运动计算机模拟软件,三体运动模拟软件ThreeBody
这是三体运动模拟软件ThreeBody,是一款三体运动模拟软件. 软件介绍三体运动模拟软件ThreeBody,可以看三体运动,程序启动后,会出现三个随机大小的球体在运动. 使用说明 1.打开已有的一 ...
三体运动python模拟（代码能直接运行）
三体运动python模拟(代码能直接运行) 原理说明程序代码(可直接运行) 代码结果示例原理说明很高兴撰写CSDN首篇文章!写了一个简单的程序对三体运动进行了模拟以验证其运动的复杂性.难以预测性 ...
Matlab模拟光粒攻击前后的三体运动
数模国赛将至,但是我的matlab水平依旧很烂,三体倒是全都看完了.我于是想到:对三体中感兴趣的情节(光粒攻击前后的三体运动)建模,然后用matlab实现出来,我不就既可以练习建模,又可以满足我贪玩的 ...
boid模型的Matlab程序,基于Boid模型以及吸引—排斥模型的沙丁鱼集群运动行为模拟...
刘培宫子涵 [摘要]本文首先为了分析了沙丁鱼集群的运动模式,引入了 Boid 模型.通过个体间的距离划分出了排斥区域.一致区域和吸引区域.当相邻个体处于排斥区域则保持足够的距离,尽量避免碰撞,当处 ...
【Python 3.7】分子运动：修改 rw_visual.py，将其中的 plt.scatter() 替换为 plt.plot() 。为模拟花粉在水滴表面的运动路径
[Python 3.7]}分子运动:修改 rw_visual.py,将其中的 plt.scatter() 替换为 plt.plot() .为模拟花粉在水滴表面的运动路径,向 plt.plot() 传 ...
如何使用matlab仿真物理实验,基于Matlab的大学物理实验模拟.pdf
维普资讯浙江科技学院学报,第 2O卷第 2期,2008年 6月 JournalofZhejiangUniversityofScienceandTechnology Vo1．20No．2,June20 ...
Algorithm之MC：基于Matlab实现通过蒙特卡洛方法模拟二维布朗运动
Algorithm之MC:基于Matlab实现通过蒙特卡洛方法模拟二维布朗运动目录输出结果设计代码输出结果基于Matlab实现通过蒙特卡洛方法模拟二维布朗运动设计代码 %基于Matlab实 ...
matlab 高斯函数_光学与matlab：高斯光的模拟
光学与matlab:高斯光的模拟高斯光束,自然何为高斯光束?百度了一下,如下所示,"通常情形,激光谐振腔发出的基模辐射场,其横截面的振幅分布遵守高斯函数,故称高斯光束."如果您对 ...
matlab中eacf函数,基于MATLAB的切比雪夫I型模拟低通滤波器设计
基于MATLAB的切比雪夫I型模拟低通滤波器设计课程设计名称:数字信号处理课程设计专业班级 : 电信0604 学生姓名 : 学号 : 20064300430 指导教师 : 课程设计时间:2009.6 ...

matlab模拟三体运动_Matlab模拟三体运动

matlab模拟三体运动_Matlab模拟三体运动相关推荐

最新文章

热门文章