基于历史数据斜坡响应的动态系统静态增益估计方法技术方案

本发明专利技术涉及基于历史数据斜坡响应的动态系统静态增益估计方法。该方法首先要完成对历史数据的分析，采用分段线性表示方法将输入和输出的时间序列分割成短数据段。其次，找出输入和输出同时呈直线且处于振幅变化较大的斜坡响应，较大幅度变化的阈值是根据振幅变化与确定系数之间的关系确定的。最后通过求解斜坡响应幅值变化较大的多个线性方程，得到估计的静态增益。本发明专利技术可以通过可视化验证该方法的有效性，克服了使用系统识别方法难以验证静态增益估计的问题。从而避免了估计的静态增益和实际增益之间偏差无法判断的情况。

【技术实现步骤摘要】
基于历史数据斜坡响应的动态系统静态增益估计方法
本公开涉及工业大数据分析
，特别是涉及基于历史数据斜坡响应的动态系统静态增益估计方法。
技术介绍
本部分的陈述仅仅是提高了与本公开相关的
技术介绍
，并不必然构成现有技术。静态增益是动态系统的重要信息，动态系统的静态增益通常用于设计反馈控制器、监测过程变化和优化运行性能，静态增益可在动态系统设计阶段获得，而设计阶段和实际操作的条件不同，因此，动态系统的静态增益在实践中往往是不可用的，必须从观察到的数据样本中估计。一种常用的方法就是从输入和输出的稳态值或某些特殊类型的测试中估计动态系统的静态增益。现有方法存在两个局限性。首先，稳态条件下的数据样本很难获得，对于某些动态系统，不允许做特殊类型的测试。因此，从日常操作收集的数据样本来估计静态增益是可取的。第二，系统辨识技术是基于一个假设条件，假设模型集足够丰富，足以包含真实模型；而且，这种假设在实践中无法证实，也无法找到估计的静态增益和实际增益之间的偏差。
技术实现思路
为了解决现有技术的不足，本公开提供了基于历史数据斜坡响应的动态系统静态增益估计方法，从隐藏在工业大数据中的相关信息估计动态系统的静态增益。此方法不依赖数据分析，而且也避免了估计的静态增益和实际增益之间偏差无法判断的情况。在模型估计方面有很好的应用价值，可以克服系统辨识的不确定性。第一方面，本公开提供了基于历史数据斜坡响应的动态系统静态增益估计方法；基于历史数据斜坡响应的动态系统静态增益估计方法，包括：采用分段线性表示方法将动态系统的输出时间序列分割成若干个输出时间段；同时，将动态系统的输入时间序列分割成若干个输入数据段；每个数据段均用一条直线表示；找出输入值和输出值同时均处于各自的直线数据段上且输入值和输出值的振幅变化均超过设定阈值的数据段，此数据段即为目标斜坡响应，将找出的所有输入数据段和所有输出数据段组成具有显著振幅变化的数据段集合；通过具有显著振幅变化的数据段集合，估计静态增益。作为一种可能的实现方式，所述将动态系统的输出时间序列分割成若干个输出时间段的具体步骤为：将时间序列分离为M个数据段对于第m个数据段nm表示第m段中第一个数据样本，m∈[1，M]，nm+1-nm-1为第m段数据段中的样本总数；用线性回归模型来描述：y(n)＝am+bmn+e(n)；其中，am表示第m段数据段的初始值，bm表示第m段数据段的斜率，e(n)是干扰；假设e(n)是具有零均值和方差为的白噪声；得到拟合线性方程：其中，和分别是y(n)、am和bm的估计值，用M表示数据段个数；M的估计值其中，L(M)为拟合误差损失函数：作为一种可能的实现方式，将动态系统的输入时间序列分割成若干个输入数据段：采用分段线性表示方法将输入的时间序列分割为M个数据段时间序列的分段个数M，M的估计值：作为一种可能的实现方式，找出输入值和输出值同时均处于各自的直线数据段上且输入值和输出值的振幅变化均超过设定阈值的数据段，此数据段即为目标斜坡响应，将找出的所有输入数据段和所有输出数据段组成具有显著振幅变化的数据段集合的具体步骤为：对于输出来说，计算第m段的振幅变化量：引入序列：其中，Ay,0为y显著幅度变化的阈值；对于输入来说，计算每个数据段的振幅变化值同样引入序列：其中，为ui显著幅度变化的阈值；整体序列：保留整体序列I0(n)＝1的数据段，当I0(n)＝1时，表示输入输出同时处于幅度变化超过设定阈值且同时处于直线数据段上，将有显著振幅变化的数据段组成一个集合作为一种可能的实现方式，通过具有显著振幅变化的数据段集合，估计静态增益的具体步骤为：步骤3.1：通过最小二乘法估计的静态增益K1,K2,……KI；K是由K1,K2,……KI组成的I维向量，K的估计是由最小二乘法得到的：从K的高斯分布，估计K的置信区间；步骤3.2：找出中Ay,l与其估计值偏差最大的数据段第l0段，如果该段幅值变化满足不等式：则将从中移除，其中δy是用户选择的参数，表示和最大偏差的可接受水平。步骤3.3：重复步骤3.1和3.2，直到找不到幅值偏差大于δy的数据段；此时，估计的静态增益矢量表示为其中S1是当前步骤结束时得到幅度显著变化的集合，即步骤3.4：剩余的具有显著变化幅值的数据段集合为：重复步骤3.1到步骤3.3，直到剩余集合中数据段数目小于I，最后，获得静态增益与现有技术相比，本公开的有益效果是：通过可视化验证该方法的有效性，克服了使用系统识别方法难以验证静态增益估计的问题。附图说明构成本申请的一部分的说明书附图用来提供对本申请的进一步理解，本申请的示意性实施例及其说明用于解释本申请，并不构成对本申请的不当限定。图1为本专利技术的基于工业大数中斜坡响应的动态系统静态增益估计流程图；图2(a)-图2(d)为本专利技术具体实施示例中的采样数据样本图；图3(a)-图3(f)为本专利技术具体实施示例中的计算数据图。具体实施方式应该指出，以下详细说明都是示例性的，旨在对本申请提供进一步的说明。除非另有指明，本文使用的所有技术和科学术语具有与本申请所属
的普通技术人员通常理解的相同含义。需要注意的是，这里所使用的术语仅是为了描述具体实施方式，而非意图限制根据本申请的示例性实施方式。如在这里所使用的，除非上下文另外明确指出，否则单数形式也意图包括复数形式，此外，还应当理解的是，当在本说明书中使用术语“包含”和/或“包括”时，其指明存在特征、步骤、操作、器件、组件和/或它们的组合。专业术语解释：所述动态系统，是指状态随时间而变化的系统。所述静态增益，是指从一个稳态到新的稳态，系统的单位变化程度。所述斜坡响应，是指一个输入量的变化斜率从零跃增到某有限值引起的时间响应。所述输出时间序列，是指输出数值按其发生的时间先后顺序排列而成的数列。所述输入时间序列，是指输入变量值按其发生的时间先后顺序排列而成的数列。如图1所示，基于历史数据斜坡响应的动态系统静态增益估计方法，包括：步骤1：将输出y和输入ui的时间序列分割成短数据段，每个数据段用一条直线表示。我们采用分段线性表示方法描述输出y的分段。将时间序列分离为M个短数据段对于第m个数据段nm表示第m段中第一个数据样本，m∈[1，M]，nm+1-nm-1为第m段数据段中的样本总数；用线性回归模型来描述：y(n)＝am+bmn+e(n)，其中，am表示第m段数据段的初始值，bm表示第m段数据段的斜率，e(n)是干扰；假设e是具有零均值和方差为的白噪声；得到拟合线性方程：其中，和分别是y(n)、am和bm的估计值，用M表示数据段个数；M的估计值其中，L(M)为拟合误差损失函数：步骤2：从数据段中找出输入和输出同时处于直线上且振幅变化较大的斜坡响应。对于第m个数据段这段数据的振幅变化量：引入指示序列：其中，Ay,0为y显著幅度变化的阈值。构造一个确定系数Dm,当Dm的值越接近1，说明拟合度越高。Dm与Ay,m可用一种关系近似表示若给定Dm的下界D0＝0.8，可得到Ay,0的值。同理，计算输入的振幅变化量和指示序列得到ui显著幅度变化的阈值那么整体序列为：保留整体序列I0(n)＝1的数据段，当I0(n)＝1时，表示输入输出同时处于幅度变化较大且为直线的数据段上，所有这些显著的振幅变化的数据段组成本文档来自技高网...

