本发明公开了一种协调控制被控对象传递函数模型辨识方法，包括DCS数据采样、应用RBF神经网络建模；判定对象的类型，选择需要辨识的参数；确定初始对象，以及初始对象的均匀随机响应；确定RBF神经网络输出；以及计算误差函数等步骤。本发明充分利用了大量的运行数据，避免了信息资源的浪费，同时节省了现场试验所花费的人力、物力以及现场条件的限制，并且辨识得到的传递函数模型精度明显优于现场试验得到的传递函数模型。
1.一种协调控制被控对象传递函数模型辨识方法，其特征在于包括以下步骤：步骤一，DCS数据采样：
①选择同一负荷下各时段的DCS数据，采样时间T为5秒，包括：发电机实发功率、锅炉主蒸汽压力、汽轮机调门开度和锅炉燃料率指令，分别用变量NE、PT、uT、uB表示；
②在各连续时段中，将变量NE、PT、uT、uB分别作时延处理，分别取各个变量前1～n个时刻的数值，形成样本数据：其中，n为3～5的自然数，N＝3000；
作为模型的输入，NE(k)；PT(k)，k＝n+1，n+2，...，N，作为模型的输出，建立RBF神经网络模型；
步骤三，判定对象的传递函数类型，确定需要辨识的参数：根据热工过程的特性，总结归纳出2×2的协调控制对象的传递函数模型为：其中，模型中的参数k1、k2、k3、k4表示对象的稳态增益，T1、T2、T3表示对象的惯性时间，τ表示滞后时间，对模型中的参数k1、k2、k3、k4、T1、T2、T3和τ进行辨识；
①根据步骤三中的传递函数模型WNB(s)，WPB(s)，选取初始参数k2、k4、T1、T2、τ，得到初始对象W′NB(s)，W′PB(s)；
②对于初始对象W′NB(s)，在其输入端加入均匀随机信号uB(i)，i＝1，2，...，M，M＝
2000，得到初始对象的输出响应y′N(j)，j＝1，2，...，M；
③同理，对于初始对象W′PB(s)，在其输入端加入步骤四②中同样的均匀随机信号uB(i)，i＝1，2，...，M，得到W′PB(s)的输出响应y′P(j)，j＝1，2，...，M；
其中，W′NB(s)，W′PB(s)分别是WNB(s)，WPB(s)的初始值；
步骤五，确定RBF神经网络输出：对于步骤二中的RBF神经网络模型，使uT(i)＝0，i＝1，
2，...，M，且uB(i)，i＝1，2，...，M为步骤四②中同样的均匀随机信号，将y′N(i)，y′P(i)，uT(i)，uB(i)，i＝1，2，...，M，经过步骤一②中同样的时延处理，得到：作为步骤二中的RBF神经网络模型的输入端，从而得到RBF神经网络模型的输出yN(j)；
将E作为目标函数，利用遗传算法进行迭代寻优，得到最佳的对象参数k2opt、k4opt、T1opt、T2opt、τopt，从而确定WNB(s)，WPB(s)；k2opt、k4opt表示对象的最佳稳态增益，T1opt、T2opt表示对象的最佳惯性时间，τopt表示最佳滞后时间；
时延处理中，使得最大的延迟时间nT为惯性时间，n取3～5的自然数；
②初始种群选择：即k1、k2、k3、k4、T1、T2、T3和τ的上下界，总结出k1∈[0，2000]，k2∈[0，
10]，k3∈[-1，0]，k4∈[0，1]，T1∈[0，500]，T2∈[0，500]，T3∈[0，50]，τ∈[0，50]；
④遗传操作：采用轮盘赌选择法、均匀交叉、均匀变异的策略；
步骤七：使uB(i)＝0，i＝1，2，...，M，uT(i)，i＝1，2，...，M均匀随机变化，重复步骤四、步骤五、步骤六，得到最佳的对象参数k1opt、k3opt、T3opt，k1opt、k3opt表示对象的最佳稳态增益，T3opt表示对象的最佳惯性时间，从而确定WNT(s)，WPT(s)。
一种协调控制被控对象传递函数模型辨识方法
技术领域
[0001] 本发明涉及热工过程控制领域，尤其涉及热工过程的辨识方法。
背景技术
[0002] 热工过程具有非线性、时延性、不确定性以及变量间的关联性等特点，难以建立精确的数学模型，这给热工控制带来了很大的难度。传统的机理建模得到的模型与现场设备在性能上有一定的差距，因为现场的绝大多数设备的性能都是非线性，不可能通过数学公式来得到精确的模型。另外由于现场中存在着较大的随机扰动，对现场试验有很大的影响。鉴于热工试验的复杂性和局限性，在现场获得多变量对象的精确模型是很困难的。
[0003] 现代化的电厂装配着庞大的DCS系统，日积月累采集了大量的运行数据，如何从蕴藏着大量信息的数据中提取各变量间的规律，成为目前的研究热点。如何找到合适的模型并求取多变量系统精确的传递函数矩阵模型，具有很现实的意义。
发明内容
[0004] 发明目的：为了克服现有技术中存在的不足，本发明提供一种协调控制被控对象传递函数模型辨识方法，利用现场的DCS数据建立神经网络模型，辨识传递函数，充分利用了大量的运行数据，避免了信息资源的浪费，得到的传递函数模型精度明显优于现场试验得到的传递函数模型，同时节省了现场试验所花费的人力、物力以及现场条件的限制。
[0005] 技术方案：为解决上述技术问题，本发明提供的一种协调控制被控对象传递函数模型辨识方法，包括以下步骤：
[0006] 步骤一，DCS数据采样；
[0007] 步骤二，应用RBF神经网络建模；
[0008] 步骤三，判定对象的类型，选择需要辨识的参数；
[0009] 步骤四，确定初始对象，以及初始对象的均匀随机响应；
[0010] 步骤五，确定RBF神经网络输出；
[0011] 步骤六，计算误差函数。
[0012] 具体地，所述步骤一的DCS数据采样方法为：
[0013] ①选择同一负荷下各时段的DCS数据，采样时间T为5秒，包括：发电机实发功率、锅炉主蒸汽压力、汽轮机调门开度和锅炉燃料率指令，分别用变量NE、PT、uT、uB表示；
[0014] ②在各连续时段中，将变量NE、PT、uT、uB分别作时延处理，分别取各个变量前1～n个时刻的数值，形成样本数据：
[0015]
[0016] 其中，n为3～5的自然数，N＝3000。
[0017] 具体地，所述步骤二的RBF神经网络建模是将：
[0018]
[0019] 作为模型的输入，NE(k)；PT(k)，k＝n+1，n+2，...，N，作为模型的输出，建立RBF神经网络模型。
[0020] 具体地，所述步骤三的判断选择过程如下：
[0021] 总结归纳出2×2的协调控制对象的传递函数模型为：
[0022]
[0023] 由该模型得出需要辨识的参数为k1、k2、k3、k4、T1、T2、T3和τ。
[0024] 具体地，所述步骤四的过程如下：
[0025] ①根据步骤三中的传递函数模型WNB(s)，WPB(s)，选取初始参数k2′、k4′、T1′、T2′、τ′，得到初始对象W′NB(s)，W′PB(s)；
[0026] ②对于初始对象W′NB(s)，在其输入端加入均匀随机信号uB(i)，i＝1，2，...，M，M＝2000，得到初始对象的输出响应y′N(j)，j＝1，2，...，M；
[0027] ③同理，对于初始对象W′PB(s)，在其输入端加入步骤四②中同样的均匀随机信号uB(i)，i＝1，2，...，M，得到W′PB(s)的输出响应y′P(j)，j＝1，2，...，M。
[0028] 具体地，所述步骤五的过程为：对于步骤二中的RBF神经网络模型，使uT(i)＝0，i＝1，2，...，M，且uB(i)，i＝1，2，...，M为步骤四②中同样的均匀随机信号，将y′N(i)，y′P(i)，uT(i)，uB(i)，i＝1，2，...，M，经过步骤二中同样的时延处理，得到：
[0029]
[0030] 作为步骤二中的RBF神经网络模型的输入端，从而得到RBF神经网络模型的输出yN(j)；yP(j)，j＝n+1，n+2，...，M。
[0031] 具体地，所述步骤六计算误差函数的过程为：
[0032]
[0033] 将E作为目标函数，利用遗传算法进行迭代寻优，得到最佳的对象参数k2opt、k4opt、T1opt、T2opt、τopt，从而确定WNB(s)，WPB(s)；
[0034] 还包括步骤七：使uB(i)＝0，i＝1，2，...，M，uT(i)，i＝1，2，...，M均匀随机变化，重复步骤四、步骤五、步骤六得到最佳的对象参数k1opt、k3opt、T3opt，从而确定WNT(s)，WPT(s)。
[0035] 具体地，所述时延处理中，使得最大的延迟时间nT为惯性时间，n取3～5的自然数。
[0036] 具体地，所述遗传算法的设置：
[0037] ①实数编码机制；
[0038] ②初始种群选择：即k1、k2、k3、k4、T1、T2、T3和τ的上下界，总结出k1∈[0，2000]，k2∈[0，10]，k3∈[-1，0]，k4∈[0，1]，T1∈[0，500]，T2∈[0，500]，T3∈[0，50]，τ∈[0，50]；③适应度函数： E为误差函数；
[0039] ④遗传操作：采用轮盘赌选择法、均匀交叉、均匀变异的策略。
[0040] 有益效果：本发明利用现场大量的DCS运行数据，建立神经网络模型，利用神经网络模型辨识单元机组协调控制对象的传递函数模型。本发明充分利用了大量的运行数据，避免了信息资源的浪费，同时节省了现场试验所花费的人力、物力以及现场条件的限制，并且辨识得到的传递函数模型精度明显优于现场试验得到的传递函数模型。
[0041] 除了上面所述的本发明解决的技术问题、构成技术方案的技术特征以及由这些技术方案的技术特征所带来的优点外，本发明的一种协调控制被控对象传递函数模型辨识方法所能解决的其他技术问题、技术方案中包含的其他技术特征以及这些技术特征带来的优点，将结合附图做出进一步详细的说明。
附图说明
[0042] 图1：遗传算法辨识神经网络模型传递函数原理图；
[0043] 图2：WPB(s)阶跃响应的真实值和预测值的对比图；
[0044] 图3：WNB(s)阶跃响应的真实值和预测值的对比图；
[0045] 图4：WPT(s)阶跃响应的真实值和预测值的对比图；
[0046] 图5：WNT(s)阶跃响应的真实值和预测值的对比图。
具体实施方式
[0047] 实施例：
[0048] 下面结合某电厂600MW亚临界燃煤单元机组的现场DCS数据作为例子，说明本发明的技术方案实施过程如下：
[0049] 步骤一：选择100％负荷下各个时段的DCS数据，采样时间为5s，包括：发电机实发功率、锅炉主蒸汽压力、汽轮机调门开度和锅炉燃料率指令，分别用变量NE、PT、uT、uB表示；
[0050] ②在各连续时段中，将变量NE、PT、uT、uB分别作时延处理，分别取各个变量前5个时刻的数值，形成样本数据：
[0051]
[0052] 步骤二：应用RBF神经网络建模，其遗传算法辨识神经网络模型传递函数原理如图1所示，将：
[0053]
[0054] 作为模型的输入，NE(k)；PT(k)，k＝n+1，n+2，...，3000，作为模型的输出，建立RBF神经网络模型；
[0055] 步骤三：总结归纳出2×2的协调控制对象的传递函数模型为：
[0056]
[0057] 由以上的模型可知，只需确定k1、k2、k3、k4、T1、T2、T3和τ，即可确定对象的四个传递函数模型；
[0058] 步骤四：确定初始对象，以及初始对象的均匀随机响应；
[0059] ①根据步骤三中的传递函数模型，如WNB(s)，WPB(s)，选取初始参数k2′、k4′、T1′、T2′、τ′，得到初始对象W′NB(s)，W′PB(s)；
[0060] ②对于初始对象W′NB(s)，在其输入端加入均匀随机信号uB(i)，i＝1，2，...，2000，，得到初始对象的输出响应y′N(j)，j＝1，2，...，M；
[0061] ③同理，对于初始对象W′PB(s)，在其输入端加入步骤四②中同样的均匀随机信号uB(i)，i＝1，2，...，2000，得到W′PB(s)的输出响应y′P(j)，j＝1，2，...，2000；
[0062] 步骤五：确定RBF神经网络输出；
[0063] 对于步骤二中的RBF神经网络模型，使uT(i)＝0，i＝1，2，...，2000，且uB(i)，i＝1，2，...，2000为步骤四②中同样的均匀随机信号，将y′N(i)，y′P(i)，uT(i)，uB(i)，i＝1，
2，...，2000，经过步骤二中同样的时延处理，得到：
[0064]
[0065] 作为步骤二中的RBF神经网络模型的输入端，从而得到RBF神经网络模型的输出yN(j)；yP(j)，j＝，n+1，n+2，...，2000；
[0066] 步骤六：计算误差函数：
[0067]
[0068] 将E作为目标函数，利用遗传算法进行迭代寻优，遗传算法的参数设置如下：
[0069] ①实数编码机制；
[0070] ②初始种群选择：即k2、k4、T1、T2和τ的上下界，其中k2∈[0，10]，k4∈[0，1]，T1∈[0，500]，T2∈[0，500]，τ∈[0，50]；
[0071] ③适应度函数： E为误差函数；
[0072] ④遗传操作：采用轮盘赌选择法、均匀交叉、均匀变异的策略。
[0073] 得到最佳参数k2opt＝8.06、k4opt＝0.30、T1opt＝140.23、T2opt＝190、τopt＝37.07，得到
[0074] 步骤七：同理，使uB(i)＝0，i＝1，2，...，M，uT(i)，i＝1，2，...，M均匀随机变化，重复步骤四、步骤五，得到最佳的对象参数k1opt＝1035、k3opt＝-0.17、T3opt＝24.95，得到[0075]
[0076] 得到对象的传递函数模型：
[0077]
[0078] 分别对四个传递函数作阶跃响应，与现场真实的阶跃响应进行比较，得到如图2的WPB(s)阶跃响应的真实值和预测值的对比图，如图3的WNB(s)阶跃响应的真实值和预测值的对比图；如图4的WPT(s)阶跃响应的真实值和预测值的对比图；以及图5的WNT(s)阶跃响应的真实值和预测值的对比图，从图中可以看出辨识得到的传递函数模型精度较高。
