Question

matlab ode45解微分方程

求状态转移矩阵，形式是W'=A*W，然后A是一个矩阵，12×12维，矩阵A是时间t的函数，W初值是一个12×12维的单位函数，积分范围t从0到3，请问这个要怎么用ode45求解，并且求出在末状态3s时的W，谢谢了！

Answer 1

ode45函数
求一阶常微分方程的初值问题
[t,y] = ode45(@(t,y)y-2*t/y,[0,4],1);
plot(t,y);
求解 y’ – y + 2*t / y且初值y(0) = 1的常微分方程初值问题，返回自变量和函数的若干个值.
若不写返回值，则会自动作出函数随自变量的变化图像.
ode45(@(t,y)y-2*t/y,[0,4],1);