X' = AX + B形式的微分方程组,其中A,B和X是矩阵。我试图通过在Matlab中使用ODE45函数来解决这个问题。我很难选择ode45参数的时间跨度,因为问题没有指定。 注意:我试过时间= [0,10];但我最终得到了一个大的矩阵,其中有一个因子乘以它。
代码编写:
%myode45function.m file
function dx_dt = myode45function (t,x)
B=[1 2 3 4 5 6 7]';
A=[2 3 4 5 6 0 7; 3 6 2 1 3 5 4; 4 2 2 4 2 7 5; 5 1 4 3 5 2 1; 6 3 2 5 4 1 2; 0 5 7 2 1 8 0; 7 4 5 1 2 0 9];
dx_dt= [(A*x)+B];
end
%main code
clear all;
clc
time=[0,4];
initial=[1,1,1,1,1,1,1];
[t,x]=ode45(@myode45function, time, initial);