CN105680071A - 基于分数阶滑模变结构sofc系统热电协同控制方法 - Google Patents

Abstract

本发明提供基于分数阶滑模变结构SOFC系统热电协同控制方法，包括以下步骤：S1、采集SOFC系统不同输入参数组合下的系统状态及输出参数，得到全负载区间稳态功率与温度、效率响应特性、旁路阀开度BP对效率优化的影响函数和指定负载切换区间及延时条件下的效率优化函数；S2、获取系统的局部最优稳态操作函数、拓展形成稳态下全局优化函数和不同切换负载区间不同延时条件下的功率跟踪函数；S3、计算得出滑模区间；S4、优化函数计算得出系列趋近律函数；S5、通过分数阶优化法对所述系列趋近律函数消除抖震，通过计算求解出趋近律。本方法控制精准、灵活、稳定，能大大的提速切换过程、克服SOFC系统大惯性的时滞特点，实现快速的负载切换。

Description

基于分数阶滑模变结构SOFC系统热电协同控制方法

技术领域

本发明涉及基于分数阶滑模变结构SOFC系统热电协同控制方法。

背景技术

燃料电池它是把氢和氧反应生成水放出的化学能转换成电能的装置。其基本原理相当于电解反应的逆向反应。燃料(H2和CO等)及氧化剂(O2)在电池的阴极和阳极上借助氧化剂作用，电离成离子，由于离子能通过在二电极中间的电介质在电极间迁移，在阴电极、阳电极间形成电压。在电极同外部负载构成回路时就可向外供电。

燃料电池种类不少，根据使用的电解质不同，主要有磷酸燃料电池(PAFC)、熔融碳酸盐型燃料电池(MCFC)、固体氧气物燃料电池(SOFC)和质子交换膜燃料电池(PEMFC)等。

燃料电池具有无污染、高效率、适用广、无噪声和能连续运转等优点。它的发电效率可达40％以上，热电联产的效率也达到80％以上。多数燃料电池正处于开发研制中，虽然磷酸燃料电池(PAFC)等技术成熟并已经推向市场，但仍较昂贵。鉴于燃料电池的独到优点，随着该项技术商业化进程的推进，必将在未来燃气采暖行业起到越来越重要的作用。

发明内容

本发明所要解决的技术问题是提供一种控制精准、灵活、稳定，能大大的提速切换过程、克服SOFC系统大惯性的时滞特点，实现快速的负载切换的基于分数阶滑模变结构SOFC系统热电协同控制方法。

本发明解决上述技术问题的技术方案如下：基于分数阶滑模变结构SOFC系统热电协同控制方法，包括以下步骤：

S1、采集SOFC系统不同输入参数组合下的系统状态及输出参数，通过系统辨识的方法得到全负载区间稳态功率与温度、效率响应特性、旁路阀开度BP对效率优化的影响函数和指定负载切换区间及延时条件下的效率优化函数；

S2、基于所述全负载区间稳态功率与温度、效率响应特性、旁路阀开度BP对效率优化的影响函数和指定负载切换区间及延时条件下的效率优化函数获取系统的局部最优稳态操作函数、拓展形成稳态下全局优化函数和不同切换负载区间不同延时条件下的功率跟踪函数；

S3、根据所述最优稳态操作函数、拓展形成稳态下全局优化函数和指定负载切换区间及延时条件下的效率优化函数计算得出滑模区间；

S4、根据所述滑模区间、全负载区间稳态功率与效率相应特性、指定负载切换区间及延时条件下的效率优化函数和不同切换负载区间不同延时条件下的改进优化函数计算得出系列趋近律函数；

S5、通过分数阶优化法对所述系列趋近律函数消除抖震，通过计算求解出趋近律。

本发明的有益效果是，燃料利用率FU是行业内公认影响系统发电效率的关键指标，且目前的研究均希望能最大化燃料利用率以提高系统效率；空气过量比AR为输入空气与参与反应的空气流量之间的直接关系表达，有利于在系统控制分析中实现协同联动。

