CN106154058A - 一种用于天线辐射等强度曲面的计算方法 - Google Patents

一种用于天线辐射等强度曲面的计算方法 Download PDF

Info

Publication number
CN106154058A
CN106154058A CN201510128647.7A CN201510128647A CN106154058A CN 106154058 A CN106154058 A CN 106154058A CN 201510128647 A CN201510128647 A CN 201510128647A CN 106154058 A CN106154058 A CN 106154058A
Authority
CN
China
Prior art keywords
point
curved surface
equal strength
radiant intensity
aerial radiation
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN201510128647.7A
Other languages
English (en)
Other versions
CN106154058B (zh
Inventor
刘磊
周渊斐
徐开诚
李南星
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Institute of Acoustics CAS
Beijing Intellix Technologies Co Ltd
Original Assignee
Institute of Acoustics CAS
Beijing Intellix Technologies Co Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Institute of Acoustics CAS, Beijing Intellix Technologies Co Ltd filed Critical Institute of Acoustics CAS
Priority to CN201510128647.7A priority Critical patent/CN106154058B/zh
Publication of CN106154058A publication Critical patent/CN106154058A/zh
Application granted granted Critical
Publication of CN106154058B publication Critical patent/CN106154058B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Abstract

本发明提供了一种用于天线辐射等强度曲面的计算方法,该计算方法考虑到了天线辐射强度沿半径射线方向的单调递减性质,采用柱坐标表示空间计算点,准确有方向的获取辐射等强度曲面点,而且计算结果便于天线辐射等强度曲面的快速绘制,从而极大的提升了计算和绘制效率。

Description

一种用于天线辐射等强度曲面的计算方法
技术领域
本发明涉及天线辐射强度仿真计算领域,尤其涉及一种天线辐射等强度曲面的计算方法。
背景技术
随着移动网络的迅速发展,信号发射天线的安装和布置也日趋频繁,尤其是在人口密度较大的城市,需要布置大量的移动信号发射天线,而大功率的天线会产生一定的辐射。我国现行的电磁辐射防护标准(GB8702-88)为40微瓦/平方厘米,比欧美各工业化国家具有更加严格的标准,如:美国1982年颁布的标准是3000微瓦/平方厘米,欧盟现行的标准是450微瓦/平方厘米。而且,在我国只有符合安全标准的设备才能入网。因此,作为网络运营商,需要严格执行国家相关标准,在移动基站的网络规划上从用户角度出发,对基站布点进行严格的覆盖优化设计,保障基站附近居民生活安全。为达到以上目的,需要计算天线的安全辐射等强度曲面,以辅助天线的选址和安装。传统的计算等强度曲面的方法需要进行大量的仿真和计算,从而使得计算速度较慢、降低了工作效率。
发明内容
本发明的目的在于,为了解决现有的天线辐射等强度曲面的计算方法由于需要大量的仿真和计算,从而导致工作效率低的技术问题,提供一种用于天线辐射等强度曲面的计算方法,通过该计算方法能够高效率的计算并显示天线辐射等强度曲面。
为实现上述目的,本发明提出一种用于天线辐射等强度曲面的计算方法,所述的计算方法包括以下步骤:
步骤1)采用柱面坐标(H,A,R)表示空间中的计算点,由高度H和角度A确定其中一层平面上某一角度的半径射线R,并设定高度H、角度A和半径射线R的初始值均为0,设定所述高度H的最大值为maxH,设定天线辐射等强度曲面的辐射强度为P,设定高度增长步长为dH,角度增长步长为dA,半径射线增长步长为dR;
步骤2)计算柱面坐标为(H,A,R)处计算点的辐射强度,若该计算点的辐射强度小于辐射强度P,则执行步骤3),若该计算点的辐射强度大于辐射强度P,则将半径射线R按照步长dR依次递增,即R=R+dR,直至某计算点的辐射强度小于辐射强度P,则确定该计算点为天线辐射等强度曲面上的点,执行步骤4);
步骤3)将半径射线R按照步长dR依次递减,即R=R-dR,直至某计算点的辐射强度大于辐射强度P或R小于0,则确定其前一计算点为天线辐射等强度曲面上的点,执行步骤4);
步骤4)停止高度H和角度A上的半径射线R的点辐射强度计算,并将角度A按步长dA增长后执行步骤2),确定下一条半径射线重复步骤2)进行计算,直到角度A的值等于360°后停止该平面上的计算,执行步骤5);
步骤5)将高度H按步长dH增长,确定下一层平面,重新初始化角度A为0后执行步骤2),直到高度H的值大于设定的高度范围maxH后,停止所有的点辐射强度计算,执行步骤6);
步骤6)将node[i][j]、node[i][j+1]、node[i+1][j]构成的三角形作为曲面微元,得到的所有曲面微元即形成天线辐射等强度曲面,其中node[i][j]表示第i层平面上的第j个计算点,按照步骤2)、3)、4)、5)计算后,得到共m=maxH/dH层平面,每层平面上有n=360/dA个等强度点。
作为上述技术方案的进一步改进,所述柱面坐标为(H,A,R)处计算点的辐射强度计算公式表示为:
S = 100 × Pt × Gt 4 × π × d 2 μW / c m 2
其中,Pt为天线发射功率,Gt为天线增益,d为计算点距离天线的距离。
本发明的一种用于天线辐射等强度曲面的计算方法优点在于:
本方法的计算过程是利用辐射等强度曲面的空间连续性,确定辐射等强度曲面点的估计位置,从而缩短获得辐射等强度曲面点的时间;同时,本方法在计算过程中考虑到了天线辐射强度沿半径射线方向的单调递减性质,采用柱坐标表示空间计算点,可以准确有方向的获取辐射等强度曲面点,而且计算结果便于天线辐射等强度曲面的快速绘制。因此,采用本发明的天线辐射等强度曲面的计算方法极大的提升了计算和绘制效率。
附图说明
图1是本发明的一种用于天线辐射等强度曲面的计算方法流程图。
具体实施方式
下面结合附图和实施例对本发明所述的一种用于天线辐射等强度曲面的计算方法进行详细说明。
如图1所示,本发明的一种用于天线辐射等强度曲面的计算方法,包括:
步骤1)采用柱面坐标(H,A,R)表示空间中的计算点,由高度H和角度A确定其中一层平面上某一角度的半径射线R,并设定天线辐射强度计算的最大高度为maxH,设定天线辐射等强度曲面的辐射强度为P;设定高度增长步长为dH、角度增长步长为dA、半径射线增长步长为dR;设定计算点高度H、角度A和半径射线R的初始值均为0,初始化计算结果集二维数组node[m][n],其中,node[i][j]表示第i层平面上的第j个计算点,所述平面的层数m=maxH/dH,每层平面上的等强度点个数n=360/dA。
根据设定的高度H得到平行于水平面的计算平面,根据设定的角度A得到高度H中计算平面上的一条半径射线R。
步骤2)计算柱面坐标为(H,A,R)的计算点处的天线辐射强度,若为某一角度上的半径曲线首次计算,则设定R为0,否则设定R为前一计算角度上的半径射线所得等强度曲面的计算点半径,计算天线辐射在计算点(H,A,R)处的辐射强度,可采用以下公式计算获得:
采用功率通量密度表示辐射强度,其计算公式为:
S = 100 × Pt × Gt 4 × π × d 2 μW / c m 2
其中,Pt为天线发射功率,假设天线发射功率为20W;Gt为天线增益,可通过查找天线射频指标表获得,假设计算点柱面坐标为(H=20m,A=30°,R=5m),天线高度为10m,则查表获得该点天线增益为14dBi;d为辐射计算点距离天线距离,可由计算公式得到利用公式可得到该点天线辐射强度为: S = 100 × 20 × 10 14 10 4 × π × 11 . 1 2 = 32.4 μW / c m 2 .
上述的天线增益是指:在输入功率相等的条件下,实际天线与理想的辐射单元在空间同一点处所产生的信号的功率密度之比。它定量地描述一个天线把输入功率集中辐射的程度。增益显然与天线方向图有密切的关系,方向图主瓣越窄,副瓣越小,其增益越高。天线增益是用来衡量天线朝一个特定方向收发信号的能力,它是选择基站天线最重要的参数之一。一般来说,增益的提高主要依靠减小垂直面向辐射的波瓣宽度,而在水平面上保持全向的辐射性能。
根据计算所得计算点的辐射强度与设定的辐射强度P的大小关系,分两种情况进行:
若该计算点的辐射强度小于辐射强度P,则执行步骤3),若该计算点的辐射强度大于辐射强度P,则将半径射线R按照步长dR依次递增,即将R设为R+dR,计算相应点的辐射强度,直到某计算点的辐射强度小于辐射强度P,则确定该计算点为天线辐射等强度曲面上的点,执行步骤4);
步骤3)将半径射线R按照步长dR依次递减,即将R设为R-dR,直到某计算点的辐射强度大于辐射强度P或R小于0,则确定其前一计算点为天线辐射等强度曲面上的点,执行步骤4);
步骤4)停止高度H和角度A上的半径射线R的点辐射强度计算;将角度A按步长dA增长,即另A=A+dA,并判断角度A是否等于360°,若是,则执行步骤5),否则,跳至步骤2)继续计算;
步骤5)将高度H按步长dH增长,即另H=H+dH,将角度A的值设定为0,并判断高度H是否大于maxH,若是,则停止所有的点辐射强度计算,执行步骤6),否则,跳至步骤2)继续计算。
步骤6)此时结果集二维数组node[m][n]全部赋值完毕,利用该结果集绘制出辐射强度为P的天线辐射等强度曲面。具体的绘制过程为:采用以上步骤计算所得结果集中,将node[i][j]、node[i][j+1]、node[i+1][j]构成的三角形作为曲面微元,得到的所有曲面微元即形成了所述天线辐射等强度曲面。其中,总共有m=maxH/dH个平面,每个平面上有n=360/dA个等强度点,即总共有m*n个等强度点。
最后所应说明的是,以上实施例仅用以说明本发明的技术方案而非限制。尽管参照实施例对本发明进行了详细说明,本领域的普通技术人员应当理解,对本发明的技术方案进行修改或者等同替换,都不脱离本发明技术方案的精神和范围,其均应涵盖在本发明的权利要求范围当中。

