李雅普诺夫方程AX-XB=C的误差界

本文对李雅普诺夫方程AX-XB=C及(A+δA)(X+δX)-(X+δX)(B+δB)=C+δC的解X及X+δX给出了相对误差‖δx‖/‖X‖的一个上界。     

对称矩阵反问题解的稳定性

考虑如下对称矩阵反问题：问题Ⅰ：给定Ｘ，Ｂ∈Ｒ（ｎ×ｎ），求Ａ∈ＳＲ（ｎ×ｎ）使得ＡＸ＝Ｂ，其中ＳＲ（ｎ×ｎ）＝｛｛Ａ｜Ａ∈Ｒ（ｎ×ｎ），ＡＴ＝Ａ｝。问题Ⅱ：给定，求使得其中是矩阵的Ｆｒｏｂｅｎｉｕｓ范数，ＳＡ表示问题Ⅰ的解集合。本文讨论问题Ⅰ解的稳定性，给出问题Ⅱ解的扰动界。     

Lyapunov矩阵方程的一种数值解法及部分并行处理

本文设计了求解Lyapunov矩阵方程的一种新方法。所考虑的矩阵方程是 AX—XB=C(1)其中A,B,C分别是m×m,n×n和m×n的已知矩阵。 该方法首先是将系数矩阵A,B初等相似约化为三对角矩阵,即存在可逆矩阵U,V,使U~(-1)AU=A,V~(-1)BV=B,其中A,B为三对角矩阵。然后设计了矩阵方程AY—YB=C的公式解法,分三步: 1)求f(λ)=det(λI—A)的λ各次幂的系数a_0,…,a_m; 2)计算sum from i=1 to m (A_(m-i)-CB~(m-i)),f(B); 3)求解Y。解方程AY—YB=C的方法称为THR算法。 最后经逆变换获得原矩阵方程(1)的解X。 求解矩阵方程(1)的方法称为R—THR算法。该方法的计算量约为m~3+4/3n~3+7m~2n+5nm~2+m~2。 本文给出了R—THR的串行计算的数值例子,并给出了THR算法的并行计算格式。最后通过几种数值方法的比较,表明该方法是可行的,也是有效的。     

线性矩阵方程AXB－CXD＝E的误差界

首先讨论线性矩阵方程ＡＸＢ－ＣＸＤ＝Ｅ在有唯一解的条件下方程之解Ｘ的一个上界，再考虑系数矩阵Ａ，Ｂ，Ｃ，Ｄ，Ｅ发生小扰动时，线性矩阵方程（Ａ＋δＡ）Ｘ（Ｂ十δＢ）－（Ｃ＋δＣ）Ｘ（Ｄ＋δＤ）＝Ｅ＋δＥ之解Ｘ作为Ｘ的近似值的相对误差的一个上界。     

第二类对称三对角矩阵逆特征值问题及其应用

本文研究如下问题: 问题Ⅰ 给定n×2实矩阵X和实对角矩阵A=diag(λ_1,λ_2),求第二类n×n实对称三对角矩阵T使得TX=XA。 问题Ⅱ 给定第二类n×n实对称三对角矩阵(?),求第二类对称三对角矩阵(?)使得,其中S_T是问题Ⅰ的解集合。 本文给出了问题Ⅱ有解的充分必要条件,研究了问题Ⅱ解的存在性和唯一性,给出了问题Ⅰ和问题Ⅱ解的表达式,描述了求解问题Ⅰ和问题Ⅱ的数值方法,讨论了数值方法的应用,并给出了一些数值例子。     

箭状矩阵的广义特征值反问题

