How can I plot the orbit of Earth over the course of a year using Julian dates?

Question

Rebekah Ericson 2020-11-13

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/647078-how-can-i-plot-the-orbit-of-earth-over-the-course-of-a-year-using-julian-dates

编辑： Cris LaPierre 2020-11-13

The error I get:

Unrecognized function or variable 't_0'.
Error in MHCP_parametricEarth (line 2)
t = linspace (t_0,(1/366),t_0+1);

What I am attempting to do:

Graph the orbit of Earth over 1 year so that its path over a day shows up every 0.01 seconds.

My code:

t_o = juliandate(2020,10,13);
t = linspace (t_0,(1/366),t_0+1);
x=150*cos(2*pi(t-t_o));
y=150*sin(2*pi(t-t_o));
tcounter=0;
while tcounter<=1
    plot([x(tcounter) x(tcounter + 1/365)], [y(tcounter) y(tcounter + 1/365)])
    tcounter = tcounter + 2/365
pause(0.01)
end

1 个评论
显示 -1更早的评论隐藏 -1更早的评论

Rebekah Ericson 2020-11-13

Thank you for the answer. I have changed that and edited my code, but it is still not plotting right. It will run, but it's not giving a graph or anything. It will just store numbers in the workspace, and t, x, and y are all 1 by 0 vectors here.

jd = juliandate(2020,10,13)

t_0 = jd

jd_2 = juliandate(2021,10,13)

t_1 = jd_2

t = linspace (jd, jd_2, 1/366)

x=150*cos(2*pi*(t-t_0));

y=150*sin(2*pi*(t-t_0));

tcounter=t_0;

while tcounter<=1

plot([x(tcounter) x(tcounter+1/365)], [y(tcounter) y(tcounter + 1/365)])

hold on

tcounter = tcounter + 2/365