Claims (2)

1.一种用于天线辐射等强度曲面的计算方法,其特征在于,包括:
步骤1)采用柱面坐标(H,A,R)表示空间中的计算点,由高度H和角度A确定其中一层平面上某一角度的半径射线R,并设定高度H、角度A和半径射线R的初始值均为0,设定所述高度H的最大值为maxH,设定天线辐射等强度曲面的辐射强度为P,设定高度增长步长为dH、角度增长步长为dA、半径射线增长步长为dR;
步骤2)计算柱面坐标为(H,A,R)处计算点的辐射强度,若该计算点的辐射强度小于辐射强度P,则执行步骤3),若该计算点的辐射强度大于辐射强度P,则将半径射线R按步长dR依次递增,直至某计算点的辐射强度小于辐射强度P,则确定该计算点为天线辐射等强度曲面上的点,执行步骤4);
步骤3)将半径射线R按步长dR依次递减,直至某计算点的辐射强度大于辐射强度P或R<0,则确定其前一计算点为天线辐射等强度曲面上的点,执行步骤4);
步骤4)停止高度H和角度A上的半径射线R的点辐射强度计算,将角度A按步长dA增长后执行步骤2),直至角度A的值等于360°执行步骤5);
步骤5)将高度H按步长dH增长,并将角度A的值设定为0后执行步骤2),直至高度H的值大于maxH停止所有的点辐射强度计算,执行步骤6);
步骤6)将node[i][j]、node[i][j+1]、node[i+1][j]构成的三角形作为曲面微元,得到的所有曲面微元即形成天线辐射等强度曲面,其中node[i][j]表示第i层平面上的第j个计算点,所述平面的层数m=maxH/dH,每层平面上的等强度点个数n=360/dA。
2.根据权利要求1所述的用于天线辐射等强度曲面的计算方法,其特征在于,所述柱面坐标为(H,A,R)处计算点的辐射强度计算公式表示为:
S = 100 &times; Pt &times; Gt 4 &times; &pi; &times; d 2 &mu;W / cm 2
其中,Pt为天线发射功率,Gt为天线增益,d为计算点距离天线的距离。
CN201510128647.7A 2015-03-23 2015-03-23 一种用于天线辐射等强度曲面的计算方法 Expired - Fee Related CN106154058B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201510128647.7A CN106154058B (zh) 2015-03-23 2015-03-23 一种用于天线辐射等强度曲面的计算方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201510128647.7A CN106154058B (zh) 2015-03-23 2015-03-23 一种用于天线辐射等强度曲面的计算方法

