Is it possible to use ode45() once instead of many times for different distinct times?

1 次查看（过去 30 天）
Majid
Majid 2021-11-12
评论： Majid 2021-11-13
Hi everyone,
I need to solve a system of differential equations and find the value of solutions in certain integer-valued times.
The system is
F(x,y,x', y')=0
G(x,y,x', y')=0
With the initial conditions
x(0) = p1, y(0) = p1
x'(0) = v1, y'(0) = v2.
I firstly used the dsolve() function but it could not find the explicit solution, and so I used the ode45() function then. Here is my code to do it.
function [Sol, Nsol] = SolvingSysEq2D(oDe, V, P, tt)
%oDe is the vector of equations, V is the velocity vector, P is the initial point, tt is the time
syms T Y t x(t) y(t) z(t) v1 v2
EqnS = [oDe(1) == 0, oDe(2) == 0];
[VF, ~] = odeToVectorField(EqnS);
ExEqSyst = matlabFunction(VF, 'Vars',{T, Y});
n = size(V,1);
Sol = zeros(1,2);
for i = 1:n
if rem(i,1)==0, disp(i), end
conD = [P(2), V(i,2), P(1), V(i,1)];
[~, Nsol] = ode45(@(T, Y)ExEqSyst(T, Y), [0, tt], conD);
Sol(i,:)= Nsol(end, [3, 1]);
end
I need to find the values of x(t) and y(t) at t = 1, 2, 3, 4, 5, …, 10. For now, I use my function ten times, each time for a t. I would like to know if there is a way that I use my function only one time and then obtain x(t) and y(t) for all the t =1, 2, …, 10.
I know that by using my function for t=10, I can draw the graphs of x(t) and y(t) on t in [0, 10], that is I can have all the data from t=0 to t=10, but I do not know how to get out the values of x(t) and y(t) at the specific times t = 1, 2, 3, 4, 5, …, 10.
Any help would be appreciated.

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

采纳的回答

Walter Roberson
Walter Roberson 2021-11-12
Instead of passing [0, tt] as the time span, pass in 0:10 . It will calculate for the intermediate times necessary to figure out all of 0 to 10, but it will only report at the times you list.
This applies any time the time vector has more than 2 elements. If you give it two elements exactly then it reports back at whatever times it feels like between those two.

更多回答（0 个）

类别

在 Help Center 和 File Exchange 中查找有关 Ordinary Differential Equations 的更多信息

标签

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by