CN105116873B - 一种火电厂多自动调节回路评价诊断方法 - Google Patents

Abstract

一种对火电厂多自动调节回路评价诊断方法。该方法以回路为单位，采集实时的数据并进行数据存储，在此数据存储的基础上，周期性的针对每一个自动调节回路的控制性能的主要指标进行计算并对自动调节回路中传感器和执行器进行故障诊断。由于控制性能评价需要大量数据而实时性不高，故障诊断需要部分数据而实时性要求高，本方法采用了不同的周期对其进行了划分。通过对自动调节回路进行排列后按周期依次进行计算，能够使多个自动调节回路控制性能评价和故障诊断的同时进行。

Description

一种火电厂多自动调节回路评价诊断方法

技术领域

[0001] 本发明属于工业控制领域，具体涉及一种火力发电机组自动调节回路控制性能的评价诊断系统。

背景技术

[0002] 随着工业自动化水平的提高，工业控制系统的规模越来越大，一个大型的工业控制系统有成百套控制回路，过程控制技术面临的重要任务之一就是对控制系统的监视和维护。在控制系统的使用过程中，有很多原因会导致控制器性能下降，比如设备的老化、非线性的作用等，为了确保控制器能在设计的性能下工作，就需要对大量的过程数据进行分析来判断控制器在工作过程中性能如何。一个大的自动化系统控制回路众多，人工来分析数据并判断控制器性能的好坏工作量相当大，这种人工监视评价控制器的性能的方法难以适应当前的高度自动化水平。研究控制系统的性能评价方法，形成与当前自动化水平相适应的控制系统性能评价体系，对实际的控制系统进行监视，确保控制器在工作过程中性能良好有很重要的意义。

[0003] 火电机组中，机组热工性能取决于各个自动调节回路的性能。执行机构状态、传感器状态、以及控制器参数选取这三个环节是一个自动调节回路关键的要素。在工业控制领域，生产过程的控制回路不断增加，生产工艺对生产过程的控制要求日益提高，而我们往往只关注采用怎样的控制器，采取怎样的控制策略，对于控制系统的性能评价虽有很多有益的探讨和研究，但还存在着一些缺点和不足。

[0004] 通常是针对个别指定的自动调节回路，针对调节过程数据的变化过程，套用理论公式得到动态调节指标。但这种方法在实际中对于调节过程要求较高，需要严格的数据趋势。另外以往的评价方法没有与故障诊断相结合，没有对故障情况下进行排除。同时，通常的评价诊断方法只是针对个别自动调节回路进行，缺少在工业控制过程中的多控制回路综合应用。

[0005] 另外，本发明的评价诊断方法都是基于在线运行系统的数据分析结果，比通常的离线分析系统具有更重要的实用意义。

发明内容

[0006] 为了解决现有技术中存在的上述技术问题，本发明公开了一种火电厂多自动调节回路评价诊断方法。

[0007] 本发明具体采用以下技术方案。

[0008] —种火电厂多自动调节回路评价诊断方法，所述方法针对火电厂多个自动调节回路进行控制性能的评价并对自动调节回路中的传感器和执行器进行故障诊断；其特征在于，所述方法包括以下步骤：

[0009] 步骤一:通过火电厂热工录波分析装置采集火电厂多个自动调节回路数据，包括机组负荷、自动调节回路中被控对象的设定值、传感器采集的测量值以及执行器开度指令、 开度反馈值；

[0010] 步骤二：以30秒为周期，由火电厂热工录波分析装置对采集数据进行处理，形成comtrade录波文件，内容包含30秒内采集的数据信息，形成录波文件之后，将其发送至作为诊断装置的电脑；

[0011] 步骤三：在诊断装置中，以30秒为周期解析由火电厂热工录波分析装置生成的comtrade录波文件，得到火电厂多个自动调节回路的数据信息，并存入历史数据库；

[0012] 步骤四：将需要进行评价诊断的火电厂自动调节回路做排列，依据采集的自动调节回路的数据信息以周期形式依次对单个自动调节回路进行诊断、性能评价，其中对传感器、执行器的诊断周期为60秒，S卩60秒内，依次对所有自动调节回路中的传感器、执行器进行诊断;对所述多个自动调节回路的控制性能评价的周期为10分钟，即10分钟内，依次对所有自动调节回路的控制性能进行评价指标计算。

[0013] 本发明进一步优选包括以下方案：

[0014] 在步骤四中，对传感器、执行器故障诊断过程如下：

