求分块周期三对角方程组的新算法

得到了求解系数矩阵为分块周期三对角矩阵线性方程组的一种新算法.     

分块带状矩阵的逆

1引言如果分块矩阵A=(A_(ij))_(n×n)满足A_(ij)=O(j-i＞p且i-j＞q),其中A_(ij)为m阶矩阵,则称A为(p,q)-分块带状矩阵．分块带状矩阵在一些实际问题中经常出现,例如在量子场论中用途很广的非线性Schr(?)dinger方程的差分离散问题,解热传导问题等,都会遇到分块带状矩阵．常见的分块三对角矩阵,分块五对角矩阵都是特殊的分块带状矩阵．采用通常的方法求解分块带状矩阵的逆矩阵时,需要进行O(n~3)次m阶矩阵的运算．本文首先将分块带状矩阵扩充成可逆的分块上(下)三角矩阵,利用其逆矩阵导出了分块带状矩阵的逆矩阵表达式；进而利用所得到的公式分别推导了分块三对角矩阵及分块五对角矩阵的逆矩阵的快速算法,所需运算量为O(n~2)次m阶矩阵的运算．本文的结果扩充了文[1]等关于分块三对角阵求逆的相关结果．     

基于角点处扭矢构造双三次Coons曲面的一种优化方法

本文针对矩形网格角点处的扭矢采用优化方法构造双三次Coons曲面,提出一种新的优化准则来确定角点处的扭矢.首先,通过变分原理,考虑曲面导矢的极小化问题转化的Euler-Lagrange偏微分方程,将该方程应用于每一个Coons块的角点上,引入一个新的极小化问题,其解是Euler-Lagrange偏微分方程的近似最优解.然后,建立一个具有块三对角系数矩阵的线性方程组来求解新的极小化问题.该系数矩阵可以表示为两个相同的形式特殊的矩阵的Kronnecker积,进而可以证明其非奇异性.最后,数值实验验证本文方法的稳定性和有效性.     

由两个右特征对构造三对角四元数矩阵

讨论由两个右特征对构造三对角四元数矩阵的数值求解问题,给出了该问题有解的充要条件,以及解的具体表达式.在已知两个特征对的条件下,进一步给出了三对角自共轭、三对角正定四元数矩阵的存在条件及计算方法.     

三对角矩阵逆的估计

1引言 三对角矩阵出现在很多应用中,例如,在求解常系数微分方程的比值问题,三次样条插值等应用中都会遇到三对角矩阵.因此这类矩阵非常重要,而且也有很多学者致力于这类矩阵的研究.在一些应用中,比如估计条件数和构造稀疏近似逆预条件子,需要计算三对角矩阵的逆,或者估计其逆元素的界.文献[1-7]给出了关于三对角矩阵逆的一些很好的结果,但是,这些结果大都建立在矩阵对角占优的条件之下,这限制了他们的应用.在本文中,我们给出一种一般三对角矩阵逆元素的估计办法.     

三对角逆M矩阵的判定

1、引言 三对角逆M矩阵是指同时为三对角矩阵和逆M矩阵的一类特殊矩阵．文用图论方法探讨三对角逆M矩阵结构,给出了三对角矩阵为逆M矩阵的充分必要条件．此条件提供了判定三对角矩阵是逆M矩阵的方法,但较复杂．文讨论了这类矩阵在Hadamard积下的封闭性．由于三对角逆M矩阵在理论和应用上都有一定价值,所以,寻求一种简单而实用的判定方法是必要的．本文通过对这类矩阵结构特点的研究找到了这样一种方法．同时,由此证明了这类矩阵在Hadamard积下的封闭性．     

三对角对称矩阵逆特征值问题存在唯一解的条件

讨论了利用给定的k(2≤k≤n)个特征对来构造相应的三对角对称矩阵的问题．在求解方法中,将已知的一些关系式等价转化成线性方程组,利用线性方程组的解存在唯一的条件,得到了所研究问题存在唯一解的充要条件,并给出了计算解的数值方法和数值实例．     

一类对称三对角矩阵的合同对角化算法的实现

从一个对称三对角矩阵的合同变换出发 ,阐述了对称三对角矩阵对应的二次型标准化的一种方法 .     

三对角矩阵计算

1 引言 在数值计算中,有许多问题最后归结为三对角矩阵的计算,因此研究它们的计算方法是有意义的。此外,有些三对角阵的计算方法可以做为带状阵计算的借鉴。 本文讨论三对角线性方程组的解耦算法,矩阵的LR~(-1)分解,求行列式,Jacobi矩阵的特征值与特征向量的关系以及三对角阵求逆等方面的问题,与现有的算法比较,本文的算法具有计算量或存贮量较少,或计算精度较高,或编程较简单等某些特点。 设A为n阶非奇实三对角阵:     

求解块三对角方程组的新算法

根据块三对角矩阵的特殊分解,给出了求解块三对角方程组的新算法.该算法含有可以选择的参数矩阵,适当选择这些参数矩阵,可以使得计算精度较著名的追赶法高,甚至当追赶法失效时,由该算法仍可得到一定精度的解.     

求解三维高次拉格朗日有限元方程的代数多重网格法

本文针对带有间断系数的三维椭圆问题，讨论任意四面体剖分下的二次拉格朗日有限元方程的代数多重网格法．通过分析线性和高次有限元空间之间的关系，我们给出了一种新的网格粗化算法和构造提升算子的代数途径．进一步，我们还对新的代数多重网格法给出了收敛性分析．数值实验表明这种代数多重网格法对求解二次拉格朗日有限元方程是健壮和有效的。     

