地月系统拉格朗日点第二版

clc;clear;clear;
%%%%%%%%%%%%%%——第二版——%%%%%%%%%%%%
%上一版是比较粗略的估计算法，以地球为地月系统中心点，实际上中心点偏离地心在地球内部
% 图设置
figure('name','地月系统拉格朗日点位置');
axis equal;
grid on;
hold on
xlabel('x/m');
ylabel('y/m');
hold on
axis([-5e8,5e8,-5e8,5e8]);
hold on
m1=5.976e24;%地球质量
m2=7.35e22;%月球质量
G=6.67*10^(-11);%万有引力常数
R=3.844e8;%地月距离
a1=m2*R/(m1+m2);
a2=m1*R/(m1+m2);
omega=sqrt(G*(m1+m2)/power(R,3));%牵连角速度
L1_x=R*0.85-a1;L1_y=0;%拉格朗日点坐标
L2_x=R*1.17-a1;L2_y=0;
L3_x=-R*0.99-a1;L3_y=0;
L4_x=R/2-a1;L4_y=sqrt(3)/2*R;
L5_x=R/2-a1;L5_y=-sqrt(3)/2*R;
x=[L1_x L2_x L3_x L4_x L5_x ];
y=[L1_y L2_y L3_y L4_y L5_y ];
plot(x,y,'ro');%绘图
hold on
line([-a1,L4_x],[0,L4_y],'linestyle',':');
line([-a1,L5_x],[0,L5_y],'linestyle',':');
line([L3_x,L2_x],[L3_y,L2_y],'linestyle',':');
hold on
theta=0:pi/30:2*pi;
x1=-a1+R*cos(theta);y1=R*sin(theta);
plot(x1,y1,'b--');
title('地月系统拉格朗日点','fontsize',18);
hold on
text(L1_x,L1_y+0.4e8,'L1','fontsize',12);
text(L2_x,L2_y+0.4e8,'L2','fontsize',12);
text(L3_x+0.5e8,L3_y,'L3','fontsize',12);
text(L4_x+0.5e8,L4_y,'L4','fontsize',12);
text(L5_x+0.5e8,L5_y,'L5','fontsize',12);