Syntax: dydt=getODEderivatives(t,y,odehandler)
This function is used to get the time course derivatives
evaluated to solve an ode system using Matlab solvers
(ode45,ode23,ode23s and so on).
Example:
Want to solve the system and plot the solution as well as how the derivatives change:
dudt=w;
dwdt=-u;
u(0)=0;
w(0)=1;
odehandler=@(t,y) [y(2);-y(1)];
[t,y]=ode45(odehandler,[0 4*pi],[0 1]);
u=y(:,1);
w=y(:,2);
(select the code below and press F9 to run it).
dydt=getODEderivatives(t,y,odehandler);
dudt=dydt(:,1);
dwdt=dydt(:,2);
%Check that dudt=w;
sum(dudt==w)==length(t)
plot(t,dudt,'r--',t,w,'g-.')
legend('dudt','w')
%Check that dwdt=-u;
sum(dwdt==-u)==length(t)
figure
plot(t,dwdt,'r--',t,-u,'g-.')
legend('dwdt','-u')
引用格式
Luca Cerone (2024). getODEderivatives() (https://www.mathworks.com/matlabcentral/fileexchange/27469-getodederivatives), MATLAB Central File Exchange. 检索来源 .
MATLAB 版本兼容性
平台兼容性
Windows macOS Linux类别
- MATLAB > Mathematics > Numerical Integration and Differential Equations > Ordinary Differential Equations >
标签
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!版本 | 已发布 | 发行说明 | |
---|---|---|---|
1.0.0.0 |