[0015] (1)以60秒为周期，从历史数据库中提取待诊断的传感器、执行器所属的自动调节回路的数据信息，包括60秒内待诊断的传感器采集的测量值历史记录、执行器的开度指令和开度反馈值历史记录，并形成CSV格式文本文件；

[0016] (2)根据生成的csv格式文本文件的内容得到60秒内待诊断的传感器采集的测量值，基于所述测量值实现对传感器的诊断；

[0017] 当传感器测量值超出了其高低限值范围，则认为传感器发生了测量值超限故障；

[0018] 计算传感器所采集的测量值的变化速率，将所述测量值的变化速率与设定的变化率阈值比较，如果超限，则判断传感器存在变化速率超限故障；

[0019] 两个传感器冗余设置时，先比较两个冗余设置的传感器在同一时刻的测量值偏差，若超过了偏差限值，则说明其中有一个发生了故障的，然后判断这两个冗余设置的传感器是否发生了测量值超限故障或变化率超限故障，如果两个冗余设置的传感器均没有发生测量值超限故障以及变化率超限故障，则认为这两个冗余传感器存在冗余故障；

[0020] 三冗余及三个以上冗余设置的传感器，比较每两个传感器间测量值的偏差，如果偏差超过了偏差限值，则认为此两个传感器中某个传感器发生了故障，再与另一个传感器的测量值相比较，判断偏差大者为发生故障的那个传感器。

[0021] 对于两个传感器冗余设置时，可以通过判断出冗余的传感器存在故障，但不能直接定位具体的故障传感器，对于三个余传感器冗余设置时，可以通过两两对比定位具体的故障传感器。

[0022] ⑶根据生成的csv文件内容，得到执行器开度指令和开度反馈历史记录，判断执行器是否存在卡死故障、增益故障和死区故障；

[0023] 当满足如下两个判据中的任一个，则判断执行器出现了卡死故障：

[0024] 判据1:在60秒内，执行器开度指令发生变化，而开度反馈不随执行器开度指令变化；

[0025] 判据2:当执行器开度指令不变化时，开度反馈也不变化，但是开度反馈和执行器开度指令之间存在偏差，偏差范围超过开度反馈量程的5%;统计60秒内执行器开度反馈的变化程度与执行器开度指令变化程度均值的比值，作为执行器增益的估计值;如果这个增 益值超过了预设的增益阈值，则认为执行器发生了增益故障，此阈值为1.05〜1.1;

[0026] 根据执行器在正常情况下的往返行程中X-Y特性图选择最大宽度为正常死区范 围。

[0027] 在步骤四中，对所述多个自动调节回路的控制性能评价过程如下：

[0028] (1)以10分钟为周期执行，对已排序的多个自动调节回路依次计算，首先要检查该控制回路是否存在执行器或者传感器故障，如果存在故障则将此回路跳过，不进行控制评价指标计算，进行下一个自动调节回路的计算过程，如无故障则进入下一个流程；

[0029] (2)将10分钟内自动调节回路中传感器测量值、执行器开度指令、自动调节回路中被控对象的设定值以及机组负荷的历史值写入CSV格式文件中；

[0030] (3)自动调节回路中被控对象的设定值发生变化，则认为该自动调节回路开始一个动态过程，否则，则认为该自动调节回路处于一个随机过程;如果自动调节回路处于动态过程，则计算该自动调节回路的确定性调节指标，包括调节时间、上升时间、超调量;如果自动调节回路处于随机过程，则进行随机性指标即最小方差指标计算。

[0031] 调节时间:单位阶跃响应首次达到稳态值的5% (或2%)的误差带并保持的时间。

[0032] 上升时间：单位阶跃响应从稳态值的10%上升到稳态值的90%，所需要的时间。对于系统单位阶跃响应单调变化的系统如此，而对有震荡的系统，可以被定义为从初始值到第一次到达稳态的时间。上升时间，反应的是系统的响应速度。

[0033] 超调量：当系统的单位阶跃有超调时，其单位阶跃响应的最大峰值减去稳态值的差与稳态值之比的百分数。

[0034] 最小方差指标:系统理论最小方差与实际方差的比值。

[0035] 本发明具有以下技术效果：

[0036] 1)可以针对火电厂多个自动调节回路进行评价诊断；

[0037] 2)实现了在线评价与诊断；

[0038] 3)在进行自动调节回路控制评价时，考虑传感器、执行器故障状态对其的影响，将故障诊断与控制评价相结合；

