采用改进的遗传优化算法精确辨识热工过程状态空间模型参数的方法技术

本发明专利技术公开了一种采用改进的遗传优化算法精确辨识热工过程状态空间模型参数的方法，其主要步骤是确定模型的结构与辨识参数、确定优化算法结构参数、求得适应度值、编码，十进制转为二进制、最优染色体高频变异、最优染色体保存机制、选择、交叉和变异、解码，二进制转为十进制和自适应空间变异。本发明专利技术引入最优染色体保存机制可避免算法后期随机搜索；最优染色体高频变异，可增强算法全局寻优能力；自适应空间变异，对全局最优解进行一定范围的实数空间变异，且空间变异的范围会随着陷入局部最优的代数增加而增加，直到跳出局部最优解，可增强算法局部寻优的能力。

An improved genetic optimization algorithm for accurate identification of the parameters of the state space model in thermal processes

The invention discloses a method of using genetic algorithm to accurately identify parameter state space model of thermal process, the main step is to determine the structure and parameter identification, to determine the model optimization of structure parameters, obtained the fitness value, encoding, decimal to two hexadecimal, optimal chromosome variation, high frequency optimal preservation mechanism, chromosome selection, crossover and mutation, decoding, binary to decimal and adaptive spatial variation. The invention introduces optimal chromosome preservation mechanism can avoid the random search algorithm later; optimal chromosome hypermuation, can enhance the ability of global optimization algorithm; adaptive spatial variability, spatial variability of real range of the global optimal solution, and the range of spatial variability will increase with the increase of the algebra into a local optimal, until out of local optimal solution that ability can enhance the local search algorithm.

