辽工大有限差分法实验报告讲解.doc

班 级： 姓 名： 学 号： 指导教师： 成 绩： 求解金属槽的电位分布 1.实验原理 matlab软件解决电位在金属槽中的分布2.有限差分法 方.1有限差分法原理 图1 有限差分法的网格划分导体槽中静电场的边值问题的拉普拉斯方程为： (1-1) 为简单起见，将场域分成足够小的正方形网格，网格线之间的距离为h，。节点0、1、2、3、4上的电位分别用、、、和表示。点1、点3在x0处可微，沿x方向在x0处的泰勒级数展开式为 (1-2) (1-3) 点2、点4在y0处可微，沿y方向在y0处的泰勒级数展开式为 (1-4) (1-5) 忽略高次项 (1-6) 稍作变化得到拉普拉斯方程的五点差分格式： (1-7) 可通过迭代法求解以上差分方程。 2.2有限差分法步骤 高斯—赛德尔迭代法 图2-1网络下标标示 (1-8) 进行迭代时可写为 (1-9) ,为行数，,为列数，为迭代次数，为前次迭代的结果，为当次迭代的结果，由于迭代从第一行、第一列开始，（）、（）点的迭代较（）点进行得早，顾可使用当次迭代的结果。直到所有的点电位满足（为所设定精度）时迭代停止。 .问题描述 设有一个直接地金属槽2-1所示，其侧壁与零，为V（）点位分布。 -1 金属槽 4.程序设计 4.1全场域问题 问题（）（）（）步距正方形化域。网格对应于矩阵中的元素中的值相邻两次的迭代值误差不应用matlab中的矩阵操作。ones(x,y)建立一个元素初值为矩阵。进行边界设置矩阵的左右为上下边界分别为0。在保证精度的情况下以差分式一级的数值计算一个迭代要求的矩阵A程序实现见附录A 4.2程序流程图 图4-1 程序设计流程图 contour(A)可以绘制出电位，这样可以观察出与理论情况的分布是否相同，分布图如2-3所示： 图 4.3问题 问题（）程序设计中ones(x,y)建立一个元素初值为矩阵。进行边界设置矩阵的左为左为-5:0产生的一个单位矩阵向量上下边界分别为0。在保证精度的情况下以差分式一级的数值计算一个迭代要求的矩阵A程序实现见附录BB中的程序就可以得到半场的数值解。contour(A)可以绘制出电位，分布图如下： 图 .4中心点的讨论 问题（）中心P()。求解电位的解析解该位置处于高度，该点所在的对称线上匀强电场。故： (4-1) 得： (4-2) 所以该点电位附录A的程序运行得出数值解为49.516V误差：0.4840V为：0.9700%.收敛因子作用 .1超松弛迭代法 加快收敛速度，采用迭代法。计算时将某点的值之差乘以一个，再加该的新老电位上，作为这一点的电位值迭代的表达式 (5-1) 式中α称为松弛因子，其值介于。其中收敛因子 (5-2) 其中为每边的节点数减去1。 其中M、N分别是xy两个方向的内节点数。本项目中41、N=21，计算： -1 最佳收敛因子 1.7516 迭代次数 89 取若干个收敛因子求得的迭代次数表-2和表-3中的数据 表4-2 收敛因子中心数值 1.1 1.2 1.3 1.4 1.5 迭代次数 740 615 508 416 335 262 中心点电位 49.5167 49.5168 49.5169 49.5170 49.5171 49.5171 表4-3 收敛因子中心数值 1.7 1.74 1.78 1.8 1.9 迭代次数 196 133 106 73 80 161 中心点电位 49.5172 49.5172 49.5173 49.5173 49.5173 49.5173 其中：中心点的解析解为V。 收敛因子到中心电位值可以得到一个利用E