因为你要建立一个圆形轨道，只需按周期缩放时间就可以得到θ。例如，因为一个周期的角度是2pi，
θ=2*pi*t/t；
那就用这个和rn来画你的图。
对一个简单的圆形卫星轨道进行时间建模
134个视图（过去30天）
显示旧的评论
因此，我基本上是试图在极坐标图上画出一个轨道，它将是一个圆，它代表一颗500kg卫星在没有外力作用的情况下的运动。我忽略了一个事实，它现在应该随着时间的推移而下降。然而，我的代码没有产生预期的结果，所以任何输入都会被估价，谢谢大家
clc；
清楚的全部的;
G=6.673e-11；%引力常数
M=5.98e24；%地球质量（单位：kg）
ra=100000；%轨道距离（米）
r=6.37e6+ra；%轨道总半径（m）
m=500；%质量卫星（kg）
a=（G*M）/（r^2）；%检查加速度
v_orb=sqrt（（G*M）/r）；%轨道速度（m/s）
T=sqrt（（（4*（pi^2））*r^3/（G*M））；%期间（s）
%让我们以150步为单位绘制周期图，以减少计算量
simt=T；
对于t=1:150:simt
v_o（t）=sqrt（（G*M）/r）；%与时间对应的速度数组，步长为150。。不应更改为1x34-但有错误
rn（t）=6.37e6+ra；%与时间对应的半径数组，步长为150。。不应更改为1x34-但有错误
分区（t）%以150的步长从1变为5101
%%%创建与t？？？%%%对应的1x34形式的数组
结束
图形
图（t，v_o）；xlabel(“模拟时间”,“字体大小”, 12);
伊拉贝尔(“速度”,“字体大小”, 12);
网格；%%%应为直线%%%的图形
图形
图（t，rn）；xlabel(“模拟时间”,“字体大小”, 12);
伊拉贝尔(“半径”,“字体大小”, 12);
网格；%%%应为直线%%%的图形
%%%将半径和时间点转换为极坐标并绘制%%%
%%%当卫星在周期%%的一段时间内围绕恒定半径的轨道运行时，预计会出现一个圆
%[theta，rho]=cart2pol（t，rn）；
θ=atan2（rn，t）；
rho=sqrt（（t.^2）+（rn.^2））；
θ=θ*（180/pi）；%到一定程度
图形
polarplot（θ，ρ）
头衔(“轨道”)
0评论
更多答案（2）
米哈尔·斯莱辛斯基
2020年2月3日
这是我的最后一段代码，它现在似乎可以工作了：
clc；
清楚的全部的;
G=6.673e-11；%引力常数
M=5.98e24；%地球质量（单位：kg）
ra=100000；%轨道距离（米）
r=6.37e6+ra；%轨道总半径（m）
m=500；%质量卫星（kg）
a=（G*M）/（r^2）；%检查加速度
v_orb=sqrt（（G*M）/r）；%轨道速度（m/s）
T=sqrt（（（4*（pi^2））*r^3/（G*M））；%期间（s）
%让我们以150步为单位绘制周期图，以减少计算量
步骤=T/35；
simt=-步；
对于i=1:1:36
simt=simt+步数；
t（i）=simt；
v_o（i）=sqrt（G*M）/r；%与时间对应的速度数组，步长为150。。不应更改为1x34-但有错误
rn（i）=6.37e6+ra；%与时间对应的半径数组，步长为150。。不应更改为1x34-但有错误
%disp（t2）%150步从1到5101
结束
图形
图（t，v_o）；xlabel(“模拟时间”,“字体大小”, 12);
伊拉贝尔(“速度”,“字体大小”, 12);
网格；%%%应为直线%%%的图形
图形
图（t，rn）；xlabel(“模拟时间”,“字体大小”, 12);
伊拉贝尔(“半径”,“字体大小”, 12);
网格；%%%应为直线%%%的图形
%%%将半径和时间点转换为极坐标并绘制%%%
%%%当卫星在周期%%的一段时间内围绕恒定半径的轨道运行时，预计会出现一个圆
%[theta，rho]=cart2pol（t，rn）；
%th=atan2（rn，t）；
th=2*pi*t/t；
rho=sqrt（（t.^2）+（rn.^2））；
%disp（th）；
图形
极性（th，rho）
头衔(“轨道”)
0评论
您还可以从以下列表中选择网站：