Skip to main content
QUICK REVIEW

[论文解读] Optimal Control of a Parabolic Distributed Parameter System Using a Barycentric Shifted Gegenbauer Pseudospectral Method

Kareem T. Elgindy|arXiv (Cornell University)|Mar 4, 2016
Advanced Numerical Methods in Computational Mathematics参考文献 2被引用 3
一句话总结

本文提出了一种新颖的重心平移勒让德-盖根鲍尔伪谱方法(BSGPM),用于求解由抛物型分布参数系统控制的最优控制问题。通过将问题重新表述为积分形式,并利用平移勒让德-盖根鲍尔-高斯积分节点的稳定重心插值,该方法在使用最少的配点数时实现了指数收敛,且精度极高,优于现有方法在效率和精度方面的表现。

ABSTRACT

In this paper, we introduce a novel pseudospectral method for the numerical solution of optimal control problems governed by a parabolic distributed parameter system. The infinite-dimensional optimal control problem is reduced into a finite-dimensional nonlinear programming problem through shifted Gegenbauer quadratures constructed using a stable barycentric representation of Lagrange interpolating polynomials and explicit barycentric weights for the shifted Gegenbauer-Gauss (SGG) points. A rigorous error analysis of the method is presented, and a numerical test example is given to show the accuracy and efficiency of the proposed pseudospectral method.

研究动机与目标

  • 开发一种高阶、稳定的数值方法,用于求解由抛物型偏微分方程控制的最优控制问题。
  • 通过积分重构和谱离散化,将无限维最优控制问题转化为条件良好的有限维非线性规划问题。
  • 通过稳定的重心插值和移位勒让德-盖根鲍尔积分,仅使用最少的配点实现指数收敛速率。
  • 相较于径向基函数等现有方法,提供更高的精度和计算效率。
  • 为一大类PDE约束的最优控制问题提供一个稳健的框架。

提出的方法

  • 该方法通过状态变量和控制变量的连续积分,将原始PDE约束的最优控制问题重构为积分形式。
  • 采用移位勒让德-盖根鲍尔(SGG)点作为配点,利用具有显式重心权重的拉格朗日插值多项式的稳定重心表示。
  • 通过移位勒让德-盖根鲍尔积分构造积分矩阵,以精确且高效地近似积分算子。
  • 将无限维问题离散化为具有线性约束的有限维非线性规划问题,可由标准优化求解器求解。
  • 通过利用勒让德-盖根鲍尔多项式的谱性质以及重心插值的良好条件性,确保高阶精度。
  • 进行了严格的误差分析,证明了随着配点数增加,误差呈谱衰减。

实验结果

研究问题

  • RQ1重心平移勒让德-盖根鲍尔伪谱方法能否在抛物型最优控制问题中实现指数收敛?
  • RQ2移位勒让德-盖根鲍尔多项式中形状参数α的选择如何影响解的精度和稳定性?
  • RQ3对于同一问题,所提方法能否在精度和计算成本方面优于径向基函数方法?
  • RQ4随着配点数的增加,该方法的收敛行为如何?
  • RQ5在离散化形式中,该方法在多大程度上保持了初始条件和边界条件等约束?

主要发现

  • BSGPM实现了指数收敛速率,随着配点数的增加,误差呈谱衰减。
  • 仅在每个空间和时间方向使用N=5个配点,该方法的精度即可与每个方向使用120个节点的径向基函数方法相媲美。
  • 对于较小的N值,SGG点处非正α值(如α = -0.2)的离散化在初始条件(ψ₁)的最大误差最小,表明在小规模问题中具有更优的精度。
  • 所有测试配置下,解的可行性以及最大边界条件违反量(ψ₂)均接近机器精度,表明约束满足性极强。
  • 所有测试的N和α值下,近似最优性能指标Jₙ,ₙ*均稳定在约15,证实了数值的一致性。
  • 即使在N=12且α=-0.2时,该方法仍能生成高度精确的状态和控制轨迹,如图2和图3所示,表明仅用极少的网格点即可实现高分辨率逼近。

更好的研究,从现在开始

从阅读论文到最终审阅,大幅缩短您的研究时间。

无需绑定信用卡

本解读由 AI 生成,并经人工编辑审核。