24小时热门版块排行榜    

查看: 1059  |  回复: 2

陈苏杭

新虫 (初入文坛)

[求助] 1stOpt 5.0计算五个自变量五个应变量,调用五组数据 已有1人参与

用1stOpt 5.0计算活化能,指前因子,反应级数,五个升温速率,分别是2k, 5k, 10k, 20k, 40k, 需要多个自变量、多个因变量调用不同的数据。x1, y1, x2, y2, x3, y3, x4, y4, x5, y5的数据依次如附录

要求:计算五个方程,共同的一个a,b, n
五个方程,共同的一个拟合度R

以下代码运行不出来

Variable x1, y1, x2, y2, x3, y3, x4, y4, x5, y5;
Parameter 80000<=[a]<=130000, 10^10<=<=10^14, 2<=[n]<=3.5;
Function y1= exp((-b/2*(x1*exp(-a/(8.314*x1))-403.15*exp(-a/(8.314*403.15))))^n); //热分解机理函数积分而来
             y2= exp((-b/5*(x2*exp(-a/(8.314*x2))-413.15*exp(-a/(8.314*413.15))))^n);
           y3= exp((-b/10*(x3*exp(-a/(8.314*x3))-429.75*exp(-a/(8.314*429.75))))^n);
            y4= exp((-b/20*(x4*exp(-a/(8.314*x4))-433.15*exp(-a/(8.314*433.15))))^n);
          y5= exp((-b/40*(x5*exp(-a/(8.314*x5))-441.15*exp(-a/(8.314*441.15))))^n);
data; //依次为 x1, y1, x2, y2, x3, y3, x4, y4, x5, y5数据
403.15        1.00004        413.15        0.99984        429.75        1.00061        433.15        0.99979        441.15        0.99942
404.15        1.00105        414.15        0.99897        433.15        0.99103        435.15        0.98674        442.15        1.00403
405.15        0.99787        415.15        1.00601        434.15        0.99645        437.15        0.98119        443.15        1.00675
406.15        0.98376        416.15        0.99449        435.15        0.98815        439.15        0.97587        444.15        1.00553
407.15        0.98109        417.15        0.98662        436.15        0.98138        445.15        0.97063        445.15        0.99903
408.15        0.97478        418.15        0.97988        437.15        0.96687        446.15        0.96606        446.15        0.99455
409.15        0.97709        419.15        0.97337        438.15        0.94818        447.15        0.95996        447.15        0.99076
410.15        0.96691        420.15        0.97106        439.15        0.93716        448.15        0.95598        448.15        0.98797
411.15        0.95392        421.15        0.96651        440.15        0.92396        449.15        0.95116        449.15        0.98879
412.15        0.94948        422.15        0.95779        441.15        0.90605        450.15        0.9421        450.15        0.98815
413.15        0.94018        424.15        0.93805        442.15        0.88589        451.15        0.92503        451.15        0.98878
414.15        0.93173        425.15        0.93542        443.15        0.86477        452.15        0.8881        452.15        0.98788
415.15        0.906        427.15        0.91393        444.15        0.83359        453.15        0.85772        453.15        0.9829
416.15        0.90224        429.15        0.89408        445.15        0.81462        454.15        0.83326        454.15        0.97599
417.15        0.89082        430.15        0.8729        446.15        0.78029        455.15        0.80258        455.15        0.96746
418.15        0.87288        431.15        0.863        447.15        0.75791        456.15        0.78046        456.15        0.96204
419.15        0.8598        432.15        0.8493        448.15        0.70797        457.15        0.75461        457.15        0.96118
420.15        0.84169        433.15        0.83477        449.15        0.66485        458.15        0.72252        458.15        0.95726
421.15        0.81976        434.15        0.81786        450.15        0.61271        459.95        0.65375        459.15        0.95386
421.35        0.81023        435.15        0.78143        450.35        0.59184        460.95        0.56292        460.15        0.94962
421.55        0.80226        436.15        0.74931        450.55        0.588        461.95        0.49363        462.15        0.93748
421.75        0.80185        437.15        0.72903        450.75        0.577        462.95        0.38784        464.15        0.89916
421.95        0.79957        438.15        0.70264        450.95        0.55669        463.95        0.30063        465.15        0.87365
422.15        0.80029        439.95        0.61315        451.15        0.54128        464.95        0.23038        466.15        0.84428
423.15        0.76824        440.95        0.55369        452.15        0.48982        465.95        0.14879        467.95        0.79174
424.15        0.73472        441.95        0.50221        453.15        0.40817        466.95        0.10936        469.95        0.72082
425.15        0.70654        442.95        0.451        454.15        0.33027        467.95        0.07257        470.15        0.71299
426.15        0.67004        443.95        0.36688        455.15        0.21011        468.95        0.0647        471.15        0.66708
427.15        0.63079        444.95        0.26829        456.15        0.12704        470.95        0.03485        471.35        0.65782
428.15        0.58119        445.95        0.17802        457.15        0.08023        471.95        0.01929        471.55        0.64801
429.15        0.53237        446.95        0.10966        458.95        0.0409        473.95        0.0083        471.75        0.63714
430.15        0.46876        447.95        0.05489        459.95        0.0207        475.95        0.00408        471.95        0.62428
431.15        0.39331        449.55        0.0338        460.95        0.01197        477.95        0.00196        473.95        0.4911
432.15        0.31175        450.55        0.00203        462.15        3.73E-4        479.95        0.04789   475.95        0.32183
433.15        0.23106        450.95        0.00141        465.15        -1.36E-4   480.95        0.01814        477.95        0.15373
434.15        0.14065        451.95        0.00174                                        481.95        0.00288
435.15        0.07112        452.15        0.00216                                        482.15        6.82982E-6
436.15        0.03465        452.35        0.0029                                        482.367        2.39643E-5
437.15        0.00577        452.55        0.00499                                        486.367        0.00572
437.55        0.00458        455.15        -7.27674E-5                               
438.15        0.00565        455.35        -7.05623E-5                               
439.15        0.00257        456.95        -9.47078E-4                                               
440.15        -0.00866        457.15        -1.35612E-4                                               
441.15        -0.00426                                                               
443.15        -0.00637
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

