CN103778477B - 监视台站布站方法和装置 - Google Patents
监视台站布站方法和装置 Download PDFInfo
- Publication number
- CN103778477B CN103778477B CN201310591118.1A CN201310591118A CN103778477B CN 103778477 B CN103778477 B CN 103778477B CN 201310591118 A CN201310591118 A CN 201310591118A CN 103778477 B CN103778477 B CN 103778477B
- Authority
- CN
- China
- Prior art keywords
- parameter
- individual
- station
- describing
- covering
- 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.)
- Active
Links
Abstract
本发明涉及空管领域，具体而言，涉及监视台站布站方法和装置。其中，监视台站布站方法，包括：对期望以监视台站覆盖服务的每个对象分别建立一个覆盖目标参数模型；对覆盖目标参数模型进行计算以设定布站区域；对布站区域进行站点优化布网分析，以确定监视台站布站方案。本发明提供的监视台站布站方法和装置，与现有技术中相比，避免了人工在野外或在地图上采集候选点，以及不断地人工选取与计算，而是采用计算机分析设置候选点，以及通过计算机分析来缺点最优站点布局，所以结果分析快速高效，大幅降低了成本。
Description
技术领域
本发明涉及空管领域，具体而言，涉及监视台站布站方法和装置。
背景技术
空管监视台站的种类和数量都在逐年增加，分布也越来越广泛。监视台站组网规划主要目的是在指定类型的空管服务区域内，科学规划一定数量的监视台站，使监视台站覆盖范围最大化覆盖服务区域。现有的组网规划方法大都采用半人工半自动化的方法，即先人工在野外或在地图上采集候选点，再利用相关方法计算指定类型的设备在该位置的覆盖范围结果，并将结果与规划区域地图进行叠加，通过不断的人工选取与计算，最终确定最佳候选点。
然而发明人发现，该方法候选点的选取对技术人员的技术和经验要求较高，因此人工成本较高，导致增加了组网规划成本。
发明内容
本发明的目的在于提供监视台站布站方法和装置，以解决上述的问题。
在本发明的实施例中提供了一种监视台站布站方法，包括：对对期望以监视台站覆盖服务的每个对象分别建立一个覆盖目标参数模型；对覆盖目标参数模型进行计算以设定布站区域；对布站区域进行站点优化布网分析，以确定监视台站布站方案。
在本发明的实施例中提供了一种监视台站布站装置，包括：模型创建模块，用于对期望以监视台站覆盖服务的每个对象分别建立一个覆盖目标参数模型；布站区域设定模块，用于对覆盖目标参数模型进行计算以设定布站区域；优化布网分析模块，用于对布站区域进行站点优化布网分析，以确定监视台站布站方案。
本发明实施例提供的监视台站布站方法和装置，与现有技术中相比，避免了人工在野外或在地图上采集候选点，以及不断地人工选取与计算，而是采用计算机分析设置候选点，以及通过计算机分析来确定最优站点布局，所以结果分析快速高效，大幅降低了成本。
附图说明
图1示出了根据本发明实施例的监视台站布站方法的流程图；
图2示出了根据本发明优选实施例的监视台站布站方法的流程图；
图3示出了根据本发明优选实施例的二次雷达信号覆盖模型示意图；
图4示出了根据本发明优选实施例的站点优化布网分析算法流程图；
图5示出了根据本发明优选实施例的邻近台站之间顶空盲区补盲示意图；
图6示出了根据本发明优选实施例的约束性判断与修正流程图；
图7示出了根据本发明实施例的监视台站布站装置的示意图。
具体实施方式
下面通过具体的实施例子并结合附图对本发明做进一步的详细描述。
图1示出了根据本发明实施例的监视台站布站方法的流程图，包括：
步骤S10，对期望以监视台站覆盖服务的每个对象分别建立一个覆盖目标参数模型；
步骤S20，对覆盖目标参数模型进行计算以设定布站区域；
步骤S30，对布站区域进行站点优化布网分析，以确定监视台站布站方案。
本实施例的监视台站布站方法，通过规范统一多源监视台站组网规划流程，组网规划所需的资源，从而为多源监视台站组网规划自动化奠定基础。通过自动化过程，本实施例避免了人工在野外或在地图上采集候选点，以及不断地人工选取与计算，而是采用计算机分析设置候选点，以及通过计算机分析来确定最优站点布局，所以结果分析快速高效，大幅降低了成本。
优选地，对期望以监视台站覆盖服务的每个对象分别建立一个覆盖目标参数模型包括设置以下至少一个参数：
覆盖目标高度参数，整型数组，用于描述监视台站规划覆盖高度；布站数量参数，整型，用于描述规划布网台站个数；终端区覆盖重数参数，整型，用于描述规划设备台站多重覆盖度要求，最低满足一重覆盖；覆盖终端区参数，浮点型二维数组，用于描述布网站点覆盖的区域，坐标设定格式为[R1,R2,..Rm],其中各元素为[x1,y1,x2,y2,...,xn,yn](n≥3)，m≥0,当m＝0时，表示当前组网规划与覆盖区域无关；覆盖航路参数，浮点型二维数组，用于描述覆盖航路坐标设定，格式为[L1,L2,..Lm],其中各元素为[x1,y1,x2,y2,...,xn,yn](n≥2)，m≥0,当m＝0时，表示当前组网规划与覆盖航路无关；航路覆盖重数参数，整型，用于描述规划台站对航路覆盖重数要求；布网区域参数，浮点型数组，用于描述布网区域坐标格式[x1,y1,x2,y2,...,xn,yn](n≥3)；航路覆盖宽度参数，整型，用于描述为航路左右延展的宽度；航路合法布站宽度参数，整型，用于描述为航路附近合法的布站区域宽度；实际覆盖搜索控制参数，布尔型，用于描述当前组网规划是否考虑现有监视台站覆盖；实际覆盖搜寻区域参数，浮点型数组，用于描述当前组网规划考虑的现有设备必须在该搜索区域内，其坐标格式为[x1,y1,x2,y2,...,xn,yn](n≥3)。
优选地，对覆盖目标参数模型进行计算以设定布站区域包括：
读取覆盖目标参数模型的配置；判断覆盖目标参数模型中的布网区域参数是否为空；如果布网区域参数为空，则判断覆盖目标参数模型中的覆盖航路参数是否为空，如果不为空，则以覆盖目标参数模型中的航路覆盖宽度、航路合法布站宽度为缓冲半径分别通过缓冲区分析得到航路覆盖面集合、航路合法布站区集合，以航路合法布站区集合进行合并，得到的合并区域作为新参数值赋给布站区域；如果覆盖航路参数为空，则将覆盖目标参数模型中的覆盖终端区参数的值赋给布站区域；如果布网区域参数不为空，则判断覆盖目标参数模型中的覆盖航路参数是否为空，如果不为空，则以覆盖目标参数模型中的航路覆盖宽度、航路合法布站宽度为缓冲半径分别通过缓冲区分析得到航路覆盖面集合、航路合法布站区集合，将航路合法布站区集合的合并区域与覆盖目标参数模型中的布网区域参数取交集，将交集结果赋给布站区域；如果覆盖航路参数为空，则将覆盖目标参数模型中的布网区域参数的值赋给布站区域。
优选地，采用遗传算法对布站区域进行站点优化布网分析。
优选地，确定监视台站布站方案包括：当站点优化布网分析满足终止条件，则根据需求输出指定个数的适应度值最高的最优个体方案，每个方案输出内容包括组网站点坐标集合，站点所在覆盖结果集。
图2示出了根据本发明优选实施例的监视台站布站方法的流程图，包括：
步骤S202，地理数据准备；
步骤S204，运行参数预处理；
步骤S206，布站区域设定；
步骤S208，运行数据预处理；
步骤S210，站点优化布网分析；
步骤S212，优化部署方案输出。
本优选实施例中引入候选站点参数模型、覆盖目标参数模型作为多约束参数条件；其中，监视设备类型包括一次雷达，二次雷达，广播式自动监视ADS-B，多点定位MLAT，每一种设备类型对应一种候选站点参数模型，该参数模型将规划台站布站约束条件进行规范化，保证了规划目标约束条件正确初始化。覆盖目标参数模型是对布网台站覆盖服务的对象进行约束。后三步是监视台站组网分析与结果输出验证步骤，先根据参考覆盖目标参数模型、候选站点参数模型设置值，确定选址布网区域，然后随机生成初始布站群方案，通过多约束条件判定，符合条件的布站群方案进入遗传算法进行适应度计算、选择、交叉和变异等，最终得到多组优化布站方案。通过对优化方案中站点地址进行微调，最终形成台站组网规划方案供决策部门参考。
优选地，对每种类型的监视台站分别建立一个候选站点参数模型包括设置以下至少一个参数：
设备类型参数，字符串型，用于描述监视台站的类型；天线高度参数，浮点型，用于描述监视台站的规划设备天线的高度；设备仰角参数，浮点型，用于描述监视台站的设备仰角；塔台限制高度参数，浮点型，用于描述监视台站的塔台设置的高度限制；地形作用距离参数，浮点型，用于描述地形对监视台站的信号覆盖影响范围的距离限制；低空盲区高度参数，整型，用于描述监视台站的低空不能覆盖的高度限制；高空盲区高度参数，整型，用于描述监视台站的高空不能覆盖的高度限制；建筑物干扰源参数，字符串类型数组，用于描述监视台站的建筑物图层数据集名称数组；建筑设施限制半径参数，浮点型，用于描述监视台站的不应有金属建筑物、密集的居民楼、高压输电线的限制半径；有源干扰源参数，字符串类型数组，用于描述监视台站的有源信号干扰源；有源干扰设施半径参数，浮点型，用于描述监视台站的指定半径范围内，不应有能产生有源干扰的电气设施。
方案详细说明如下：
步骤S202，地理数据准备
地理数据是指带空间数据和属性数据的各类地图或图集，是台站组网规划必须条件，主要包括等高线数据(大比例尺，如1:25万)、数字高程DEM数据(大比例尺，如1:25万)、地面交通数据、居民区数据、各类建筑物干扰源数据集，各类有源干扰源分布数据、现有四类监视设备台站分布数据，其中建筑物干扰源数据有居民区、金属建筑物、高压输电线等数据，有源干扰源数据包括气象雷达、高频炉等分布数据。以上为常见地理数据，也可以更加实际需要增加其他地理数据，这些添加的地理数据与覆盖目标参数模型和候选站点参数模型中的参数相关。
步骤S204，运行参数预处理
参数预处理是对候选站点参数模型、覆盖目标参数模型和遗传算法参数模型进行参数配置。三个参数模型中设置了当前规划监视站点的台站建设约束条件，只有满足这些条件的监视台站，才能符合规划组网要求。
候选站点参数模型中参数为：设备类型，天线高度，设备仰角，建筑设施限制半径，有源干扰设施半径等。具体说明参考表1。该模型为组网设备内部约束条件，模型参数值为该类型设备中国民用航空行业标准参考值，具有通用性，也可以根据实际行业应用做对应调整。
表1候选站点参数模型
覆盖目标参数模型中参数为：覆盖目标高度，布站数量，终端区覆盖重数，覆盖终端区，覆盖航路，航路覆盖重数，布网区域，航路覆盖宽度，航路合法布站宽度，实际覆盖搜索控制，实际覆盖搜索区域，具体说明参见表2，该模型为组网设备外部约束条件，其中参数中覆盖终端区与覆盖航路2个参数配置不能同时为空；模型参数值为当前规划设备覆盖服务对象的限制条件，配置值只对当前规划组网有效。
表2覆盖目标参数模型
遗传算法参数模型是遗传算法的运行参数，参数为：种群规模，个体编码长度，编码精度，终止迭代数，适应度变化阈，交叉率，变异率，最优个体保存率，具体参数说明参见表3，该参数模型参数值一般采用经验值。
表3遗传算法参数模型
步骤S206，布网区域设定
布网区域是规划台站布站有效范围限定，该区域范围受到上述中的覆盖目标参数模型配置和航路覆盖预处理结果影响，通过该步骤，最终确定实际布网区域RealPlanRegion。设定方法如下：
(1)首先读取覆盖目标参数模型配置；判断布网区域参数值是否为空，如果该参数为空，则执行第(2)步；否则执行第(3)步。
(2)判断覆盖航路是否为空，如果不为空，则将步骤S208得到的航路合法布站区集合进行合并，得到的合并区域作为新参数值赋给RealPlanRegion。如果为空，则将覆盖终端区参数值赋给RealPlanRegion。
(3)判断覆盖航路是否为空，如果不为空，则将步骤S208得到的航路合法布站区集合进行合并，得到的合并区域与布网区域参数取交集，交集结果作为新参数值赋给RealPlanRegion。如果为空，则RealPlanRegion为配置参数布网区域的值。
布网区域的重新设定，为台站布网缩小搜索范围。
步骤S208，运行数据预处理：航路覆盖预处理和设备覆盖预处理。
步骤S206与步骤S208可以交换执行秩序。
航路覆盖预处理和设备覆盖预处理是站点优化布网分析前的数据准备，该两个处理过程是否执行，受到覆盖目标参数模型中参数配置值的控制。
(1)航路覆盖预处理：
根据配置的覆盖目标参数模型，如果覆盖航路参数配置不为空，则以配置的覆盖航路为对象，然后分别以航路覆盖宽度为缓冲半径通过缓冲区分析得到航路覆盖面集合，以航路合法布站宽度为缓冲半径通过缓冲分析得到航路合法布站区集合。
(2)设备覆盖预处理：
设备覆盖预处理实现在指定监视设备类型的信号覆盖模型基础上计算监视设备信号在某高度上的覆盖。以二次雷达为例，根据2.2中配置的模型参数，如果实际覆盖搜索控制参数配置为TRUE，则查询地理数据中二次雷达数据集中所有落在指定实际覆盖搜索区域内二次雷达站点，并计算这些站点在覆盖目标高度上的信号覆盖。因此该步骤关键在于信号覆盖模型的建立。以下为监视设备通用信号覆盖模型建立的说明：
图3示出了根据本发明优选实施例的二次雷达信号覆盖模型示意图。如图3，假设当前监视设备站点其在一个方位上发现目标的仰角范围是0°～45°(该仰角范围根据设备出厂设置为准)，考虑遮蔽物的影响，对于某一高度层的空域被划分为顶空盲区、可视区和遮蔽区。考虑一阶球面近似可以得到遮蔽角θs计算公式：
其中ae为等效地球半径，hs为设备视线上最高遮蔽点的海拔高度，ls为遮蔽点到天线的斜距，ha为天线海拔高度(即设备台站的海拔标高+天线高度)。在图中所示方位，海拔高度hr的空域可视区在监视设备仰角[θs，45°]范围内。
在遮蔽角θs下，监视设备天线到高度hr的视在距离：
在遮蔽角θs下，监视设备的实际作用距离为：
其中Rmax是监视设备的最大作用距离，典型的最大作用距离约为360km；G(θ)是监视设备天线的归一化方向函数，由监视设备特性决定，不同的类型，G(θ)是不同的。
监视设备不同高度层上可视区(覆盖范围)计算步骤：
(1)选择空域高度hr；
(2)计算不同方位扇区(0～360°)的遮蔽角，初始记方位ai＝0(0≦i<360)；
(3)计算以设备为起点，地形作用距离为长度，方位ai上相交的等高线，其中该长度与等高线的交点为地形遮蔽点；
(4)根据上步求出的等高线的地理信息，找出高程值最大且距站点最近的等高线L，并获得L上的遮蔽点坐标；
(5)计算该设备站点到上述遮蔽点的视在距离Rs(θs)；
(6)如果Rs(θs)≤Rr(θs)，则设置高度hr上可视区范围[θs，45°]；
(7)如果Rs(θs)＞Rr(θs)，令Rs(θn)＝Rr(θn)，求得θn，则设置高度hr上可视区范围[θn，45°]；
(8)以设备位置为中心，由球面坐标公式和直角坐标公式可以计算在ai方位角度上可视区角度θn上的作用距离在平面上相对设备位置的偏移量，从而得到ai方位上水平覆盖坐标点；
(9)i＝i+1，ai＝ai+1，循环执行(a.2)～(a.8)，直到i＝360；
(10)分别将0～360°方位上获取的水平覆盖坐标点组成多边形区域，该区域则是设备信号水平覆盖区域。
步骤S210，站点优化布网分析
站点优化布网分析是在前面步骤的基础上，采用遗传算法思想，解决布网区域内布站最优化的问题。遗传算法是一种基于生物自然选择与遗传机理的随机搜索算法，一种有效的解最优化问题的方法。
本方案结合遗传算法和GIS(Geography Information System，地理信息系统)来解决复杂的空间优化配置问题，具有智能的搜索算法可以大大提高空间的搜索能力。在基于进化的优化过程中，根据GIS的空间数据来计算不同解决方案(染色体)的适应度。
图4示出了根据本发明优选实施例的站点优化布网分析算法流程图。如图4所示，以监视设备二次雷达为例，站点优化布网分析的步骤是：
步骤S302，根据设定的布网区域，初始布站方案群体；
群体方案中的每个个体可以表示为A：(F,Pts,Code,P),F为个体适应度，Pts表示初始化个体中站点的坐标数组[(x1,y1),(x2,y2),...,(xn,yn)]，Code为染色体串(经纬度坐标二进制编码)，P为个体适应度概率。
通过步骤S206中得到的RealPlanRegion区域可能是规则或不规则区域得到该区域的最小外接矩形，从而得到布网站点的地理坐标范围，经度[Xmin,Xmax],纬度[Ymin,Ymax]。
设雷达网由N部雷达单元组成，以网内各雷达单元的地理位置坐标为染色体，采用二进制编码方式，即各雷达单元的地理位置分别用Lx,Ly位二进制编码表示，N部雷达单元的染色体串可表示为：
Code＝N(Lx+Ly)位编码形成染色体串对应着一种雷达网部署方案。
Lx或Ly位长度可以根据遗传算法公式：
x＝xLxL-1xL-2...x2x1 公式2
计算得到，其中公式1为二进制编码精度公式，公式2为一个个体的编码，公式3为解码公式，Umax，Umin分别是随机数值的最大值、最小值，L为二进制编码长度，δ为二进制编码的编码精度。根据上述三个公式可以很容易进行编码和解码。
步骤S304，进行多约束条件判定和群体修正；
多约束条件判定是为了保障初始化站点方案群或当前方案群体中的所有个体都要满足布站基本要求。根据候选站点参数模型中配置条件进行站点过滤。多约束条件判定条件可以分为环境约束条件和性能约束条件。
环境约束条件包括布网点RealPlanRegion实际布网区域；候选站点参数模型中的建筑设施限制半径，有源干扰设施半径，即以候选站点为中心，在配置半径范围内如果有建筑设施或有源干扰设施，则该站点不符合条件。
性能约束条件是指监视设备台站建设应满足相邻单元之间最小部署距离dmin，相邻单元最大部署距离dmax等条件。
dmin是为了避免设备之间近距离的信号干扰而要求的最小部署距离，可以根据现有公式得到：PT+GT-PL+GR＝Pmin，PL＝K+20lg(f)+20lg(dmin)
其中：PT＝发射机功率(dBm)，GT＝发射机天线增益(dB)，PL＝路径损耗，GR＝接收机天线旁瓣增益(dB)，f＝频率(MHz)，K＝32.45(距离以km为单位)
则：
dmax是为了防止由于相邻单元设备之间的距离太远而产生顶空盲区而要求的最大部署距离，这样设备之间才具有一定的高空和顶空补盲作用。邻近站间最大距离由以下的两个参数决定：低空盲区高度Hi和顶空盲区高度Hl，如图5所示。具体的计算方法可根据监视设备天线的方向函数和地球曲率半径求得不同Hl和Hh要求下相邻设备站点的最大距离dmax(Hi)和dmax(Hh)，则最终相邻雷达的最大部署距离表示为：dmax＝min[dmax(Hi),dmax(Hh)]，如要求不超过3km和Hh不超过6km的条件下，计算得到相邻雷达间的距离应当小于230km。
如图6所示，种群方案多约束条件判定与群体修正具体详细操作如下：
步骤S602，输入种群方案；
步骤S604，对种群中每个个体的单元进行约束条件(1)判断；
步骤S606，对不满足的个体单元进行局部修正，即在原有坐标的基础上加入随机项将该单元移出干扰区域，直到所有个体单元都满足约束条件(1)；
步骤S608，对种群中每个个体进行约束条件(2)判断，如果群体中存在不满足约束条件(2)的个体，执行以下的步骤610-616；如果所有个体都满足约束条件(2)，则执行步骤S618；
步骤S610，以个体中某个单元坐标(x1,y1)为基准；
步骤S612，随机产生[0,2π]上均匀分布的随机数θ和[dmin,dmax]上均匀分布随机数k；
步骤S614，对不满足约束条件(2)的个体，相邻单元的位置调整为(x2,y2)＝(x1,y1)+k(cosθ,sinθ)；
步骤S616，对修正后的个体单元进行约束条件(1)判断，如果不满足，则返回执行步骤S612；直到所有的个体都满足约束条件(1)，则执行步骤S618；
步骤S618，输出修正后的种群；修正后的种群保证了所有个体都同时满足了约束条件(1)和(2)。
其中，条件(1)为环境约束条件判定，条件(2)为性能约束条件判定。
步骤S306，进行覆盖适应度计算和规划站点覆盖计算；
遗传算法对一个个体(解)的好坏用适应度函数值来评价，适应度函数值越大，解的质量越好。适应度函数是遗传算法进化过程的驱动力，也是进行自然选择的唯一标准，它的设计应结合求解问题本身的要求而定。
本方案根据监视台站布站条件，设计适合本方案的覆盖适应度函数，具体说明如下：
假设系统由N部监视设备(D1D2…DN)组成，航路覆盖区为Sh，终端区Sz，M为高度层数(一般分为4个高度层分别是3600m、4800m、6000m、8400m)，Sij为监视台站i在第j高度层的覆盖区域(与不同设备的覆盖能力及地形因素相关)。ωj为加权系数，ωj取决于对各个高度层的关心程度，取值在[0，1]上，组网覆盖性能达到最优要求满足以下两个条件：
1)管制区覆盖范围最大化
管制空域内覆盖范围越大，则覆盖连续性和严密性越好。可用覆盖系数ρ描述，取值范围为[0,1]，即：
其中Sg为管制空域面积，包括终端区Sz和航路覆盖区Sh。
2)不同管制空域覆盖层数要求
如终端区域要求覆盖重数为Nzc＝3，航路要求覆盖层Nhc＝2，用冗余系数μ描述，取值[0，1]，μz为终端区覆盖冗余系数，μh为航路覆盖冗余系数：
μ＝μz+μh 公式7
其中为覆盖冗余因子；
设置f＝k1ρ+k2μ 公式8
其中k1,k2，是权重系数，表示终端区覆盖与航路覆盖的重要程度，可以通过层次分析方法得到。这些都是统计值，与经验值有关
在上述公式4-公式8的基础上，具体覆盖适应度计算操作过程如下：
(I)令当前迭代次数为t，第一次初始化t＝1，t<T；
(II)根据候选站点设备类型，分别计算Q个群体方案中所有候选站点在指定的地理位置上的信号覆盖：
(III)判定实际覆盖搜索控制，如果该参数为TRUE，覆盖适应度需要考虑现有台站覆盖，执行(IV)；否则执行(V)；
(IV)通过步骤S206可以得知覆盖条件的现有台站个数为N1，规划台站数配置值为N2，则公式5中的设备个数为N＝N1+N2；
(V)如果没有现有台站参与规划，则N＝N2；
(VI)判定覆盖目标模型参数配置，确定此次规划是正对终端区、航路或二者。如果正对终端区(即覆盖终端区域参数不为空，覆盖航路参数为空)，则执行第(ⅥI)；如果正对航路，则执行(ⅦI)；如果正对二者，则执行(IX)。
(VII)结合公式4～8，令管制空域面积Sg＝Sz，μh＝0，μ＝μz；
(VIII)结合公式4～8，令管制空域面积Sg＝Sh，μz＝0，μ＝μh；
(IX)结合公式4～8，令管制空域面积Sg＝Sz∪Sh，μ＝μz+μh；
(X)以i＝0，循环计算群体方案Q个个体的适应度Fi＝fi＝k1ρi+k2μi；其中当Nzc>＝N,冗余系数公式5和6中的冗余因子为0.
其中，方案中单个候选站点信号覆盖计算方法如下：
(i)利用DEM高程数据获取候选点Pt所在位置的高程值；
(ii)计算候选点Pt的天线海拔高度＝高程值+塔台限制高度；
(iii)根据当前配置的设备类型，利用该设备类型归一化方向函数计算该候选站点在指定高度上的信号覆盖。
上面的地形数据DEM和等高线可以通过预先准备地理数据来获取。
地理数据是指带空间数据的各类地图或图集，是台站组网规划必须条件，主要包括等高线(大比例尺，如1:25万)、数字高程DEM图(大比例尺，如1:25万)、各类空管区分布图集、航路图、地面交通图、居民区图、干扰源设备分布图、现有四类监视设备台站分布图；其中空管区分布图集包括管制区、情报区、管制扇区、进近管制区、进近管制扇区、终端区等区域，这些区域都可能成为监视台站布站区域。以上为常见地理数据，也可以更加实际需要增加其他地理数据，这些添加的地理数据与覆盖目标参数模型和候选站点参数模型中的参数相关。
步骤S308，终止判断；
GA遗传算法是一种反复迭代的随机搜索算法，它通过多次进化逐步逼近最优解，因此需要确定终止条件。这里采用规定遗传代数和个体适应度判断相结合的方法作为终止条件。当遗传代数大于设定值T或群体适应度平均值变化小于门限β时，则终止搜索。操作如下：
(I)设全局变量Favg为群体适应度平均值，Fcavg为当前群体适应度平均值，计算
(II)判断当前迭代数t，当t＝1时，Favg＝Fcavg；执行第S310步；否则执行第Ⅲ步。
(III)当t<T且Favg-Fcavg<＝β，或t＝T时，终止迭代，执行第S212步输出最优方案。否则执行第S310步。
步骤S310，存储最佳个体
(I)对群体的所有个体按适应度大小Fi进行降序排序。
(II)如果当前迭代数t＝1时，令全局变量ColonyGroup表示迄今为止的适应度最高最好个体，大小为G·Q；否则通过适应度将ColonyGroup数组中最优个体与当前群体中的G·Q个适应度最高的个体进行比较，将适应度最高的G·Q个个体保存到ColonyGroup数组中。也就是当前群体中最佳个体的适应度比总的迄今为止的最好个体的适应度还要高，则以当前群体中的最佳个体作为新的迄今为止的最好个体。
(III)将迄今为止的ColonyGroup数组中最好的个体替换当前群体中G·Q个适应度最差的个体。
(IV)重新按照适应度对群体个体排序。
(V)计算当前群体中每个个体适应度概率各个个体分配到的概率作为其能够被遗传到下一代的概率，基于这些概率值用比例选择的方法产生下一代群体。
步骤S312，选择交叉变异；
执行经典遗传算法产生下一代新个体。经典遗传算法包括选择、交叉、变异。
(I)选择：即从当前群体中按照一定的概率选出优良的个体，使它们有机会作为父代繁殖下一代子孙形成新群体，新群体个体个数为Q。
(II)交叉：将选择得到的新群体个体进行两两相互配对，配对的染色体依据交叉概率Pc按照某种方式交互其部分基因，从而形成两个新的个体，将新个体保存到新的群体中。
(III)变异：针对新群体个体通过一个变异率(pn)改变染色体位串上的某一个或某一些基因，将变异后的个体添入新群体中。
(IV)分别将新群体中的Q个个体染色体进行解码，将解码后的地理坐标数组赋给Pts。
(V)至此一次站点优化布网分析结束，t＝t+1，进入第S304步，执行多约束条件判定。
步骤S212，优化部署方案输出
当站点优化布网分析满足终止条件，则根据需求输出指定个数的适应度值最高的最优个体方案。每个方案输出内容包括组网站点坐标集合，站点所在覆盖结果集。
图7示出了根据本发明实施例的监视台站布站装置，包括：
模型创建模块10，用于对期望以监视台站覆盖服务的每个对象分别建立一个覆盖目标参数模型；
布站区域设定模块20，用于对覆盖目标参数模型进行计算以设定布站区域；
优化布网分析模块30，用于对布站区域进行站点优化布网分析，以确定监视台站布站方案。
本实施例避免了人工在野外或在地图上采集候选点，以及不断地人工选取与计算，而是采用计算机分析设置候选点，以及通过计算机分析来确定最优站点布局，所以结果分析快速高效，大幅降低了成本。
优选地，优化布网分析模块采用遗传算法对布站区域进行站点优化布网分析。
从以上的描述中可以看出，本发明具有以下优点：
1、结果分析快速高效
2、充分考虑了民用航空监视台站选址规范和标准，使分析结果根据符合实际；
3、方法灵活性大，适用于各类监视设备类型。
显然，本领域的技术人员应该明白，上述的本发明的各模块或各步骤可以用通用的计算装置来实现，它们可以集中在单个的计算装置上，或者分布在多个计算装置所组成的网络上，可选地，它们可以用计算装置可执行的程序代码来实现，从而，可以将它们存储在存储装置中由计算装置来执行，或者将它们分别制作成各个集成电路模块，或者将它们中的多个模块或步骤制作成单个集成电路模块来实现。这样，本发明不限制于任何特定的硬件和软件结合。
以上所述仅为本发明的优选实施例而已，并不用于限制本发明，对于本领域的技术人员来说，本发明可以有各种更改和变化。凡在本发明的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。
Claims (6)
1.一种监视台站布站方法，其特征在于，包括：
对期望以所述监视台站覆盖服务的每个对象分别建立一个覆盖目标参数模型；
对所述覆盖目标参数模型进行计算以设定布站区域；
对所述布站区域进行站点优化布网分析，以确定监视台站布站方案，其中，采用遗传算法对所述布站区域进行站点优化布网分析；
所述遗传算法采用了遗传算法参数模型，所述遗传算法参数模型包括以下至少一个参数：
种群规模参数，整型，用于描述群体中所含个体的数量，一般取20～100；
个体编码长度参数，整型，用于描述个体编码长度，同时决定了编码精度；
编码精度参数，浮点型，用于描述求解台站位置坐标小数精度；
终止迭代数参数，整型，用于描述传算法终止进化代数；
适应度变化阈参数，浮点型，用于描述根据实际问题的求解过程得到的经验值；
交叉率参数，浮点型，用于描述进化过程中交叉操作的概率；
变异率参数，浮点型，用于描述进化过程中变异操作的概率；
最优个体保存率参数，浮点型，用于描述产生新一代群体时保留上一代最佳个体的百分率；
所述遗传算法采用了遗传算法参数模型，所述遗传算法参数模型包括以下至少一个参数：
种群规模参数，整型，用于描述群体中所含个体的数量，一般取20～100；
个体编码长度参数，整型，用于描述个体编码长度，同时决定了编码精度；
编码精度参数，浮点型，用于描述求解台站位置坐标小数精度；
终止迭代数参数，整型，用于描述传算法终止进化代数；
适应度变化阈参数，浮点型，用于描述根据实际问题的求解过程得到的经验值；
交叉率参数，浮点型，用于描述进化过程中交叉操作的概率；
变异率参数，浮点型，用于描述进化过程中变异操作的概率；
最优个体保存率参数，浮点型，用于描述产生新一代群体时保留上一代最佳个体的百分率；
采用遗传算法对所述布站区域进行站点优化布网分析包括：
(a)确定所述监视台站的信号覆盖区域，具体包括：
设置遮蔽角
其中ae为等效地球半径，hs为设备视线上最高遮蔽点的海拔高度，ls为遮蔽点到天线的斜距，ha为天线海拔高度；
设置在遮蔽角θs下，所述监视台站天线到高度hr的视在距离
设置在遮蔽角θs下，所述监视台站的实际作用距离
其中Rmax是所述监视台站的最大作用距离，G(θ)是所述监视台站天线的归一化方向函数；
设置所述监视台站不同高度层上的可视区，步骤如下：
(a.1)选择空域高度hr；
(a.2)计算不同方位扇区的遮蔽角，初始记方位ai＝0，其中0≦i<360；
(a.3)计算以所述监视台站为起点，地形作用距离为长度，方位ai上相交的等高线，其中该长度与等高线的交点为地形遮蔽点；
(a.4)根据上步求出的等高线的地理信息，找出高程值最大且距站点最近的等高线L，并获得L上的遮蔽点坐标；
(a.5)计算所述监视台站到上述遮蔽点的视在距离Rs(θs)；
(a.6)如果Rs(θs)≤Rr(θs)，则设置高度hr上可视区范围[θs，45°]；
(a.7)如果Rs(θs)＞Rr(θs)，令Rs(θn)＝Rr(θn)，求得θn，则设置高度hr上可视区范围[θn，45°]；
(a.8)以所述监视台站为中心，由球面坐标公式和直角坐标公式计算在ai方位角度上可视区角度θn上的作用距离在平面上相对设备位置的偏移量，从而得到ai方位上水平覆盖坐标点；
(a.9)i＝i+1，ai＝ai+1，循环执行(a.2)～(a.8)，直到ai＝360；
(a.10)分别将0～360°方位上获取的水平覆盖坐标点组成多边形区域，设置该区域为所述监视台站的信号覆盖区域；
(b)设置初始布站方案群体，具体包括：
布站方案群体中的每个个体表示为A：(F,Pts,Code,P)，F为个体适应度，Pts表示初始化个体中站点的坐标数组[(x1,y1),(x2,y2),...,(xn,yn)]，Code为对应经纬度坐标二进制编码的染色体串，P为个体适应度概率；
确定所述布站区域的最小外接矩形的经度[Xmin,Xmax]和纬度[Ymin,Ymax]；
设雷达网由N部雷达单元组成，以网内各雷达单元的地理位置坐标为染色体，各雷达单元的地理位置分别用Lx,Ly位二进制编码表示，N部雷达单元的染色体串表示为：
Code＝N(Lx+Ly)位编码形成染色体串对应着一种雷达网部署方案；
Lx或Ly位长度根据遗传算法公式计算得到，公式如下：
x＝xLxL-1xL-2...x2x1 公式2
其中公式1为二进制编码精度公式，公式2为一个个体的编码，公式3为解码公式，Umax，Umin分别是随机数值的最大值、最小值，L为二进制编码长度，δ为二进制编码的编码精度；
(c)进行多约束条件判定和群体修正，具体包括：
多约束条件判定条件分为环境约束条件和性能约束条件，环境约束条件包括实际布网区域；候选站点参数模型中的建筑设施限制半径，有源干扰设施半径；性能约束条件包括所述监视台站应满足的相邻单元之间最小部署距离dmin和相邻单元最大部署距离dmax，设置：
dmax＝min[dmax(Hi),dmax(Hh)]；
其中：PT＝发射机功率(dBm)，GT＝发射机天线增益(dB)，PL＝路径损耗，GR＝接收机天线旁瓣增益(dB)，f＝频率(MHz)，K＝32.45(距离以km为单位)；
(c.1)输入种群方案；
(c.2)对种群中每个个体的单元进行环境约束条件判断，对不满足的个体单元，在原有坐标的基础上加入随机项将该单元移出干扰区域，直到所有个体单元都满足环境约束条件；
(c.3)对种群中每个个体进行性能约束条件判断，如果群体中存在不满足性能约束条件的个体，执行步骤(c.4)-(c.6)；如果所有个体都满足性能约束条件，则执行步骤(c.7)；
(c.4)随机产生[0,2π]上均匀分布的随机数θ和[dmin,dmax]上均匀分布随机数k；
(c.5)对不满足性能约束条件的个体，以个体中某个单元坐标(x1,y1)为基准，相邻单元的位置调整为(x2,y2)＝(x1,y1)+k(cosθ,sinθ)；
(c.6)对修正后的个体单元进行环境约束条件判断，如果不满足，则返回执行步骤(c.3)；直到所有的个体都满足环境约束条件，则执行步骤(c.7)；
(c.7)输出修正后的种群；
(d)确定覆盖适应度和规划站点覆盖，具体包括：
假设系统由N部所述监视台站(D1D2...DN)组成，航路覆盖区为Sh，终端区Sz，M为高度层数，Sij为监视台站i在第j高度层的所述信号覆盖区域，ωj为加权系数，ωj取决于对各个高度层的关心程度，取值在[0，1]上，组网覆盖性能达到最优要求满足以下两个条件：
管制区覆盖范围最大化
设置覆盖系数ρ，取值范围为[0,1]，计算如下：
其中Sg为管制空域面积，包括终端区Sz和航路覆盖区Sh；
不同管制空域覆盖层数要求
设置冗余系数μ，取值[0，1]，μz为终端区覆盖冗余系数，μh为航路覆盖冗余系数，计算如下：
μ＝μz+μh 公式7
其中为覆盖冗余因子；
设置f＝k1ρ+k2μ 公式8
其中k1，k2，是权重系数，表示终端区覆盖与航路覆盖的重要程度；
具体覆盖适应度计算操作过程如下：
(d.1)令当前迭代次数为t，第一次初始化t＝1，t<T；
(d.2)根据候选站点设备类型，分别计算Q个群体方案中所有候选站点在指定的地理位置上的信号覆盖：
(d.3)判定实际覆盖搜索控制，如果该参数为TRUE，执行步骤(d.4)；否则执行步骤(d.5)；
(d.4)如果已知覆盖条件的现有台站个数为N1，规划台站数配置值为N2，则设置公式5中的设备个数为N＝N1+N2；
(d.5)如果没有现有台站参与规划，则设置N＝N2；
(d.6)判定覆盖目标模型参数配置，确定此次规划是正对终端区、航路或二者，如果正对终端区，则执行步骤(d.7)；如果正对航路，则执行步骤(d.8)；如果正对二者，则执行步骤(d.9)；
(d.7)结合公式4～8，令管制空域面积Sg＝Sz，μh＝0，μ＝μz；
(d.8)结合公式4～8，令管制空域面积Sg＝Sh，μz＝0，μ＝μh；
(d.9)结合公式4～8，令管制空域面积Sg＝Sz∪Sh，μ＝μz+μh；
(d.10)以i＝0，循环计算群体方案Q个个体的适应度Fi＝fi＝k1ρi+k2μi；其中当Nzc>＝N，冗余系数公式5和6中的冗余因子为0，Nzc为终端区域要求覆盖重数；
其中，方案中单个候选站点信号覆盖计算方法如下：
(1)利用DEM高程数据获取候选点Pt所在位置的高程值；
(2)计算候选点Pt的天线海拔高度＝高程值+塔台限制高度；
(3)根据当前配置的设备类型，利用该设备类型归一化方向函数计算该候选站点在指定高度上的信号覆盖；
(e)终止判断，具体包括：
当遗传代数大于设定值T或群体适应度平均值变化小于门限β时，则终止搜索，操作如下：
设全局变量Favg为群体适应度平均值，Fcavg为当前群体适应度平均值，计算
(e.1)判断当前迭代数t，当t＝1时，Favg＝Fcavg，执行第f步；否则执行步骤(e.2)；
(e.2)当t<T且Favg-Fcavg<＝β，或t＝T时，终止迭代，输出最优方案；否则执行第f步；
(f)存储最佳个体，具体包括：
(f.1)对群体的所有个体按适应度大小进行降序排序；
(f.2)如果当前迭代数t＝1时，令全局变量ColonyGroup表示迄今为止的适应度最高最好个体，大小为G·Q；否则通过适应度将ColonyGroup数组中最优个体与当前群体中的G·Q个适应度最高的个体进行比较，将适应度最高的G·Q个个体保存到ColonyGroup数组中，其中G为最优个体保存率；
(f.3)将迄今为止的ColonyGroup数组中最好的个体替换当前群体中G·Q个适应度最差的个体；
(f.4)重新按照适应度对群体个体排序；
(f.5)计算当前群体中每个个体适应度概率各个个体分配到的概率作为其能够被遗传到下一代的概率，基于这些概率值用比例选择的方法产生下一代群体；
(g)选择交叉变异；
执行经典遗传算法产生下一代新个体，包括：
(g.1)选择：从当前群体中按照一定的概率选出优良的个体，使它们有机会作为父代繁殖下一代子孙形成新群体，新群体个体个数为Q；
(g.2)交叉：将选择得到的新群体个体进行两两相互配对，配对的染色体依据交叉概率Pc按照某种方式交互其部分基因，从而形成两个新的个体，将新个体保存到新的群体中；
(g.3)变异：针对新群体个体通过一个变异率(pn)改变染色体位串上的一个或多个基因，将变异后的个体添入新群体中；
(g.4)分别将新群体中的Q个个体染色体进行解码，将解码后的地理坐标数组赋给Pts；
(g.5)至此一次站点优化布网分析结束，t＝t+1，进入第(c)步，执行多约束条件判定。
2.根据权利要求1所述的方法，其特征在于，对期望以所述监视台站覆盖服务的每个对象分别建立一个覆盖目标参数模型包括设置以下至少一个参数：
覆盖目标高度参数，整型数组，用于描述监视台站规划覆盖高度；
布站数量参数，整型，用于描述规划布网台站个数；
终端区覆盖重数参数，整型，用于描述规划设备台站多重覆盖度要求，最低满足一重覆盖；
覆盖终端区参数，浮点型二维数组，用于描述布网站点覆盖的区域，坐标设定格式为[R1,R2,..Rm],其中各元素为[x1,y1,x2,y2,...,xn,yn](n≥3)，m≥0,当m＝0时，表示当前组网规划与覆盖区域无关；
覆盖航路参数，浮点型二维数组，用于描述覆盖航路坐标设定，格式为[L1,L2,..Lm],其中各元素为[x1,y1,x2,y2,...,xn,yn](n≥2)，m≥0,当m＝0时，表示当前组网规划与覆盖航路无关；
航路覆盖重数参数，整型，用于描述规划台站对航路覆盖重数要求；
布网区域参数，浮点型数组，用于描述布网区域坐标格式[x1,y1,x2,y2,...,xn,yn](n≥3)；
航路覆盖宽度参数，整型，用于描述为航路左右延展的宽度；
航路合法布站宽度参数，整型，用于描述为航路附近合法的布站区域宽度；
实际覆盖搜索控制参数，布尔型，用于描述当前组网规划是否考虑现有监视台站覆盖；
实际覆盖搜寻区域参数，浮点型数组，用于描述当前组网规划考虑的现有设备必须在该搜索区域内，其坐标格式为[x1,y1,x2,y2,...,xn,yn](n≥3)。
3.根据权利要求1所述的方法，其特征在于，对所述覆盖目标参数模型进行计算以设定布站区域包括：
读取所述覆盖目标参数模型的配置；
判断所述覆盖目标参数模型中的布网区域参数是否为空；
如果所述布网区域参数为空，则判断所述覆盖目标参数模型中的覆盖航路参数是否为空，如果不为空，则以所述覆盖目标参数模型中的航路覆盖宽度、航路合法布站宽度为缓冲半径分别通过缓冲区分析得到航路覆盖面集合、航路合法布站区集合，以所述航路合法布站区集合进行合并，得到的合并区域作为新参数值赋给所述布站区域；如果所述覆盖航路参数为空，则将所述覆盖目标参数模型中的覆盖终端区参数的值赋给所述布站区域；
如果所述布网区域参数不为空，则判断所述覆盖目标参数模型中的覆盖航路参数是否为空，如果不为空，则以所述覆盖目标参数模型中的航路覆盖宽度、航路合法布站宽度为缓冲半径分别通过缓冲区分析得到航路覆盖面集合、航路合法布站区集合，将所述航路合法布站区集合的合并区域与所述覆盖目标参数模型中的布网区域参数取交集，将所述交集结果赋给所述布站区域；如果所述覆盖航路参数为空，则将所述覆盖目标参数模型中的布网区域参数的值赋给所述布站区域。
4.根据权利要求1所述的方法，其特征在于，对每种类型的监视台站分别建立一个所述候选站点参数模型，其包括设置以下至少一个参数：
设备类型参数，字符串型，用于描述所述监视台站的类型；
天线高度参数，浮点型，用于描述所述监视台站的规划设备天线的高度；
设备仰角参数，浮点型，用于描述所述监视台站的设备仰角；
塔台限制高度参数，浮点型，用于描述所述监视台站的塔台设置的高度限制；
地形作用距离参数，浮点型，用于描述地形对所述监视台站的信号覆盖影响范围的距离限制；
低空盲区高度参数，整型，用于描述所述监视台站的低空不能覆盖的高度限制；
高空盲区高度参数，整型，用于描述所述监视台站的高空不能覆盖的高度限制；
建筑物干扰源参数，字符串类型数组，用于描述所述监视台站的建筑物图层数据集名称数组；
建筑设施限制半径参数，浮点型，用于描述所述监视台站的不应有金属建筑物、密集的居民楼、高压输电线的限制半径；
有源干扰源参数，字符串类型数组，用于描述所述监视台站的有源信号干扰源；
有源干扰设施半径参数，浮点型，用于描述所述监视台站的指定半径范围内，不应有能产生有源干扰的电气设施。
5.根据权利要求1所述的方法，其特征在于，确定监视台站布站方案包括：
当所述站点优化布网分析满足终止条件，则根据需求输出指定个数的适应度值最高的最优个体方案，每个方案输出内容包括组网站点坐标集合，站点所在覆盖结果集。
6.一种监视台站布站装置，其特征在于，包括：
模型创建模块，用于对期望以所述监视台站覆盖服务的每个对象分别建立一个覆盖目标参数模型；
布站区域设定模块，用于对所述覆盖目标参数模型进行计算以设定布站区域；
优化布网分析模块，用于对所述布站区域进行站点优化布网分析，以确定监视台站布站方案；
所述优化布网分析模块采用遗传算法对所述布站区域进行站点优化布网分析；
所述遗传算法采用了遗传算法参数模型，所述遗传算法参数模型包括以下至少一个参数：
种群规模参数，整型，用于描述群体中所含个体的数量，一般取20～100；
个体编码长度参数，整型，用于描述个体编码长度，同时决定了编码精度；
编码精度参数，浮点型，用于描述求解台站位置坐标小数精度；
终止迭代数参数，整型，用于描述传算法终止进化代数；
适应度变化阈参数，浮点型，用于描述根据实际问题的求解过程得到的经验值；
交叉率参数，浮点型，用于描述进化过程中交叉操作的概率；
变异率参数，浮点型，用于描述进化过程中变异操作的概率；
最优个体保存率参数，浮点型，用于描述产生新一代群体时保留上一代最佳个体的百分率；
所述遗传算法采用了遗传算法参数模型，所述遗传算法参数模型包括以下至少一个参数：
种群规模参数，整型，用于描述群体中所含个体的数量，一般取20～100；
个体编码长度参数，整型，用于描述个体编码长度，同时决定了编码精度；
编码精度参数，浮点型，用于描述求解台站位置坐标小数精度；
终止迭代数参数，整型，用于描述传算法终止进化代数；
适应度变化阈参数，浮点型，用于描述根据实际问题的求解过程得到的经验值；
交叉率参数，浮点型，用于描述进化过程中交叉操作的概率；
变异率参数，浮点型，用于描述进化过程中变异操作的概率；
最优个体保存率参数，浮点型，用于描述产生新一代群体时保留上一代最佳个体的百分率；
采用遗传算法对所述布站区域进行站点优化布网分析包括：
(a)确定所述监视台站的信号覆盖区域，具体包括：
设置遮蔽角
其中ae为等效地球半径，hs为设备视线上最高遮蔽点的海拔高度，ls为遮蔽点到天线的斜距，ha为天线海拔高度；
设置在遮蔽角θs下，所述监视台站天线到高度hr的视在距离
设置在遮蔽角θs下，所述监视台站的实际作用距离
其中Rmax是所述监视台站的最大作用距离，G(θ)是所述监视台站天线的归一化方向函数；
设置所述监视台站不同高度层上的可视区，步骤如下：
(a.1)选择空域高度hr；
(a.2)计算不同方位扇区的遮蔽角，初始记方位ai＝0，其中0≦i<360；
(a.3)计算以所述监视台站为起点，地形作用距离为长度，方位ai上相交的等高线，其中该长度与等高线的交点为地形遮蔽点；
(a.4)根据上步求出的等高线的地理信息，找出高程值最大且距站点最近的等高线L，并获得L上的遮蔽点坐标；
(a.5)计算所述监视台站到上述遮蔽点的视在距离Rs(θs)；
(a.6)如果Rs(θs)≤Rr(θs)，则设置高度hr上可视区范围[θs，45°]；
(a.7)如果Rs(θs)＞Rr(θs)，令Rs(θn)＝Rr(θn)，求得θn，则设置高度hr上可视区范围[θn，45°]；
(a.8)以所述监视台站为中心，由球面坐标公式和直角坐标公式计算在ai方位角度上可视区角度θn上的作用距离在平面上相对设备位置的偏移量，从而得到ai方位上水平覆盖坐标点；
(a.9)i＝i+1，ai＝ai+1，循环执行(a.2)～(a.8)，直到ai＝360；
(a.10)分别将0～360°方位上获取的水平覆盖坐标点组成多边形区域，设置该区域为所述监视台站的信号覆盖区域；
(b)设置初始布站方案群体，具体包括：
布站方案群体中的每个个体表示为A：(F,Pts,Code,P)，F为个体适应度，Pts表示初始化个体中站点的坐标数组[(x1,y1),(x2,y2),...,(xn,yn)]，Code为对应经纬度坐标二进制编码的染色体串，P为个体适应度概率；
确定所述布站区域的最小外接矩形的经度[Xmin,Xmax]和纬度[Ymin,Ymax]；
设雷达网由N部雷达单元组成，以网内各雷达单元的地理位置坐标为染色体，各雷达单元的地理位置分别用Lx,Ly位二进制编码表示，N部雷达单元的染色体串表示为：
Code＝N(Lx+Ly)位编码形成染色体串对应着一种雷达网部署方案；
Lx或Ly位长度根据遗传算法公式计算得到，公式如下：
x＝xLxL-1xL-2...x2x1 公式2
其中公式1为二进制编码精度公式，公式2为一个个体的编码，公式3为解码公式，Umax，Umin分别是随机数值的最大值、最小值，L为二进制编码长度，δ为二进制编码的编码精度；
(c)进行多约束条件判定和群体修正，具体包括：
多约束条件判定条件分为环境约束条件和性能约束条件，环境约束条件包括实际布网区域；候选站点参数模型中的建筑设施限制半径，有源干扰设施半径；性能约束条件包括所述监视台站应满足的相邻单元之间最小部署距离dmin和相邻单元最大部署距离dmax，设置：
dmax＝min[dmax(Hi),dmax(Hh)]；
其中：PT＝发射机功率(dBm)，GT＝发射机天线增益(dB)，PL＝路径损耗，GR＝接收机天线旁瓣增益(dB)，f＝频率(MHz)，K＝32.45(距离以km为单位)；
(c.1)输入种群方案；
(c.2)对种群中每个个体的单元进行环境约束条件判断，对不满足的个体单元，在原有坐标的基础上加入随机项将该单元移出干扰区域，直到所有个体单元都满足环境约束条件；
(c.3)对种群中每个个体进行性能约束条件判断，如果群体中存在不满足性能约束条件的个体，执行步骤(c.4)-(c.6)；如果所有个体都满足性能约束条件，则执行步骤(c.7)；
(c.4)随机产生[0,2π]上均匀分布的随机数θ和[dmin,dmax]上均匀分布随机数k；
(c.5)对不满足性能约束条件的个体，以个体中某个单元坐标(x1,y1)为基准，相邻单元的位置调整为(x2,y2)＝(x1,y1)+k(cosθ,sinθ)；
(c.6)对修正后的个体单元进行环境约束条件判断，如果不满足，则返回执行步骤(c.3)；直到所有的个体都满足环境约束条件，则执行步骤(c.7)；
(c.7)输出修正后的种群；
(d)确定覆盖适应度和规划站点覆盖，具体包括：
假设系统由N部所述监视台站(D1D2...DN)组成，航路覆盖区为Sh，终端区Sz，M为高度层数，Sij为监视台站i在第j高度层的所述信号覆盖区域，ωj为加权系数，ωj取决于对各个高度层的关心程度，取值在[0，1]上，组网覆盖性能达到最优要求满足以下两个条件：
管制区覆盖范围最大化
设置覆盖系数ρ，取值范围为[0,1]，计算如下：
其中Sg为管制空域面积，包括终端区Sz和航路覆盖区Sh；
不同管制空域覆盖层数要求
设置冗余系数μ，取值[0，1]，μz为终端区覆盖冗余系数，μh为航路覆盖冗余系数，计算如下：
μ＝μz+μh 公式7
其中为覆盖冗余因子；
设置f＝k1ρ+k2μ 公式8
其中k1，k2，是权重系数，表示终端区覆盖与航路覆盖的重要程度；
具体覆盖适应度计算操作过程如下：
(d.1)令当前迭代次数为t，第一次初始化t＝1，t<T；
(d.2)根据候选站点设备类型，分别计算Q个群体方案中所有候选站点在指定的地理位置上的信号覆盖：
(d.3)判定实际覆盖搜索控制，如果该参数为TRUE，执行步骤(d.4)；否则执行步骤(d.5)；
(d.4)如果已知覆盖条件的现有台站个数为N1，规划台站数配置值为N2，则设置公式5中的设备个数为N＝N1+N2；
(d.5)如果没有现有台站参与规划，则设置N＝N2；
(d.6)判定覆盖目标模型参数配置，确定此次规划是正对终端区、航路或二者，如果正对终端区，则执行步骤(d.7)；如果正对航路，则执行步骤(d.8)；如果正对二者，则执行步骤(d.9)；
(d.7)结合公式4～8，令管制空域面积Sg＝Sz，μh＝0，μ＝μz；
(d.8)结合公式4～8，令管制空域面积Sg＝Sh，μz＝0，μ＝μh；
(d.9)结合公式4～8，令管制空域面积Sg＝Sz∪Sh，μ＝μz+μh；
(d.10)以i＝0，循环计算群体方案Q个个体的适应度Fi＝fi＝k1ρi+k2μi；其中当Nzc>＝N，冗余系数公式5和6中的冗余因子为0，Nzc为终端区域要求覆盖重数；
其中，方案中单个候选站点信号覆盖计算方法如下：
(1)利用DEM高程数据获取候选点Pt所在位置的高程值；
(2)计算候选点Pt的天线海拔高度＝高程值+塔台限制高度；
(3)根据当前配置的设备类型，利用该设备类型归一化方向函数计算该候选站点在指定高度上的信号覆盖；
(e)终止判断，具体包括：
当遗传代数大于设定值T或群体适应度平均值变化小于门限β时，则终止搜索，操作如下：
设全局变量Favg为群体适应度平均值，Fcavg为当前群体适应度平均值，计算
(e.1)判断当前迭代数t，当t＝1时，Favg＝Fcavg，执行第f步；否则执行步骤(e.2)；
(e.2)当t<T且Favg-Fcavg<＝β，或t＝T时，终止迭代，输出最优方案；否则执行第f步；
(f)存储最佳个体，具体包括：
(f.1)对群体的所有个体按适应度大小进行降序排序；
(f.2)如果当前迭代数t＝1时，令全局变量ColonyGroup表示迄今为止的适应度最高最好个体，大小为G·Q；否则通过适应度将ColonyGroup数组中最优个体与当前群体中的G·Q个适应度最高的个体进行比较，将适应度最高的G·Q个个体保存到ColonyGroup数组中，其中G为最优个体保存率；
(f.3)将迄今为止的ColonyGroup数组中最好的个体替换当前群体中G·Q个适应度最差的个体；
(f.4)重新按照适应度对群体个体排序；
(f.5)计算当前群体中每个个体适应度概率各个个体分配到的概率作为其能够被遗传到下一代的概率，基于这些概率值用比例选择的方法产生下一代群体；
(g)选择交叉变异；
执行经典遗传算法产生下一代新个体，包括：
(g.1)选择：从当前群体中按照一定的概率选出优良的个体，使它们有机会作为父代繁殖下一代子孙形成新群体，新群体个体个数为Q；
(g.2)交叉：将选择得到的新群体个体进行两两相互配对，配对的染色体依据交叉概率Pc按照某种方式交互其部分基因，从而形成两个新的个体，将新个体保存到新的群体中；
(g.3)变异：针对新群体个体通过一个变异率(pn)改变染色体位串上的一个或多个基因，将变异后的个体添入新群体中；
(g.4)分别将新群体中的Q个个体染色体进行解码，将解码后的地理坐标数组赋给Pts；
(g.5)至此一次站点优化布网分析结束，t＝t+1，进入第(c)步，执行多约束条件判定。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201310591118.1A CN103778477B (zh) | 2013-11-20 | 2013-11-20 | 监视台站布站方法和装置 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201310591118.1A CN103778477B (zh) | 2013-11-20 | 2013-11-20 | 监视台站布站方法和装置 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN103778477A CN103778477A (zh) | 2014-05-07 |
| CN103778477B true CN103778477B (zh) | 2017-07-18 |
Family
ID=50570685
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201310591118.1A Active CN103778477B (zh) | 2013-11-20 | 2013-11-20 | 监视台站布站方法和装置 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN103778477B (zh) |
Families Citing this family (12)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105046368A (zh) * | 2015-08-11 | 2015-11-11 | 国家电网公司 | 一种变电站主变容量和最优站址选取的方法 |
| CN105978642B (zh) * | 2016-05-18 | 2019-02-01 | 南京新兴空间电子系统有限公司 | 基于干扰大数据的无线监测站分析选址方法及系统 |
| CN106251080A (zh) * | 2016-08-05 | 2016-12-21 | 成都酉辰科技有限公司 | 边境电磁环境监测台情管理系统 |
| CN106371083B (zh) * | 2016-08-30 | 2018-11-27 | 长沙深之瞳信息科技有限公司 | 雷达监测点选取方法及装置以及回波信号处理方法及装置 |
| CN107830838A (zh) * | 2017-11-03 | 2018-03-23 | 中国银行股份有限公司 | 电磁态势地图绘制方法和装置 |
| CN108537361A (zh) * | 2018-03-06 | 2018-09-14 | 西安大衡天成信息科技有限公司 | 一种机动雷达和通信台站的选址评估方法 |
| CN109031309B (zh) * | 2018-06-20 | 2020-11-27 | 清华大学 | 一种组网气象雷达系统的拓扑结构优化方法及装置 |
| CN109581349B (zh) * | 2018-11-08 | 2020-06-26 | 中国电子科技集团公司第三十八研究所 | 一种雷达地形遮蔽分析及显示方法 |
| CN111091266A (zh) * | 2019-11-12 | 2020-05-01 | 青岛民航空管实业发展有限公司 | 空管台站的等级确定方法、装置及存储介质 |
| CN112652195B (zh) * | 2020-12-10 | 2022-01-14 | 北京航空航天大学杭州创新研究院 | 一种基于太阳能供电和移动通信传输的大规模ads-b飞行监视系统 |
| CN114950805B (zh) * | 2022-06-20 | 2023-11-10 | 成都飞机工业（集团）有限责任公司 | 一种飞机机翼机器人喷涂系统的站位优化设计方法 |
| CN115655282B (zh) * | 2022-12-14 | 2023-03-21 | 湖南高至科技有限公司 | 基于飞行器航迹和高程数据的遥测布站方法、装置及设备 |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2000333239A (ja) * | 1999-05-24 | 2000-11-30 | Mitsubishi Electric Corp | 屋内基地局の最適配置決定方法および最適配置決定システム |
| JP2009163511A (ja) * | 2008-01-07 | 2009-07-23 | Mazda Motor Corp | 部品搬送ダイヤの設定方法 |
| CN101650879A (zh) * | 2009-07-06 | 2010-02-17 | 民航数据通信有限责任公司 | 空管通信导航监视设备的布网优化的方法 |
| CN102522752A (zh) * | 2011-12-14 | 2012-06-27 | 天津大学 | 面向供电能力提高的配电网联络结构优化方法 |
-
2013
- 2013-11-20 CN CN201310591118.1A patent/CN103778477B/zh active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2000333239A (ja) * | 1999-05-24 | 2000-11-30 | Mitsubishi Electric Corp | 屋内基地局の最適配置決定方法および最適配置決定システム |
| JP2009163511A (ja) * | 2008-01-07 | 2009-07-23 | Mazda Motor Corp | 部品搬送ダイヤの設定方法 |
| CN101650879A (zh) * | 2009-07-06 | 2010-02-17 | 民航数据通信有限责任公司 | 空管通信导航监视设备的布网优化的方法 |
| CN102522752A (zh) * | 2011-12-14 | 2012-06-27 | 天津大学 | 面向供电能力提高的配电网联络结构优化方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN103778477A (zh) | 2014-05-07 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN103778477B (zh) | 监视台站布站方法和装置 | |
| CN101047937B (zh) | 基于遗传算法的移动通信频率规划方法 | |
| CN110418354A (zh) | 一种基于机器学习的免传播模型无线网络规划方法 | |
| CN108920481B (zh) | 一种基于手机定位数据的道路网重建方法及系统 | |
| CN102404752B (zh) | 一种gsm网络的频率规划方法 | |
| CN105430664B (zh) | 一种基于分类拟合预测传播路损的方法和装置 | |
| CN112902969B (zh) | 一种无人机在数据收集过程中的路径规划方法 | |
| CN112016812A (zh) | 多无人机任务调度方法、系统及存储介质 | |
| CN102170697A (zh) | 一种室内定位方法及装置 | |
| CN115630978B (zh) | 一种基于人口就业分布的多交通设施选址方法 | |
| CN111417066B (zh) | 定位方法、装置、电子设备和计算机可读存储介质 | |
| CN105530651A (zh) | 无线定位系统中最佳基站布设方法 | |
| CN114501530B (zh) | 基于深度强化学习的天线参数的确定方法和装置 | |
| CN110912627B (zh) | 一种基于数据驱动的小区接收功率预测方法 | |
| CN108833144A (zh) | 一种基于聚类统计模型的智能电表集中器选址方法 | |
| CN110366188B (zh) | 干扰测量点部署方法、干扰测量路径规划方法及系统 | |
| CN101650879B (zh) | 空管通信导航监视设备的布网优化的方法 | |
| CN114828026A (zh) | 基站规划方法、装置、设备、存储介质及程序产品 | |
| CN112835997B (zh) | 一种基于GeoSOT网格的电网巡检空间网格编码方法 | |
| CN115412926A (zh) | 基站站址规划方法、基站站址规划系统、设备及介质 | |
| CN112114386B (zh) | 高时空分辨率微波降雨监测方法 | |
| CN114513805A (zh) | 无线建模方法和系统 | |
| CN104700154A (zh) | 含不确定性的井下射频识别阅读器的高维多目标优化布局 | |
| CN112488343B (zh) | 一种基于超启发式算法的智慧城市智能感知终端选址方法 | |
| CN112434857B (zh) | 一种基于ga-pso优化的无人机机场网络选址布置方法 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| CB03 | Change of inventor or designer information |
Inventor after: Li Rui Inventor after: Zhu Liang Inventor after: Li Jian Inventor after: Zhou Zili Inventor after: Xue Kang Inventor after: Hou Changbo Inventor after: Fan Lijuan Inventor before: Li Rui Inventor before: Zhu Liang Inventor before: Li Jian |
|
| CB03 | Change of inventor or designer information | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |