基于文法的舰船目标运动建模与仿真方法与流程
本发明涉及一种基于文法的舰船目标运动建模与仿真方法，特别涉及基于文法的舰船目标的运动建模与仿真方法，属于目标识别与仿真技术领域。
背景技术：
海上集成监视系统主要包括以各类雷达(岸基、船载、机载)组成的传感系统、通信系统和数据融合处理中心。主要是对近海区域、江河、湖泊、港口以及海峡等航道的船舶进行连续跟踪处理，其最终目的是对海上舰船的真实轨迹进行跟踪，并显示处于被监视海洋区域内各种类型不同的船舶目标。国家海上监视系统的成熟完善与否决定着其保卫领海安全、维护国家海洋权利、发展区域海洋经济的能力，也是提升海军作战能力、促进国家海军向科技化、信息化升级的关键环节。因此对海上舰船目标跟踪监视成为海上集成监视系统的首要任务。
已有的海上舰船目标的运动建模方法难以产生多样化的运动轨迹，难以反映舰船航行时的机动性，且产生的数据维度高。大多采用蒙特卡罗仿真方法，采用单一的运动模型，包括直线运动、马尔科夫运动模型。已有方法难以将不同的运动模型灵活的融合，生成的舰船轨迹难以反映真实舰船的航行情况，并且反映不出舰船目标运动的机动性。海上监视区域内舰船在一段时间内运动变化幅度大，且受到的干扰因素较多，这些复杂的态势运动构成了一个情景。已有的舰船运动模型产生的情景存在难以真实描述舰船运动的复杂性和多样性、数据量大、维度高等问题。
因此，对海上目标舰船运动建模与仿真，来产生逼真的实验情景、大量的测试数据对海上集成监视系统中舰船目标的跟踪监视非常重要，这对诸如目标跟踪算法的选型和改进、使用风险的合理规避并最终提高国家海上安全都意义重大。
技术实现要素：
针对现有技术存在的上述缺陷，本发明提出了一种基于文法的舰船目标运动建模与仿真方法，能够自动产生形式简洁、包含多艘舰船目标且运动态势多变的海上舰船目标运动的情景；且生成的舰船多为高机动目标，舰船航行运动状态改变丰富，能够模拟出一艘乃至数百艘移动目标在海域的不同运动状态下的情景，可以很好的体现出多艘舰船航行时，彼此之间产生的运动关系。
本发明所述的基于文法的舰船目标运动建模与仿真方法，其特征在于，包括如下步骤：
s1：构建舰船目标运动模型：构建参数化上下文无关文法的情景模型，用有限个数的文法规则，来描述舰船的运动轨迹；
s2：构建文法模型：将舰船航行中的不确定因素，包括初始状态、实体、事件、关系事件模型，用上下文无关文法表示；
s3：推演文法模型：将文法模型进行推演，首先经过推演转换规则，产生舰船情景的推演树，再将推演树中包含的情景信息进一步展开和丰富，产生派生树作为情景框架；
s4：读取和赋值：对情景框架进行读取和赋值，将情景框架中每一艘舰船目标的运动态势字串读取出来，获得情景字串，对其中各个用终止符表示的参数进行赋值。
优选地，所述s1中，文法规则包括：
规则1：整个模型的建立从此语法开始，公式如下：
s-＞t(0)i0(x)a(x)
其中，终止符t(0)为实体x的初始时间，即舰船x在海域内出现的时间；非终止符a(x)表示实体x还可以继续变化；引入了第一艘船x并且初始化此船的各项初始参数、包括初始时间t(0)和初始速度、航向、经纬度，用i0(x)统一表示；
规则2：
a(x)-＞tδ(x)vδ(x)a(x)
其中，终止符tδ(x)和vδ(x)分别表示实体x发生改变的时刻和改变量，具体情况为：船x在tδ(x)时刻发生了速度改变，大小为vδ(x)；a(x)表示实体x还能继续改变；
规则3：
a(x)-＞tδ(x)dδ(x)a(x)
其中，终止符tδ(x)和dδ(x)分别表示实体x发生改变的时刻和改变量，具体情况为：船x在tδ(x)时刻发生了转向，转向角的大小为dδ(x)；a(x)表示实体x还能继续改变；
规则4：
a(x)-＞tδ(x)aδ(x)tk(x)a(x)
其中，经过tδ(x)时间后加速度改变(默认所有模式的a＝0)，加速度改变的大小为aδ(x)即进入到加速度模型(ca)，且保持tk(x)的时间；
规则5：
a(x)-＞tδ(x)aδ(x)tk(x)b(x)
其中，经过tδ(x)时间后进入加速度模型，与规则4不同的是，维持时间结束后，舰船仍以当前的加速度继续维持上一时刻的状态；
规则6：
a(x)-＞tδ(x)wδ(x)tk(x)a(x)
其中，经过tδ(x)时间后引入角速度w，即转向角的改变为(θ＝w*t+θδ)，进入恒定转率和速度模型，并且还要维持tk(x)的时间；
规则7：
a(x1)-＞a(x1)+t(x1)[i0(x2)a(x2)+i0(x3)a(x3)
+i0(x4)a(x4)+i0(x5)a(x5)+i0(x6)a(x6)]
其中，船x1引入五艘新的舰船，分别用a(x2)、a(x3)、a(x4)、a(x5)、a(x6)表示这些实体还可继续变化，同时初始化这些新实体的初始属性分别用i0(x2)、i0(x3)、i0(x4)、i0(x5)、i0(x6)表示，而a(x1)表示实体x1还可以继续变化；
规则8：
a(x1)-＞a(x1)+[i0(x2)a(x2)+i0(x3)a(x3)]
其中，船x1引入两艘新的舰船，分别用a(x2)、a(x3)表示这些实体还可继续变化，同时初始化这些新实体的初始属性分别用i0(x2)、i0(x3)表示，而a(x1)表示实体x1还可以继续变化；
规则9：
a(x1)-＞a(x1)+i0(x2)a(x2)
其中，船x1引入1艘新的舰船实体，用a(x2)表示这个新实体还可继续变化，同时初始化这个实体的初始属性用i0(x2)来表示，新舰船a(x2)的初始时间为a(x1)的初始时间，而a(x1)表示实体x1也还可以继续变化；
规则10：
a(x)-＞ε
其中，此规则意为实体x结束变化，非终止符变为终止符；
规则11：
b(x)-＞*a(x)
其中，此规则意为维持时间结束后,找到参照的运动模式；当上一规则结束后,按照上一时刻的运动模型继续运动,直到下一次规则选取的到来；若没有选择规则11,则按照上一时刻的角度和速度进行匀速直线运动。
优选地，所述s1中，文法规则包括：舰船从雷达探测边缘区域进入，在节点引入新舰船时有不同：
规则7：
a(x1)-＞a(x1)+t(x1)[i0b(x2)a(x2)+i0b(x3)a(x3)
+i0b(x4)a(x4)+i0b(x5)a(x5)+i0b(x6)a(x6)]
其中，终止符t(x1)表示过了t(x1)段时间后，船x1从雷达探测的边界点引入五艘新的舰船实体，分别用a(x2)、a(x3)、a(x4)、a(x5)、a(x6)表示这些实体还可继续变化，i0b可以从i0bn、i0bs、i0be、i0bw的初始属性中选取，而a(x1)表示实体x1还可以继续变化；
规则8：
a(x1)-＞a(x1)+t(x1)[i0b(x2)a(x2)+i0b(x3)a(x3)]
其中，终止符t(x1)表示过了t(x1)段时间后，船x1从雷达探测的边界点引入两艘新的舰船实体，分别用a(x2)、a(x3)表示这些实体还可继续变化，i0b可以从i0bn、i0bs、i0be、i0bw的初始属性中选取，而a(x1)表示实体x1还可以继续变化；
规则9：
a(x1)-＞a(x1)+t(x1)i0b(x2)a(x2)
其中，终止符t(x1)表示过了t(x1)段时间后，船x1从雷达探测的边界点引入1艘新的舰船实体，用a(x2)表示这个新实体还可继续变化，同时初始化这个实体的初始属性用i0b(x2)来表示，新舰船a(x2)的初始时间为a(x1)的初始时间加上t(x1)，而a(x1)表示实体x1也还继续变化。
优选地，所述s2中，舰船航行中的不确定因素，包括：
s21：初始状态：初始状态为舰船目标的初始运动状态信息、雷达系统的目标检测误差；
s22：实体：实体包含不同的舰船目标；
s23：事件：事件用实体状态的改变来描述，包括增加新船、改变速度、改变角度、改变弧度；
s24：关系：关系可分为实体与初始状态间的依存关系、实体与实体间的时间关系、实体与事件间的依存关系、事件间的关系、实体与参数间的归属关系、以及实体与事件间的归属关系。
优选地，所述s3中，推演转换规则包括如下小步：
s31：推演从起始节点1开始(200)，且除了起始时不能再调用规则1，推演过程是自顶向下，按层推演的；
s32：每次推演从规则2-9随机选择非终止符来调用(201)；
s33：根据新调用的规则增删非终止符的数量(202)，每当调用规则7、8、9时，会引入新的非终止符，此节点可拥有的子节点数量将会增加；
s34：将规则号作为子节点，继续推演(203)。
优选地，所述s4中，对情景框架进行读取和赋值，将情景框架中每一艘舰船目标的运动态势字串分别读取出来，获得情景字串，对其中各个用终止符表示的参数进行赋值，即完成了一个测试情景的产生，多次推演情景模型即可产生大量各种各样的情景。
本发明的有益效果是：(1)能够自动产生形式简洁、包含多艘舰船目标且运动态势多变的海上舰船目标运动的情景；(2)生成的舰船多为高机动目标，舰船航行运动状态改变丰富，包括运动方向的改变、速度改变、加速度大小改变、角速度改变、从边界点进入等；能够模拟出一艘乃至数百艘移动目标在海域的不同运动状态下的情景，可以很好的体现出多艘舰船航行时，彼此之间产生的运动关系；(3)产生的情景数据量小、维度低，方便数据处理，可以作为后续工作的处理；(4)舰船运动模型丰富，且产生的舰船转角范围、速度改变区间、舰船出现的经纬度等范围可人为调节；(5)事件性情景模型可以用一连串的事件字串来表达多目标复杂的运动轨迹，且产生的数据更加灵活多变，便于重新组合。
附图说明
图1是本发明的情景生成总示意图。
图2(a)是本发明的推演树示意图。
图2(b)是本发明的情景推演过程流程图。
图3是本发明的具体推演过程。
图4是本发明的推演树产生派生树示意图。
图5是本发明的遍历派生树示意图。
图6(a)是本发明无噪声的实验结果图。
图6(b)是本发明有噪声的实验结果图。
图6(c)是本发明有噪声和海杂波的实验结果图。
具体实施方式
为了使本发明目的、技术方案更加清楚明白，下面结合实施例，对本发明作进一步详细说明。
实施例1：
如图1至图6(c)所示，本发明所述的基于文法的舰船目标运动建模与仿真方法，其特征在于，包括如下步骤：
s1：构建舰船目标运动模型：构建参数化上下文无关文法的情景模型，用有限个数的文法规则，来描述舰船的运动轨迹；
s2：构建文法模型：将舰船航行中的不确定因素，包括初始状态、实体、事件、关系事件模型，用上下文无关文法表示；
s3：推演文法模型：将文法模型进行推演，首先经过推演转换规则，产生舰船情景的推演树，再将推演树中包含的情景信息进一步展开和丰富，产生派生树作为情景框架；
s4：读取和赋值：对情景框架进行读取和赋值，将情景框架中每一艘舰船目标的运动态势字串读取出来，获得情景字串，对其中各个用终止符表示的参数进行赋值。
优选地，所述s1中，文法规则包括：
规则1：整个模型的建立从此语法开始，公式如下：
s-＞t(0)i0(x)a(x)
其中，终止符t(0)为实体x的初始时间，即舰船x在海域内出现的时间；非终止符a(x)表示实体x还可以继续变化；引入了第一艘船x并且初始化此船的各项初始参数、包括初始时间t(0)和初始速度、航向、经纬度，用i0(x)统一表示；
规则2：
a(x)-＞tδ(x)vδ(x)a(x)
其中，终止符tδ(x)和vδ(x)分别表示实体x发生改变的时刻和改变量，具体情况为：船x在tδ(x)时刻发生了速度改变，大小为vδ(x)；a(x)表示实体x还能继续改变；
规则3：
a(x)-＞tδ(x)dδ(x)a(x)
其中，终止符tδ(x)和dδ(x)分别表示实体x发生改变的时刻和改变量，具体情况为：船x在tδ(x)时刻发生了转向，转向角的大小为dδ(x)；a(x)表示实体x还能继续改变；
规则4：
a(x)-＞tδ(x)aδ(x)tk(x)a(x)
其中，经过tδ(x)时间后加速度改变(默认所有模式的a＝0)，加速度改变的大小为aδ(x)即进入到加速度模型(ca)，且保持tk(x)的时间；
规则5：
a(x)-＞tδ(x)aδ(x)tk(x)b(x)
其中，经过tδ(x)时间后进入加速度模型，与规则4不同的是，维持时间结束后，舰船仍以当前的加速度继续维持上一时刻的状态，即仍保持匀加速状态继续航行；
规则6：
a(x)-＞tδ(x)wδ(x)tk(x)a(x)
其中，经过tδ(x)时间后引入角速度w，即转向角的改变为(θ＝w*t+θδ)，进入恒定转率和速度模型，并且还要维持tk(x)的时间；
规则7：
a(x1)-＞a(x1)+t(x1)[i0(x2)a(x2)+i0(x3)a(x3)
+i0(x4)a(x4)+i0(x5)a(x5)+i0(x6)a(x6)]
其中，船x1引入五艘新的舰船，分别用a(x2)、a(x3)、a(x4)、a(x5)、a(x6)表示这些实体还可继续变化，同时初始化这些新实体的初始属性分别用i0(x2)、i0(x3)、i0(x4)、i0(x5)、i0(x6)表示，而a(x1)表示实体x1还可以继续变化；
规则8：
a(x1)-＞a(x1)+[i0(x2)a(x2)+i0(x3)a(x3)]
其中，船x1引入两艘新的舰船，分别用a(x2)、a(x3)表示这些实体还可继续变化，同时初始化这些新实体的初始属性分别用i0(x2)、i0(x3)表示，而a(x1)表示实体x1还可以继续变化；
规则9：
a(x1)-＞a(x1)+i0(x2)a(x2)
其中，船x1引入1艘新的舰船实体，用a(x2)表示这个新实体还可继续变化，同时初始化这个实体的初始属性用i0(x2)来表示，新舰船a(x2)的初始时间为a(x1)的初始时间，而a(x1)表示实体x1也还可以继续变化；
规则10：
a(x)-＞ε
其中，此规则意为实体x结束变化，非终止符变为终止符；
规则11：
b(x)-＞*a(x)
其中，此规则意为维持时间结束后,找到参照的运动模式；当上一规则结束后,按照上一时刻的运动模型继续运动,直到下一次规则选取的到来；若没有选择规则11,则按照上一时刻的角度和速度进行匀速直线运动。
优选地，所述s1中，文法规则包括：舰船从雷达探测边缘区域进入，在节点引入新舰船时有不同：
规则7：
a(x1)-＞a(x1)+t(x1)[i0b(x2)a(x2)+i0b(x3)a(x3)
+i0b(x4)a(x4)+i0b(x5)a(x5)+i0b(x6)a(x6)]
其中，终止符t(x1)表示过了t(x1)段时间后，船x1从雷达探测的边界点引入五艘新的舰船实体，分别用a(x2)、a(x3)、a(x4)、a(x5)、a(x6)表示这些实体还可继续变化，i0b可以从i0bn、i0bs、i0be、i0bw的初始属性中选取，而a(x1)表示实体x1还可以继续变化；
规则8：
a(x1)-＞a(x1)+t(x1)[i0b(x2)a(x2)+i0b(x3)a(x3)]
其中，终止符t(x1)表示过了t(x1)段时间后，船x1从雷达探测的边界点引入两艘新的舰船实体，分别用a(x2)、a(x3)表示这些实体还可继续变化，i0b可以从i0bn、i0bs、i0be、i0bw的初始属性中选取，而a(x1)表示实体x1还可以继续变化；
规则9：
a(x1)-＞a(x1)+t(x1)i0b(x2)a(x2)
其中，终止符t(x1)表示过了t(x1)段时间后，船x1从雷达探测的边界点引入1艘新的舰船实体，用a(x2)表示这个新实体还可继续变化，同时初始化这个实体的初始属性用i0b(x2)来表示，新舰船a(x2)的初始时间为a(x1)的初始时间加上t(x1)，而a(x1)表示实体x1也还继续变化。
优选地，所述s2中，舰船航行中的不确定因素，包括：
s21：初始状态：初始状态为舰船目标的初始运动状态信息、雷达系统的目标检测误差；
s22：实体：实体包含不同的舰船目标；
s23：事件：事件用实体状态的改变来描述，包括增加新船、改变速度、改变角度、改变弧度；
s24：关系：关系可分为实体与初始状态间的依存关系、实体与实体间的时间关系、实体与事件间的依存关系、事件间的关系、实体与参数间的归属关系、以及实体与事件间的归属关系。
优选地，所述s3中，推演转换规则包括如下小步：
s31：推演从起始节点1开始(200)，且除了起始时不能再调用规则1，推演过程是自顶向下，按层推演的；
s32：每次推演从规则2-9随机选择非终止符来调用(201)；
s33：根据新调用的规则增删非终止符的数量(202)，每当调用规则7、8、9时，会引入新的非终止符，此节点可拥有的子节点数量将会增加；
s34：将规则号作为子节点，继续推演(203)。
优选地，所述s4中，对情景框架进行读取和赋值，将情景框架中每一艘舰船目标的运动态势字串分别读取出来，获得情景字串，对其中各个用终止符表示的参数进行赋值，即完成了一个测试情景的产生，多次推演情景模型即可产生大量各种各样的情景。
实施例2：
如图1示出本发明的情景生成示意图，情景生成过程主要分为三步：
第一步：推演定义的规则，生成情景的文法“推演树”(100)；
第二步：使用规则的内容对推演树进行替换，将规则号表示的事件用文法中的终止和非终止字符表示，得到“情景框架”(101)；
第三步：将情景框架中每一艘舰船目标的运动态势字串分别读取出来，获得“情景字串”(102)，对其中各个参数(用终止符表示)进行赋值，多次推演情景模型即可产生大量各种各样的情景。
在具体实施例中：需要确定舰船航行中的不确定因素，将其分类为初始状态、实体与事件，并建立因素间的相互关系。
舰船航行中的不确定因素，包括以下四个方面：
1)初始状态：初始状态为舰船目标的初始运动状态信息(如目标的初始速度、相对雷达的径向距离、径向速度、初始运动模式，以及目标出现的初始时刻等)、雷达系统的目标检测误差(如径向距离误差、角度误差、速度误差等)。
2)实体：实体包含不同的舰船目标。
3)事件：事件用实体状态的改变来描述(如目标速度大小的改变、运动方向的改变、航迹形状的改变等)；本发明目前定义的事件类型包括：增加新船、改变速度、改变角度、改变弧度等。
4)关系：关系可分为实体与初始状态间的依存关系、实体与实体间的时间关系、实体与事件间的依存关系、事件间的关系(如时间关系、运动状态上的量值关系等)、实体与参数间的归属关系、以及实体与事件间的归属关系。
将上述等因素的事件型模型用上下文无关文法表示，本发明在上下文无关文法的基础之上，进一步为各文法增添了终止符号、非终止符号和起始符号，并引入参数(如x、x1、x2等)来表达他们与舰船目标的归属关系。
文法的各类符号定义如下：
(1)n：非终结符集合，在这里a(x)是除了初始符s以外唯一的一种非终止符，是模拟情景中的主要变量，表示船x还可以继续变化。b(x)表示当舰船结束当前运动模型的时候，找到参照的运动模型，例如继续按照上一时刻的运动模型继续运动。
(2)σ：终止符集，它包括t(x),tδ(x),i0(x),i0bn(x)、i0bs(x)、i0be(x)、i0bw(x)，ε,[,]等符号。终止符i0(x)表示实体x的其他初始数据，这些初始数据包含舰船在海域出现时的速度、朝向、位置等，其中i0bn(x)、i0bs(x)、i0be(x)、i0bw(x)分别表示为舰船从雷达探测范围边界的北、南、东、西方向进入雷达探测范围中的初始数据。
(3)s：是一个指定的文法开始的符号，s属于n。
(4)p：转换规则，共有七条，“->”表示左边的符号被替换为右边的符号，有两套转换体系，规则如下，第一套是初始化情景：
s-＞t(0)i0(x)a(x)
起始规则：公式如上，整个模型的建立从此语法开始。其中终止符t(0)为实体x的初始时间，即舰船x在海域内出现的时间；；非终止符a(x)表示实体x还可以继续变化。它引入了第一艘船x并且初始化此船的各项初始参数、包括初始时间t(0)和初始速度、航向、经纬度(用i0(x)统一表示)。
a(x)-＞tδ(x)vδ(x)a(x)
规则2：公式如上，终止符tδ(x)和vδ(x)分别表示实体x发生改变的时刻和改变量，具体情况为：船x在tδ(x)时刻发生了速度改变，大小为vδ(x)。a(x)表示实体x还能继续改变。
a(x)-＞tδ(x)dδ(x)a(x)
规则3：公式如上，终止符tδ(x)和dδ(x)分别表示实体x发生改变的时刻和改变量，具体情况为：船x在tδ(x)时刻发生了转向，转向角的大小为dδ(x)。a(x)表示实体x还能继续改变。
a(x)-＞tδ(x)aδ(x)tk(x)a(x)
规则4：公式如上，经过tδ(x)时间后加速度改变(默认所有模式的a＝0)，加速度改变的大小为aδ(x)即进入到加速度模型(ca)，且保持tδ(x)的时间。
a(x)-＞tδ(x)aδ(x)tk(x)b(x)
规则5：公式如上，经过tδ(x)时间后进入加速度模型，与规则4不同的是，维持时间结束后，舰船仍以当前的加速度继续维持上一时刻的状态(即仍保持匀加速状态继续航行)。
a(x)-＞tδ(x)wδ(x)tk(x)a(x)
规则6：公式如上，经过tδ(x)时间后引入角速度w，即转向角的改变为(θ＝w*t+θδ)，进入恒定转率和速度模型，并且还要维持tk(x)的时间。
a(x1)-＞a(x1)+t(x1)[i0(x2)a(x2)+i0(x3)a(x3)
+i0(x4)a(x4)+i0(x5)a(x5)+i0(x6)a(x6)]
规则7：公式如上，船x1引入五艘新的舰船，分别用用a(x2)、a(x3)、a(x4)、a(x5)、a(x6)表示这些实体还可继续变化，同时初始化这些新实体的初始属性分别用i0(x2)、i0(x3)、i0(x4)、i0(x5)、i0(x6)表示，而a(x1)表示实体x1还可以继续变化。
a(x1)-＞a(x1)+[i0(x2)a(x2)+i0(x3)a(x3)]
规则8：公式如上，船x1引入两艘新的舰船，分别用a(x2)、a(x3)表示这些实体还可继续变化，同时初始化这些新实体的初始属性分别用i0(x2)、i0(x3)表示，而a(x1)表示实体x1还可以继续变化。
a(x1)-＞a(x1)+i0(x2)a(x2)
规则9：公式如上，船x1引入1艘新的舰船实体，用a(x2)表示这个新实体还可继续变化，同时初始化这个实体的初始属性用i0(x2)来表示，新舰船a(x2)的初始时间为a(x1)的初始时间，而a(x1)表示实体x1也还可以继续变化。
a(x)-＞ε
规则10：公式如上，此规则意为实体x结束变化，非终止符变为终止符。
b(x)-＞*a(x)
规则11：公式如上，此规则意为维持时间结束后,找到参照的运动模式。当上一规则结束后,按照上一时刻的运动模型继续运动,直到下一次规则选取的到来。若没有选择规则11,则按照上一时刻的角度和速度进行匀速直线运动。
第二套是舰船从雷达探测边缘区域进入，与第一套规则相比只是在节点引入新舰船时有不同，这里只表现与第一套的不同：
a(x1)-＞a(x1)+t(x1)[i0b(x2)a(x2)+i0b(x3)a(x3)
+i0b(x4)a(x4)+i0b(x5)a(x5)+i0b(x6)a(x6)]
规则7：公式如上，终止符t(x1)表示过了t(x1)段时间后，船x1从雷达探测的边界点引入五艘新的舰船实体，分别用a(x2)、a(x3)、a(x4)、a(x5)、a(x6)表示这些实体还可继续变化，i0b可以从i0bn、i0bs、i0be、i0bw的初始属性中选取，而a(x1)表示实体x1还可以继续变化。
a(x1)-＞a(x1)+t(x1)[i0b(x2)a(x2)+i0b(x3)a(x3)]
规则8：公式如上，终止符t(x1)表示过了t(x1)段时间后，船x1从雷达探测的边界点引入两艘新的舰船实体，分别用a(x2)、a(x3)表示这些实体还可继续变化，i0b可以从i0bn、i0bs、i0be、i0bw的初始属性中选取，而a(x1)表示实体x1还可以继续变化。
a(x1)-＞a(x1)+t(x1)i0b(x2)a(x2)
规则9：公式如上，终止符t(x1)表示过了t(x1)段时间后，船x1从雷达探测的边界点引入1艘新的舰船实体，用a(x2)表示这个新实体还可继续变化，同时初始化这个实体的初始属性用i0b(x2)来表示，新舰船a(x2)的初始时间为a(x1)的初始时间加上t(x1)，而a(x1)表示实体x1也还可以继续变化。
如图2示出本发明的情景推演过程流程图，具体实现步骤如下：
步骤1：推演从起始节点1开始(200)，且除了起始时不能再调用规则1，推演过程是自顶向下，按层推演的。
步骤2：每次推演从规则2-9随机选择非终止符来调用(201)。
步骤3：根据新调用的规则增删非终止符的数量(202)，每当调用规则7、8、9时，会引入新的非终止符，此节点可拥有的子节点数量将会增加。
步骤4：将规则号作为子节点，继续推演(203)。
如图3示出本发明的具体推演过程图，图示是在推演过程中产生的一个例子。
在具体实施例中：
推演开始，推演树以规则1(300)为根节点，非终止符调用规则9(301)，并引入一艘新的舰船，因为规则9引入新的非终止符，所以节点9有两个节点(302)(303),第一个非终止符(302)调用规则10，转换为终止符(304)，此路停止推演，非终止符(303)调用规则2(305)，舰船发生速度的变化，并继续推演…
图4示出本发明的推演树产生派生树示意图。如图所示，将推演树中包含的情景信息进一步展开和丰富，产生了派生树作为情景框架。
在具体实施例中：
图示右边派生树中的
图示右边的派生树其实就是将左边推演树的规则展开来表示，派生完毕后情景中表达内容更加充实和逼真，能直接表达出每艘舰船目标经历了哪些运动状态改变的事件，以及不同目标间运动状态的相互关系。
图5示出本发明的遍历派生树示意图，将情景框架中每一艘舰船目标的运动态势字串分别读取出来，获得情景字串，对其中各个用终止符表示的参数进行赋值，即完成了一个测试情景的产生，多次推演情景模型即可产生大量各种各样的情景。
在具体实施例中：
遍历时前序遍历，遇到白色底色的非终止符节点时跳过，只从左向右依次记录遍历到的绿色底色的终止符节点。
遍历如图派生树的结果是一串字符串：“t(0)i0(1)tδ(1)dδ(1)εt(1)[i0(2)tδ(2)vδ(2)tδ(2)dδ(2)εt(2)[i0(3)tδ(3)dδ(3)ε]]”。这个字符串共表示三艘舰船目标的运动态势，分别为船1、船2、船3。由于“[”和“]”符号表示在此符号中间的舰船具有相同的初始时刻，故应将带“[”和“]”符号的字符串展开，展开后，三艘舰船目标的运动态势字符串为：
(1)t(0)i0(1)tδ(1)dδ(1)ε(500)
此字符串序列表示：船1在t(0)时刻在海域出现在其初始位置并以其初始速度和航向航行，经过tδ(1)段时间后发生了一次转向，转向角为dδ(1)。此后保持运动状态不变继续航行。
(2)t(0)t(1)i0(2)tδ(2)vδ(2)tδ(2)dδ(2)ε(501)
此字符串序列表示：船1出现后再经过t(1)时刻后(即t(0)+t(1)段时间后)，船2在海域出现在其初始位置并以其初始速度和航向开始航行，又经过tδ(2)时刻后，船2发生速度改变，改变值为vδ(2)，又经过tδ(2)时刻后，船2发生转向，转向角为dδ(2)。此后保持运动状态不变继续航行。
(3)t(0)t(1)t(2)i0(3)tδ(3)dδ(3)ε(502)
此字符串序列表示：船2出现后再经过t(2)时刻后(即t(0)+t(1)+t(2)段时间后)，船3在海域出现在其初始位置并以其初始速度和航向开始航行，又经过tδ(3)时刻后，船3发生转向，转向角为dδ(3)。此后保持运动状态不变继续航行。
图6示出本发明的实验结果图,(a)(b)(c)分别为无噪声、有噪声、有噪声和海杂波的情景图。由图可见模拟出真实复杂海域中多舰船目标环境，且生成舰船多为高机动目标，舰船航行运动状态改变丰富，可以很好的体现出多艘舰船航行时，彼此之间产生的运动关系。需要说明的是：本发明将真实ais数据使用数据挖掘软件进行等频分箱处理，将得出的位置、速度改变量、角度改变量、加速度改变量等参数来进行赋值。
本发明可广泛运用于目标识别与仿真场合。
以上所述仅为本发明的较佳实施例而己，并不以本发明为限制，凡在本发明的精神和原则之内所作的均等修改、等同替换和改进等，均应包含在本发明的专利涵盖范围内。
- 还没有人留言评论。精彩留言会获得点赞！