基于站间距及方位角加权的邻区添加工具实现

1、基于站间距及方位角加权的邻区添加工具实现龙颖中国电信CHINATELECOM2016/8/5摘要：本课题设计一个针对小区间方位角进行距离加权的算法，从而更加贴近网络现状进行添加邻区。基于站间距及方位角加权的邻区添加工具实现新余电信2022-4-27第 页,共16页目录TOC o 1-5 h z HYPERLINK l bookmark10 o Current Document 一、设计背景2 HYPERLINK l bookmark12 o Current Document 二、设计目标及方法2 HYPERLINK l bookmark14 o Current Document 三、设计流程图

2、3 HYPERLINK l bookmark16 o Current Document 3.1设计流程图3 HYPERLINK l bookmark22 o Current Document 四、设计过程及结果44.1核心算法4 HYPERLINK l bookmark24 o Current Document 4.2检索优化算法6 HYPERLINK l bookmark30 o Current Document 4.3经纬度算两点间距离算法104.4经纬度算两点间方位角算法11 HYPERLINK l bookmark46 o Current Document 五、设计结论与收获13 HY

3、PERLINK l bookmark48 o Current Document 六、参考资料14 HYPERLINK l bookmark54 o Current Document 七、使用说明书及工具14一、设计背景目前在众多软件中有很多邻区添加工具，但是较多工具中邻区添加算法大多没有对站点间小区的覆盖范围纳入考虑，故在小区对边缘站点添加邻区时，他们将邻站小区与本站小区正朝向的优先级与小区负朝向的优先级相同，显然按照此类邻区添加算法与现网网络模型不匹配。现本课题设计一个针对小区间方位角进行距离加权的算法，从而更加贴近网络现状进行添加邻区。二、设计目标及方法主要设计目标如下：设计一个算法，使得

4、在邻区关系添加过程中，能根据方位角和地理位置区分对原小区的邻区优先级，从而更合理的添加邻区。完成设计共需要实现4个功能：小区间方位角加权的距离实现优化检索算法计算本小区与邻区间的站间距计算本小区与邻区间方位角差值本设计将采用VBA程序语言，结合EXCEL界面做成一个工具，主要在EXCEL工具内部实现上述功能。三、设计流程图3.1设计流程图本工具主要设计流程图如下：四、设计过程及结果4.1核心算法方位角加权，可以把扇区考虑成一个扇形的区域，加权系数则是从扇柄到扇形外部的距离，当加权系数大时，距离扇柄的距离大，当加权系数小时，则到扇柄的距离小。如下图所示：加权情况示意图则为了有效的将三个小区区分，

5、可以对经纬度点根据方位角进行加权，这样加权越大三个扇区的位置差异就越明显，从而能够对方位角与距离进行判断，添加更合理的网络邻区加权情况在实际网络中的应用启区方fit甬系数。m权后小区距iBfi厉小核心算法由于小区距离的方位角DAL是由A向B的，所以夹角a=方位角A-(DAL+180)方向夹角a小区距离,方夹角bAjJE駆方位角系如转换后sin(a)*D+sin(b)*D小区柜离小区距S+cos(a)*D+cos(b)*DJO权厉小区距离二KKSosa)+D+cos(b)*D)(sin(a)*D+sin(b)*D)夹角b=方位角B-DAL，算得夹角后。在垂直于小区连线方向通过SIN函数乘以加权系

6、数D得到D_SHORT，小区连线方向通过COS函数乘以加权系数D加上小区距离，得到加权后连线的方向的长度D_LONG。上述D_LONG与D_SHORT得到可以通过勾股定理，算得加权后的小区距离(见以下代码实现部分)核心代码实现部分DoWhileifield_O,tiers1-field_1,tiers2-field_2,tiers19-field_19（般20层即可）。每个对应笛卡尔层次的域将根据当前这条记录的经纬度通过笛卡尔算法计算出归属于当前层的网格，然后将gridld（网格唯一标示）以term的方式存入索引。这样每条记录关于笛卡尔0-19的域将都会有一个gridld对应起来。目前索引的实