【技术实现步骤摘要】
采用改进的遗传优化算法精确辨识热工过程状态空间模型参数的方法
本专利技术属于工程
，尤其是一种采用改进的遗传优化算法精确辨识热工过程状态空间模型参数的方法。
技术介绍
现实热工过程模型多为非线性，大迟延，大惯性，非线性状态空间模型。当前一般采用基本遗传算法来实现。基本遗传算法(SGA)是一种基于适者生存的自然选择机制和生物遗传的进化搜索算法，其主要特点是群体搜索策略，具有较强的全局搜索能力且对被优化的数学模型无先验要求，被广泛用于自动控制、图像识别、机器学习和故障诊断等领域，但同时也暴露出其理论和技术上的许多不足和缺陷。SGA中的三个基本算子，选择算子、交叉算子和变异算子。交叉与变异算子可以通过交叉和变异概率来产生新的抗体，但其搜索盲目性和随机性，故求出全局最优解的概率不高，且对最优个体无保存机制，在进化的中后期如随机搜索。因此，基本遗传算法(SGA)对高维、非线性的复杂系统求解问题，往往产生早熟现象。
技术实现思路
有鉴于此，为解决现有技术存在的上述问题，申请人提供了一种采用改进的遗传优化算法精确辨识热工过程状态空间模型参数的方法。具体地，该技术方案为：一种采用改进的遗传优化算法精确辨识热工过程状态空间模型参数的方法，包括如下步骤：步骤一、确定模型的结构与辨识参数热工过程状态空间模型写为：为n维状态变量的导数，X为n维状态变量，n为状态变量个数，A＝A(X)，A(X)为关于状态变量X的函数，n×n阵A为系统矩阵，B＝B(X)，B(X)为关于状态变量X的函数，n×p阵B为输入矩阵，p为输入变量个数，C＝C(X)，C(X)为关于状态变量X的函数，q×n阵C为输出矩阵，q为输出变量个数，D＝D(X)，D(X)关于状态变量X的函数,q×p阵D为传输矩阵，U为p维输入，Y为q维输出，n,q,p为自然数；确定辨识参数(c1,c2,...,cz),其中，c为辨识参数，z为辨识参数个数；步骤二、确定优化算法结构参数确定染色体个数popsize，辨识参数个数z，二进制码长[L1,L2,...,Lz]，迭代次数，G，交叉概率Pc,变异概率Pm，高频变异概率Pmax，高频变异染色体组数n1，自适应变异染色体组数n2，a自适应变化因子，辨识参数变化范围，最小范围MinX＝[xmin1,xmin2,...,xminz]；最大范围MaxX＝[xmax1,xmax2,...,xmaxz],xminj和xmaxj为第j个辨识参数对应的最小值和最大值，1≤j≤z，z为辨识参数个数；步骤三、求得适应度值J在算法寻优开始时，初始化染色体，即获得初始代的辨识参数(c10,c20,...,cz0)，将当前输入量U(k)＝(u1(k),u2(k),...,up(k))，p为输入变量个数，k为第k时刻的运行数据点，1≤k≤N，N为实测数据总数，k、N为自然数；在超临界协调控制系统模型，共3个系统输入变量，3个系统状态变量，3个系统输出变量；u1为燃料量指令；u2为给水流量；u3为汽轮机调门开度；x1为入炉煤量；x2为汽水分离器焓值；x3为汽水分离器压力；y1为主蒸汽压力；y2为汽水分离器焓值；y3为机组功率；对状态变量导数进行离散化，即n为状态变量个数，T为采样周期；已知当前输入量U(k)和当前状态变量X(k)，即可获得下一时刻状态变量，X(k+1)，将X(k+1)和U(k)代入Y＝C(X)+D(X)U，即可获得下一时刻模型输出即q为输出变量个数；再将X(k+1)作为当前状态变量X(k)；同理，令下一时刻U(k+1)作为U(k)，迭代往复，即可获得模型输出u(k)为系统实际测得输入量，y(k)为系统的实际输出，为模型计算输出，k为第k时刻的运行数据点，N为实测数据总数，再结合实测运行数据Y(k)＝(y1(k),y2(k),...,yq(k)),1≤k≤N，故可得适应度值J：上式中，αj为第j个输出变量对应的权值；步骤四、采用二进制进行编码；步骤五、最优染色体高频变异；通过对适应度值进行排序，可获得最佳适应度值及最优染色体；选取一定数量的染色体组，n1,对最优染色体上的二进制基因，进行高频变异，其概率为Pmax,使得0→1,1→0(0转换成1,1转换成0)，并将变异过的染色体n1组,来取代适应度最差的前n1个染色体；步骤六、最优染色体保存机制将初代最佳适应度值保留下来作为全局最佳适应度值，并选其对应的染色体为全局最佳染色体，在以后的进化过程中，不断将本代的最佳适应度值与全局的最佳适应度值进行比较，如果全局的适应度值较小，则将全局最佳染色体注入当代种群中；如果当代最佳适应度值较小，将当代最佳染色体和适应度值更新全局的最佳染色体和适应度值。从而，最优染色体可以保存下来；步骤七、选择算子本算法中选择算子选用轮盘赌策略；令PP0＝0，其中PPi为累计概率，ppi为个体的选择概率，其计算公式为：其中fitness(xi)为个体的适应度值；共转轮popsize次，popsize为染色体个数，每次转轮时，随机产生0到1之间的随机数r，当PPi-1≤r≤PPi时选择个体i；故从选择概率的计算公式可以看出，个体的适应值越大，其选择概率越大；步骤八、交叉和变异算子步骤九、解码：将二进制转为十进制步骤十、自适应空间变异在全局最优解陷入停滞时，对最优解进行一定范围大小的实数空间变异，并将所变异的染色体组n2注入当前种群中，且实数变异的空间范围会随着停滞的代数增加而增加，即自适应变异因子,a增大，并对所变异的最优解进行限幅操作，来满足可行性，提高算法局部搜索能力。当跳出局部最优解时，a恢复到初始值；然后，进化下一代，达到算法最大进化代数时，停止计算，输出最佳染色体和适应度值。实施本专利技术，可获得的有益效果是：1、本技术将热工过程状态空间模型参数的辨识问题转化为相关参数的组合寻优问题，利用改进的遗传优化算法对整个参数空间进行多组参数全局优化同步搜索，并将所求的下一时刻的状态量与输入量代入模型中，计算下一时刻模型输出，反复更新、迭代，与实际输出值计算适应度值，多组个体不断寻优进化，从而获得系统模型参数的最优估计。2、与基本遗传算法相比，本专利技术将最优个体保存机制引入改进的算法中，可避免后期随机搜索现象。3、与基本遗传算法相比，本专利技术将最优个体高频变异优化算子引入算法中，大大提高染色体变异的成功率，增强了全局寻优的能力。4、本专利技术引入自适应空间变异，对全局最优解进行一定的实数空间范围变异，且空间变异的范围会随着陷入局部最优的代数增加而增加，直到跳出局部最优解，可增强局部寻优的能力。附图说明图1是本专利技术的流程图。图2展示了遗传算法中适应度值的变化。图3展示了改进的遗传算法适应度值的变化。图4是主蒸汽压力算法比较图。图5是汽水分离器焓值算法比较图。图6是机组功率算法比较图。具体实施方式本专利技术的主要步骤如下：定义适应度值为：上式中，q,输出变量个数，αj第j个输出变量对应的权值,y(k)为系统的实际输出，为模型计算输出,N为实测数据总数。改进的遗传算法寻优，即寻找一组辨识参数,使得适应度值J为最小。热工过程状态空间模型辨识方案流程设计如下：1、确定模型的结构与辨识参数状态空间模型一般可写为其中，为n维状态变量的导数，X为n维状态变量，n为状态变量个数，A＝A(X)，A(X)为本文档来自技高网...

