24小时热门版块排行榜    

查看: 1080  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 1uaZt2kAZj2S 2026-10-30 16/800 2026-11-01 07:01 by YXM1FDCSE8Bn
[硕博家园] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +3 X5bHxOULxFqb 2026-10-30 17/850 2026-11-01 06:48 by YXM1FDCSE8Bn
[基金申请] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +5 X5bHxOULxFqb 2026-10-30 21/1050 2026-11-01 06:46 by YXM1FDCSE8Bn
[基金申请] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 xvJS1BYDxvBS 2026-10-31 8/400 2026-11-01 05:48 by YXM1FDCSE8Bn
[硕博家园] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +3 xvJS1BYDxvBS 2026-10-31 8/400 2026-11-01 05:43 by YXM1FDCSE8Bn
[找工作] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 1uaZt2kAZj2S 2026-10-31 6/300 2026-11-01 05:41 by YXM1FDCSE8Bn
[公派出国] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 1uaZt2kAZj2S 2026-10-31 6/300 2026-11-01 05:41 by YXM1FDCSE8Bn
[博后之家] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 1uaZt2kAZj2S 2026-10-30 15/750 2026-11-01 05:22 by YXM1FDCSE8Bn
[论文投稿] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 1uaZt2kAZj2S 2026-10-31 8/400 2026-11-01 03:59 by YXM1FDCSE8Bn
[找工作] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 1uaZt2kAZj2S 2026-10-31 6/300 2026-11-01 03:56 by YXM1FDCSE8Bn
[论文投稿] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +4 X5bHxOULxFqb 2026-10-30 25/1250 2026-11-01 03:30 by YXM1FDCSE8Bn
[考博] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 1uaZt2kAZj2S 2026-10-31 5/250 2026-11-01 02:25 by YXM1FDCSE8Bn
[考研] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 1uaZt2kAZj2S 2026-10-30 12/600 2026-11-01 02:19 by YXM1FDCSE8Bn
[公派出国] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 1uaZt2kAZj2S 2026-10-31 4/200 2026-11-01 00:41 by YXM1FDCSE8Bn
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +4 X5bHxOULxFqb 2026-10-30 12/600 2026-11-01 00:15 by YXM1FDCSE8Bn
[教师之家] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 X5bHxOULxFqb 2026-10-30 13/650 2026-10-31 17:12 by hKe7U1kad5Sg
[硕博家园] 越来越担心硕士毕业后没有竞争力了 +8 zuesli 2026-10-27 14/700 2026-10-30 22:44 by 好运连连公主
[基金申请] 论文的挂名怎么评判贡献?基金申请中,排名第一的挂名文章 +3 biily 2026-10-25 3/150 2026-10-30 18:33 by 旁观2020
[硕博家园] 申博自荐 2篇一区TOP+1篇EI+2公开发明专利 2027博士申请 +3 小小阿光 2026-10-27 3/150 2026-10-29 23:56 by 起风了的快乐
[考博] 申博求助 20+3 Katrinae 2026-10-26 4/200 2026-10-27 09:56 by 北京莱茵编辑
信息提示
请填处理意见