Программа построения временных диаграмм. clc g r pi 180 r g 180 pi load m y.dat t m y ,1 x m y ,2 y m y ,3 z m y ,4 Vx m y ,5 Vy m y ,6 Vz m y ,7 clear m y s tmp size t s m s tmp 1 clear s tmp load m f.dat Fz m f ,2 Fs m f ,3 Fl m f ,4 Fa m f ,5 U20 m f ,6 clear m f load m s.dat xs m s ,2 ys m s ,3 zs m s ,4 clear m s load m par.dat p m par ,2 e m par ,3 a m par ,4 Om m par ,5 i m par ,6 omg m par ,7 T m par ,8 u m par ,9 clear m par p n 6952137. e n 0 a n 6952137. Om n0 28.1 g r i n 97.6 g r omg n 346.725 g r T n 5765 ws 2 pi 365.2422 24 3600 for j 1 s m, tmp j Om n0 ws t j end Om n tmp clear tmp map 1,1,1 colormap map plot t, p, y min t max t, p n p n, r title Фокальный параметр, grid on print -dwin pause plot t, p-p n, y title dp, grid on print -dwin pause plot t, e, y min t max t, e n e n, r title Эксцентриситет, grid on print -dwin pause plot t, e-e n, y title de, grid on print -dwin pause plot t, a, y min t max t, a n a n, r title Большая полуось орбиты, grid on print -dwin pause plot t, a-a n, y title da, grid on print -dwin pause plot t,Om r g, y- ,t,Om n r g, r title Долгота восходящего узла, grid on print -dwin pause plot t,Om r g-Om n r g, y title dOm, grid on print -dwin pause plot t, i r g, y min t max t, i n r g i n r g, r title Наклонение, grid on print -dwin pause plot t, i r g-i n r g, y title di, grid on print -dwin pause plot t,T, y min t max t , T n T n, r title Период, grid on print -dwin pause plot t,T-T n, y title dT , grid on print -dwin pause plot3 x, y,z, b axis min x max x min y max y min z max z set gca, box, on title Положение МКА hold on plt plot3 0,0,0 erasemode, xor, markersize ,24 dk ceil length y 2500 for k 1 dk length y set plt, xdata, x k, ydata, y k, zdata, z k drawnow end hold off pause plot t,Fz, y title Гравитация Земли, grid on print -dwin pause plot t,Fs, y title Гравитация Солнца и солнечное давление, grid on print -dwin pause plot t, Fl, y title Гравитация Луны, grid on print -dwin pause plot t,Fa, y title Сопротивление атмосферы, grid on print -dwin pause plot t,U20, y title Нецентральность гравитационного поля Земли, grid on print -dwin pause plot t,Fz Fs Fl Fa U20, y title Суммарное возмущающее ускорение, grid on print -dwin pause clear all clc g r pi 180 r g 180 pi p n 6952137. e n 0 a n 6952137. Om n0 28.1 g r i n 97.6 g r omg n 346.725 g r T n 5765 load u par.dat t u u par ,1 p u u par ,2 e u u par ,3 a u u par ,4 Om u u par ,5 i u u par ,6 omg u u par ,7 T u u par ,8 u u u par ,9 clear u par load u f.dat Fz u u f ,2 Fs u u f ,3 Fl u u f ,4 Fa u u f ,5 U20 u u f ,6 clear u f s tmp size t u s m u s tmp 1 clear s tmp ws 2 pi 365.2422 24 3600 for j 1 s m u, tmp j Om n0 ws t u j end Om n u tmp clear tmp plot t u, p u, y min t u max t u, p n p n, r title Фокальный параметр, grid on print -dwin pause plot t u, p u-p n, y title dp, grid on print -dwin pause plot t u, e u, y min t u max t u, e n e n, r title Эксцентриситет, grid on print -dwin pause plot t u, e u-e n, y title de, grid on print -dwin pause plot t u, a u, y min t u max t u, a n a n, r title Большая полуось орбиты, grid on print -dwin pause plot t u, a u-a n, y title da, grid on print -dwin pause plot t u,Om u r g, y- ,t u,Om n u r g, r title Долгота восходящего узла, grid on print -dwin pause plot t u,Om u r g-Om n u r g, y title dOm, grid on print -dwin pause plot t u, i u r g, y min t u max t u, i n r g i n r g, r title Наклонение, grid on print -dwin pause plot t u, i u r g-i n r g, y title di, grid on print -dwin pause plot t u,T u, y min t u max t u , T n T n, r title Период, grid on print -dwin pause plot t u,T u-T n, y title dT , grid on print -dwin pause clear all.