基于四次矩阵样条的矩阵微分方程近似解.docxVIP

基于四次矩阵样条的矩阵微分方程近似解 第30卷第11期 2007年11月 合肥工业大学(自然科学版) JOURNALOFHEFEIUNIVERSITYOFTECHNOLOGY Vol.3ONo.ll Nov.2007 基手四次矩阵样条的矩阵微分方程近似解 郭清伟，王芳 (合SLT.业大学理学院，安徽合肥230009) 摘要:矩阵微分方程经常出现在物理模型和工程技术模型中.文章给 出了用四次矩阵样条构造形如1,,= A(x)Y+.B⑵,Y(O)=Yo,z G ,hi,A(z),?B⑵丘,6］的一阶矩阵线性微分方 程初值问题近似解的方 法，研究了该方法的逼近误差并编制了实现该方法的一个算法，最后 给出一些数值实例;比较结果表明，用四 次矩阵样条所构造的近似解的逼近效果要比用三次矩阵样条所构 造的近似解的逼近效果好. 关键词:一阶矩阵线性微分方程;三次矩阵样条;四次矩阵样条;逼近 中图分类号:0221.73文献标识码：A文章编 号:1003—5060(2007) 11—1537—05 Approximatesolutionofmatrixdifferential equationsbythequarticmatrixspline GUOQing-wei.WANGFang9 (SchoolofSciences,HefeiUniversityofTechnology,Hefei230009,China) Abstract:Matrixdifferentialequationsappearfrequentlyinawidevarietyof modelsinphysicsandan一 gineering.Thispaperdealswiththeconstructionoftheapproximatesolutiono ffirst——ordermatrixlin ? eardifferentialequationsgivenby=A(x)Y+B(z)andy(O) 一 II0,wherez 三 a,6],andA(z), B(z) e a,6],usingthequarticmatrixspline(QMS).Anestimationoftheapproximation error andanalgorithmfortheimplementationoftheconstructionmethodaredescri bed.Twoillustrative examplesareincluded,andforthesamefirst― rdermatrixlineardifferential equation,theresultsob— tainedbyusingtheQMSarecomparedwiththosebyusingthecubicmatrixspli ne,whichshowsthat theerrorofapproximationbyusingtheQMSissmallerthanthatbyusingthecu bicmatrixspline. Keywords :first-ordermatrixlineardifferentialequation;cubicmatrixspline ;quarticmatrix spline;approximation 矩阵微分方程经常出现在许多物理模型和工 程技术模型中Ei?4］?除了许多数学模型中的问题 经常写成矩阵的形式外，许多特殊的方法为了解 决数量或向量问题也经常以矩阵形式出现?例 如:用嵌入法求解具有边值条件的常微分方程的 问题［5;用打靶法求解具有边值条件的标量和向 量微分方程的问题［6;用线性方法研究偏微分方 程的数值积分问题［7］;用同伦方法求解非线性方 程组的问题_8］等. 向量化方法将一个矩阵方程化为几个孤立的 标量方程或向量方程有许多缺点［?用向量化方 法会使方程中的量失去物理意义;计算的成本会 增加很多;向量化方法不具有利用一些符号语言 处理矩阵表达式的优点. 本文在原有三次矩阵样条算法的基础上，利 用四次矩阵样条函数构造一阶矩阵线性微分方程 的近似解，并给出实现算法和近似解的逼近误差 估计以及一些数值实例?主要目的在于验证用三 次矩阵样条和四次矩阵样条分别构造一阶矩阵线 性微分方程的近似解时,如果所取的步长相同，则 四次矩阵样条所构造的近似解的逼近效果会不会 比三次矩阵样条所构造的近似解的逼近效果要 收稿日期:2006—11—10;修改日期:2007—05—10 基金项目:合肥工业大学基金资助项目(061007F) 作者简介:郭清伟(1968-),男，安徽滩溪人，合肥工业大学副教授，硕士 生导师 1538合肥工业大学(自然科学版)第30卷 好?即在相同的计