讨论实对称箭状矩阵（除对角元及最后一行、最后一列元素外，其余位置元素全为零）的广义特征值反问题，它可以用来描述星形弹簧质量系统的振动问题，即给出系统的振动频率如何来确定质点的质量或弹簧的刚度。通过对箭状矩阵特征多项式性质的研究，运用部分分式理论，证明了给定正定箭状矩阵B，实数{λi}i=1^n,{μi}i=1^n-1，满足λ1＜μ1＜…＜μn-1＜λn，存在箭状矩阵A，使广义特征值问题Ax-λBx有解{λi}i=1^n，而广义特征值问题A(n-1)x=λB(n-1)x有解{μi}i=1^n-1，其中A（n-1),B(n-1)分别表示A，B的n-1级主子矩阵。     

带比例矩阵逆特征值问题

研究讨论了一类带比例矩阵的特征值反问题:任意给定2n-1(n≥2)个实数λ(n)1…λ(2)1λ(1)1λ(2)2…λ(n)n,求一个带比例矩阵A,使得λ(j)1和λ(j)j分别是其顺序主子阵Aj(1!j!n)的最小和最大特征值。文中给出了此问题有唯一解的充要条件以及有解的充分条件,并给出了解的表达式,最后用数值算例验证了结论的正确性。     

具有圆形极点和方差约束的不确定离散系统鲁棒控制

对一类具有不确定参数和扰动的线性离散系统，研究了使闭环系统的稳态状态方差小于某个给定的上界，同时闭环极点位于一给定圆盘中的状态反馈鲁棒方差控制器设计问题。导出了控制器存在的条件，并证明了该条件等价于一个线性矩阵不等式系统的可行解问题，进而用这组矩阵系统的可行解给出了一组所求控制器的一个参数化表示。最后通过一个算例说明了该方法的可行性和有效性。     

一类矩阵方程的迭代解法

利用Kronecker积、矩阵拉直运算以及线性方程组的迭代解法,给出了解矩阵方程A1X1B1 A2X2B2 … AnXnBn=D的一种新的迭代方法,并给出了它的几种特殊形式.     

非对称线性方程组的块广义最小向后误差算法

许多科学领域都需要求多个右边值的大型非对称线性方程组，使用块方法同时计算所有的方程组比分别计算每一个方程要有效得多。因此，能同时计算所有方程的块迭代方法比单独计算每一个方程的迭代法要有效得多。本提出了一个块GMBACK方法求解有多个右边值的大型非对称线性方程组，该方法利用块Arnoldi过程构造Krylov子空间来求解Xm∈X0 Km(A,R0)使得矩阵A的扰动范数最小。     

关于一类非线性连续Leontief模型的条件投入产出方程的可解性结果

引入了一类非线性连续型Leontief模型及对应的条件投入产出方程．提出在正向或负向边界条件下的两个基本问题。通过应用非线性分析方法．得到了该方程的可解性结果与对应解的扰动性质，并通过注解说明了其经济意义。     

非确定结构系统区间分析的直接优化法

工程中的非确定性问题可以用区间分析、随机理论或模糊集理论进行求解。本文采用区间分析法来处理结构静力分析和设计中的不确定性问题。将结构系统中的不确定性参数用区间数来表示，用有限元法建立系统的控制方程。该控制方程是线性区间方程组。本文讨论了一些求解区间线性方程组的方法并提出了一种直接优化法。将方程组中的所有区间数都作为设计变量，区间量的变化区间作为相应的设计变量的边界约束，运用约束优化法求出方程组解的各元素的最大值和最小值。文中给出了两个算例，列出了本文算法与其他算法的结果比较     

奇异摄动法解最优拦截问题

本文应用奇异摄动法探讨水平面内飞机的中远程拦截问题。在介绍航迹优化的一般方法之后,通过采用三重时间尺度分离法建立奇异摄动的数学模型,求出了外层解、边界层解和一致有效的组合解,导出了以推力参数和滚转角为控制变量的最短时间拦截的最优控制律。并以F-4飞机为例,计算了拦截成功的最短时间和最优航迹。其结果可作为空中拦截实时在线优化的参考。     

输出跟踪问题中内部动态有界解的存在性