此外，面向负载跟踪，针对不同功率的稳态输出需求，可以快速获取最优操作点(Is，BP，AR，FU)，使其在保障系统热特性安全前提下满足负载需求，并能实现在需求功率输出下实现系统最大效率。

在面向负载跟踪全过程观测器设计中，进一步与系统过程的响应速度进行综合匹配分析，设计出满足系统发电全过程需求的电堆内部温度观测器，以提高动静态控制策略及整体控制架构的精确性与稳定性。

向分数阶滑模变结构SOFC热电管控架构进行映射，此优化稳态操作策略集可将目标滑模面拓展为滑模区间，从而将为以安全、快速、低抖震抵达滑模面为目标的趋近律(即动态切换策略)设计，考虑不同气体延时前提下，提供更大的灵活性与品质保证。

在功率跟踪过程中，功率变化导致系统温度变化的过程需要一段时间才表现出来，所以，只要在进行快速功率跟踪的同时，施加能够抑制功率波动带来的温度变化的控制量即可实现温度约束管控的目的。而本项目分析获得的优化操作集(Is，BP，AR，FU)能够保证效率优化的情况下抑制特定功率下的温度约束；同时，拟采用的包括两个阶段的组合切换策略将可大大的提速切换过程、克服SOFC系统大惯性的时滞特点，实现快速的负载切换

进一步，所述步骤S1中，通过运行电堆的分布式节点模型获取电堆内部温度分布参数。

采用上述进一步的有益效果是，基于已有指定稳态下的电堆内部温度梯度观测方法、进行热电特性响应时间尺度、气固热特性响应时间尺度两层分解，实现负载跟踪全过程中的温度梯度观测的实时性，从而原理上保障动静态控制策略及分数阶滑模变结构控制架构的物理可实现性。在面向负载跟踪全过程观测器设计中，进一步与系统过程的响应速度进行综合匹配分析，以提高动静态控制策略及整体控制架构的精确性与稳定性。

进一步，所述电堆的维数为一维。

采用上述进一步的有益效果是，一维化是为了简化实际电堆模型，这个简化的过程保证了所需要参数是准确的，例如电堆内的温度，温度梯度等。

进一步，所述步骤S1中局部最优稳态操作函数的确定方法,包括以下步骤：

固定工作电流Is和旁路阀开度BP，得到以燃料利用率FU为横坐标、以空气过量比AR为纵坐标的输出性能图谱；

对于所述输出性能图谱，固定工作电流Is，得到不同旁路阀开度下以燃料利用率为横坐标、以空气过量比为纵坐标的输出性能图谱；

遍历不同工作电流Is的值，获取任一操作点(Is，BP，AR，FU)下的系统输出性能；

面向负载跟踪，针对任意功率的稳态输出需求，确定与该需求对应的最优稳态操作函数。

采用上述进一步的有益效果是，由于实际系统物理实现过程中，包括温度、功率、效率等热电特性均与电流直接相关，本发明中基于模型的稳态与动态分析将均基于电流分析模式展开。面向SOFC热电管控需求，为实时调节电堆内部温度分布，需要通入过量空气；为实时跟踪负载需求，需调节燃料利用率；为快速调节电堆内部温度，且优化控制系统发电效率，需实时调节旁路冷空气的输出量，即调节旁路阀开度BP。燃料利用率FU是行业内公认影响系统发电效率的关键指标，且目前该领域的研究均希望能最大化燃料利用率以提高系统效率；空气过量比AR为输入空气与参与反应的空气流量之间的直接关系表达，有利于在系统控制分析中实现协同联动。本发明采用电流控制模式，选取空气过量比(AR)、燃料利用率(FU)，旁路阀开度(BP)、电流(Is)为组合调节量。

进一步，所述步骤S1中电堆内部温度分布参数是通过电堆内部温度梯度观测器来获取的，所述电堆内部温度梯度观测器的构建包括以下步骤：

建立基于线性化状态空间方程涉及降维龙伯格观测器；

