Is it possible to get a time-dependent function as an output of an ODE solver?

Question

Nang Su Lwin 2022-11-23

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1861058-is-it-possible-to-get-a-time-dependent-function-as-an-output-of-an-ode-solver

回答： Steven Lord 2022-11-24

My code is below, it plots as I expect it to. I would like to know if it is possible to get an output where the ODE solver gives the solution as a function in time rather than 2 column vectors [T,Y]. I wanted to use the output to further analyze something, but need a function in time rather than the column vectors.

clear all

clc

I0= 500; % maximum flow

Tc=60/72; % heart period

Ts=(2/5*Tc); % time in systole

P_ss=80; % diastolic pressure

R= 1;

C=1;

R1=0.05;

L=0.005;

I=@(t)I0*sin((pi*t)/Ts).^2.*(t<=Ts); %input current flow

Idot = @(t)I0*2*sin(pi*t/Ts).*cos(pi*t/Ts)*pi/Ts.*(t<=Ts);

Idotdot = @(t) I0*2*(cos(pi*t/Ts).^2 - sin(pi*t/Ts).^2)*(pi/Ts)^2.*(t<=Ts);

fun = @(t,y)[y(2);(Idotdot(t)*(R*L*C*R1)+Idot(t)*(L*(R+R1))+I(t)*(R*R1) - (y(2)*(C*R*R1+L)+y(1)*R1))/(L*C*R)];

y0 = [80 ;0];

tspan = [0 60/72];

[T,Y] = ode45(fun,tspan,y0);

plot(T,Y(:,1),'g')

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

请先登录，再进行评论。

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

Answer 1

Torsten 2022-11-23

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1861058-is-it-possible-to-get-a-time-dependent-function-as-an-output-of-an-ode-solver#answer_1110028

在 MATLAB Online 中打开