常微分方程组的数值计算一直是备受人们关注的领域,对此已发展了丰富的数值方法。近年来,精细积分方法得到广泛关注,已扩展到时变、非线性微分方程、偏微分方程的求解,并成功地应用到结构动力响应、随机振动、波导、热传导以及最优控制等领域,为不同领域的数值计算提供了一个高精度、高稳定性的算法平台,值得深入研究。另一方面,控制领域对数值计算的关注度和重要性意识正在加强,而合适的理论框架对于构造高性能算法有重要意义。现代控制论所奠基的状态空间法的起点至少也应回溯到Hamilton正则方程体系,表明经典力学与现代控制论有共同的数学形式和理论基础,两个学科的问题是相互对应的。因此,借鉴力学中成熟的有限元、子结构分析等方法,展开对最优控制领域数值方法和控制系统设计的研究是有意义的。本论文以发展高效、可靠的数值算法为主线,改进了精细积分算法平台的性能,研究了时滞、时变、非线性系统最优控制的数值计算和控制器设计等问题,开发了最优控制系统设计工具箱并将其应用于卫星编队飞行控制的研究。主要工作如下:(1)采用矩阵函数逼近理论,提出了基于Pade级数逼近的矩阵指数精细积分方法中加权参数N和级数展开项数q的递推自适应选择算法,提高了精细积分方法的计算效率。并与MATLAB内置函数expm()进行了比较,表明本文方法在达到相同的效率的同时具有更高的精度和稳定性。(2)提出了动力初值问题中非齐次项产生的Duhamel积分响应矩阵的扩展精细积分方法(EPIM),该方法不需对系统矩阵(或相关动力矩阵)求逆。当非齐次项为多项式函数、指数函数、正/余弦函数及其组合函数的形式时,可以得到计算机意义上的精确解。并推广应用于:1)与虚拟激励法结合,应用于随机振动响应的计算;2)结合传统数值积分技术(如Taylor级数单步法和Adams多步法),构造了求解非线性微分方程的显式/隐式算法;3)利用系数周期性变化的特点,导出了周期时变Floquet转移矩阵和一类非线性周期系统响应的计算格式;等。算例表明,基于扩展精细积分方法构造的算法提高了数值稳定性和适用范围,具有高效、高精度、高稳定性的优点。(3)提出了两点边值问题中非齐次项产生的区段响应矩阵的扩展精细积分方法(EPIM),当非齐次项为多项式函数、指数函数、正/余弦函数及其组合函数的形式时,可以得到计算机意义上的精确解。在此基础上,研究了一般非齐次项的处理方法以及在无限长区段和变系数两点边值问题中的应用。还结合周期时变Floquet转移矩阵的扩展精细积分方法,导出了周期变系数Riccati、Lyapunov、Sylvester等矩阵微分方程的保结构算法,数值算例验证了算法的有效性。(4)对时滞系统的H_∞最优控制和滤波进行了研究。首先采用扩展精细积分方法对连续时滞系统方程和性能指标离散化,以最大程度地保证与原系统的等价性。然后引入合适的增维向量,化为不显含时滞的标准离散形式,采用区段混合能方法和扩展W-W算法进行计算分析,增强了增维方法的可行性,从而为时滞H_∞最优控制和滤波系统的分析和设计提供了一套精确、稳定的算法。并导出了含输入时滞的H_∞全信息控制器,应用于建筑结构的减振控制,仿真显示对于不同的时滞量和地震激励形式,结构的振动响应都得到了有效抑制,验证了控制器的有效性。(5)时变、非线性最优控制系统设计导出Hamilton系统两点边值问题,其数值算法应该保辛。本文在区段分析的框架下,提出了时变Hamilton两点边值问题基于常值精细积分的保辛摄动方法,导出了零阶、摄动系统分别基于区段混合能矩阵和区段传递矩阵的组合公式以及对应关系,指出前者具有内在的稳定性从而是更好的选择。进一步提出了时变非齐次Hamilton两点边值问题的保辛摄动方法,并应用于非线性最优控制问题的迭代计算,结果表明,迭代过程中关键算法的改进显著地提高了收敛速度,降低了对初始迭代值的敏感性,说明保辛摄动方法是一种高精度和稳定的算法。(6)传统终端控制器往往存在终端高增益或奇异现象,只好在靠近终端区段采用开环控制。本文引入终端“软约束项”改进了性能指标,并利用Lagrange乘子的常数本质,构造了非奇异的、两个区段都具有反馈-前馈控制结构的终端控制器。分析了引入的“软约束项”对构造反馈结构控制器的重要影响,对于最小能量控制问题尤为重要。进一步利用区段混合能矩阵构造了反馈增益矩阵和控制系统方程的闭合解,导出了保结构递推算法,方便了控制器的设计与实现。并将该方法推广应用于离散时间系统的终端控制器设计。(7)针对当前主流商业控制系统设计软件MATLAB缺乏有限长时间时变最优控制器设计功能的现状开发了PIMCSD Toolbox;在此基础上研究了典型双星编队重构的时变最优控制方案,研究成果为航天器编队控制系统的工程设计和应用提供了重要参考。