24小时热门版块排行榜    

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

月只蓝

主管区长 (职业作家)

【答案】应助回帖

★ ★ ★ ★ ★
感谢参与,应助指数 +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的回帖

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的回帖
相关版块跳转 我要订阅楼主 bg0east 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 我不理解! +13 Edward_pc 2026-08-26 19/950 2026-08-26 19:46 by jurkat.1640
[基金申请] 2026年的国家社科基金项目通讯评审的新规则与新动向、新挑战 +7 process2012 2026-08-23 10/500 2026-08-26 19:23 by hmhminy
[基金申请] 国自然面上复盘~欢迎讨论 +3 晴天加油 2026-08-26 4/200 2026-08-26 17:59 by Tide man
[基金申请] 系统查不到 +10 董八千 2026-08-26 10/500 2026-08-26 16:30 by Equinoxhua
[基金申请] 怎么查啊 +5 huang1991js 2026-08-26 5/250 2026-08-26 16:09 by Equinoxhua
[基金申请] 范进中举一文的中心思想 +9 炎黄贵胄 2026-08-22 10/500 2026-08-26 15:40 by semaglutide
[基金申请] 能否退出参与的面上项目解除限项 +23 koalala 2026-08-24 26/1300 2026-08-26 14:29 by 宝贝虫子
[基金申请] 出来了 +9 trojank 2026-08-26 9/450 2026-08-26 14:25 by 宝贝虫子
[基金申请] 为什么国自然不能直接公布 +4 bjdxyxy 2026-08-26 4/200 2026-08-26 13:12 by qingmu1201
[基金申请] 科研孤儿太难了 +21 我4大白菜 2026-08-20 22/1100 2026-08-26 11:29 by nlgza
[基金申请] 项目信息和经费信息在系统里都可以看到了 +6 wittyboy 2026-08-26 14/700 2026-08-26 10:55 by wittyboy
[基金申请] 怎么看青基中了没有啊 +4 叶九微 2026-08-26 4/200 2026-08-26 10:52 by xiacongxu
[基金申请] 今天务委会开完了,明天出结果吗 +19 angus9576 2026-08-25 23/1150 2026-08-26 10:03 by zp519
[基金申请] 没有任何消息-是不是就凉了 +9 图啦图啦 2026-08-24 10/500 2026-08-25 11:59 by 南海小哥
[基金申请] 人气不行了 +11 fansofjerry 2026-08-21 11/550 2026-08-25 11:04 by 孤独的英雄6
[基金申请] 建议基金发布提前给出明确的时间点 +13 kulium 2026-08-21 16/800 2026-08-24 16:27 by superceng
[基金申请] 让我中一个面上吧! +13 大萍1987 2026-08-20 16/800 2026-08-24 10:23 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 紫杉醇
[基金申请] 应该是下周三26日公布了吧? +4 哈哈蛤? 2026-08-21 4/200 2026-08-21 10:58 by Vivilian
信息提示
请填处理意见