一种连续时间系统的并行仿真及误差补偿方法技术方案

本发明专利技术公开的一种连续时间系统的并行仿真及误差补偿方法，属于复杂系统仿真技术领域。本发明专利技术首先建立连续时间系统的雅可比矩阵，然后将雅可比矩阵转化为超图形式表达，利用超图分割策略将所述超图分割成相互间耦合度尽可能低的多个部分，再将超图分割结果映射到连续时间系统中以实现将连续时间系统分割成多个子连续时间系统且子连续时间系统间耦合低。最后各子连续时间系统被分布在不同处理机上并行仿真，每个子连续时间系统在仿真过程中采用动态选取的插值方式进行输入量的插值以补偿误差。本发明专利技术具有误差补偿策略通用性好，能够缩短连续时间系统求解耗时，能降低单个计算机的负载的优点，适用于工程领域连续时间系统仿真。

【技术实现步骤摘要】
一种连续时间系统的并行仿真及误差补偿方法
本专利技术涉及一种针对于连续时间系统的并行仿真及误差补偿方法，属于复杂系统仿真

技术介绍
20世纪以来，随着计算机、互联网、可视化等技术的飞速发展与融合，仿真技术在工业产品的研制过程中获得了大规模应用，尤其在设计、制造和综合保障等多个关键研制环节发挥着举足轻重的作用。连续时间系统广泛见于工业各领域的建模与仿真中，通常以常微分方程组(OrdinaryDifferentialEquations，ODEs)或微分代数方程(DifferentialAlgebraicEquation，DAE)的形式表达。随着系统越来越复杂，系统求解的耗时问题日益突出，尤其在实时仿真领域这一问题尤其致命。如何提高仿真的效率成为迫切需要解决的问题。并行仿真是解决这一问题的有效方法之一。并行仿真的主要任务是将仿真目标分解为多个子目标并将其分布在不同处理机上同时仿真，从而提高仿真的效率。将一个串行的仿真问题转化成并行仿真问题的过程称为仿真并行化，针对连续时间系统的仿真并行化研究尚有不足。此外，仿真并行化过程中由于子目标间数据交换不频繁会产生额外的时延误差。针对于此，现有的误差补偿方法普遍采用单一的插值策略：或者是常量插值、或者是线性插值等，这一单一策略无法适应系统不同时段的变化要求，更无法适应不同系统间的变化。
技术实现思路
针对现有技术存在的下述不足：首先是没有用于连续时间系统并行仿真的系统性方法，其次是误差补偿策略不够通用。本专利技术公开的一种连续时间系统的并行仿真及误差补偿方法要解决的技术问题是：(1)将串行连续时间系统仿真问题转化为并行仿真问题，缩短连续时间系统求解耗时；(2)通过更加通用且合理的误差补偿策略尽量减小由串行问题转化成并行问题导致的额外误差，使误差补偿策略更通用。本专利技术的目的是通过下述技术方案实现的。本专利技术公开的一种连续时间系统的并行仿真及误差补偿方法，首先建立连续时间系统的雅可比矩阵，然后将雅可比矩阵转化为超图形式表达，利用超图分割策略将所述超图分割成相互间耦合度尽可能低的多个部分，再将超图分割结果映射到连续时间系统中以实现将连续时间系统分割成多个子连续时间系统且子连续时间系统间耦合低。最后各子连续时间系统被分布在不同处理机上并行仿真，每个子连续时间系统在仿真过程中采用动态选取的插值方式进行输入量的插值以补偿误差。所述方法具有误差补偿策略通用性好，能够缩短连续时间系统求解耗时，能降低单个计算机的负载的优点，适用于工程领域连续时间系统仿真。本专利技术公开的一种连续时间系统的并行仿真及误差补偿方法，包括以下步骤：步骤1：建立连续时间系统的雅可比矩阵。步骤1具体实现方法如下：连续时间系统用式(1)表示：式(1)中fi(i＝1,2,...,n)为状态转移函数，xi(i＝1,2,...,n)为状态变量。那么连续时间系统的雅可比矩阵表示为：矩阵中为0的元素(即)表示fi不依赖于状态变量xj，也即是的求解不依赖于状态变量xj，反之表示依赖于xj。即实现直观表达连续时间系统中的数据依赖关系。步骤2：将步骤1中雅可比矩阵转换为超图形式表示。步骤2具体实现方法如下：式(2)所示的雅可比矩阵的每列对应超图的顶点，每行对应超图的超边结点。若则将对应的顶点i与对应的超边结点j相连接从而构成超边的一条边，顶点集合和超边集合共同构成完整的超图。步骤3：利用超图划分策略将所有的超边结点划分为多个超边结点子集合。步骤3具体实现步骤如下：步骤3.1：利用超图划分策略将所有超边结点划分为多个完全独立的超边结点子集合以及耦合超边结点集合。所述完全独立的超边结点子集合是指与子集合中所有超边结点相连接的顶点的集合不会发生交叉，即不存在一个顶点能够同时连接到属于两个完全独立的超边结点子集合的两个超边结点。所述耦合超边结点是指该超边结点不能归属于任一独立的超边结点子集合，否则会造成独立性不成立。耦合超边结点集合是原系统内部的最基本耦合关系。完全独立的超边结点子集合允许有多个，而耦合超边结点集合只有一个。对于无法直接计算的复杂的连续时间系统，步骤3.1优选借助PaToH工具实现利用超图划分策略将所有超边结点划分为多个完全独立的超边结点子集合以及耦合超边结点集合。步骤3.2：将耦合超边结点分配到独立的各超边结点子集合中，无需特别的分配规则，但应尽量使得各子集合中超边结点数量均匀。步骤4：将步骤3划分的结果映射到原连续时间系统中，使得原连续时间系统分为多个子连续时间系统。步骤5：对子连续时间系统分别建模以得到各子模型并将其分布到不同处理机上。步骤6：确定子模型间的输入输出依赖关系也即数据流。步骤7：各子模型在不同处理机上并行仿真。对每个子模型仿真处理，具体实现方法包括如下步骤：步骤7.1：模型初始化。所述初始化内容包括：当前仿真时钟、状态变量及输入量初始值、当前所选取插值函数。步骤7.2：判断当前仿真时刻是否为仿真结束时刻，如果是，则仿真结束；如果不是则继续。步骤7.3：判断当前仿真时刻是否为通信点。如果是，则继续；如果不是，跳至步骤7.6。所述通信点是子模型间数据交换的时间点。在并行仿真中，各子模型的依赖关系被限制在离散的通信点上，子模型的数据交换只发生在通信点上，两个通信点之间(称之为“通信区间”)各子模型并行独立求解。通信区间的长度为通信步长，通信步长通常为子模型仿真步长的n倍。步骤7.4：将各输出量的值及其当前时刻的导数传给下游模型以更新下游模型的输入；同时检查自身输入是否已经更新，若没有则等待上游模型更新，若已更新则继续。所述输出量当前时刻导数通过最后两次计算值的差除以仿真步长得出。步骤7.5：根据当前通信时刻和上一通信时刻各输入量的导数，确定当前各输入量采用的插值方式。步骤7.5为误差补偿策略，步骤7.5具体实现方法如下：子模型的每个输入ui(i＝1,2,...,m)当前通信时刻的导数为上一通信时刻的导数为规定：当时，用符号“\”表示其变化性；时，用符号“＿”表示其变化性；时，用符号“/”表示其变化性。则该通信区间内各输入量的插值方式通过表1选取。所述的a和b根据要求设定，优选a＝-0.2，b＝0.2。表1插值方式选取表步骤7.6：各输入量采用其当前插值方式插值获取当前时刻的值。步骤7.7：模型仿真一个仿真步长，同时仿真时钟推进一个仿真步长。然后跳至步骤7.2。步骤8：根据步骤1到步骤7将连续时间系统分割成若干子连续时间系统后并行仿真，将串行连续时间系统仿真问题转化为并行仿真问题，缩短连续时间系统求解耗时长；根据步骤7.5所述的误差补偿策略尽量减小由串行问题转化成并行问题导致的额外误差，使误差补偿策略更通用。还包括步骤9：将步骤1至8所述一种连续时间系统的并行仿真及误差补偿方法应用于工程领域连续时间系统仿真，解决相应工程问题，所述的工程领域连续时间系统仿真包括：弹道仿真、结构模态仿真等。有益效果：1.本专利技术公开的一种连续时间系统的并行仿真及误差补偿方法，将现有技术中单一的插值策略转变为步骤7.5所示的动态选择的插值策略，即实现更加合理的误差补偿策略，尽量减小由串行问题转化成并行问题导致的额外误差，使误差补偿策略更通用。2.本专利技术公开的一种连续时间系统的并行仿真及误差补偿方法，采用雅可比矩阵和超图划分策略本文档来自技高网...

