orbit formula graphing orbits

15 次查看（过去 30 天）
spaceman45
spaceman45 2022-10-20
评论： Steven Lord 2022-10-22
%Orbit One
thetadeg = 0:1:360; h = 51593.15284; mu = 398600; e = 0; r = (h^2/mu)./(1-e.*cosd(thetadeg)); polarplot(thetadeg,r)
%Orbit Two
thetadeg = 0:1:360; h = 56935.55284; mu = 398600; e = 0.217819; r = (h^2/mu)./(1-e.*cosd(thetadeg)); polarplot(thetadeg,r)
%Orbit Three
thetadeg = 0:1:360; h = 62277.95284; mu = 398600; e = 0.45708; r = (h^2/mu)./(1-e.*cosd(thetadeg)); polarplot(thetadeg,r)
%Orbit Four
thetadeg = 0:1:360; h = 67620.35284; mu = 398600; e = 0.7177924; r = (h^2/mu)./(1-e.*cosd(thetadeg)); polarplot(thetadeg,r)
%Escape thetadeg = 0:1:120; h = 72962.75284; mu = 398600; e = 1.00; r = (h^2/mu)./(1-e.*cosd(thetadeg)); polarplot(thetadeg,r)
  1 个评论
Steven Lord
Steven Lord 2022-10-22
Rather than flagging your own question as Unclear, you should edit the question or add a comment with more information clarifying what you've asked. I've removed the flag.

请先登录，再进行评论。

请先登录，再回答此问题。

回答（1 个）

David Hill
David Hill 2022-10-21
Ran in:
%Orbit One
thetadeg = 0:.1:360;
h = 51593.15284;
mu = 398600;
e = 0;
r = (h^2/mu)./(1-e.*cosd(thetadeg));
polarplot(deg2rad(thetadeg),r)
hold on; %need hold on
%Orbit Two
h = 56935.55284;
e = 0.217819;
r = (h^2/mu)./(1-e.*cosd(thetadeg));
polarplot(deg2rad(thetadeg),r)
%Orbit Three
h = 62277.95284;
e = 0.45708;
r = (h^2/mu)./(1-e.*cosd(thetadeg));
polarplot(deg2rad(thetadeg),r)
%Orbit Four
h = 67620.35284;
e = 0.7177924;
r = (h^2/mu)./(1-e.*cosd(thetadeg));
polarplot(deg2rad(thetadeg),r)
%Escape need to reduce the range so it fits on the plot
%
h = 72962.75284;
e = 1.00;
r = (h^2/mu)./(1-e.*cosd(thetadeg));
Rad=deg2rad(thetadeg);
polarplot(Rad(484:3118),r(484:3118))

类别

在 Help Center 和 File Exchange 中查找有关 Satellite and Orbital Mechanics 的更多信息

标签

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by