Я использую ode45 для решения системы уравнений. Моя программа работает отлично. Поскольку размер шага, который решает решатель, составляет автоматический c, я хочу получить доступ к вектору решения на каждом шаге времени моделирования. Как это сделать? Мой код:
%% Solve
opts = odeset('RelTol',1e-3);
[T,X]=ode45(@(T,X)sys(T,X),Tspan,X(:,1),opts);
%% defining the function
function Xdot = sys(T,X)
Xdot= A*X + B*U
end
Мне нужен доступ к X-вектору на каждом временном шаге.