[求助] matlab解微分方程

» 猜你喜欢

Springer期刊投稿求助已经有4人回复
读博已经有3人回复
博士读完未来一定会好吗已经有10人回复
小论文投稿已经有3人回复
Bioresource Technology期刊，第一次返修的时候被退回好几次了已经有9人回复
心脉受损已经有3人回复
到新单位后，换了新的研究方向，没有团队，持续积累2区以上论文，能申请到面上吗已经有8人回复
申请2026年博士已经有6人回复
请问哪里可以有青B申请的本子可以借鉴一下。已经有5人回复

» 本主题相关价值贴推荐，对您同样有帮助:

用matlab求解一偏微分方程，求助！已经有1人回复
流体力学一维定常不可压缩流体ns方程，已知压差变化，怎么求解速度u关于x的常微分方程已经有0人回复
matlab书籍已经有103人回复
求大神帮忙解题啊！！！已经有4人回复
Matlab求解微分方程里面有参数带入怎么求解，求救！！！已经有0人回复
求助CFD软件可以自定义模拟算法已经有0人回复
【紧急求助】带有微分平方项的ODE求解已经有11人回复
在控制器设计中关于大家如何给定？如果给定是方波怎么给？已经有2人回复
请教高手：PDE工具箱能解两个偏微分方程相互迭代求解吗已经有2人回复
【求助】求解形如y''+a*y'^2+by'+c=0 二阶微分方程已经有14人回复
【资料】一些数学实验,有程序源代码已经有7人回复

1楼2012-03-06 15:25:15

dbb627

荣誉版主 (著名写手)

专家经验: +4
仿真EPI: 6
应助: 289 (大学生)
贵宾: 0.589
金币: 24640.4
散金: 551
红花: 61
沙发: 1
帖子: 1246
在线: 1794.8小时
虫号: 149791
注册: 2005-12-29
性别: GG
专业: 污染控制化学
管辖: 计算模拟

【答案】应助回帖

★ ★ ★ ★
感谢参与，应助指数 +1
御剑江湖(金币+3): 呵呵，谢谢 2012-03-07 20:05:33
xiegangmai(金币+1): 谢谢参与！ 2012-03-08 18:40:39

当然有，解什么样的微分方程？
常微分？偏微分？
常微分一帮用ode解算器有ode15 ode23 ode45还有解刚性和延迟微分方程等等
偏微分pde，就比较复杂些，只能解某些特定类型的椭圆，抛物线，双曲型等

The more you learn, the more you know, the more you know, and the more you forget. The more you forget, the less you know. So why bother to learn.

2楼2012-03-06 18:47:20

hanvs2009

银虫 (小有名气)

应助: 1 (幼儿园)
金币: 1079.8
帖子: 263
在线: 9小时
虫号: 787494
注册: 2009-06-05
性别: GG
专业: 计算机软件

【答案】应助回帖

★ ★ ★
感谢参与，应助指数 +1
御剑江湖(金币+2): 呵呵，谢谢 2012-03-07 20:05:41
xiegangmai(金币+1): 谢谢参与！ 2012-03-08 18:40:46

去百度文库下载一个文章，就叫MATLABA解微分方程，照着做应该就会了。

3楼2012-03-07 07:30:14

zhufayan

金虫 (正式写手)

西北狼

应助: 0 (幼儿园)
金币: 406.1
散金: 555
红花: 4
帖子: 587
在线: 134.1小时
虫号: 798203
注册: 2009-06-24
性别: GG
专业: 配位化学

-dc/dt=k(c0-c)^1/3(c-c∞)
-dc/dt=k(c0-c)^2/3(c-c∞)p
-dc/dt=k(c0-c)^4/3(c-c∞)p
-dc/dt=k (c-c∞)n
现在我知道上述每一微分方程的c,t和c0（其中c为浓度，t为时间，c0为初始浓度）要用拟合的方法给出c∞和p，书上说用runge-kutta计算方法，但是我看不懂，请指教，不胜感激！
t c
0.25 0.0873
0.5 0.1172
0.7 0.136
1 0.1486
1.5 0.1581
2.5 0.1597
4 0.1601
7 0.166

4楼2012-03-07 10:09:00

dbb627

荣誉版主 (著名写手)

