介紹了對複雜網格的變換
標準的有限差分法要求在均勻網格上進行，但是在實際當中，比如速度邊界層（越靠近邊界，速度變化越大），還用均勻的網格顯然會造成較大的誤差（如果網格劃分的很細的話也可以減小誤差，但計算量太大）因此對於這種情況應該構造非均勻的網格，針對非均勻網格構造差分主要採用座標變換法
如上圖的一維非均勻網格，求解的思路是：嘗試將物理平面上的非均勻網格轉移到計算平面的均勻網格中，如下圖所示：
在轉換過程中我們可以得到關係式
根據該關係式可以得到
稱作jacobi矩陣這種變換要求相鄰網格差距不能過大，否則會帶來較大誤差，通常要求在10%以內。
二維情況仍是採用座標變換的方式，如下圖
根據變換情況可知：
從而 由於是可以利用差分法計算出的
那麼進行下列變換得到
其中 ，物理含義是網格單元面積
此外也可以得到一些性質：（證明略）
對原方程左右兩端同乘
，利用上述性質可以得到
所以原方程可以寫成如下形式：
其中類似於二維情況，可得
其中
有限差分法MATLAB程式
設有乙個長直接地金屬矩形槽，長a 40，寬b 20，其側壁與底面電位均為零，頂蓋電位為100v 相對值 求槽內電位分布。利用高斯迭代求解 如下 相鄰兩次迭代值最大允許誤差為0.001 a zeros 21,41 a 1,100 b zeros 19,39 c eye 19,39 count 1 d ...
有限差分法Eluer演算法（求解常微分方程）
有限差分法c 實現 幾種eluer演算法 在matlab 裡面寫乙個程式 要求用 隱式尤拉法 backward euler 去解決常微分方程。下面是兩個例題，1.x 2x 還給出準確值是 x e 2t 要求 求出t 0，這是乙個初始值，然後算出在區間 0,5 的值。還給出 步長 step size ...
有限元法 有限差分法 有限體積法
有限元法也叫有限單元法 finite element method,fem 是隨著電子計算機的發展而迅速發展起來的一種彈性力學問題的數值求解方法。五十年代初，它首先應用於連續體力學領域 飛機結構靜 動態特性分析中，用以求得結構的變形 應力 固有頻率以及振型。由於這種方法的有效性，有限單元法的應用已從...