一种面向航空器自主运行的空中交通态势复杂性评价方法与流程

1.本发明涉及空中交通管理技术领域，特别是涉及一种面向航空器自主运行的空中交通态势复杂性评价方法。


背景技术：

2.航空器自主保持间隔是未来空中交通管理的发展趋势。在航空器自主保持间隔运行模式下，航空器将部分或全部承担当前由管制员承担的间隔保障责任，也将具备选择最优航路、最优高度的可能性。然而航空器所具有的更为丰富的机动行为之间会产生冲突，这将极大影响空中交通态势复杂程度。现有的空中交通复杂性评价方法大多是面向当前管制模式，是以管制员工作负荷为基准，往往需要依靠管制员主观评价才能最终确定复杂性值。并没有考虑到航空器更为丰富的自主运行行为，并不适合未来航空器自主运行下的空中交通态势。


技术实现要素：

3.为了解决上述现有技术存在的问题，本发明的目的在于提供一种面向航空器自主运行的空中交通态势复杂性评价方法。该评价方法简单易用，可实现航空器自主保持间隔运行模式下空中交通态势复杂性的评价。本方法在接入航空器的航迹数据后，基于航空器机动行为建立轨迹空间，进一步生成空中交通态势网络模型，最后进行空中交通态势复杂性评价，无需人工干预、性能可靠，是一种通过数据对航空器自主运行模式下空中交通态势复杂性进行自动化评价的方法。可以弥补当前技术和方法的不足，为未来空中交通管理运行提供基础。
4.为了达到上述目的，本发明采取的技术方案是：一种面向航空器自主运行的空中交通态势复杂性评价方法，其特征在于，所述方法包括按顺序进行的下列步骤：
5.步骤1、生成航空器机动行为：接入航空器实时数据，对于每一航空器i，先读入其横坐标x0、纵坐标y0、高度h0、航向d0、速度s0；再依次生成5个机动行为：航向、高度保持不变的机动行为记为m1、左转航向15度的机动行为记为m2、右转航向15度的机动行为记为m3、爬升600米的机动行为记为m4、下降600米的机动行为记为m5。
6.步骤2、计算航空器轨迹空间：分别计算t分钟后航空器i在5类机动行为下的空间位置，并与当前位置连线生成每一航空器轨迹空间；当航空器i采取航向、高度保持不变的机动行为m1时，航空器i在t时刻的横坐标x
t,1
、纵坐标y
t,1
、高度h
t,1
、航向d
t,1
、速度s
t,1
的计算方法如式(1)所示：
[0007][0008]
将航空器i初始时刻空间位置点n0(x0,y0,h0)和采取机动行为m1后t时刻的空间位置点n1(x
t,1
,y
t,1
,h
t,1
)连线，生成航空器i采取机动行为m1的轨迹线段n0n1(i)；当航空器i采取左转航向15度的机动行为m2时，航空器i在t时刻的横坐标x
t,2
、纵坐标y
t,2
、高度h
t,2
、航向d
t,2
、速度s
t,2
的计算方法如式(2)所示：
[0009][0010]
将航空器i初始时刻空间位置点n0(x0,y0,h0)和采取机动行为m2后t时刻的空间位置点n1(x
t,2
,y
t,2
,h
t,2
)连线，生成航空器i采取机动行为m2的轨迹线段n0n2(i)；当航空器i采取右转航向15度的机动行为m3时，航空器i在t时刻的横坐标x
t,3
、纵坐标y
t,3
、高度h
t,3
、航向d
t,3
、速度s
t,3
的计算方法如式(3)所示：
[0011][0012]
将航空器i初始时刻空间位置点n0(x0,y0,h0)和采取机动行为m3后t时刻的空间位置点n1(x
t,3
,y
t,3
,h
t,3
)连线，生成航空器i采取机动行为m3的轨迹线段n0n3(i)；当航空器i采取爬升600米的机动行为m4时，航空器i在t时刻的横坐标x
t,4
、纵坐标y
t,4
、高度h
t,4
、航向d
t,4
、速度s
t,4
的计算方法如式(4)所示：
[0013][0014]
将航空器i初始时刻空间位置点n0(x0,y0,h0)和采取机动行为m4后t时刻的空间位置点n1(x
t,4
,y
t,4
,h
t,4
)连线，生成航空器i采取机动行为m4的轨迹线段n0n4(i)；当航空器i采取下降600米的机动行为m5时，航空器i在t时刻的横坐标x
t,5
、纵坐标y
t,5
、高度h
t,5
、航向d
t,5
、速度s
t,5
的计算方法如式(5)所示：
[0015][0016]
将航空器i初始时刻空间位置点n0(x0,y0,h0)和采取机动行为m5后t时刻的空间位置点n1(x
t,5
,y
t,5
,h
t,5
)连线，生成航空器i采取机动行为m5的轨迹线段n0n5(i)。
[0017]
步骤3、判断航空器轨迹空间中是否存在冲突：将所有航空器轨迹空间合并，并遍历轨迹空间集合，两两计算轨迹间的最小水平距离、最小垂直距离，如最小水平距离小于水平间隔标准10千米且最小垂直距离小于垂直间隔标准300米，则说明相应的两个机动行为存在冲突，标记为1，否则为零；设i和j分别是两个航空器，m1(i)、m2(i)、m3(i)、m4(i)、m5(i)是第i个航空器的5个机动行为，m1(j)、m2(j)、m3(j)、m4(j)、m5(j)是第j个航空器的5个机动行为,n0n1(i)、n0n2(i)、n0n3(i)、n0n4(i)、n0n5(i)是第i个航空器的轨迹空间，n0n1(j)、n0n2(j)、n0n3(j)、n0n4(j)、n0n5(j)是第j个航空器的轨迹空间；遍历两个航空器轨迹空间集合，两两计算轨迹空间集合中轨迹线段间的最小水平距离、最小垂直距离，如n0n1(i)和n0n2(j)的最小水平距离小于10千米且最小垂直距离小于300米，则航空器i的第1个机动行为m1(i)与航空器j的第2个机动行为m2(j)存在冲突，记为c(m1(i),m2(j))＝1，否则c(m1(i),m2(j))＝0。
[0018]
步骤4、基于轨迹空间的冲突关系生成空中交通态势网络：以空中交通态势中所有航空器的5种机动行为作为节点，将航空器机动行为间的冲突关系作为边，建立空中交通态势网络模型；此时空中交通态势网络中边的数量即为空中交通态势中所有航空器机动行为间冲突关系的总和，记为en。
[0019]
步骤5、评价空中交通态势复杂性：设空中交通态势中航空器数量为n，则n个航空器的机动行为数量为5n，这5n个节点中冲突关系最多有(5n
×
(5n-1)/2)个，再剔除同一个航空器不同机动行为轨迹间的冲突关系10n，即空中交通态势网络中冲突关系最多有(5n
×
(5n-1)/2-10n)个；最终，综合空中交通态势中的航空器数量、航空器机动行为轨迹间的冲突关系数量按照如下公式计算空中交通态势复杂性c：
[0020]
c＝n
×
(1 en/(5n
×
(5n-1)/2)-10n)
ꢀꢀꢀ
(6)。
[0021]
本发明所产生的有益效果是：采用面向航空器自主运行的空中交通态势复杂性评价方法，可以准确、快速地通过实际数据在充分考虑航空器机动行为的基础上，量化评价航空器自主保持间隔运行下空中交通态势的复杂性程度，本方法不受人为因素影响，评估方法简单易用。
附图说明
[0022]
图1为本发明提供的面向航空器自主运行的空中交通态势复杂性评价方法流程图；
[0023]
图2为本发明实施例1中空中交通态势网络示意图；
[0024]
图3为本发明实施例2中空中交通态势网络示意图。
具体实施方式
[0025]
以下结合附图及实施例对本发明做进一步说明。
[0026]
如图1所示，本发明提供的面向航空器自主运行的空中交通态势复杂性评价方法包括按顺序进行的下列步骤：
[0027]
步骤1、生成航空器机动行为：接入航空器实时数据，本实施例1的空中交通态势中有两个航空器，分别是航空器1、航空器2。航空器1的横坐标x0(1)＝600千米、纵坐标y0(1)＝300千米、高度h0(1)＝8800米、速度s0(1)＝900千米/小时、航向d0(1)＝90度，航空器1的初始空间位置点n0的坐标为(600,300,8800)。航空器2的横坐标x0(2)＝900千米、纵坐标y0(2)＝260千米、高度h0(2)＝8900米、速度s0(2)＝900千米/小时、航向d0(2)＝270度，航空器2的初始空间位置点n0的坐标为(900,260,8900)。
[0028]
本实施例2的空中交通态势中有两个航空器，分别是航空器3、航空器4。航空器3的横坐标x0(3)＝600千米、纵坐标y0(3)＝300千米、高度h0(3)＝8800米、速度s0(3)＝900千米/小时、航向d0(3)＝90度，航空器3的初始空间位置点n0的坐标为(600,300,8800)。航空器4的横坐标x0(4)＝900千米、纵坐标y0(4)＝260千米、高度h0(4)＝8900米、速度s0(4)＝900千米/小时、航向d0(4)＝270度，航空器4的初始空间位置点n0的坐标为(900,260,8900)。
[0029]
再对航空器1、航空器2、航空器3、航空器4依次生成5个机动行为，包括航向、高度保持不变的机动行为记为m1、左转航向15度的机动行为记为m2、右转航向15度的机动行为记为m3、爬升600米的机动行为记为m4、下降600米的机动行为记为m5。
[0030]
步骤2、计算航空器轨迹空间：设t为10分钟，分别计算10分钟后航空器i在5类机动行为下的空间位置，并与当前位置连线生成每一航空器轨迹空间。
[0031]
本实施例1中，当航空器1采取机动行为m1时，根据公式(1)计算航空器1在10分钟后的横坐标x
t,1
(1)＝750千米、纵坐标y
t,1
(1)＝300千米、高度h
t,1
(1)＝8800米、航向d
t,1
(1)＝90度、速度s
t,1
(1)＝900千米/小时。此时航空器1的空间位置点n1的坐标为(750,300,8800)，航空器1的轨迹n0n1(1)为从初始时刻空间位置点(600,300,8800)到空间位置点(750,300,8800)的一条直线段。
[0032]
当航空器1采取机动行为m2时，根据公式(2)计算航空器1在10分钟后的横坐标x
t,2
(1)＝744.9千米、纵坐标y
t,2
(1)＝338.8千米、高度h
t,2
(1)＝8800米、航向d
t,2
(1)＝75度、速度s
t,2
(1)＝900千米/小时。此时航空器1的空间位置点n2的坐标为(744.9,338.8,8800)，航空器1的轨迹n0n2(1)为从初始时刻空间位置点n0(600,300,8800)到空间位置点n2(744.9,338.8,8800)的一条直线段。
[0033]
当航空器1采取机动行为m3时，根据公式(3)计算航空器1在10分钟后的横坐标x
t,3
(1)＝744.9千米、纵坐标y
t,3
(1)＝261.2千米、高度h
t,3
(1)＝8800米、航向d
t,3
(1)＝105度、速度s
t,3
(1)＝900千米/小时。此时航空器1的空间位置点n3的坐标为(744.9,261.2,8800)，航空器1的轨迹n0n3(1)为从初始时刻空间位置点n0(600,300,8800)到空间位置点n3(744.9,261.2,8800)的一条直线段。
[0034]
当航空器1采取机动行为m4时，根据公式(4)计算航空器1在10分钟后的横坐标x
t,4
(1)＝750千米、纵坐标y
t,4
(1)＝300千米、高度h
t,4
(1)＝9400米、航向d
t,4
(1)＝90度、速度s
t,4
(1)＝900千米/小时。此时航空器1的空间位置点n4的坐标为(750,300,9400)，航空器1的轨迹n0n4(1)为从初始时刻空间位置点n0(600,300,8800)到空间位置点n4(750,300,
9400)的一条直线段。
[0035]
当航空器1采取机动行为m5时，根据公式(5)计算航空器1在10分钟后的横坐标x
t,5
(1)＝750千米、纵坐标y
t,5
(1)＝300千米、高度h
t,5
(1)＝8200米、航向d
t,5
(1)＝90度、速度s
t,5
(1)＝900千米/小时。此时航空器1的空间位置点n5的坐标为(750,300,8200)，航空器1的轨迹n0n5(1)为从初始时刻空间位置点n0(600,300,8800)到空间位置点n5(750,300,8200)的一条直线段。
[0036]
本实施例1中，当航空器2采取机动行为m1时，根据公式(1)计算航空器2在10分钟后的横坐标x
t,1
(2)＝750千米、纵坐标y
t,1
(2)＝260千米、高度h
t,1
(2)＝8900米、航向d
t,1
(2)＝270度、速度s
t,1
(2)＝900千米/小时。此时航空器2的空间位置点n1的坐标为(750,260,8900)，航空器2的轨迹n0n1(2)为从初始时刻空间位置点n0(600,300,8800)到空间位置点n1(750,260,8900)的一条直线段。
[0037]
当航空器2采取机动行为m2时，根据公式(2)计算航空器2在10分钟后的横坐标x
t,2
(2)＝755.1千米、纵坐标y
t,2
(2)＝221.2千米、高度h
t,2
(2)＝8900米、航向d
t,2
(2)＝255度、速度s
t,2
(2)＝900千米/小时。此时航空器2的空间位置点n2的坐标为(755.1,221.2,8900)，航空器2的轨迹n0n2(2)为从初始时刻空间位置点n0(600,300,8800)到空间位置点n2(755.1,221.2,8900)的一条直线段。
[0038]
当航空器2采取机动行为m3时，根据公式(3)计算航空器2在10分钟后的横坐标x
t,3
(2)＝755.1千米、纵坐标y
t,3
(2)＝298.8千米、高度h
t,3
(2)＝8900米、航向d
t,3
(2)＝285度、速度s
t,3
(2)＝900千米/小时。此时航空器2的空间位置点n3的坐标为(755.1,298.8,8900)，航空器2的轨迹n0n3(2)为从初始时刻空间位置点n0(600,300,8800)到空间位置点n3(755.1,298.8,8900)的一条直线段。
[0039]
当航空器2采取机动行为m4时，根据公式(4)计算航空器2在10分钟后的横坐标x
t,4
(2)＝750千米、纵坐标y
t,4
(2)＝260千米、高度h
t,4
(2)＝9500米、航向d
t,4
(2)＝270度、速度s
t,4
(2)＝900千米/小时。此时航空器2的空间位置点n4的坐标为(750,260,9500)，航空器2的轨迹n0n4(2)为从初始时刻空间位置点n0(600,300,8800)到空间位置点n4(750,260,9500)的一条直线段。
[0040]
当航空器2采取机动行为m5时，根据公式(5)计算航空器2在10分钟后的横坐标x
t,5
(2)＝750千米、纵坐标y
t,5
(2)＝260千米、高度h
t,5
(2)＝8300米、航向d
t,5
(2)＝270度、速度s
t,5
(2)＝900千米/小时。此时航空器2的空间位置点n5的坐标为(750,260,8300)，航空器2的轨迹n0n5(2)为从初始时刻空间位置点n0(600,300,8800)到空间位置点n5(750,260,8300)的一条直线段。
[0041]
本实施例2中，当航空器3采取机动行为m1时，根据公式(1)计算航空器3在10分钟后的横坐标x
t,1
(3)＝750千米、纵坐标y
t,1
(3)＝300千米、高度h
t,1
(3)＝8800米、航向d
t,1
(3)＝90度、速度s
t,1
(3)＝900千米/小时。此时航空器3的空间位置点n1的坐标为(750,260,8800)，航空器3的轨迹n0n1(3)为从初始时刻空间位置点n0(600,300,8800)到空间位置点n1(750,260,8300)的一条直线段。
[0042]
当航空器3采取机动行为m2时，根据公式(2)计算航空器3在10分钟后的横坐标x
t,2
(3)＝744.9千米、纵坐标y
t,2
(3)＝338.8千米、高度h
t,2
(3)＝8800米、航向d
t,2
(3)＝75度、速度s
t,2
(3)＝900千米/小时。此时航空器3的空间位置点n2的坐标为(744.9,338.8,8800)，航
空器3的轨迹n0n2(3)为从初始时刻空间位置点n0(600,300,8800)到空间位置点n2(744.9,338.8,8800)的一条直线段。
[0043]
当航空器3采取机动行为m3时，根据公式(3)计算航空器3在10分钟后的横坐标x
t,3
(3)＝744.9千米、纵坐标y
t,3
(3)＝261.2千米、高度h
t,3
(3)＝8800米、航向d
t,3
(3)＝105度、速度s
t,3
(3)＝900千米/小时。此时航空器3的空间位置点n3的坐标为(744.9,261.2,8800)，航空器3的轨迹n0n3(3)为从初始时刻空间位置点n0(600,300,8800)到空间位置点n3(744.9,261.2,8800)的一条直线段。
[0044]
当航空器3采取机动行为m4时，根据公式(4)计算航空器3在10分钟后的横坐标x
t,4
(3)＝750千米、纵坐标y
t,4
(3)＝300千米、高度h
t,4
(3)＝9400米、航向d
t,4
(3)＝90度、速度s
t,4
(3)＝900千米/小时。此时航空器3的空间位置点n4的坐标为(750,300,9400)，航空器3的轨迹n0n4(3)为从初始时刻空间位置点n0(600,300,8800)到空间位置点n4(750,300,9400)的一条直线段。
[0045]
当航空器3采取机动行为m5时，根据公式(5)计算航空器3在10分钟后的横坐标x
t,5
(3)＝750千米、纵坐标y
t,5
(3)＝300千米、高度h
t,5
(3)＝8200米、航向d
t,5
(3)＝90度、速度s
t,5
(3)＝900千米/小时。此时航空器3的空间位置点n5的坐标为(750,300,8200)，航空器3的轨迹n0n5(3)为从初始时刻空间位置点n0(600,300,8800)到空间位置点n5(750,300,8200)的一条直线段。
[0046]
本实施例2中，当航空器4采取机动行为m1时，根据公式(1)计算航空器4在10分钟后的横坐标x
t,1
(4)＝740千米、纵坐标y
t,1
(4)＝300千米、高度h
t,1
(4)＝8900米、航向d
t,1
(4)＝270度、速度s
t,1
(4)＝900千米/小时。此时航空器4的空间位置点n1的坐标为(740,300,8900)，航空器4的轨迹n0n1(4)为从初始时刻空间位置点n0(900,260,8900)到空间位置点n1(740,300,8900)的一条直线段。
[0047]
当航空器4采取机动行为m2时，根据公式(2)计算航空器4在10分钟后的横坐标x
t,2
(4)＝745.1千米、纵坐标y
t,2
(4)＝261.2千米、高度h
t,2
(4)＝8900米、航向d
t,2
(4)＝255度、速度s
t,2
(4)＝900千米/小时。此时航空器4的空间位置点n2的坐标为(745.1,261.2,8900)，航空器4的轨迹n0n2(4)为从初始时刻空间位置点n0(900,260,8900)到空间位置点n2(745.1,261.2,8900)的一条直线段。
[0048]
当航空器4采取机动行为m3时，根据公式(3)计算航空器4在10分钟后的横坐标x
t,3
(4)＝745.1千米、纵坐标y
t,3
(4)＝338.8千米、高度h
t,3
(4)＝8900米、航向d
t,3
(4)＝285度、速度s
t,3
(4)＝900千米/小时。此时航空器4的空间位置点n3的坐标为(745.1,338.8,8900)，航空器4的轨迹n0n3(4)为从初始时刻空间位置点n0(900,260,8900)到空间位置点n3(745.1,338.8,8900)的一条直线段。
[0049]
当航空器4采取机动行为m4时，根据公式(4)计算航空器4在10分钟后的横坐标x
t,4
(4)＝740千米、纵坐标y
t,4
(4)＝300千米、高度h
t,4
(4)＝9500米、航向d
t,4
(4)＝270度、速度s
t,4
(4)＝900千米/小时。此时航空器4的空间位置点n4的坐标为(740,300,9500)，航空器4的轨迹n0n4(4)为从初始时刻空间位置点n0(900,260,8900)到空间位置点n4(740,300,9500)的一条直线段。
[0050]
当航空器4采取机动行为m5时，根据公式(5)计算航空器4在10分钟后的横坐标x
t,5
(4)＝740千米、纵坐标y
t,5
(4)＝300千米、高度h
t,5
(4)＝8300米、航向d
t,5
(4)＝270度、速度st,5
(4)＝900千米/小时。此时航空器4的空间位置点n5的坐标为(740,300,8300)，航空器4的轨迹n0n5(4)为从初始时刻空间位置点n0(900,260,8900)到空间位置点n5(740,300,8300)的一条直线段。
[0051]
步骤3、判断航空器轨迹空间中是否存在冲突：将实施例1和实施例2中所有航空器轨迹空间合并，并遍历轨迹空间集合，两两判断机动行为对应的轨迹是否有冲突。
[0052]
实施例1的轨迹空间中有10条轨迹，其中航空器1机动行为的5条轨迹为n0n1(1)、n0n2(1)、n0n3(1)、n0n4(1)、n0n5(1)，航空器2机动行为的5条轨迹为n0n1(2)、n0n2(2)、n0n3(2)、n0n4(2)、n0n5(2)。如表1所示，在这10条轨迹中，轨迹n0n1(1)与轨迹n0n3(2)的最小水平距离5.2千米小于等于水平间隔标准10千米、最小垂直距离100米小于等于垂直间隔标准300米，即c(m1(1),m3(2))＝1；轨迹n0n1(2)与轨迹n0n3(1)的最小水平距离5.2千米小于等于水平间隔标准10千米、最小垂直距离100米小于等于垂直间隔标准300米，即c(m3(1),m1(2))＝1。
[0053]
表1
[0054][0055]
实施例2的轨迹空间中有10条轨迹，其中航空器3机动行为的5条轨迹为n0n1(3)、n0n2(3)、n0n3(3)、n0n4(3)、n0n5(3)，航空器4机动行为的5条轨迹为n0n1(4)、n0n2(4)、n0n3(4)、n0n4(4)、n0n5(4)。如表2所示，在这10条轨迹中，轨迹n0n1(3)与轨迹n0n1(4)的最小水平距离10千米小于等于水平间隔标准10千米、最小垂直距离100米小于等于垂直间隔标准300米，即c(m1(3),m1(4))＝1；轨迹n0n2(3)与轨迹n0n3(4)的最小水平距离0.2千米小于等于水平间隔标准10千米、最小垂直距离100米小于等于垂直间隔标准300米，即c(m2(3),m3(4))＝1；轨迹n0n2(4)与轨迹n0n3(3)的最小水平距离0.2千米小于等于水平间隔标准10千米、最小垂直距离100米小于等于垂直间隔标准300米，即c(m2(4),m3(4))＝1；轨迹n0n4(3)与轨迹n0n4(4)的最小水平距离10千米小于等于水平间隔标准10千米、最小垂直距离100米小于等于垂直间隔标准300米，即c(m4(3),m4(4))＝1；轨迹n0n5(3)与轨迹n0n5(4)的最小水平距离10千米小于等于水平间隔标准10千米、最小垂直距离100米小于等于垂直间隔标准300米，即c(m5(3),m5(4))＝1。
[0056]
表2
[0057][0058]
步骤4、基于轨迹空间的冲突关系生成网络：在实施例1中，航空器数量n＝2，以航空器1对应的机动行为m1(1)、m2(1)、m3(1)、m4(1)、m5(1)以及航空器2对应的机动行为m1(2)、m2(2)、m3(2)、m4(2)、m5(2)作为节点，即节点数量为5n。以航空器机动行为对应轨迹间的冲
突关系作为边。步骤3计算得知，c(m1(1),m3(2))＝1、c(m3(1),m1(2))＝1，故实施例1对应的网络模型中，只有m1(1)和m3(2)相连、m3(1)和m1(2)相连，如图2所示。
[0059]
在实施例2中，以航空器3对应的机动行为m1(3)、m2(3)、m3(3)、m4(3)、m5(3)以及航空器4对应的机动行为m1(4)、m2(4)、m3(4)、m4(4)、m5(4)作为节点，即节点数量为5n。以航空器机动行为对应轨迹间的冲突关系作为边。步骤3计算得知，c(m1(3),m1(4))＝1、c(m2(3),m3(4))＝1、c(m2(4),m3(3))＝1、c(m4(3),m4(4))＝1、c(m5(3),m5(4))＝1，故实施例2对应的网络模型中，只有m1(3)和m1(4)相连、m2(3)和m3(4)相连、m2(4)和m3(4)相连、m4(3)和m4(4)相连、m5(3)和m5(4)相连，如图3所示。
[0060]
步骤5、评价空中交通态势复杂性：在实施例1的空中交通态势网络中，节点之间只有两条连边，故航空器机动行为间的冲突关系之和en＝2，根据公式(6)可计算出当前空中交通态势复杂性c＝2.16。在实施例2的空中交通态势网络中，节点之间只有5条连边，故航空器机动行为间的冲突关系之和en＝5，根据公式(6)可计算出当前空中交通态势复杂性c＝2.40。
[0061]
对比实施例1和实施例2，虽然两个实施例中的航班数相同，但从航空器机动行为之间的冲突关系看，实施例1中航空器机动行为间的冲突关系较少，所反映的空中交通态势复杂性则相对实施例2较小。
[0062]
本方法相对现有只考虑航空器数量的评价方法而言，充分考虑了航空器的机动行为，并在描述机动行为轨迹冲突关系的基础上，建立了空中交通态势网络模型，形成了更为精细化的空中交通态势复杂性评价方法，可以更全面反映航空器自主运行下空中交通态势复杂程度。