CN113433819A - 一种筛选pid控制回路中数据的方法和系统辨识方法 - Google Patents

一种筛选pid控制回路中数据的方法和系统辨识方法 Download PDF

Info

Publication number
CN113433819A
CN113433819A CN202110642209.8A CN202110642209A CN113433819A CN 113433819 A CN113433819 A CN 113433819A CN 202110642209 A CN202110642209 A CN 202110642209A CN 113433819 A CN113433819 A CN 113433819A
Authority
CN
China
Prior art keywords
data
loop
segment
value
excitation
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN202110642209.8A
Other languages
English (en)
Other versions
CN113433819B (zh
Inventor
刘志勇
吴庆尉
吴洁芸
裘坤
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Zhongkong Technology Co ltd
Original Assignee
Zhejiang Supcon Technology Co Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Zhejiang Supcon Technology Co Ltd filed Critical Zhejiang Supcon Technology Co Ltd
Priority to CN202110642209.8A priority Critical patent/CN113433819B/zh
Publication of CN113433819A publication Critical patent/CN113433819A/zh
Application granted granted Critical
Publication of CN113433819B publication Critical patent/CN113433819B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G05—CONTROLLING; REGULATING
    • G05B—CONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
    • G05B11/00—Automatic controllers
    • G05B11/01—Automatic controllers electric
    • G05B11/36—Automatic controllers electric with provision for obtaining particular characteristics, e.g. proportional, integral, differential
    • G05B11/42—Automatic controllers electric with provision for obtaining particular characteristics, e.g. proportional, integral, differential for obtaining a characteristic which is both proportional and time-dependent, e.g. P. I., P. I. D.

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Automation & Control Theory (AREA)
  • Testing And Monitoring For Control Systems (AREA)

Abstract

本发明提供一种筛选PID控制回路中数据的方法，包括：获取PID控制回路的历史运行数据；对历史运行数据进行预处理，获得预处理后的运行数据；根据PID控制回路的开闭环状态，对预处理后的运行数据进行分割，获得多个开环数据段和多个闭环数据段；基于开环数据激励段检测策略，根据每一个开环数据段中输入值的变化和输入值变化后引起的输出变化，从每一个开环数据段中筛选出激励数据段；基于闭环数据激励段检测策略，根据每一个闭环数据段中设定值的变化和设定值变化后引起的输出变化，从每一个闭环数据段中筛选出激励数据段。能够从PID控制回路的历史运行数据中自动筛选出充分激励的数据段。本发明还提供一种系统辨识方法和计算机设备。

Description

一种筛选PID控制回路中数据的方法和系统辨识方法
技术领域
本发明涉及工业自动化控制技术领域，尤其涉及一种筛选PID控制回路中数据的方法和系统辨识方法。
背景技术
在我国，实际工业现场还存在大量的回路采用的是人工控制。其中PID控制(比例积分微分控制)是一种很常见的回路控制，降低回路PID参数的整定难度，提供经济实惠的模型预测控制(MPC)方案，是提高我国过程工业自动化的重要方法，也是当前我国工业自动化公司的重要使命。系统辨识是MPC的关键步骤，同时也是PID参数整定的一种重要方法。
系统辨识要求数据是充分激励的。在以往的系统辨识中，常常需要进行激励实验，导致回路的正常运行受到干扰，降低了工厂的效益。而实际上，工厂储存了大量的回路运行历史数据，虽然历史数据中包含各种异常段、平稳段等，但是仍然存在能反映过程真实动态变化的数据段，而这些数据段可用于系统辨识。如果能够从历史数据中，自动筛选出充分激励的数据段，那么将显著降低模型预测控制、PID参数整定等的成本，提高工厂经济效益。
因此，亟需一种筛选PID控制回路中数据的方法和系统辨识方法，能够从回路历史运行数据中自动筛选出适合于系统辨识的数据段。
发明内容
(一)要解决的技术问题
鉴于上述技术中存在的问题，本发明至少从一定程度上进行解决。为此，本发明的一个目的在于提出了一种筛选PID控制回路中数据的方法，能够从PID控制回路的历史运行数据中自动筛选出充分激励的数据段。
本发明的第二个目的在于提出一种系统辨识方法，使辨识获得的模型精度高，数值计算可靠。
本发明的第三个目的在于提出一种计算机设备。
(二)技术方案
为了达到上述目的，本发明一方面提供一种筛选PID控制回路中数据的方法，包括：
获取PID控制回路的历史运行数据；运行数据包括每一时刻PID控制回路的输入，输出，设定值，开闭环状态。
对历史运行数据进行预处理，获得预处理后的运行数据。
根据PID控制回路的开闭环状态，对预处理后的运行数据进行分割，获得多个开环数据段和多个闭环数据段。
基于开环数据激励段检测策略，根据每一个开环数据段中输入值的变化和输入值变化后引起的输出变化，从每一个开环数据段中筛选出激励数据段；基于闭环数据激励段检测策略，根据每一个闭环数据段中设定值的变化和设定值变化后引起的输出变化，从每一个闭环数据段中筛选出激励数据段。
可选地，基于开环数据激励段检测策略，根据每一个开环数据段中输入值的变化和输入值变化后引起的输出变化，从每一个开环数据段中筛选出激励段数据，包括：根据每一个开环数据段中输入值的变化梯度Δu(t)，筛选出满足|Δu(t)|＞cuRu的目标输入值变化梯度；根据稳态段判断标准检测每一个开环数据段中发生目标输入值变化梯度后引起的输出变化，若所述输出变化中存在稳态段，则将所述输出变化到达稳态段之前对应的运行数据作为激励段数据，若所述输出变化中不存在稳态段，则将所述输出变化对应的运行数据作为激励段数据。
基于闭环数据激励段检测策略，根据每一个闭环数据段中设定值的变化和设定值变化后引起的输出变化，从每一个闭环数据段中筛选激励段数据，包括：根据每一个闭环数据段中设定值的变化梯度Δr(t)，筛选出满足|Δr(t)|＞crRy的目标设定值变化梯度；根据稳态段判断标准检测每一个闭环数据段中发生目标设定值变化梯度后引起的输出变化，若输出变化中存在稳态段，则将输出变化到达稳态段之前对应的运行数据作为激励段数据，若输出变化中不存在稳态段，则将输出变化对应的运行数据作为激励段数据。
其中，稳态段判断标准为|Δy(t)|＜cyRy，Δy(t)为输出值的变化梯度，cy为第四阈值，Ry为r时刻输出值y(t)的量程；cu为第三阈值，Ru为t时刻输入值u(t)的量程；cr为第六阈值。
可选地，对历史运行数据进行预处理，包括：对历史运行数据中的缺失数据进行剔除处理；对历史运行数据中的异常峰值进行处理。
可选地，对历史运行数据中的缺失数据进行剔除处理，包括：基于数据剔除策略S(t)＝Su(t)∩Sy(t)∩Sr(t)，检测历史运行数据中每一时刻的运行数据，剔除S(t)＝0的运行数据；根据历史运行数据剔除处理后形成的数据段，生成多个有效数据段。
其中，Su(t)为t时刻输入值u(t)的状态码，如果u(t)缺失，则Su(t)＝0，否则Su(t)＝1；Sy(t)为t时刻输出值y(t)的状态码，如果y(t)缺失，则Sy(t)＝0，否则Sy(t)＝1；Sr(t)为r时刻设定值r(r)的状态码，如果r(t)缺失，则Sr(t)＝0，否则Sr(t)＝1；S(t)为t时刻运行数据的状态码，符号∩表示且运算，只有当Su(t)、Sy(t)、Sr(t)都等于1时，S(t)＝1，否则S(t)＝0。
可选地，根据历史运行数据剔除处理后形成的数据段，生成多个有效数据段，包括：根据历史运行数据剔除处理后形成的数据段，将间距小于第一阈值的相邻数据段间的缺失数据进行线性插值，形成新的数据段；检测每一个数据段的长度，剔除长度小于第二阈值的数据段，生成多个有效数据段。
可选地，对历史运行数据中的异常峰值进行处理，包括：对每一个有效数据段进行异常峰值初处理，获得每一个初处理数据段。
其中异常峰值初处理包括：计算每一个有效数据段中u(t)的平均值
Figure BDA0003107989160000031
和标准差σu，以及y(t)的平均值
Figure BDA0003107989160000041
和标准差σy；根据
Figure BDA0003107989160000042
筛选出每一个有效数据段中u(t)的异常峰值；根据
Figure BDA0003107989160000043
筛选出每一个有效数据段中y(t)的异常峰值；根据每一个有效数据段中异常峰值邻近的数据对该异常峰值进行插值处理。
对每一个初处理数据段进行异常峰值再处理，获得每一个再处理数据段。
其中异常峰值再处理包括：采用半窗口大小为p的中值滤波算法对每个初处理数据段进行滤波，获得每一个滤波数据段；根据每个初处理数据段及其对应的滤波数据段，计算两者之间的u(t)差值和y(t)差值，获得每一个差值Δu(t)数据段和每一个差值Δy(t)数据段；计算每一个差值Δu(t)数据段的平均值
Figure BDA0003107989160000044
和标准差σΔu，以及计算每一个差值Δy(r)数据段的平均值
Figure BDA0003107989160000045
和标准差σΔy；根据
Figure BDA0003107989160000046
筛选出每一个差值Δu(t)数据段中Δu(t)的异常峰值；根据
Figure BDA0003107989160000047
筛选出每一个差值Δy(t)数据段的异常峰值；根据每一个差值Δu(t)数据段中异常峰值邻近的数据对该异常峰值进行插值处理，根据每一个差值Δy(t)数据段中异常峰值邻近的数据对该异常峰值进行插值处理。
本发明第二方面提供一种系统辨识方法，包括：
采用如上所述的方法获取激励数据段。
采用高阶ARX模型对激励数据段进行辨识，获得高阶过程模型。
采用MORSM方法对高阶过程模型进行降阶处理，获得第一低阶过程模型。
根据阻尼高斯牛顿法和第一低阶过程模型求解高阶过程模型的渐近负对数似然函数最小化时对应的低阶过程模型，获得第二低阶过程模型，将第二低阶过程模型作为PID回路模型。
可选地，高阶辨识模型的渐近负对数似然函数为：
Figure BDA0003107989160000048
其中，m为考虑的频率离散点数；ns为激励数据段的总数；
Figure BDA0003107989160000051
Figure BDA0003107989160000052
Φu(ω)、Φv(ω)分别为输入u(t)、扰动项H(q)e(t)的自谱，Φue(ω)为输入u(t)与白噪声e(t)的互谱，R表示ARX模型的方程误差，即A(q)y(t)-B(q)u(t)的估计方差；
Figure BDA0003107989160000053
是高阶过程模型的频率响应，n表示高阶过程模型的阶数，符号^表示估计值，
Figure BDA0003107989160000054
上标l是要求解的第二低阶过程模型；j表示虚数。
在高阶辨识模型是对单激励数据段的辨识时，渐近负对数似然函数表示为：
Figure BDA0003107989160000055
根据阻尼高斯牛顿法和所述第一低阶过程模型求解高阶过程模型的渐近负对数似然函数最小化时对应的低阶过程模型，获得第二低阶过程模型，包括：
ε(ω)的雅可比矩阵J为：
Figure BDA0003107989160000056
其中，ε＝[ε(ω1)，…，ε(ωm)]T；θ为第一低阶过程模型，
Figure BDA0003107989160000057
Figure BDA0003107989160000058
下标nl表示低阶过程模型的阶数。
由于
Figure BDA0003107989160000059
则
Figure BDA00031079891600000510
其中，
Figure BDA0003107989160000061
令
Figure BDA0003107989160000062
Figure BDA0003107989160000063
ε(ω)的偏导表示为：
Figure BDA0003107989160000064
Figure BDA0003107989160000065
其中，
Figure BDA0003107989160000066
Figure BDA0003107989160000067
Figure BDA0003107989160000068
Figure BDA0003107989160000069
根据JTJ·Δθ＝-JTε和θ(k+1)＝θ(k)+α·Δθ，获得第二低阶过程模型；其中α为阻尼因子，初始值θ(0)是第一低阶过程模型。
本发明第三方面提出一种计算机设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的PID控制回路中的数据筛选程序，处理器执行PID控制回路中的数据筛选程序时，实现如上所述的筛选PID控制回路中数据的方法。
本发明第四方面提出一种计算机设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的系统辨识程序，处理器执行系统辨识程序时，实现如上所述的系统辨识方法。
(三)有益效果
本发明的有益效果是：
1、通过预先定义的运行数据状态码和基于状态码构建的数据剔除策略来检测历史运行数据并剔除缺失数据，简化了缺失数据剔除处理的流程。以及通过预先定义的开环状态码和闭环状态码用于激励数据段的筛选，简化了激励数据段筛选的流程。
2、基于3σ原则，通过检测有效数据段输入与输入平均值的偏差，和输出与输出平均值的偏差来确定异常峰值，以及结合检测数据段滤波前后输入差值与输入差值平均值的偏差，和输出差值和输出差值平均值的偏差来确定异常峰值；并根据插值法对异常峰值进行处理。能够全面地检测出异常峰值，使异常峰值处理的结果更好。
3、本发明提出的筛选PID控制回路中数据的方法，实现了通过检测开环数据中u(t)的变化，来筛选开环数据中的激励数据段，通过检测闭环数据中r(r)的变化，来筛选闭环数据中的激励数据段，能够从PID控制回路的历史运行数据中自动筛选出充分激励的数据段。
4、本发明提出的系统辨识方法，可以同时辨识多个激励数据段，适用于开环、闭环PID控制过程的辨识，普适性高，并且使辨识获得的模型精度高，数值计算可靠。
附图说明
本发明借助于以下附图进行描述：
图1为根据本发明一个实施例的PID闭环控制回路示意图；
图2为根据本发明一个实施例的PID开环控制回路示意图；
图3为根据本发明一个实施例的筛选PID控制回路中数据的方法的流程示意图；
图4为根据本发明一个实施例的系统辨识方法的流程示意图。
具体实施方式
为了更好的解释本发明，以便于理解，下面结合附图，通过具体实施方式，对本发明作详细描述。
PID控制在工业现场应用十分普遍。图1为本发明一个实施例示出的PID闭环控制回路示意图，图2为本发明一个实施例示出的PID开环控制回路示意图。图中，G(q)为被控对象的传递函数模型，H(q)为干扰噪声的传递函数模型，Gc(q)为PID控制器的传递函数，y(t)为t时刻被控对象的输出值，u(t)为t时刻被控对象的输入值，r(t)为t时刻的设定值，e(t)为高斯白噪声，均值为0，方差为λ，v(t)为干扰噪声。使用线性函数来描述被控对象，如下：
y(t)＝G(q)u(r)+v(r)
v(t)＝H(q)e(r)
式中，G(q)为被控对象模型，t为时间，q为延迟算子，有q-1u(t)＝u(t-1)。
被控对象模型G(q)的准确辨识要求被控对象的输入和输出数据是充分激励的。为获取充分激励的数据，在以往的系统辨识中，常常需要进行激励实验。而实际上，工厂的PID回路存在大量的历史运行数据，如果能从中筛选出可用的充分激励数据段，则可以减少激励实验的时间，甚至完全避免做激励实验。
在PID控制回路中主要有两类情况可能包含充分激励段：
a)回路运行于如图2所示的开环状态下，输入u(t)是人为给定的，只要u(t)的变化足以激励过程，回路的数据就可用于系统辨识。
b)回路运行于如图1所示的闭环状态下，只要设定值r(t)的变化足以激励过程，回路的数据就可用于系统辨识。
由此，本发明将历史运行数据中充分激励数据段的筛选问题转化为开环状态下u(t)的变化检测问题和闭环状态下r(t)的变化检测问题。
为此，本发明实施例提出一种筛选PID控制回路中数据的方法，通过获取包括PID控制回路的输入，输出，设定值，开闭环状态的历史运行数据，根据PID控制回路的开闭环状态，对历史运行数据进行分割，获得多个开环数据段和多个闭环数据段；并基于开环数据激励段检测策略，根据每一个开环数据段中输入值的变化和输入值变化后引起的输出变化，从每一个开环数据段中筛选出激励数据段，基于闭环数据激励段检测策略，根据每一个闭环数据段中设定值的变化和设定值变化后引起的输出变化，从每一个闭环数据段中筛选出激励数据段。实现了通过检测开环数据中u(t)的变化，来筛选开环数据中的激励数据段，通过检测闭环数据中r(t)的变化，来筛选闭环数据中的激励数据段，能够从PID控制回路的历史运行数据中自动筛选出充分激励的数据段。
为了更好的理解上述技术方案，下面将参照附图更详细地描述本发明的示例性实施例。虽然附图中显示了本发明的示例性实施例，然而应当理解，可以以各种形式实现本发明而不应被这里阐述的实施例所限制。相反，提供这些实施例是为了能够更清楚、透彻地理解本发明，并且能够将本发明的范围完整的传达给本领域的技术人员。
图3为本发明一个实施例的筛选PID控制回路中数据的方法的流程示意图。
如图3所示，该筛选PID控制回路中数据的方法包括以下步骤：
步骤101、获取PID控制回路的历史运行数据。
具体地，运行数据包括每一时刻PID控制回路的输入，输出，设定值，开闭环状态。
具体的，作为一个实施例，步骤101还包括：检测历史运行数据，如果运行数据t时刻缺失开闭环状态数据，则根据r时刻PID控制回路的输出值和设定值，确定t时刻的开闭环状态数据。这是由于在t时刻，开环状态下设定值和输出值是相等的，闭环状态下设定值和输出值是不同的，且设定值呈现阶梯变化。
步骤102、对历史运行数据进行预处理，获得预处理后的运行数据。
PID控制回路运行不可避免会存在各种类型的异常值，如停机导致的数据缺失，仪表失灵或上下游回路干扰导致的异常峰值等。在系统辨识前，需要处理这些异常值，因此需要对历史运行数据进行预处理。
具体地，作为一个实施例，对历史运行数据进行预处理，包括：对历史运行数据中的缺失数据进行剔除处理；对历史运行数据中的异常峰值进行处理。
优选地，作为一个实施例，对历史运行数据中的缺失数据进行剔除处理，包括：基于数据剔除策略S(t)＝Su(t)∩Sy(t)∩Sr(t)，检测历史运行数据中每一时刻的运行数据，剔除S(t)＝0的运行数据；根据历史运行数据剔除处理后形成的数据段，生成多个有效数据段。
其中，Su(t)为t时刻输入值u(t)的状态码，如果u(t)缺失，则Su(t)＝0，否则Su(t)＝1；Sy(t)为t时刻输出值y(t)的状态码，如果y(t)缺失，则Sy(r)＝0，否则Sy(t)＝1；Sr(t)为t时刻设定值r(t)的状态码，如果r(t)缺失，则Sr(t)＝0，否则Sr(t)＝1；S(t)为t时刻运行数据的状态码，符号∩表示且运算，只有当Su(t)、Sy(t)、Sr(t)都等于1时，S(t)＝1，否则S(t)＝0。
通过预先定义的运行数据状态码和基于状态码构建的数据剔除策略来检测历史运行数据并剔除缺失数据，简化了缺失数据剔除处理的流程。
进一步地，作为一个实施例，根据历史运行数据剔除处理后形成的数据段，生成多个有效数据段，包括：根据历史运行数据剔除处理后形成的数据段，将间距小于第四阈值的相邻数据段间的缺失数据进行线性插值，形成新的数据段；检测每一个数据段的长度，剔除长度小于第五阈值的数据段，生成多个有效数据段。避免数据段碎片化，同时舍弃很短的数据段。
具体地，检测历史运行数据剔除处理后形成的数据段的位置和长度，第i个数据段表示为：
[xstart，i，xend，i]
第i个数据段的长度表示为：
Li＝xend，i-xstart，i+1
式中，xstart，i和xend，i分别表示第i个数据段在历史运行数据中的起始位置和终止位置。
根据数据段的位置，筛选出间距小于第一阈值的相邻数据段，公式表示为：
Di＝xstart，i+1-xend，i
Di＜d1
式中，Di为相邻数据段的间距，d1为第一阈值。
将间距小于第一阈值的相邻数据段间的缺失数据区间(xend，i，xstart，i+1)进行线性插值，并将该缺失数据区间的全部S(t)设置为1，形成新的数据段。
根据数据段的长度，剔除长度小于第二阈值的数据段，生成多个有效数据段。即当Li＜L1时，L1为第二阈值，将该数据段的S(t)设置为0，也就是说舍弃很短的数据段。
优选地，作为一个实施例，对历史运行数据中的异常峰值进行处理，包括：
步骤102-1、对每一个有效数据段进行异常峰值初处理，获得每一个初处理数据段。
其中异常峰值初处理包括：计算每一个有效数据段中u(r)的平均值
Figure BDA0003107989160000111
和标准差σu，以及y(t)的平均值
Figure BDA0003107989160000112
和标准差σy；根据
Figure BDA0003107989160000113
筛选出每一个有效数据段中u(t)的异常峰值；根据
Figure BDA0003107989160000114
筛选出每一个有效数据段中y(t)的异常峰值；根据每一个有效数据段中异常峰值邻近的数据对该异常峰值进行插值处理。
步骤102-2、对每一个初处理数据段进行异常峰值再处理，获得每一个再处理数据段。
其中异常峰值再处理包括：采用半窗口大小为p的中值滤波算法对每个初处理数据段进行滤波，获得每一个滤波数据段。根据每个初处理数据段及其对应的滤波数据段，计算两者之间的u(t)差值和y(t)差值，获得每一个差值Δu(t)数据段和每一个差值Δy(t)数据段。计算每一个差值Δu(t)数据段的平均值
Figure BDA0003107989160000115
和标准差σΔu，以及计算每一个差值Δy(t)数据段的平均值
Figure BDA0003107989160000116
和标准差σΔy；根据
Figure BDA0003107989160000117
筛选出每一个差值Δu(t)数据段中Δu(t)的异常峰值；根据
Figure BDA0003107989160000121
筛选出每一个差值Δy(t)数据段的异常峰值。根据每一个差值Δu(t)数据段中异常峰值邻近的数据对该异常峰值进行插值处理，根据每一个差值Δy(t)数据段中异常峰值邻近的数据对该异常峰值进行插值处理。
基于3σ原则，通过检测有效数据段输入与输入平均值的偏差，和输出与输出平均值的偏差来确定异常峰值，以及结合检测数据段滤波前后输入差值与输入差值平均值的偏差，和输出差值和输出差值平均值的偏差来确定异常峰值；并根据插值法对异常峰值进行处理。能够全面地检测出异常峰值，使异常峰值处理的结果更好。
步骤103、根据PID控制回路的开闭环状态，对预处理后的运行数据进行分割，获得多个开环数据段和多个闭环数据段。
回路运行的历史数据包含了开环和闭环状态，而判断开环数据和闭环数据充分激励段的标准是不同的，故需要根据PID控制回路的开闭环状态，将预处理后的运行数据中的开环段和闭环段分割出来。
具体地，作为一个实施例，开环数据段为t时刻开环状态码Sopen(t)均为1的数据段，闭环数据段为t时刻闭环状态码Sclosed(t)均为1的数据段。
步骤104、基于开环数据激励段检测策略，根据每一个开环数据段中输入值的变化和输入值变化后引起的输出变化，从每一个开环数据段中筛选出激励数据段；基于闭环数据激励段检测策略，根据每一个闭环数据段中设定值的变化和设定值变化后引起的输出变化，从每一个闭环数据段中筛选出激励数据段。
优选地，作为一个实施例，基于开环数据激励段检测策略，根据每一个开环数据段中输入值的变化和输入值变化后引起的输出变化，从每一个开环数据段中筛选出激励数据段，包括：根据每一个开环数据段中输入值的变化梯度Δu(t)，筛选出满足|Δu(t)|＞cuRu的目标输入值变化梯度；根据稳态段判断标准检测每一个开环数据段中发生目标输入值变化梯度后引起的输出变化，若输出变化中存在稳态段，则将输出变化到达稳态段之前对应的运行数据作为激励段数据，若输出变化中不存在稳态段，则将输出变化对应的运行数据作为激励段数据。
其中，cu为第三阈值，要求输入的改变引起的输出变化要明显大于扰动和噪声引起的变化，可选取5％；Ru为t时刻输入值u(t)的量程；稳态段判断标准为|Δy(t)|＜cyRy，Δy(t)为输出值的变化梯度，cy为第四阈值，Ry为t时刻输出值y(t)的量程。
进一步地，作为一个实施例，根据开环数据段中筛选出的激励数据段，将间距小于第五阈值且之间不存在闭环数据的相邻激励数据段进行合并。避免激励数据段碎片化。
优选地，作为一个实施例，基于闭环数据激励段检测策略，根据每一个闭环数据段中设定值的变化和设定值变化后引起的输出变化，从每一个闭环数据段中筛选激励段数据，包括：根据每一个闭环数据段中设定值的变化梯度Δr(t)，筛选出满足|Δr(t)|＞crRy的目标设定值变化梯度；根据稳态段判断标准检测每一个闭环数据段中发生目标设定值变化梯度后引起的输出变化，若输出变化中存在稳态段，则将输出变化到达稳态段之前对应的运行数据作为激励段数据，若输出变化中不存在稳态段，则将输出变化对应的运行数据作为激励段数据。
其中，cr为第六阈值，要求设定值的改变引起的输出变化要明显大于扰动噪声引起的变化，可选取5％；稳态段判断标准为|Δy(t)|＜cyRy。
进一步地，作为一个实施例，根据闭环数据段中筛选出的激励数据段，将间距小于第五阈值且之间不存在开环数据的相邻数据段进行合并。
在激励数据段的筛选中，通过梯度来判断激励段，且阈值的选取具有很大的普适应，方法简单有效。
在以往的系统辨识方法中，不能对多个激励数据段进行辨识，并且还存在仅能对开环PID控制过程进行辨识，不能对闭环PID控制过程进行辨识的方法。为对上述筛选PID控制回路中数据方法获取的激励数据段进行系统辨识，本发明实施例还提供一种系统辨识方法。
下面就参照附图来描述根据本发明实施例提出的系统辨识方法。
图4为本发明一个实施例的系统辨识方法的流程示意图。
如图4所示，该系统辨识方法包括以下步骤：
步骤201、获取激励数据段。
步骤202、采用高阶ARX(Auto Regressive with eXogeneous，带外部输出的自回归)模型对激励数据段进行辨识，获得高阶过程模型。
具体地，作为一个实施例，高阶ARX模型结构为：
A(q)y(t)＝B(q)u(t)+e(t)
其中，A(q)＝1+a1q-1+a2q-2+…+anq-n，B(q)＝b1q-1+b2q-2+…+bnq-n，a1…an和b1…bn是要求解的模型的参数，n表示高阶ARX模型的阶数。
根据高阶ARX模型，辨识出的高阶过程模型
Figure BDA0003107989160000141
扰动模型
Figure BDA0003107989160000142
和扰动信号谱
Figure BDA0003107989160000143
表示为：
Figure BDA0003107989160000144
其中，上标n表示高阶ARX模型的阶数；下标N表示激励数据段；符号^表示估计值；ω表示频率；R表示ARX模型的方程误差，即A(q)y(t)-B(q)u(t)的估计方差；j表示虚数。
步骤203、采用MORSM(Model Order Reduction Steiglitz-McBride，斯蒂格里茨·麦克布赖德模型降阶)方法对高阶过程模型进行降阶处理，获得第一低阶过程模型。
具体地，获得的第一低阶过程模型表示为：
Figure BDA0003107989160000145
式中，q是算符，nl表示低阶过程模型的阶数。
步骤204、根据阻尼高斯牛顿法和第一低阶过程模型求解高阶过程模型的渐近负对数似然函数最小化时对应的低阶过程模型，获得第二低阶过程模型，将第二低阶过程模型作为PID回路模型。
根据极大似然准则，当N→∞时，高阶过程模型的渐近负对数似然函数为：
Figure BDA0003107989160000151
式中，
Figure BDA0003107989160000152
是高阶过程模型的频率响应；j表示虚数；上标l是要求解的第二低阶过程模型；Φu(ω)、Φv(ω)分别为输入u(t)、扰动项H(q)e(t)的自谱；Φue(ω)为输入u(t)与白噪声e(t)的互谱；R表示ARX模型的方程误差，即A(q)y(t)-B(q)u(t)的估计方差。
在高阶辨识模型是对单激励数据段的辨识时，令dω≈2πΔf，Δf为常数，单位Hz，为分析方便，取Δf＝1，渐近负对数似然函数可离散化为：
Figure BDA0003107989160000153
令
Figure BDA0003107989160000154
Figure BDA0003107989160000155
将ε(ωp)带入离散渐近负对数似然函数，得
Figure BDA0003107989160000156
式中，m为考虑的频率离散点数。
上式为平方准则，故可根据阻尼高斯牛顿法和第一低阶过程模型求解高阶过程模型的渐近负对数似然函数最小化时对应的低阶过程模型，包括：
ε(ω)的雅可比矩阵J为：
Figure BDA0003107989160000161
其中，ε＝[ε(ω1)，…，ε(ωm)]T；θ为第一低阶过程模型，
Figure BDA0003107989160000162
Figure BDA0003107989160000163
下标nl表示低阶过程模型的阶数。
由于
Figure BDA0003107989160000164
则
Figure BDA0003107989160000165
其中，
Figure BDA0003107989160000166
令
Figure BDA0003107989160000167
Figure BDA0003107989160000168
ε(ω)的偏导表示为：
Figure BDA0003107989160000169
Figure BDA00031079891600001610
其中，
Figure BDA00031079891600001611
Figure BDA00031079891600001612
Figure BDA00031079891600001613
Figure BDA00031079891600001614
根据JTJ·Δθ＝-JTε和θ(k+1)＝θ(k)+α·Δθ，获得第二低阶过程模型；其中α为阻尼因子，初始值θ(0)是第一低阶过程模型。
在高阶辨识模型是对多个激励数据段的辨识时，高阶辨识模型的渐近负对数似然函数变为：
Figure BDA0003107989160000171
式中，ns为激励数据段的总数。
求解过程与单激励数据段的辨识类似，此处不再赘述。
本发明实施例提出的系统辨识方法，可以同时辨识多个激励数据段，适用于开环、闭环PID控制过程的辨识，普适性高，并且使辨识获得的模型精度高，数值计算可靠。
本发明实施例提供一种计算机设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的PID控制回路中的数据筛选程序，处理器执行PID控制回路中的数据筛选程序时，实现如上所述的筛选PID控制回路中数据的方法。
本发明实施例提供一种计算机设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的系统辨识程序，处理器执行系统辨识程序时，实现如上所述的系统辨识方法。
需要理解的是，以上对本发明的具体实施例进行的描述只是为了说明本发明的技术路线和特点，其目的在于让本领域内的技术人员能够了解本发明的内容并据以实施，但本发明并不限于上述特定实施方式。凡是在本发明权利要求的范围内做出的各种变化或修饰，都应涵盖在本发明的保护范围内。

