电磁场矩量法解金属圆柱导体散射资料

《电磁场矩量法解金属圆柱导体散射资料》由会员分享，可在线阅读，更多相关《电磁场矩量法解金属圆柱导体散射资料（8页珍藏版）》请在金锄头文库上搜索。

1、 矩量法解无限长金属圆柱导体散射场矩量法解无限长金属圆柱导体散射场 摘要摘要 矩量法（MOM）在研究电磁散射的问题中十分有效，能帮我们解决很多电 磁场的问题。本文主要对电磁场的矩量法进行了学习和研究，并利用 Matlab 进 行编程仿真，对特定问题进行求解，求解出无限长圆柱导体的雷达散射截面， 并对结果进行分析。 0 引言引言 电磁场在我们的生活中扮演者越来越重要的作用。 我们手机的信号要靠电磁场来进行传播，很多探测类的 工作都要靠电磁场的探测来进行实现。特别是在军事中， 探测飞机的雷达散射信号来获取敌方军情和设计己方 飞机不被探测到，即检测和设计隐形飞机有着十分重要 的意义。所以我们要对物体

2、对电磁场的散射情况进行分 析，主要是通过麦克斯韦方程这一基本公式来进行各种 计算。而由于实际问题十分复杂，导致解析方法难以实 现，所以需要数值计算的方法来进行实现。利用计算机 强大的计算能力，可以迅速的算出电磁场对物体的散射 情况，再对其进行分析，以得到我们所需要的场，从而 探测到物体的形状位置等重要信息。其中重要的一个环 节，就是对问题的分析和方法的应用，而矩量法在电磁 散射中有着广泛的应用，本文就将对应用矩量法来求解 无限长圆柱导体的雷达散射问题。 1 矩量法概述矩量法概述 矩量法（The Method of Moments,简称 MOM） ，在 天线、微波技术和电磁波散射等方面应用十分广

3、泛。矩 量法最早被 Richmand 和 Harrington 用于求解电磁场的 问题，而后在 Harrington 的著作中得到了系统的论述， 从此成为求解电磁场问题数值解的主要方法，并成功应 用于天线和电磁散射问题， 至今已有50年的发展历史。 矩量法是一种将连续方程离散化为代数方程组的方法。 其基本原理是：先选定基函数对未知函数进行近似展开， 带入算子方程，再选取适当的权函数，使在加权平均的 意义下方程的余量等于零，由此将连续的算子方程转换 为代数方程。剩下的问题就是利用计算机进行大量的数 字计算。原则上，矩量法可用于求解微分方程和积分方 程，但用于微分方程时所得到的的系数矩阵往往是病态

4、 的，故在电磁场中主要用于求解积分方程。矩量法是一 种严格的数值方法，求解精度高，加之格林函数直接满 足辐射条件，无需设置吸收边界条件，因而可以灵活解 决边界比较复杂的一些问题，在电磁辐射和散射、天线 电流分布、天线设计、微波网络、生物电磁学、辐射效 应研究、微带线分析、电磁兼容等方面得到广泛应用。 下面本文就将对 TM 波入射无限长金属圆柱导体的雷达 散射截面采用矩量法做进一步的计算分析。 2 问题问题分析分析 矩量法把泛函方程转化为矩阵方程，然后通过矩 阵方程求解。一般表达式会是如下的泛函方程： （1 ） 是一个函数，求解一个未知函数比较困难，我们熟悉 的是求解方程中的未知数而不是未知函数

5、，所以我们可 以将其展开为级数的形式，如下： （ 2） 这样 （1） 变为： （3） 两边乘权函数并积分得： （ 4） 可以写为： （ 5） 这样可以表示为一个矩阵形式： （ 6） 其中， Amn= m,n=1,2,N （7） Bm= m,n=1,2,N （8） 取基函数为： （ 9） 权函数为： （10 ） 可得： （11） bm= (12) 这样就可以求出c： c=A-1b （13） 则各点的场就能求解。 已知一个无限长的导体圆柱， 有 TM 波 Ei射向圆柱 体，求散射后的场分布，只考虑二维情况。 图一 示意图 TM 极化波电场满足亥姆霍兹方程： （14） 在导体圆柱表面，由边界条件有：

6、 （15） （16） 还有方程： （17） （18） 由上述方程我们可以得到： （19） （20） 即可得散射场在圆柱体表面和远处的分布。 对其进行数值化： （21） （22） （23） 可得： （24） 用计算机进行编程运算，即可求得解。 3 Matlab 编程仿真编程仿真 主要流程如下： 具体做法可见附录，有详细的文字说明。这里简要说 明一下圆柱圆周的分割问题。示意图如图所示： 主程序 几何分割 添矩阵 解矩阵方程 计算并绘图 End 图二 圆柱划分示意图 先将其分为 N 等分，求出各点的坐标，再取每段的中 点，求出新的坐标，避免端点的奇异性： （25） 每段长度间隔： （26） 参数设置

7、： 波长为 1，入射波表达式为： （27） 其中入射角为。 4 仿真结果及分析仿真结果及分析 运行程序得到了导体圆柱表面电流密度分布图 （图三） ，散射场图（图四）和总场图（图五） ，还有 散射动态图，无法展示可运行查看结果。其参数设置 都为入射角为 ，半径 r=，图片所显示的分割精度为 100 次，可自行调节分割次数，次数越多精度越高，但 速度也会越慢。 图三 TM 波入射的金属圆柱表面等效电流分布 图四 TM 波入射金属圆柱散射场分布 图五 TM 波入射金属圆柱总场分布图 图六 课本上的例图 对比仿真图与课本上的从图中我们可以看到 仿真结果基本正确，与书上的结果大致吻合。电流密 度在入射的

