| 查看: 3238 | 回复: 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
» 猜你喜欢
27申博推荐：南京航空航天大学国际前沿院光学方向招收博士生！
已经有19人回复
-
华南师范大学+双一流学科物理学+广东省杰青团队2027年博士招生+机器学习驱动材料设计
已经有1人回复
-
物理学I论文润色/翻译怎么收费?
已经有148人回复
-
apl投稿状态
已经有2人回复
-
华南师范大学+双一流学科物理学+广东省杰青团队2027年博士招生+机器学习驱动材料设计
已经有1人回复
-
日内瓦国际发明展申报完整流程（科研人员通用指南）
已经有11人回复
-
3月 中国石油和化学工业联合会科学技术奖 马上开始，申获奖要点
已经有15人回复
-
中国石油和化学工业联合会科学技术奖的评审标准是什么？
已经有5人回复
-
经费如何使用，花不完怎么办？
已经有5人回复
-
兼职
已经有12人回复
-
海外Gustau Catalan课题组情况了解
已经有1人回复
17楼2013-06-17 22:25:23
2楼2013-06-13 16:33:58
3楼2013-06-13 20:22:16
4楼2013-06-13 21:22:31
