《线天线辐射及散射(共12页).doc》由会员分享，可在线阅读，更多相关《线天线辐射及散射(共12页).doc（12页珍藏版）》请在taowenge.com淘文阁网|工程机械CAD图纸|机械工程制图|CAD装配图下载|SolidWorks_CaTia_CAD_UG_PROE_设计图分享下载上搜索。
1、精选优质文档-倾情为你奉上摘要：矩量法是将连续方程离散为代数方程组的方法，此方法对于求解微分方程和积分方程均适用，本文以半波振子天线为例，系统的阐述了半波振子天线的海伦积分方程的建立，利用矩量法求解海伦积分方程而得半波振子天线上的电流分布，并进一步根据电流肺部，求解半波振子天线的方向图。关键字：半波振子天线；海伦积分方程；矩量法1 引言电磁辐射和散射问题的分析方法一般可以分为两大类，即解析方法和数值方法，而实际中只有极少数集合形状特别的电磁问题才能通过解析方法求解，大部分只能用数值方法获得近似解。矩量法则是求解微分方程和积分方程的一种重要的数值分析方法，它从函数空间和线性算子的观点来处理问题，
2、具有计算效率高、处理灵活、快速而精准、不限定物体几何形状、理论基础较健全等诸多优点，因而在电磁场数值计算方面得到了广泛的应用。同时，由于线天线在实际中应用广泛，而且是分析其它天线的基础。在矩量法分析过程中，有多种不同的积分方程可供选择，如双位积分方程、Hallen积分方程、Pocklington积分方程、Schelkunoff积分方程、响应积分方程等等，而前3种应用最为广泛，因此本文采用Hallen积分方程对半波振子天线进行深入的数值分析，并用matlab编程仿真。2 矩量法的基本原理2.1 矩量法的概念矩量法是将一个算子方程化为矩阵方程，然后求解该矩阵方程的方法。在历史上把采用基函数和检验函
3、数离散化的积分方程的数值方法称为矩量法。矩量法是一种基于泛函分析理论的积分形式的数值方法，这种方法具有计算结果准确且误差小、处理过程灵活、分析目标不限定物体几何形状和理论基础健壮等优点。如果非齐次方程为 （2-1） 式中是线性算子，为已知函数，为未知函数。令： （2-2）式中是系数。被称为展开函数或基函数。用算子的线性便可以得到: （2-3）规定了一个适当的内积，那么在的值域内定义一个权函数或检验函数的集合，并对每个取式(2-3)的内积，则（2-4）式中1，2，3。此方程组可以写成如下的矩阵形式 （2-5）式中（2-6） （2-7）如果矩阵是非奇异性的，于是可写成： （2-8）2.2基函数和权
4、函数的选择 矩量法求解算子方程的关键问题是基函数和权函数的选择。在特定的问题中,主要任务是选择基函数和权函数,它们必须是线性无关的。基函数可以分为整域基和分域基;权函数一般有点匹配法,伽略金法,最小二乘法等。应用矩量法需注意:误差分析;方程收敛性;积分奇异点处理等。下面分别介绍基函数和权函数的选择：本文的基函数选择的是正弦函数，选取这样的基函数是考虑到它可以满足细导体末端电流为零的边界条件。文中用点匹配法选择权函数。下面将简要介绍点匹配法。若选取狄拉克函数为权函数，即令 （2-9）式中，狄拉克函数的定义为（2-10）由于，因而矩量法方程(2-6)( 2-7)相应矩阵元素的计算结果为（2-11）
5、 （2-12）由此表明和的计算归结为只需计算所在点处的对应值，因此称这种方法为点匹配法。3 半波振子天线的矩量法求解假设半波振子天线放在Z轴线上，原点位于中点。如下图所示：图1半波振子天线示意图图1为的半波振子天下，可做如下假设：（1）电流沿导线轴流动，体电流密度可以用线电流I来近似，体电流密度用线电荷密度近似；（2）忽略天线端面的周向电流径向电流；（3）天线上电流仅为长度变量的函数，即当半波振子天线上在加上馈电信号时，或在接收电场作用时，其上将产生电流，其电流分布按基函数展开，在此，我们将选用正弦函数作为基函数。 通过以上假设可以知道，半波振子天线的矢量磁位，而 （3-1）其中为格林函数，通
6、过麦克劳林级数展开，舍弃高次项，这样其稳定度也较高。的表达式为： （3-2）为波数。于是z分量的电场表达式为： （3-3）由边界条件可知：（由边界条件可知，在导体表面切向电场为零），其中为电源所产生的场，令，为电压，为狄拉克函数，将上述条件和（3-1），（3-2）式代入（3-3）于是便得到： （3-4）对上式进行求解为： （3-5）公式中B为常数，为自由空间的波阻抗。该方程为半波振子天线的海伦积分方程（Hallen）。海伦积分方程简单变型为： （3-6）选择分布电流展开表达式： （3-7） （3-8）为待定系数，这里为简单起见，只取到N=2。将（3-7）式代入（3-6）式，并整理，可得： （3
7、-9）式中 ， ， 式中为天线的一般长度，为导体半径，为导体表面上场点的坐标，为导体表面上源点的坐标，在轴线上。（3-9）式中有三个未知量，C，因此应选三个加权函数作为三个方程式。课题中我们求解的是=/4，为波长，为了求半波振子天线上电流分布，采用点配法，权函数选择狄拉克函数，即，其中m=1,2N，其中选择的点为：变换，可以选择的值为=0，=/8，=/4，即以这三点为观察点，分别对分成99段的天线进行观察。运用选择的的值对每一对（3-9）式两边求内积，这样可以将（3-9）式转变成矩阵的形式 (4)根据的取值，等号右边的内积计算如下：并根据函数的特性，上式可化简为 (5)式中其中除了，三个未知量
8、，其余各个元素都可以通过计算机编程来求解，然后再求解出，这样根据公式（2）就可以求得天线上的电流分布。4应用举例这里采用半波振子天线为例：波长为，长度为。把天线分成99（奇数）段，基函数采用三角函数，权函数采用函数。其他条件为，则半波振子部分电流值，电流分布和方向图分别如图所示：图2 部分分段点的电流值图3 电流分布示意图图4 E面方向性图图5 H面方向性图图6 电流分布计算数据比较本文用的基函数是正弦函数，与基函数为狄拉克函数的程序相比，推算过程有其特殊之处，故不能用普通方法验证。但通过图3与图6相比可以看出，其电流分布，并与king的三项理论及Mack实验数据曲线比较基本符合，仿真结果正确
9、且较好。附录：程序如下所示：clear all;close all;clf;tic; % 计时N=99;% 分段数e0=8.854e-012;u0=4*pi*10(-7);c=3e+008;w=2*pi*c;%光速,角频率lamda=1; % 波长measurement =7.022*(10)(-3); % 物体电尺寸v0=1; % 电压常数point=3; % 待求点数range=1/4; % 求解的天线归一化长度D=(N-1)/2+1;% D指中间段k=(2*pi)/lamda; % 波数a=measurement*lamda; % 导体半径l=range*lamda; % 天线长度len
10、= l /N;% 将线分成奇数段,注意首末两端的电流为0step=l/(point-1); % 点匹配间距i_point=1:point; % 测试点赋值LL(i_point)=step*(i_point-1); z=linspace(-l,l,N+1); % 积分离散% 求解A向量for i_point=1:point; r=(LL(i_point)-z).2 + a2).(1/2); % 场源距离离散 g=exp(-j*k*r)./r; % 格林函数离散 A(i_point)=trapz(z,cos(k*z).*g); % A元素确定end% 求解B向量for i_point=1:poin
11、t; r=(LL(i_point)-z).2 + a2).(1/2); % 场源距离离散 g=exp(-j*k*r)./r; % 格林函数离散 B(i_point)=trapz(z,sin(2*k*abs(z).*g); % B元素确定end% 求解C向量for i_point=1:point; C(i_point)=cos(k*LL(i_point); % B元素确定end% 阻抗矩阵确定ZZ=A.,B.,C.;% 电压矩阵确定for i_point=1:point; VV(i_point)=(-j*v0/60)*sin(k*abs(LL(i_point); % B元素确定 VV= 0;-j
12、/(60*sqrt(2);-j/60;end% 求解a1,a2,C_contant;current=inv(ZZ)*VV;% 图形表示z_distribute=linspace(0,l, N+1);current_function=current(1,1)*sin(k*(l-abs(z_distribute) +. current(2,1)*sin(2*k*(l-abs(z_distribute); % 离散化电流分布% 电流实虚部current_re=real(current_function);current_im=imag(current_function);% current_valu
13、e=double(current_function)此处加这句可查看具体电流值% 绘图figure(1)plot(current_re,z_distribute,r);hold on;plot(current_im,z_distribute,g);xlabel(current distribution);ylabel(unitary distance)title(antenna current distribution plot);legend(real current,imag current,2);% E面方向图figure(2)tt=1;for theta=0:pi/100:2*pi;f
14、or n=1: point F(n)=current_function(n+1)*len*exp(i*k*len*n*cos(theta)*sin(theta);end;FF(tt)=sum(F);tt=tt+1;end;FFmax=max(abs(FF);for ii=1:tt-1 FF(ii)=FF(ii)/FFmax;end;theta=0:pi/100:2*pi;polar(theta,abs(FF),r);title(E-plane pattern,F(theta,phi);title(E-plane pattern,F(theta,phi)legend(E面方向性图,天线长L=0.
15、25,天线分段N=99,天线半径a= 0.)%view(90,78)view(90,-68)% H面方向图figure(3) theta=pi/2;tt=1;for fai=0:pi/180:2*pifor ii=1: point FH(ii)=current_function(ii)*len*exp(i*k*len*ii*cos(theta)*sin(theta);end;FHA(tt)=sum(FH);tt=tt+1;end;FHAMAX=FHA/max(FHA);fai=0:pi/180:2*pi;polar(fai,abs(FHAMAX),r);title(H-plane patter
16、n,F(theta,phi),theta=90);legend(归一化的H面方向性图,天线长L=0.25,天线分段N=99,天线半径a= 0.)参考文献1 R.F.哈林顿. 计算电磁学的矩量法.北京：国防工业出版社,1981.2 李世智. 电磁辐射与散射问题的矩量法.北京：电子工业出版社,1985.3 伍刚.半波振子天线矩量法的研究.微计算机信息.20064 梁荣江.基于Matlab分析矩量法在线天线中的应用.西安邮电学院学报.2009.015 姜光兴，曹伟，朱洪波.基于不同积分方程的线天线矩量法分析.南京邮电学院学报.2005.026 高宝春，彭戈，胡一峰.基于矩量法的线天线电磁散射和电磁辐射分析.中国舰船研究.2011.067 吴晓英，刘学观.线天线辐射特性的矩量法分析.苏州大学学报（工科版）.2005.068 RECollin，FJZucker”Antenna Theory”(1969)专心-专注-专业