解浅水方程离散化线性方程组的PEk方法

定义于球面的浅水方程能够很好地描述浅齐次的不可压缩非黏滞流体层的性状,它在全球大气模型、海洋数字模型和天气预报的数值计算中都有广泛的应用,浅水方程的一般形式如下：     

On spectral decompositions of solutions to discrete Lyapunov equations

A new approach to solving discrete Lyapunov matrix algebraic equations is based on methods for spectral decomposition of their solutions. Assuming that all eigenvalues of the matrices on the left-hand side of the equation lie inside the unit disk, it is shown that the matrix of the solution to the equation can be calculated as a finite sum of matrix bilinear quadratic forms made up by products of Faddeev matrices obtained by decomposing the resolvents of the matrices of the Lyapunov equation. For a linear autonomous stochastic discrete dynamic system, analytical expressions are obtained for the decomposition of the asymptotic variance matrix of system’s states.     

Chebyshev polynomial approximation for high‐order partial differential equations with complicated conditions

In this article, a new method is presented for the solution of high‐order linear partial differential equations (PDEs) with variable coefficients under the most general conditions. The method is based on the approximation by the truncated double Chebyshev series. PDE and conditions are transformed into the matrix equations, which corresponds to a system of linear algebraic equations with the unknown Chebyshev coefficients, via Chebyshev collocation points. Combining these matrix equations and then solving the system yields the Chebyshev coefficients of the solution function. Some numerical results are included to demonstrate the validity and applicability of the method. © 2008 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2009     

A new collocation method for solution of mixed linear integro-differential-difference equations

Numerical solution of mixed linear integro-differential-difference equation is presented using Chebyshev collocation method. The aim of this article is to present an efficient numerical procedure for solving mixed linear integro-differential-difference equations. Our method depends mainly on a Chebyshev expansion approach. This method transforms mixed linear integro-differential-difference equations and the given conditions into matrix equation which corresponds to a system of linear algebraic equation. The reliability and efficiency of the proposed scheme are demonstrated by some numerical experiments and performed on the computer algebraic system Maple10.     

On the numerical solution of stochastic quadratic integral equations via operational matrix method

In this paper, stochastic operational matrix of integration based on delta functions is applied to obtain the numerical solution of linear and nonlinear stochastic quadratic integral equations (SQIEs) that appear in modelling of many real problems. An important advantage of this method is that it dose not need any integration to compute the constant coefficients. Also, this method can be utilized to solve both linear and nonlinear problems. By using stochastic operational matrix of integration together collocation points, solving linear and nonlinear SQIEs converts to solve a nonlinear system of algebraic equations, which can be solved by using Newton's numerical method. Moreover, the error analysis is established by using some theorems. Also, it is proved that the rate of convergence of the suggested method is O(h²). Finally, this method is applied to solve some illustrative examples including linear and nonlinear SQIEs. Numerical experiments confirm the good accuracy and efficiency of the proposed method.     

Sequential synthesis of time-optimal control for linear systems with disturbances

A method of sequential synthesis of time-optimal control for a linear system with unknown disturbances is considered. A system of linear algebraic equations is obtained which relates the increments of phase coordinates to the increments of initial conditions of a normalized adjoint system and to the increment of control completion time. Evaluations consist in solving repeatedly a system of linear algebraic equations and integrating a matrix differential equation on the displacement intervals of control switching times and on the displacement interval of final control time. A procedure of correcting the switching times and the completion time in moving along the phase trajectory of a controllable object is examined. Simple and constructive conditions are specified for a discontinuous mode to occur, for a representation point to move along the switching manifolds, and for the optimal control structure to transform in moving along the phase trajectory of a system with uncontrollable disturbance. A computational algorithm is presented. It is proved that a sequence of controls converges locally at a quadratic rate and globally to a time-optimal control.     

一类特征值反问题(IEP)的基于矩阵方程的Ulm型算法

应用求解算子方程的Ulm方法构造了求解一类矩阵特征值反问题（IEP）的新算法.所给算法避免了文献[Aishima K.，A quadratically convergent algorithm based on matrix equations for inverse eigenvalue problems，Linear Algebra and its Applications，2018，542：310-33]中算法在每次迭代中要求解一个线性方程组的不足，证明了在给定谱数据互不相同的条件下所给算法具有根收敛意义下的二次收敛性.数值实验表明本文所给算法在矩阵阶数较大时计算效果优于上文所给算法.     

Stochastic version of the degenerate kernel method for solving large systems of linear equations with sparse matrix

A stochastic method for solving systems of linear algebraic equations with sparse matrix is proposed. The method is based on the determination of the principal part of a matrix operator in the form of a degenerate kernel. Results of a numerical experiment witnessing the efficiency of the method are presented.     

On the convergence of inexact Newton methods for discrete-time algebraic Riccati equations

In this paper, we present a convergence analysis of the inexact Newton method for solving Discrete-time algebraic Riccati equations (DAREs) for large and sparse systems. The inexact Newton method requires, at each iteration, the solution of a symmetric Stein matrix equation. These linear matrix equations are solved approximatively by the alternating directions implicit (ADI) or Smith?s methods. We give some new matrix identities that will allow us to derive new theoretical convergence results for the obtained inexact Newton sequences. We show that under some necessary conditions the approximate solutions satisfy some desired properties such as the d-stability. The theoretical results developed in this paper are an extension to the discrete case of the analysis performed by Feitzinger et al. (2009) [8] for the continuous-time algebraic Riccati equations. In the last section, we give some numerical experiments.