24小时热门版块排行榜    

查看: 811  |  回复: 2
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

bg0east

新虫 (初入文坛)

[求助] 1stopt 微分方程数据拟合求助 已有1人参与

代码如下:
Parameter b1, b2, b3, b4;
Variable x, y;
ODEFunction y'=0.001*b3-(y/b1)^b2 *b3-(y/b1)^(b2-1)*b4/b1*b3
DataFile "CodeSheet1[A1:B27]"; //Press F2 to view data

数据见附件,多谢先!
回复此楼

» 本帖附件资源列表

  • 欢迎监督和反馈:小木虫仅提供交流平台,不对该内容负责。
    本内容由用户自主发布,如果其内容涉及到知识产权问题,其责任在于用户本人,如对版权有异议,请联系邮箱:xiaomuchong@tal.com
  • 附件 1 : data1.xls
  • 2015-10-07 19:09:29, 9 K

» 猜你喜欢

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

bg0east

新虫 (初入文坛)

引用回帖:
2楼: Originally posted by 月只蓝 at 2015-10-07 21:52:20
Root of Mean Square Error (RMSE): 8.08398907127184
Residual Sum of Squares (SSE): 2091.22813774216
Correlation Coef. (R): 0.998652058700909
R-Square: 0.997305934347565
Determination Coef. (DC): 0 ...

非常感谢!能否用同样的数据再次代跑一下这个代码?非线性微分方程的参数拟合,用matlab做初值貌似很难。

Parameter c1, c2, c3, N;
Variable x, y;
ODEFunction y'=c1-c2*y^(N-1)-c3*y^N
DataFile "CodeSheet1[A1:B32]"; //Press F2 to view data
3楼2015-10-18 16:03:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 3 个回答

月只蓝

主管区长 (职业作家)

【答案】应助回帖

★ ★ ★ ★ ★
感谢参与,应助指数 +1
bg0east: 金币+5, ★★★★★最佳答案 2015-10-18 15:44:29
Root of Mean Square Error (RMSE): 8.08398907127184
Residual Sum of Squares (SSE): 2091.22813774216
Correlation Coef. (R): 0.998652058700909
R-Square: 0.997305934347565
Determination Coef. (DC): 0.997004327699384
F-Statistic: 2908.86039638179
Parameters Best Estimate
-------------------- -------------
b1 65726.5347698722
b2 -2.86969203064131
b3 -0.00574926776698071
b4 559.843140932705
====== Output Results =====
File: Data file - 1
No Obs. y Cal. y
1 513 499.966124334547
2 558 532.063700692123
3 607 590.036631538515
4 635 636.355901922696
5 670 676.764333806462
6 695 703.106016472069
7 723 733.137193055944
8 747 755.565582488483
9 765 776.329652868411
10 793 798.195338713362
11 810 814.256131562792
12 838 835.766754254898
13 845 848.153789626046
14 873 862.612409855305
15 884 877.225267353007
16 895 891.299648257148
17 905 902.430510407337
18 916 916.479655827076
19 933 929.862016863361
20 944 940.854808625361
21 954 949.707418242262
22 969 961.668817185868
23 976 971.55123175092
24 986 981.118733642333
25 993 991.912613342909
26 1000 1000.87160792104
27 1010 1009.57894001842
28 1010 1018.05050222379
29 1030 1027.65510194984
30 1030 1035.66336683093
31 1030 1038.28899536714
MATLAB、MS小问题、普通问题请发帖求助!时间精力有限,恕不接受无偿私信求助。
2楼2015-10-07 21:52:20
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 7K1CJE38xLG4 2026-08-25 4/200 2026-08-26 02:04 by cNXvBfCpiZOM
[论文投稿] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 7K1CJE38xLG4 2026-08-25 4/200 2026-08-26 02:04 by cNXvBfCpiZOM
[基金申请] 有没有大神帮我看看基金代码 28+4 1234567wang 2026-08-24 10/500 2026-08-25 19:15 by lfy8008
[基金申请] 2026国自然函评费到账 +20 羊腰板 2026-08-21 23/1150 2026-08-25 16:34 by zsna
[基金申请] 今日不放榜?网传国自然预计 8 月 27 日可查结果 +17 医学老男孩 2026-08-20 22/1100 2026-08-25 15:36 by 医学老男孩
[基金申请] 如果此刻你正在为国基感到焦虑,不妨来听听这首《基金之外》 +8 scalable 2026-08-24 8/400 2026-08-25 12:52 by jnhyjjm
[基金申请] 人气不行了 +11 fansofjerry 2026-08-21 11/550 2026-08-25 11:04 by 孤独的英雄6
[基金申请] 2026年的国家社科基金项目通讯评审的新规则与新动向、新挑战 +5 process2012 2026-08-23 7/350 2026-08-25 09:42 by huixian257
[基金申请] 今天基金会出结果吗?20260819 +17 kkkl_v 2026-08-19 18/900 2026-08-25 09:41 by windflowerwy
[基金申请] 我面上完蛋了 +13 且听虎啸 2026-08-20 14/700 2026-08-25 09:10 by mrkang
[基金申请] 只有每年这种时候来逛逛小木虫 +25 yaoyewhu2008 2026-08-20 27/1350 2026-08-25 08:01 by Equinoxhua
[基金申请] 能否退出参与的面上项目解除限项 +21 koalala 2026-08-24 24/1200 2026-08-24 19:25 by 家与远方
[基金申请] filecode,4个jtjc了 +14 ziyangfang 2026-08-19 17/850 2026-08-24 18:37 by 哈哈蛤?
[基金申请] 建议基金发布提前给出明确的时间点 +13 kulium 2026-08-21 16/800 2026-08-24 16:27 by superceng
[教师之家] 跳槽后在研项目怎么办? +5 简单化xn 2026-08-22 10/500 2026-08-23 12:38 by 简单化xn
[基金申请] 今天放榜吗? +15 布布和一二 2026-08-19 16/800 2026-08-23 09:55 by 张春生
[基金申请] 看来今天不会放榜了? +8 chengyan1220 2026-08-21 11/550 2026-08-21 17:52 by dcqxinyang
[基金申请] 时间戳又变了 +13 wuchongjun 2026-08-20 19/950 2026-08-21 17:21 by 紫杉醇
[基金申请] 今天放榜没戏了吧 +9 yuleib84 2026-08-19 11/550 2026-08-21 10:06 by gltch
[基金申请] 基金啊基金 +4 longfie172 2026-08-20 4/200 2026-08-21 08:58 by mark mao
信息提示
请填处理意见