Run a code for 100 years

Question

Adam Kelly 2020-4-11

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/517032-run-a-code-for-100-years

回答： Les Beckham 2020-4-11

I have a program that I need to run for 100 years, the program is calculating the angle of a circular orbit at 451 days. I want to calculate the angle changing over 100 years, I have the angle calculated for 451 days just want to repeat the same process but add 451 to the last calculation. T=451 days my idea was to have a for loop for the 100 years starting at T=451 days for the first year and then adding 451 to T for each year till 100.

clear all;
clc;
ro=1.49e11;             %Sun orbital radius (1 AU)            
G=6.67408e-11;
MS=2e+30;
%synodic period
perp = 365;             %planet period (days)
perd = 1896.59;         %red dwarf period (days)
np = 2*pi / perp;       % mean motion of planet (rad/day)
nd = 2*pi / perd;       % mean motion of red dwarf (rad/day)
T = 2*pi / (np - nd)    %This is the average time between stellar eclipses.
v=sqrt((G*MS)/ro);
theta=((v*T)*86400)/ro  %Angle in radians 
     

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

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

Answer 1

Les Beckham 2020-4-11

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/517032-run-a-code-for-100-years#answer_425305

在 MATLAB Online 中打开

I don't think the results are going to be very interesting but you could do what you are asking as follows:

ro=1.49e11;             %Sun orbital radius (1 AU)            
G=6.67408e-11;
MS=2e+30;
%synodic period
perp = 365;             % planet period (days)
perd = 1896.59;         % red dwarf period (days)
np = 2*pi / perp;       % mean motion of planet (rad/day)
nd = 2*pi / perd;       % mean motion of red dwarf (rad/day)
dT = 2*pi / (np - nd);  % This is the average time between stellar eclipses.
v=sqrt((G*MS)/ro);
T = dT:dT:dT*ceil(100*perp/dT); % up to next multiple of dT after 100 years (multiples of perp)
theta=((v*T)*86400)/ro  % Angle in radians at each multiple of dT

This is just going to add a constant increment to theta every dT. If you plot(T, theta) you will get a boring straight line. You could, of course, apply a mod by 2*pi but this won't be much more useful (IMHO).

Anyway, I hope this answers your question. If not, please clarify.