 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 查看: 2541 | 回復(fù): 31 | ||
| 本帖產(chǎn)生 1 個 計算強帖 ，點擊這里進(jìn)行查看 | ||
[求助]
將PDE轉(zhuǎn)化為狀態(tài)空間后用matlab求解，結(jié)果竟然發(fā)散了
|
||
|
系統(tǒng)的對象是一個傳熱方程，瞬態(tài)，并且包含物體的速度。具體就是帶鋼層流冷卻的溫度分布。說白了，就是往一塊運動著的鋼板上噴水，對其溫度進(jìn)行冷卻。 我先講物理方程進(jìn)行離散，得出來其狀態(tài)空間形式，既xdot=AX+BU，然后再利用輸出方程（既：y=CX的形式）求解其溫度分布。 理論上來說，各個取樣點的溫度應(yīng)該是逐漸下降的（因為進(jìn)行噴水冷卻了）。但是，我的求解結(jié)果卻不全是如此。 四個附件的程序運行順序為：先運行b_1.m求出來狀態(tài)空間中的稀疏矩陣A和B，再運行l(wèi)_ode.m 求出狀態(tài)空間的解，既x的值。然后運行b_1_c.m 求出來輸出矩陣C的值，并得出計算結(jié)果。exf.m為 l_ode.m 的調(diào)用函數(shù)。 圖片：物理方程為原始溫度方程。離散結(jié)果為將物理方程離散（用有限元離散）后的狀態(tài)空間形式。 PS：太窮了，沒有金幣。  我為此挨了不少罵了，搞的老板現(xiàn)在罵我的力氣都沒了。唉，全是淚水啊。求壇子里的隱藏的大�？炜旖o點指點！物理方程.png 插值函數(shù).png 離散結(jié)果.png |
回復(fù)此樓» 本帖附件資源列表
-
歡迎監(jiān)督和反饋：小木蟲僅提供交流平臺，不對該內(nèi)容負(fù)責(zé)。
本內(nèi)容由用戶自主發(fā)布，如果其內(nèi)容涉及到知識產(chǎn)權(quán)問題，其責(zé)任在于用戶本人，如對版權(quán)有異議，請聯(lián)系郵箱：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
» 猜你喜歡
2026年循環(huán)經(jīng)濟功能材料國際會議（ICFMCE 2026）
已經(jīng)有0人回復(fù)
-
2026年第五屆電氣、電子與信息工程國際會議（ISEEIE 2026）
已經(jīng)有0人回復(fù)
-
物理學(xué)I論文潤色/翻譯怎么收費?
已經(jīng)有281人回復(fù)
-
哈爾濱理工大學(xué)物理系招收物理學(xué)考研調(diào)劑
已經(jīng)有5人回復(fù)
-
0702一志愿吉大B區(qū)求調(diào)劑
已經(jīng)有5人回復(fù)
-
求調(diào)劑
已經(jīng)有0人回復(fù)
-
0702一志愿吉大B區(qū)求調(diào)劑有論文
已經(jīng)有0人回復(fù)
-
請問還有沒有用Latex寫文章的小伙伴們？
已經(jīng)有0人回復(fù)
-
光學(xué)工程學(xué)碩調(diào)劑信息
已經(jīng)有26人回復(fù)
-
歡迎加入課題組
已經(jīng)有0人回復(fù)
-
散金幣，求好運，祝面上順利！
已經(jīng)有34人回復(fù)
» 本主題相關(guān)價值貼推薦，對您同樣有幫助:
- 【討論】CFD軟件是求解PDE的，matlab中也可以求解PDE
已經(jīng)有6人回復(fù)
，下班后參考你的模型隨便弄的二維，一維的我覺得你壓力山大，二維的考驗?zāi)愕木幊�，雖然在comsol里很容易搞，你可以去試試。但是我覺得了解方程本身比同軟件更值得花時間
21樓2013-06-18 18:25:59
2樓2013-06-13 16:33:58
3樓2013-06-13 20:22:16
4樓2013-06-13 21:22:31
5樓2013-06-13 21:34:31
6樓2013-06-14 09:13:22
7樓2013-06-14 09:38:15
8樓2013-06-14 12:27:48
9樓2013-06-14 21:27:22
10樓2013-06-15 08:05:43
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 yuahua 的主題更新
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 311求調(diào)劑 +9 | 冬十三 2026-03-15 | 10/500 |
|
|---|---|---|---|---|
|
|
[考研] 085601專碩，總分342求調(diào)劑，地區(qū)不限 +4 | share_joy 2026-03-16 | 4/200 |
|
|
[考研] 268求調(diào)劑 +7 | 好運連綿不絕 2026-03-12 | 8/400 |
|
|
|
[考研] 308求調(diào)劑 +4 | 是Lupa啊 2026-03-16 | 4/200 |
|
|
|
[考研] 274求調(diào)劑0856材料化工 +13 | z2839474511 2026-03-11 | 14/700 |
|
|
|
[考研] 211本，11408一志愿中科院277分，曾在中科院自動化所實習(xí) +6 | Losir 2026-03-12 | 7/350 |
|
|
|
[考研] 278求調(diào)劑 +3 | Yy7400 2026-03-13 | 3/150 |
|
|
|
[考研] 0856求調(diào)劑 +3 | 劉夢微 2026-03-15 | 3/150 |
|
|
|
[考研] 085601材料工程315分求調(diào)劑 +3 | yang_0104 2026-03-15 | 3/150 |
|
|
|
[考研] 288求調(diào)劑 +4 | 奇點0314 2026-03-14 | 4/200 |
|
|
|
[考研] 289求調(diào)劑 +4 | 這么名字咋樣 2026-03-14 | 6/300 |
|
|
|
[基金申請] 現(xiàn)在如何回避去年的某一個專家，不知道名字 +3 | zk200107 2026-03-12 | 6/300 |
|
|
|
[考研] 308 085701 四六級已過求調(diào)劑 +7 | 溫喬喬喬喬 2026-03-12 | 14/700 |
|
|
|
[考研] 0703，333分求調(diào)劑 一志愿鄭州大學(xué)-物理化學(xué) +3 | 李魔女斗篷 2026-03-11 | 3/150 |
|
|
|
[考研] 26調(diào)劑/材料/英一數(shù)二/總分289/已過A區(qū)線 +6 | 步川酷紫123 2026-03-13 | 6/300 |
|
|
|
[考研] 281求調(diào)劑 +9 | Koxui 2026-03-12 | 11/550 |
|
|
|
[考研] 【考研調(diào)劑求收留】 +3 | Ceciilia 2026-03-11 | 3/150 |
|
|
|
[考研] 307求調(diào)劑 +5 | 超級伊昂大王 2026-03-12 | 5/250 |
|
|
|
[考研] 求調(diào)劑 +3 | 程雨杭 2026-03-12 | 3/150 |
|
|
|
[考研] 289求調(diào)劑 +3 | 李政瑩 2026-03-12 | 3/150 |
|