基于矩量法的二维金属体散射(内含matlab程序)资料.doc

基于矩量法的二维金属体散射计算 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