本文關(guān)鍵詞：計(jì)算物理，由筆耕文化傳播整理發(fā)布。

第四章《有限差分方法》習(xí)題參考答案
? ? 2? ( x, y ) ? 0, 1. 解：依題，對(duì)于拉普拉斯方程 ? ?? ( x, 0) ? ? ( x,1) ? 0, ? (0, y ) ? ? (1, y ) ? 1.

有限差分計(jì)算步驟如下： i).首先用正方形網(wǎng)格劃分區(qū)域 D： 0 ? x ? 1, 0 ? y ? 1. 其中， xi ? h

i, y j ? hj, i, j ? 1,? , N , h ? ii).寫(xiě)出有限差分方程： 由微分的中心差值表示，
? ? 2? ? ?i ?1, j ? 2?i , j ? ?i ?1, j ? ? 2? ? ?i , j ?1 ? 2?i , j ? ?i , j ?1 , ? 2? ? . ? 2? ? h2 h2 ? ?x ?ij ? ?y ?ij

1 . N

有所求的拉普拉斯方程對(duì)應(yīng)的差分方程表達(dá)式
1 ? ??i , j ? (?i ?1, j ? ?i ?1, j ? ?i , j ?1 ? ?i , j ?1 ) ? 0, (i, j ) ? D; 4 ? ? ?i , j ? gij , (i, j ) ? ?D. ?

引入層向量，
? ?1, j ? ? ?1 ? ? ? ? ? ? ? ? ? ? , ? j ? ? ? ? , j ? 1,? , N ? 1; ?? ? ?? ? ? N ?1 ? ? N ?1, j ? ? ? g 0, j ? ? ? 0 ? ? b1 ? ? 1 ? ? B ? ? ? ? , b j ? ? ? ? ? , j ? 2,? , N ; 4? ? ?b ? ? N ?1 ? ? 0 ? ? ?g ? ? N, j ?

? ? G ? ?? I ? 4 K ?? ? ? ? ? ? 0 ?

?

I 4

? ? I 4

G ? ?

? ? I 4

? ? 0 ? ? 1 ? ? ? ?? 1 ? ? ? 4 ?, G ? ? I ? ? ? ? ? 4? ? ? ? 0 G ? ? ? K? ? B .

?

1 4

? ? ? ? 1 4

1 ? ?

? 0 ? ? ? ? ? ? 1 ? ? 4? ? 1 ? ?

有差分方程：

iii).用迭代法解差分方程：
由 ? ? R? ? B, R ? I ? K .

有 ? ( k ) ? R? ( k ?1) ? B, ?i ( k ) ? ? Rij? j
j

( k ?1)

? Bi .

引入松弛因子
??
2 ? [1, 2], 其中(l ? 1)為每邊節(jié)點(diǎn)數(shù). 1 ? sin(? / l )

? ( k ?1) ? ?? GS ( k ?1) ? (1 ? ? )? ( k ) ? R? ? ( k ) ? ? ( I ? ? L) ?1 B, R? ? ? ( I ? ? L) ?1[?U ? (1 ? ? ) L].
?0 ? 0? ?0 ? ? ? ? 其中， L ? U ? R, L ? ? ? ? ? , U ? ? ? ? ? ?. ? ? ? ? 0? ? ?0 ? 0?

先賦初始值： ? (0) ? {?i(0) }, i ? 1,? , N . 再運(yùn)用超松弛迭代法
( k ?1) (k ) ?ij( k ?1) ? ?? GS ij ? (1 ? ? )?ij

=?ij( k ) ?
當(dāng)

