| 查看: 3094 | 回复: 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
» 猜你喜欢
散金币祈福
已经有187人回复
-
求助，如何提取ELK的rt-TDDFT在某一时刻的自旋密度分布和电子密度分布
已经有1人回复
-
物理学I论文润色/翻译怎么收费?
已经有277人回复
-
澳门大学IAPME李海峰教授常年招收博士研究生
已经有0人回复
-
澳门大学IAPME李海峰教授常年招收博士研究生
已经有0人回复
-
澳门大学IAPME李海峰教授常年招收博士研究生
已经有0人回复
-
深圳大学应用技术学院招聘凝聚态物理博士后
已经有118人回复
-
深圳大学2027级光电信息工程/物理学研究生推免
已经有70人回复
-
求助中医院临床中医中标的青年标书学习，谢谢
已经有4人回复
-
浙大物院表面物理课题组长期招收普博及直博学生
已经有257人回复
-
【征稿】Springer Nature热电专题 Discover Materials / Discover Chemistry
已经有47人回复
15楼2013-06-17 22:10:03
2楼2013-06-13 16:33:58
3楼2013-06-13 20:22:16
4楼2013-06-13 21:22:31
