实验七 无约束优化

K=[0.0910 0.2543 0.3121 0.3792 0.4754 0.4410 0.4517 0.5595 0.8080 1.307 1.7042 2.0019 2.2914 2.4941 2.8406 2.9854 3.2918 3.7314 4.3500]';

L=[4.8179 4.9873 5.1282 5.2783 5.4334 5.5329 6.4749 6.5491 6.6152 6.6808 6.7455 6.8065 6.8950 6.9820 7.0637 7.1394 7.2085 7.3025 7.3740]';

y=log(Q);

f=[ones(19,1),log(K),log(L)]; aa=f\\y

aa(1)=exp(aa(1)); aa a11 aa =

-1.1556 0.6600 1.2677 aa =

0.3149 0.6600 1.2677 用非线性方法

function f=exam01fun(x,t) f=x(1)*t(1,:).^(x(2)).*t(2,:).^(x(3)); m文件

c=[0.7171 0.8964 1.0202 1.1962 1.4928 1.6909 1.8548 2.1618 2.6638 3.4634 4.6759 5.8478 6.7885 7.4463 7.8345 8.2068 9.9468 9.7315 10.4791];

K=[0.0910 0.2543 0.3121 0.3792 0.4754 0.4410 0.4517 0.5595 0.8080 1.307 1.7042 2.0019 2.2914 2.4941 2.8406 2.9854 3.2918 3.7314 4.3500];

L=[4.8179 4.9873 5.1282 5.2783 5.4334 5.5329 6.4749 6.5491 6.6152 6.6808 6.7455 6.8065 6.8950 6.9820 7.0637 7.1394 7.2085 7.3025 7.3740];

x0=[1,0.8,0.8];

[x,norm]=lsqcurvefit(@exam01fun,x0,[K;L],c) x =

0.9666 0.7812 0.6521

norm =

2.7492