请教matlab解常微分方程组

mathtype · 发表于 2005-12-20 21:04

马上注册，结交更多好友，享用更多功能，让你轻松玩转社区。

您需要登录才可以下载或查看，没有账号？我要加入

x

为什么结果中会出现1.0e-006 *的字样，程序见最后，运行结果如下 Results by using ode45()： Columns 1 through 4 0 5.01000000000000 7.51500000000000 0 0.01000000000000 4.13485131470842 6.43984788572855 1.39799207811633 0.02000000000000 3.44758152320393 5.52114353135189 2.56308685158067 0.03000000000000 2.88227350501698 4.74679740666201 3.50341358451941 0.04000000000000 2.41653435229744 4.09216256729740 4.26205575708555 0.05000000000000 2.03217524208578 3.53707621966476 4.87385280064125 0.06000000000000 1.71435628646564 3.06495078243123 5.36707396630125 0.07000000000000 1.45097171346299 2.66212545408824 5.76464463014267 0.08000000000000 1.23227208332134 2.31748114244378 6.08490640206789 0.09000000000000 1.05019512891767 2.02170960733801 6.34292946270182 0.10000000000000 0.89816160057097 1.76711088582958 6.55091023207195 Columns 5 through 8 0 0 0 0 1.85417181101065 0.24242050252632 0.00000000127553 0.02755016531156 3.47326108268142 0.46693935334930 0.00000000134597 0.02361905620216 4.82434258441260 0.67032127023475 0.00000000141331 0.02027676642426 5.95468857583128 0.85454973825349 0.00000000147659 0.01742564578963 6.90285331290088 1.02140383872519 0.00000000153439 0.01498924068954 7.70042088738478 1.17248184606668 0.00000000158497 0.01290519252663 8.37328400050171 1.30922201148748 0.00000000162644 0.01112137975657 8.94241237043581 1.43290354260834 0.00000000165653 0.00959604442601 9.42527170416793 1.54469575551820 0.00000000167425 0.00828965815828 9.83622305409098 1.64565833968169 0.00000000167814 0.00717051296770 Columns 9 through 12 0 0 0 0 0.01487537005217 0.00000635950258 0.00000159951088 0.00004449791139 0.01188996848660 0.00000619685408 0.00000148163931 0.00004654861490 0.00957538000765 0.00000603525221 0.00000132989826 0.00004844465404 0.00775862261725 0.00000587728523 0.00000119685800 0.00005022986351 0.00632119422608 0.00000572427936 0.00000108178341 0.00005193644928 0.00517615801544 0.00000557772102 0.00000098246945 0.00005358686151 0.00425878093285 0.00000543879846 0.00000089690102 0.00005519692867 0.00352040597361 0.00000530984469 0.00000082333614 0.00005676360925 0.00292380589974 0.00000518995721 0.00000076018048 0.00005830197757 0.00243978006073 0.00000507966164 0.00000070602171 0.00005981314719 Columns 13 through 16 0 0 0 0 0.00001432058256 0.00000208776181 0.00007554196860 0.05361993867842 0.00001346707926 0.00000187341242 0.00006446725765 0.04743270307420 0.00001263867168 0.00000168826725 0.00005485176875 0.04167209530931 0.00001183812644 0.00000152627217 0.00004657371352 0.03637666942422 0.00001106451339 0.00000138282258 0.00003948062110 0.03153849678097 0.00001031728713 0.00000125453517 0.00003342976593 0.02715070521960 0.00000959597228 0.00000113888270 0.00002828763617 0.02320405525700 0.00000890289827 0.00000103414809 0.00002393339278 0.01968787828191 0.00000823565632 0.00000093871262 0.00002025565213 0.01658389853409 0.00000759495746 0.00000085143871 0.00001715646937 0.01387123157608 Columns 17 through 18 0 0 0.00000000000000 0.00000000000016 0.00000000000000 0.00000000000015 0.00000000000000 0.00000000000013 0.00000000000000 0.00000000000012 0.00000000000000 0.00000000000011 0.00000000000000 0.00000000000009 0.00000000000000 0.00000000000008 0.00000000000000 0.00000000000007 0.00000000000000 0.00000000000006 0.00000000000000 0.00000000000005 Cbanlance,Hbanlance,Obanlance: 1.0e-006 * -0.00000000088818 0 0 0.01352816670742 -0.02065940068556 0.02346317154819 -0.01075123884248 -0.10054292322770 -0.01225026438725 0.04524795755838 0.00057467034584 0.05518473977143 -0.08111397253933 -0.12845423767213 -0.05916933099101 0.00696812014700 0.12466885183850 -0.01876261282997 -0.00199972709785 0.07646786177685 -0.00829187888496 -0.00733501590228 -0.01100774752679 -0.00116685988822 0.00854392779299 0.03457823715796 0.05003402113966 -0.00880554723832 -0.02175283199375 -0.00626595662938 -0.00724876952105 -0.00868484242837 -0.00322983494963 >> 程序如下： function Simulate
clear all;clc
 Y0=[5.010 7.515 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0]; [t,Y]=ode23s(@f,[0:0.01:0.1],Y0); options.AbsTol=1e-30; % 绝对误差为1e-30; disp('Results by using ode45()：') disp([t Y]) plot(t,Y),xlabel('t'),ylabel('Y') legend('CH4','CO2','H2','CO','H2O','O2','CH3','CH2','CH','C','H','O','OH','HCO','CH2O','CH3O','HO2') Cbanlance=Y0(1)+Y0(2)-Y(:,1)-Y(:,2)-Y(:,4)-Y(:,7)-Y(:,8)-Y(:,9)-Y(:,10)-Y(:,14)-Y(:,15)-Y(:,16); Hbanlance=4*Y0(1)-4*Y(:,1)-2*Y(:,3)-2*Y(:,5)-3*Y(:,7)-2*Y(:,8)-Y(:,9)-Y(:,11)-Y(:,13)-Y(:,14)-2*Y(:,15)-3*Y(:,16)-Y(:,17); Obanlance=2*Y0(2)-2*Y(:,2)-Y(:,4)-Y(:,5)-2*Y(:,6)-Y(:,12)-Y(:,13)-Y(:,14)-Y(:,15)-Y(:,16)-2*Y(:,17); disp('Cbanlance,Hbanlance,Obanlance:') disp([Cbanlance Hbanlance Obanlance]) %-------------------------------------------------------------------------- function dYdt=f(t,Y)