独孤神宇

版主 (知名作家)

【答案】应助回帖

感谢参与,应助指数 +1
Variable x1, y1, x2, y2, x3, y3, x4, y4, x5, y5;
Parameter 80000<=[a]<=130000, 10^10<=<=10^14, 2<=[n]<=3.5;
SharedModel;
Function y1= exp((-b/2*(x1*exp(-a/(8.314*x1))-403.15*exp(-a/(8.314*403.15))))^n); //热分解机理函数积分而来
             y2= exp((-b/5*(x2*exp(-a/(8.314*x2))-413.15*exp(-a/(8.314*413.15))))^n);
           y3= exp((-b/10*(x3*exp(-a/(8.314*x3))-429.75*exp(-a/(8.314*429.75))))^n);
            y4= exp((-b/20*(x4*exp(-a/(8.314*x4))-433.15*exp(-a/(8.314*433.15))))^n);
          y5= exp((-b/40*(x5*exp(-a/(8.314*x5))-441.15*exp(-a/(8.314*441.15))))^n);
Data;
x1, y1, x2, y2, x3, y3, x4, y4, x5, y5
数值计算
2楼2018-08-06 09:27:08
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

陈苏杭

新虫 (初入文坛)

引用回帖:
2楼: Originally posted by 独孤神宇 at 2018-08-06 09:27:08
Variable x1, y1, x2, y2, x3, y3, x4, y4, x5, y5;
Parameter 80000&lt;=&lt;=130000, 10^10&lt;=&lt;=10^14, 2&lt;=&lt;=3.5;
SharedModel;
Function y1= exp((-b/2*(x1*exp(-a/(8.31 ...

谢谢,但还是运行不了啊
3楼2018-08-09 10:32:18
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 陈苏杭 的主题更新
信息提示
请填处理意见