計算

1 次查看（过去 30 天）

显示更早的评论

高倉 2023-3-19

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1931480-

移动：Atsushi Ueno 2023-3-21

在 MATLAB Online 中打开

ISSとFENGYUN 1C DEB が最も接近する直前， ISSは日本上空を通過していることがわかった.

東京(N35.68944, E139.69167)の基地局から ISSとの通信リンクは確立していたか

(通信できる状態であったか)確認しなさい. 通信可能な開始時間と終了時間を求めなさい.

startTime = datetime(2022,12,08,00,00,0);
stopTime = startTime + hours(12);
sampleTime = 5;
sc = satelliteScenario(startTime,stopTime,sampleTime);
sat = satellite(sc,"iss_orbit.tle");
rc=255/255;
gc=1/255;
bc=1/255;
sat(1).MarkerColor = [rc gc bc];
show(sat);
names = sat.Name + " Camera";
cam = conicalSensor(sat,"Name",names,"MaxViewAngle",90);
% groundTrack(sat, "LeadTime", 200);
groundTrack(sat(1), LeadTime=200, ...
    TrailTime=100, ...
    LineWidth=3, LeadLineColor=[0 1 0], TrailLineColor=[1 0 1]);
ele1 = orbitalElements(sat(1));
ele2 = orbitalElements(sat(2));
ele3 = orbitalElements(sat(3));
disp(ele3);
time = datetime(2022,12,08,00,00,0);
[pos, vel] = states(sat(1),time,"CoordinateFrame","geographic");
disp(pos);
disp(vel);
name = ["Madrid Deep Space Communications Complex", ...
    "Canberra Deep Space Communications Complex"];
lat = [40.43139, -35.40139];
lon = [-4.24806, 148.98167];
gs = groundStation(sc,"Name",name,"Latitude",lat, ...
    "Longitude", lon, ...
    "minElevationAngle", 0);
%access setting
ac1 = access(gs(1), sat);
ac2 = access(gs(2), sat);
% ac = access(gs(2), sat(1));
%interval status
intvls = accessIntervals(ac1);
disp(intvls);
%display FOV of camera
fieldOfView(cam);
time = datetime(2020,6,02,12,30,0);
[az,elev,r] = aer(gs(1),sat(1),time);
disp(az),
disp(elev);
disp(r);
play(sc);

このコードを活用して，軌道解析する

2 个评论
显示无隐藏无

Atsushi Ueno 2023-3-21

移动：Atsushi Ueno 2023-3-21

在 MATLAB Online 中打开

urlwrite("https://celestrak.org/NORAD/elements/stations.txt","iss_orbit.tle");
startTime = datetime(2023,3,5,22,49,45); % startTime = datetime(2022,12,08,00,00,0);
stopTime = startTime + minutes(3);       % stopTime = startTime + hours(12);
sampleTime = 5;                          % sampleTime = 5;
sc = satelliteScenario(startTime,stopTime,sampleTime);
sat = satellite(sc,"iss_orbit.tle");
sat(1).MarkerColor = [255/255 1/255 1/255];
names = sat.Name + " Camera";
cam = conicalSensor(sat,"Name",names,"MaxViewAngle",90);
groundTrack(sat(1), LeadTime=200, TrailTime=100, ...
    LineWidth=3, LeadLineColor=[0 1 0], TrailLineColor=[1 0 1]);
ele1 = orbitalElements(sat(1));
ele2 = orbitalElements(sat(2));
ele3 = orbitalElements(sat(3));
[pos, vel] = states(sat(1),"CoordinateFrame","geographic");
name = ["Madrid Deep Space Communications Complex", ...
    "Canberra Deep Space Communications Complex"];
lat = [40.43139, -35.40139];
lon = [-4.24806, 148.98167];
gs = groundStation(sc,"Name",name,"Latitude",lat, ...
    "Longitude", lon, "minElevationAngle", 0);
%access setting
ac1 = access(gs(1), sat);
ac2 = access(gs(2), sat);
% ac = access(gs(2), sat(1));
%interval status
intvls = accessIntervals(ac1);
disp(intvls);
%display FOV of camera
fieldOfView(cam);
time = datetime(2020,6,02,12,30,0);
[az,elev,r] = aer(gs(1),sat(1),time);