[0039] 4)可以对调节过程进行确定性指标评价和随机性指标评价。

附图说明

[0040] 图1为本发明火电厂多回路自动调节回路评价诊断方法实现过程示意图。

具体实施方式

[0041] 下面结合说明书附图对本发明的技术方案做进一步详细介绍。如附图1所示为本发明公开的火电厂多回路自动调节回路评价诊断方法实现过程方案示意图。

[0042] 附图1中的设备代表火电机组所配置的数据转发装置。

[0043] 附图1中的热工数据采集装置通过接收转发装置送来的数据，并形成comtrade格式录波文件。

[0044] 附图1中的实时库用来存储所有需要实时显示的数据信息，包括所有的诊断、评价 结果。

[0045] 附图1中的历史库用来存储所有需要存储历史的数据信息，包括诊断、评价过程中 需要的测量值信息、执行器开度指令信息、执行器开度反馈信息。

[0046] 附图1中的平台扩展应用APP功能接口主要提供实现本发明的方法的应用接口。

[0047] 下面对本发明的技术方法实现过程方案做进一步详细介绍。步骤1:通过火电厂热工录波分析装置通过通讯规约，完成从火电厂DCS系统中相关数据的采集。

[0048] 主要采集的数据包括参与进行评价指标计算和诊断的机组负荷、自动调节回路中被控对象的设定值、传感器采集的测量值以及执行器开度指令、开度反馈。

[0049] 步骤2:以30秒为周期，由火电厂热工录波分析装置对采集数据进行处理，形成comtrade录波文件，内容包含30秒内采集的数据信息，形成录波文件之后，将其发送至进行评价诊断的电脑。

[0050] 步骤3:以30秒为周期解析由火电厂热工录波分析装置生成的comtrade录波文件，得到数据信息，并存储历史。存储内容为各个数据点的时间、值。

[0051] 步骤1一一步骤3完成了数据的采集和存储工作。

[0052] 步骤4:将需要进行评价诊断的火电厂自动调节回路在文件中做排列，之后以周期形式依次对单个自动调节回路进行诊断、性能评价。其中传感器、执行器的诊断周期为60秒，S卩60秒内，依次对所有自动调节回路中的传感器、执行器进行诊断。控制性能评价的周期为10分钟，即10分钟内，依次对所有自动调节回路的控制性能进行评价指标计算。具体如下：

[0053] 在步骤4中，对传感器、执行器故障诊断过程如下：

[0054] (1)以60秒为周期，从历史数据库中提取待诊断的传感器、执行器所属的自动调节回路的数据信息，包括60秒内待诊断的传感器采集的测量值历史记录、执行器的开度指令和开度反馈值历史记录，并形成CSV格式文本文件；

[0055] (2)根据生成的csv格式文本文件的内容，得到60秒内待诊断的传感器采集的测量值，基于所述测量值实现对传感器的诊断；

[0056] 当传感器测量值超出了其高低限制范围，则认为传感器发生了测量值超限故障，此高低限制范围为传感器量程，具体的需要参考被诊断传感器的铭牌参数；

[0057] 计算传感器所采集的测量值的变化速率，将所述测量值的变化速率与设定的变化率阈值比较，如果超限，则判断传感器存在变化速率超限故障，所述的变化率阈值可以参考被诊断传感器的铭牌参数；

[0058] 两个传感器冗余设置时，先比较两个冗余设置的传感器在同一时刻的测量值偏差，若超过了偏差限值，则说明其中有一个发生了故障的，此偏差限制为传感器量程的5%，然后判断这两个冗余设置的传感器是否发生了测量值超限故障或变化率超限故障，如果两个冗余设置的传感器均没有发生测量值超限故障以及变化率超限故障，则认为这两个冗余传感器存在冗余故障。

[0059] 三冗余及三个以上冗余设置的传感器，比较每两个传感器间测量值的偏差，如果偏差超过了偏差限值，则认为此两个传感器中某个传感器发生了故障，再与另一个传感器的测量值相比较，判断偏差大者为发生故障的那个传感器，所述的偏差限制为传感器量程的5 %。

[0060] 对于两个传感器冗余设置时，可以通过判断出冗余的传感器存在故障，但不能直接定位具体的故障传感器，对于三个余传感器冗余设置时，可以通过两两对比定位具体的 故障传感器。

