1stopt7.0跑出来的结果,有问题站内我。请采纳。
==================
模型公式: y1 = (50/(a-k1))*(a*(1-exp(0-k1*t))-k1*(1-exp(0-a*t)))
y1 = (50/(0.0084033309045955-322363191667758))*(0.0084033309045955*(1-exp(0-322363191667758*t))-322363191667758*(1-exp(0-0.0084033309045955*t)))
模型公式: y2 = 50*k1*(k3-a)*50/(k3*(a-k1))*(exp(0-k1*t)-exp(0-a*t))
y2 = 50*322363191667758*1.064782396E-5*50/(0.00841397872855583*(0.0084033309045955-322363191667758))*(exp(0-322363191667758*t)-exp(0-0.0084033309045955*t))
迭代数: 60
计算用时(时:分:秒:微秒): 00:00:06:411
优化算法: 通用全局优化算法(UGO1)
计算结束原因: 达到收敛判断标准
均方差(RMSE): 2.59320635116301
残差平方和(SSE): 134.494383594244
相关系数(R): 0.949540727213403
相关系数之平方(R^2): 0.901627592636958
修正R平方(Adj. R^2): 0.85408615799387
确定系数(DC): 0.76042613935766
F统计(F-Statistic): 15.4019567906444
参数 最佳估算
-------------------- -------------
a 0.0084033309045955
k1 322363191667758
k3 0.00841397872855583
====== 结果输出 ======
文件: 数据文件-1
No 目标 y1 计算 y1 目标 y2 计算 y2
1 0 0 0 0
2 5.177 1.24474372221455 3.205 3.0849696839899
3 8.335 2.85633276190374 2.953 2.9829970207315
4 10.575 4.7961149954351 2.764 2.86025806208689
5 11.441 9.13217560988145 2.212 2.58589553043911
6 14.645 11.1416462946358 2.445 2.45874706242545
7 15.705 13.0523113176466 2.415 2.33785048383576
8 15.677 16.5964185923686 2.167 2.11359848858356
9 17.529 19.8005669461762 2.144 1.91085726047507
10 17.725 22.6973660206378 1.876 1.72756343725306,