matlab 二元微分方程组,求助，matlab求解二元二阶的常微分方程组

边值问题Matlab可用BVP4C命令，但感觉比较麻烦，下面用1stOpt求解，很简单快捷：CODE:

Constant Pey=9.73, Nox=8.05, uxuy=3, bd=1, cx1e=2.3,

b1=309.7,b2=2832.5, a1=27.8,a2=2.15,a3=-0.84,a4=0.935;

Variable t=[0:0.025:1],x=1,y'=[0,-Pey*bd*y];

Plot x,y,y';

ODEFunction x'=-(Nox/b1)*(b1*x+cx1e-(a1*y+a2)/(a3*y+a4));

y''=-Pey*bd*y'-(Nox*Pey*bd*uxuy/b2)*(b1*x+cx1e-(a1*y+a2)/(a3*y+a4));

边值估算:

y(t=0): 0.327786532200411

算法: 龙格-库塔-费尔博格法(Runge-Kutta-Fehlberg Method)

步长值: 0.025

步长数: 40

种群数: 5

结果:

t             x(t)                         y'(t)                     y(t)                         x'(t)                        y''(t)                  y'(t)

0             1                            0                            0.327786532200411            -7.66600291732416            -24.4666416985093          0

1             0.000693840858838046        -0.00586485658562458        0.000602760183379987        -0.00510084725252375        0.0407853057400083        -0.00586485658562458

结果过程:

No.        t             x(t)                         y'(t)                     y(t)                         x'(t)                        y''(t)                  y'(t)

0        0             1                            0                            0.327786532200411            -7.66600291732416            -24.4666416985093          0

1        0.025        0.826313919267495            -0.490075360338193           0.321182502343127            -6.27870464940025            -15.2705386363518          -0.490075360338193

2        0.05          0.683848068843495