[0061] (3)根据生成的csv文件内容，得到执行器开度指令和开度反馈历史记录，判断执行器是否存在卡死故障、增益故障和死区故障；

[0062] 当满足如下两个判据中的任一个，则判断执行器出现了卡死故障：

[0063] 判据1:在60秒内，执行器开度指令发生变化，而开度反馈不随执行器开度指令变化；

[0064] 判据2:当执行器开度指令不变化时，开度反馈也不变化，但是开度反馈和执行器开度指令之间存在偏差，偏差范围超过开度反馈量程的5% ;

[0065] 统计60秒内执行器开度反馈的变化程度与执行器开度指令变化程度均值的比值，作为执行器增益的估计值;如果这个增益值超过了预设的增益阈值，则认为执行器发生了增益故障；

[0066] 根据执行器在正常情况下的往返行程中X-Y特性图选择最大宽度为正常死区范 围。

[0067] 在步骤4中，对所述多个自动调节回路的控制性能评价过程如下：

[0068] (1)以10分钟为周期执行，按步骤四中排列的自动调节回路依次计算，首先要检查该控制回路是否存在执行器或者传感器故障，如果存在故障则将此回路跳过，不进行控制评价指标计算，进行下一个自动调节回路的计算过程，如无故障则进入下一个流程；

[0069] (2)将10分钟内自动调节回路中传感器测量值、执行器开度指令、自动调节回路中被控对象的设定值以及机组负荷的历史值写入csv格式文件中；

[0070] (3)自动调节回路中被控对象的设定值发生变化，则认为该自动调节回路开始一个动态过程，否则，则认为该自动调节回路处于一个随机过程;如果自动调节回路处于动态过程，则计算该自动调节回路的确定性调节指标，包括调节时间、测量值上升时间、超调量；如果自动调节回路处于随机过程，则进行随机性指标即最小方差指标计算。

[0071] 计算自动调节回路的确定性调节指标的过程如下：

[0072] ①建立自动调节回路中被控对象的数学模型：

[0073] 自动调节回路的输入一输出模型的差分方程为：

[0074]

[0075] 其中：e⑴为噪声，沖)'为输入一输出观测向量，y⑴为自动调节回路的测量值，u ⑴为自动调节回路中的阀位反馈，Θ为未知参数向量，共有na，个输出、nb个输入，m……ana，bi......bnb为其对应参数，d为迟延参数；

[0076] 且

[0079] 令根据N次观测{y (i)，u(i) :i = l，2,…，N，N彡na+nb+l}对Θ的估计为表，引入符号：

[0080]

[0081] y⑴，y⑵一y (N)为连续N次采样得到的测量值；

[0082] 可以得到对参数的最小二乘估计为：

[0083]

[0084] 其中

为参数的最小二乘估计，得到的矩阵依次为ai、a2……ana的值；

[0085] ②根据前述的自动调节回路的数学模型数学模型，施加单位阶跃扰动信号，形成理想的输出曲线，根据调节曲线，得到每个时刻对应的输出值，在此基础上，计算得到上升时间、调节时间和超调量。

[0086] 计算自动调节回路的随机性调节指标的过程如下：

[0087] ①通过火电厂热工录波分析装置采集火电厂多个自动调节回路数据信息，包括机组负荷、自动调节回路中被控对象的设定值、传感器采集的测量值以及执行器开度指令、开度反馈值；

[0088] ②在采集的数据中，选择数据长度至少为n = 500且包含系统典型的扰动特性的一段数据进行评价，η即表示数据长度；

[0089] ③进行开环测试或用闭环估计技术估计采样数据的时间延迟k;

[0090] ④自动调节回路的输出y (t)表示为下式

[0091] y (t) =F| (t) +aiy (t_k) +a2y (t-k-l) + · · ·+amy (t-k-m+1)

[0092]

[0093] Q1……知为模型的参数；

[0094] k为迟延因子，由自动调节回路的迟延时间决定；

[0095] m为自回归模型的长度，m取5〜30,模型长度取决于采样时间，并且模型长度的选择要使最终模型能全部捕捉系统的脉冲响应；

[0096] Fl (t)为系统输出在t时刻的响应，其余项为扰动信号对系统输出的响应；

[0097] ⑤由④可以获得数据长度η范围内的自动调节回路中传感器采集的测量值，叠加写成以下形式

[0098] y=Xa+F|

[0099] 其中，

[0100] 由获得数据根据上式构建上述向量y和矩阵X;

