电磁场实验报告 孙颖哲.docxVIP

基于有限差分法的二维边值问题的数值分析班级：1402502学号：140250218姓名：舒 荣一、实验目的1.掌握简单二维边值问题的分离变量求解方法；2.通过有限差分法的实现来熟悉数值法的求解过程。二、预习内容1.预习如下边值问题（课本149页例3.6.1）的分离变量求解过程；具体参数为：盖板电位U=100 V，其余三面电位=0 V，尺寸a=10 m，b=10 m。2.预习课本161页“3.7有限差分法”的内容。三、实验内容及步骤1.在matlab中分析基于分离变量法的解析解；2.利用简单迭代法求解，与解析法结论对比，分析求解结果的精确度。分析过程至少包括：在网格尺寸为0.1 m和1 m两种条件下，两次迭代差值最大为10-10时的分析结论； 3.利用超松弛迭代法分析，选择松弛因子，分析其对收敛速度（即迭代次数）的影响，并确定最优值。分析过程至少包括：在网格尺寸为0.1和1两种条件下，两次迭代差值最大为10-10时，迭代次数随松弛因子的变化，得到对应的最优松弛因子，与经验值（课本165页式子3.7.15）进行对比。四、实验过程1、解析解根据直角坐标系中的分离变量法求得解析解:本题中，，a=10，b=10。带入Matlab中画图计算得到解析解的图像和各点的数值作为基准：三维图像解析x方向的电势解析y方向的电势2、简单迭代法计算公式：式中：i,j=0,1,2…,k=0,1,2…,并