假定电化学反应是瞬间完成的，进行准静态假设后，模型为：

$\left\{\begin{array}{c}\frac{{\mathrm{dx}}_1}{dt}=f_1(x_1,y,u,w)\\ \frac{dy}{dt}=f_2\left(x_1,y,u,w\right)\end{array}\right.---\left(1\right)$

其中x1为前四个节点的固体层温度和空气层温度，y为第五个节点的固体层温度和空气层温度；

w＝I_tot，

对以上模型，设计如下形式的观测器：

$\left\{\begin{array}{c}\frac{d{\hat{x}}_1}{dt}=f_1\left({\hat{x}}_1,\hat{y},u,w\right)+Lv\\ \frac{d\hat{y}}{dt}=f_2\left({\hat{x}}_1,\hat{y},u,w\right)-v\end{array}\right.---\left(2\right)$

其中和分别是x₁和y的估计；L是待求增益，v是估计值与实际值之间的误差，u为输入的阳极燃料流量和阴极燃料流量，W为电堆的输出电流；

为求解增益L，先将模型进行线性化成如下形式：

$\left\{\begin{array}{c}\frac{{\mathrm{dx}}_1}{dt}=A_{11}{x}_1+A_{12}y+B_{1}u+C_{1}w\\ \frac{dy}{dt}=A_{21}{x}_1+A_{22}y+B_{2}u+C_{2}w\end{array}\right.---\left(3\right)$

对比非线性模型，验证线性化后模型的准确度。然后，根据SOFC气、固体间热传递响应速度差异很大的特性，对模型进行如下分解：

气体热传递模型：

$\left\{\begin{array}{c}\frac{{\mathrm{dx}}_{air}^1}{dt}=A_{air}^{11}{x}_1+A_{air}^{12}y+B_{air}^1{u}_{air}+C_{air}^{1}w\\ \frac{{\mathrm{dy}}_{air}}{dt}=A_{air}^{21}{x}_1+A_{air}^{22}y+B_{air}^2{u}_{air}+C_{air}^{2}w\end{array}\right.---\left(4\right)$

固体热传递模型：

$\left\{\begin{array}{c}\frac{{\mathrm{dx}}_{sol}^1}{dt}=A_{sol}^{11}{x}_1+A_{sol}^{12}y+B_{sol}^1{u}_{sol}+C_{sol}^{1}w\\ \frac{{\mathrm{dy}}_{sol}}{dt}=A_{sol}^{21}{x}_1+A_{sol}^{22}y+B_{sol}^2{u}_{sol}+C_{sol}^{2}w\end{array}\right.---\left(5\right)$

对以上的气体热传递模型和固体热传递模型，分别用线性理论设计观测器，求解出各自对应的观测器增益L₁和L₂。然后，我们用如下方法构造出非线性观测器增益矩阵L：

$L={\left[\begin{array}{cccccccc}{L}_1\left(1,1\right)& 0& L_1\left(2,1\right)& 0& L_1\left(3,1\right)& 0& L_1\left(4,1\right)& 0\\ 0& L_2\left(1,1\right)& 0& L_2\left(2,1\right)& 0& L_2\left(3,1\right)& 0& L_2\left(4,1\right)\end{array}\right]}^T---\left(6\right)$

其中：

[L₁(1,1)L₁(2,1)L₁(3,1)L₁(4,1)]^T＝L₁

[L₂(1,1)L₂(2,1)L₂(3,1)L₂(4,1)]^T＝L₂。

采用上述进一步的有益效果是，_SOFC的电特性由电化学反应决定，响应速度在毫秒级；SOFC的热特性由电堆之间的热传递决定，响应速度在秒甚至分钟级。另，由于气体和固体热传递系数的差异，气、固体之间热传递响应速度也存在时间量级上的差距。针对不同负载变换时的切换策略，需要在已有电堆内部观测器设计方法基础上、进行与负载跟踪过程中时间尺度相匹配的观测器设计优化。根据已有的“全可测最小维状态空间”方法，可以确定出低成本、快速、高精度的观测变量选择方案然后，基于该变量组合，为了进一步实现全负载跟踪过程电堆温度分布快速观测的需要，拟基于线性化状态空间方程设计降维龙伯格观测器。

与直接进行系统极点配置求解L相比，该构造方法有效的避免了因气体和固体热传递响应速度的数量级不同；且该方法可以扩展在不同的操作点进行，获取在该工作点的观测器。

进一步，所述步骤S4中，计算得出系列趋近律函数的具体方法是：

根据所述效率优化函数得到系统全局的系列最优滑模面；

根据不同滑模面之间的功率区间的功率与温度、效率特性确定趋近率函数；

根据功率的最优切换函数确定趋近率函数的分数阶阶次减少抖震；

用改进的恒权重计算方法求解趋近率，获得系统功率跟踪切换过程中的系统输入参数。

所述全局的最优函数是局部最优稳态操作函数遍历不同电流的值得到，即由每一个电流对应的局部最优函数组成全局最优函数。

采用上述进一步的有益效果是，通过这种方法可以得到SOFC系统在所有允许输入参数下的输出参数，方便计算得出滑模区间。

进一步，所述SOFC系统不同输入参数组合下的系统状态及输出参数包括：电堆内部温度分布参数、系统效率、系统功率。

采用上述进一步的有益效果是，SOFC系统安全运行必须保证电堆内部温度分布均匀且不超过温度约束，较大的温度梯度或者较高的温度将会对SOFC电堆造成严重的不可逆转的损害，影响SOFC系统寿命和运行安全。为了让SOFC系统安全，高效，长时间的运行，通过分数阶滑模变结构热电协同控制保证功率输出和内部温度参数要求的同时消耗较少的燃料，达到高效率运行的目标。

附图说明

图1为本发明的步骤流程图。

具体实施方式

以下结合附图对本发明的原理和特征进行描述，所举实例只用于解释本发明，并非用于限定本发明的范围。

如图1所示，基于分数阶滑模变结构SOFC系统热电协同控制方法，包括以下步骤：

S1、采集SOFC系统不同输入参数组合下的系统状态及输出参数，通过系统辨识的方法得到全负载区间稳态功率与温度、效率响应特性、旁路阀开度BP对效率优化的影响函数和指定负载切换区间及延时条件下的效率优化函数；

S2、基于所述全负载区间稳态功率与温度、效率响应特性、旁路阀开度BP对效率优化的影响函数和指定负载切换区间及延时条件下的效率优化函数获取系统的局部最优稳态操作函数、拓展形成稳态下全局优化函数和不同切换负载区间不同延时条件下的功率跟踪函数；

S3、根据所述最优稳态操作函数、拓展形成稳态下全局优化函数和指定负载切换区间及延时条件下的效率优化函数计算得出滑模区间；

S4、根据所述滑模区间、全负载区间稳态功率与效率相应特性、指定负载切换区间及延时条件下的效率优化函数和不同切换负载区间不同延时条件下的改进优化函数计算得出系列趋近律函数；

S5、通过分数阶优化法对所述系列趋近律函数消除抖震，通过计算求解出趋近律。

所述步骤S1中，通过运行电堆的分布式节点模型获取电堆内部温度分布参数。

所述电堆的维数为一维。

获得温度分布的方法可以理解为：所建立的电堆模型实际上是将电堆分为五个节点，即所说的分布式模型，每个节点都有温度参数，运行这个模型就可以得到对应节点的温度，这是一组温度参数，它可以反应电堆内部温度的分布。

所述步骤S1中局部最优稳态操作函数的确定方法,包括以下步骤：

固定工作电流Is和旁路阀开度BP，得到以燃料利用率FU为横坐标、以空气过量比AR为纵坐标的输出性能图谱；

对于所述输出性能图谱，固定工作电流Is，得到不同旁路阀开度下以燃料利用率为横坐标、以空气过量比为纵坐标的输出性能图谱；

遍历不同工作电流Is的值，获取任一操作点(Is，BP，AR，FU)下的系统输出性能；

面向负载跟踪，针对任意功率的稳态输出需求，确定与该需求对应的最优稳态操作函数。

所述步骤S1中电堆内部温度分布参数是通过电堆内部温度梯度观测器来获取的，所述电堆内部温度梯度观测器的构建包括以下步骤：

建立基于线性化状态空间方程涉及降维龙伯格观测器；

假定电化学反应是瞬间完成的，进行准静态假设后，模型为：

$\left\{\begin{array}{c}\frac{{\mathrm{dx}}_1}{dt}=f_1(x_1,y,u,w)\\ \frac{dy}{dt}=f_2\left(x_1,y,u,w\right)\end{array}\right.---\left(1\right)$

其中x1为前四个节点的固体层温度和空气层温度，y为第五个节点的固体层温度和空气层温度；

w＝I_tot，

对以上模型，设计如下形式的观测器：

$\left\{\begin{array}{c}\frac{d{\hat{x}}_1}{dt}=f_1\left({\hat{x}}_1,\hat{y},u,w\right)+Lv\\ \frac{d\hat{y}}{dt}=f_2\left({\hat{x}}_1,\hat{y},u,w\right)-v\end{array}\right.---\left(2\right)$

其中和分别是x₁和y的估计；L是待求增益，v是估计值与实际值之间的误差，u为输入的阳极燃料流量和阴极燃料流量，W为电堆的输出电流；

为求解增益L，先将模型进行线性化成如下形式：

$\left\{\begin{array}{c}\frac{{\mathrm{dx}}_1}{dt}=A_{11}{x}_1+A_{12}y+B_{1}u+C_{1}w\\ \frac{dy}{dt}=A_{21}{x}_1+A_{22}y+B_{2}u+C_{2}w\end{array}\right.---\left(3\right)$

对比非线性模型，验证线性化后模型的准确度。然后，根据SOFC气、固体间热传递响应速度差异很大的特性，对模型进行如下分解：

气体热传递模型：

$\left\{\begin{array}{c}\frac{{\mathrm{dx}}_{air}^1}{dt}=A_{air}^{11}{x}_1+A_{air}^{12}y+B_{air}^1{u}_{air}+C_{air}^{1}w\\ \frac{{\mathrm{dy}}_{air}}{dt}=A_{air}^{21}{x}_1+A_{air}^{22}y+B_{air}^2{u}_{air}+C_{air}^{2}w\end{array}\right.---\left(4\right)$

固体热传递模型：

$\left\{\begin{array}{c}\frac{{\mathrm{dx}}_{sol}^1}{dt}=A_{sol}^{11}{x}_1+A_{sol}^{12}y+B_{sol}^1{u}_{sol}+C_{sol}^{1}w\\ \frac{{\mathrm{dy}}_{sol}}{dt}=A_{sol}^{21}{x}_1+A_{sol}^{22}y+B_{sol}^2{u}_{sol}+C_{sol}^{2}w\end{array}\right.---\left(5\right)$

对以上的气体热传递模型和固体热传递模型，分别用线性理论设计观测器，求解出各自对应的观测器增益L₁和L₂。然后，我们用如下方法构造出非线性观测器增益矩阵L：

$L={\left[\begin{array}{cccccccc}{L}_1\left(1,1\right)& 0& L_1\left(2,1\right)& 0& L_1\left(3,1\right)& 0& L_1\left(4,1\right)& 0\\ 0& L_2\left(1,1\right)& 0& L_2\left(2,1\right)& 0& L_2\left(3,1\right)& 0& L_2\left(4,1\right)\end{array}\right]}^T---\left(6\right)$

其中：

[L₁(1,1)L₁(2,1)L₁(3,1)L₁(4,1)]^T＝L₁

[L₂(1,1)L₂(2,1)L₂(3,1)L₂(4,1)]^T＝L₂。

所述步骤S4中，计算得出系列趋近律函数的具体方法是：

根据所述效率优化函数得到系统全局的系列最优滑模面；

根据不同滑模面之间的功率区间的功率与温度、效率特性确定趋近率函数；

根据功率的最优切换函数确定趋近率函数的分数阶阶次减少抖震；

用改进的恒权重计算方法求解趋近率，获得系统功率跟踪切换过程中的系统输入参数。

所述全局的最优函数是局部最优稳态操作函数遍历不同电流的值得到，即由每一个电流对应的局部最优函数组成全局最优函数。

所述SOFC系统不同输入参数组合下的系统状态及输出参数包括：电堆内部温度分布参数、系统效率、系统功率。

以上所述仅为本发明的较佳实施例，并不用以限制本发明，凡在本发明的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。

Claims (8)

1.基于分数阶滑模变结构SOFC系统热电协同控制方法，其特征在于，包括以下步骤：

S1、采集SOFC系统不同输入参数组合下的系统状态及输出参数，通过系统辨识的方法得到全负载区间稳态功率与温度、效率响应特性、旁路阀开度BP对效率优化的影响函数和指定负载切换区间及延时条件下的效率优化函数；

S2、基于所述全负载区间稳态功率与温度、效率响应特性、旁路阀开度BP对效率优化的影响函数和指定负载切换区间及延时条件下的效率优化函数获取系统的局部最优稳态操作函数、拓展形成稳态下全局优化函数和不同切换负载区间不同延时条件下的功率跟踪函数；

S3、根据所述最优稳态操作函数、拓展形成稳态下全局优化函数和指定负载切换区间及延时条件下的效率优化函数计算得出滑模区间；

S4、根据所述滑模区间、全负载区间稳态功率与效率相应特性、指定负载切换区间及延时条件下的效率优化函数和不同切换负载区间不同延时条件下的改进优化函数计算得出系列趋近律函数；

S5、通过分数阶优化法对所述系列趋近律函数消除抖震，通过计算求解出趋近律。

2.根据权利要求1所述基于分数阶滑模变结构SOFC系统热电协同控制方法，其特征在于，所述步骤S1中，通过运行电堆的分布式节点模型获取电堆内部温度分布参数。

3.根据权利要求2所述基于分数阶滑模变结构SOFC系统热电协同控制方法，其特征在于，所述电堆的维数为一维。

4.根据权利要求1所述基于分数阶滑模变结构SOFC系统热电协同控制方法，其特征在于，所述步骤S1中局部最优稳态操作函数的确定方法,包括以下步骤：

固定工作电流Is和旁路阀开度BP，得到以燃料利用率FU为横坐标、以空气过量比AR为纵坐标的输出性能图谱；

对于所述输出性能图谱，固定工作电流Is，得到不同旁路阀开度下以燃料利用率为横坐标、以空气过量比为纵坐标的输出性能图谱；

遍历不同工作电流Is的值，获取任一操作点(Is，BP，AR，FU)下的系统输出性能；

面向负载跟踪，针对任意功率的稳态输出需求，确定与该需求对应的最优稳态操作函数。

5.根据权利要求1所述基于分数阶滑模变结构SOFC系统热电协同控制方法，其特征在于，所述步骤S1中电堆内部温度分布参数是通过电堆内部温度梯度观测器来获取的，所述电堆内部温度梯度观测器的构建包括以下步骤：

建立基于线性化状态空间方程涉及降维龙伯格观测器；

假定电化学反应是瞬间完成的，进行准静态假设后，模型为：

$\left\{\begin{array}{c}\frac{{\mathrm{dx}}_1}{dt}=f_1(x_1,y,u,w)\\ \frac{dy}{dt}=f_2(x_1,y,u,w)\end{array}\right.---\left(1\right)$

其中为前四个节点的固体层温度和空气层温度，y为第五个节点的固体层温度和空气层温度；

$u={\left[\begin{array}{cc}{F}_{ca}^{in}& F_{an}^{in}\end{array}\right]}^T,w=I_{tot},$

对以上模型，设计如下形式的观测器：

$\left\{\begin{array}{c}\frac{d{\hat{x}}_1}{dt}=f_1({\hat{x}}_1,\hat{y},u,w)+Lv\\ \frac{d\hat{y}}{dt}=f_2({\hat{x}}_1,\hat{y},u,w)-v\end{array}\right.---\left(2\right)$

其中和分别是x₁和y的估计；L是待求增益，v是估计值与实际值之间的误差，u为输入的阳极燃料流量和阴极燃料流量，W为电堆的输出电流；

为求解增益L，先将模型进行线性化成如下形式：

$\left\{\begin{array}{c}\frac{{\mathrm{dx}}_1}{dt}=A_{11}{x}_1+A_{12}y+B_{1}u+C_{1}w\\ \frac{dy}{dt}=A_{21}{x}_1+A_{22}y+B_{2}u+C_{2}w\end{array}\right.---\left(3\right)$

对比非线性模型，验证线性化后模型的准确度，然后根据SOFC气、固体间热传递响应速度差异很大的特性，对模型进行如下分解：

气体热传递模型：

$\left\{\begin{array}{c}\frac{{\mathrm{dx}}_{air}^1}{dt}=A_{air}^{11}{x}_1+A_{air}^{12}y+B_{air}^1{u}_{air}+C_{air}^{1}w\\ \frac{{\mathrm{dy}}_{air}}{dt}=A_{air}^{21}{x}_1+A_{air}^{22}y+B_{air}^2{u}_{air}+C_{air}^{2}w\end{array}\right.---\left(4\right)$

固体热传递模型：

$\left\{\begin{array}{c}\frac{{\mathrm{dx}}_{sol}^1}{dt}=A_{sol}^{11}{x}_1+A_{sol}^{12}y+B_{sol}^1{u}_{sol}+C_{sol}^{1}w\\ \frac{{\mathrm{dy}}_{sol}}{dt}=A_{sol}^{21}{x}_1+A_{sol}^{22}y+B_{sol}^2{u}_{sol}+C_{sol}^{2}w\end{array}\right.---\left(5\right)$

对以上的气体热传递模型和固体热传递模型，分别用线性理论设计观测器，求解出各自对应的观测器增益L₁和L₂，然后用如下方法构造出非线性观测器增益矩阵L：

$L={\left[\begin{array}{cccccccc}{L}_1(1,1)& 0& L_1(2,1)& 0& L_1(3,1)& 0& L_1(4,1)& 0\\ 0& L_2(1,1)& 0& L_2(2,1)& 0& L_2(3,1)& 0& L_2(4,1)\end{array}\right]}^T---\left(6\right)$

其中：

[L₁(1,1)L₁(2,1)L₁(3,1)L₁(4,1)]^T＝L₁

[L₂(1,1)L₂(2,1)L₂(3,1)L₂(4,1)]^T＝L₂。

6.根据权利要求1所述基于分数阶滑模变结构SOFC系统热电协同控制方法，其特征在于，所述步骤S4中，计算得出系列趋近律函数的具体方法是：

根据所述效率优化函数得到系统全局的系列最优滑模面；

根据不同滑模面之间的功率区间的功率与温度、效率特性确定趋近率函数；

根据功率的最优切换函数确定趋近率函数的分数阶阶次减少抖震；

用改进的恒权重计算方法求解趋近率，获得系统功率跟踪切换过程中的系统输入参数。

7.根据权利要求6所述基于分数阶滑模变结构SOFC系统热电协同控制方法，其特征在于，所述全局的最优函数是局部最优稳态操作函数遍历不同电流的值得到，即由每一个电流对应的局部最优函数组成全局最优函数。

8.根据权利要求1-7任一项所述基于分数阶滑模变结构SOFC系统热电协同控制方法，其特征在于，所述SOFC系统不同输入参数组合下的系统状态及输出参数包括：电堆内部温度分布参数、系统效率、系统功率。