[0101] y、X矩阵中值为对应下标时刻的自动调节回路中传感器采集的测量值；

[0102] ⑥选择合适的模型获得<^的参数；

[0103] 用线性回归方法估计自回归模型参数{αι}，最小二乘结果为：

[0104]

[0105] ⑦计算自动调节回路的最小方差，最小方差的估计值可以由残余方差代替：

[0106]

[0107]

:是理论的最小方差；

[0108] ⑧计算系统实际的输出方差：

[0109]

[0110]

是实际的输出方差，

是输出方差中出去理论最小方差的部分。

[0111] ⑨计算随机性调节指标：

[0112]

[0113] ;)(/c)是最小方差指标，由上式可知，其取值范围为0-1，此值越接近1，则表示实际 方差与理论最小方差越接近，说明自动调节回路抗随机扰动能力更强。

[0114] 申请人结合说明书附图对本发明的实施例做了详细的说明与描述，但是本领域技术人员应该理解，以上实施例仅为本发明的优选实施方案，详尽的说明只是为了帮助读者更好地理解本发明精神，而并非对本发明保护范围的限制，相反，任何基于本发明的发明精神所作的任何改进或修饰都应当落在本发明的保护范围之内。

Claims (4)

1. 一种火电厂多自动调节回路评价诊断方法，所述方法针对火电厂多个自动调节回路 进行控制性能的评价并对自动调节回路中的传感器和执行器进行故障诊断；其特征在于， 所述方法包括以下步骤： 步骤一:通过火电厂热工录波分析装置采集火电厂多个自动调节回路数据信息，包括 机组负荷、自动调节回路中被控对象的设定值、传感器采集的测量值以及执行器开度指令、 开度反馈值； 步骤二：以30秒为周期，由火电厂热工录波分析装置对采集数据信息进行处理，形成 comtrade录波文件，内容包含30秒内采集的数据信息，形成录波文件之后，将其发送至作为 诊断装置的电脑； 步骤三：在诊断装置中，以30秒为周期解析由火电厂热工录波分析装置生成的 comtrade录波文件，得到火电厂多个自动调节回路的数据信息，并存入历史数据库； 步骤四：将需要进行评价诊断的火电厂自动调节回路做排列，依据采集的自动调节回 路的数据信息以周期形式依次对单个自动调节回路进行诊断、性能评价，其中对传感器、执 行器的诊断周期为60秒，S卩60秒内，依次对所有自动调节回路中的传感器、执行器进行诊 断;对所述多个自动调节回路的控制性能评价的周期为10分钟，即10分钟内，依次对所有自 动调节回路的控制性能进行评价指标计算： 其中，对传感器、执行器故障诊断过程如下： (1) 以60秒为周期，从历史数据库中提取待诊断的传感器、执行器所属的自动调节回路 的数据信息，包括60秒内待诊断的传感器采集的测量值历史记录、执行器的开度指令和开 度反馈值历史记录，并形成CSV格式文本文件； (2) 根据生成的csV格式文本文件的内容得到60秒内待诊断的传感器采集的测量值，基 于所述测量值实现对传感器的诊断； 当传感器测量值超出了其高低限值范围，则认为传感器发生了测量值超限故障； 计算传感器所采集的测量值的变化速率，将所述测量值的变化速率与设定的变化率阈 值比较，如果超限，则判断传感器存在变化速率超限故障； 两个传感器冗余设置时，先比较两个冗余设置的传感器在同一时刻的测量值偏差，若 超过了偏差限值，则说明其中有一个发生了故障的，然后判断这两个冗余设置的传感器是 否发生了测量值超限故障或变化率超限故障，如果两个冗余设置的传感器均没有发生测量 值超限故障以及变化率超限故障，则认为这两个冗余传感器存在冗余故障； 三冗余及三个以上冗余设置的传感器，比较每两个传感器间测量值的偏差，如果偏差 超过了偏差限值，则认为此两个传感器中某个传感器发生了故障，再与另一个传感器的测 量值相比较，判断偏差大者为发生故障的那个传感器； (3) 根据生成的CSV文件内容，得到执行器开度指令和开度反馈历史记录，判断执行器 是否存在卡死故障、增益故障和死区故障； 当满足如下两个判据中的任一个，则判断执行器出现了卡死故障： 判据1:在60秒内，执行器开度指令发生变化，而开度反馈不随执行器开度指令变化； 判据2:当执行器开度指令不变化时，开度反馈也不变化，但是开度反馈和执行器开度 指令之间存在偏差，偏差范围超过开度反馈量程的5% ; 统计60秒内执行器开度反馈的变化程度与执行器开度指令变化程度均值的比值，作为 执行器增益的估计值;如果这个增益值超过了预设的增益阈值，则认为执行器发生了增益 故障； 根据执行器在正常情况下的往返行程中X-Y特性图，选择最大宽度为正常死区范围来 判断执行器的死区故障。

