| 查看: 1419 | 回复: 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程序呢？我遇到了这个问题，目前束手无策，想请哪位大神回复一下。多谢！ |
回复此楼» 猜你喜欢
2026年的国家社科基金项目通讯评审的新规则与新动向、新挑战
已经有7人回复
-
今天基金会出结果吗？20260819
已经有18人回复
-
我面上完蛋了
已经有14人回复
-
范进中举一文的中心思想
已经有8人回复
-
只有每年这种时候来逛逛小木虫
已经有27人回复
-
明天应该可查了！？
已经有5人回复
-
如果此刻你正在为国基感到焦虑，不妨来听听这首《基金之外》
已经有7人回复
-
没有任何消息-是不是就凉了
已经有8人回复
-
人气不行了
已经有10人回复
-
能否退出参与的面上项目解除限项
已经有24人回复
chendequan
铁虫 (小有名气)
- 应助: 152 (高中生)
- 金币: 918.6
- 红花: 17
- 帖子: 262
- 在线: 208.7小时
- 虫号: 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
20