大型线性方程组(集)的分段分片直接解法及其与子结构法的偶合

本文提出了一种求解大型线代数方程(正定)组集的直接解法,称为分段、分片解法(D.S.P)。它是基于LDL~T分解的一种改进算法,可用于含多组右端向量的方程集,并适用于大带宽矩阵,可兼顾节省计算机容量与时间。与文[1]、[2]相比,本解法用分片处理克服了带宽过大时引起的所谓“超段”问题;并提出与子结构法偶合,以有效地处理多叉式结构。     

非线性半定规划若干进展

讨论非线性半定规划的四个专题, 包括半正定矩阵锥的变分分析、非凸半定规划问题的最优性条件、非凸半定规划问题的扰动分析和非凸半定规划问题的增广Lagrange方法.     

非线性半定规划若干进展

讨论非线性半定规划的四个专题,包括半正定矩阵锥的变分分析、非凸半定规划问题的最优性条件、非凸半定规划问题的扰动分析和非凸半定规划问题的增广Lagrange方法．     

半定规划的近似中心投影法

1．引言半定规划问题标准形的数学形式是这里C，AIEIR”””及变量XEIRn“”为对称矩阵，Tr（·）表示矩阵的迹，用符号＞0和三0分别表示矩阵正定和半正定．由于半定规划在控制论，结构优化，组合优化方面有重要应用［1，3，16，17］以及线性规划内点法取得的巨大成就［7］，将线性规划的内点法推广到半定规划上，是数学规划领域内近年来受到重视的一个研究课题．线性规划内点法中的势函数下降法［10，16］原始对偶中心路径跟踪法［2，4，8，9，11。15］已经先后被推广到半定规划上．ROOS－Visl近似中心法则是求解线性规划的另一类内…     

拟线性微分代数方程的标准形

讨论了拟线性微分代数方程在一类特殊的奇点-拟障碍点附近的标准形.通过矩阵广义逆理论,拟线性微分代数方程可化为半显式形式.然后运用标准形理论,在微分同胚变换下,给出了拟线性微分代数方程在拟障碍点附近的标准形.在此基础上进一步讨论了这类标准形的去奇异化性质.     

有限元线代数方程组的杂交算法及相应的节点编号方法

在有限元问题的计算中,大家知道耗时最多的是求解线代数方程组的部分。其计算方法一般可分为直接法和迭代法两大类。这两类方法从有限元发展历史看,基本上是以直接法占优势,特别是大型计算机的发展更加强了这个优势。但是从七十年代中后期起,有限元的发展由线性问题向非线性问题过渡,由二维问题向三维问题过渡,特别是微机     

求解最大割问题的半定规划松驰的序列线性规划方法

本文基于最大割问题的半定规划松弛,利用矩阵分解的方法给出了与半定规划松弛等价的非线性规划模型,提出一种序列线性规划方法求解该模型.并在适当的条件下,证明了算法的全局收敛性.数值实验表明:序列线性规划方法在时间上要优于半定规划的内点算法.所以序列线性规划方法能更有效地求解大规模的最大割问题的半定规划松弛.     

分片Bernstein多项式的样条配点法求解四阶微分方程

本文对四阶微分方程边值问题,给出一种基于分片Bernstein多项式的样条配点法求解,该格式构造过程容易理解,形成的线代数方程组系数矩阵稀疏,可用迭代法求解.数值实验表明,该方法可有效求解一般四阶线性微分方程边值问题,结合非等距配置点亦可用于求解含小参数的扰动问题.     

矩阵位移法在连系梁计算中的应用

矩阵位移法是结构矩阵分析法的一种,矩阵位移法的基本体系和节点位移未知量的选择一般是唯一的,它容易编制通用的计算机程序,矩阵位移法又分为刚度法和直接刚度法,二者的基本原理并无本质区别,只在形成总刚度矩阵时使用方法不同,直接刚度法相对简便,应用也较广泛,本文研究直接刚度法在连续梁内力计算的应用.     

计算应力强度因子的无限相似单元法

为了计算应力强度因子,将平面弹性体所占的区域作三角形单元剖分。在裂纹尖端附近,三角形单元为可数无穷多个,且具有相似性,它们组成组合单元体。用矩阵方法求出组合单元体的组合刚度矩阵,从而得到问题的解。与有限单元法比较,无限相似单元法更好地反映了解的奇性,因此提高了精确度。同时,由于组合刚度矩阵的求得,使代数方程组的阶数大为降低。     

半定优化与半光滑牛顿算法

