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

mello

程序如下：
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

程序如下：
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

我是初学者 不懂 请见谅啊

dingd

你的代码等同于下：

1. 
   
2. Title "myfunction";
   
3. Parameter k1，n;
   
4. Variable x,y;
   
5. Function y=(exp(k1*x))^(1/n);
   
6. data;
   
7. //x,y=
   
8. 0 0.96020507
   
9. 25 0.835247809
   
10. 50 0.531616364
    
11. 75 0.279758509
    
12. 100 0.124089619
    
13. 180 0.017483159
    
14. 300 0
    
15. 460 0
    
16. 600 0
    

复制代码

均方差(RMSE): 0.0657333451876741
残差平方和(SSE): 0.0388878540260573
相关系数(R): 0.983670555487176
相关系数之平方(R^2): 0.967607761732448

参数        最佳估算
----------        -------------
k1        -7.44034606813273
n        506.002270003345