How to check the complementary time-dependent term of differential equations.

2 views (last 30 days)
CT on 7 Jun 2021
Commented: CT on 7 Jun 2021
Is there any way to check how the time-dependent term ft in the following code is completed by interp1?
function [T1,Y1] = ode_param_sample1()
TSPAN = [0 pi/2];
IC = 1;
tv = linspace(0,pi/2,10);
ftv = sin(tv);
[T1 Y1] = ode45(@myODE1, TSPAN, IC);
function dydt1 = myODE1(t1,y1)
ft = interp1(tv, ftv, t1,'spline');
dydt1 = - ft *y1;
end
end
I want to check the ft in this code.
By checking, I mean I want to compare it with the original time-variant parameters, such as plotting it on a graph.

SALAH ALRABEEI on 7 Jun 2021
% Try this
function [T1,Y1] = ode_param_sample1()
TSPAN = [0 pi/2];
IC = 1;
tv = linspace(0,pi/2,10);
ftv = sin(tv);
ft = interp1(tv, ftv, T1,'spline');
plot(tv,ftv,'or',T1,ft,'b')
CT on 7 Jun 2021
I tried, but I couldn't.
I get the message "The function or variable 'T1' is not recognized."

R2021a

Community Treasure Hunt

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

Start Hunting!