matlab有限差分法编程波导_有限差分法及matlab实现.pdf

有限差分法及matlab实现

有限差分法解静电场的边值问题的算法实现及相关问题讨论：

王宁远

中国科学技术大学 09 级物理 2 班

E-mail wny@mail.ustc.edu.cn

摘要：

本文用 MATLAB 实现了有限差分法解静电场边值问题的算法，将偏微分方程的问题化为线性方程

组问题，并使用了迭代法进行线性方程组的数值解。讨论了从几个角度去优化迭代法的措施。并运用这

样的方法解决了文 [1]的闪电模拟问题，，使用了更优化的算法对重新进行了计算，并一定程度上改进了

模型,讨论了几个与文[1]所持的不同的观点。

正文：

经典场的边值问题在数学上表达为泊松方程和拉普拉斯方程，但解偏微分方程往往是困难的。幸而

很多时候对于具体问题我们需要的不是解析解，而是数值解，所以可以考虑用连续变量离散化的方法求

出数值解,在足够的精度上进行逼进,这就引出了有限差分法。

 1.1 有限差分法:

有限差分法：

微分： f 'x =f x h−f x  h 0=dy

h dx

用有限的 h 代替，使得

f 'x ≈ △ y

△ x

差分的种类：

f x h−f x f x −f x −h  f x h−f x −h

一阶差分： 或者 或者

h h 2h

f x h−f x  f x −f x −h

−