基于样条的高阶矩阵微分方程近似解

矩阵微分方程经常出现在许多物理模型和工程技术模型中.利用矩阵样条构造形如{y(p)(x)=Ap-1(x)y(p-1)(x)+Ap-2(x)y(p-2)(x)+…+A1(x)y(1)(x)+A0(x)y(x)+B0(x),y(a)=ya,…,y(p-1)(a)=y(p-1)a,x∈[a,b];Ai(x),B0(x)∈C4[a,b],0≤i≤p-烅烄烆1的高阶矩阵线性微分方程初值问题的数值解.给出实现算法和数值解的近似误差估计以及数值实例.先将高阶矩阵微分方程转化为一阶矩阵微分方程,然后利用三次矩阵样条求出一阶矩阵线性微分方程的数值解,从而解决高阶微分方程问题.

On the Matrix Equations of Higher Degrees

A generalm atrix equation ∑m s= 0 ∑s t= 0 XtAstXs- t= 0 in unknow n X ofdegree m is con- sidered. The solvability of the equation w ithin som e restricted subsets of n×n m atrices is dis- cussed. The m ethods offinding the set ofsolutions in the subsetfor the equation are also provid- ed.