|
|
【答案】应助回帖
lihongru(金币+30): 回去试一下，不知是否会因为初值不好而不能拟合。 2011-10-27 13:02:00
clear,clc
F=@(q,x) q(1)*log(x(:,1)-q(1)*x(:,2))+q(2)*x(:,3)+q(3)-log(x(:,2))
x=;
y=;
T=;
q0=[];%初值
warning off
q=nlinfit([x y T],zeros(size(x,1),1),F,q0);
disp('k a b分别为')
disp(num2str(q)); |
|
回复此楼
