声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

查看: 2704|回复: 3

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

[复制链接]
发表于 2011-1-8 19:24 | 显示全部楼层 |阅读模式

马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。

您需要 登录 才可以下载或查看,没有账号?我要加入

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的值相差很大,不知道取哪个值?才是正确的
请高手指点!
回复
分享到:

使用道具 举报

 楼主| 发表于 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的值相差很大,不知道取哪个值?才是正确的
请高手指点!
 楼主| 发表于 2011-1-8 19:25 | 显示全部楼层
我是初学者 不懂 请见谅啊
发表于 2011-1-9 13:51 | 显示全部楼层
你的代码等同于下:

  1. Title "myfunction";
  2. Parameter k1,n;
  3. Variable x,y;
  4. Function y=(exp(k1*x))^(1/n);
  5. data;
  6. //x,y=
  7. 0 0.96020507
  8. 25 0.835247809
  9. 50 0.531616364
  10. 75 0.279758509
  11. 100 0.124089619
  12. 180 0.017483159
  13. 300 0
  14. 460 0
  15. 600 0
复制代码
均方差(RMSE): 0.0657333451876741
残差平方和(SSE): 0.0388878540260573
相关系数(R): 0.983670555487176
相关系数之平方(R^2): 0.967607761732448

参数        最佳估算
----------        -------------
k1        -7.44034606813273
n        506.002270003345
您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

QQ|小黑屋|Archiver|手机版|联系我们|声振论坛

GMT+8, 2024-12-23 11:49 , Processed in 0.065339 second(s), 17 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

快速回复 返回顶部 返回列表