求助：如何确定自定义函数拟合时的参数初值？感谢！

[azure_blue]

我做一个人口增长的拟合，在确定参数初值时遇到了严重的麻烦，求高手给与帮助！！！万分感谢！！！
以下是我写的程序：

>> t=[1790:10:2000];

>> population=[3.9 5.3 7.2 9.6 12.9 17.1 23.2 31.4 38.6 50.2 62.9 76.0 92.0 106.5 123.2 131.7 150.7 179.3 204.0 226.5 251.4 281.4];

>> ft=inline('a(1).*a(2)./(a(1)+exp(-a(3).*t).*a(2)-exp(-a(3).*t).*a(1))','a','t');

>> a=nlinfit(t,population,ft,[1 1 1])

Warning: Rank deficient, rank = 1, tol = 2.6167e-014.
> In nlinfit>LMfit at 296
In nlinfit at 168

Warning: Rank deficient, rank = 1, tol = 2.6050e-014.
> In nlinfit>LMfit at 296
In nlinfit at 168

Warning: Rank deficient, rank = 1, tol = 2.6038e-014.
> In nlinfit>LMfit at 296
In nlinfit at 168

Warning: The Jacobian at the solution is ill-conditioned, and some
model parameters may not be estimated well (they are not identifiable).
Use caution in making predictions.
> In nlinfit at 225
a =
1.0000 94.7727 1.0000

当然[1 1 1]是我瞎写的，显然最后的结果不对啊！但到底应该是什么呢？怎么确定？:handshake

[slgu]

a=[1.391057891*10^(-11), 44657285.16, 0.02154729647]

[azure_blue]

哥，请问你是如何算出的？能给与解释吗？非常感谢啊！！！

[slgu]

本例的非线性拟合有较大难度，原因在于对初值的要求太高。随机初值与最优值之间相差较大，用matlab自带的拟合命令(lsqcurvefit, nlinfit) 很难实现全局最优拟合。本人自编了一个小程序，能在任意初值条件下稳定实现一般曲线曲面的拟合问题。

[azure_blue]

功力深厚啊！能不能将您那个程序发给我用用呢？我邮箱[email protected] QQ是349041471，如果可以的话非常感谢！！:)

[qq437544897]

本例的非线性拟合有较大难度，原因在于对初值的要求太高。随机初值与最优值之间相差较大，用matlab自带的拟合命令(lsqcurvefit, nlinfit) 很难实现全局最优拟合。本人自编了一个小程序，能在任意初值条件下稳定实现一般曲线曲面的拟合问题。


slgu,您好。我也遇到一个对初值非常敏感的非线性拟合问题，麻烦您也传我一份您的程序吧！还请您出手相救啊！！ [email protected] 或者 [email protected] 叩谢！

[liaol09]

您好！ 刚进实验室，导师让我编个程序对数据进行高斯拟合，也遇到了这种问题。能把您的程序发给我试试么？导师交的第一个任务就不完成以后再实验室都混不下去了。。。急求，谢谢啦!!我的邮箱是[email protected]

[slgu]

可将数据给出，或以文件的形式置上，并将可能的模型给出，否则我代为寻找模型需花费较长时间。我将尽力给出合适的结果。很抱歉，源程序尚不想公示。