2. 根据权利要求1所述的火电厂多自动调节回路评价诊断方法，其特征在于： 在步骤四中，对所述多个自动调节回路的控制性能评价过程如下： (1) 以10分钟为周期执行，对已排序的多个自动调节回路依次计算，首先要检查该控制 回路是否存在执行器或者传感器故障，如果存在故障则将此回路跳过，不进行控制评价指 标计算，进行下一个自动调节回路的计算过程，如无故障则进入下一个流程； (2) 将10分钟内自动调节回路中传感器测量值、执行器开度指令、自动调节回路中被控 对象的设定值以及机组负荷的历史值写入csv格式文件中； (3) 自动调节回路中被控对象的设定值发生变化，则认为该自动调节回路开始一个动 态过程，否则，则认为该自动调节回路处于一个随机过程；如果自动调节回路处于动态过 程，则计算该自动调节回路的确定性调节指标，包括调节时间、测量值上升时间、超调量;如 果自动调节回路处于随机过程，则进行随机性指标即最小方差指标计算。

3. 根据权利要求2所述的火电厂多自动调节回路评价诊断方法，其特征在于： 计算自动调节回路的确定性调节指标过程如下： (1)建立自动调节回路中被控对象的数学模型： 自动调节回路的输入一输出模型的差分方程为：

其中：e⑴为噪声，p的为输入一输出观测向量，y⑴为自动调节回路的测量值，u⑴为 自动调节回路中的阀位反馈，Θ为未知参数向量，共有na，个输出、nb个输入，ai......ana， bi bnb为其对应参数，d为迟延参数； 且

令根据N次观测{y(i)，u(i) :i = l，2,…具

对Θ的估计为

，引入符号：

y⑴，y⑵一y (N)为连续N次采样得到的测量值； 得到对参数的最小二乘估计为：

其中：

为参数的最小二乘估计，得到的矩阵依次为ai、a2……ana的值； (2)根据前述的自动调节回路的数学模型数学模型，施加单位阶跃扰动信号，形成理想 的输出曲线，根据调节曲线，得到每个时刻对应的输出值，在此基础上，计算得到上升时间、 调节时间和超调量。

4.根据权利要求2所述的火电厂多自动调节回路评价诊断方法，其特征在于： 计算自动调节回路的随机性调节指标的过程如下： (1) 通过火电厂热工录波分析装置采集火电厂多个自动调节回路数据信息，包括机组 负荷、自动调节回路中被控对象的设定值、传感器采集的测量值以及执行器开度指令、开度 反馈值； (2) 在采集的数据中，选择数据长度至少为n = 500且包含系统典型的扰动特性的一段 数据进行评价，η即表示数据长度； ⑶进行开环测试或用闭环估计技术估计采样数据的时间延迟k; ⑷自动调节回路的输出y⑴表示为下式

αι......<^为模型的参数； k为迟延因子，由自动调节回路的迟延时间决定； m为自回归模型的长度，m取5〜30,模型长度取决于采样时间，并且模型长度的选择要 使最终模型能全部捕捉系统的脉冲响应； Fl (t)为自动调节回路输出在t时刻的响应，其余项为扰动信号对输出的响应； (5) 由（4)可以获得数据长度η范围内的自动调节回路中传感器采集的测量值，叠加写 成以下形式

由获得数据根据上式构建上述向量y和矩阵X; y、X矩阵中值为对应下标时刻的自动调节回路中传感器采集的测量值； (6) 选择合适的模型获得〜的参数； 用线性回归方法估计自回归模型参数1^}，最小二乘结果为：

⑵计算自动调节回路的最小方差，最小方差的估计值可以由残余方差代替：

是理论的最小方差； ⑶计算系统实际的输出方差：

是实际的输出方差

是输出方差中出去理论最小方差的部分； (9)计算随机性调节指标：

是最小方差指标，由上式可知，其取值范围为〇_1，此值越接近1，则表不实际方差与 理论最小方差越接近，说明自动调节回路抗随机扰动能力更强。