Claims (10)

1.一种筛选PID控制回路中数据的方法，其特征在于，包括：
获取PID控制回路的历史运行数据；所述运行数据包括每一时刻PID控制回路的输入，输出，设定值，开闭环状态；
对所述历史运行数据进行预处理，获得预处理后的运行数据；
根据PID控制回路的开闭环状态，对所述预处理后的运行数据进行分割，获得多个开环数据段和多个闭环数据段；
基于开环数据激励段检测策略，根据每一个开环数据段中输入值的变化和输入值变化后引起的输出变化，从每一个开环数据段中筛选出激励数据段；
基于闭环数据激励段检测策略，根据每一个闭环数据段中设定值的变化和设定值变化后引起的输出变化，从每一个闭环数据段中筛选出激励数据段。
2.根据权利要求1所述的筛选PID控制回路中数据的方法，其特征在于，所述基于开环数据激励段检测策略，根据每一个开环数据段中输入值的变化和输入值变化后引起的输出变化，从每一个开环数据段中筛选出激励段数据，包括：
根据每一个开环数据段中输入值的变化梯度Δu(r)，筛选出满足|Δu(t)|＞cuRu的目标输入值变化梯度；根据稳态段判断标准检测每一个开环数据段中发生目标输入值变化梯度后引起的输出变化，若所述输出变化中存在稳态段，则将所述输出变化到达稳态段之前对应的运行数据作为激励段数据，若所述输出变化中不存在稳态段，则将所述输出变化对应的运行数据作为激励段数据；
所述基于闭环数据激励段检测策略，根据每一个闭环数据段中设定值的变化和设定值变化后引起的输出变化，从每一个闭环数据段中筛选激励段数据，包括：
根据每一个闭环数据段中设定值的变化梯度Δr(t)，筛选出满足|Δr(t)|＞crRy的目标设定值变化梯度；根据稳态段判断标准检测每一个闭环数据段中发生目标设定值变化梯度后引起的输出变化，若所述输出变化中存在稳态段，则将所述输出变化到达稳态段之前对应的运行数据作为激励段数据，若所述输出变化中不存在稳态段，则将所述输出变化对应的运行数据作为激励段数据；
其中，稳态段判断标准为|Δy(t)|＜cyRy，Δy(t)为输出值的变化梯度，cy为第四阈值，Ry为r时刻输出值y(t)的量程；cu为第三阈值，Ru为t时刻输入值u(t)的量程；cr为第六阈值。
3.根据权利要求1所述的筛选PID控制回路中数据的方法，其特征在于，所述对历史运行数据进行预处理，包括：
对历史运行数据中的缺失数据进行剔除处理；
对历史运行数据中的异常峰值进行处理。
4.根据权利要求3所述的筛选PID控制回路中数据的方法，其特征在于，所述对历史运行数据中的缺失数据进行剔除处理，包括：
基于数据剔除策略S(t)＝Su(t)∩Sy(t)∩Sr(t)，检测历史运行数据中每一时刻的运行数据，剔除S(t)＝0的运行数据；
根据历史运行数据剔除处理后形成的数据段，生成多个有效数据段；
其中，Su(t)为t时刻输入值u(t)的状态码，如果u(t)缺失，则Su(t)＝0，否则Su(t)＝1；Sy(t)为t时刻输出值y(t)的状态码，如果y(t)缺失，则Sy(t)＝0，否则Sy(t)＝1；Sr(t)为t时刻设定值r(t)的状态码，如果r(r)缺失，则Sr(r)＝0，否则Sr(t)＝1；S(t)为t时刻运行数据的状态码，符号∩表示且运算，只有当Su(t)、Sy(t)、Sr(t)都等于1时，S(t)＝1，否则S(t)＝0。
5.根据权利要求4所述的筛选PID控制回路中数据的方法，其特征在于，所述根据历史运行数据剔除处理后形成的数据段，生成多个有效数据段，包括：
根据历史运行数据剔除处理后形成的数据段，将间距小于第一阈值的相邻数据段间的缺失数据进行线性插值，形成新的数据段；检测每一个数据段的长度，剔除长度小于第二阈值的数据段，生成多个有效数据段。
6.根据权利要求4或5所述的筛选PID控制回路中数据的方法，其特征在于，所述对历史运行数据中的异常峰值进行处理，包括：
对每一个有效数据段进行异常峰值初处理，获得每一个初处理数据段；
其中异常峰值初处理包括：计算每一个有效数据段中u(t)的平均值
Figure FDA0003107989150000031
和标准差σu，以及y(t)的平均值
Figure FDA0003107989150000032
和标准差σy；根据
Figure FDA0003107989150000033
筛选出每一个有效数据段中u(t)的异常峰值；根据
Figure FDA0003107989150000034
筛选出每一个有效数据段中y(t)的异常峰值；根据每一个有效数据段中异常峰值邻近的数据对该异常峰值进行插值处理；
对每一个初处理数据段进行异常峰值再处理，获得每一个再处理数据段；
其中异常峰值再处理包括：采用半窗口大小为p的中值滤波算法对每个初处理数据段进行滤波，获得每一个滤波数据段；根据每个初处理数据段及其对应的滤波数据段，计算两者之间的u(r)差值和y(t)差值，获得每一个差值Δu(t)数据段和每一个差值Δy(t)数据段；计算每一个差值Δu(t)数据段的平均值
Figure FDA0003107989150000035
和标准差σΔu，以及计算每一个差值Δy(t)数据段的平均值
Figure FDA0003107989150000036
和标准差σΔy；根据
Figure FDA0003107989150000037
筛选出每一个差值Δu(t)数据段中Δu(t)的异常峰值；根据
Figure FDA0003107989150000038
筛选出每一个差值Δy(t)数据段的异常峰值；根据每一个差值Δu(t)数据段中异常峰值邻近的数据对该异常峰值进行插值处理，根据每一个差值Δy(t)数据段中异常峰值邻近的数据对该异常峰值进行插值处理。
7.一种系统辨识方法，其特征在于，包括：
采用如权利要求1至6任一项所述的方法获取激励数据段；
采用高阶ARX模型对所述激励数据段进行辨识，获得高阶过程模型；
采用MORSM方法对所述高阶过程模型进行降阶处理，获得第一低阶过程模型；
根据阻尼高斯牛顿法和所述第一低阶过程模型求解所述高阶过程模型的渐近负对数似然函数最小化时对应的低阶过程模型，获得第二低阶过程模型，将所述第二低阶过程模型作为PID回路模型。
8.根据权利要求7所述的系统辨识方法，其特征在于，
所述高阶辨识模型的渐近负对数似然函数为：
Figure FDA0003107989150000041
其中，m为考虑的频率离散点数；ns为激励数据段的总数；
Figure FDA0003107989150000042
Figure FDA0003107989150000043
Φu(ω)、Φv(ω)分别为输入u(t)、扰动项H(q)e(t)的自谱，Φue(ω)为输入u(t)与白噪声e(t)的互谱，R表示ARX模型的方程误差，即A(q)y(t)-B(q)u(t)的估计方差；
Figure FDA0003107989150000044
是高阶过程模型的频率响应，n表示高阶过程模型的阶数，符号^表示估计值，
Figure FDA0003107989150000045
上标l是要求解的第二低阶过程模型；j表示虚数；
在高阶辨识模型是对单激励数据段的辨识时，渐近负对数似然函数表示为：
Figure FDA0003107989150000046
所述根据阻尼高斯牛顿法和所述第一低阶过程模型求解高阶过程模型的渐近负对数似然函数最小化时对应的低阶过程模型，获得第二低阶过程模型，包括：
ε(ω)的雅可比矩阵J为：
Figure FDA0003107989150000051
其中，ε＝[ε(ω1)，…，ε(ωm)]T；θ为第一低阶过程模型，
Figure FDA0003107989150000052
Figure FDA0003107989150000053
下标nl表示低阶过程模型的阶数。
由于
Figure FDA0003107989150000054
则
Figure FDA0003107989150000055
其中，
Figure FDA0003107989150000056
令
Figure FDA0003107989150000057
Figure FDA0003107989150000058
ε(ω)的偏导表示为：
Figure FDA0003107989150000059
Figure FDA00031079891500000510
其中，
Figure FDA00031079891500000511
Figure FDA00031079891500000512
Figure FDA00031079891500000513
Figure FDA00031079891500000514
根据JTJ·Δθ＝-JTε和θ(k+1)＝θ(k)+α·Δθ，获得第二低阶过程模型；其中α为阻尼因子，初始值θ(0)是第一低阶过程模型。
9.一种计算机设备，其特征在于，包括存储器、处理器及存储在存储器上并可在处理器上运行的PID控制回路中的数据筛选程序，所述处理器执行所述PID控制回路中的数据筛选程序时，实现如权利要求1-6中任一项所述的筛选PID控制回路中数据的方法。
10.一种计算机设备，其特征在于，包括存储器、处理器及存储在存储器上并可在处理器上运行的系统辨识程序，所述处理器执行所述系统辨识程序时，实现如权利要求7或8中所述的系统辨识方法。
CN202110642209.8A 2021-06-09 2021-06-09 一种系统辨识方法和计算机设备 Active CN113433819B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110642209.8A CN113433819B (zh) 2021-06-09 2021-06-09 一种系统辨识方法和计算机设备

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110642209.8A CN113433819B (zh) 2021-06-09 2021-06-09 一种系统辨识方法和计算机设备

