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.
|
|
|
|
|
% 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>
|
|
|
|
|
|
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))].';
|
|
|
|
|
|
Rp_b=[0;0;1];
|
|
|
|
|
|
[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,Lat);
|
