二元常系数线性微分方程matlab算法

如题所述

求解常微分方程可用ode*系列函数(典型如ode45)。

参考代码:
dx = @(t,x) [ ...
0.097*x(1) - 0.1289*x(2) + 465.2078; ...
0.2477*x(1) - 0.3503*x(2) + 61.6214 ...
];
x0 = [481.8 166.2];
t = 1 : 11;
[t, y] = ode45(dx, t, x0);
disp(num2str([t y]))

输出(三列依次为t、x1、x2):
1 481.8 166.2
2 987.09858 325.12652
3 1518.0848 545.61211
4 2069.9726 814.28125
5 2639.0849 1120.8389
6 3222.5984 1457.3624
7 3818.348 1817.7561
8 4424.675 2197.3316
9 5040.3111 2592.4845
10 5664.2885 3000.445
11 6295.8714 3419.0863
温馨提示:答案为网友推荐,仅供参考