How to check the complementary time-dependent term of differential equations.
1 visualizzazione (ultimi 30 giorni)
Mostra commenti meno recenti
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.
0 Commenti
Risposte (1)
SALAH ALRABEEI
il 7 Giu 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')
Vedere anche
Categorie
Scopri di più su Ordinary Differential Equations in Help Center e File Exchange
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!