You cannot select more than 25 topics Topics must start with a letter or number, can include dashes ('-') and can be up to 35 characters long.

25 lines
750 B
Plaintext

% function check_main()
clc
clear
close all
[AOC_t0,AOC_fai,AOC_thita,AOC_pusai,AOC_a0,AOC_e0,AOC_i0,AOC_OMG0,AOC_os0,AOC_u0,AOC_Tstep] = readInput();
AOC_dTime=86400;
AOC_Tin=0:AOC_Tstep:AOC_dTime;
Mu=3.986005e5; %<25><><EFBFBD><EFBFBD><EFBFBD><EFBFBD><EFBFBD><EFBFBD><EFBFBD><EFBFBD><EFBFBD><EFBFBD>
w0=(Mu/AOC_a0^3)^0.5; %<25><><EFBFBD>ǹ<EFBFBD><C7B9><EFBFBD><EFBFBD><EFBFBD><EFBFBD>ٶ<EFBFBD>
we=7.29211585530487e-5; %<25><><EFBFBD><EFBFBD><EFBFBD><EFBFBD>ת<EFBFBD><D7AA><EFBFBD>ٶ<EFBFBD>
Re=6371.23; %<25><><EFBFBD><EFBFBD>ƽ<EFBFBD><C6BD><EFBFBD>
Len=length(AOC_Tin);
Thita=linspace(-8.6,8.6,Len);
for num=1:Len
Tin=AOC_Tin(num);
AOC_u=AOC_u0+we*Tin;
Rp_b=[sind(Thita(num)),0,cosd(Thita(num))].';
[Tmp_Lon,Tmp_Lat]=SightPoint(Tin,AOC_fai,AOC_thita,AOC_pusai,AOC_a0,AOC_e0,AOC_i0,AOC_OMG0,AOC_os0,AOC_u,Rp_b);
Lon(num,1)=Tmp_Lon;
Lat(num,1)=Tmp_Lat;
end
saveData(Lon,)