有限差分法解偏微分方程式.PDF

有限差分法解偏微分方程式 第一節 偏微分方程 已知二階偏微分方程 2 2 2 u u u  u u  A B C f x , t,u, , 0 2 2   x xt t  x t  若B 2 4AC 0 ，則稱上式為拋物線型偏微分方程。 若B 2 4AC 0 ，則稱上式為雙曲線型偏微分方程。 若B 2 4AC 0 ，則稱上式為橢圓型偏微分方程。 第二節 Heatconduction偏微分方程之數值解 HeatConductionEquation ： 2 u 1 u ，0 x L ，t 0 2 2 t x  邊界條件與初始條件必須給定 。 L 將 x分成n段 x h ，並取t方向增量k t n 【方法一】顯性近似法 ：       考慮節點 i, j 處之差分公式 ：gridpoints ：x , t ，u u x , t i j ij i j 二階偏微分為中央差分近似 2 u 2u u u 2u u u i1,j i ,j i1,j i1,j i ,j i1,j   x 2 x 2 h2 一階偏微分為前向差分近似 u ui ,j 1 ui ,j ui ,j 1 ui ,j   t t k 2 u 1 u 代入熱傳方程 2 2 ，得 x  t ui1,j 2ui ,j ui1,j 1 ui ,j 1 ui ,j 2 2 k h  1 移項整理得 2 k   u 2u u u u h2 i1,j i ,j i1,j i ,j 1 i ,j 2 k   令* h2 代入，得節點 i, j 之有限差分近似公式   ui ,j 1 ui ,j * ui1,j 2ui,j ui1,j ，i 0,1,2, , n ，j 0,1,2, 