有关热力学的吗?还是别的?
哦!解析法是不行,那就用数值方法吧。
function odemethod
[T,X]=ode23t(@fun,0:0。1:10,[0 0 0 0])
plot(T,X)
legend('x','y','z','t'),xlabel('T')
function dX=fun(t,X)
x=X(1);y=X(2);z=X(3);w=X(4);
dX1=-0。
0558*x-0。9968*y 0。0802*z 0。0415*w 0。00729;
dX2=0。598*x-0。115*y-0。0318*z-0。475 0。00775;
dX3=-0。05*x 0。
388*y-0。465*z 0。296;
dX4=0。805*x z;
dX=[dX1;dX2;dX3;dX4];
运行结果:
T =
0
0。
1000
0。2000
0。3000
0。4000
0。5000
0。6000
0。
7000
0。8000
0。9000
1。0000
1。1000
1。2000
1。
3000
1。4000
1。5000
1。6000
1。7000
1。8000
1。
9000
2。0000
2。1000
2。2000
2。3000
2。4000
2。
5000
2。6000
2。7000
2。8000
2。9000
3。0000
3。
1000
3。2000
3。3000
3。4000
3。5000
3。6000
3。