系统监测模型生成方法、处理器芯片以及工业系统与流程

1.本发明涉及监测技术领域，具体涉及一种系统监测模型生成方法、处理器芯片以及工业系统。


背景技术：

2.在工业制造系统中，一般可以通过传感器来采集系统的输入输出信号；进一步的，由于传感器一般是间隔进行数据采集的，为了能够实时获取系统的输入输出，一般会基于系统的真实模型构建对应的虚拟系统模型，该虚拟系统模型即为通过系统的输入输出时间函数来模拟系统行为的数学模型。
3.在构建虚拟系统模型时，会对虚拟系统模型中的系统参数进行辨识，常用的辨识方法如自回归滑动平均模型(autoregressive moving average model，简称arma模型)、状态空间模型以及后向arma模型等。
4.然而，在系统参数求解过程中，存在矩阵不可逆、计算过程复杂等问题，并且会因为非凸问题导致得到虚拟系统模型的精度不高。另外，若采用神经网络通过多线性拟合非线性的方法进行系统参数辨识，也会存在计算消耗过大而无法对系统参数进行实时判别的问题。


技术实现要素：

5.本发明的目的是提供了一种系统监测模型生成方法、处理器芯片以及工业系统，在系统参数辨识过程中，不会受到不可逆矩阵对算力的影响，从而能够快速地迭代收敛；降低了对设备算力的要求，应用于算力较低的处理器芯片，能够实时的对系统进行监测。
6.为实现上述目的，本发明提供了一种系统监测模型生成方法，包括：根据待监测系统的真实数学模型，构建所述待监测系统对应虚拟系统的数字孪生模型，所述数字孪生模型中包括所述虚拟系统的多个虚拟系统参数；构建所述真实数学模型与所述数字孪生模型之间的误差反馈模型；基于所述待监测系统的参考历史数据与所述误差反馈模型进行回归分析，得到各所述虚拟系统参数的值。
7.本发明还提供了一种处理器芯片，处理器芯片用于执行上述系统监测模型生成方法生成待监测系统对应的虚拟系统的数字孪生模型。
8.本发明还提供了一种工业系统，包括上述的处理器芯片与连接于所述处理器芯片的工业设备；所述处理器芯片用于生成所述工业设备对应的虚拟系统的数字孪生模型。
9.本发明实施例中，在生成待监测系统所对应的虚拟系统的数字孪生模型时，首先根据待监测系统的真实数学模型，构建所述待监测系统对应虚拟系统的数字孪生模型，所述数字孪生模型中包括所述虚拟系统的多个虚拟系统参数，然后构建所述真实数学模型与所述数字孪生模型之间的误差反馈模型，再基于所述待监测系统的参考历史数据与所述误差反馈模型进行回归分析，得到各所述虚拟系统参数的值，从而能够得到待监测系统的数字孪生模型；采用本发明的系统监测模型生成方法，在系统参数辨识过程中，不会受到不可
逆矩阵对算力的影响，从而能够快速地迭代收敛；降低了对设备算力的要求，应用于算力较低的处理器芯片，能够实时的对系统进行监测。
10.在一个实施例中，所述待监测系统的真实数学模型的表达式为：
11.x(t k)＝f[x(t k
‑
1)，x(t k
‑
2)，
…
，x(t)，u(t)]，其中，x(t)表示t时刻所述待监测系统的输出测量值，k表示预设的第一滑窗长度，u(t)表示t时刻所述待监测系统的输入值；
[0012]
所述待监测系统对应的数字孪生模型的表达式为：
[0013][0014]
其中，表示n时刻所述待监测系统的输出预测值，x(n)表示n时刻所述待监测系统的输出测量值，ω
i
、ω
i
′
以及ω0为所述虚拟系统参数，l表示预设的第二滑窗长度，u(n)表示n时刻所述待监测系统的输入值，k0表示第一预设常数。
[0015]
在一个实施例中，所述误差反馈模型的表达式为：
[0016][0017]
其中，e(n)表示n时刻所述真实数学模型与所述数字孪生模型之间的误差，表示n时刻所述待监测系统的输出预测值，x(n)表示n时刻所述待监测系统的输出测量值，k表示所述第一滑窗长度，l表示所述第二滑窗长度。
[0018]
在一个实施例中，所述基于所述待监测系统的参考历史数据与所述误差反馈模型进行回归分析，得到各所述虚拟系统参数的值，包括：
[0019]
利用所述待监测系统的参考历史数据，采用单次测量天牛须搜索算法对所述误差反馈模型进行回归分析，得到各所述虚拟系统参数的值。
[0020]
在一个实施例中，所述利用所述待监测系统的参考历史数据，采用单次测量天牛须搜索算法对所述误差反馈模型进行回归分析，得到各所述虚拟系统参数的值，包括：
[0021]
利用所述多个虚拟系统参数构建天牛的位置，将所述误差反馈模型作为所述单次测量天牛须搜索算法的适应度函数，并对天牛进行初始化；
[0022]
基于所述历史参考数据、天牛的当前位置以及所述适应度函数，得到天牛在当前位置的适应度；
[0023]
计算天牛左右两个触须的适应度，并更新天牛的位置；
[0024]
更新天牛的步长、须长以及迭代次数，并判断是否满足预设的迭代条件；
[0025]
若满足预设的迭代条件，终止所述单次测量天牛须搜索算法，并输出当前的天牛位置作为虚拟系统参数的值；
[0026]
若不满足预设的迭代条件，将更新后的天牛的位置代入到所述适应度函数中，重复所述单次测量天牛须搜索算法。
[0027]
在一个实施例中，所述计算天牛左右两个触须的适应度，包括：
[0028]
基于天牛的当前位置、天牛左右触须的长度，得到天牛左右触须的位置；
[0029]
将所述历史参考数据中最近一个测量数据与天牛左触须的位置代入所述适应度函数，得到天牛左触须的适应度；
[0030]
将所述历史参考数据中所述最近一个测量数据的前一个测量数据与天牛右触须的位置代入所述适应度函数，得到天牛右触须的适应度。
[0031]
在一个实施例中，天牛的位置更新公式为：
[0032]
t
i 1
＝t
i
‑
step*g/σ；
[0033]
其中，t
i
表示第i次迭代天牛的位置，step表示天牛的当前步长，g表示单次测量天牛须搜索算法的梯度，σ表示第二预设常数。
[0034]
在一个实施例中，单次测量天牛须搜索算法的梯度的公式为：
[0035]
g＝sign(e(t
i
)
‑
e(t
i
‑1))*normal(t
i
‑
t
i
‑1)；
[0036]
其中，g表示单次测量天牛须搜索算法的梯度，sign表示天牛下一步的方向，t
i
表示第i次迭代天牛的位置，e(t
i
)表示天牛在t
i
位置时的适应度。
附图说明
[0037]
图1是根据本发明第一实施例的系统监测模型生成方法的框架图；
[0038]
图2是根据本发明第一实施例的系统监测模型生成方法的具体流程图；
[0039]
图3是图2的系统监测模型生成方法的步骤103的具体流程图；
[0040]
图4是根据本发明第一实施例的数据孪生模型在基于单次测量天牛须搜索算法进行回归分析的示意图；
[0041]
图5是根据本发明第三实施例的工业系统的示意图。
具体实施方式
[0042]
以下将结合附图对本发明的各实施例进行详细说明，以便更清楚理解本发明的目的、特点和优点。应理解的是，附图所示的实施例并不是对本发明范围的限制，而只是为了说明本发明技术方案的实质精神。
[0043]
在下文的描述中，出于说明各种公开的实施例的目的阐述了某些具体细节以提供对各种公开实施例的透彻理解。但是，相关领域技术人员将认识到可在无这些具体细节中的一个或多个细节的情况来实践实施例。在其它情形下，与本技术相关联的熟知的装置、结构和技术可能并未详细地示出或描述从而避免不必要地混淆实施例的描述。
[0044]
除非语境有其它需要，在整个说明书和权利要求中，词语“包括”和其变型，诸如“包含”和“具有”应被理解为开放的、包含的含义，即应解释为“包括，但不限于”。
[0045]
在整个说明书中对“一个实施例”或“一实施例”的提及表示结合实施例所描述的特定特点、结构或特征包括于至少一个实施例中。因此，在整个说明书的各个位置“在一个实施例中”或“在一实施例”中的出现无需全都指相同实施例。另外，特定特点、结构或特征可在一个或多个实施例中以任何方式组合。
[0046]
如该说明书和所附权利要求中所用的单数形式“一”和
“”
包括复数指代物，除非文中清楚地另外规定。应当指出的是术语“或”通常以其包括“或/和”的含义使用，除非文中清楚地另外规定。
[0047]
在以下描述中，为了清楚展示本发明的结构及工作方式，将借助诸多方向性词语进行描述，但是应当将“前”、“后”、“左”、“右”、“外”、“内”、“向外”、“向内”、“上”、“下”等词语理解为方便用语，而不应当理解为限定性词语。
[0048]
本发明第一实施方式涉及一种系统监测模型生成方法，应用于处理器芯片，该处理器芯片能够执行本实施例中的系统监测模型生成方法生成待监测系统对应的虚拟系统
的数字孪生模型，该数字孪生模型能够基于待监测系统的输入输出时间函数来模拟系统的行为；其中，待监测系统可以为工业系统，例如变压器系统、风电系统等。
[0049]
本实施例的系统监测模型生成方法，包括：根据待监测系统的真实数学模型，构建待监测系统对应虚拟系统的数字孪生模型，数字孪生模型中包括虚拟系统的多个虚拟系统参数；构建真实数学模型与数字孪生模型之间的误差反馈模型；基于待监测系统的参考历史数据与误差反馈模型进行回归分析，得到各虚拟系统参数的值。
[0050]
采用本实施例中的系统监测模型生成方法，在系统参数辨识过程中，不会受到不可逆矩阵对算力的影响，从而能够快速地迭代收敛；降低了对设备算力的要求，应用于算力较低的处理器芯片，能够实时的对系统进行监测。
[0051]
下面结合待监测系统的具体数学模型进行详细说明，仅为示例性说明，并非实现本发明的必须。其中，图1为系统监测模型生成方法的框架图。
[0052]
本实施方式的系统监测模型生成方法的具体流程如图2所示。
[0053]
步骤101，根据待监测系统的真实数学模型，构建待监测系统对应虚拟系统的数字孪生模型，数字孪生模型中包括虚拟系统的多个虚拟系统参数。
[0054]
待监测系统的真实数学模型的表达式为：
[0055]
x(t k)＝f[x(t k
‑
1)，x(t k
‑
2)，
…
，x(t)，u(t)]，其中，x(t)表示t时刻待监测系统的输出测量值，k表示预设的第一滑窗长度，u(t)表示t时刻待监测系统的输入值。以待监测系统为变压器系统为例，变压器系统的输入u为输入到变压器系统中的电压，变压器系统的输出x为变压器系统中的传感器测量的变压器系统输出的电压。
[0056]
待监测系统对应的数字孪生模型的表达式为：
[0057][0058]
其中，表示n时刻待监测系统的输出预测值，x(n)表示n时刻待监测系统的输出测量值，ω
i
、ω
i
′
以及ω0为虚拟系统参数，l表示预设的第二滑窗长度，u(n)表示n时刻待监测系统的输入值，k0表示第一预设常数。其中，第一预设常数k0一般设置为远大于0的常数，例如100、1000、2000等。
[0059]
本实施例中，在待监测系统的真实数学模型以及与待监测系统对应的虚拟系统的数字孪生模型中加入了滑动窗口，能够增加待监测系统的真实数学模型以及数字孪生模型对系统监测的完整性，使得模型能够更加准确的反应系统的输入输出。
[0060]
步骤102，构建真实数学模型与数字孪生模型之间的误差反馈模型。
[0061]
本实施例中采用滑动平均误差来构建误差反馈机制，构建得到的真实数学模型与数字孪生模型之间的误差反馈模型的表达式为：
[0062][0063]
其中，e(n)表示n时刻真实数学模型与数字孪生模型之间的误差，表示n时刻待监测系统的输出预测值，x(n)表示n时刻待监测系统的输出测量值，k表示第一滑窗长度，l表示第二滑窗长度。
[0064]
步骤103，利用待监测系统的参考历史数据，采用单次测量天牛须搜索算法对误差反馈模型进行回归分析，得到各虚拟系统参数的值。
[0065]
请参考图3，步骤103包括以下子步骤：
[0066]
子步骤1031，利用多个虚拟系统参数构建天牛的位置，将误差反馈模型作为单次测量天牛须搜索算法的适应度函数，并对天牛进行初始化。
[0067]
具体而言，本实施例中选择步骤102中构建的误差反馈模型e(n)作为单次测量天牛须搜索算法(one measure beetle antennae search algorithm，简称ombas)的适应度函数(即待优化的目标函数)，然后利用数字孪生模型中所包含的虚拟系统参数ω
i
、ω
i
′
以及ω0构建ombas算法中天牛的位置，即天牛在第i次迭代中的位置t
i
＝(ω
i
，ω
i
′
，ω0)，并对天牛进行初始化，初始化过程包括：设置ω
i
、ω
i
′
以及ω0的范围，然后可以从这个范围中随机产生天牛位置，同时还包括对迭代次数的初始化，即设置迭代次数i＝1，以及天牛左右触须的长度d、步长step。
[0068]
子步骤1032，基于历史参考数据、天牛的当前位置以及适应度函数，得到天牛在当前位置的适应度。
[0069]
具体而言，将天牛的当前位置t
i
带入到数字孪生模型和适应度函数e(n)中，然后将历史参考数据中各个时刻的输入信号代入到数字孪生模型中，得到数字孪生模型在各个时刻的输出预测值，然后基于适应度函数e(n)的表达式与数字孪生模型的输出预测值，便能够计算到天牛在当前位置t
i
的适应度。
[0070]
子步骤1033，计算天牛左右两个触须的适应度，并更新天牛的位置。
[0071]
具体而言，计算天牛左右两个触须的适应度，包括：基于天牛的当前位置、天牛左右触须的长度，得到天牛左右触须的位置；将历史参考数据中最近一个测量数据与天牛左触须的位置代入适应度函数，得到天牛左触须的适应度；将历史参考数据中最近一个测量数据的前一个测量数据与天牛右触须的位置代入适应度函数，得到天牛右触须的适应度。即只需一个传感器提供的历史参考数据便能够运行单次测量天牛须搜索算法，减少了传感器的测量次数，进一步提高了历史参考数据的利用率及算法效率。
[0072]
在单次测量天牛须搜索算法中，天牛左触须与右触须之间的距离d0是固定的，可以用xl表示天牛左触须的坐标，xr表示右须坐标，x表示质心坐标(即为天牛当前的位置)，并生成随机向量dir＝rands(n，1)来表示天牛进行下一步飞行后的天牛头的朝向，n为大于或等于1的正整数，即天牛右触须指向天牛左触须的向量的朝向是任意的，归一化后得到dir＝dir/norm(dir)；
[0073]
xl
‑
xr＝d0*dir；
[0074]
因此，xl、xr还可以表示成关于质心x的表达式：
[0075]
xl＝x d0*dir/2；
[0076]
xr＝x
‑
d0*dir/2。
[0077]
随后将天牛左触须的坐标xl带入到适应度函数e(n)中，利用历史参考数据最近一个测量数据得到天牛左触须的适应度fleft；同理，将天牛右触须的坐标xr带入到适应度函数e(n)中，利用历史参考数据最近一个测量数据的前一个测量数据得到天牛右触须的适应度fright。其中，天牛左右触须的适应度的计算方式与天牛的适应度的计算方式大致相同，在此不再赘述。
[0078]
将天牛左触须的适应度fleft与天牛右触须的适应度ffright进行比较，能够得到天牛下一步的前进方向，如果fleft＜fright，为了探寻天牛的适应度的最小值，则天牛向着左触须方向行进，如果fleft＞fright，为了探寻天牛的适应度的最小值，则天牛向着右
触须方向行进，可以采用如下的符号函数sign来表示：
[0079][0080]
其中，sign＝1时，表示下一步天牛向着右触须方向行进；sign＝
‑
1时，表示下一步天牛向着左触须方向行进。
[0081]
在确定天牛下一步的行进方向之后，便可以根据天牛当前的步长step来更新天牛的位置，天牛的位置更新公式为
[0082]
t
i 1
＝t
i
‑
step*g/σ；
[0083]
其中，t
i
表示第i次迭代天牛的位置，step表示天牛的当前步长，g表示单次测量天牛须搜索算法的梯度，σ表示第二预设常数。
[0084]
单次测量天牛须搜索算法的梯度的公式为：
[0085]
g＝sign(e(t
i
)
‑
e(t
i
‑1))*normal(t
i
‑
t
i
‑1)；
[0086]
g表示单次测量天牛须搜索算法的梯度，sign表示天牛下一步的方向，t
i
表示第i次迭代天牛的位置，e(t
i
)表示天牛在t
i
位置时的适应度。
[0087]
子步骤1034，更新天牛的步长、须长以及迭代次数，并判断是否满足预设的迭代条件。若是，则进入子步骤1034；若否，则进入子步骤1035。
[0088]
子步骤1035，终止单次测量天牛须搜索算法，并输出当前的天牛位置作为虚拟系统参数的值。
[0089]
子步骤1036，将更新后的天牛的位置代入到适应度函数中，重复单次测量天牛须搜索算法。
[0090]
天牛步长step的更新公式为：
[0091]
step＝step*eta；
[0092]
其中，eta表示天牛的步长调整比例，0＜eta＜1，例如eta取0.95。
[0093]
天牛须长p的更新公式为：
[0094]
p
i
＝r*p
i
‑
1；
[0095]
其中，r表示天牛的须长衰减系数，i表示迭代次数。
[0096]
迭代次数的更新方式为i＝i 1。
[0097]
在完成上述的天牛的步长、须长以及迭代次数的更新之后，判断本次迭代是否满足预设的迭代条件，迭代条件可以为天牛的适应度是否满足预设的适应度阈值，和/或天牛的迭代次数i是否达到预设的迭代次数阈值，当本次迭代满足迭代条件时，终止本次的单次测量天牛须搜索算法，并将当前天牛的位置t
i
中的ω
i
，ω
i
′
以及ω0的值输入到数字孪生模型中，该数字孪生模型便能够基于时间来模拟待监测系统的输出。当本次迭代不满足迭代条件时，将更新后的天牛的位置代入到适应度函数中，回到步骤103，重复上述过程，直至迭代至满足预设的迭代条件。
[0098]
在一个实施例中，还可以在得到待监测系统所对应的虚拟系统的数字孪生模型后，虚拟系统的误差e
n
为：
[0099][0100]
虚拟系统的误差更新方式如下：
[0101]
[0102]
由此可以得到虚拟系统的误差判别公式：
[0103]
|e
n
|＜k1*e
′
n
；
[0104]
其中，∈与k1为待定参数，此时可以利用参数∈与k1构成天牛的位置，并采用天牛须搜索算法来整定参数∈与k1，便能够得到虚拟系统的误差判别公式，具体过程与上述类似，在此不再赘述。当虚拟系统的误差符合误差判别公式时，说明虚拟系统的误差较小，该虚拟系统符合要求；当虚拟系统的误差不符合误差判别公式时，说明虚拟系统的误差较大，该虚拟系统不符合要求。
[0105]
请参考图4，包括了一种数据孪生模型在基于单次测量天牛须搜索算法进行回归分析的收敛速度图、精度误差图以及损失图，由图4可见，单次测量天牛须搜索算法在迭代了112次后收敛，损失值为0.085，精度误差在0.2左右。
[0106]
本实施例中，处理器芯片在利用待监测系统所对应的虚拟系统的数字孪生模型后，能够结合待监测系统中的传感器测量的数据与数字孪生模型反馈的数据进行对比，从而在虚拟系统与真实系统之间的误差较大时，利用最近的待监测系统中的传感器测量的数据重新执行上述的系统监测模型生成方法，生成一个新的数字孪生模型，从而能够适用于待监测系统在各个阶段的实时状态。
[0107]
由上可知，本实施例中提出了基于单次测量天牛须搜索算法来生成数字孪生模型，采用单次测量天牛须搜索算法回归分析方法，对数字孪生模型中的虚拟系统参数进行辨识，实现了算力消耗小、孪生精度高、收敛速度快的效果。
[0108]
本发明的第二实施例涉及一种处理器芯片，处理器芯片用于执行第一实施方式中的系统监测模型生成方法生成待监测系统对应的虚拟系统的数字孪生模型，该数字孪生模型能够基于待监测系统的输入输出时间函数来模拟系统的行为；其中，待监测系统可以为工业系统，例如变压器系统、风电系统等。
[0109]
本发明第三实施例涉及一种工业系统，请参考图5，包括第二实施例中的处理器芯片1与连接于处理器芯片1的工业设备2；工业设备例如为变压器系统、风电系统等。
[0110]
处理器芯片1用于生成工业设备2对应的虚拟系统的数字孪生模型。
[0111]
以上已详细描述了本发明的较佳实施例，但应理解到，若需要，能修改实施例的方面来采用各种专利、申请和出版物的方面、特征和构思来提供另外的实施例。
[0112]
考虑到上文的详细描述，能对实施例做出这些和其它变化。一般而言，在权利要求中，所用的术语不应被认为限制在说明书和权利要求中公开的具体实施例，而是应被理解为包括所有可能的实施例连同这些权利要求所享有的全部等同范围。