怎样用matlab求解复杂的微分方程组vx'=-(vx+1.4)*((vx+1.4)^2+vy^2+(vz-1.4)^2)^0.5; vy'=-vy*((vx+1.4)^2+vy^2+(vz-1.4)^2)^0.5; vz'=-(vz-1.4)*((vx+1.4)^2+vy^2+(vz-1.4)^2)^0.5; 其中： vx(0)=1500; vy(0)=-500; vz(0)=0;

怎样用matlab求解复杂的微分方程组vx'=-(vx+1.4)*((vx+1.4)^2+vy^2+(vz-1.4)^2)^0.5; vy'=-vy*((vx+1.4)^2+vy^2+(vz-1.4)^2)^0.5; vz'=-(vz-1.4)*((vx+1.4)^2+vy^2+(vz-1.4)^2)^0.5; 其中： vx(0)=1500; vy(0)=-500; vz(0)=0;
怎样用matlab求解复杂的微分方程组
vx'=-(vx+1.4)*((vx+1.4)^2+vy^2+(vz-1.4)^2)^0.5;
vy'=-vy*((vx+1.4)^2+vy^2+(vz-1.4)^2)^0.5;
vz'=-(vz-1.4)*((vx+1.4)^2+vy^2+(vz-1.4)^2)^0.5;
其中：
vx(0)=1500;
vy(0)=-500;
vz(0)=0;

怎样用matlab求解复杂的微分方程组vx'=-(vx+1.4)*((vx+1.4)^2+vy^2+(vz-1.4)^2)^0.5; vy'=-vy*((vx+1.4)^2+vy^2+(vz-1.4)^2)^0.5; vz'=-(vz-1.4)*((vx+1.4)^2+vy^2+(vz-1.4)^2)^0.5; 其中： vx(0)=1500; vy(0)=-500; vz(0)=0;
用数值解法吧.matlab程序为：
function hhh
clear;clc;
[TOUT,X] = ode23t(@odefun,[0,0.01],[1500 -500 0])
%t=TOUT,vx=X(:,1),vy=X(:,2),vz=X(:,3)
plotyy(TOUT,X(:,1:2),TOUT,X(:,3))
function yhat=odefun(t,x)
vx=x(1);vy=x(2);vz=x(3);
dx=-(vx+1.4)*((vx+1.4)^2+vy^2+(vz-1.4)^2)^0.5;
dy=-vy*((vx+1.4)^2+vy^2+(vz-1.4)^2)^0.5;
dz=-(vz-1.4)*((vx+1.4)^2+vy^2+(vz-1.4)^2)^0.5;
yhat=[dx;dy;dz];
结果：
TOUT =
0
0.0000
0.0000
0.0000
0.0000
0.0000
0.0000
0.0000
0.0001
0.0001
0.0001
0.0002
0.0002
0.0002
0.0003
0.0003
0.0004
0.0005
0.0006
0.0006
0.0007
0.0008
0.0009
0.0010
0.0012
0.0013
0.0014
0.0016
0.0017
0.0019
0.0021
0.0023
0.0025
0.0028
0.0030
0.0033
0.0035
0.0038
0.0042
0.0045
0.0049
0.0054
0.0058
0.0063
0.0068
0.0073
0.0078
0.0085
0.0091
0.0098
0.0100
X =
1.0e+003 *
1.5000 -0.5000 0
1.4991 -0.4997 0.0000
1.4973 -0.4991 0.0000
1.4947 -0.4982 0.0000
1.4907 -0.4969 0.0000
1.4710 -0.4903 0.0000
1.4397 -0.4799 0.0001
1.3976 -0.4659 0.0001
1.3499 -0.4500 0.0001
1.3054 -0.4352 0.0002
1.2500 -0.4167 0.0002
1.1990 -0.3998 0.0003
1.1393 -0.3799 0.0003
1.0853 -0.3619 0.0004
1.0248 -0.3418 0.0004
0.9708 -0.3238 0.0005
0.9110 -0.3039 0.0005
0.8582 -0.2863 0.0006
0.8004 -0.2670 0.0007
0.7500 -0.2502 0.0007
0.6954 -0.2321 0.0008
0.6482 -0.2163 0.0008
0.6070 -0.2026 0.0008
0.5629 -0.1879 0.0009
0.5248 -0.1752 0.0009
0.4915 -0.1641 0.0009
0.4557 -0.1522 0.0010
0.4249 -0.1420 0.0010
0.3979 -0.1330 0.0010
0.3690 -0.1234 0.0011
0.3441 -0.1151 0.0011
0.3223 -0.1078 0.0011
0.2989 -0.1000 0.0011
0.2787 -0.0933 0.0011
0.2611 -0.0874 0.0012
0.2422 -0.0811 0.0012
0.2259 -0.0757 0.0012
0.2116 -0.0709 0.0012
0.1963 -0.0658 0.0012
0.1831 -0.0614 0.0012
0.1692 -0.0568 0.0012
0.1572 -0.0528 0.0013
0.1468 -0.0494 0.0013
0.1359 -0.0457 0.0013
0.1265 -0.0426 0.0013
0.1183 -0.0399 0.0013
0.1111 -0.0375 0.0013
0.1029 -0.0347 0.0013
0.0959 -0.0324 0.0013
0.0897 -0.0303 0.0013
0.0880 -0.0298 0.0013
积分区间自己定,这里暂定为[0,0.01].