本发明属于数据驱动故障诊断和控制技术领域，具体涉及一种数据驱动的闭环系统稳定性监测方法。
背景技术
近年来，在计算机技术、电子、信息和通信技术快速发展的带动下，当今工业系统，例如化工生产、机械制造、能源系统等，在规模不断扩大的同时，集成化和复杂程度也变得越来越高。在复杂工业系统中，一个局部异常事件甚至都有可能导致整个工业系统性能下降或者导致重大的工业事故并造成巨大的经济损失。为了提高经济效益和保持行业竞争力，现代工业过程的安全性和可靠性成为了最关键的因素，并获得了学术界和工业领域的广泛关注。
在现有的基于模型的闭环反馈系统稳定性分析工具中，稳定裕度和间隙度量技术起到了至关重要的作用。但是由于这两项技术对系统模型的依赖，闭环反馈系统的稳定性往往只能离线进行分析，这使得系统故障对闭环系统稳定性的影响没有办法在线评估和监测。在现有的数据驱动故障诊断技术中：一方面，大多数方法仅仅考虑系统故障的监测和诊断问题，缺乏系统故障对闭环反馈系统稳定性影响的分析；另一方面，很少考虑和分析闭环数据中反馈控制器带来的数据间的耦合关系，严重影响了数据驱动故障诊断方法的实际应用效果。
技术实现要素：
本发明为解决现有技术缺乏对数据驱动的闭环系统稳定性监测的问题，提供了一种数据驱动的闭环系统稳定性监测方法。
本发明所述一种数据驱动的闭环系统稳定性监测方法，通过以下技术方案实现：
步骤一、采集过去某时刻闭环系统的闭环数据，包括闭环系统的输入信号、输出信号以及参考输入信号；
步骤二、利用采集的数据构造汉克尔矩阵；
步骤三、对所构造的汉克尔矩阵进行cholesky分解，构造系统的稳定象描述、系统的稳定核描述，并求取归一化的稳定象描述和归一化的稳定核描述；
步骤四、利用步骤三构造得到的稳定象描述计算稳定裕度；
步骤五、根据稳定裕度设置监测阈值，采集系统当前闭环数据，重复步骤二至步骤四，得到系统当前归一化的稳定象描述和系统当前归一化的稳定核描述；并计算系统当前的间隙度量；
步骤六、结合监测阈值和系统当前的间隙度量对系统进行实时监测。
作为对上述技术方案的进一步阐述：
进一步的，步骤二具体包括以下步骤：
步骤二一、利用反馈控制器k(z)＝(a,b,c,d)参数构造稳定滤波器
步骤二二、通过下述公式计算得到滤波后的参考输入信号w(z)：
其中，ω(z)为参考输入信号；
步骤二三、选取维度参数sp，sf和n，构造关于系统滤波后的参考输入信号w(z)、输入信号u(z)、输出信号y(z)的汉克尔矩阵：
其中，
进一步的，步骤三具体包括以下步骤：
步骤三一、对所构造的汉克尔矩阵做如下cholesky分解：
其中，
步骤三二、构建系统的稳定象描述
其中，
步骤三三、求取如下左零空间：
其中，
步骤三四、求取归一化的稳定象描述
其中，
进一步的，步骤四中所述稳定裕度的具体计算过程包括：
其中，
进一步的，步骤五具体包括以下步骤：
步骤五一、根据步骤四中所得到的稳定裕度
步骤五二、采集系统当前闭环数据，重复步骤二至步骤四得到系统当前归一化的稳定象描述
步骤五三、计算系统当前的间隙度量
进一步的，步骤六具体为：
将步骤五中得到的间隙度量
进一步的，步骤二一中所述稳定滤波器具体为：
其中，l表示观测器增益矩阵，其选取需要使得a-lc极点在单位圆之内。
本发明最为突出的特点和显著的有益效果是：
本发明所涉及的一种数据驱动的闭环系统稳定性监测方法，具有以下几方面优势：
1、本发明利用系统的输入输出闭环数据，即可对闭环系统的稳定性进行评估和监测，在闭环系统的稳定性遭到破坏时能够产生报警，能够监测到99％的系统故障。
2、本发明基于闭环过程数据，不依赖于系统模型，适用于绝大多数(90％以上)闭环控制系统。
附图说明
图1为本发明逻辑控制图；
图2为本发明方法流程图；
图3为实施例中闭环系统稳定性监测结果曲线图。
具体实施方式
具体实施方式一：结合图1、图2对本实施方式进行说明，本实施方式给出的一种数据驱动的闭环系统稳定性监测方法，具体包括以下步骤：
步骤一、采集过去某时刻闭环系统的闭环数据，包括闭环系统的输入信号、输出信号以及参考输入信号；
步骤二、利用采集的数据构造汉克尔矩阵；
步骤三、对所构造的汉克尔矩阵进行cholesky分解，构造系统的稳定象描述、系统的稳定核描述，并求取归一化的稳定象描述和归一化的稳定核描述；
步骤四、利用步骤三构造得到的稳定象描述计算稳定裕度；
步骤五、根据稳定裕度设置监测阈值，采集系统当前闭环数据，重复步骤二至步骤四，得到系统当前归一化的稳定象描述和系统当前归一化的稳定核描述；并计算系统当前的间隙度量；
步骤六、结合监测阈值和系统当前的间隙度量对系统进行实时监测。
具体实施方式二：本实施方式与具体实施方式一不同的是，步骤二具体包括以下步骤：
步骤二一、利用所述闭环系统中的反馈控制器k(z)＝(a,b,c,d)参数构造稳定滤波器
步骤二二、通过下述公式计算得到滤波后的参考输入信号w(z)：
其中，ω(z)为参考输入信号；
步骤二三、选取适当的维度参数sp，sf和n，构造关于系统滤波后的参考输入信号w(z)、输入信号u(z)、输出信号y(z)的汉克尔矩阵：
其中，
其他步骤及参数与具体实施方式一相同。
具体实施方式三：本实施方式与具体实施方式二不同的是，步骤三具体包括以下步骤：
步骤三一、对所构造的汉克尔矩阵做如下cholesky分解：
其中，
步骤三二、构建系统的稳定象描述
其中，
步骤三三、求取如下左零空间：
其中，
步骤三四、求取归一化的稳定象描述
其中，
其他步骤及参数与具体实施方式二相同。
具体实施方式四：本实施方式与具体实施方式三不同的是，步骤三具体包括：步骤四中所述稳定裕度的具体计算过程包括：
其中，
其他步骤及参数与具体实施方式三相同。
具体实施方式五：本实施方式与具体实施方式四不同的是，步骤五具体包括以下步骤：
步骤五一、根据步骤四中所得到的稳定裕度
步骤五二、采集系统当前(在线)闭环数据，重复步骤二至步骤四得到系统当前归一化的稳定象描述
步骤五三、计算系统当前的间隙度量
其他步骤及参数与具体实施方式四相同。
具体实施方式六：本实施方式与具体实施方式五不同的是，步骤六具体为：
将步骤五中得到的间隙度量
其他步骤及参数与具体实施方式五相同。
具体实施方式七：本实施方式与具体实施方式二、三、四、五或六不同的是，步骤二一中所述稳定滤波器具体为：
其中，l表示观测器增益矩阵，其选取需要使得a-lc极点在单位圆之内，也就是控制学里说的使得a-lc稳定。
其他步骤及参数与具体实施方式二、三、四、五或六相同。
实施例
利用以下仿真过程来验证本发明的有益效果：
选取正常系统g0(z)＝(a0,b0,c0,d0)为：
反馈控制器k(z)＝(a,b,c,d)为：
假设正常系统g0(z)将在运行过程中慢慢转变为故障系统gf(z)＝(af,bf,cf,df)：
本实施例所述一种数据驱动的闭环系统稳定性监测方法按照以下步骤进行：
步骤一、采集过去某时刻闭环系统的闭环数据，包括闭环系统的输入信号、输出信号以及参考输入信号；
步骤二、利用采集的数据构造汉克尔矩阵；
步骤二一：选取
步骤二二、通过下述公式计算得到滤波后的参考输入信号w(z)：
步骤二三、选取维度参数sp＝20，sf＝60和n＝5000，构造关于系统滤波后的参考输入信号w(z)、输入信号u(z)、输出信号y(z)的汉克尔矩阵：
步骤三、对所构造的汉克尔矩阵进行cholesky分解，构造系统的稳定象描述、系统的稳定核描述，并求取归一化的稳定象描述和归一化的稳定核描述；
步骤四、利用步骤三构造得到的稳定象描述计算稳定裕度
步骤五、根据稳定裕度设置监测阈值jth＝0.12，采集系统当前闭环数据，重复步骤二至步骤四，得到系统当前归一化的稳定象描述和系统当前归一化的稳定核描述；并计算系统当前的间隙度量；
步骤六、结合监测阈值和系统当前的间隙度量对系统进行实时监测。
闭环系统的稳定性监测结果如图3所示。由图3可看出，本发明方法能够通过系统闭环数据有效地监测故障对闭环系统稳定性的影响。
本发明还可有其它多种实施例，在不背离本发明精神及其实质的情况下，本领域技术人员当可根据本发明作出各种相应的改变和变形，但这些相应的改变和变形都应属于本发明所附的权利要求的保护范围。