【技术保护点】
1.一种连续时间系统的并行仿真及误差补偿方法，其特征在于：包括以下步骤：步骤1：建立连续时间系统的雅可比矩阵；步骤2：将步骤1中雅可比矩阵转换为超图形式表示；步骤3：利用超图划分策略将所有的超边结点划分为多个超边结点子集合；步骤4：将步骤3划分的结果映射到原连续时间系统中，使得原连续时间系统分为多个子连续时间系统；步骤5：对子连续时间系统分别建模以得到各子模型并将其分布到不同处理机上；步骤6：确定子模型间的输入输出依赖关系也即数据流；步骤7：各子模型在不同处理机上并行仿真；步骤8：根据步骤1到步骤7将连续时间系统分割成若干子连续时间系统后并行仿真，将串行连续时间系统仿真问题转化为并行仿真问题，缩短连续时间系统求解耗时长。

【技术特征摘要】
1.一种连续时间系统的并行仿真及误差补偿方法，其特征在于：包括以下步骤：步骤1：建立连续时间系统的雅可比矩阵；步骤2：将步骤1中雅可比矩阵转换为超图形式表示；步骤3：利用超图划分策略将所有的超边结点划分为多个超边结点子集合；步骤4：将步骤3划分的结果映射到原连续时间系统中，使得原连续时间系统分为多个子连续时间系统；步骤5：对子连续时间系统分别建模以得到各子模型并将其分布到不同处理机上；步骤6：确定子模型间的输入输出依赖关系也即数据流；步骤7：各子模型在不同处理机上并行仿真；步骤8：根据步骤1到步骤7将连续时间系统分割成若干子连续时间系统后并行仿真，将串行连续时间系统仿真问题转化为并行仿真问题，缩短连续时间系统求解耗时长。2.如权利要求1所述的一种连续时间系统的并行仿真及误差补偿方法，其特征在于：还包括步骤9：将步骤1至8所述一种连续时间系统的并行仿真及误差补偿方法应用于工程领域连续时间系统仿真，解决相应工程问题，所述的工程领域连续时间系统仿真包括弹道仿真、结构模态仿真。3.如权利要求1或2所述的一种连续时间系统的并行仿真及误差补偿方法，其特征在于：步骤1具体实现方法如下，连续时间系统用式(1)表示：式(1)中fi(i＝1,2,...,n)为状态转移函数，xi(i＝1,2,...,n)为状态变量；那么连续时间系统的雅可比矩阵表示为：矩阵中为0的元素(即)表示fi不依赖于状态变量xj，也即是的求解不依赖于状态变量xj，反之表示依赖于xj；即实现直观表达连续时间系统中的数据依赖关系。4.如权利要求3所述的一种连续时间系统的并行仿真及误差补偿方法，其特征在于：步骤2具体实现方法如下，式(2)所示的雅可比矩阵的每列对应超图的顶点，每行对应超图的超边结点；若则将对应的顶点i与对应的超边结点j相连接从而构成超边的一条边，顶点集合和超边集合共同构成完整的超图。5.如权利要求4所述的一种连续时间系统的并行仿真及误差补偿方法，其特征在于：步骤3具体实现步骤如下，步骤3.1：利用超图划分策略将所有超边结点划分为多个完全独立的超边结点子集合以及耦合超边结点集合；所述完全独立的超边结点子集合是指与子集合中所有超边结点相连接的顶点的集合不会发生交叉，即不存在一个顶点能够同时连接到属于两个完全独立的超边结点子集合的两个超边结点；所述耦合超边结点是指该超边结点不能归属于任一独立的超边结点子集合，否则会造成独立性不成立；耦合超边结点集合是原系统内部的最基本耦合关系；完全独立的超边结点子集合...

【专利技术属性】
技术研发人员：李文光，王强，刘莉，
申请(专利权)人：北京理工大学，
类型：发明
国别省市：北京,11