7、现核心代码：ActiveCell.FormulaR1C1=lF(MOD(RC-5,0.1)-MOD(R-1C-5,0.1)=&(xnow-1),_Operator:=xlAnd,Criteria2:=&(ynow-1),_Operator:=xlAnd,Criteria2:=&(xnow-1),_Operator:=xlAnd,Criteria2:=&(ynow-1),_Operator:=xlAnd,Criteria2:=和两向量的夹角K。用公式*vOB=IOAI*IOBI*cosK可以得到；还有对相同点进行处理等。代码实现部分：PublicFunctionCal_Long_Lat(ByVa

8、llong1AsDouble,ByVallat1AsDouble,ByVallong2AsDouble,ByVallat2AsDouble)AsDouble单位公里ConstpiAsDouble=3.1415926535DimAngleLong1,AngleLat1,AngleLong2,AngleLat2AsDoubleAngleLong1=long1*pi/180AngleLat1=lat1*pi/180AngleLong2=long2*pi/180AngleLat2=lat2*pi/180DimsinX,cosXAsDoublesinX=Sin(AngleLat1)*Sin(AngleL

9、at2)cosX=Cos(AngleLat1)*Cos(AngleLat2)*Cos(AngleLong2-AngleLong1)X=sinX+cosXOnErrorResumeNextax=Atn(-X/Sqr(-X*X+1)+2*Atn(1)Cal_Long_Lat=6368.16*axEndFunctionA,B,C表示球面上的三个点及球面上“弧线”在该点处所夹的角a,b,c表示A,B,C三点的对“弧”两端点与地心连线所夹的角（其实这里解释成ABC三点对弧的弧度更方便）0为球心L为AB两点间球面距离首先：在知道AB点经纬度后，要用到第一个公式，三面角余弦公式，cos(c)=cos(a)X

10、cos(b)+sin(a)Xsin(b)Xcos(AOCMB)AOCB是面AOC与面BOC的二面角，为了方便，写成这个样子这里将已知数据代入，公式便写成：cos(c)-cos(90-Bw)Xcos(90-Aw)-i-sin(90-Bw)Xsin(90-Aw)Xcos(Bj-Aj)二面角A0CB的度数就是两点经度之差其次：知道了角c的余弦值后要求得它的正弦值，所用的公式就是三角函数公式里最基本的“扣方加赛方等于1”的一个变形sin(c)=J1cos2(c)+jVrVrVr-rY*第三步：求得正弦后，接下来要用一个不太常用的公式，球面正弦公式abcsiii(A)sin.(Bsin(C)将已知数据代

11、入并稍微变形一下，公式写为：涮A)二刎匹3)用反正弦函数求角度，于是上式可直接写成Er働严賊闷sin这里需要注意一点，一开始的假设便是求B点相对于A点的方位角，因此这里是Bj-Aj，不要写反，否则得不到正确结果。在此得到的结果并不总符合对方位角的定义，因此要根据B相对于A的位置在四个象限两个轴上进行讨论，依据不同情况对计算结果进行不同处理。假设A点固定于原点，则:B点在第一象限，Azimuth=A；B在第二象限，Azimuth=360+A;B在第三四象限，Azimuth=180-A。经过一系列计算，最后，就得到了最终结果。代码实现部分：PublicFunctionqiufangweijiao(

12、ByVallon1AsDouble,ByVallat1AsDouble,ByVallon2AsDouble,ByVallat2AsDouble)AsIntegerConstpiAsDouble=3.1415926535DimavlatAsDoubleDimqiufangweijiaosAsDoubleavlat=(lat1+lat2)/2Iflat1-lat2=0Thenqiufangweijiaos=90Elseqiufangweijiaos=180-Atn(lon1-lon2)*Cos(avlat)/(lat1-lat2)*180/piEndIfIflat1lat2Thenqiufangw

13、eijiaos=qiufangweijiaos+180EndIfIfqiufangweijiaos=360Thenqiufangweijiaos=-360+qiufangweijiaosEndIfqiufangweijiao=Round(qiufangweijiaos,0)EndFunction五、设计结论与收获目前该工具实现了基于方位角加权邻区的添加功能，并邻区代码实现。较之前的邻区添加工具相比，更加贴近网络情况。在设计过程中熟悉了VBA些循环算法及三角函数命令，最重要的是加强了设计队伍对于扇区结构的认识，把扇区的结构量化具体化。本次的设计从数量、距离两个方面对邻区添加功能进行实现，现网网络

14、复杂，不同场景下无线环境不一，后续可以结合网络指标情况及地理环境情况对邻区进行优化，希望本工具能对后续一些工具或软件的开发提供灵感及借鉴，为网络优化及现代化建设出功出力。感谢新余电信及华为公司对本次设计的大力支持！六、参考资料1、基于Solr的LBS（地理位置搜索）实现原理 HYPERLINK /a221133/article/details/14525197 /a221133/article/details/145251972、根据两点的经纬度求方位角和距离，等 HYPERLINK /s/blog_658a93570101hynw.html /s/blog_658a93570101hynw.

15、html七、使用说明书及工具7.1使用说明书1首先在现网工参SHEET中贴入工参，包括：eNodeBID、CelllD、SectorlD、Longitude、Latitude、站名、Azimuth。以上所有字段都是必填项，请选择性粘贴为数值。如下图所示：ABcHE|rG1aKodaBTI)CfiLlIDSfictorlBLangituiaLfl.titud.6可占名A.zimutk491115.093T5ZT.9ZT153XYHFDI渝朮河T龙ti140:i47144650115.093T5ST.92T153KTHFD：痂水疽卞东j耘肚HIti24DL471446513115.0937527.

