[slgu]

先建立函数的m文件：
function fy=qst(a,x)
fy=a(1)*x(:,1).*(1-x(:,2)/a(2))-a(3)*x(:,2);
存成qst.m
然后输入X=[x;z]';
及y=[...]';
用：
b0=[8.69 355 0.39]
再用：
b=nlinfit(X,y,'qst',b0)
b =
8.6925 355.3126 0.3902