24小时热门版块排行榜

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

beina

新虫 (正式写手)

应助: 3 (幼儿园)
金币: 55
散金: 254
红花: 2
帖子: 570
在线: 163.4小时
虫号: 1377290
注册: 2011-08-23
性别: MM
专业: 敏感电子学与传感器

[求助] 需要用matlab拟合数据的外行人急需帮助

我有一函数y=f(x),里面有个未知数k，还有个x,y数值的txt文件，用这些条件求k
首先建立了一个M文件，
function F=myfun(x,xdata)
F=0.6518 + 0.7713./atan(xdata) - 0.0911*atan(xdata) + ((0.9121* xdata + 1.572)./(0.9121* xdata + 1.59 + 2.3026./ xdata + 0.0637*log(1 + 15.7./ xdata)) - 1)./((1 + 5.04*x* xdata).*(1 + (x^0.674)* (xdata.^0.173)));
然后就运行下面的命令：
>> load data.txt
>> xdata=data(:,1);
>> ydata=data(:,2);
>> x=lsqcurvefit(@myfun,0.05,xdata,ydata)
然后就有错误了
??? Error using ==> qr
Complex sparse QR is not yet available.

Error in ==> aprecon at 57
RPCMTX = qr(TM(:,p));

Error in ==> trdog at 47
[R,permR] = feval(pcmtx,H,pcoptions,DM,DG,varargin{:});

Error in ==> snls at 346
[sx,snod,qp,posdef,pcgit,Z] = trdog(x,g,A,D,delta,dv,...

Error in ==> lsqncommon at 285
[x,FVAL,LAMBDA,JACOB,EXITFLAG,OUTPUT,msg]=...

Error in ==> lsqcurvefit at 186
[x,Resnorm,FVAL,EXITFLAG,OUTPUT,LAMBDA,JACOB] = ...

这是哪里的错误呢， lsqcurvefit 的运用中有木有语法错误，哪位大侠可以帮我解答一下，感激涕零，如果能帮我跑一下程序就最好不过了……