1.一种固体氧化物燃料电池的量子模型预测控制方法，其特征在于，该方法将量子计算和模型预测控制应用到固体氧化物燃料电池系统控制中，实现固体氧化物燃料电池系统的快速负载跟踪；固体氧化物燃料电池两极发生的电化学反应方程为：+ ‑

其中，H2为氢气；H为氢正离子；e为电子；O2为氧气；H2O为水蒸气；固体氧化物燃料电池的输出电压VSOFC表达式为：

VSOFC＝ENernst‑Vactivation‑Vohmic‑Vconcentration        (2)其中，VSOFC为固体氧化物燃料电池的输出电压；ENernst为固体氧化物燃料电池的开路电压；Vactivation为固体氧化物燃料电池的反应物活化的损耗电压降；Vohmic为固体氧化物燃料电池的内部欧姆电阻引起的损耗电压降；Vconcentration为固体氧化物燃料电池发生化学反应后由反应物浓度降低引起的损耗电压降；固体氧化物燃料电池的开路电压ENernst表示为：其中，ENernst,0为固体氧化物燃料电池的参考电动势；kE为固体氧化物燃料电池的经验常数；TSOFC为固体氧化物燃料电池的工作温度；RSOFC为固体氧化物燃料电池的内部电阻；F为法拉第常数； 为固体氧化物燃料电池发生化学反应的氢气的实际压强； 为固体氧化物燃料电池发生化学反应的氧气的实际压强； 为固体氧化物燃料电池发生化学反应的水蒸气的实际压强；ln()为以常数e为底数的对数；固体氧化物燃料电池的反应物活化的损耗电压降Vactivation为：Vactivation＝a+bln(ISOFC)             (4)其中，a和b为塔菲尔参数；ISOFC为固体氧化物燃料电池的输出电流；固体氧化物燃料电池的内部欧姆电阻引起的损耗电压降Vohmic由以下计算：其中，aelectrolyte和belectrolyte为固体氧化物燃料电池的电解质电阻常数；exp为常数e的指数函数；ASOFC为固体氧化物燃料电池反应区的面积；δelectrolyte为固体氧化物燃料电池的电解质厚度； 与 为多个固体氧化物燃料电池互连的材料电阻； 为多个固体氧化物燃料电池之间的互连厚度固体氧化物燃料电池发生化学反应后由反应物浓度降低引起的损耗电压降Vconcentration通过以下计算：其中， 为固体氧化物燃料电池的供气管道内氢气的压强； 为固体氧化物燃料电池的供气管道内水蒸气的压强；模型预测控制方法的优化问题为：其中，J为模型预测控制中滚动优化的最优解；N为模型预测控制中的预测次数；y(k+i)T

为第k+i时刻的固体氧化物燃料电池的输出，且y(k+i)＝[VSOFC(k+i)Fu(k+i)] ；VSOFC(k+i)为第k+i时刻固体氧化物燃料电池的输出电压；Fu(k+i)为第k+i时刻固体氧化物燃料电池的输出燃料；yref为系统给定的参考输出； u(k+i‑1)为第k+i‑1时刻的固体氧化物燃料电池的输入，且 为第k+i‑1时刻的固体氧化物燃料电池的氢气输入量；Mair(k+i‑1)第k+i‑1时刻的固体氧化物燃料电池的氧气输入量；μref为与参考输出相对应的参考输入； 与参考输出相对应的参考输入μref表示为：

‑1 ‑1

μref＝(C(I‑A) B)             (8)其中，

模型预测控制第k+i时

刻的输出y(k+i)表示为：

y(k+i)＝Cx(k+i)             (9)T

其中，x(k)＝[VSOFC(k) Fu(k) ISOFC(k)] ，ISOFC(k)为第k时刻的固体氧化物燃料电池输出的电流；x(k+i)为第k+i时刻的固体氧化物燃料电池系统状态；量子模型预测控制第k+1时刻的预测状态 为：

其中， 为量子模型预测控制第k+1时刻固体氧化物燃料电池系统的预测状态；

为量子模型预测控制第k时刻的固体氧化物燃料电池系统状态；u(k)为量子模型预测控制第k时刻的固体氧化物燃料电池系统输入；y(k)为量子模型预测控制第k时刻的固体氧化物燃料电池系统输出； 三维量子旋转门η为：其中，为量子旋转角；θ为量子旋转角；固体氧化物燃料电池系统的输入约束条件为：umin≤u(k+i)≤umax            (12)其中，umin为固体氧化物燃料电池系统的输入下限；u(k+i)为第k+i时刻的固体氧化物燃料电池的输入；umax为固体氧化物燃料电池系统的输入上限；固体氧化物燃料电池系统第k+i时刻的输出约束条件为：

ymin≤y(k+i)≤ymax         (13)其中，ymin为固体氧化物燃料电池系统的输出下限；ymax为固体氧化物燃料电池系统的输出上限。