function y= myfun(beta,x)
b1 = beta(1);
b2 = beta(2);
b3 = beta(3);
b4 = beta(4);
b5 = beta(5);

x1 = x(:,1);
x2 = x(:,2);
x3 = x(:,3);

y = (b1*x2 - x3/b5)./(1+b2*x1+b3*x2+b4*x3); %模型方程

%待拟合的自变量
x0=[
   470   300    10
   285    80    10
   470   300   120
   470    80   120
   470    80    10
   100   190    10
   100    80    65
   470   190    65
   100   300    54
   100   300   120
   100    80   120
   285   300    10
   285   190   120];
%待拟合的函数值
y0=[
8.55
3.79
4.82
0.02
2.75
14.39
2.54
4.35
13
8.5
0.05
11.32
3.13];
%参数初值
beta0=[1 0.05 0.02 0.1 2]';
%调用
beta = nlinfit(x0,y0,@myfun,beta0);
%相关系数
x1=x0(:,1);
x2=x0(:,2);
x3=x0(:,3);
y=(beta(1)*x2 - x3/beta(5))./(1+beta(2)*x1+beta(3)*x2+beta(4)*x3);%估计值
R=corrcoef(y0,y)