电磁场实验报告(一)

电磁场实验报告(一)

ID:46825454

大小:311.25 KB

页数:6页

时间:2019-11-28

电磁场实验报告(一)_第1页
电磁场实验报告(一)_第2页
电磁场实验报告(一)_第3页
电磁场实验报告(一)_第4页
电磁场实验报告(一)_第5页
资源描述:

《电磁场实验报告(一)》由会员上传分享,免费在线阅读,更多相关内容在行业资料-天天文库

1、电磁场实验报告(一)实验目的(1)掌握有限差分法的原理与步骤。(2)理解并掌握求解差分方程组的高斯迭代法和超松弛迭代法。(3)分析超松弛迭代法中加速收敛因子α的作用。(4)学习应用有限差分法求解接地金属槽问题,并编制计算程序。实验原理有限差分法是以差分原理为基础的一种数值计算方法。应用有限差分法通常所采取的的步骤是:(1)采用一定的网格分割方式对求解场域离散化。(2)进行差分离散化处理。用离散的、只含有限个未知数的差分方程组,来近似代替场域内具有连续变量的偏微分方程以及边界条件。(3)结合选定的代数方程组的解法,编制计算机程序,求解由上面所得对应于待

2、求边值问题的差分方程组,所得解答即为该边值问题的数值解。启动给定边值给定场域内φ的初值迭代次数n=0n=n+1按超松弛迭代法进行一次迭代所有内点相邻两次迭代值是否小于W?是否打印迭代次数停机程序框图如下:实验内容MATLAB程序如下:a=zeros(21,41);a(1,:)=100;a(21,:)=0;a(:,1)=0;a(:,41)=0;prea=zeros(21,41);n=0;pi=1.74;while(sum(sum(abs(a-prea)<10^-3))~=861)prea=a;n=n+1;fori=2:20forj=2:40a(i,j)

3、=prea(i,j)+pi/4*(prea(i+1,j)+prea(i,j+1)+a(i-1,j)+a(i,j-1)-4*prea(i,j));endendenda(:,20)contour3(a)mesh(a)结果如下:100.000094.082988.191882.351976.585870.913965.352859.915454.610849.443944.416039.524834.764630.127325.602121.176416.835812.56488.34684.16440n=72绘图如下:仿真程序如下:0014.104372

4、1350766928.28422081581828312.5395460422248416.8703478142963521.2766261320327625.7583809954341730.3156124045004834.9483203592317939.65650485962791044.54404046901121149.63079358441521254.81951177486261360.11019504035331465.50284338088731570.99745679646461676.59403528708521782.292

5、57885274921888.09308749345651993.99556120920720100.000000000001实验分析从上述实验过程中我们可以发现通过matlab的数值计算得到的数据与仿真实验得到的数据误差较小,因此可以认为实验过程是正确的。并且通过此次实验我们可以直观的体会到有限差分法在求解电磁场域的数值解中的应用、matlab和maxwell等工具在求解数值解和仿真实验中的应用,并且对进一步优化实验过程有了更深刻的体会。

当前文档最多预览五页,下载文档查看全文

此文档下载收益归作者所有

当前文档最多预览五页,下载文档查看全文
温馨提示:
1. 部分包含数学公式或PPT动画的文件,查看预览时可能会显示错乱或异常,文件下载后无此问题,请放心下载。
2. 本文档由用户上传,版权归属用户,天天文库负责整理代发布。如果您对本文档版权有争议请及时联系客服。
3. 下载前请仔细阅读文档内容,确认文档内容符合您的需求后进行下载,若出现内容与标题不符可向本站投诉处理。
4. 下载文档时可能由于网络波动等原因无法下载或下载错误,付费完成后未能成功下载的用户请联系客服处理。