?
4

) (k ) (k ) (k ) (k ) (?i(?k1, j ? ?i , j ?1 ? ?i ?1, j ? ?i , j ?1 ? 4?ij )

?(k ) ?
(k )

? ? 時(shí)，stop ! 其中，，移位矢量 ? ( k ) ? {? i( k ) } ? {?i( k ) ? ?i( k ?1) }.

4. 證：依題，中子擴(kuò)散方程的形式為 ? 2? ? f ( x, y )? ? q( x, y ).

其中， f ( x, y) ? ?a 2 , q( x, y ) ? ? sin ? ?

?? y ? ? sin ? ?. ? L ? ? L ?

?x?

則依“五點(diǎn)差分”格式（正方形網(wǎng)格劃分） ，
?ij ?
? 2 ? ? ? ih ? ? ? jh ? ? h sin ? L ? sin ? L ? ? (?i ?1, j ? ?i ?1, j ? ?i , j ?1 ? ?i , j ?1 ) ? ? ? ? ? ? ? 2 ? A? ? h S (i, j ) ? (?i ?1, j ? ?i ?1, j ? ?i , j ?1 ? ?i , j ?1 ) ? ? 1 4 ? a 2h2 1 ? ? ih ? ? ? jh ? , S (i, j ) ? sin ? ? sin ? ?. 2 2 4?a h ? L ? ? L ?

（1）

其中， A ?

由已知差分方程
? 8 ? ? h ?? ?ij ? S (i. j ) B, B ? ? a 2 ? 2 sin 2 ? ? ? . h ? 2L ?? ?
?1

有
?i ?1, j ? B sin ?
? (i ? 1)? h ? ? ? jh ? ? sin ? ? L ? ? ? L ? ? ? i? h ? ??h ? ? i? h ? ? ? h ? ? ? ? jh ? ? B ?sin ? ? cos ? ? ? cos ? ? sin ? ? ? sin ? ?, ? L ? ? L ? ? L ?? ? L ? ? ? L ? ? (i ? 1)? h ? ? ? jh ? ? sin ? ? L ? ? ? L ? ? ? i? h ? ??h ? ? i? h ? ? ? h ? ? ? ? jh ? ? B ?sin ? ? cos ? ? ? cos ? ? sin ? ? ? sin ? ?, ? L ? ? L ? ? L ?? ? L ? ? ? L ?

?i ?1, j ? B sin ?

?i , j ?1 ? B sin ?

? i? h ? ? ( j ? 1)? h ? ? sin ? ? L ? L ? ? ? ? i? h ? ? ? j? h ? ??h ? ? j? h ? ? ? h ? ? ? B sin ? ? ?sin ? ? cos ? ? ? cos ? ? sin ? ?? , ? L ?? ? L ? ? L ? ? L ? ? L ?? ? i? h ? ? ( j ? 1)? h ? ?i , j ?1 ? B sin ? ? sin ? ? L ? L ? ? ? ? i? h ? ? ? j? h ? ??h ? ? j? h ? ? ? h ? ? ? B sin ? ? ?sin ? ? cos ? ? ? cos ? ? sin ? ?? . ? L ?? ? L ? ? L ? ? L ? ? L ??

把上面四個(gè)式子代入方程（1）中，有
? ? ? ? h ?? ??h ? ? ?ij ? A ? h 2 S (i, j ) ? BS (i. j )4 ? ? ? =S (i, j ) ? Ah 2 ? AB[4 ? ?]? . ? L ?? ? L ?? ? ?

下面證明： B ? Ah 2 ? AB[4 cos ? ?

?] ， ? L ?

?h ?

（2）

h2 Ah 2 4 ? a 2h2 ? ??h ? ??h ? 2 2 1 ? 4 A cos ? ? 1 ? 4 cos ? ? / (4 ? a h ) ? L ? ? L ? 事實(shí)上，由于 1 1 ?B ? ? 4? ? ? h ? ? a 2 ? 4 sin 2 ? ? h ? 2 a ? 2 ? 1 ? cos ? ? ? ?? h2 h ? ? L ? ? L ??

顯見(jiàn)， （2）式得證。也即
? ? ih ? ? ? jh ? ? 2 8 2 ? ? h ?? ?ij ? S (i, j ) B ? sin ? ? sin ? ? ? a ? 2 sin ? ?? . h ? L ? ? L ?? ? 2L ??
?1



更多相關(guān)文檔:

原子物理學(xué)第四章習(xí)題解答

原子物理學(xué)第四章習(xí)題解答_理學(xué)_高等教育_教育專(zhuān)區(qū)...4-11 試計(jì)算在B為2.5T的磁場(chǎng)中,鈉原子的D雙線(xiàn)所...2014小學(xué)教師資格考試《... 2014年幼兒園教師資格考...

計(jì)算物理(彭芳麟)課后部分答案

《計(jì)算物理》第四章習(xí)題... 4頁(yè) 免費(fèi) 《計(jì)算物理》第三章習(xí)題... 4頁(yè) 免費(fèi) 《計(jì)算物理》第五章習(xí)題... 4頁(yè) 免費(fèi) 計(jì)算物理參考答案(修正版... 暫無(wú)評(píng)價(jià)...

計(jì)算物理考試題及答案

計(jì)算物理考試題及答案_理學(xué)_高等教育_教育專(zhuān)區(qū)。:page...是二十世紀(jì)四十年代中期由于 科學(xué)技術(shù)的發(fā)展和電子...

計(jì)算物理課后部分習(xí)題答案2

《計(jì)算物理》第四章習(xí)題參... 4頁(yè) 免費(fèi)喜歡此文檔的還喜歡 計(jì)算物理學(xué)(科學(xué)...計(jì)算物理課后部分習(xí)題答案2 隱藏>> 1.1 Show that the error in the nth-orde...

大學(xué)物理學(xué)-計(jì)算物理-編程練習(xí)題答案

大學(xué)物理學(xué)-計(jì)算物理-編程練習(xí)題答案_理學(xué)_高等教育_教育專(zhuān)區(qū)。clear all close...2014小學(xué)教師資格考試《... 2014年幼兒園教師資格考... 2014教師資格中學(xué)教育知...

物理計(jì)算練習(xí)及答案

物理計(jì)算練習(xí)及答案_理化生_初中教育_教育專(zhuān)區(qū)。14 周聯(lián)考復(fù)習(xí)資料——計(jì)算專(zhuān)題...2014小學(xué)教師資格考試《... 2014年幼兒園教師資格考... 2014教師資格中學(xué)教育知...

《操作系統(tǒng)》第4章教材習(xí)題解答

《操作系統(tǒng)》第4章教材習(xí)題解答_其它考試_資格考試/...10 4 7 計(jì)算邏輯地址0A5C(H)所對(duì)應(yīng)的物理地址。...操作系統(tǒng)課后習(xí)題的答案... 77頁(yè) 免費(fèi) 第...

物理計(jì)算題習(xí)題解答

大​學(xué)​物​理​期​末​計(jì)​算​題​習(xí)​題​解​答習(xí)題解答第三章 3-3 解: (1)設(shè)兩物體的加速度為 a,滑輪的受力如圖,設(shè)滑輪的...

操作系統(tǒng)習(xí)題及答案四

TLB 是是否否否是是 (1)計(jì)算兩個(gè)邏輯地址:0AC5H、1AC5H 對(duì)應(yīng)的物理地址...操作系統(tǒng)習(xí)題第四章答案 5頁(yè) 免費(fèi) 操作系統(tǒng)練習(xí)題四參考答... 4頁(yè) 免費(fèi) ©...

物理計(jì)算題及答案

物理計(jì)算題及答案_調(diào)查/報(bào)告_表格/模板_實(shí)用文檔。...所以并聯(lián)部分分得的電壓是第三只串聯(lián) 聯(lián)電路有中...2014小學(xué)教師資格考試《... 2014年幼兒園教師資格考...

更多相關(guān)標(biāo)簽:

  本文關(guān)鍵詞：計(jì)算物理，由筆耕文化傳播整理發(fā)布。