24小时热门版块排行榜    

查看: 970  |  回复: 4

东景辉

新虫 (正式写手)

[求助] 1stopt微分方程组参数拟合(各位大哥求帮忙,没有高版本的软件) 已有1人参与

初始金币 5
Parameter  a11=0.01[0,10000];
Parameter  a12=0.01[0,10000];
Parameter  a13=0.01[0,10000];
Parameter  b1=0.01[0,10000];
Parameter  a21=0.01[0,10000];
Parameter  a22=0.01[0,10000];
Parameter  a23=0.01[0,10000];
Parameter  b2=0.01[0,10000];
Parameter  a31=0.01[0,10000];
Parameter  a32=0.01[0,10000];
Parameter  a33=0.01[0,10000];
Parameter  b3=0.01[0,10000];
Variable t,x,y;
ODEFunction x'=a11* x+a12*y+a13*z-b1*x*x;
         y'= a21* x+a22*y+a23*z-b2*y*y;
z'= a31* x+a32*y+a33*z-b2*y*y;
Data;
1     0  0      0
33    15 12     6
43    39 30     30
72   84   60     63
92    174  102    119
102   334  194    188
109   414  258    284
123  608  452    404 @beefly
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dingd

铁杆木虫 (职业作家)

【答案】应助回帖

感谢参与,应助指数 +1
内容已删除
2楼2017-05-11 13:48:49
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

东景辉

新虫 (正式写手)

非常感谢!!参数a11为0有点问题啊。 @beefly

发自小木虫Android客户端
3楼2017-05-11 15:16:45
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

东景辉

新虫 (正式写手)

能否再算一下Parameter??a11=0.01[0,10000];
Parameter??a12=0.01[0,10000];
Parameter??b1=0.01[0,10000];
Parameter??a21=0.01[0,10000];
Parameter??a22=0.01[0,10000];
Parameter??b2=0.01[0,10000];
Variable t,x,y;
ODEFunction x'=a11* x+a12*y-b1*x*x);
? ?? ?? ?y'= a21* x+a22*y-b2*y*y);
Data;
1? ???0? ?? ?0
33? ? 12? ???6
43? ? 30? ???30
72? ? 60? ???63
92? ? 102? ? 119
102? ?194? ? 188
109? ?258? ? 284
123? ?452? ? 404 @beefly

发自小木虫Android客户端
4楼2017-05-11 15:18:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

东景辉

新虫 (正式写手)

引用回帖:
2楼: Originally posted by dingd at 2017-05-11 13:48:49
参考下:
均方差(RMSE): 22.4162459568179
残差平方和(SSE): 7034.83315915174
相关系数(R): 0.99124454515698
相关系数之平方(R^2): 0.982565748303469
修正R平方(Adj. R^2): 0.965161343617583
确定系数(DC ...

非常感谢楼主,不过参数a11理论上不能是0,能否再重新算一下,谢谢
Parameter a11=0.01[0,10000];
Parameter a12=0.01[0,10000];
Parameter b1=0.01[0,10000];
Parameter a21=0.01[0,10000];
Parameter a22=0.01[0,10000];
Parameter b2=0.01[0,10000];
Variable t,x,y;
ODEFunction x'=a11* x+a12*y-b1*x*x);
y'= a21* x+a22*y-b2*y*y);
Data;
1 0 0
33 12 6
43 30 30
72 60 63
92 102  119
102 194 188
109 258  284
123 452 404 @beefly
5楼2017-05-11 15:32:08
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 东景辉 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 朋友圈看到的 +4 wangzilk 2026-08-18 5/250 2026-08-18 19:25 by wangzilk
[基金申请] 欢迎发来filecode的Mz6后的代码验证其规律 +51 医学老男孩 2026-08-13 110/5500 2026-08-18 19:07 by icerianxia
[基金申请] 明天放榜? +5 Shxjjxjkx 2026-08-18 5/250 2026-08-18 18:14 by -大大大大大-
[基金申请] 时间戳变了,能看出什么问题? +8 基诺咪客 2026-08-17 11/550 2026-08-18 17:13 by 悠然熊猫
[基金申请] 重要消息,中午系统在维护 +9 yuleib84 2026-08-18 10/500 2026-08-18 16:52 by archvillain
[教师之家] 两姐妹,一个辛辛苦苦30岁工作当大学教师,另一个23岁就工作当开车教练 +3 瞬息宇宙 2026-08-11 7/350 2026-08-18 16:42 by LNP@mRNA
[论文投稿] 投稿咨询 +4 wwm09 2026-08-17 6/300 2026-08-18 15:36 by wwm09
[基金申请] 时间戳又变了8-15 +14 archvillain 2026-08-15 26/1300 2026-08-18 13:37 by phantomgost
[基金申请] 今天维护系统维护 祝所有人 高中 +8 gjjjzhong 2026-08-18 9/450 2026-08-18 13:01 by 家与远方
[基金申请] 哪位老哥知道今年的国自然具体哪一天放榜? +13 Ldrop2023 2026-08-13 16/800 2026-08-18 12:25 by 淀粉搬运工
[基金申请] 2027广东省杰青 +3 奶牛小黑 2026-08-15 8/400 2026-08-18 11:47 by gltch
[基金申请] 93BebMhtakh前后11位开头都是大写 +4 且听虎啸 2026-08-17 5/250 2026-08-18 00:49 by 蔡棒棒菂
[基金申请] 感觉是下周放榜了 +6 angus9576 2026-08-17 11/550 2026-08-17 23:57 by angus9576
[基金申请] 今天系统多次维护,明天很可能放榜! +8 zju2000 2026-08-16 9/450 2026-08-17 12:20 by lmz0216
[精细化工] 招聘 金属平磨液,抛光液研发工程师 +3 小天0311 2026-08-14 3/150 2026-08-16 07:31 by H9PLUS
[基金申请] 各位道友,我要去昆明玩几天,回来见。 +7 Tide man 2026-08-14 8/400 2026-08-15 01:11 by arzu_hma
[基金申请] 重要来源:本周末出结果 +10 瞬息宇宙 2026-08-12 10/500 2026-08-13 15:46 by likettle
[基金申请] 不应该看fileCode +7 且听虎啸 2026-08-12 9/450 2026-08-13 14:27 by flydreamws
[基金申请] Filecode 又变了,巨变 +3 WH3796 2026-08-12 4/200 2026-08-13 14:13 by 小木虫6752397
[基金申请] 结合人工智能,周易传统文化,filecode打分制来了,3分以上希望很大。 +3 Tide man 2026-08-12 4/200 2026-08-13 08:35 by ZJTJZ
信息提示
请填处理意见