静电场问题的离化模型.pptVIP

概述 1.模型背景 2.有限差分法 二维泊松方程的差分格式 Error as a function of the mesh size. 算例 5.误差分析和改进措施 1、误差分析：误差来源分为三种：替代误差、跳跃误差、舍入误差。 替代误差：差商代替微商引入的误差。  由h高阶舍入引起！ 替代误差是本次模拟的误差产生的主要因素！ 2、改进措施： 1）、减小h 2)、合理设置网格线。 附录 源程序 e=ones(36,42); a=zeros(36,42); for i=24:26 for j=1:12 a(i,j)=1; end end x1=a ; %进行迭代 while norm(e)1e-6%差量限制在10^(-6) for i=2:23 for j=2:12 a(i,j)=0.25*[a(i-1,j)+a(i+1,j)+a(i,j-1)+a(i,j+1)]; end end for i=27:35 for j=2:12 a(i,j)=0.25*[a(i-1,j)+a(i+1,j)+a(i,j-1)+a(i,j+1)]; end end 程序设计原理 启动 ↓ 给定边值Φ1，Φ2 ↓ 填写场域内Φ的初值 ↓ 迭代次数计数N =0 ↓ ------→ N=N+1 ↓ ↑ 以普通迭代法进行一次迭代，求Φi,j(N+1) 否 ↓ ←---所有内点相邻两次迭代的方均根值是否小于W ↓是 输出结果 注：1、普通迭代法即用a(i,j)=0.25*[a(i-1,j)+a(i+1,j)+a(i,j-1)+a(i,j+1)]对矩阵顺序运算最终得到方程组的解。 2、齐次诺曼条件的使用：在原有矩阵多加一列作为平衡列！即在迭代过程中始终使a(i,j-1)=a(i,j+1) 其中第j列所在的位置为对称面！ * 静电场问题离散化模型 及求解 混合班第20组：冀军 李甲子 吴联芳 1.模型背景 2.对于同种介质中静电场的求解（二维泊松方程的差分格式） 3.对于不同种介质静电场的求解（基于通量守恒的差分格式） 4.屏蔽微带传输线的静电场的模拟应用 5.误差分析和改进措施 静电场的求解是解决电磁问题的关键，通常我们希望建立基于麦克斯韦方程的连续模型来求解，但实际上静电场的性质并不象想象中那样易于描述，即使可以用连续模型建立的式子来描述，求解可能非常困难。我们希望通过建立静电场的离散模型，把微分变成差分，来近似描述静电场的电位分布情况，更重要的是这样将使得求解变得容易，并且通过数字化的软件处理可以用图象来模拟静电场的实际电位分布。我们本模型的目的正是要寻求建立这种离散模型的方法以及如何求解这种离散模型。 有限差分法（ FDM ——Finite Difference Method）是基于差分原理的一种数值计算法。其基本思想：将场域离散为许多小网格，应用差分原理，将连续函数的微分采用差分近似，从而将求解连续函数的泊松方程的问题转换为求解网格节点上的差分方程组的问题。 Notice : finite [’fainait]，infinite [’infinit] 有限差分的网格分割