编了个程序,得到了预期的多条x-y曲线,但是想得到各时间点的具体值。研究了好几天,还是做不出。哪位高手能演示下?随便输出哪一条曲线都可以。程序如下:
function IV
tspan=[0, 30];
x0 = [0 0 0 0 0 100 0 0 0 0 0 0 0 0 0 0 100];
[T, X] = ode45(@lun_fun, tspan, x0);
plot(T, X)
function dx = lun_fun(t, x)
dx = zeros(17, 1);
k12=0.69;
k13=4.6;
k21=0.35;
k26=0.33;
k27=0.67;
k34=0.66;
k30=6.6;
k43=0.33;
k46=0.17;
k47=0.33;
k56=0.01;
k62=0.67;
k64=0.33;
k65=0.005;
k68=0.17;
k69=1.7;
k610=10;
k611=1.7;
k612=2.2;
k76=0.023;
k86=0.083;
k96=0.023;
k100=0.0046;
k103=0.023;
k106=0.023;
k110=0.023;
k116=0.83;
k120=0.0046;
k126=1.1;
dx(1)=k21*x(2)-k13*x(1)-k12*x(1);
dx(2)=k12*x(1)+k62*x(6)-k26*x(2)-k27*x(2)-k21*x(2);
dx(3)=k13*x(1)-k34*x(3)-k30*x(3)+k103*x(10)+k43*x(4);
dx(4)=k34*x(3)+k64*x(6)-k47*x(4)-k46*x(4)-k43*x(4);
dx(5)=k65*x(6)-k56*x(5);
dx(6)=k26*x(2)+k46*x(4)+k76*x(7)+k56*x(5)+k86*x(8)+k96*x(9)+k106*x(10)+k116*x(11)+k126*x(12)-k62*x(6)-k64*x(6)-k65*x(6)-k68*x(6)-k69*x(6)-k610*x(6)-k611*x(6)-k612*x(6);
dx(7)=k47*x(4)+k27*x(2)-k76*x(7);
dx(8)=k68*x(6)-k86*x(8);
dx(9)=k69*x(6)-k96*x(9);
dx(10)=k610*x(6)-k100*x(10)-k103*x(10)-k106*x(10);
dx(11)=k611*x(6)-k110*x(11)-k116*x(11);
dx(12)=k612*x(6)-k120*x(12)-k126*x(12);
dx(13)=k30*x(3);
dx(14)=k100*x(10);
dx(15)=k110*x(11);
dx(16)=k120*x(12);
dx(17)=-k13*x(1)+k34*x(3)+k30*x(3)-k103*x(10)-k43*x(4)-k21*x(2)+k13*x(1)+k12*x(1)-k30*x(3)-k100*x(10)-k110*x(11)-k120*x(12);