声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

查看: 2869|回复: 2

[1stopt] 这是曲面拟合吗 为什么每次算出的结果都不同呢

[复制链接]
发表于 2010-9-17 14:56 | 显示全部楼层 |阅读模式

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

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

x
Title 1;
Variable A,B,N;
Parameters m,n,k,c;
Function 0.014442395*c^2*A^(2*m)*B^(2*n)+1.531243101*c*A^m*B^n*N-0.9842*N^2+k=0;
data;
2.739E-02        2.543E+01        0.007297312
2.815E-02        2.539E+01        0.005969001
2.920E-02        2.537E+01        0.005838403
3.101E-02        2.531E+01        0.004894325
3.331E-02        2.527E+01        0.004772042


2.559E-02        2.768E+01        0.00735714
2.414E-02        2.769E+01        0.007010803
2.704E-02        2.764E+01        0.006601412
2.801E-02        2.763E+01        0.006364402


2.323E-02        3.117E+01        0.006169685
2.406E-02        3.114E+01        0.005547479
2.761E-02        3.107E+01        0.004825175
2.925E-02        3.102E+01        0.004372091
2.985E-02        3.100E+01        0.00414964

输出结果:
代数: 55
计算用时(时:分:秒:微秒): 00:00:15:297
优化算法: 准牛顿法(BFGS) + 通用全局优化法
计算结束原因: 达到收敛判断标准
均方差(RMSE): 0.000441861661807771
残差平方和(SSE): 2.73338419445734E-6
相关系数(R): 0.903725744395038
相关系数之平方(R^2): 0.816720221082365
决定系数(DC): 0.816720219656846
卡方系数(Chi-Square): 0.000220479940498009
F统计(F-Statistic): 53.4736712956903

参数        最佳估算
----------        -------------
m                  -0.115112154309556
n                  -0.11145076950358
k                  0.630225348291101
c                  6.12366603218399

====== 结果输出 =====

No        实测值N        计算值N
1        0.007297312        0.0068027
2        0.005969001        0.0065021
3        0.005838403        0.0061041
4        0.004894325        0.0055079
5        0.004772042        0.0048609
6        0.00735714        0.0066372
7        0.007010803        0.0073395
8        0.006601412        0.0060384
9        0.006364402        0.0056756
10        0.006169685        0.0064286
11        0.005547479        0.0060526
12        0.004825175        0.0047549
13        0.004372091        0.0043050
14        0.00414964        0.0041593
我继续点“执行” 又得到这个结果:
1

迭代数: 361
计算用时(时:分:秒:微秒): 00:01:35:797
优化算法: 准牛顿法(BFGS) + 通用全局优化法
计算结束原因: 达到收敛判断标准
均方差(RMSE): 0.000442233966328597
残差平方和(SSE): 2.73799233364612E-6
相关系数(R): 0.903554775973902
相关系数之平方(R^2): 0.816411233185247
决定系数(DC): 0.81641123318505
卡方系数(Chi-Square): 0.000220983294584277
F统计(F-Statistic): 53.3634762529257

参数        最佳估算
----------        -------------
m                  -1.43223869633129
n                  -1.37881619892022
k                  3.24859882199665
c                  -0.00192382061892242

====== 结果输出 =====

No        实测值N        计算值N
1        0.007297312        0.0067914
2        0.005969001        0.0064968
3        0.005838403        0.0061057
4        0.004894325        0.0055157
5        0.004772042        0.0048640
6        0.00735714        0.0066346
7        0.007010803        0.0073229
8        0.006601412        0.0060460
9        0.006364402        0.0056872
10        0.006169685        0.0064372
11        0.005547479        0.0060669
12        0.004825175        0.0047653
13        0.004372091        0.0042956
14        0.00414964        0.0041396

每次算出的结果都不同 我改怎么办呢
回复
分享到:

使用道具 举报

发表于 2010-9-17 22:32 | 显示全部楼层
本帖最后由 dingd 于 2010-9-17 22:33 编辑

1stOpt不区分大小写:
Variable A,B,N;
Parameters m,n,k,c;
上面都出现N,是不正确的。

  1. Title 1;
  2. Variable A,B,N1;
  3. Parameters m,n,k,c;
  4. Function 0.014442395*c^2*A^(2*m)*B^(2*n)+1.531243101*c*A^m*B^n*N1-0.9842*N1^2+k=0;
  5. data;
  6. 2.739E-02        2.543E+01        0.007297312
  7. 2.815E-02        2.539E+01        0.005969001
  8. 2.920E-02        2.537E+01        0.005838403
  9. 3.101E-02        2.531E+01        0.004894325
  10. 3.331E-02        2.527E+01        0.004772042
  11. 2.559E-02        2.768E+01        0.00735714
  12. 2.414E-02        2.769E+01        0.007010803
  13. 2.704E-02        2.764E+01        0.006601412
  14. 2.801E-02        2.763E+01        0.006364402
  15. 2.323E-02        3.117E+01        0.006169685
  16. 2.406E-02        3.114E+01        0.005547479
  17. 2.761E-02        3.107E+01        0.004825175
  18. 2.925E-02        3.102E+01        0.004372091
  19. 2.985E-02        3.100E+01        0.00414964
复制代码
1stOpt 4.0可得唯一稳定结果:

(RMSE): 2.1429602047289E-9
(RSS): 6.4291898146724E-17
(R): 0.999999999998941
(DC): 0.999999999995689

m        0.116487877315394
n        0.114643154090158
k        -3.65008026534162E-5
c        0.00816136216148614
 楼主| 发表于 2010-9-18 07:58 | 显示全部楼层
非常感谢您!您说的太对了!回复 dingd 的帖子


   
您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

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

GMT+8, 2024-12-23 11:40 , Processed in 0.063663 second(s), 18 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

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