Publications (2)

Publication Number Publication Date
CN106154058A true CN106154058A (zh) 2016-11-23
CN106154058B CN106154058B (zh) 2018-10-02

Family

ID=58064054

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201510128647.7A Expired - Fee Related CN106154058B (zh) 2015-03-23 2015-03-23 一种用于天线辐射等强度曲面的计算方法

Country Status (1)

Country Link
CN (1) CN106154058B (zh)

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20050171706A1 (en) * 2002-06-06 2005-08-04 Mark Leckenby Method for determining field radiation levels for a radiating device
EP1870718A1 (fr) * 2006-06-23 2007-12-26 The Swatch Group Research and Development Ltd. Système de mesure du diagramme de rayonnement d'une antenne d'émission
CN101354411A (zh) * 2008-08-28 2009-01-28 深圳华为通信技术有限公司 一种测试电磁辐射强度的方法及装置
WO2009038388A1 (en) * 2007-09-19 2009-03-26 Electronics And Telecommunications Research Institute Apparatus and method for measuring antenna radiation patterns
CN103076503A (zh) * 2012-12-27 2013-05-01 广东省辐射防护协会 一种gsm移动通信基站环境电磁辐射三维空间预测方法
CN104237651A (zh) * 2014-09-18 2014-12-24 国家电网公司 移动通信基站天线辐射功率密度计算方法

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20050171706A1 (en) * 2002-06-06 2005-08-04 Mark Leckenby Method for determining field radiation levels for a radiating device
EP1870718A1 (fr) * 2006-06-23 2007-12-26 The Swatch Group Research and Development Ltd. Système de mesure du diagramme de rayonnement d'une antenne d'émission
WO2009038388A1 (en) * 2007-09-19 2009-03-26 Electronics And Telecommunications Research Institute Apparatus and method for measuring antenna radiation patterns
CN101354411A (zh) * 2008-08-28 2009-01-28 深圳华为通信技术有限公司 一种测试电磁辐射强度的方法及装置
CN103076503A (zh) * 2012-12-27 2013-05-01 广东省辐射防护协会 一种gsm移动通信基站环境电磁辐射三维空间预测方法
CN104237651A (zh) * 2014-09-18 2014-12-24 国家电网公司 移动通信基站天线辐射功率密度计算方法