由于系统的内部动态不稳定，非最小相位系统的输出跟踪问题是一个具有挑战性的非常困难的问题，非最小相位系统在跟踪期望的输出轨这时，内部动态的不稳定变量可能发散。研究表明，不稳定内部动态有界解的存在性是非最小相位系统能够稳定跟踪期望输出轨迹的关键。对于系统内部动态具有双曲型形式的情形，Ｄｅｖａｓｉａ，Ｐａｄｅｎ和Ｃｈｅｎ利用对于不稳内部动态变量反时间方向积分的方法，证明了系统的内部地部存在有界解。文中讨     

白噪声参激余维2分岔系统奇点分析及最大Lyapunov指数

在随机动力系统中,最大Lyapunov指数是定义随机分岔系统概率1意义分岔的重要指标.为了分析白噪声参数激励对三维中心流形上一类余维2分岔系统的影响,本文基于一维扩散过程的奇点理论,通过L Arnold的摄动方法,针对不同的噪声作用项系数矩阵讨论了相应的奇异边界分类情况和最大Lyapunov指数渐近分析式.     

一类有阻尼的强非线性自由振动问题的一种解法

本文用插值摄动法及改进的Ｌ－Ｐ法相结合的方法求解了一类有限尼的强非线性自由振动问题，得到了一级近似解。     

逐次逼近摄动法及其在力学中的应用

提出逐次逼近摄动法，将预估的零级近似解逐渐加以修正，使之逐渐逼近于其真实解。有非线性问题的算例。算例表明，解法的精度和其它奇异摄动解法（如L—P法）的基本相同，但方法允许进行一次选代，使其精度更加提高。按照方法，可求得任一级近似解，其精度逐级提高，计算过程也很简便。     

数值模拟具有四个腔室的多分支管道二维非定常流动

应用基于Ｇｏｄｕｎｏｖ格式的二维非定常流数学模型，模拟具有四个腔室的多分支管道内的瞬变气体运动，以便对该气体运动的发展过程有一更好的了解。根据所给出的模型编制了相应的计算及绘图软件。数值模拟定量给出了多分支管道内详细的速度矢量场，压力样梯度场及各分支管道的质量流量变化曲线。并对多分支管道内非定常气体运动过程的某些特征做出了分析和讨论。计算结果表明：基于Ｇｏｄｕｎｏｖ格式的非线性波作用计算方法能成功地处理合压力波（包括激波）的传播发展和含复杂几何形状边界的多分支管道流动问题。因此该方法可以作为一个有效的数值分析工具，用来精确预测多分支管道内详细的压力变化、速度变化和质量流量变化的过程。     

不确定离散系统二次稳定的充分必要条件

研究了具有快时变参数变化离散系统的鲁棒稳定性问题，给出了该系统二次稳定（ｑｕａｄｒａｔｉｃｓｔａｂｉｌｉｔｙ）的充要条件，并且进一步定义了稳定鲁棒性指标来度量系统的稳定鲁棒性。对于参数具有已知变化范围的系统，其充要条件的检验和稳定鲁棒指标的计算转化为一个极小极大优化问题。当参数变化范围为一个凸多面体时，该极大化问题退化成其端点值间的相互比较，并证明了极小化问题中的任意局部最小值就是全局最小值     

Exact 3D Thermoelastic Solutions for a Penny-Shaped Crack in an Infinite Magnetoelectric Medium

Exact solutions of three-dimensional （3D） crack problems are much less in number than those of two-di- mensional ones, especially for multi-field coupling media exhibiting a certain kind of material anisotropy. An exact 3D thermoelastic solution has been reported for a uniformly heated penny-shaped crack in an infinite magnetoelec- tric space, with impermeable electromagnetic conditions assumed on the crack faces. Exact 3D solutions for the penny-shaped crack subjected to uniform or point temperature load are further presented here when the crack faces are electrically and magnetically permeable. The solutions, obtained by the potential theory method, are exact in the sense that all field variables are explicitly derived and expressed in terms of elementary functions. Along with the previously reported solution, the limits or bounds of the stress intensity factor at the crack-tip for a practical crack can be identified.