用nlinfit多元函数拟合提示秩不足怎么解决

程序如下
主程序
clc;clear;
PA=[6.18 6.32 11.9 12.7 6.31 6.95 12.7 13.1 9.58 9.81 4.65 14.6 9.67 9.98 9.66 9.78 9.60 9.63 5.65 14.3 5.65 20.0 10.1 19.8];
PB=[2.52 2.45 3.36 3.05 7.91 8.77 9.14 9.18 6.64 5.91 5.56 6.92 0.714 10.7 5.85 6.14 5.95 5.92 4.52 2.24 2.24 2.82 1.00 3.14];
PH=[14.0 27.5 14.3 27.3 15.3 27.4 14.5 27.6 9.07 32.0 20.4 20.5 27.0 21.4 20.8 20.4 20.8 20.8 20.4 7.25 28.4 17.1 17.0 9.87];
ra=[3.541 2.397 1.021 4.722 0.593 0.268 2.797 2.451 3.196 2.021 0.896 5.084 5.686 1.193 2.648 3.303 3.054 3.302 1.271 11.648 2.002 9.604 7.754 11.59];
x=[PA;PB;PH];
beta=nlinfit(x,ra,@myfuna,[1 1 1 1 1]);
KA=beta(1);
KB=beta(2);
KH=beta(3);
Ks=beta(4);
Kp=beta(5);

调用程序
function ra=myfuna(const,x)
KA=const(1);
KB=const(2);
KH=const(3);
Ks=const(4);
Kp=const(5);
ra=Ks*KA.*(x(1,:)-x(2,:)/Kp)./((1+KA*x(1,:)+KB*x(2,:)+KH*x(3,:)).^2);
end

这个时候就提示秩不足，5个参数中Kp的数值错误，怎么解决啊....