【技术保护点】
一种采用改进的遗传优化算法精确辨识热工过程状态空间模型参数的方法，其特征在于，包括如下步骤：步骤一、确定模型的结构与辨识参数热工过程状态空间模型写为：

【技术特征摘要】
1.一种采用改进的遗传优化算法精确辨识热工过程状态空间模型参数的方法，其特征在于，包括如下步骤：步骤一、确定模型的结构与辨识参数热工过程状态空间模型写为：为n维状态变量的导数，X为n维状态变量，n为状态变量个数，A＝A(X)，A(X)为关于状态变量X的函数，n×n阵A为系统矩阵，B＝B(X)，B(X)为关于状态变量X的函数，n×p阵B为输入矩阵，p为输入变量个数，C＝C(X)，C(X)为关于状态变量X的函数，q×n阵C为输出矩阵，q为输出变量个数，D＝D(X)，D(X)关于状态变量X的函数,q×p阵D为传输矩阵，U为p维输入，Y为q维输出，n,q,p为自然数；确定辨识参数(c1,c2,...,cz),其中，c为辨识参数，z为辨识参数个数；步骤二、确定优化算法结构参数确定染色体个数popsize，辨识参数个数z，二进制码长[L1,L2,...,Lz]，迭代次数，G，交叉概率Pc,变异概率Pm，高频变异概率Pmax，高频变异染色体组数n1，自适应变异染色体组数n2，a自适应变化因子，辨识参数变化范围，最小范围MinX＝[xmin1,xmin2,...,xminz]；最大范围MaxX＝[xmax1,xmax2,...,xmaxz],xminj和xmaxj为第j个辨识参数对应的最小值和最大值，1≤j≤z，z为辨识参数个数；步骤三、求得适应度值J在算法寻优开始时，初始化染色体，即获得初始代的辨识参数(c10,c20,...,cz0)，将当前输入量U(k)＝(u1(k),u2(k),...,up(k))，p为输入变量个数，k为第k时刻的运行数据点，1≤k≤N，N为实测数据总数，k、N为自然数；在超临界协调控制系统模型，共3个系统输入变量，3个系统状态变量，3个系统输出变量；u1为燃料量指令；u2为给水流量；u3为汽轮机调门开度；x1为入炉煤量；x2为汽水分离器焓值；x3为汽水分离器压力；y1为主蒸汽压力；y2为汽水分离器焓值；y3为机组功率；对状态变量导数进行离散化，即n为状态变量个数，T为采样周期；已知当前输入量U(k)和当前状态变量X(k)，即可获得下一时刻状态变量，X(k+1)，将X(k+1)和U(k)代入Y＝C(X)+D(X)U，即可获得下一时刻模型输出即q为输出变量个数；再将X(k+1)作为当前状态变量X(k)；同...

【专利技术属性】
技术研发人员：范赫，张雨飞，
申请(专利权)人：东南大学，
类型：发明
国别省市：江苏,32