| 查看: 2612 | 回复: 31 | ||
| 本帖产生 1 个 计算强帖 ，点击这里进行查看 | ||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||
[求助]
将PDE转化为状态空间后用matlab求解，结果竟然发散了
|
||
|
系统的对象是一个传热方程，瞬态，并且包含物体的速度。具体就是带钢层流冷却的温度分布。说白了，就是往一块运动着的钢板上喷水，对其温度进行冷却。 我先讲物理方程进行离散，得出来其状态空间形式，既xdot=AX+BU，然后再利用输出方程（既：y=CX的形式）求解其温度分布。 理论上来说，各个取样点的温度应该是逐渐下降的（因为进行喷水冷却了）。但是，我的求解结果却不全是如此。 四个附件的程序运行顺序为：先运行b_1.m求出来状态空间中的稀疏矩阵A和B，再运行l_ode.m 求出状态空间的解，既x的值。然后运行b_1_c.m 求出来输出矩阵C的值，并得出计算结果。exf.m为 l_ode.m 的调用函数。 图片：物理方程为原始温度方程。离散结果为将物理方程离散（用有限元离散）后的状态空间形式。 PS：太穷了，没有金币。  我为此挨了不少骂了，搞的老板现在骂我的力气都没了。唉，全是泪水啊。求坛子里的隐藏的大牛快快给点指点！ 物理方程.png  插值函数.png  离散结果.png |
回复此楼» 本帖附件资源列表
-
欢迎监督和反馈：小木虫仅提供交流平台，不对该内容负责。
本内容由用户自主发布，如果其内容涉及到知识产权问题，其责任在于用户本人，如对版权有异议，请联系邮箱：xiaomuchong@tal.com - 附件 1 : b_1.m
- 附件 2 : b_1_C.m
- 附件 3 : l_ode.m
- 附件 4 : exf.m
2013-06-13 16:25:54, 6.47 K
2013-06-13 16:26:02, 1.23 K
2013-06-13 16:26:13, 418 bytes
2013-06-13 16:26:19, 222 bytes
» 猜你喜欢
法国博士后职位
已经有0人回复
-
重庆交大26年硕士生招生拟调剂通知已出！欢迎加入机器视觉与3D光学成像课题组。
已经有0人回复
-
物理学I论文润色/翻译怎么收费?
已经有127人回复
-
**
已经有1人回复
-
0702一志愿吉大B区求调剂 本科期间发表一篇Sci
已经有3人回复
-
基底STO,薄膜SRO，XRD里面的振荡，是laue震荡还是kiessig振荡? 怎么判断？
已经有2人回复
-
B区学生调剂-兰州交通大学材料科学与工程学院
已经有8人回复
-
山西大同大学物理学专业还有调剂名额，欢迎调剂！
已经有12人回复
-
桂林理工大学物理学专业招收调剂
已经有11人回复
-
VASP 的一组 GPU / CPU 基准测试记录
已经有0人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 【讨论】CFD软件是求解PDE的，matlab中也可以求解PDE
已经有6人回复
5楼2013-06-13 21:34:31
2楼2013-06-13 16:33:58
3楼2013-06-13 20:22:16
4楼2013-06-13 21:22:31
