matlab 无量纲,[MATLAB数学相关] 求助 （已经-无量纲化）的非线性方程组（已经仔细检查，能算，但exitflag=0求高手…

其实是后边要用matcont做常微分方程组延拓分析的。但是那个延拓软件必须有一个稳态解(充分发展的，dx/dt=0时的解)作为初值，才能开始延拓。

代码如下：已经仔细检查了一周了。能算，但是exitflag=0，估计初值太敏感。菜鸟求助啊。。。(%后边是多余的项，不用管，直接贴代码进编辑器就行。。)

小女再次感谢！

steady.mX0=[0.278212173322273,0.912568632338739,0.085310950018308,0.912568731077048,0.087431259625095,0.087431260178122,0.499999721583934,0.000000000015842,0.087431268922953,0.053882265858466,0.087432024167792,0.087431276521867,0.500000000000000,0.087431268922956,0.087431268922956,0.500000000000000,0.087431273571452,0.500000000000001,0.078613189530066];

options=optimset(‘Display’,’iter’,’MaxFunEvals’,100000,’MaxIter’,800)

%’MaxFunEvals’,100000000);

[x,fval,exitflag] =fsolve(‘steadys’,X0,options)

steadys.m

function pp=steadys(x)

pp(1)=(1-x(1))+0.00001*48.432471*(0.0132913*x(3)*x(2)*x(7)+0.0132913*1.6306e+04*2.71828^(-9.6993/x(1))*x(2)*x(5)+0.0132913*1.6306e+04*2.71828^(-9.6993/x(1))*1.6037e+08*2.71828^(-18.8930/x(1))*x(6)+5.781275*1.6306e+04*2.71828^(-9.6993/x(1))*x(2)*x(10))-2.782248*(x(1)-x(19));

pp(2)=1-x(2)-8.522*10^(-3)*1.3028e+12*2.71828^(-27.8955/x(1))*x(2)^2-5.6228*10^(-3)*1.6037e+08*2.71828^(-18.8930/x(1))*x(2)*x(3)-0.0132913*1.6306e+04*2.71828^(-9.6993/x(1))*x(2)*x(7)-0.0132913*1.6306e+04*2.71828^(-9.6993/x(1))*x(2)*x(6)-5.781275*1.6306e+04*2.71828^(-9.6993/x(1))*x(2)*x(10)-2.493604*10^(-4)*6.6735e+19*2.71828^(-45.6473/x(1))*x(9);

%-y7*c4*x(2)*x(10)

pp(3)=-x(3)+1.1722874*1.3028e+12*2.71828^(-27.8955/x(1))*x(2)^2-1.546946*1.6037e+08*2.71828^(-18.8930/x(1))*x(2)*x(3)-3.198*10^(-4)*0.01*x(3)*x(8);

%-y10*c5*x(3)*x(10);

pp(4)=1-x(4)-65.3447*1.1663e+16*2.71828^(-36.9952/x(1))*x(4)-0.01923587*17.6895*2.71828^(-2.8730/x(1))*x(8)*x(6);

pp(5)=-x(5)+2.445938*1.6037e+08*2.71828^(-18.8930/x(1))*x(2)*x(3)-5.781275*1.6306e+04*2.71828^(-9.6993/x(1))*x(2)*x(5)-8.496798*17.6895*2.71828^(-2.8730/x(1))*x(5)*x(8)+3.2672*10^10*1.1663e+16*2.71828^(-36.9952/x(1))*x(9)-(4.54117*10^(-3)*8.0952*2.71828^(-2.0910/x(1)))*x(5)*x(10);

%+y16*c4*x(2)*x(10);+y18*c11

pp(6)=-x(6)-5.781275*1.6306e+04*2.71828^(-9.6993/x(1))*x(2)*x(6)-2.8425*10^10*1.1663e+16*2.71828^(-36.9952/x(1))*x(4)-8.496798*10^6*17.6895*2.71828^(-2.8730/x(1))*x(6)*x(8);

pp(7)=-x(7)*2.44595*10^8*1.6037e+08*2.71828^(-18.8930/x(1))*x(2)*x(4)-5.781275*1.6306e+04*2.71828^(-9.6993/x(1))*x(2)*x(7)-(4.54117*10^(-3)*8.0952*2.71828^(-2.0910/x(1)))*x(7)*x(10);

%-y23*c5*x(3)*x(10);+y18*c11

