救命啊！！谁能帮看看这个是什么原因，到底是哪错了？！！急啊！

林偏殇

原始代码
syms x
x0=[1978;1979;1980;1981;1982;1983;1984;1985;1986;1987;1988;1989;1990;1991;1992;1993;1994;1995;1996;1997;1998;1999;2000;2001;2002;2003;2004;2005;2006;2007;2008;2009;2010];
y0=[566;632;745;755;769;789;985;1110;1313;1428;1782;1920;2150;2292;2601;3149;4338;5145;5809;6241;6854;7656;8772;10007;11374;12567;14332;16614;19228;22844;26404;29688;32074];
f=fittype('a*exp(b*x) + c*exp(d*x)','coeff',{'a','b','c','d'});
cfun=fit(x0,y0,f)
xi=1978:2010;
yi=cfun(xi);
plot(x0,y0,'r*',xi,yi,'b-');
运行之后，出现很多错误。。
我想通过上面的数据，根据y='a*exp(b*x) + c*exp(d*x)拟合出来
但是 一直错误，求教！！！
错误的提示
Error using ==> eig
NaN or Inf prevents convergence.
Error in ==> fit at 436
                errstr = handleerr(es.identifier, es.message, suppresserr); return;
Error in ==> aiyo at 5
cfun=fit(x0,y0,f)

林偏殇

求救啊！！！！

zhouyang664

感觉像是矩阵自身的性质不行，从错误提示来看，不收敛吧