24小时热门版块排行榜    

查看: 1081  |  回复: 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的回帖

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的回帖
查看全部 7 个回答

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的回帖

lipenggg

专家顾问 (正式写手)

一团和气
7楼2016-08-18 09:45:37
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +4 xvJS1BYDxvBS 2026-10-31 10/500 2026-11-01 14:50 by BCdvBAEpt3zm
[论文投稿] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 xvJS1BYDxvBS 2026-10-31 7/350 2026-11-01 14:50 by BCdvBAEpt3zm
[找工作] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +4 1uaZt2kAZj2S 2026-10-31 8/400 2026-11-01 14:48 by BCdvBAEpt3zm
[考研] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +4 1uaZt2kAZj2S 2026-10-30 13/650 2026-11-01 14:37 by BCdvBAEpt3zm
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +5 X5bHxOULxFqb 2026-10-30 13/650 2026-11-01 14:16 by BCdvBAEpt3zm
[公派出国] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 1uaZt2kAZj2S 2026-10-30 7/350 2026-11-01 12:27 by BCdvBAEpt3zm
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +4 1uaZt2kAZj2S 2026-10-30 17/850 2026-11-01 12:24 by BCdvBAEpt3zm
[教师之家] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +5 X5bHxOULxFqb 2026-10-30 15/750 2026-11-01 12:11 by BCdvBAEpt3zm
[硕博家园] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +4 X5bHxOULxFqb 2026-10-30 18/900 2026-11-01 12:11 by BCdvBAEpt3zm
[公派出国] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 hKe7U1kad5Sg 2026-10-31 5/250 2026-11-01 11:14 by BCdvBAEpt3zm
[考博] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 xvJS1BYDxvBS 2026-10-31 4/200 2026-11-01 11:01 by BCdvBAEpt3zm
[硕博家园] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +4 xvJS1BYDxvBS 2026-10-31 9/450 2026-11-01 10:53 by BCdvBAEpt3zm
[基金申请] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 1uaZt2kAZj2S 2026-10-30 14/700 2026-11-01 10:43 by BCdvBAEpt3zm
[考博] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +4 1uaZt2kAZj2S 2026-10-30 16/800 2026-11-01 10:27 by BCdvBAEpt3zm
[论文投稿] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +5 X5bHxOULxFqb 2026-10-30 26/1300 2026-11-01 10:22 by BCdvBAEpt3zm
[基金申请] 售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
[考博] 2027博士申请-临床检验诊断学、病原生物学、微生物方向 50+3 yyh@ 2026-10-28 4/200 2026-10-31 12:47 by 可淡不可忘
[硕博家园] 申博自荐 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 北京莱茵编辑
信息提示
请填处理意见