跪求助：关于lsqcurvefit命令拟合的问题

lysgreat

求助：关于lsqcurvefit命令拟合的问题
已知自变量x和对应函数y的一组数据值，利用函数拟合命令lsqcurvefit拟合出表达式中的两个参数a和b的大小。
其中，函数表达式是：

y=a*(1-erf(x/2/sqrt(bt)))

其中 erf(x)=2/sqrt(pi)*integral from 0 to x of exp(-t^2) dt.
t为常数 t=4*365*24*3600
也可见链接：
http://home.simwe.com/space-419822-do-album-picid-35278.html


自变量：x=[10 20 30 40 50 60 70]*10^(-3)
函数值：y=[5.1 3.9 3.3 2.2 1.7 1.6 1.3]
自己编的程序如下：


xdata=[10  20 30 40 50 60 70]*10^(-3);
Ccl=[5.1 3.9 3.3 2.2 1.7 1.6 1.3];

x=xdata;
y=Ccl;

a0= [1.5 ,  4];
% 初识揣测

[a, resnorm] = lsqcurvefit(@myfun , a0 , x , y)

子程序myfun.m  为
function  F = myfun(a, data)
x = data;
t =4*365*24*3600;
F= a(1)*(1-erf(x./2/sqrt(a(2)*t)));

但是不知道为何，总是有错误，请各位大侠帮帮忙啊，有着急啊，谢谢！！

messenger

什么错误呀？贴出来呀，这么大段程序，难道让别人替你运行一下？

lysgreat

谢谢！！错误提示为：

??? Undefined function or method 'erf' for input arguments of type 'double' and attributes 'full row complex'.

Error in ==> myfun at 6
F= a(1)*(1-erf(x./2/sqrt(a(2)*t)));

Error in ==> lsqcurvefit>objective at 266
         F = feval(funfcn_x_xdata{3},x,XDATA,varargin{:});

Error in ==> snls at 362
            newfvec = feval(funfcn{3},xcurr,varargin{:});

Error in ==> lsqncommon at 149
    [xC,FVAL,LAMBDA,JACOB,EXITFLAG,OUTPUT,msgData]=...

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

Error in ==> runme at 10
[a, resnorm] = lsqcurvefit(@myfun,a0,data,zdata)

这里lsqcurvefit拟合命令运行出错了。。。
运行中出现复数吗？erf命令用不了复数吗？

messenger

是的，erf命令用不了复数。你可以令a(2)=/sqrt(a(2)*t))./2，得到拟合结果后再反算原来的a(2)。即采用F= a(1)*(1-erf(x.*a(2)));来拟合，可以得到比较理想的拟合结果：
a =
    5.809752094073831    13.84893396425852
resnorm =
    0.3.412381193673270

lysgreat

好主意啊！！
谢谢 messenger 总版主！！！

的确可以成功运行啊！！