M=12.525;NA=6.02217e+23; electrondensity=2e+11;Ne=electrondensity*1e+06/NA;
K01=3.234e+07;K02=1.960e+07;K03=1.251e+07;K04=3.788e+02;K05=2.673e+07; K1=3.137e+05;K2=2.984e+05;K3=2.165e+06;K4=5.443e+07;K5=8.000e+06;K6=8.000e+07;K7=3.337e+03; K8=2.000e+07;K9=8.131e+05;K10=4.000e+07;K11=1.475e+07;K12=1.119e+05;K13=9.479e+05;K14=2.500e+07; K15=2.000e+07;K16=1.000e+07;K17=4.000e+07;K18=1.688e+07;K19=4.490e+04;K20=5.697e+06; K21=3.175e+06;K22=6.896e+06;K23=1.401e+04;K24=6.009e+04;K25=5.000e+07;K26=2.000e+08; K27=3.000e+07;K28=3.000e+07;K29=3.000e+06;K30=1.118e+06;K31=9.358e+06;K32=2.324e+04; K33=1.000e+04;K34=6.846e+02;K35=1.881e+03;K36=9.359e+03;K37=3.824e+06;K38=1.000e+07; K39=2.219e+05;K40=1.746e+05;K41=1.958e+05;K43=1.407e+06;K44=2.894e+10;K45=1.741e+07; K46=2.000e+07;K47=1.500e+07;K48=1.965e+07;K49=1.965e+07;K50=1.581e+06;K51=1.066e+06;
f1=K4*Y(7)*Y(11)+K9*Y(7)*Y(17)+K23*Y(15)*Y(7)+K31*Y(14)*Y(7)-(K01+K02+K03+K04)*Y(1)*Ne-K1*Y(1)*Y(11)-K2*Y(1)*Y(12)-K3*Y(1)*Y(13); f2=K13*Y(8)*Y(6)+K28*Y(14)*Y(12)+K41*Y(4)*Y(13)-K05*Y(2)*Ne-K51*Y(2)*Y(9); f3=K02*Y(1)*Ne+K03*Y(1)*Ne+K1*Y(1)*Y(11)+K5*Y(7)*Y(13)+K10*Y(8)*Y(11)+K21*Y(15)*Y(11)+K26*Y(14)*Y(11)+K35*Y(11)*Y(11)*M+K39*Y(11)*Y(13)+K45*Y(17)*Y(11)-K40*Y(3)*Y(12)-K43*Y(13)*Y(3); f4=K05*Y(2)*Ne+K12*Y(8)*Y(6)+K15*Y(9)*Y(6)+K17*Y(9)*Y(12)+K24*Y(14)*M+K25*Y(14)*Y(13)+K26*Y(14)*Y(11)+K27*Y(14)*Y(12)+K29*Y(14)*Y(6)+K31*Y(14)*Y(7)+K48*Y(10)*Y(6)+K49*Y(10)*Y(13)+K50*Y(10)*Y(12)*M+K51*Y(2)*Y(9)-K41*Y(4)*Y(13); f5=K3*Y(1)*Y(13)+K22*Y(15)*Y(13)+K25*Y(14)*Y(13)+K32*Y(13)*Y(11)*M+K37*Y(13)*Y(13)+K43*Y(13)*Y(3)+K47*Y(17)*Y(13); f6=K9*Y(7)*Y(17)+K30*Y(14)*Y(17)+K34*Y(12)*Y(12)*M+K38*Y(12)*Y(13)+K45*Y(17)*Y(11)+K46*Y(17)*Y(12)+K47*Y(17)*Y(13)-K7*Y(7)*Y(6)-K11*Y(8)*Y(6)-K12*Y(8)*Y(6)-K13*Y(8)*Y(6)-K15*Y(9)*Y(6)-K16*Y(9)*Y(6)-K19*Y(16)*Y(6)-K29*Y(14)*Y(6)-K36*Y(11)*Y(6)*M-K48*Y(10)*Y(6); f7=K01*Y(1)*Ne+K1*Y(1)*Y(11)+K2*Y(1)*Y(12)+K3*Y(1)*Y(13)-K4*Y(7)*Y(11)-K5*Y(7)*Y(13)-K6*Y(7)*Y(12)-K7*Y(7)*Y(6)-K8*Y(7)*Y(17)-K9*Y(7)*Y(17)-K23*Y(15)*Y(7)-K31*Y(14)*Y(7); f8=K02*Y(1)*Ne-K10*Y(8)*Y(11)-K11*Y(8)*Y(6)-K12*Y(8)*Y(6)-K13*Y(8)*Y(6)-K14*Y(8)*Y(13); f9=K03*Y(1)*Ne+K10*Y(8)*Y(11)-K15*Y(9)*Y(6)-K16*Y(9)*Y(6)-K17*Y(9)*Y(12)-K51*Y(2)*Y(9); f10=K04*Y(1)*Ne-K48*Y(10)*Y(6)-K49*Y(10)*Y(13)-K50*Y(10)*Y(12)*M; f11=(K01+K03+4*K04)*Y(1)*Ne+K6*Y(7)*Y(12)+K12*Y(8)*Y(6)+2*K13*Y(8)*Y(6)+K14*Y(8)*Y(13)+K17*Y(9)*Y(12)+K18*Y(16)*M+K24*Y(14)*M+K28*Y(14)*Y(12)+K38*Y(12)*Y(13)+K40*Y(3)*Y(12)+K41*Y(4)*Y(13)+K43*Y(13)*Y(3)+K49*Y(10)*Y(13)-K1*Y(1)*Y(11)-K4*Y(7)*Y(11)-K10*Y(8)*Y(11)-K21*Y(15)*Y(11)-K26*Y(14)*Y(11)-K32*Y(13)*Y(11)*M-K33*Y(12)*Y(11)*M-2*K35*Y(11)*Y(11)*M-K36*Y(11)*Y(6)*M-K39*Y(11)*Y(13)-K44*Y(17)*Y(11)-K45*Y(17)*Y(11); f12=K05*Y(2)*Ne+K16*Y(9)*Y(6)+K37*Y(13)*Y(13)+K39*Y(11)*Y(13)+K48*Y(10)*Y(6)-K2*Y(1)*Y(12)-K6*Y(7)*Y(12)-K17*Y(9)*Y(12)-K20*Y(15)*Y(12)-K27*Y(14)*Y(12)-K28*Y(14)*Y(12)-K33*Y(12)*Y(11)*M-2*K34*Y(12)*Y(12)*M-K38*Y(12)*Y(13)-K40*Y(3)*Y(12)-K46*Y(17)*Y(12)-K50*Y(10)*Y(12)*M; f13=K2*Y(1)*Y(12)+K7*Y(7)*Y(6)+K8*Y(7)*Y(17)+K11*Y(8)*Y(6)+K12*Y(8)*Y(6)+K15*Y(9)*Y(6)+K20*Y(15)*Y(12)+K27*Y(14)*Y(12)+K33*Y(12)*Y(11)*M+K40*Y(3)*Y(12)+2*K44*Y(17)*Y(11)+K46*Y(17)*Y(12)-K3*Y(1)*Y(13)-K5*Y(7)*Y(13)-K14*Y(8)*Y(13)-K22*Y(15)*Y(13)-K25*Y(14)*Y(13)-K32*Y(13)*Y(11)*M-2*K37*Y(13)*Y(13)-K38*Y(12)*Y(13)-K39*Y(11)*Y(13)-K41*Y(4)*Y(13)-K43*Y(13)*Y(3)-K47*Y(17)*Y(13)-K49*Y(10)*Y(13); f14=K11*Y(8)*Y(6)+K16*Y(9)*Y(6)+K20*Y(15)*Y(12)+K21*Y(15)*Y(11)+K22*Y(15)*Y(13)+K23*Y(15)*Y(7)+K51*Y(2)*Y(9)-K24*Y(14)*M-K25*Y(14)*Y(13)-K26*Y(14)*Y(11)-K27*Y(14)*Y(12)-K28*Y(14)*Y(12)-K29*Y(14)*Y(6)-K30*Y(14)*Y(17)-K31*Y(14)*Y(7); f15=K5*Y(7)*Y(13)+K6*Y(7)*Y(12)+K7*Y(7)*Y(6)+K14*Y(8)*Y(13)+K18*Y(16)*M+K19*Y(16)*Y(6)+K30*Y(14)*Y(17)-K20*Y(15)*Y(12)-K21*Y(15)*Y(11)-K22*Y(15)*Y(13)-K23*Y(15)*Y(7); f16=K8*Y(7)*Y(17)-K18*Y(16)*M-K19*Y(16)*Y(6); f17=K19*Y(16)*Y(6)+K29*Y(14)*Y(6)+K36*Y(11)*Y(6)*M-K8*Y(7)*Y(17)-K9*Y(7)*Y(17)-K30*Y(14)*Y(17)-K44*Y(17)*Y(11)-K45*Y(17)*Y(11)-K46*Y(17)*Y(12)-K47*Y(17)*Y(13); dYdt=[f1;f2;f3;f4;f5;f6;f7;f8;f9;f10;f11;f12;f13;f14;f15;f16;f17];

tammy · 发表于 2005-12-20 21:14

表示其解的结果是显示的数值乘以1.0e-006 你输入a=[0.0000001 0.0000002]看看结果是什么样的就知道这表示什么意思了

mathtype · 发表于 2005-12-21 11:30

谢谢

账号		自动登录	找回密码
密码			我要加入

请教matlab解常微分方程组

马上注册，结交更多好友，享用更多功能，让你轻松玩转社区。

回复：（mathtype）请教matlab解常微分方程组

回复：（tammy）回复：（mathtype）请教matlab解常微分...