求助：一自变量对应多个变量的函数拟合该怎么做啊？

albb

RT，现有一组数据x,y1,y2，请改该用什么方法拟合啊？

shamohu

用1stopt试试，很容易搞定。

albb

shamohu 发表于 2013-3-14 10:04
用1stopt试试，很容易搞定。

多谢版主！我用了下，但是效果不是很好，
我用的是1stopt1.0，sharedmodel，不知道跟版本有没有关系？

shamohu

1.0版好像根本就不支持SharedModel命令，试一下高版本。

albb

shamohu 发表于 2013-3-27 08:44
1.0版好像根本就不支持SharedModel命令，试一下高版本。

网上的那个1.5的版本我也试了一下，您看一下我的程序

Title "Type your title here";
Parameters g0,g1,g2,g3,g4,g5,t1=[0.01,0.05],t2=[0.1,0.5],t3=[1,5],t4=[10,20],t5=[90,100];
Variable x,y1,y2;
SharedModel;
Function  y2=g1*x*t1/(1+(x*t1)^2)+g2*x*t2/(1+(x*t2)^2)+g3*x*t3/(1+(x*t3)^2)+g4*x*t4/(1+(x*t4)^2)+g5*x*t5/(1+(x*t5)^2);
y1=g0+g1*(x*t1)^2/(1+(x*t1)^2)+g2*(x*t2)^2/(1+(x*t2)^2)+g3*(x*t3)^2/(1+(x*t3)^2)+g4*(x*t4)^2/(1+(x*t4)^2)+g5*(x*t5)^2/(1+(x*t5)^2);


Data;
1.00E-02        3.22E+03        1.21E+04
1.29E-02        5.09E+03        1.45E+04
1.65E-02        7.23E+03        1.72E+04
2.12E-02        9.40E+03        2.13E+04
2.73E-02        1.23E+04        2.55E+04
3.50E-02        1.53E+04        3.06E+04
4.50E-02        2.09E+04        3.62E+04
5.78E-02        2.63E+04        4.13E+04
7.42E-02        3.28E+04        4.60E+04
9.54E-02        4.04E+04        5.22E+04
1.23E-01        4.90E+04        5.72E+04
1.58E-01        5.90E+04        6.32E+04
2.03E-01        7.19E+04        6.80E+04
2.60E-01        8.20E+04        7.43E+04
3.34E-01        9.59E+04        7.76E+04
4.29E-01        1.07E+05        8.09E+04
5.51E-01        1.21E+05        8.83E+04
7.08E-01        1.37E+05        9.06E+04
9.09E-01        1.54E+05        9.53E+04
1.17E+00        1.63E+05        1.04E+05
1.50E+00        1.75E+05        1.16E+05
1.93E+00        2.02E+05        1.14E+05
2.47E+00        2.00E+05        8.86E+04
3.18E+00        2.36E+05        1.06E+05
4.08E+00        2.61E+05        1.30E+05
5.25E+00        2.72E+05        1.54E+05
6.74E+00        2.85E+05        1.34E+05
8.66E+00        3.17E+05        1.34E+05
1.11E+01        3.19E+05        1.36E+05
1.43E+01        3.53E+05        1.55E+05
1.84E+01        3.61E+05        1.50E+05
2.36E+01        3.84E+05        1.53E+05
3.03E+01        3.95E+05        1.84E+05
3.89E+01        4.11E+05        1.97E+05
5.00E+01        4.25E+05        2.01E+05

shamohu

5.0版计算没任何问题啊：

均方差(RMSE): 8363.49977005659
残差平方和(SSE): 4896368988.26156
相关系数(R): 0.998791583684114
相关系数之平方(R^2): 0.997584627638221
决定系数(DC): 0.997584627638221
F统计(F-Statistic): 775.071431079709

参数                  最佳估算
--------------------        -------------
g0        -796.487145792553
g1        486617.223767348
g2        170138.786741471
g3        117258.651634488
g4        67167.9648846951
g5        7990.15466824903
t1        0.0100000307298376
t2        0.172817920505244
t3        1.40251547105646
t4        10.0000306431163
t5        99.3860627372653

albb

shamohu 发表于 2013-4-6 11:12
5.0版计算没任何问题啊：

均方差(RMSE): 8363.49977005659

非常感谢！
看来就是版本问题，我用1.0拟合之后结果根本没有g0，估计sharedmodel没有被调用？