1stopt 数据分析

1274

收藏 2015-03-12

悬赏 50 个论坛币未解决

以1stopt sharedmodel 拟合方程组得不出结果，请各位大神看看代码哪里有错，如果方便的话再帮忙跑一遍程序
Parameters c1,c2,c3,c4,c5,c6,c7,c8,c9,c10;
SharedModel;
Variable t,total,one,two,three,four,five;
function
one=(1-exp(a1))*(-a2/a1)*exp(-a1*(t-1));
two=(1-exp(a3))*(-a4/a3)*exp(-a3*(t-1));
three=(1-exp(a5))*(-a6/a5)*exp(-a5*(t-1));
four=(1-exp(a7))*(-a8/a7)*exp(-a7*(t-1));
five=(1-exp(a9))*(-a10/a9)*exp(-a9*(t-1));
total=(1-exp(a1))*(-a2/a1)*exp(-a1*(t-1))+(1-exp(a3))*(-a4/a3)*exp(-a3*(t-1))+(1-exp(a5))*(-a6/a5)*exp(-a5*(t-1))+(1-exp(a7))*(-a8/a7)*exp(-a7*(t-1))+(1-exp(a9))*(-a10/a9)*exp(-a9*(t-1));
data;
1 45235 11761 22654 4578 3489 2753
2 47138 8922 24193 6120 4443 3460
3 55625 9045 25012 9604 8549 3415
4 53360 7714 20191 13303 9742 4410
5 56015 6610 15503 14736 13895 5271

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

全部回复

gooxin

2016-8-27 20:03:57

模型公式: one = (1-exp(a1))*(-a2/a1)*exp(0-a1*(t-1))
      one = (1-exp(0.098742301576816))*(-9987.68471207633/0.098742301576816)*exp(0-0.098742301576816*(t-1))
模型公式: two = (1-exp(a3))*(-a4/a3)*exp(0-a3*(t-1))
      two = (1-exp(0.0678622270641146))*(-23581.2539541204/0.0678622270641146)*exp(0-0.0678622270641146*(t-1))
模型公式: three = (1-exp(a5))*(-a6/a5)*exp(0-a5*(t-1))
      three = 0.2158367758*26020.2124*exp(0-(-0.243138086123822)*(t-1))
模型公式: four = (1-exp(a7))*(-a8/a7)*exp(0-a7*(t-1))
      four = 0.2436871109*17503.14212*exp(0-(-0.279300113992489)*(t-1))
模型公式: five = (1-exp(a9))*(-a10/a9)*exp(0-a9*(t-1))
      five = 0.1219824744*23402.26073*exp(0-(-0.130088724722968)*(t-1))
模型公式: total = (1-exp(a1))*(-a2/a1)*exp(0-a1*(t-1))+(1-exp(a3))*(-a4/a3)*exp(0-a3*(t-1))+(1-exp(a5))*(-a6/a5)*exp(0-a5*(t-1))+(1-exp(a7))*(-a8/a7)*exp(0-a7*(t-1))+(1-exp(a9))*(-a10/a9)*exp(0-a9*(t-1))
      total = (1-exp(0.098742301576816))*(-9987.68471207633/0.098742301576816)*exp(0-0.098742301576816*(t-1))+(1-exp(0.0678622270641146))*(-23581.2539541204/0.0678622270641146)*exp(0-0.0678622270641146*(t-1))+0.2158367758*26020.2124*exp(0-(-0.243138086123822)*(t-1))+0.2436871109*17503.14212*exp(0-(-0.279300113992489)*(t-1))+0.1219824744*23402.26073*exp(0-(-0.130088724722968)*(t-1))

迭代数: 53
计算用时(时:分:秒:微秒): 00:00:07:278
优化算法: 通用全局优化算法(UGO1)
计算结束原因: 达到收敛判断标准
均方差(RMSE): 1656.79523577114
残差平方和(SSE): 82349113.5982186
相关系数(R): 0.897426042984151
相关系数之平方(R^2): 0.805373502626191
修正R平方(Adj. R^2): 1.24454267368637
确定系数(DC): 0.784743980423066
F统计(F-Statistic): -1.45128662992522

参数       最佳估算
-------------------- -------------
a1 0.098742301576816
a2 9987.68471207633
a3 0.0678622270641146
a4 23581.2539541204
a5 -0.243138086123822
a6 6326.50464236963
a7 -0.279300113992489
a8 4888.62958875922
a9 -0.130088724722968
a10 3044.37025394851
c1 1.5672286208719
c2 1.430365646258
c3 -2.63553559267893
c4 1.66500383219682
c5 1.7622335809283
c6 -3.15033708768897
c7 -0.178150566294789
c8 0.50601847236976
c9 0
c10 -1.55967446928844

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

文件: 数据文件-1
No 目标 total 计算 total 目标 one 计算 one 目标 two 计算 two 目标 three 计算 three 目标 four 计算 four 目标 five 计算 five
1 45235 47633.3046703029 11761 10497.4269503682 22654 24399.8031665979 4578 5616.11874815824 3489 4265.29013489372 2753 2854.66567028479
2 47138 48362.1034681167 8922 9510.41841737679 24193 22798.9125183617 6120 7161.92569938108 4443 5639.58408818472 3460 3251.26274481234
3 55625 50212.1176564183 9045 8616.21223002719 25012 21303.057589066 9604 9133.20782974476 8549 7456.68118272072 3415 3702.95882485955
4 53360 53435.1677731245 7714 7806.08275417466 20191 19905.346900979 13303 11647.0749296547 9742 9859.25440445707 4410 4217.40878385905
5 56015 58363.6140867622 6610 7072.1247734204 15503 18599.3411317482 14736 14852.8706392945 13895 13035.9465598527 5271 4803.33098244646

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

扫码加我 拉你入群

扫码加我 拉你入群

分享

扫码加好友，拉您进群

扫码加我拉你入群

扫码加我拉你入群