1. 程式人生 > >【原始碼】用有限差分法求解二維拉普拉斯方程

【原始碼】用有限差分法求解二維拉普拉斯方程

拉普拉斯方程(Laplace’sequation)又稱調和方程、位勢方程,是一種偏微分方程,因由法國數學家拉普拉斯首先提出而得名,涉及領域包括電磁學、天體物理學、力學、數學等。拉普拉斯方程表示液麵曲率與液體表面壓強之間的關係的公式。

使用隱式矩陣求逆技術和顯式迭代求解,使用5點有限差分模板求解二維拉普拉斯方程。

Laplace’s equation is solved in 2d using the 5-point finite differencestencil using both implicit matrix inversion techniques and explicit iterativesolutions.

所使用的邊界條件包括Dirichlet和Neumann條件。

The boundary conditions used include bothDirichlet and Neumann type conditions.

% Solving the 2-D Laplace’sequation by the Finite Difference Method

% 通過有限差分方法求解二維拉普拉斯方程

% Numerical scheme used is an implicitsecond order central difference in space (5-point difference)

% 所使用的數值求解方案是空間中的隱式二階中心差分(5點差分)

% Specifying parameters

% 引數設定

nx=100; %Number of steps inspace(x) 空間x的步數

ny=100; %Number of steps inspace(y)

dx=2/(nx-1); %Width of space step(x) 空間x的寬度

dy=2/(ny-1); %Width of space step(y)

……

下載相關原始碼地址:

更多精彩文章請關注微訊號:在這裡插入圖片描述