第22卷第l期2004年3月飞行力学FLIGHTDYNAMICSVOI．22NO．1Mar．2004文章编号：1002—0853(2004)01—0014-04改进直接多重打靶算法及其应用胡朝江，陈士橹(1．北京航空工程技术研究所11室．北京100076；2．西北工业大学航天工程学院．陕西西安710072)摘要：直接多重打靶算法是求解最优控制问题很有效的方法之一。通过仅假设出节点处的控制变量值．使该算法在求解最优控制问题时更方便，收敛更快。利用改进算法求解了飞机在风切变中着陆和敏捷性管理系统的优化设计问题。结果表明．改进算法能较好地求解一类受约束最优控制问题，对奇异最优控制问题也能较好求解。关键词：直接多重打靶算法；最优控制；风切变；敏捷性管理系统中图分类号：V212．1文献标识码：A引言对飞行力学逆问题的求解一直是飞行力学研究的一个重要领域之一，飞行力学逆问题通常都可视为一个最优控制问题。长期以来，人们已经为这类问题的求解探索了一系列行之有效的方法，归纳起来可分为直接求解法和间接求解法两大类。相比之下，直接法无需象间接法那样，把最优控制问题转化为比较麻烦的两点边值问题，而是将时间历程分段，使状态和控制参数化，构成参数最优控制问题直接求解。直接多重打靶算法则是直接法中比较有名的，该算法能较好地求解一类飞行力学逆问题Ll]。但从适时控制的角度而言，除提高计算机硬件的性能外，也希望改进计算方法。从这个角度考虑，通过仅假设出节点处的控制变量值，对该算法作了适当的改进，从而使该算法所面临的非线性规划子问题维数降低很多，因而使得该算法在求解具体问题上更方便，更节省时间，更有利于适时最优控制问题的实现。1方法简介最优控制问题的形式主要有Lagrange，Mayer及Bolza3种，这3种形式在数学上是等价的。为不失一般性，考虑到本文将要求解的问题，设最优控制问题的形式如下：广7’J：==minI厂f，(，)，lI(f)]dtⅣ‘，)JO式中，状态变量(￡)及控制变量“(f)满足常微分方程dx／dt=f[t，(f)，“(f)]。控制变量满足上下限约束，。≤lI(f)≤g，状态变量满足上下限约束，≤(f)≤及边界约束(0)=。，(丁)~--X7，。用直接多重打靶法求解最优控制问题通常的做法是把时间历程分为，段，通过假设出节点处的状态变量、控制变量及时间，把微分方程约束化为，个初值问题，然后求解使节点处的匹配条件得到满足的控制规律lI(f)，使所取泛函指标最小。由于采用把节点处的状态变量、控制变量及时间都同时假设出来的做法，使得所转化而成的非线性规划问题变量个数会很多，从而导致求解比较困难。如设所要求解的最优控制问题的阶数为4阶，控制变量为1个，时间历程被l0等分，则对应的非线性规划问题的变量将多达50余个，这无疑会使问题变得非常复杂。为此，经过分析认为：仅假设出节点处的控制变量值即可，具体如下：(1)将原问题化为控制终端时间给定的Mayer问题收稿日期：2002—12—09；修订日期：2003一¨一10作者简介：胡朝江(1968一)，男．贵州开阳人，副教授，博士。主要从事飞行动力学及控制方面的研究；陈士橹(1920一)．男．浙江东阳人．俄罗斯宇航科学院外籍院士．中国工程院院士，西北工业大学教授．博士生导师，主要从事飞行动力学及控制方面的研究。维普资讯http://www.cqvip.com第1期胡朝江等．改进直接多重打靶算法及其应用15①引进新的时间变量r∈[0，1]，设终端时间丁为可变参数，定义f=丁r；②引进增广状态变量及状态方程一(，U，f)，-z(0)一0，则目标泛函化为：J—J((1))其中：(1)一[．rl(1)，2(1)，⋯，，r十l(1)]‘将t=Tr代入状态方程及约束式中，问题中的自变量t换为r，积分区间变为Eo，1]，则上述终时不给定的最优控制问题就化为了终时给定的标准Mayer问题。(2)将问题化为有限维非线性规划问题①将时间区间r∈I-o，13,n等分，得到+1个节点；②引入一组向量，作为节点(—o，1，⋯，一1)处控制变量的估计值。节点之间的控制变量值由相邻两点之间的值线性插值求得。若当节点处控制变量(一0，1，⋯，，一1)已知时，则在已知状态变量初值的情况下，各节点处的状态变量可顺次积分求得，因而可求得(1)，并进而求得泛函指标值。因此，可以认为，微分方程的解及泛函指标都仅...
