【原始碼】用有限差分法求解二維拉普拉斯方程
拉普拉斯方程(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)
……
下載相關原始碼地址:
更多精彩文章請關注微訊號: