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 ...

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