专家经验: +4
仿真EPI: 6
应助: 289 (大学生)
贵宾: 0.589
金币: 24640.4
散金: 551
红花: 61
沙发: 1
帖子: 1246
在线: 1794.8小时
虫号: 149791
注册: 2005-12-29
性别: GG
专业: 污染控制化学
管辖: 计算模拟

★ ★ ★
御剑江湖(金币+3): 呵呵，谢谢 2012-03-07 20:05:54

引用回帖:

4楼: Originally posted by zhufayan at 2012-03-07 10:09:00:
-dc/dt=k(c0-c)^1/3(c-c∞)
-dc/dt=k(c0-c)^2/3(c-c∞)p
-dc/dt=k(c0-c)^4/3(c-c∞)p
-dc/dt=k (c-c∞)n
现在我知道上述每一微分方程的c,t和c0（其中c为浓度，t为时间，c0为初始浓度）要用拟合的方法给出c∞和 ...

你的公式-dc/dt=k(c0-c)^1/3(c-c∞)这个公式的意思是什么
等价于-dc/dt=k*(c-c∞)*(c0-c)^1/3？还是-dc/dt=k*(c-c∞)*(c0-c)^1/[3*(c-c∞)]?

你的数据是给那个公式用，按照你的模型浓度是降低的，怎么c的数据反而增大了呢

The more you learn, the more you know, the more you know, and the more you forget. The more you forget, the less you know. So why bother to learn.

5楼2012-03-07 12:06:42

zhufayan

金虫 (正式写手)

西北狼

应助: 0 (幼儿园)
金币: 406.1
散金: 555
红花: 4
帖子: 587
在线: 134.1小时
虫号: 798203
注册: 2009-06-24
性别: GG
专业: 配位化学

引用回帖:

楼: Originally posted by dbb627 at 2012-03-07 12:06:42:
你的公式-dc/dt=k(c0-c)^1/3(c-c∞)这个公式的意思是什么
等价于-dc/dt=k*(c-c∞)*(c0-c)^1/3？还是-dc/dt=k*(c-c∞)*(c0-c)^1/?

你的数据是给那个公式用，按照你的模型浓度是降低的，怎么c的数据反而增大了呢

是这个公式-dc/dt=k*(c0-c)^0.66*(c-c∞)^p
，它是溶解动力学模型，上述把数据给错了，不好意思，新数据如下：现在就是要通过拟合给出c∞和p
t c
11 0.166
24 0.1678
31.5 0.161
38.5 0.1554
39 0.146
40 0.1386
42 0.1256
46 0.1184
51.5 0.1178
72 0.1181

6楼2012-03-07 15:36:59

zhufayan

金虫 (正式写手)

西北狼

应助: 0 (幼儿园)
金币: 406.1
散金: 555
红花: 4
帖子: 587
在线: 134.1小时
虫号: 798203
注册: 2009-06-24
性别: GG
专业: 配位化学

是这个公式：-dc/dt=k*(c0-c)^0.66*(c-c∞)^p
上述把数据给错了，不好意思，新数据如下：
11 0.166
24 0.1678
31.5 0.161
38.5 0.1554
39 0.146
40 0.1386
42 0.1256
46 0.1184
51.5 0.1178
72 0.1181
我的目的就是通过上述动力学公式通，过拟合给出c∞和p两个参数，十分感谢你的帮助

7楼2012-03-07 15:40:24

dbb627

荣誉版主 (著名写手)

专家经验: +4
仿真EPI: 6
应助: 289 (大学生)
贵宾: 0.589
金币: 24640.4
散金: 551
红花: 61
沙发: 1
帖子: 1246
在线: 1794.8小时
虫号: 149791
注册: 2005-12-29
性别: GG
专业: 污染控制化学
管辖: 计算模拟

引用回帖:

7楼: Originally posted by zhufayan at 2012-03-07 15:40:24:
是这个公式：-dc/dt=k*(c0-c)^0.66*(c-c∞)^p
上述把数据给错了，不好意思，新数据如下：
11 0.166
24 0.1678
31.5 0.161
38.5 0.1554
39 0.146
40 0.1386
42 0.1256
46 0.1184
51.5 0.1178
72 0.1181 ...

t为时间，c0为初始浓度，c0呢？

The more you learn, the more you know, the more you know, and the more you forget. The more you forget, the less you know. So why bother to learn.

8楼2012-03-08 09:53:55