1.一种基于空域网格的赋值方法,其特征在于,具体步骤如下:首先、针对目标空域A进行网格递归划分,形成m行n列均匀分布的子网格,视为m×n阶矩阵;然后、在t=t a 时刻针对目标空域A的单个航空器i,根据实时的飞行动态、气象信息和飞行计划,对该航空器在每个子网格的位置分别赋予参数,构成确定型网格赋值矩阵为 确定型网格赋值矩阵为 表示为:u=1,2,...,m,v=1,2,...,n; 表示在t=t a 时刻航空器i所在的序号为u×v的子网格被赋予的参数值;根据实时的飞行动态、气象信息和飞行计划,将航空器i实际所在的序号为u×v的子网格赋值为 其余位置均为0,由此得到确定的航空器位置信息的网格赋值矩阵;同理,得到该目标空域A内每个航空器在t=t a 时刻各自对应的确定网格赋值矩阵;接着、根据各个航空器在飞行过程中存在的误差,建立概率型网格赋值矩阵,并对各航空器所在的实际子网格的概率进行赋值;所述的概率型网格赋值矩阵的建立过程为:首先,航空器i经过子网格x u×v 的概率为P(X=x u×v )=p u×v ,则将子网格x u×v 在网格赋值矩阵 的对应位置赋值为p u×v ;将该子网格的编码序号x u×v 转化为矩阵行列号;矩阵行号为 矩阵列号为n u×v ≡[u×v-1](modn)+1;其中 为上取整函数,mod为取余函数;同理,对目标空域A内在t=t a 时刻每个航空器各自所在的实际子网格赋值概率值;得到各自对应的矩阵行列号;最后、利用各航空器的确定型网格赋值矩阵和概率型网格赋值矩阵分别计算目标空域A的复杂度、网格内飞行流量和碰撞概率,为空管提供参考数据;所述目标空域A的各子网格的交通复杂度,计算过程如下:首先,针对子网格x l ,计算交通复杂度的各个指标:①交通密度N(x l ):N(x l )=W k (x l )/(f s (x l )×H);其中W k (x l )表示网格x l 在第k个时间段内的交通流量;H表示该飞行空域管理高度;f s (x l )表示空域网格x l 的球面投影面积;②飞行状态cond(x l ):统计每个网格单元在第k个时间段内飞行器正在进行爬升和下降的次数和时间;③飞机机型混合影响v(x l ):即单元网格x l 内的不同航空器速度不同引起的复杂度;④天气影响交通指数:WITI(x l )=N k (x l )×W k (x l )其中N k (x l )表示第k个时间段内,交通流经过网格x l 的国家对流天气预报报告次数;W k (x l )表示网格x l 在第k个时间段内的飞行流量;然后,利用各个指标计算子网格x l 的交通复杂度,公式为:TC(x l )=N(x l )×θ(x l )其中θ(x l )为空域网格x l 的复杂度因子,由复杂度指标因子复合得到,其计算公式为:θ(x l )=cond(x l )×v(x l )×WITI(x l );所述的飞行流量是指:在时段T=[t 1 ,t N ]内进入空域A的航空器总数;计算过程为:首先,针对t=t a 时刻,将空域A内所有航空器的确定型网格赋值矩阵内每个位置的数值累计求和,得到 并将 转换为二进制形式; 然后,计算二进制字符串 的Hamming重量,即空域A内t a 时刻的飞行流量 Hamming重量是二进制字符串中1的个数;同理,分别计算各个时刻下的二进制字符串 得到各时刻对应的飞行流量;遍历各个时刻,将当前时刻t j 的飞行流量记为 根据当前时刻与相邻时刻的关系 计算相邻两时刻的Hamming距离 并更新下一时刻的飞行流量值: 直至j=N;j初始化取1;Hamming距离 是指对二进制字符串 和 的元素进行异或运算并统计结果为1的个数值;所述的空域A内的碰撞概率,计算过程如下:首先,针对目标空域A,利用网格层级为r 0 的空域网格单元,对其进行网格递归划分;网格层级r 0 的选择根据航空器安全包络的大小选取;然后,计算M架飞机Y={y 0 ,y 1 ,...,y j ,...,y M-1 }飞往不同空域网格单元中的概率;编号为y j 的飞机在编码为x i 的空域网格单元的概率p ij ,表示如下:p ij =P(X=x i ,Y=y i )进一步,根据各飞机在空域网格单元中的概率分布,计算每个子网格的模糊碰撞概率;编码为x i 的空域网格单元,模糊碰撞概率P i '公式为:所有其余的每个空域网格单元的模糊碰撞概率最终组成该单元的模糊碰撞概率集合P';设计空域网格单元x i 的三元组数据G i (x i ,P i ',S i );S i 为空域网格单元x i 中满足概率p ij ≠0的所有飞机组成的集合;映射关系如下: 且 使得f 1 -1 (P i ')=x i ; 且 使得f 2 -1 (S i )=x i ;利用设定的碰撞概率算子P A 对有碰撞风险的空域网格单元进行筛选,依次选择每个单元的三元组数据,当P i '≥P A 则在该单元存在碰撞风险,返回其对应的三元组数据保存至数组集合G*,否则,碰撞风险忽略不计,继续判段下一个空域网格单元;直至所有单元判断完毕得到新的集合G*;集合G*即为所有存在碰撞风险的空域A内的网格单元、碰撞概率以及其对应的航空器编号。
2.如权利要求1所述的一种基于空域网格的赋值方法,其特征在于,所述的对目标空域A进行网格递归划分的过程为:从网格坐标的原点开始,按矩阵行号从小到大、同行号的按列号从小到大,建立“Z”字排序法,依次对子网格进行数值编码;得到网格的编码序号l=1,2,...,m×n;用离散随机变量X={x 1 ,x 2 ,...,x l ,...,x m×n }表示。