[发明专利]采用改进的遗传优化算法精确辨识热工过程状态空间模型参数的方法

摘要

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

主权项

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为自然数；确定辨识参数(c₁,c₂,...,c_z),其中，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)＝(y₁(k),y₂(k),...,y_q(k)),1≤k≤N，故可得适应度值J：上式中，αj为第j个输出变量对应的权值；步骤四、采用二进制进行编码；步骤五、最优染色体高频变异；通过对适应度值进行排序，可获得最佳适应度值及最优染色体；选取一定数量的染色体n1,对最优染色体上的二进制基因，进行高频变异，其概率为Pmax,使得0→1,1→0，并将变异过的n1个染色体，来取代适应度最差的前n1个染色体；步骤六、最优染色体保存机制将初代最佳适应度值保留下来作为全局最佳适应度值，并选其对应的染色体为全局最佳染色体，在以后的进化过程中，不断将本代的最佳适应度值与全局的最佳适应度值进行比较，如果全局的适应度值较小，则将全局最佳染色体注入当代种群中；如果当代最佳适应度值较小，将当代最佳染色体和适应度值更新全局的最佳染色体和适应度值。从而，最优染色体可以保存下来；步骤七、选择算子本算法中选择算子选用轮盘赌策略；令其中PP_i为累计概率，pp_i为个体的选择概率，其计算公式为：其中fitness(xi)为个体的适应度值；共转轮popsize次，popsize为染色体个数，每次转轮时，随机产生0到1之间的随机数r，当PPi‑1≤r≤PPi时选择个体i；故从选择概率的计算公式可以看出，个体的适应值越大，其选择概率越大；步骤八、交叉和变异算子步骤九、解码：将二进制转为十进制步骤十、自适应空间变异在全局最优解陷入停滞时，对最优解进行一定范围大小的实数空间变异，并将所变异的染色体组n2注入当前种群中，且实数变异的空间范围会随着停滞的代数增加而增加，即自适应变异因子，a增大，并对所变异的最优解进行限幅操作，来满足可行性，提高算法局部搜索能力， 当跳出局部最优解时，a恢复到初始值；然后，进化下一代，达到算法最大进化代数时，停止计算，输出最佳染色体和适应度值。