banana = @(x)100*(x(2)-x(1)^2)^2+(1-x(1))^2;
[x,fval] = fminsearch(banana,[-1.2, 1])
x =
    1.0000    1.0000
fval =
    8.1777e-010

4楼: Originally posted by dingd at 2016-08-29 16:35:33
有数据吗？包括三次不同的T值，一起放上来看看

5楼: Originally posted by duu123 at 2016-09-02 09:37:47
麻烦帮忙看一下，我试了很长时间都没弄出来。...

6楼: Originally posted by dingd at 2016-09-02 10:38:15
公式里的Rg是什么？是待求参数还是已知常数？...

VarConstant T=[303,323,343];
Constant Rg=8.314;
ConstStr k=k0*exp(-H/(Rg*T));
Variable P,q;
Function q=qm*k*p/(1+(k*p)^n)^(1/n);
Data;
//P        q
0.89        0.0319
4.15        0.0511
8.77        0.0866
16.23        0.2242
39.31        0.396
57.84        0.6076
79.32        0.7489
103.26        0.8667
195.4        1.1721
306.1        1.3996
403.4        1.5746
509.68        1.7066
610.4        1.8212
706.67        1.9082
Data;
//P        q
10.18        0.0966
20.17        0.1815
38.72        0.3208
64.27        0.517
79.6        0.538
106.98        0.6583
197.4        0.9437
302.85        1.1369
386.16        1.2879
501.96        1.4217
578.7        1.5101
707.29        1.6153
Data;
//P        q
10.84        0.0694
21.58        0.1328
40        0.2219
62.76        0.3243
81.27        0.4026
102.2        0.481
150.66        0.6438
218.15        0.7933
318.51        0.9567
437.4        1.1178
518.7        1.2097
594.17        1.2903
716.44        1.3903

8楼: Originally posted by dingd at 2016-09-02 11:38:16
1stOpt计算：

VarConstant T=;
Constant Rg=8.314;
ConstStr k=k0*exp(-H/(Rg*T));
Variable P,q;
Function q=qm*k*p/(1+(k*p)^n)^(1/n);
Data;
//P        q
0.89        0.0319
4.15        0.0511
8.77        0.0866
16.23        0.2 ...