mello 发表于 2011-1-8 19:24

我这个方程ln(y^n)=-(n+1)*k1*x应该用什么方法来拟合呢?

程序如下:
Title "myfunction";
Parameter k1,n;
Variable x,y;
Function ln(y^n)=-(n+1)*k1*x;
data;
//x,y=
0 0.96020507
25 0.835247809
50 0.531616364
75 0.279758509
100 0.124089619
180 0.017483159
300 0
460 0
600 0
我用标准简面体爬山法+通用全局优化法
结果我没戏运行出来的K1、n的值相差很大,不知道取哪个值?才是正确的
请高手指点!

mello 发表于 2011-1-8 19:25

程序如下:
Title "myfunction";
Parameter k1,n;
Variable x,y;
Function ln(y^n)=-(n+1)*k1*x;
data;
//x,y=
0 0.96020507
25 0.835247809
50 0.531616364
75 0.279758509
100 0.124089619
180 0.017483159
300 0
460 0
600 0
我用标准简面体爬山法+通用全局优化法
结果我每次运行出来的K1、n的值相差很大,不知道取哪个值?才是正确的
请高手指点!

mello 发表于 2011-1-8 19:25

我是初学者 不懂 请见谅啊

dingd 发表于 2011-1-9 13:51

你的代码等同于下:
Title "myfunction";
Parameter k1,n;
Variable x,y;
Function y=(exp(k1*x))^(1/n);
data;
//x,y=
0 0.96020507
25 0.835247809
50 0.531616364
75 0.279758509
100 0.124089619
180 0.017483159
300 0
460 0
600 0
均方差(RMSE): 0.0657333451876741
残差平方和(SSE): 0.0388878540260573
相关系数(R): 0.983670555487176
相关系数之平方(R^2): 0.967607761732448

参数        最佳估算
----------        -------------
k1        -7.44034606813273
n        506.002270003345
页: [1]
查看完整版本: 我这个方程ln(y^n)=-(n+1)*k1*x应该用什么方法来拟合呢?