电磁场及电磁波演示验.docVIP

实验一 基于有限差分法的二维边值问题的数值分析 一、实验目的 1.掌握简单二维边值问题的分离变量求解方法； 2.通过有限差分法的实现来熟悉数值法的求解过程。 二、预习内容 1.预习如下边值问题（课本149页例3.6.1）的分离变量求解过程； 具体参数为：盖板电位U=100V，其余三面电位=0，尺寸a=10，b=10； 2.预习课本161页“3.7有限差分法”的内容。 三、实验内容及步骤 1.在matlab中分析基于分离变量法的解析解； 2.利用简单迭代法求解，与解析法结论对比，分析求解结果的精确度。分析过程至少包括：在网格尺寸为0.1和1两种条件下，两次迭代差值最大为10-10时的分析结论； 3.利用超松弛迭代法分析，选择松弛因子，分析其对收敛速度（即迭代次数）的影响，并确定最优值。分析过程至少包括：在网格尺寸为0.1和1两种条件下，两次迭代差值最大为10-10时，松弛因子随迭代次数的变化，得到对应的最优松弛因子。 四、实验内容及步骤分析 1. 在matlab中分析基于分离变量法的解析解； 分析：利用解析法绘制图像，将二维图像分别转换为x=5和y=4的两个一维图像进行分析可知，观察y=4处的电位变化，呈左右对称，说明空间内电位是以x=5左右对称的。观察x=5处的电位变化，从y=0至y=10由低到高单调递增，说明空间内的电位也是以y=0至y=10单调递增的。 2.利用简单迭代法求解，与解析法结论对比，分析求解结果的精确度。 分析：当网格尺寸=1时，简单迭代法在y=4处，x的电位与解析法差距很大，而在x=5，y的电位与解析法近似。 分析：当网格尺寸=0.1时，迭代法一共迭代38728次，远远大于网格尺寸1时的迭代次数，而且从图像来看，与解析法的曲线基本重合，说明网格尺寸越小简单跌单越接近解析法的解。 3.利用超松弛迭代法分析，选择松弛因子，分析其对收敛速度（即迭代次数）的影响，并确定最优值。 网格尺寸=0.1 a= 1.9391 486 网格尺=1 a=1.5279 47 分析：当网格尺寸=0.1时，通过公式得出超松弛因子a=1.9391，此时迭代次数为486次。当a取值离最佳迭代因子越远时，迭代次数越多。当网格尺寸=1时，通过公式得出超松弛因子a=1.5279，此时迭代次数为47次。当a取值离最佳迭代因子越远时，迭代次数越多。五、实验结论： 1：实验1得到了解析法的图像，说明空间内的电位是关于x=5左右对称的，而且是从y=10到y=0单调递减的。 2：实验2通过网格尺寸分别为1和0.1的简单迭代发与解析法进行对比，说明网格尺寸0.1得到的电位解更接近解析解，网格尺寸越小，迭代越准确。 3：当网格尺寸=0.1时，通过公式得出超松弛因子a=1.9391，此时迭代次数为486次。当a取值离最佳迭代因子越远时，迭代次数越多。当网格尺寸=1时，通过公式得出超松弛因子a=1.5279，此时迭代次数为47次。当a取值离最佳迭代因子越远时，迭代次数越多。