24小时热门版块排行榜    

查看: 1629  |  回复: 17
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

小婧子1993

新虫 (初入文坛)

[求助] 20171217——方程组求解显示函数表达式有误1stOpt 已有2人参与

求大神出手 小白感激不尽

程序如下:
Variable x,y;
Function y=(1.0906*10^(-5)*1.036744442^a*0.65^b*0.016^c*0.00207^(-c-d)*4.5^(d-1))*x^d;
         a>0;
         b>0;
         c>0;
         d>1;
Data;
   0        0
120        0.092200061
240        0.155061401
360        0.210170689
480        0.260781153
600        0.308289484
720        0.353463558
840        0.396784536
960        0.438579873
1080        0.479085569
1200        0.518479044
1320        0.556898081
1380        0.575777365
1440        0.594452478
1500        0.612934006

拟合后结果如下:
迭代数: 16
计算用时(时:分:秒:微秒): 00:00:00:80
优化算法: 麦夸特法(Levenberg-Marquardt) + 通用全局优化法
计算结束原因: 达到收敛判断标准
均方差(RMSE): 0.415014295003779
残差平方和(SSE): 2.58355297586225
相关系数(R): 0
相关系数之平方(R^2): 0
决定系数(DC): -3.89435349548653
卡方系数(Chi-Square): 5.552957298
F统计(F-Statistic): 0

参数        最佳估算
----------        -------------
d                  0.676846150308847

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

No        实测值y        计算值y
1        0        0.0000000
2        0.092200061        0.0000000
3        0.155061401        0.0000000
4        0.210170689        0.0000000
5        0.260781153        0.0000000
6        0.308289484        0.0000000
7        0.353463558        0.0000000
8        0.396784536        0.0000000
9        0.438579873        0.0000000
10        0.479085569        0.0000000
11        0.518479044        0.0000000
12        0.556898081        0.0000000
13        0.575777365        0.0000000
14        0.594452478        0.0000000
15        0.612934006        0.0000000

Q : 我的约束条件为a>0,b>0,c>0,d>1,但是结果不对呀 ,是不是我的约束条件有问题 还是程序写错了呢 如果写错了  麻烦大神帮忙修改下
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lipenggg

专家顾问 (正式写手)

【答案】应助回帖

a: 7.13388623539761E-14
b: 2.99999999999994
c: 5.5992331458212E-14
d: 1
一团和气
8楼2017-12-18 16:04:42
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 18 个回答

dingd

铁杆木虫 (职业作家)

【答案】应助回帖

感谢参与,应助指数 +1
估计是你的版本有问题吧,7.0运行如下,注意楼主的拟合方程是过拟合的,参数组值不唯一:
1:
均方差(RMSE): 0.0352143054057536
残差平方和(SSE): 0.0186007095781453
相关系数(R): 0.994135874720372
相关系数之平方(R^2): 0.988306137406039
修正R平方(Adj. R^2): 0.986357160307046
确定系数(DC): 0.964762306485339
卡方系数(Chi-Square): 0.045371714138712
F统计(F-Statistic): 132.579053952861
约束条件: a-0 = 0
          b-0 = 8.29933096981543
          c-0 = 0.529347831510118
          d-1 = 0

参数        最佳估算
----------        -------------
a        0
b        8.29933096981543
c        0.529347831510118
d        1

2:
均方差(RMSE): 0.0352143054057536
残差平方和(SSE): 0.0186007095781453
相关系数(R): 0.994135874720372
相关系数之平方(R^2): 0.98830613740604
修正R平方(Adj. R^2): 0.986357160307046
确定系数(DC): 0.964762306485339
卡方系数(Chi-Square): 0.045371714154723
F统计(F-Statistic): 132.579053933663
约束条件: a-0 = 34.7152563389921
          b-0 = 8.74470347086842
          c-0 = 0.0106014793294766
          d-1 = 0

参数        最佳估算
----------        -------------
a        34.7152563389921
b        8.74470347086842
c        0.0106014793294766
d        1

3:
均方差(RMSE): 0.0352143054057536
残差平方和(SSE): 0.0186007095781453
相关系数(R): 0.994135874720372
相关系数之平方(R^2): 0.98830613740604
修正R平方(Adj. R^2): 0.986357160307047
确定系数(DC): 0.964762306485339
卡方系数(Chi-Square): 0.0453717141184233
F统计(F-Statistic): 132.579053977187
约束条件: a-0 = 2.06640356342782
          b-0 = 6.70284020175292
          c-0 = 0.156588361284413
          d-1 = 0

参数        最佳估算
----------        -------------
a        2.06640356342782
b        6.70284020175292
c        0.156588361284413
d        1
2楼2017-12-17 22:26:49
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

小婧子1993

新虫 (初入文坛)

引用回帖:
2楼: Originally posted by dingd at 2017-12-17 22:26:49
估计是你的版本有问题吧,7.0运行如下,注意楼主的拟合方程是过拟合的,参数组值不唯一:
1:
均方差(RMSE): 0.0352143054057536
残差平方和(SSE): 0.0186007095781453
相关系数(R): 0.994135874720372
相关系 ...

请问下代码是用的我的代码吗?约束条件好像有变化?
3楼2017-12-18 09:37:02
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dingd

铁杆木虫 (职业作家)

【答案】应助回帖

是用楼主的代码,没做任何改动。多解,约束条件也都满足。
4楼2017-12-18 10:19:08
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
信息提示
请填处理意见