| 查看: 1439 | 回复: 3 | |||
[交流]
【求助】matlab微分方程参数拟合问题 已有2人参与
|
|
之前有一个帖子是这样的： 本人小硕士一枚，到了毕业关键的时候，写论文发现需要对数据进行参数拟合，看类似文献上使用matlab进行参数拟合的，于是开始借相关书籍自学，但是学了一段时间发现了一个结果，越学越难，再这样下去论文就完不成了，这个问题又无法绕开，所以只能在网上求助各位大哥了。 动力学模型：dc/dt=(-50000*k*(125*c-1)^g)*(3-66.67*c)^0.667 需要拟合的参数：k，g 实验数据：t(min): 0 100 200 300 400 c(kg/kg): 0.030 0.0224 0.0149 0.0124 0.0111 希望能有贵人相助。详见 http://muchong.com/t-6727903-1 我的问题是： 请问一下这个例子中的待拟合参数k或g其中一个如果改变成与c一样也是一组采样值，也就是数据为三组t,c,g（以g为例）， 实验数据：t(min): 0 100 200 300 400 c(kg/kg): 0.030 0.0224 0.0149 0.0124 0.0111 g: 1 2 3 4 5 也就是待辨识的微分方程里面不仅有t,c还需要有g这个变量，那么如何编写这种问题的matlab程序呢？我遇到了这个问题，目前束手无策，想请哪位大神回复一下。多谢！ |
回复此楼» 猜你喜欢
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有5人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有7人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有6人回复
-
售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急
已经有8人回复
-
上海工程技术大学激光智能制造课题组｜2027级博士研究生招生公告
已经有8人回复
-
上海工程技术大学激光智能制造课题组招收博士研究生
已经有8人回复
-
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有3人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有3人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有4人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急
已经有3人回复
chendequan
铁虫 (小有名气)
- 应助: 152 (高中生)
- 金币: 918.6
- 红花: 17
- 帖子: 262
- 在线: 209小时
- 虫号: 994474
- 注册: 2010-04-11
- 专业: 结构工程
2楼2018-01-21 17:45:17
|
你好，感谢回复！如果这样的话微分方程里面的程序是不是也要响应的改变呢，例如： 将原来的： function dydt = ModelEqs(t,y,beta) % Model equations dydt =((-50000*beta(1)*(125*y-1)^beta(2))*(3-66.67*y)^0.667); 变为： function dydt = ModelEqs(t,y,beta) % Model equations dydt1 =((-50000*beta(1)*(125*y(1)-1)^y(2))*(3-66.67*y(1))^0.667); dydy=[dydt1, y(2)]'; 请问这么改合理吗？ |
3楼2018-01-22 00:02:06
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
nono2009: 金币-50, 屏蔽内容, 违规存档, 违规广告 2018-04-19 10:51:16
nono2009: 金币-50, 屏蔽内容, 违规存档, 违规广告 2018-04-19 10:51:16
|
本帖内容被屏蔽 |
4楼2018-04-06 06:31:32
10