pp(8)=-x(8)-5.871963*10^13*17.6895*2.71828^(-2.8730/x(1))*x(8)*x(10)+5.05856*1.1663e+16*2.71828^(-36.9952/x(1))*x(13)+4.515847*10^5*1.1663e+16*2.71828^(-36.9952/x(1))*x(4)-134.9877*17.6895*2.71828^(-2.8730/x(1))*x(5)*x(8)+5.871963*10^13*1.1663e+16*2.71828^(-36.9952/x(1))*x(9)-8.033211*0.01*x(3)*x(8);

pp(9)=-x(9)+1.699359668*10^(-2)*17.6895*2.71828^(-2.8730/x(1))*x(8)*x(5)-65.3447*1.1663e+16*2.71828^(-36.9952/x(1))*x(9)-0.216944*6.6735e+19*2.71828^(-45.6473/x(1))*x(9);

pp(10)=-x(10)+0.0132913*1.6306e+04*2.71828^(-9.6993/x(1))*x(2)*x(7)+0.0132913*1.6306e+04*2.71828^(-9.6993/x(1))*x(2)*x(5)+0.0132913*1.6306e+04*2.71828^(-9.6993/x(1))*x(2)*x(6)*65.3447*1.1663e+16*2.71828^(-36.9952/x(1))*x(13)-x(10)*(8.496798*17.6895*2.71828^(-2.8730/x(1))*x(8)+4.54117*10^(-3)*8.0952*2.71828^(-2.0910/x(1))*x(10));

%+y7*c4*x(2);+y18*c11

pp(11)=-x(11)+0.001*(0.0132913*x(3)*x(2)*x(7)+0.0132913*1.6306e+04*2.71828^(-9.6993/x(1))*x(2)*x(5)+0.0132913*1.6306e+04*2.71828^(-9.6993/x(1))*x(2)*x(6)+5.781275*1.6306e+04*2.71828^(-9.6993/x(1))*x(2)*x(10))+65.3447*1.1663e+16*2.71828^(-36.9952/x(1))*x(14)-x(11)*(8.496798*17.6895*2.71828^(-2.8730/x(1))*x(8)+(4.54117*10^(-3)*8.0952*2.71828^(-2.0910/x(1)))*x(10));

%+y7*c4*x(2);+y32*c5*x(3);+y18*c11

pp(12)=-x(12)+0.00001*(0.0132913*x(3)*x(2)*x(7)+0.0132913*1.6306e+04*2.71828^(-9.6993/x(1))*x(2)*x(5)+0.0132913*1.6306e+04*2.71828^(-9.6993/x(1))*x(2)*x(6)+5.781275*1.6306e+04*2.71828^(-9.6993/x(1))*x(2)*(200*x(11)+x(10))+65.3447*1.1663e+16*2.71828^(-36.9952/x(1))*x(15)-x(12)*(150.3041*2.71828^(-2.8730/x(1))*x(8)+0.036761679*2.71828^(-2.0910/x(1)))*x(10));

%+y7*c4*x(2);+y32*c5*x(3);+y18*c11

pp(13)=-x(13)*8.496798*17.6895*2.71828^(-2.8730/x(1))*x(8)*x(10)-65.3447*1.1663e+16*2.71828^(-36.9952/x(1))*x(13);

pp(14)=-x(14)+8.496798*17.6895*2.71828^(-2.8730/x(1))*x(8)*x(11)-65.3447*1.1663e+16*2.71828^(-36.9952/x(1))*x(14);

pp(15)=-x(15)+8.496798*17.6895*2.71828^(-2.8730/x(1))*x(8)*x(12)-65.3447*1.1663e+16*2.71828^(-36.9952/x(1))*x(14);

pp(16)=-x(16)*(0.5*4.54117*10^(-3)*8.0952*2.71828^(-2.0910/x(1)))*x(10)^2;

%+y18*c11;+x(10)*(y7*c4*x(2)+y32*c5*x(3));

pp(17)=-x(17)+(4.54117*10^(-3)*8.0952*2.71828^(-2.0910/x(1)))*x(10)*x(11);

%+y18*c11;+x(11)*(y7*c4*x(2)+y32*c5*x(3));

pp(18)=-x(18)*(4.54117*10^(-3)*8.0952*2.71828^(-2.0910/x(1)))*x(10)*x(12)+4.54117*10^(-3)*8.0952*2.71828^(-2.0910/x(1))*x(11);

%+y18*c11;+x(12)*(y7*c4*x(2)+y32*c5*x(3))

pp(19)=(0.727148704-12.7016)*x(19)+4.897673*(x(1)-x(19));

相关资源：本草纲目下载李时珍本草纲目查询软件版v1.4_本草纲目pdf彩图版…

来源：weixin_39795419