Разработка и анализ математических моделей с использованием MATLAB и MAPLE. Коробейников А.Г. - 136 стр.

UptoLike

Составители:

Рубрика:

 135

% (маркеры - точки, линия - сплошная)

plot(T,X(:,1),'r.-')

%Вывод графика производной от решения исходного дифференциального

уравнения

%маркеры - точки, линия - сплошная)

hold on

plot(T,X(:,2),'k.:')

% Вывод пояснений на график

xlabel('\itt')

legend('x', 'y',4)

grid on

hold off

%Подпрограмма вычисления правых частей уравнений

function F = osci1(t,x)

F=[x(3); x(4); -x(1)/(x(1)^2+x(2)^2)^(3/2); -x(2)/(x(1)^2+x(2)^2)^(3/2)];

0 2 4 6 8 10 12 14

-1

-0.8

-0.6

-0.4

-0.2

0.2

0.4

0.6

0.8

Рис.18.5. Движение спутника по эллиптической орбите (а).