5 (@ fck, [0 1500], [1500 0 0 0]);(t,c,«LineWidth»,2);=legend(«c(1)»,«c(2)»,«c(3)»,«c(4)»,40);(h,«Interpreter»,«none»);;
function dc=fck (t, c)
K1=9.6207e +010;=6.8719e +010;=9.6207e +010;=2.0616e +011;=0; _R=9232.7;
w=25;=300 + w * sin (2 * pi * t/10); =zeros(4,1);(1)=K1*exp(-E_R/T1)*c(2)-K4*exp(-E_R/T1)*c(1);(2)=K2*exp(-E_R/T1)*c(3)-K1*exp(-E_R/T1)*c(2)(3)=K4*exp(-E_R/T1)*c(1)-K3*exp(-E_R/T1)*c(3)-K2*exp(-E_R/T1)*c(3);(4)=K3*exp(-E_R/T1)*c(3);
Лістінг3ky3
[T, C]=ode15s (@ gif, [0:1000], [1500 0 0 0300]); (T,C,«LineWidth»,2);=legend(«c(1)»,«c(2)»,«c(3)»,«c(4)»,«T»);(h,«Interpreter»,«none»);on=C(end,1)=C(end,2)=C(end,3)=C(end,4)=C(end,5)dc=gif(t,c)=9232.7;=9.6207e+010;=6.8719e+010;=9.6207e+010;=2.0616e+011;=0;=0.5/60; =1.1; =2000 * 1300; =12000;=10000;=10000;=20000;=0;=200; =7; =1400 * 1000; =V=275; =zeros(5,1);(1)=(v/V)*(1500-c(1))+K1*c(2)*exp(-Er/c(5))-K4*c(1)*exp(-Er/c(5));(2)=(v/V)*(0-c(2))+K2*c(3)*exp(-Er/c(5))-K1*c(2)*exp(-Er/c(5));(3)=(v/V)*(0-c(3))+K4*c(1)*exp(-Er/c(5))-K3*c(3)*exp(-Er/c(5))-K2*c(3)*exp(-Er/c(5));(4)=(v/V)*(0-c(4))+K3*c(3)*exp(-Er/c(5));(5)=(v/V)*(300-c(5))+(1/Cpcm)*((H1)*K1*c(2)*exp(-Er/c(5))+(H2)*K2*exp(-Er/c(5))*c(3)+(H3)*K3*c(3)*exp(-Er/c(5))+(H4)*K4*exp(-Er/c(5))*c(1))+(Kt*F/(V*Cpcm))*((F*Kt*c(5)+2*Cpx*vx*Txo)/(F*Kt+2*Cpx*vx)-c(5));
Лістинг 4ky4
[T, C]=ode15s (@ gif, [0:2000], [252.5062 144.2706 266.2115 837.0118 318.0719]); (T, C (:, 5)) on
T0=T, C (:, 5) dc=gif (t, c)=9232.7;=9.6207e +010;=6.8719e +010;=9.6207e +010;=2.0616e +011;=0 ;=0.5/60; =1.1; =2000 * 1300; =12000;=10000;=10000;=20000;=0;
Kt=200; =7; =1400 * 1000; =V * 1.05; =275; =Zeros (5,1);
dc(1)=(v/V)*(1500-c(1))+K1*c(2)*exp(-Er/c(5))-K4*c(1)*exp(-Er/c(5));(2)=(v/V)*(0-c(2))+K2*c(3)*exp(-Er/c(5))-K1*c(2)*exp(-Er/c(5));(3)=(v/V)*(0-c(3))+K4*c(1)*exp(-Er/c(5))-K3*c(3)*exp(-Er/c(5))-K2*c(3)*exp(-Er/c(5));(4)=(v/V)*(0-c(4))+K3*c(3)*exp(-Er/c(5));(5)=(v/V)*(300-c(5))+(1/Cpcm)*((H1)*K1*c(2)*exp(-Er/c(5))+(H2)*K2*exp(-Er/c(5))*c(3)+(H3)*K3*c(3)*exp(-Er/c(5))+(H4)*K4*exp(-Er/c(5))*c(1))+(Kt*F/(V*Cpcm))*((F*Kt*c(5)+2*Cpx*vx*Txo)/(F*Kt+2*Cpx*vx)-c(5));
Лістінг5ky55
[T, C]=ode15s (@ gif, [0:2000], [252.5062 144.2706 266.2115 837.0118 318.0719 100]); (T, C (:, 5)) on dc=gif (t, c)=9232.7;=9.6207e +010;=6.8719e +010;=9.6207e +010;=2.0616e +011;=0;= 0.5/60; =1.1; =2000 * 1300; =12000;=10000;=10000;=20000;=0;=200; =7; =1400 * 1000; =V;=275; =7.63;=104.44;=320;=zeros(5,1);(1)=(v/V)*(1500-c(1))+K1*c(2)*exp(-Er/c(5))-K4*c(1)*exp(-Er/c(5));(2)=(v/V)*(0-c(2))+K2*c(3)*exp(-Er/c(5))-K1*c(2)*exp(-Er/c(5));(3)=(v/V)*(0-c(3))+K4*c(1)*exp(-Er/c(5))-K3*c(3)*exp(-Er/c(5))-K2*c(3)*exp(-Er/c(5));(4)=(v/V)*(0-c(4))+K3*c(3)*exp(-Er/c(5));(5)=(v/V)*(300-c(5))+(1/Cpcm)*((H1)*K1*c(2)*exp(-Er/c(5))+(H2)*K2*exp(-Er/c(5))*c(3)+(H3)*K3*c(3)*exp(-Er/c(5))+(H4)*K4*exp(-Er/c(5))*c(1))+(Kt*F/(V*Cpcm))*((F*Kt*c(5)+2*Cpx*(vx*c(6)/100)*Txo)/(F*Kt+2*Cpx*(vx*c(6)/100))-c(5));(6)=-Kp*((v/V)*(300-c(5))+(1/Cpcm)*((H1)*K1*c(2)*exp(-Er/c(5))+(H2)*K2*exp(-Er/c(5))*c(3)+(H3)*K3*c(3)*exp(-Er/c(5))+(H4)*K4*exp(-Er/c(5))*c(1))+(Kt*F/(V*Cpcm))*((F*Kt*c(5)+2*Cpx*(vx*c(6)/100)*Txo)/(F*Kt+2*Cpx*(...
