x(:,1)=[20 20 20 20 40 40 40 40 60 60 60 60]';
x(:,2)=[15 30 45 60 15 30 45 60 15 30 45 60]';
beta0=[1 1 1 1];
y=[0.06 0.076 0.183 0.284 0.37 0.19 0.371 0.556 0.812 0.938 0.296 0.644]';
f=@(beta,x)((x(:,1).*sin(x(:,2))-beta(1)).^2/(2*beta(2))+(x(:,1).^2.*cos(x(:,2)).^2.*sin(beta(3).*x(:,2)))/(2*beta(4)));
[beta,r,j]=nlinfit(x,y,f,beta0);
beta =
1.0e+03 *
0.0094 2.8425 0.0010 1.2828
你看看结果是不是符合,另外你给的Y的维数和x不一样,我去掉了最后三个