本文讨论半光滑牛顿算法的基本概念与其在求解半定优化问题中的应用．特别地，该算法可用于求解线性或非线性半定互补问题．本文同时综述最近在矩阵方程，增广拉格朗日公式和半定优化稳定性方面的、源于半光滑牛顿算法的理论成果．     

循环矩阵及其在结构计算中的应用(Ⅱ)

在[1]中提出了利用循环矩阵进行结构计算的方法.如果区域和剖分都是规则的,在周期边界条件的情况下,形成的代数方程系数矩阵为循环矩阵,在第一类或第二类边界条件的情况下,系数矩阵为准循环矩阵.[1]中对这种类型的方程进行了讨论,利用快速富氏变换的工具,得到一个速度快而且存贮量节省的计算方法,所述的方法对某些非规则区     

解强刚性块线代数方程组的两类L-收敛迭代法

对解强刚性块线代数方程组X=(A(?)J)X φ,本文提出了L-收敛的最佳单参数迭代法(L-OOPI)和L-收敛的多参数迭代直接法(L-MPID),并给出了数值例子．数例表明,对于强刚性块线代数方程组,该二迭代法是有效的．     

求方程近似解的三种方法

我们知道,在许多实际问题中,经常会遇到超越方程、高次代数方程或其他类型的方程．要想求得这类方程的精确值,几乎是不可能的,因此需要寻求方程的近似解．对于求方程近似解的方法,吴老师给同学们上了一次这样的习题课：首先,她出了一道题：如图1所示,设有一半径为R的半圆,试过原点O做一直线L,将半圆面积二等分,求L的倾角a,要求精确到小数点第二位．同学们经过思考后,认为应根据问题的要求建立关于a的方程,他们采用了两种比较简单的方法．－（1）利用初等数学知识,如图12）利用定积分知识：,均可得到方程（其中x=2a）接着,…     

特殊矩阵乘积在非线性数值计算中的应用

本文将Hadamard矩阵乘积引入到非线性数值计算,获得了简单的矩阵形式的非线性代数模拟方程,利用Hadamard矩阵乘积和Hadamard矩阵函数的方法,我们能够容易地构造快速收敛的简单迭代法解非线性代数方程组的迭代公式,使该法成为与Newton-Raphson法相比有竞争力的方法,我们也首次定义了一种新的特殊矩阵乘积—SJT矩阵乘积。运用SJT积,我们能够方便高效的计算Newton-Raphson法中Jacobi导数矩阵的精确解,利用Hadamard矩阵乘积的范数性质,我们也导出了非线性计算摄动误差的分析公式,此外,Hadamard积和SJT积能够被用于非线性数值解耦计算,这极大地减少了求解耦合的非线性偏微分方程组的计算工作量和内存需要量。     

多参数结构特征二阶灵敏度

提出了一种有效计算多参数结构特征值与特征向量二阶灵敏度矩阵--Hessian矩阵的方法.将特征值和特征向量二阶摄动法转变为多参数形式,推导出二阶摄动灵敏度矩阵,由此得到特征值和特征向量的二阶估计式.该法解决了无法用直接求导法计算特征值和特征向量二阶灵敏度矩阵的问题.数值算例说明了该算法的应用和计算精度.     

多目标半定规划的互补弱鞍点和G-鞍点最优性条件

对于含矩阵函数半定约束和多个目标函数的多目标半定规划问题，给出Lagrange函数在弱有效意义下的互补弱鞍点和Geofrrion恰当有效意义下的G-鞍点的定义及其等价定义．然后，在较弱的凸性条件下，利用含矩阵和向量约束的择一性定理，建立多目标半定规划的互补弱鞍点和G-鞍点充分必要条件．     

利用第二类Chebyshev小波求二阶常微分方程组边值问题的数值解

给出了一个求二阶常微分方程组边值问题数值解的第二类Chebyshev小波配点法.利用第二类Chebyshev小波积分算子矩阵,将问题转化成代数方程组的运算.数值例子说明了方法的准确性及易操作性.另外,为了表明方法的高精度性和有效性,数值算例结果与解析解,以及运用变分迭代法,B样条配点法,连续遗传算法等得到的结果进行了比较.     

非凸半定规划的广义Fakars引理及最优性条件

1引言在本文中,我们用(?),S~n,S_ ~n分别表示有限维向量空间,n阶对称矩阵空间及n阶半正定矩阵锥.我们考虑如下形式的非凸半定规划问题:     

连续型线性定常系统的Riccati代数方程的小摄动问题

本文讨论了连续型线性定常系统摄动的Riccati代数方程所对应的稳定性问题.通过矩阵范数分析建立了摄动的Riccati代数方程的解的摄动界估计(以系统参数摄动界表出),从而提供了一种方便的实用计算方法.