Also Published As

Publication number Publication date
CN106154058B (zh) 2018-10-02

Similar Documents

Publication Publication Date Title
CN107729640A (zh) 一种使用最少阵元的稀疏天线阵列综合布阵方法
CN103323845B (zh) 一种非均匀采样综合孔径辐射计的图像反演方法
CN106772260A (zh) 基于凸优化算法的雷达阵列和差波束方向图优化方法
Huang et al. A novel 3D GBSM for mmWave MIMO channels
CN101706839B (zh) 一种基于时间反演的共形阵列天线激励相位确定方法
US11962097B2 (en) Method and apparatus for selecting antenna elements of an antenna array
CN106771673B (zh) 一种gps天线方向性测试方法及系统
CN105929388A (zh) 一种基于Wi-Fi网络的新型室内定位方法
JP6696859B2 (ja) 品質推定装置及び品質推定方法
CN107087286A (zh) 一种窗户面向通信基站的室内电磁辐射预测方法
CN107333290A (zh) 一种基站天线电性能参数监测方法
CN106842191B (zh) 一种电离层参数的获取方法
CN104038959B (zh) 一种信道模拟器及其建模方法
CN110133588A (zh) 一种天线定位方法、装置及设备
Hou et al. A hybrid method of FDTD and vector parabolic equation for radio wave propagation prediction in tunnels
CN106154225A (zh) 基于波达方向歧义消除的定位方法及装置
Jain et al. Efficient time domain HF geolocation using multiple distributed receivers
Bai et al. A 3-D wideband multi-confocal ellipsoid model for wireless MIMO communication channels
CN104683949A (zh) 一种应用于无线Mesh网中基于天线阵列的混合自定位方法
Inomata et al. Radio propagation prediction for high frequency bands using hybrid method of ray-tracing and ER model with point cloud of urban environments
CN102986152B (zh) 电磁波传播特性的分析方法和装置
CN106154058A (zh) 一种用于天线辐射等强度曲面的计算方法
CN104181529A (zh) 一种Ka波段合成孔径雷达SAR信号处理方法和设备
CN106066904A (zh) 一种基于多目标进化算法的高维多目标天线设计方法
Yin et al. Study on calculation and verification of radiowave propagation using parabolic equation for the antenna near the ground

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20181002

Termination date: 20200323

CF01 Termination of patent right due to non-payment of annual fee