8、方向最大，背对入射方向最小，符合基本 原理。但是在背对入射场的方向，即角度为 0 度附 近，仿真的等效电流不均匀，与解析解有一定的偏 差，可能是精度不够，或者是奇异点问题处理的不够 好，有待进一步改进。 此题采用 Mie 级数进行求解，得到解析解为： n n jnn z KRH ej R E J 2 0 2 （28） 利用 matlab 对其进行绘图，程序见附录二，得到的等 效电流分布如图七所示： 图七 Mie 级数近似解析解 与之前采用矩量法得到的等效电流分布图出了刚刚所 说的端点处，基本完全吻合，再一次证明了所求解的 场分布的正确性。 但是当波长变短，即频率增加时，电流分布会变 得极其不稳

9、定，但改变剖分精度，又能基本稳定，如 图八和图九所示。参数分别表在图下方。 图八 =1/3.5, N=100 图九 =1/3.5, N=500 个人认为有可能是剖分精度不同，导致有些点在 计算中边为奇异点，导致解发生震荡。其中原理还没 研究太明白。 5 总结总结 通过此次矩量法的研究和学习，成功地仿真 出无限长金属圆柱导体对 TM 波散射情况的分布 图，收获很大。在刚开始接触矩量法时，完全没 有思路。后来通过大量的文献查阅和学习，逐渐 理解了矩量法的基本思路和编程方法。在仿真的 过程中，也遇到过各种问题，挑战很大，后来通 过反复看文献理解矩量法的思想， 学习 Matlab 的 编程方法才慢慢有

10、了一定的体会，最终才得到了 较为正确的结果，也有了对矩量法更为深刻的认 识。这次的仿真，不仅仅是一次完成作业，更是 一次挑战，不仅学到了课本上的知识，还学到了 Matlab 的编程技巧， 收获远不止一次完成任务这 么简单。 虽然仿真结果基本吻合，但是此次仿真还有 很多不足之处。首先是程序计算速度还不够快， 未能结合数值方法对其进行加速处理，程序有待 进一步优化。其次就是未能对导体内部共振对散 射场的影响做出讨论和分析，结果可能不是特别 准确。最后，由于时间比较紧，理论推导中的公 式比较复杂，所以没有花时间进行公示的编辑， 只是在书上进行截图，有些地方可能不清晰，有 待改善和提高。最后感谢老师这

11、一学期认真负责 的为我们讲解电磁场数值计算的知识，经常在一 些我们困惑的地方知道我们，学到了很多，感谢 李老师！ 参考文献参考文献 1Jin J. Theory and Computation of Electromagnetic FieldsM/ Theory and computation of electromagnetic fields. Wiley :, 2010. 2 佚名. 电磁场数值计算M. 高等教育出版社, 1996. 3 何红雨. 电磁场数值计算法与 MATLAB 实现J. 2004. 4 Jin J. 电磁场矩量法M. 西安电子科技大学出版社, 1998. 5 王长清.

12、现代计算电磁学基础J. 2005. 6 陈涌频, 孟敏, 方宙奇. 电磁场数值方法M. 科学出版社, 2016. 7 胡光华. 计算电磁学要论J. 2008(10):3-4. 附录一附录一 矩量法程矩量法程序序 function MOM wavelength=1; %设置入射波长 rho=1; %导体圆柱半径 Z0=377; %真空波阻抗 k0=2*pi/wavelength; %波矢 b=(Z0*k0/4); %计算要用到的常量 N=100; %导体圆柱所分段数 Z_mn=zeros(N,N); %定义阻抗矩阵 V_m=zeros(N,1); %定义入射场矩阵 width=10*rho; %

13、作图区域宽度 Nx=100; %作图划分精度 Ny=100; %*划分圆柱导体* for n=1:N+1 %分成 N 段，求每点的坐标 x(n)=rho*cos(2.0*pi*(n-1.0)/N); y(n)=rho*sin(2.0*pi*(n-1.0)/N); end for n=1:N %求两点中点坐标，避免端点奇异性 xc(n)=(x(n+1)+x(n)/2; yc(n)=(y(n+1)+y(n)/2; L(n)=sqrt(x(n+1)-x(n)2+(y(n+1)-y(n)2); %分段长度 end %*求解矩阵* for m=1:N V_m(m)=Einc(xc(m),yc(m),k0

14、); %不同点的入射场赋值 for n=1:N %带入分析后的式子 if (m=n) %考虑奇异点情况 Z_mn(m,n)=b*L(n)*(1-1j*(2/pi)*(log(0.455*k0*L(n)-1); else r=sqrt(xc(m)-xc(n)2+(yc(m)-yc(n)2); Z_mn(m,n)=b*L(n)*besselh(0,2,k0*r); end end end J=Z_mnV_m; %求解表面电流分布 figure(1) %画出电流密度在 0360 度的分布图 w=(1:1:m); plot(w*360/m,abs(J) grid on;xlabel(degrees);ylabel(Js/H0); %添加网格线,设置坐标轴 title(TM 波入射的金属圆柱表面等效电流密度分布图) %*计算并画散射场场图* A=zeros(Nx,Ny); %定义散射场矩阵 B=zeros(Nx,Ny); %定义总场矩阵 for ix=1:Nx %画一定区域内的场图 xo=(2*(ix-1)/(Nx-1)-1)*width; for iy=1:Ny yo=