在求解一组线性微分方程时,如何为Matlab选择ODE时间跨度?

时间:2015-03-21 04:11:48

标签: matlab linear ode differential-equations

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);

0 个答案:

没有答案