1stopt的初值设定
最近在解一组方程，可是用1stopt程序运行后,结果很奇怪，各参数甚至不在设定的初始范围内，请问是为什么呢？使用的是1stop1.5PJ版，会是版本的问题么？以下是程序内容，能否请大家帮忙看一下呢：Parameter x(1:5) , x(11:17) , x23 , x(29:31), x(35:37), x(40:42) ;
Y(1:6) ;
x6 , x7 ,x8 ,x9 ,x10=500e-6 ;
x18,x19 ,x20 ,x21 ,x22=500e-6 ;
x24 , x25 , x26 , x27 , x28 ;
x32 ,x33 ,x34 ;
x38=3500 ,x39=3500,x43 , x44;
Constant P=68 , Ramb=1000e-06,
Rsurf=100e-06 , Tamb=2550,TCZ=4000 ,T=298.15,
rod=0.089, bet=0.02, la=0.6, hsCO=-3950 ,
hsCO2=-8950 , hsH2=0 ,sH2O=-13443.9,
hsAl2O= -1954.3 , hsAl= 12005.2, hsAlxOx = -2848.6,
CpsCO=1.04, CpsCO2=0.84, CpsH2=14.43 , CpsH2O=1.87 , CpsAl2O=0.62,
CpsAl=0.79,
CpsAlxOx = 0.63,
CpCO=1.39,CpCO2=1.53, CpH2=17.6 , CpH2O=2.98 , CpAl=0.77,
CpAlxOx=0.95, hAl2O3TCZ= -10447 , hAlTsurf = 2765.2,
aco=-4260.076, aco2=-9200.45, ah2=-4302.3, ah2o=-14001.44,
aalxox=-3036.43, aal2o=-2139.15,aal=11775.62, aalo=1812.5,
hAl2O3Tsurf=-11153.2, CpAl2O=0.83, CpAlO=0.9;
Function x1*(1/x10-1/Ramb)/(4*pi*rod*(1-bet)) = Ln ((x1*0.18-x2)/(x1*x6-x2)); // 1
x1*(1/x10-1/Ramb)/(4*pi*rod*(1-bet)) = Ln ((x1*0.174-x3)/(x1*x7-x3)); // 2
x1*(1/x10-1/Ramb)/(4*pi*rod*(1-bet)) = Ln ((x1*0.008-x4)/(x1*x8-x4)); // 3
x1*(1/x10-1/Ramb)/(4*pi*rod*(1-bet)) = Ln ((x1*0.27-x5)/(x1*x9-x5)); // 4
x11*(1/x22-1/x10)/(4*pi*rod*(1-bet)) = Ln ((x11*x6-x12)/(x11*x18-x12)); // 5
x11*(1/x22-1/x10)/(4*pi*rod*(1-bet)) = Ln ((x11*x7-x13)/(-x13)); // 6
x11*(1/x22-1/x10)/(4*pi*rod*(1-bet)) = Ln ((x11*x8-x14)/(x11*x19-x14)); // 7
x11*(1/x22-1/x10)/(4*pi*rod*(1-bet)) = Ln ((x11*x9-x15)/(-x15)); // 8
x11*(1/x22-1/x10)/(4*pi*rod*(1-bet)) = Ln (-x16/(x11*x20-x16)); // 9
x11*(1/x22-1/x10)/(4*pi*rod*(1-bet)) = Ln (-x17/(x11*x21-x17)); //10
x23*(1/Rsurf-1/x22)/(4*pi*rod*(1-bet)) = Ln (x18/x24); // 11
x23*(1/Rsurf-1/x22)/(4*pi*rod*(1-bet)) = Ln (x19/x25); // 12
x23*(1/Rsurf-1/x22)/(4*pi*rod*(1-bet)) = Ln ((x23*x20-x29)/(x23*x26-x29)); // 13
x23*(1/Rsurf-1/x22)/(4*pi*rod*(1-bet)) = Ln ((x23*x21-x30)/(x23*x27-x30)); // 14
x23*(1/Rsurf-1/x22)/(4*pi*rod*(1-bet)) = Ln (-x31/(x23*x28-x31)); // 15
x1=x2+x3+x4+x5; // 16
x11=x12+x13+x14+x15+x16+x17; // 17
x23=x29+x30+x31; // 18
1=x6+x7+x8+x9+x32; // 19
1=x18+x19+x20+x21+x33; // 20
1=x24+x25+x26+x27+x28+x34; // 21
Y1=x2*CpCO+x3*CpCO2+x4*CpH2+x5*CpH2O; // 44
Y3=x12*CpCO+x13*CpCO2+x14*CpH2+x15*CpH2O+x16*CpAlxOx+x17*CpAl2O; // 45
Y5=x29*CpAlO+x30*CpAl2O+x31*CpAl; // 46
Y2=x2*aco+x3*aco2+x4*ah2+x5*ah2o; // 47
Y4=x12*aco+x13*aco2+x14*ah2+x15*ah2o+x16*aalxox+x17*aal2o; // 48
Y6=x29*aalo+x30*aal2o+x31*aal; // 49
Y1*(1/x10-1/Ramb)/(4*pi*(1-bet)*la) = Ln ((Y1*Tamb+Y2-x35)/(Y1*TCZ+Y2-x35)); // 22
Y3*(1/x22-1/x10)/(4*pi*(1-bet)*la) = Ln ((Y3*TCZ+Y4-x36)/(Y3*x38+Y4-x36)); // 23
Y5*(1/Rsurf-1/x22)/(4*pi*(1-bet)*la) = Ln ((Y5*x38+Y6-x37)/(Y5*x39+Y6-x37)); // 24
x37=0.593*(x29*hAl2O3Tsurf)+x31*hAlTsurf; // 25
x36=x37; // 26
x35=x36-x40*hAl2O3TCZ; // 27
x2/28.011 = -x3/44.011; // 28
x4/2.016 = -x5/18.016; // 29
(x13-x3)/44.011 = (x2-x12)/28.011; // 30
(x15-x5)/18.016 = (x4-x14)/2.016; // 31
x16/86+x17/35=(x13-x3)/44.011 +(x15-x5)/18.016; // 32
x40=(102/70)*x17+(102/86)*x16; // 33
x30=-(70/172)*x29; // 34
x41=-1.6e-4; //50
x41=-x31; // 35
x42=(-102/172)*x29; // 36
x31/27=(x16-x29)/43; // 37
x30=x17; // 38
x31/27 = -x13/44.011-x15/18.016 ; // 39
x43/(x43+x44) = 0.122+0.000577*(x38-3500); // 40
-x29=(3653*(Rsurf^2)*P*(1-bet)*x26/43)/((x39^0.5）*（x24/28.011+x25/2.016+x26/43+x27/70+x28/27+x34/33.28)); // 41
P*(x28/27)/(x24/28.011+x25/2.016+x26/43+x27/70+x28/27+x34/33.28)= exp(12.43-34680/x39); // 42
x43+x44=x20; // 43
本帖最后由 dingd 于 2011-6-10 22:45 编辑
meconglin 发表于 2011-6-10 17:03 http://www.chinavib.com/static/image/common/back.gif
最近在解一组方程，可是用1stopt程序运行后,结果很奇怪，各参数甚至不在设定的初始范围内，请问是为 ...
Parameter x(1:5) , x(11:17) , x23 , x(29:31), x(35:37), x(40:42) ;
Y(1:6) ;
x6 , x7 ,x8 ,x9 ,x10=500e-6 ;
x18,x19 ,x20 ,x21 ,x22=500e-6 ;
x24 , x25 , x26 , x27 , x28 ;
x32 ,x33 ,x34 ;
x38=3500 ,x39=3500,x43 , x44;
改为：
Parameter x(1:5) , x(11:17) , x23 , x(29:31), x(35:37), x(40:42),
Y(1:6) ,
x6 , x7 ,x8 ,x9 ,x10=500e-6 ,
x18,x19 ,x20 ,x21 ,x22=500e-6,
x24 , x25 , x26 , x27 , x28 ,
x32 ,x33 ,x34 ,
x38=3500 ,x39=3500,x43 , x44;
回复 2 # dingd 的帖子
谢谢dingd主任！
页:
[1]