| 查看: 1428 | 回复: 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程序呢？我遇到了这个问题，目前束手无策，想请哪位大神回复一下。多谢！ |
回复此楼» 猜你喜欢
我就是申请一个面上项目而已，这评审意见是按照杰青的条件评的吧？
已经有3人回复
-
有多少人是今天查系统知道结果的？
已经有18人回复
-
国自然评审意见
已经有15人回复
-
为什么资助数各大高校都创新高，自己申请怎么就这么难
已经有13人回复
-
国自然面上复盘~欢迎讨论
已经有14人回复
-
系统查不到
已经有11人回复
-
有没有仍没收到信息的
已经有3人回复
-
导师吐槽：我怎么摊上了这么个极品研究生！
已经有8人回复
-
基金不中，共勉
已经有11人回复
-
麻烦专家们看看评委们的意见（F口面上）
已经有9人回复
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