16、92713XHFD：谕办诃下工代毂2H：【h34D471391491IIS.1082T.8576XYHfTL.分宜凤阳乡玦2FY3J:.4T1391502115.1002T.0576XYHFTL艸育風阳各$FZFY130471301513115.IOS27.8576XYHFTL分亩凤汨岭V240:：47142249111亍10939河RAH3XYHFTI.分亘启飯机榕乐畠岚羊展喘趋帀FP；T.T.04T14Z250Z115.10309ZT.06193XYHFTL分宜高岚机房眾高岚乡展鹏超市陀GLL1Z00471422513115.1093927.0C133KITTL分宜高岚机尽馮岚乡辰鵬起床胆

17、九L27014713QQ4Q1115.112B2?7的4诃XYHFT1.分官国11弋商山村炸円2Wa:4T139950Z115.uzazr.Qb454XIHfTL分山*壮比3tiITO34T1399513115.1130227.054B4KTHF_TL分甘铃山买岗lllHtFEQS.ft24D14713Q5401115.OSQIE-27.8E805XYHFTL分肓铃1|新申由$F2QW10b471395SO2115.0591627.85895KYHTTL.分宜铃Ll新祉tF2QSti_一LI-*ri1I14.-,13J2在规划首页贴入需要规划的数据，包括：eNodeBID、CellID、Sec

18、torID、Longitude、Latitude、站名、Azimuth。以上所有字段都是必填项，请选择性粘贴为数值。如下图所示：ABDE7-eNodeEIDCellIDSetj.yilJLongitudeLatltiide站名AzLuuth471295502115.0782127.B7092KY_HF_U新余市二医K_M5AtKlK(；_W120471295513115.078212?.B?092XY_1F_TL_新余市二医院_共联餘1X（：_/24=0+71359491113.0732127.B7052K_HF_U新余市芬直电厂?ILCI_L30471359502I:.27.B5312H_-

19、1_新帝十匸口丄!一：；亲：；：烏：.604L35?513LL5.临二27.B53L2KYHFTI新余币分宜电J$际1XSLISO471230491115.06845327.B5312XYHFTL新余H工行城东办共联期1閉下280471230SO2115-135327.05924:l3!【二=m1:141工？G1:：2：X?H?j5ADDEUTEABINnimEQIICELL:LOC虹CELLIh1:、【二J门1471296工：?H?jJ丄二|znr.dm.femiu二让応口豆二二匕1：!【二=i11!XY_HF_DY.1-JThU.A-h：!(.；|II：L|.H.izl/l：=1|z：l1

20、471420二3阳H?DrADDEUTEAJIIHJWWMCELL:LQCALCEL1ID=1.IKC=41工占J1471397：；旺芒ADDEUTRiJIINBAfEEQHCELL:LOCALCE二二匸1.：=TJ1二：:TF：i5匸JLUTF.kli丁二.fUIU二:Li：CAL：Ii二二i-1:I：-=1门1zr14i：ADDEBIEAJIIHTRmEQUCELL:L0CrtLCEL1II1:HOIt汇加tj1汽忆2XY_HF_DY.ADDBUTEABIN1WKB3IICELL:LOCALCELLII1:IKOit;71曲1J二F二:?HF.DYADDZITEhl：丁二丄邛Z：QU二：L疋灶了二二lIJ【二J;7詁;：i1尸；:；：.HT工;a?jTrTE.i,r7YATET；:ir-:.:：iii-.um；i=i：r-=y-11-：1：；H-!1-：IThM.-b-：-l(-.:ln：t|.H.il.l：.=、lt7144bJ1J二-