| 7 | 1/1 | 返回列表 |
| 查看: 1114 | 回復(fù): 6 | ||
[求助]
用Matlab編程序來解決一個(gè)小問題,但結(jié)果不太滿意,哪位蟲友大神能給些建議哈 已有1人參與
|
原題是傳熱學(xué)穩(wěn)態(tài)傳熱的特別簡單的一個(gè)題,雖然之前試過用mathematica的Dsolve函數(shù),一步就可以算出來符號(hào)解析解,但最近剛開始自學(xué)Matlab,所以想用Matlab來做一下數(shù)值計(jì)算試試,,但果然還是剛開始學(xué),也是第一次用Matlab來編程,編的程序問題很多,,自己感覺最大的兩個(gè)問題1.對(duì)矩陣的賦值,全用了循環(huán)語句,明明Matlab的循環(huán)效率很差的,,但不知道怎樣寫比較簡便。2,循環(huán)達(dá)到穩(wěn)定的判斷條件太麻煩,不夠準(zhǔn)確而且?隙ㄟ有其它各種問題,得到的結(jié)果在邊界的有些地方好像還不太對(duì),,哪位大神來幫忙改一改,謝謝謝謝謝了 |
至尊木蟲 (文壇精英)
你編寫的邊界條件確實(shí)不簡潔。一般來講這類問題可以用有限差分法很好的解決,而有限差分一般使用迭代方法求值。你的邊界條件完全可以使用向量。如你的勢函數(shù)是dot=zeros(100,100)當(dāng)然只是舉個(gè)例子,你的四個(gè)邊界就是dot(1, ,dot(:,100),dot(:,1),dot(:,100)。你完全可以dot(1, =tw1,類似的定義其他幾個(gè)邊界。迭代的話設(shè)置迭代次數(shù),然后比較上一遍迭代值和這一遍迭代值之差,當(dāng)達(dá)到你要的精度后,迭代終止。 |

|
enen, ,邊界賦初值那個(gè)我明白了,,但迭代那個(gè)還是不太清楚,,之前自己試過讓循環(huán)的判斷條件是上一遍與這一遍各點(diǎn)之間的差值小于某個(gè)精度,比如0.01,但最后結(jié)果表明,顯然沒有達(dá)到穩(wěn)態(tài),,節(jié)點(diǎn)稍多一些的話,區(qū)域中間那塊兒都是凹下去的。。。。是不是我那個(gè)循環(huán)編的不對(duì),,自己看不少書但里面關(guān)于有限差分的講的都比較少,也沒講怎么在Matlab上實(shí)現(xiàn),所以我這點(diǎn)還是不清楚,,能麻煩再大致解釋一下么? 謝謝了 |
至尊木蟲 (文壇精英)

| 7 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 一志愿天大材料與化工(085600)總分338 +3 | 蔡大美女 2026-03-13 | 3/150 |
|
|---|---|---|---|---|
|
[考研] 298-一志愿中國農(nóng)業(yè)大學(xué)-求調(diào)劑 +7 | 手機(jī)用戶 2026-03-17 | 7/350 |
|
|
[考研] 0703化學(xué)調(diào)劑 +3 | 妮妮ninicgb 2026-03-17 | 3/150 |
|
|
[考研] 材料工程專碩調(diào)劑 +5 | 204818@lcx 2026-03-17 | 5/250 |
|
|
[考研] 302求調(diào)劑 +9 | 負(fù)心者當(dāng)誅 2026-03-11 | 9/450 |
|
|
[考研] 材料與化工專碩調(diào)劑 +5 | heming3743 2026-03-16 | 5/250 |
|
|
[考博] 26申博 +4 | 八旬速覽 2026-03-16 | 4/200 |
|
|
[基金申請(qǐng)] 國自科面上基金字體 +6 | iwuli 2026-03-12 | 7/350 |
|
|
[考研] 中科院材料273求調(diào)劑 +4 | yzydy 2026-03-15 | 4/200 |
|
|
[考研] 0703化學(xué)調(diào)劑,求各位老師收留 +8 | 秋有木北 2026-03-14 | 8/400 |
|
|
[考研] 0703 物理化學(xué)調(diào)劑 +3 | 我可以上岸的對(duì)?/a> 2026-03-13 | 5/250 |
|
|
[考研] 326求調(diào)劑 +3 | mlpqaz03 2026-03-15 | 3/150 |
|
|
[考研] 0856專碩279求調(diào)劑 +5 | 加油加油!? 2026-03-15 | 5/250 |
|
|
[考研] 265求調(diào)劑 +4 | 威化餅07 2026-03-12 | 4/200 |
|
|
[考研] 材料080500調(diào)劑求收留 +3 | 一顆meteor 2026-03-13 | 3/150 |
|
|
[考研] 330求調(diào)劑 +3 | ?醬給調(diào)劑跪了 2026-03-13 | 3/150 |
|
|
[考研] 0856材料與化工301求調(diào)劑 +5 | 奕束光 2026-03-13 | 5/250 |
|
|
[考研] 26調(diào)劑/材料/英一數(shù)二/總分289/已過A區(qū)線 +6 | 步川酷紫123 2026-03-13 | 6/300 |
|
|
[考研] 0703化學(xué)求調(diào)劑 +7 | 綠豆芹菜湯 2026-03-12 | 7/350 |
|
|
[考研] 308求調(diào)劑 +3 | 是Lupa啊 2026-03-12 | 3/150 |
|