二维矩量法的matlab仿真,基于HALLEN方程法的半波天线的参数计算

摘  要： 在现代电磁工程中，对于边界不复杂的问题，可用解析法得到精准解，但较复杂的边值问题用解析法不能得到解答，需要数值法。在此导出HALLEN方程法，用MATLAB软件实现，它能够非常接近解析值，说明其正确性。矩量法是求解电磁场边界值问题中一种行之有效的数值方法，利用矩量法将HALLEN法的积分方程化为差分方程，将积分方程中积分化为有限求和，从而建立代数方程组，用MATLAB编程进行数值计算和仿真。最后得出结论，矩量法分析能够接近其解析解，说明其正确性。

关键词： 矩量法；HALLEN方法；半波天线；MATLAB

0 引言

自赫兹和马可尼发明了天线以来，天线在社会生活中的重要性与日俱增，如今已成不可或缺之势。天线是人们见闻世界的耳目，是人类与太空联系的纽带，是文明社会的组成要素[1]。

矩量法是求解电磁场边界值问题中一种行之有效的数值方法，其数学本质是一种求解线性方程的方法[2]。矩量法[3-4]所做的工作是将积分方程化为差分方程，或将积分方程中的积分化为有限求和，从而建立代数方程组，故它的主要工作是用计算机求解代数方程组，因而在电磁场数值计算方面得到了广泛的应用，并已得到了国内外学者的公认。HALLEN E在1938年提出积分方程法[5-6]，首先从假定的电流分布用矢量磁位写出天线表面外部的电场公式，然后令此电场等于天线表面内部因克服集肤效应引起的阻抗所需要的电场强度，由此求得表示电流分布的积分方程式。国内外对半波天线的研究主要集中于激励原理、特性、实验等方面的讨论，而关于天线方向图和电流分布等参数计算和分析很少。张黎明[7]、姜光兴[8]等人利用矩量法对半波天线的参数进行了计算。HAFIANE A等人[9]利用天线阵模型对天线的参数进行了计算。基于以上的研究，本文利用HALLEN方法计算并仿真出半波天线的电流分布、功率增益和方向图，计算过程中不但简化了数值分析的复杂性，而且减小了计算量并提高了其精度。

1 天线的关键参数推导

1.1 矩量法的数学表述[7，10]

设一个非其次方程为：

L(f)=g(1)

其中，L为微分算子或积分算子，f为待求的场量或响应，g为已知的源或激励。

设空间为线性的，在算子L的定义域内选择一组函数{fn}，将待求函数f展开为它们的组合：