24小时热门版块排行榜    

查看: 861  |  回复: 6

zyugiec

金虫 (正式写手)

[求助] 关于共参数的多因变量拟合 已有1人参与

大家好,以下程序是用1stopt写的,用到了sharemodel模块,但由于是盗版非正版,无法拟合,请问有其他软件或其他办法吗?
盗版的1stopt无法用sharemodel模块。求赐教,非常感激!

两个因变量y1和y2是共用参数k1, k2和k3的。

Parameters a[0,],k1[0,],k3[0,];
Variable t,y1,y2;
SharedModel;
Function y1=(50/(a-k1))*(a*(1-exp(-k1*t))-k1*(1-exp(-a*t))); y2=50*k1*(k3-a)*50/(k3*(a-k1))*(exp(-k1*t)-exp(-a*t));
Data;
0                0                0
3                5.177        3.205
7                8.335        2.953
12                10.575        2.764
24                11.441        2.212
30                14.645        2.445
36                15.705        2.415
48                15.677        2.167
60                17.529        2.144
72                17.725        1.876
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lipenggg

专家顾问 (正式写手)

【答案】应助回帖

感谢参与,应助指数 +1
迭代数: 40
计算用时(时:分:秒:微秒): 00:00:03:237
优化算法: 麦夸特法(Levenberg-Marquardt) + 通用全局优化法
计算结束原因: 达到收敛判定标准
均方差(RMSE): 2.59320636093186
残差平方和(SSE): 134.494384607549
相关系数(R): 0.960355233351647
相关系数之平方(R^2): 0.922282174225897
决定系数(DC): 0.921963506716319
F统计(F-Statistic): 94.9364874829745

参数                  最佳估算
--------------------        -------------
a                  0.00840333088924141
k1                 13346110.5010209
k3                 0.0084139787133626

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

No        实测y1        计算y1        实测y2        计算y2
1        0        0        0        0
2        5.177        1.2447436892702        3.205        3.08496973825087
3        8.335        2.85633272715297        2.953        2.9829970733821
4        10.575        4.79611495864386        2.764        2.8602581127907
5        11.441        9.13217556908946        2.212        2.58589557675574
6        14.645        11.1416462522698        2.445        2.45874710669121
7        15.705        13.0523112739599        2.415        2.33785052614034
8        15.677        16.5964185467179        2.167        2.11359852721962
9        17.529        19.8005668993401        2.144        1.91085729575714
10        17.725        22.6973659732639        1.876        1.72756346946909
关于共参数的多因变量拟合
y1.png


关于共参数的多因变量拟合-1
y2.png

» 本帖已获得的红花(最新10朵)

一团和气
2楼2016-08-17 13:59:38
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zyugiec

金虫 (正式写手)

送红花一朵
内容已删除
3楼2016-08-17 14:18:27
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zyugiec

金虫 (正式写手)

送红花一朵
引用回帖:
2楼: Originally posted by lipenggg at 2016-08-17 13:59:38
迭代数: 40
计算用时(时:分:秒:微秒): 00:00:03:237
优化算法: 麦夸特法(Levenberg-Marquardt) + 通用全局优化法
计算结束原因: 达到收敛判定标准
均方差(RMSE): 2.59320636093186
残差平方和(SSE): 134.49438 ...

能否把代码发给我一下
4楼2016-08-17 14:31:19
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lipenggg

专家顾问 (正式写手)

内容已删除
一团和气
5楼2016-08-18 07:32:53
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zyugiec

金虫 (正式写手)

内容已删除
6楼2016-08-18 09:35:28
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lipenggg

专家顾问 (正式写手)

一团和气
7楼2016-08-18 09:45:37
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 zyugiec 的主题更新
信息提示
请填处理意见