【技术保护点】
1.基于历史数据斜坡响应的动态系统静态增益估计方法，其特征是，包括：采用分段线性表示方法将动态系统的输出时间序列分割成若干个输出时间段；同时，将动态系统的输入时间序列分割成若干个输入数据段；每个数据段均用一条直线表示；找出输入值和输出值同时均处于各自的直线数据段上且输入值和输出值的振幅变化均超过设定阈值的数据段，此数据段即为目标斜坡响应，将找出的所有输入数据段和所有输出数据段组成具有显著振幅变化的数据段集合；通过具有显著振幅变化的数据段集合，估计静态增益。

【技术特征摘要】
1.基于历史数据斜坡响应的动态系统静态增益估计方法，其特征是，包括：采用分段线性表示方法将动态系统的输出时间序列分割成若干个输出时间段；同时，将动态系统的输入时间序列分割成若干个输入数据段；每个数据段均用一条直线表示；找出输入值和输出值同时均处于各自的直线数据段上且输入值和输出值的振幅变化均超过设定阈值的数据段，此数据段即为目标斜坡响应，将找出的所有输入数据段和所有输出数据段组成具有显著振幅变化的数据段集合；通过具有显著振幅变化的数据段集合，估计静态增益。2.如权利要求1所述的方法，其特征是，所述将动态系统的输出时间序列分割成若干个输出时间段的具体步骤为：将时间序列分离为M个数据段对于第m个数据段nm表示第m段中第一个数据样本，m∈[1，M]，nm+1-nm-1为第m段数据段中的样本总数；用线性回归模型来描述：y(n)＝am+bmn+e(n)；其中，am表示第m段数据段的初始值，bm表示第m段数据段的斜率，e(n)是干扰；假设e(n)是具有零均值和方差为的白噪声；得到拟合线性方程：其中，和分别是y(n)、am和bm的估计值，用M表示数据段个数；M的估计值其中，L(M)为拟合误差损失函数：3.如权利要求1所述的方法，其特征是，将动态系统的输入时间序列分割成若干个输入数据段：采用分段线性表示方法将输入的时间序列分割为M个数据段时间序列的分段个数M，M的估计值：4.如权利要求1所述的方法，其特征是，找出输入值和输出值同时均处...

【专利技术属性】
技术研发人员：王建东，周东华，
申请(专利权)人：山东科技大学，
类型：发明
国别省市：山东,37