Publications (2)

Publication Number Publication Date
CN113433819A true CN113433819A (zh) 2021-09-24
CN113433819B CN113433819B (zh) 2022-05-10

Family

ID=77755566

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110642209.8A Active CN113433819B (zh) 2021-06-09 2021-06-09 一种系统辨识方法和计算机设备

Country Status (1)

Country Link
CN (1) CN113433819B (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114488780A (zh) * 2022-03-21 2022-05-13 浙江中控技术股份有限公司 一种pid控制方法以及pid控制器

Citations (27)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPS62224809A (ja) * 1986-03-26 1987-10-02 Toshiba Corp 同定装置
CN101493677A (zh) * 2008-06-30 2009-07-29 李华嵩 一种神经网络的系统控制器结构及系统辨识结构
CN102385307A (zh) * 2011-09-26 2012-03-21 北京空间飞行器总体设计部 一种星地时间系统的钟差测量方法
CN103529698A (zh) * 2013-10-17 2014-01-22 广东电网公司电力科学研究院 发电机调速系统参数辨识方法
CN103916339A (zh) * 2014-01-19 2014-07-09 中国人民解放军后勤工程学院 梯度最速下降anf频率无偏估计方法
CN104614985A (zh) * 2014-11-27 2015-05-13 北京航空航天大学 一种基于非线性规划的高阶系统最优降阶方法
CN105373094A (zh) * 2015-10-22 2016-03-02 太极控制有限公司 基于数据快采样的工业过程动态模型辨识装置
CN105426341A (zh) * 2015-10-29 2016-03-23 中国科学院力学研究所 复杂对象的参数辨识方法和装置
CN105589334A (zh) * 2016-02-26 2016-05-18 北京控制工程研究所 一种基于输入平滑的特征模型参数辨识方法
CN105868165A (zh) * 2016-04-15 2016-08-17 华中科技大学 一种电站锅炉运行数据清洗方法
CN107168060A (zh) * 2017-05-31 2017-09-15 苏州博众精工科技有限公司 基于继电反馈技术的带弹簧伺服系统辨识方法
CN107370806A (zh) * 2017-07-12 2017-11-21 北京京东尚科信息技术有限公司 Http状态码监控方法、装置、存储介质和电子设备
JP2018063637A (ja) * 2016-10-14 2018-04-19 トヨタ自動車株式会社 パラメータ同定装置
CN108279567A (zh) * 2017-12-29 2018-07-13 浙江中控软件技术有限公司 用于鲁棒控制的系统辨识方法
CN108388108A (zh) * 2018-02-27 2018-08-10 浙江中控技术股份有限公司 一种多重冗余控制系统中同步数据的方法及装置
CN109863105A (zh) * 2016-09-13 2019-06-07 因温特奥股份公司 用于监控电梯设备的方法
CN109947088A (zh) * 2019-04-17 2019-06-28 北京天泽智云科技有限公司 基于模型全生命周期管理的设备故障预警系统
CN110083526A (zh) * 2019-03-15 2019-08-02 深圳壹账通智能科技有限公司 应用程序测试方法、装置、计算机装置及存储介质
EP3550390A1 (en) * 2016-12-02 2019-10-09 Kao Corporation Manufacturing line control system, manufacturing line control device, and manufacturing line control method
CN110389530A (zh) * 2019-07-18 2019-10-29 西北工业大学 基于数据筛选的mems陀螺仪参数辨识驱动控制方法
JP2019191836A (ja) * 2018-04-23 2019-10-31 東洋電機製造株式会社 異常検知システムおよびプログラム
CN110713087A (zh) * 2019-10-21 2020-01-21 北京猎户星空科技有限公司 一种电梯门状态检测方法及装置
CN111030758A (zh) * 2019-12-16 2020-04-17 哈尔滨工业大学（威海） 带有稀疏约束的自适应零吸引因子盲判决反馈均衡算法
CN111725840A (zh) * 2020-06-29 2020-09-29 浙江大学 一种直驱风电机组控制器参数辨识方法
CN112130542A (zh) * 2020-07-01 2020-12-25 浙江浙能台州第二发电有限责任公司 基于正常运行数据与系统辨识的控制回路性能评价方法
CN112861318A (zh) * 2021-01-12 2021-05-28 华中科技大学 一种基于切削成像的材料本构模型参数辨识方法及系统
US20210164676A1 (en) * 2017-12-22 2021-06-03 Watrix Technology Corporation Limited Systems and methods for controlling an air-conditioning system based on gait recognition

Patent Citations (27)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPS62224809A (ja) * 1986-03-26 1987-10-02 Toshiba Corp 同定装置
CN101493677A (zh) * 2008-06-30 2009-07-29 李华嵩 一种神经网络的系统控制器结构及系统辨识结构
CN102385307A (zh) * 2011-09-26 2012-03-21 北京空间飞行器总体设计部 一种星地时间系统的钟差测量方法
CN103529698A (zh) * 2013-10-17 2014-01-22 广东电网公司电力科学研究院 发电机调速系统参数辨识方法
CN103916339A (zh) * 2014-01-19 2014-07-09 中国人民解放军后勤工程学院 梯度最速下降anf频率无偏估计方法
CN104614985A (zh) * 2014-11-27 2015-05-13 北京航空航天大学 一种基于非线性规划的高阶系统最优降阶方法
CN105373094A (zh) * 2015-10-22 2016-03-02 太极控制有限公司 基于数据快采样的工业过程动态模型辨识装置
CN105426341A (zh) * 2015-10-29 2016-03-23 中国科学院力学研究所 复杂对象的参数辨识方法和装置
CN105589334A (zh) * 2016-02-26 2016-05-18 北京控制工程研究所 一种基于输入平滑的特征模型参数辨识方法
CN105868165A (zh) * 2016-04-15 2016-08-17 华中科技大学 一种电站锅炉运行数据清洗方法
CN109863105A (zh) * 2016-09-13 2019-06-07 因温特奥股份公司 用于监控电梯设备的方法
JP2018063637A (ja) * 2016-10-14 2018-04-19 トヨタ自動車株式会社 パラメータ同定装置
EP3550390A1 (en) * 2016-12-02 2019-10-09 Kao Corporation Manufacturing line control system, manufacturing line control device, and manufacturing line control method
CN107168060A (zh) * 2017-05-31 2017-09-15 苏州博众精工科技有限公司 基于继电反馈技术的带弹簧伺服系统辨识方法
CN107370806A (zh) * 2017-07-12 2017-11-21 北京京东尚科信息技术有限公司 Http状态码监控方法、装置、存储介质和电子设备
US20210164676A1 (en) * 2017-12-22 2021-06-03 Watrix Technology Corporation Limited Systems and methods for controlling an air-conditioning system based on gait recognition
CN108279567A (zh) * 2017-12-29 2018-07-13 浙江中控软件技术有限公司 用于鲁棒控制的系统辨识方法
CN108388108A (zh) * 2018-02-27 2018-08-10 浙江中控技术股份有限公司 一种多重冗余控制系统中同步数据的方法及装置
JP2019191836A (ja) * 2018-04-23 2019-10-31 東洋電機製造株式会社 異常検知システムおよびプログラム
CN110083526A (zh) * 2019-03-15 2019-08-02 深圳壹账通智能科技有限公司 应用程序测试方法、装置、计算机装置及存储介质
CN109947088A (zh) * 2019-04-17 2019-06-28 北京天泽智云科技有限公司 基于模型全生命周期管理的设备故障预警系统
CN110389530A (zh) * 2019-07-18 2019-10-29 西北工业大学 基于数据筛选的mems陀螺仪参数辨识驱动控制方法
CN110713087A (zh) * 2019-10-21 2020-01-21 北京猎户星空科技有限公司 一种电梯门状态检测方法及装置
CN111030758A (zh) * 2019-12-16 2020-04-17 哈尔滨工业大学（威海） 带有稀疏约束的自适应零吸引因子盲判决反馈均衡算法
CN111725840A (zh) * 2020-06-29 2020-09-29 浙江大学 一种直驱风电机组控制器参数辨识方法
CN112130542A (zh) * 2020-07-01 2020-12-25 浙江浙能台州第二发电有限责任公司 基于正常运行数据与系统辨识的控制回路性能评价方法
CN112861318A (zh) * 2021-01-12 2021-05-28 华中科技大学 一种基于切削成像的材料本构模型参数辨识方法及系统

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
孟杰: "电力系统阻尼控制设计中的降阶状态空间模型辨识", 《中国电力》 *
彭煊等: "改进的格型Steiglitz-McBride 算法*", 《电路与系统学报》 *
涂亚庆等: "自适应陷波器的科氏流量计信号频率跟踪方法", 《重庆大学学报》 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114488780A (zh) * 2022-03-21 2022-05-13 浙江中控技术股份有限公司 一种pid控制方法以及pid控制器
CN114488780B (zh) * 2022-03-21 2022-07-26 浙江中控技术股份有限公司 一种pid控制方法以及pid控制器

Also Published As

Publication number Publication date
CN113433819B (zh) 2022-05-10

Similar Documents

Publication Publication Date Title
US11604442B2 (en) Predictive monitoring and diagnostics systems and methods
JP5696171B2 (ja) 制御パラメータ調整方法、制御パラメータ調整システム及び制御パラメータ設定装置
Piovoso et al. Monitoring process performance in real-time
US20070112442A1 (en) Apparatus and method for analyzing model quality in a process control environment
SE529454C2 (sv) Förfarande och anordning för trimning och styrning
WO2008094833A1 (en) Apparatus and method for automated closed-loop identification of an industrial process in a process control system.
Shardt et al. Closed-loop identification with routine operating data: Effect of time delay and sampling time
CN113433819B (zh) 一种系统辨识方法和计算机设备
EP3296822B1 (en) Model-plant mismatch detection using model parameter data clustering for paper machines or other systems
GB2547137A (en) Rapid process model identification and generation
WO2013061126A1 (en) A method and a system for tuning multivariable pid controller
CN114861522A (zh) 一种基于人工智能元学习技术的精密制造质量监测方法和装置
Massimiliano Performance improvement of Smith predictor through automatic computation of dead time
Huang Performance assessment of processes with abrupt changes of disturbances
Xu et al. Auto-tuning of PID controller parameters with supervised receding horizon optimization
Aufderheide et al. A variably tuned multiple model predictive controller based on minimal process knowledge
Petersson et al. A comparison of two feedforward control structure assessment methods
CN111830833B (zh) 基于模范系统的对象辨识及控制系统优化整定方法
Kelly et al. Overcoming encoder quantisation noise in an adaptive position controller
Petruck et al. Online parameter estimation for model-based force control in milling processes
WO2022077180A1 (zh) 确定无味卡尔曼滤波器的模型参数的方法、装置和系统
CN117076260B (zh) 一种参数及设备异常检测方法及装置
EP3715963B1 (en) Dead time estimation device and test device including the same
Bobál et al. Self-tuning control of nonlinear servo system: Comparison of LQ and predictive approach
Tham et al. Application of on-line estimation techniques to fermentation processes

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant
CP01 Change in the name or title of a patent holder

Address after: 310053 No. 309 Liuhe Road, Binjiang District, Hangzhou City, Zhejiang Province

Patentee after: Zhongkong Technology Co.,Ltd.

Address before: 310053 No. 309 Liuhe Road, Binjiang District, Hangzhou City, Zhejiang Province

Patentee before: ZHEJIANG SUPCON TECHNOLOGY Co.,Ltd.

CP01 Change in the name or title of a patent holder