funcofmu.m

太阳—地月系平衡解,5个平衡解

close all;clear all;clc;
format long;
%% --------------------地月系平动解
r=1;
global mu
mu=0.01215057;


s=[mu*r,0,0];
e=[-(1-mu)*r,0,0];
%------------------------------------%------------------------------------
%------------------------------------%------------------------------------

plot3(mu*r,0,0,'rO');
axis([-1.5,1.5,-1.5,1.5]);axis equal;
hold on;
%------------------------------------%------------------------------------
%------------------------------------%------------------------------------
plot3(-(1-mu)*r,0,0,'kh');

L1=fsolve('myxyz',[-1.0 0 0]);
plot3(L1(1),L1(2),L1(3),'m*');
L2=fsolve('myxyz',[-0 0 0]);
plot3(L2(1),L2(2),L2(3),'c*');
L3=fsolve('myxyz',[1.0 0 0]);
plot3(L3(1),L3(2),L3(3),'y*');
L4=fsolve('myxyz',[-0.5+mu 0.866 0]);
plot3(L4(1),L4(2),L4(3),'b*');
L5=fsolve('myxyz',[-0.5+mu -0.866 0]);
plot3(L5(1),L5(2),L5(3),'g*');


grid on;
legend('Earth','Moon','L1','L2','L3','L4','L5');
title('地月系5个平动点（mu=0.01215057）');

text(mu*r-0.1,0+0.2,0,' Earth');
text(-(1-mu)*r-0.1,0+0.2,0,'Moon');
text(L1(1)-0.2,L1(2),L1(2),'L1');
text(L2(1)+0.1,L2(2),L2(2),'L2');
text(L3(1)+0.1,L3(2),L3(2),'L3');
text(L4(1)+0.1,L4(2),L4(2),'L4');
text(L5(1)+0.1,L5(2),L5(2),'L5');