lsqcurvefit函数拟合

用lsqcurvefit函数拟合下沉参数时遇到的问题,求解答

X=xlsread('C:\Users\YDH\Desktop\高程.xlsx','sheet1','A2:A15');
X1=xlsread('C:\Users\YDH\Desktop\高程.xlsx','sheet2','A2:A15');
Y3=1000*(X1-X);
fun= inline ('-3400/2*beta(1)*(erf((X/1000-beta(3)-0).*(sqrt(pi)/beta(2)))+1)+3400/2*beta(1)*(erf(((X/1000-(275+beta(4))+0)*sin((beta(5)+1)*pi/180)/sin((beta(5))*pi/180)).*(sqrt(pi)/beta(2)))+1)','beta','X');
x0=[0.7,50,5,5,50];
parm=lsqcurvefit(fun,x0,X,Y3);
H=120;
beta(1)=parm(1);
beta(2)=parm(2);
beta(3)=parm(3);
beta(4)=parm(4);
beta(5)=parm(5);
q=beta(1);
tanb=H/(beta(2)*1000);
S1=beta(3);
S2=parm(4);
the=parm(5);
parameters1=[q,tanb1,S3,S4,tanb2];
W=-3400/2* beta(1)*(erf((X/1000-beta(3)-0).*(sqrt(pi)/beta(2)))+1)+3400/2*beta(1)*(erf(((X/1000-(275+beta(4))+0)*sin((beta(5)+1)*pi/180)/sin((beta(5))*pi/180)).*(sqrt(pi)/beta(2)))+1);
plot(X,Y3,'o',X,W,'r');
xlabel('倾向 L (m)');
ylabel('下沉 W (mm)');
set(gca,'XTick',-250:50:450);set(gca,'YTick',-3000:500:500);%坐标轴刻度控制
legend('原始数据','拟合曲线');

运行结果及详细报错内容

出错 nh (line 7)
parm=lsqcurvefit(fun,x0,X,Y3);

原因:
Failure in initial objective function evaluation. LSQCURVEFIT cannot continue.