有限差分法及matlab实现,有限差分法与matlab实现

《有限差分法与matlab实现》由会员分享，可在线阅读，更多相关《有限差分法与matlab实现(13页珍藏版)》请在金锄头文库上搜索。

1、 有限差分法解静电场的边值问题的算法实现及相关问题讨论： 王宁远 中国科学技术大学 09 级物理 2班 E-mail wnymail.ustc.edu.cn 摘要： 本文用MATLAB实现了有限差分法解静电场边值问题的算法，将偏微分方程的问题化为线性方程 组问题，并使用了迭代法进行线性方程组的数值解。讨论了从几个角度去优化迭代法的措施。并运用这 样的方法解决了文1的闪电模拟问题，使用了更优化的算法对重新进行了计算，并一定程度上改进了 模型,讨论了几个与文1所持的不同的观点。 正文： 经典场的边值问题在数学上表达为泊松方程和拉普拉斯方程，但解偏微分方程往往是困难的。幸而 很多时候对于具体问题我们。

2、需要的不是解析解，而是数值解，所以可以考虑用连续变量离散化的方法求 出数值解,在足够的精度上进行逼进,这就引出了有限差分法。 1.1 有限差分法: 有限差分法： 微分： f x= f xhf x h h0= dy dx 用有限的h代替，使得 f x y x 差分的种类： 一阶差分： f xhf x h 或者 f xf xh h 或者 f xhf xh 2h 二阶差分: f xhf x h f xf xh h h = f xhf xh2fx h2 设Ux,y,z为空间电势的函数。 泊松方程： 2U= 使用二阶差分代替泊松方程中的拉普拉斯算符，有： 2U= 2U x2 2U y2 2U z2 = 。

3、f xh,y,zf xh,y ,z2fx,y,z h2 表示分别对三个变元求差分之和，以下相同 矩阵(数组)是计算机中重要的数据结构，为了方便用矩阵去存储数据&#x