大家好，又见面了，我是你们的朋友全栈君。
matlab解矩阵微分方程
MATLAB
2018 06 25
—
MATLAB
1
Definition1 t
′ ( )
Definition2
∫ ( )
∫
2
4 MATLAB
ode45
1: —
1
1–MATLAB
⩾
(1)
2.1 4
MATLAB
1 % RUNGEKUTTA4 4
2 % hyowinner(MATALB )
3 %
4 % ( )
5 % 2018 624
6 % http ://www. ilovematlab . cn/ thread 4280811.html
7 % 2.0
8
9 function [ x, y ] =rungeKutta4 ( ufunc , a , b, y0 , h )
10 % ufunc a
11 % b y0 h ode45
12 n =floor ((b a) / h);%
13 x =ones(n ,1);
14 x(1) =a;%
15 y(: ,1) =y0;%
16 for ii =1:n
17 x( ii + 1) =x( ii ) + h;
18 k1 =ufunc (x( ii ) , y(: , ii ));
19 k2 =ufunc (x( ii ) + h / 2, y(: , ii ) + h * k1 / 2);
20 k3 =ufunc (x( ii ) + h / 2, y(: , ii ) + h * k2 / 2);
21 k4 =ufunc (x( ii ) + h, y(: , ii ) + h * k3 );
22 y(: , ii +1) =y(: , ii ) + h * (k1 +
发布者：全栈程序员-站长，转载请注明出处：https://javaforall.net/127324.html原文链接：https://javaforall.net
