原创力文档- 1、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。。
- 2、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击 这里二次下载。
- 3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
- 4、该文档为VIP文档，如果想要下载，成为VIP会员后，下载免费。
- 5、成为VIP后，下载本文档将扣除1次下载权益。下载后，不支持退款、换文档。如有疑问请联系我们。
- 6、成为VIP后，您将拥有八大权益，权益包括：VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
- 7、VIP文档为合作方或网友上传，每下载1次， 网站将根据用户上传文档的质量评分、类型等，对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
基于矩量法地二维金属体散射(内含matlab程序)
基于矩量法的二维金属体散射计算
1 问题的描述
本题是用矩量法计算二维金属圆柱体的散射场，如图所示为一圆柱体和一个椭圆柱的截面，为了计算简单，选入射波为垂直z轴入射的TM或TE平面波
y
2
x
2 矩量法求解过程
2.1 电场积分方程
2.1.1问题的分析
由麦克斯韦方程组
(1)
(2)
可得电场积分方程为
(3)
表示在圆柱表面的面电流在远处产生的总场。设入射场为E，散射场为E，由金属表面的边界条件
=0 (4)
得 (5)
2.1.2 离散化
设入射波为,将散射体截面C分为N份△C，用点匹配法对上述积分式子进行离散化， 即基函数可取
(6)
可得下列离散方程：
[P]{J}={b} (7)
其中： (8)
(9)
当m≠n时，
(10)
当m=n时 解析积分为
(11)
其中=1.781,e=2.718
2.1.3方程组的求解
可用LU分解求解方程组，即P=LU ，其中P为可逆矩阵，L为上三角矩阵，U为下三角矩阵，则可利用这两个基本的三角矩阵进行求解J，求出J之后，就可求散射场
(12)
(13)
与二维场中的散射截面
(14)
2.1.4输出结果的验证
此散射问题也可用模式展开法进行求解，可用此结果对本问题进行验证。所得J为
(15)
2.2 磁场积分方程
对于TE波垂直与z方向入射时的金属体的散射。对于一般的TE波而言只有场分量，电流密度方程只有横向分量。则MFIE为： (16)
其中
(17) y
y t
n
x
x
(18)
其中t表示边界上的一点，是和X的夹角。
根根据前面的过程，圆柱边界分成N分。等效电流密度可以近似为一些脉冲函数的迭加：
(19)
其中 (20)
则得到
(21)
矩阵非对角元
(22)
(23)
在上认为是常量故
(24)
对角元 (25)
其中
回波宽度的近似公式为：
(26)
3 计算机数值实验及分析
本论文通过数值计算验证前面理论分析的结果，并对数值计算结果进行分析。分别以金属圆柱体和金属椭圆体为计算例子，做数值实验和分析。所使用的计算机程序是商业软件MATLAB6.5，数值实验在本人机子（celeron4 1.8
文档评论（0）