共查询到20条相似文献,搜索用时 0 毫秒
1.
2.
Peter Benner Grece El Khoury Miloud Sadkane 《Numerical Linear Algebra with Applications》2014,21(5):645-665
A squared Smith type algorithm for solving large‐scale discrete‐time Stein equations is developed. The algorithm uses restarted Krylov spaces to compute approximations of the squared Smith iterations in low‐rank factored form. Fast convergence results when very few iterations of the alternating direction implicit method are applied to the Stein equation beforehand. The convergence of the algorithm is discussed and its performance is demonstrated by several test examples. Copyright © 2013 John Wiley & Sons, Ltd. 相似文献
3.
K. Jbilou 《Linear algebra and its applications》2010,432(10):2473-1398
In the present paper, we propose preconditioned Krylov methods for solving large Lyapunov matrix equations AX+XAT+BBT=0. Such problems appear in control theory, model reduction, circuit simulation and others. Using the Alternating Direction Implicit (ADI) iteration method, we transform the original Lyapunov equation to an equivalent symmetric Stein equation depending on some ADI parameters. We then define the Smith and the low rank ADI preconditioners. To solve the obtained Stein matrix equation, we apply the global Arnoldi method and get low rank approximate solutions. We give some theoretical results and report numerical tests to show the effectiveness of the proposed approaches. 相似文献
4.
5.
Tiexiang Li Peter Chang-Yi Weng Eric King-wah Chu Wen-Wei Lin 《Numerical Algorithms》2013,63(4):727-752
We consider the solution of large-scale Lyapunov and Stein equations. For Stein equations, the well-known Smith method will be adapted, with $A_k = A^{2^k}$ not explicitly computed but in the recursive form $A_k = A_{k-1}^{2}$ , and the fast growing but diminishing components in the approximate solutions truncated. Lyapunov equations will be first treated with the Cayley transform before the Smith method is applied. For algebraic equations with numerically low-ranked solutions of dimension n, the resulting algorithms are of an efficient O(n) computational complexity and memory requirement per iteration and converge essentially quadratically. An application in the estimation of a lower bound of the condition number for continuous-time algebraic Riccati equations is presented, as well as some numerical results. 相似文献
6.
BIT Numerical Mathematics - Linear and nonlinear evolution equations have been formulated to address problems in various fields of science and technology. Recently, methods using an exponential... 相似文献
7.
In this paper, we introduce an extension of multiple set split variational inequality problem (Censor et al. Numer. Algor. 59, 301–323 2012) to multiple set split equilibrium problem (MSSEP) and propose two new parallel extragradient algorithms for solving MSSEP when the equilibrium bifunctions are Lipschitz-type continuous and pseudo-monotone with respect to their solution sets. By using extragradient method combining with cutting techniques, we obtain algorithms for these problems without using any product space. Under certain conditions on parameters, the iteration sequences generated by the proposed algorithms are proved to be weakly and strongly convergent to a solution of MSSEP. An application to multiple set split variational inequality problems and a numerical example and preliminary computational results are also provided. 相似文献
8.
The extended Krylov subspace method has recently arisen as a competitive method for solving large-scale Lyapunov equations.
Using the theoretical framework of orthogonal rational functions, in this paper we provide a general a priori error estimate
when the known term has rank-one. Special cases, such as symmetric coefficient matrix, are also treated. Numerical experiments
confirm the proved theoretical assertions. 相似文献
9.
In the present paper, we propose a hierarchical identification method (SSHI) for solving Lyapunov matrix equations, which is based on the symmetry and skew-symmetry splitting of the coefficient matrix. We prove that the iterative algorithm consistently converges to the true solution for any initial values with some conditions, and illustrate that the rate of convergence of the iterative solution can be enhanced by choosing the convergence factors appropriately. Furthermore, we show that the method adopted can be easily extended to study iterative solutions of other matrix equations, such as Sylvester matrix equations. Finally, we test the algorithms and show their effectiveness using numerical examples. 相似文献
10.
We discuss the numerical solution of large-scale discrete-time algebraic Riccati equations (DAREs) as they arise, e.g., in fully discretized linear-quadratic optimal control problems for parabolic partial differential equations (PDEs). We employ variants of Newton??s method that allow to compute an approximate low-rank factor of the solution of the DARE. The principal computation in the Newton iteration is the numerical solution of a Stein (aka discrete Lyapunov) equation in each step. For this purpose, we present a low-rank Smith method as well as a low-rank alternating-direction-implicit (ADI) iteration to compute low-rank approximations to solutions of Stein equations arising in this context. Numerical results are given to verify the efficiency and accuracy of the proposed algorithms. 相似文献
11.
In this note we present a modified cyclic low-rank Smith method to compute low-rank approximations to solutions of Lyapunov equations arising from large-scale dynamical systems. Unlike the original cyclic low-rank Smith method introduced by Penzl in [20], the number of columns required by the modified method in the approximate solution does not necessarily increase at each step and is usually much lower than in the original cyclic low-rank Smith method. The modified method never requires more columns than the original one. Upper bounds are established for the errors of the low-rank approximate solutions and also for the errors in the resulting approximate Hankel singular values. Numerical results are given to verify the efficiency and accuracy of the new algorithm. 相似文献
12.
Zeng Zhijun 《Journal of Mathematical Sciences》2009,161(2):337-343
This paper presents converse Lyapunov theorems for exponential stability of nonautonomous discrete-time systems with disturbances
and free of disturbances, respectively. It is shown that Lyapunov functions exist for discrete-time systems if the systems
are exponentially stable. Moreover, in the periodic case, we explicitly construct a Lyapunov function for systems with disturbances.
Translated from Sovremennaya Matematika i Ee Prilozheniya (Contemporary Mathematics and Its Applications), Vol. 61, Optimal
Control, 2008. 相似文献
13.
A new version of the Smith method for solving Sylvester equation and discrete-time Sylvester equation
下载免费PDF全文
![点击此处可从《Journal of Applied Analysis & Computation》网站下载免费的PDF全文](/ch/ext_images/free.gif)
Recently, Xue etc. \cite{28}
discussed the Smith method for solving Sylvester equation $AX+XB=C$,
where one of the matrices $A$ and $B$ is at least a nonsingular
$M$-matrix and the other is an (singular or nonsingular) $M$-matrix.
Furthermore, in order to find the minimal non-negative solution of a
certain class of non-symmetric algebraic Riccati equations, Gao and
Bai \cite{gao-2010} considered a doubling iteration scheme to
inexactly solve the Sylvester equations. This paper discusses the
iterative error of the standard Smith method used in \cite{gao-2010}
and presents the prior estimations of the accurate solution $X$ for
the Sylvester equation. Furthermore, we give a new version of the
Smith method for solving discrete-time Sylvester equation or Stein
equation $AXB+X=C$, while the new version of the Smith method can
also be used to solve Sylvester equation $AX+XB=C$,
where both $A$ and $B$ are positive definite. % matrices.
We also study the convergence rate of the new Smith method. At last, numerical examples are given to illustrate
the effectiveness of our methods 相似文献
14.
In this paper, we propose two derivative-free iterative methods for solving nonlinear monotone equations, which combines two modified HS methods with the projection method in Solodov and Svaiter (1998) [5]. The proposed methods can be applied to solve nonsmooth equations. They are suitable to large-scale equations due to their lower storage requirement. Under mild conditions, we show that the proposed methods are globally convergent. The reported numerical results show that the methods are efficient. 相似文献
15.
In the present paper, we propose a preconditioned Newton–Block Arnoldi method for solving large continuous time algebraic Riccati equations. Such equations appear in control theory, model reduction, circuit simulation amongst other problems. At each step of the Newton process, we solve a large Lyapunov matrix equation with a low rank right hand side. These equations are solved by using the block Arnoldi process associated with a preconditioner based on the alternating direction implicit iteration method. We give some theoretical results and report numerical tests to show the effectiveness of the proposed approach. 相似文献
16.
In this paper, by the use of the residual vector and an approximation to the steepest descent direction of the norm function, we develop a norm descent spectral method for solving symmetric nonlinear equations. The method based on the nomonotone line search techniques is showed to be globally convergent. A specific implementation of the method is given which exploits the recently developed cyclic Barzilai–Borwein (CBB) for unconstrained optimization. Preliminary numerical results indicate that the method is promising. 相似文献
17.
用迭代Lyapunov函数方法对微分代数系统稳定域进行了研究,根据所研究的微分代数系统形式,构造一个Lyapunov函数,然后对这个Lyapunov函数进行逐次迭代,给出了微分代数系统稳定域逐次扩大的迭代算法,数值实验表明迭代Lyapunov函数方法应用于微分代数系统稳定域的估计比单个Lyapunov函数具有良好的优越性。 相似文献
18.
《Journal of Computational and Applied Mathematics》2005,176(2):463-466
We characterize the stability of discrete-time Lyapunov equations with periodic coefficients. The characterization can be seen as the analog of the classical stability theorem of Lyapunov equations with constant coefficients. It involves quantities readily computable with good accuracy. 相似文献
19.
Numerical Algorithms - A fully derivative-free conjugate residual method, using secant condition, is introduced to solve general large-scale nonlinear equations. Under some conditions, global and... 相似文献
20.
This paper is devoted to the study of Lp Lyapunov-type inequalities (1?p?+∞) for linear partial differential equations. More precisely, we treat the case of Neumann boundary conditions on bounded and regular domains in RN. It is proved that the relation between the quantities p and N/2 plays a crucial role. This fact shows a deep difference with respect to the ordinary case. The linear study is combined with Schauder fixed point theorem to provide new conditions about the existence and uniqueness of solutions for resonant nonlinear problems. 相似文献