ALBERT

All Library Books, journals and Electronic Records Telegrafenberg

Your email was sent successfully. Check your inbox.

An error occurred while sending the email. Please try again.

Proceed reservation?

Export
Filter
  • Articles  (73)
  • Articles: DFG German National Licenses  (73)
  • CR: G1.3  (42)
  • Optimal control  (31)
  • 2015-2019
  • 1985-1989  (73)
  • 1945-1949
  • Mathematics  (73)
  • Architecture, Civil Engineering, Surveying
Collection
  • Articles  (73)
Source
  • Articles: DFG German National Licenses  (73)
Publisher
Years
Year
Topic
  • 1
    Electronic Resource
    Electronic Resource
    Springer
    Numerische Mathematik 50 (1986), S. 263-271 
    ISSN: 0945-3245
    Keywords: AMS(MOS): 65F10 ; CR: G1.3
    Source: Springer Online Journal Archives 1860-2000
    Topics: Mathematics
    Notes: Summary We study semi-iterative (two step interative) methods of the form $$x_{n + 1} = T^* T(\alpha x_n + \gamma x_{n - 1} ) + \beta x_n + (1 - \beta )x_{n - 1} - (\alpha + \gamma )T^* y$$ for the approximate solution of ill-posed or ill-conditioned linear equationsTx=y in (infinite or finite dimensional) Hilbert spaces. We present results on convergence, convergence rates, the influence of perturbed data, and on the comparison of different methods.
    Type of Medium: Electronic Resource
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 2
    Electronic Resource
    Electronic Resource
    Springer
    Numerische Mathematik 52 (1987), S. 241-250 
    ISSN: 0945-3245
    Keywords: AMS (MOS): 15A12 ; 49D15 ; 65F35 ; CR: G1.3 ; G1.6
    Source: Springer Online Journal Archives 1860-2000
    Topics: Mathematics
    Notes: Summary We derive lower bounds for the ∞-condition number of then×n-Vandermonde matrixV n(x) in the cases where the node vectorx T=[x1, x2,...,xn] has positive elements or real elements located symmetrically with respect to the origin. The bounds obtained grow exponentially inn. withO(2n) andO(2n/2), respectively. We also compute the optimal spectral condition numbers ofV n(x) for the two node configurations (including the optimal nodes) and compare them with the bounds obtained.
    Type of Medium: Electronic Resource
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 3
    Electronic Resource
    Electronic Resource
    Springer
    Numerische Mathematik 54 (1989), S. 591-599 
    ISSN: 0945-3245
    Keywords: AMS(MOS): 65F05, 65G05, 15A06 ; CR: G1.3
    Source: Springer Online Journal Archives 1860-2000
    Topics: Mathematics
    Notes: Summary In this paper a Gauss-Jordan algorithm with column interchanges is presented and analysed. We show that, in contrast with Gaussian elimination, the Gauss-Jordan algorithm has essentially differing properties when using column interchanges instead of row interchanges for improving the numerical stability. For solutions obtained by Gauss-Jordan with column interchanges, a more satisfactory bound for the residual norm can be given. The analysis gives theoretical evidence that the algorithm yields numerical solutions as good as those obtained by Gaussian elimination and that, in most practical situations, the residuals are equally small. This is confirmed by numerical experiments. Moreover, timing experiments on a Cyber 205 vector computer show that the algorithm presented has good vectorisation properties.
    Type of Medium: Electronic Resource
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 4
    Electronic Resource
    Electronic Resource
    Springer
    Numerische Mathematik 55 (1989), S. 265-280 
    ISSN: 0945-3245
    Keywords: AMS(MOS) ; 65H10 ; 58C99 ; 55M25 ; CR: G1.3
    Source: Springer Online Journal Archives 1860-2000
    Topics: Mathematics
    Notes: Summary The homotopy method can be used to solve eigenvalue-eigenvector problems. The purpose of this paper is to report the numerical experience of the homotopy method of computing eigenpairs for real symmetric tridiagonal matrices together with a couple of new theoretical results. In practice, it is rerely of any interest to compute all the eigenvalues. The homotopy method, having the order preserving property, can provide any specific eigenvalue without calculating any other eigenvalues. Besides this advantage, we note that the homotopy algorithm is to a large degree a parallel algorithm. Numerical experimentation shows that the homotopy method can be very efficient especially for graded matrices.
    Type of Medium: Electronic Resource
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 5
    Electronic Resource
    Electronic Resource
    Springer
    Numerische Mathematik 52 (1987), S. 279-300 
    ISSN: 0945-3245
    Keywords: AMS(MOS): 65F15 ; CR: G1.3
    Source: Springer Online Journal Archives 1860-2000
    Topics: Mathematics
    Notes: Summary In this paper we compare several implementations of Kogbetliantz's algorithm for computing the SVD on sequential as well as on parallel machines. Comparisons are based on timings and on operation counts. The numerical accuracy of the different methods is also analyzed.
    Type of Medium: Electronic Resource
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 6
    Electronic Resource
    Electronic Resource
    Springer
    Numerische Mathematik 53 (1988), S. 571-593 
    ISSN: 0945-3245
    Keywords: AMS(MOS): 65F10 ; CR: G1.3
    Source: Springer Online Journal Archives 1860-2000
    Topics: Mathematics
    Notes: Summary The Chebyshev and second-order Richardson methods are classical iterative schemes for solving linear systems. We consider the convergence analysis of these methods when each step of the iteration is carried out inexactly. This has many applications, since a preconditioned iteration requires, at each step, the solution of a linear system which may be solved inexactly using an “inner” iteration. We derive an error bound which applies to the general nonsymmetric inexact Chebyshev iteration. We show how this simplifies slightly in the case of a symmetric or skew-symmetric iteration, and we consider both the cases of underestimating and overestimating the spectrum. We show that in the symmetric case, it is actually advantageous to underestimate the spectrum when the spectral radius and the degree of inexactness are both large. This is not true in the case of the skew-symmetric iteration. We show how similar results apply to the Richardson iteration. Finally, we describe numerical experiments which illustrate the results and suggest that the Chebyshev and Richardson methods, with reasonable parameter choices, may be more effective than the conjugate gradient method in the presence of inexactness.
    Type of Medium: Electronic Resource
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 7
    Electronic Resource
    Electronic Resource
    Springer
    Numerische Mathematik 54 (1988), S. 19-32 
    ISSN: 0945-3245
    Keywords: AMS(MOS): 65F20 ; CR: G1.3
    Source: Springer Online Journal Archives 1860-2000
    Topics: Mathematics
    Notes: Summary A direct method is developed for solving linear least squares problems $$\mathop {\min \left\| {Ax - b} \right\|_2 }\limits_x $$ , whereA is large and sparse and the solution is subject to lower and upper boundsl≦x≦u. The problem is initially transformed to upper triangular form by a sparseQR-factorization. An active set algorithm is then used. The key step is the stable updating of theR-factor associated with the columns ofA corresponding to the free variables, when theQ-factor is not available. For this a new method is developed, which uses the semi-normal equations and iterative refinement.
    Type of Medium: Electronic Resource
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 8
    Electronic Resource
    Electronic Resource
    Springer
    Numerische Mathematik 56 (1989), S. 627-633 
    ISSN: 0945-3245
    Keywords: AMS(MOS): 65F15 ; CR: G1.3
    Source: Springer Online Journal Archives 1860-2000
    Topics: Mathematics
    Notes: Summary We propose a “one-sided” or “implicit” variant of the Jacobi diagonalization algorithm for positive definite matrices. The variant is based on a previous Cholesky decomposition and currently uses essentially one square array which, on output, contains the matrix of eigenvectors thus reaching the storage economy of the symmetric QL algorithm. The current array is accessed only columnwise which makes the algorithm attractive for various parallelized and/or vectorized implementations. Even on a serial computer our algorithm shows improved efficiency, in particular if the Cholesky step is made with diagonal pivoting. On matrices of ordern=25–200 our algorithm is about 2–3.5 times slower than QL thus being almost on the halfway between the standard Jacobi and QL algorithms. The previous Cholesky decomposition can be performed with higher precision without extra time or storage costs thus offering considerable gains in accuracy with highly conditioned input matrices.
    Type of Medium: Electronic Resource
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 9
    Electronic Resource
    Electronic Resource
    Springer
    Numerische Mathematik 56 (1989), S. 721-734 
    ISSN: 0945-3245
    Keywords: AMS(MOS): 65F10 ; CR: G1.3
    Source: Springer Online Journal Archives 1860-2000
    Topics: Mathematics
    Notes: Summary The acceleration by Tchebychev iteration for solving nonsymmetric eigenvalue problems is dicussed. A simple algorithm is derived to obtain the optimal ellipse which passes through two eigenvalues in a complex plane relative to a reference complex eigenvalue. New criteria are established to identify the optimal ellipse of the eigenspectrum. The algorithm is fast, reliable and does not require a search for all possible ellipses which enclose the spectrum. The procedure is applicable to nonsymmetric linear systems as well.
    Type of Medium: Electronic Resource
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 10
    Electronic Resource
    Electronic Resource
    Springer
    Numerische Mathematik 54 (1989), S. 639-654 
    ISSN: 0945-3245
    Keywords: AMS(MOS): 35P15, 47A55, 65M15 ; CR: G1.3
    Source: Springer Online Journal Archives 1860-2000
    Topics: Mathematics
    Notes: Summary LetA, B be essentially self-adjoint and positive definite differential operators defined inL 2(G). Using Svirskij's construction of the base operator and some results from the analytic perturbation theory of linear operators a formula providing eigenvalue lower bounds of the problemAu=λBu is derived. In this formula a rough lower bound of some higher eigenvalue and the residual convergence of the Rayleigh-Ritz eigenfunction approximations are needed. Some numerical results are presented.
    Type of Medium: Electronic Resource
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 11
    Electronic Resource
    Electronic Resource
    Springer
    Numerische Mathematik 55 (1989), S. 463-476 
    ISSN: 0945-3245
    Keywords: AMS(MOS): 65F05, 65N30 ; CR: G1.3
    Source: Springer Online Journal Archives 1860-2000
    Topics: Mathematics
    Description / Table of Contents: Résumé On considère la méthode de dissections emboîtées basée sur des théorèmes de séparation introduite par Gilbert-Tarjan et Roman utilisée pour la résolution par élimination de Gauss de grands systèmes linéaires creux. Plus précisemment, on étudie une structure de données par blocs similaire à celle proposée par George dans le cadre des graphes en grille, et on démontre les propriétés suivantes: d'une part, pour des familles de graphes à degré borné admettant unn σ-théorème de séparation, 1/2≦σ〈1, le stockage secondaire de la structure par blocs contenant la matrice factorisée est linéaire par rapport à la taille du système; d'autre part, en rajoutant une hypothèse non restrictive sur la manière d'effectuer la séparation, la structure peut être construite en temps linéaire par une factorisation logique par blocs. Des exemples numériques illustrent ces résultats théoriques.
    Notes: Summary We consider the nested dissection method based on separator theorems introduced by Gilbert-Tarjan and Roman used for solving large sparse systems of linear equations. More precisely, we study a block storage scheme such as proposed by George for regular square grids and we prove the following results: first, for families of graphs of bounded degree withn σ-separator theorem, 1/2≦σ〈1, the overhead storage of the block data structure for the factored matrix is linear in system dimension; on the other hand, by adding a non restrictive assumption on the separation, this structure can be constructed by a block symbolic factorization which runs in time linear in matrix dimension. Numerical experiments illustrating these theoretical results are provided.
    Type of Medium: Electronic Resource
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 12
    Electronic Resource
    Electronic Resource
    Springer
    Numerische Mathematik 55 (1989), S. 667-684 
    ISSN: 0945-3245
    Keywords: AMS(MOS): 65F20 ; CR: G1.3
    Source: Springer Online Journal Archives 1860-2000
    Topics: Mathematics
    Notes: Summary We study the augmented system approach for the solution of sparse linear least-squares problems. It is well known that this method has better numerical properties than the method based on the normal equations. We use recent work by Arioli et al. (1988) to introduce error bounds and estimates for the components of the solution of the augmented system. In particular, we find that, using iterative refinement, we obtain a very robust algorithm and our estimates of the error are accurate and cheap to compute. The final error and all our error estimates are much better than the classical or Skeel's error analysis (1979) indicates. Moreover, we prove that our error estimates are independent of the row scaling of the augmented system and we analyze the influence of the Björck scaling (1967) on these estimates. We illustrate this with runs both on large-scale practical problems and contrived examples, comparing the numerical behaviour of the augmented systems approach with a code using the normal equations. These experiments show that while the augmented system approach with iterative refinement can sometimes be less efficient than the normal equations approach, it is comparable or better when the least-squares matrix has a full row, and is, in any case, much more stable and robust.
    Type of Medium: Electronic Resource
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 13
    Electronic Resource
    Electronic Resource
    Springer
    Numerische Mathematik 56 (1989), S. 215-227 
    ISSN: 0945-3245
    Keywords: AMS(MOS): 65F05 ; 65N20 ; CR: G1.3
    Source: Springer Online Journal Archives 1860-2000
    Topics: Mathematics
    Notes: Summary Discretization of the Poisson equation on a rectangle by finite differences using the standard five-point stencil yields a linear system of algebraic equations, which can be solved rapidly by the cyclic reduction method. In this method a sequence of tridiagonal linear systems is solved. The matrices of these systems commute, and we investigate numerical aspects of their ordering. In particular, we present two new ordering schemes that avoid overflow and loss of accuracy due to underflow. These ordering schemes improve the numerical performance of the subroutine HWSCRT of FISHPAK. Our orderings are also applicable to the solution of Helmholtz's equation by cyclic reduction, and to related numerical schemes, such as FACR methods.
    Type of Medium: Electronic Resource
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 14
    Electronic Resource
    Electronic Resource
    Springer
    Numerische Mathematik 56 (1989), S. 179-213 
    ISSN: 0945-3245
    Keywords: AMS(MOS): 65F10 ; 65H10 ; CR: G1.3 ; G1.5
    Source: Springer Online Journal Archives 1860-2000
    Topics: Mathematics
    Notes: Summary We first discuss the solution of a fixed point equationxΦx of a Fréchet differentiable self-mapping Φ by iterative methods of the general form $$x_m : = \left[ {\Phi \left( {\sum\limits_{i = 0}^{m - 1} {x_i \gamma _{i,m - 1} } } \right) - \sum\limits_{j = 0}^{m - 1} {x_j \beta _{j,m} } } \right]\frac{1}{{\beta _{m,m} }},m = 1,2, \ldots ,$$ defined by two infinite nonsingular upper triangular matricesB=(β j, m ) andC=(γ j, m ) with column sums 1. We show that two such methods, defined byB, C and $$\tilde B,\tilde C$$ , respectively, are in a certain sense equivalent if and only if $$C^{ - 1} B = \tilde C^{ - 1} \tilde B$$ . In particular, $$\hat B: = C^{ - 1} $$ and the unit matrix (replacingC) define an equivalent so-called semiiterative method. We introduce (k, l)-step methods as those whereB andC have upper bandwidthk andl−1, respectively. They require storing max {k, l} previous iterates only. For stationary methodsB andC have Toeplitz structure except for their first row, which is chosen such that the column sum condition holds. An Euler method, which may require to store all iterates, is equivalent to a (stationary) (k, l)-step method if and only if the underlying conformal mapg is a rational function of the formg(w)=w(μ 0+...+μ k −k )/ (v 0+...+v l−1 w −+1). By choosingg withg(1)=1 such that for some ρ〈1 it maps |w|〉ρ onto the exterior of some continuumS known to contain the eigen values of the Fréchet derivative of Φ, one obtains a feasible procedure for designing a locally converging stationary (k, l)-step method custom-made for a set of problems. In the case Φx≔Tx+d, with a linear operatorT, where one wants to solve a linear system of equations, we show that the residual polynomials of a stationary semiiterative method are generalized Faber polynomials with respect to a particular weight function. Using another weight function leads to what we call almost stationary methods. (The classical Chebyshev iteration is an example of such a method.) We define equivalent almost stationary (k, l)-step methods and give a corresponding convergence result.
    Type of Medium: Electronic Resource
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 15
    Electronic Resource
    Electronic Resource
    Springer
    Numerische Mathematik 56 (1989), S. 625-626 
    ISSN: 0945-3245
    Keywords: AMS (MOS): 65F15 ; CR: G1.3
    Source: Springer Online Journal Archives 1860-2000
    Topics: Mathematics
    Notes: Summary A-posteriori bound is given for the computed eigenpair ( $$\hat \lambda ,\hat x$$ ), of the eigenvalue problemAx=λx, which is shown to be more realistic than the available one. A simple expression is further presented for calculating the backward error.
    Type of Medium: Electronic Resource
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 16
    Electronic Resource
    Electronic Resource
    Springer
    Numerische Mathematik 52 (1988), S. 365-375 
    ISSN: 0945-3245
    Keywords: AMS(MOS): 65F15 ; CR: G1.3
    Source: Springer Online Journal Archives 1860-2000
    Topics: Mathematics
    Notes: Summary For each λ in some domainD in the complex plane, letF(λ) be a linear, compact operator on a Banach spaceX and letF be holomorphic in λ. Assuming that there is a ξ so thatI−F(ξ) is not one-to-one, we examine two local methods for approximating the nonlinear eigenvalue ξ. In the Newton method the smallest eigenvalue of the operator pencil [I−F(λ),F′(λ)] is used as increment. We show that under suitable hypotheses the sequence of Newton iterates is locally, quadratically convergent. Second, suppose 0 is an eigenvalue of the operator pencil [I−F(ξ),I] with algebraic multiplicitym. For fixed λ leth(λ) denote the arithmetic mean of them eigenvalues of the pencil [I−F(λ),I] which are closest to 0. Thenh is holomorphic in a neighborhood of ξ andh(ξ)=0. Under suitable hypotheses the classical Muller's method applied toh converges locally with order approximately 1.84.
    Type of Medium: Electronic Resource
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 17
    Electronic Resource
    Electronic Resource
    Springer
    Numerische Mathematik 52 (1988), S. 665-682 
    ISSN: 0945-3245
    Keywords: AMS(MOS): 65F05 ; CR: G1.3
    Source: Springer Online Journal Archives 1860-2000
    Topics: Mathematics
    Notes: Summary The paper deals with the problems of fast inversion of matricesA=T+H, whereT is Toeplitz andH is Hankel. Several algorithms are presented and compared, among them algorithms working for arbitrary strongly nonsingular matricesA=T+H.
    Type of Medium: Electronic Resource
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 18
    Electronic Resource
    Electronic Resource
    Springer
    Numerische Mathematik 55 (1989), S. 501-519 
    ISSN: 0945-3245
    Keywords: AMS(MOS): 65F50 ; 65K05 ; CR: G1.3
    Source: Springer Online Journal Archives 1860-2000
    Topics: Mathematics
    Notes: Summary We study the computation of sparse null bases of equilibrium matrices in the context of the force method in structural optimization. Two classes of structural problems are considered. For the class of rigid jointed skeletal structures, we use a partitioning method suggested by Henderson and Maunder to partition the problem into a set of independent null basis computations. For the class of structures represented by a continuum, we compute a sizable fraction of the null vectors in a basis from a consideration of the finite element formulation and the bipartite graph of the equilibrium matrix. The remaining null vectors are computed by the triangular algorithm in [6]. The new algorithms find sparser bases than the triangular algorithm and are also faster.
    Type of Medium: Electronic Resource
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 19
    Electronic Resource
    Electronic Resource
    Springer
    Numerische Mathematik 51 (1987), S. 143-180 
    ISSN: 0945-3245
    Keywords: AMS(MOS): 65N20 ; 65F10 ; 60K25 ; CR: G1.3
    Source: Springer Online Journal Archives 1860-2000
    Topics: Mathematics
    Notes: Summary Markovian queueing networks having overflow capacity are discussed. The Kolmogorov balance equations result in a linear homogeneous system, where the right null-vector is the steady-state probability distribution for the network. Preconditioned conjugate gradient methods are employed to find the null-vector. The preconditioner is a singular matrix which can be handled by separation of variables. The resulting preconditioned system is nonsingular. Numerical results show that the number of iterations required for convergence is roughly constant independent of the queue sizes. Analytic results are given to explain this fast convergence.
    Type of Medium: Electronic Resource
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 20
    Electronic Resource
    Electronic Resource
    Springer
    Numerische Mathematik 51 (1987), S. 429-439 
    ISSN: 0945-3245
    Keywords: AMS(MOS): 65F05 ; CR: G1.3
    Source: Springer Online Journal Archives 1860-2000
    Topics: Mathematics
    Notes: Summary In this paper we consider an extension to nonlinear algebraic systems of the class of algorithms recently proposed by Abaffy, Broyden and Spedicato for general linear systems. We analyze the convergence properties, showing that under the usual assumptions on the function and some mild assumptions on the free parameters available in the class, the algorithm is locally convergent and has a superlinear rate of convergence (per major iteration, which is computationally comparable to a single Newton's step). Some particular algorithms satisfying the conditions on the free parameters are considered.
    Type of Medium: Electronic Resource
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 21
    Electronic Resource
    Electronic Resource
    Springer
    Numerische Mathematik 53 (1988), S. 351-366 
    ISSN: 0945-3245
    Keywords: AMS(MOS)65F05 ; CR: G1.3
    Source: Springer Online Journal Archives 1860-2000
    Topics: Mathematics
    Notes: Summary New techniques for fast Toeplitz QR decomposition are presented. The methods are based on the shift invariance property of a Toeplitz matrix. The numerical properties of the algorithms are discussed and some comparisons are made with two other fast Toeplitz orthogonalization methods.
    Type of Medium: Electronic Resource
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 22
    Electronic Resource
    Electronic Resource
    Springer
    Numerische Mathematik 53 (1988), S. 485-498 
    ISSN: 0945-3245
    Keywords: AMS (MOS): 65F10 ; CR: G1.3
    Source: Springer Online Journal Archives 1860-2000
    Topics: Mathematics
    Notes: Summary Recently, special attention has been given in the literature, to the problems of accurately computing the least-squares solution of very largescale over-determined systems of linear equations which occur in geodetic applications. In particular, it has been suggested that one can solve such problems iteratively by applying the block SOR (Successive Overrelaxation) and EGS1 (Extrapolated Gauss Seidel 1) plus semi-iterative methods to a linear system with coefficient matrix 2-cyclic or 3-cyclic. The comparison of 2-block SOR and 3-block SOR was made in [1] and showed that the 2-block SOR is better. In [6], the authors also proved that 3-block EGS1-SI is better than 3-block SOR. Here, we first show that the 2-block DJ (Double Jacobi)-SI, GS-SI and EGS1-SI methods are equivalent and all of them are equivalent to the 3-block EGS1-SI method; then, we prove that the composite methods and 2-block SOR have the same asymptotic rate of convergence, but the former has a better average rate of convergence; finally, numerical experiments are reported, and confirm that the 3-block EGS1-SI is better than the 2-block SOR.
    Type of Medium: Electronic Resource
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 23
    Electronic Resource
    Electronic Resource
    Springer
    Numerische Mathematik 47 (1985), S. 483-504 
    ISSN: 0945-3245
    Keywords: AMS(MOS) ; 65F30 ; CR: G1.3
    Source: Springer Online Journal Archives 1860-2000
    Topics: Mathematics
    Notes: Summary LetA be a realm×n matrix with full row rankm. In many algorithms in engineering and science, such as the force method in structural analysis, the dual variable method for the Navier-Stokes equations or more generally null space methods in quadratic programming, it is necessary to compute a basis matrixB for the null space ofA. HereB isn×r, r=n−m, of rankr, withAB=0. In many instancesA is large and sparse and often banded. The purpose of this paper is to describe and test a variation of a method originally suggested by Topcu and called the turnback algorithm for computing a banded basis matrixB. Two implementations of the algorithm are given, one using Gaussian elimination and the other using orthogonal factorization by Givens rotations. The FORTRAN software was executed on an IBM 3081 computer with an FPS-164 attached array processor at the Triangle Universities Computing Center and on a CYBER 205 vector computer. Test results on a variety of structural analysis problems including two- and three-dimensional frames, plane stress, plate bending and mixed finite element problems are discussed. These results indicate that both implementations of the algorithm yielded a well-conditioned, banded, basis matrixB whenA is well-conditioned. However, the orthogonal implementation yielded a better conditionedB for large, illconditioned problems.
    Type of Medium: Electronic Resource
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 24
    Electronic Resource
    Electronic Resource
    Springer
    Numerische Mathematik 50 (1987), S. 613-632 
    ISSN: 0945-3245
    Keywords: AMS(MOS): 65G05, 65F05 ; CR: G1.3
    Source: Springer Online Journal Archives 1860-2000
    Topics: Mathematics
    Notes: Summary A forward error analysis is presented for the Björck-Pereyra algorithms used for solving Vandermonde systems of equations. This analysis applies to the case where the points defining the Vandermonde matrix are nonnegative and are arranged in increasing order. It is shown that for a particular class of Vandermonde problems the error bound obtained depends on the dimensionn and on the machine precision only, being independent of the condition number of the coefficient matrix. By comparing appropriate condition numbers for the Vandermonde problem with the forward error bounds it is shown that the Björck-Pereyra algorithms introduce no more uncertainty into the numerical solution than is caused simply by storing the right-hand side vector on the computer. A technique for computing “running” a posteriori error bounds is derived. Several numerical experiments are presented, and it is observed that the ordering of the points can greatly affect the solution accuracy.
    Type of Medium: Electronic Resource
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 25
    Electronic Resource
    Electronic Resource
    Springer
    Numerische Mathematik 56 (1989), S. 73-91 
    ISSN: 0945-3245
    Keywords: AMS(MOS): 65F15 ; CR: G1.3
    Source: Springer Online Journal Archives 1860-2000
    Topics: Mathematics
    Notes: Summary Linear convergence of the row cyclic Jacobi and Kogbetliantz methods can be guaranteed if certain constraints concerning the angles of rotations are implemented. Unlike the results of Forsythe and Henrici, convergence can be achieved without under-rotations.
    Type of Medium: Electronic Resource
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 26
    Electronic Resource
    Electronic Resource
    Springer
    Numerische Mathematik 56 (1989), S. 157-177 
    ISSN: 0945-3245
    Keywords: AMS(MOS): 65F10 ; 65N20 ; 65N30 ; CR: G1.3
    Source: Springer Online Journal Archives 1860-2000
    Topics: Mathematics
    Notes: Summary A recursive way of constructing preconditioning matrices for the stiffness matrix in the discretization of selfadjoint second order elliptic boundary value problems is proposed. It is based on a sequence of nested finite element spaces with the usual nodal basis functions. Using a nodeordering corresponding to the nested meshes, the finite element stiffness matrix is recursively split up into two-level block structures and is factored approximately in such a way that any successive Schur complement is replaced (approximated) by a matrix defined recursively and thereform only implicitely given. To solve a system with this matrix we need to perform a fixed number (v) of iterations on the preceding level using as an iteration matrix the preconditioning matrix already defined on that level. It is shown that by a proper choice of iteration parameters it suffices to use $$v 〉 \left( {1 - \gamma ^2 } \right)^{ - \tfrac{1}{2}} $$ iterations for the so constructedv-foldV-cycle (wherev=2 corresponds to aW-cycle) preconditioning matrices to be spectrally equivalent to the stiffness matrix. The conditions involve only the constant λ in the strengthened C.-B.-S. inequality for the corresponding two-level hierarchical basis function spaces and are therefore independent of the regularity of the solution for instance. If we use successive uniform refinements of the meshes the method is of optimal order of computational complexity, if $$\gamma ^2〈 \tfrac{8}{9}$$ . Under reasonable assumptions of the finite element mesh, the condition numbers turn out to be so small that there are in practice few reasons to use an accelerated iterative method like the conjugate gradient method, for instance.
    Type of Medium: Electronic Resource
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 27
    Electronic Resource
    Electronic Resource
    Springer
    Numerische Mathematik 56 (1989), S. 283-289 
    ISSN: 0945-3245
    Keywords: AMS(MOS): 65F10 ; CR: G1.3
    Source: Springer Online Journal Archives 1860-2000
    Topics: Mathematics
    Notes: Summary Comparison results for weak regular splittings of monotone matrices are derived. As an application we get upper and lower bounds for the convergence rate of iterative procedures based on multisplittings. This yields a very simple proof of results of Neumann-Plemmons on upper bounds, and establishes lower bounds, which has in special cases been conjectured by these authors.
    Type of Medium: Electronic Resource
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 28
    Electronic Resource
    Electronic Resource
    Springer
    Numerische Mathematik 56 (1989), S. 269-282 
    ISSN: 0945-3245
    Keywords: AMS(MOS): 65H10 ; 65W05 ; CR: G1.3
    Source: Springer Online Journal Archives 1860-2000
    Topics: Mathematics
    Notes: Summary Linear multisplitting methods are known as parallel iterative methods for solving a linear systemAx=b. We extend the idea of multisplittings to the problem of solving a nonlinear system of equationsF(x)=0. Our nonlinear multisplittings are based on several nonlinear splittings of the functionF. In a parallel computing environment, each processor would have to calculate the exact solution of an individual nonlinear system belonging to ‘his’ nonlinear multisplitting and these solutions are combined to yield the next iterate. Although the individual systems are usually much less involved than the original system, the exact solutions will in general not be available. Therefore, we consider important variants where the exact solutions of the individual systems are approximated by some standard method such as Newton's method. Several methods proposed in literature may be regarded as special nonlinear multisplitting methods. As an application of our systematic approach we present a local convergence analysis of the nonlinear multisplitting methods and their variants. One result is that the local convergence of these methods is determined by an induced linear multisplitting of the Jacobian ofF.
    Type of Medium: Electronic Resource
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 29
    Electronic Resource
    Electronic Resource
    Springer
    Numerische Mathematik 46 (1985), S. 365-395 
    ISSN: 0945-3245
    Keywords: AMS(MOS) 65F, 65G ; CR: G1.3
    Source: Springer Online Journal Archives 1860-2000
    Topics: Mathematics
    Notes: Summary Part I of this work deals with the forward error analysis of Gaussian elimination for general linear algebraic systems. The error analysis is based on a linearization method which determines first order approximations of the absolute errors exactly. Superposition and cancellation of error effects, structure and sparsity of the coefficient matrices are completely taken into account by this method. The most important results of the paper are new condition numbers and associated optimal component-wise error and residual estimates for the solutions of linear algebraic systems under data perturbations and perturbations by rounding erros in the arithmetic floating-point operations. The estimates do not use vector or matrix norms. The relative data and rounding condition numbers as well as the associated backward and residual stability constants are scaling-invariant. The condition numbers can be computed approximately from the input data, the intermediate results, and the solution of the linear system. Numerical examples show that by these means realistic bounds of the errors and the residuals of approximate solutions can be obtained. Using the forward error analysis, also typical results of backward error analysis are deduced. Stability theorems and a priori error estimates for special classes of linear systems are proved in Part II of this work.
    Type of Medium: Electronic Resource
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 30
    Electronic Resource
    Electronic Resource
    Springer
    Numerische Mathematik 46 (1985), S. 69-83 
    ISSN: 0945-3245
    Keywords: AMS(MOS): 65F30 ; CR: G1.3
    Source: Springer Online Journal Archives 1860-2000
    Topics: Mathematics
    Notes: Summary Let $$\mathfrak{A}$$ be a real irreduciblen×n interval matrix. Then a necessary and sufficient condition is given for the sequence $$\{ \mathfrak{A}^k \} $$ of the powers of an interval matrix $$\mathfrak{A}$$ to converge to a matrix $$\mathfrak{A}^\infty $$ which is not the null matrix. In addition a criterion for $$\mathfrak{A}$$ is proved to decide whether the limit matrix $$\mathfrak{A}^\infty $$ satisfies the condition of symmetry $$\mathfrak{A}^\infty = - \mathfrak{A}^\infty $$ .
    Type of Medium: Electronic Resource
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 31
    Electronic Resource
    Electronic Resource
    Springer
    Numerische Mathematik 48 (1986), S. 349-368 
    ISSN: 0945-3245
    Keywords: AMS(MOS): 65F15 ; CR: G1.3
    Source: Springer Online Journal Archives 1860-2000
    Topics: Mathematics
    Notes: Summary The Symmetric Tridiagonal Eigenproblem has been the topic of some recent work. Many methods have been advanced for the computation of the eigenvalues of such a matrix. In this paper, we present a divide-and-conquer approach to the computation of the eigenvalues of a symmetric tridiagonal matrix via the evaluation of the characteristic polynomial. The problem of evaluation of the characteristic polynomial is partitioned into smaller parts which are solved and these solutions are then combined to form the solution to the original problem. We give the update equations for the characteristic polynomial and certain auxiliary polynomials used in the computation. Furthermore, this set of recursions can be implemented on a regulartree structure. If the concurrency exhibited by this algorithm is exploited, it can be shown that thetime for computation of all the eigenvalues becomesO(nlogn) instead ofO(n 2) as is the case for the approach where the order is increased by only one at every step. We address the numerical problems associated with the use of the characteristic polynomial and present a numerically stable technique for the eigenvalue computation.
    Type of Medium: Electronic Resource
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 32
    ISSN: 0945-3245
    Keywords: AMS(MOS): 65F10 ; CR: G1.3
    Source: Springer Online Journal Archives 1860-2000
    Topics: Mathematics
    Notes: Summary Ann×n complex matrixB is calledparacontracting if ‖B‖2≦1 and 0≠x∈[N(I-B)]⊥⇒‖Bx‖2〈‖x‖2. We show that a productB=B k B k−1 ...B 1 ofk paracontracting matrices is semiconvergent and give upper bounds on the subdominant eigenvalue ofB in terms of the subdominant singular values of theB i 's and in terms of the angles between certain subspaces. Our results here extend earlier results due to Halperin and due to Smith, Solomon and Wagner. We also determine necessary and sufficient conditions forn numbers in the interval [0, 1] to form the spectrum of a product of two orthogonal projections and hence characterize the subdominant eigenvalue of such a product. In the final part of the paper we apply the upper bounds mentioned earlier to provide an estimate on the subdominant eigenvalue of the SOR iteration matrix ℒω associated with ann×n hermitian positive semidefinite matrixA none of whose diagonal entries vanish.
    Type of Medium: Electronic Resource
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 33
    Electronic Resource
    Electronic Resource
    Springer
    Numerische Mathematik 46 (1985), S. 339-349 
    ISSN: 0945-3245
    Keywords: AMS(MOS): 65F05, 65B05 ; CR: G1.3
    Source: Springer Online Journal Archives 1860-2000
    Topics: Mathematics
    Description / Table of Contents: Zusammenfassung In dieser Arbeit wird die folgende Fragestellung behandelt: Wenn man weiß, daß gewisse TeilsystemeA 1·x 1=b 1, ...,A m ·x m =b m eines linearen GleichungssystemsA·x=b eindeutig lösbar sind, was kann man über die Lösbarkeit des SystemsA·x=b sagen und wie kann man gegebenenfallsx ausx 1,...,x m berechnen? Die Lösung dieses Problems erlaubt die rekursive Berechnung gewisser linearer oder quasilinearer Folgentransformationen [7, 13, 15].
    Notes: Summary This note is concerned with the following problem: Given a systemA·x=b of linear equations and knowing that certains of its subsystemsA 1·x 1=b 1, ...,A m ·x m =b m can be solved uniquely what can be said about the regularity ofA and how to find the solutionx fromx 1, ...,x m ? This question is of particular interest for establishing methods computing certain linear or quasilinear sequence transformations recursively [7, 13, 15].
    Type of Medium: Electronic Resource
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 34
    Electronic Resource
    Electronic Resource
    Springer
    Numerische Mathematik 46 (1985), S. 255-268 
    ISSN: 0945-3245
    Keywords: AMS(MOS) 65F05 ; 65F25 ; 65G99 ; CR: G1.3
    Source: Springer Online Journal Archives 1860-2000
    Topics: Mathematics
    Notes: Summary The accumulation of rounding errors in a method used to compute the solution of an underdetermined system of linear equations at the least distance from a given point is being studied. The method used is based on orthogonal matrix decompositions.
    Type of Medium: Electronic Resource
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 35
    Electronic Resource
    Electronic Resource
    Springer
    Numerische Mathematik 46 (1985), S. 443-454 
    ISSN: 0945-3245
    Keywords: AMS(MOS): 65F15 ; CR: G1.3
    Source: Springer Online Journal Archives 1860-2000
    Topics: Mathematics
    Notes: Summary The DAP architecture brings into consideration the Jacobi method where several non-interacting rotations can be performed in parallel. However the design of the algorithm is crucial in a parallel environment. In this paper we shall consider two techniques specifically designed to reduce the organisation and the arithmetic components of the parallel Jacobi method.
    Type of Medium: Electronic Resource
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 36
    Electronic Resource
    Electronic Resource
    Springer
    Numerische Mathematik 46 (1985), S. 479-491 
    ISSN: 0945-3245
    Keywords: AMS(MOS):65F30 ; CR: G1.3
    Source: Springer Online Journal Archives 1860-2000
    Topics: Mathematics
    Notes: Summary If the columns of a matrix are orthonormal and it is partitioned into a 2-by-1 block matrix, then the singular value decompositions of the blocks are related. This is the essence of the “CS decomposition”. The computation of these related SVD's requires some care. Stewart has given an algorithm that uses the LINPACK SVD algorithm together with a Jacobitype “clean-up” operation on a cross-product matrix. Our technique is equally stable and fast but avoids the cross product matrix. The simplicity of our technique makes it more amenable to parallel computation on systolic-type computer architectures. These developments are of interest because a good way to compute the generalized singular value decomposition of a matrix pair (A, B) is to compute the CS decomposition of a certain orthogonal column matrix related toA andB.
    Type of Medium: Electronic Resource
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 37
    Electronic Resource
    Electronic Resource
    Springer
    Numerische Mathematik 48 (1986), S. 239-249 
    ISSN: 0945-3245
    Keywords: AMS (MOS): 65F05 ; CR: G1.3
    Source: Springer Online Journal Archives 1860-2000
    Topics: Mathematics
    Notes: Summary We show that the greedy algorithm introduced in [1] and [5] to perform the parallel QR decomposition of a dense rectangular matrix of sizem×n is optimal. Then we assume thatm/n 2 tends to zero asm andn go to infinity, and prove that the complexity of such a decomposition is asymptotically2n, when an unlimited number of processors is available.
    Type of Medium: Electronic Resource
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 38
    Electronic Resource
    Electronic Resource
    Springer
    Numerische Mathematik 48 (1986), S. 525-542 
    ISSN: 0945-3245
    Keywords: AMS(MOS): 65F10, 41A50, 41A17 ; CR: G1.3
    Source: Springer Online Journal Archives 1860-2000
    Topics: Mathematics
    Notes: Summary We investigate approximation problems of the form (P) $$\mathop {\min }\limits_{p:p(0) = 1} \mathop {\max |p(\lambda )|}\limits_{\lambda \in S} $$ withp ranging over all polynomials of degree ≦k andS any line segment of the complex plane. Problems of this type appear in connection with error bounds for the minimum residual algorithm for solving linear systems with coefficient matrixA=e iθ(dI−N), whereN=−N H andd∈ℝ. Besides a few general results on (P), we present the explicit solution of (P) for segmentsS which are parallel to the imaginary and symmetric to the real axis. This result leads to an optimal semi-iterative method for the class of matricesA=I−N, whereN=−N T.
    Type of Medium: Electronic Resource
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 39
    Electronic Resource
    Electronic Resource
    Springer
    Numerische Mathematik 50 (1986), S. 17-26 
    ISSN: 0945-3245
    Keywords: AMS(MOS): 65F10 ; 65G10 ; CR: G1.3
    Source: Springer Online Journal Archives 1860-2000
    Topics: Mathematics
    Notes: Summary In classical numerical analysis the asymptotic convergence factor (R 1-factor) of an iterative processx m+1=Axm+b coincides with the spectral radius of then×n iteration matrixA. Thus the famous Theorem of Stein and Rosenberg can at least be partly reformulated in terms of asymptotic convergence factor. Forn×n interval matricesA with irreducible upper bound and nonnegative lower bound we compare the asymptotic convergence factor (α T ) of the total step method in interval analysis with the factorα S of the corresponding single step method. We derive a result similar to that of the Theorem of Stein and Rosenberg. Furthermore we show thatα S can be less than the spectral radius of the real single step matrix corresponding to the total step matrix |A| where |A| is the absolute value ofA. This answers an old question in interval analysis.
    Type of Medium: Electronic Resource
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 40
    Electronic Resource
    Electronic Resource
    Springer
    Numerische Mathematik 47 (1985), S. 175-190 
    ISSN: 0945-3245
    Keywords: AMS (MOS): 65F05 ; 65N30 ; CR: G1.3
    Source: Springer Online Journal Archives 1860-2000
    Topics: Mathematics
    Description / Table of Contents: Summary A ≪nested dissection≫ ordering is given for solving any system of linear equationsAX=B, for the family ofsparse symmetric positive definite matrices corresponding to the class of undirected graphs withbounded degree, and satisfyinga $$\sqrt n $$ -separator theorem. If we compare the methods and results presented by Rose [9], our algorithm is more simple, but the complexity results are less general because of the restriction of bounded degree. Besides, our proofs use continually the arborescent structure on the family of separators, which is, by another way, a partition of the set of vertices for the initial graph. Then, the general implementation scheme in the finite element package MODULEF, fortwo-dimensional finite element problems, is presented, and numerical comparisons between our ordering and the standard envelope method are given.
    Notes: Résumé Nous présentons une numérotation de type ≪Nested Dissection≫ des inconnues d'un système linéaireAX=B pour des ensembles de matricescreuses symétriques définies positives correspondant à des famille de graphes non orientés,à degré borné, et admettant un $$\sqrt n $$ -thérème de séparation. Comparativement aux méthodes et résultats de Rose [9], l'algorithme présenté est plus simple, mais les théorèmes de complexité moins généraux, en raison de l'hypothèse restrictive de degré borné. En outre, les démonstrations font appel en permanence à la structure d'arbre sur la famille des séparateurs qui constitue, par ailleurs, une partition de l'ensemble des sommets du graphe initial. Nous présentons ensuite le schéma général d'implémentation dans le cadre du code d'éléments finis MODULEF pourdes problèmes plans d'éléments finis, et nous donnons quelques mesures comparatives avec la numérotation plus classique qui tend à minimiser le profil de la matrice.
    Type of Medium: Electronic Resource
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 41
    Electronic Resource
    Electronic Resource
    Springer
    Numerische Mathematik 49 (1986), S. 81-94 
    ISSN: 0945-3245
    Keywords: AMS(MOS): 65F05 ; CR: G1.3
    Source: Springer Online Journal Archives 1860-2000
    Topics: Mathematics
    Notes: Summary This paper presents a new algorithm for computing theQR factorization of anm×n Toeplitz matrix inO(mn) operations. The algorithm exploits the procedure for the rank-1 modification and the fact that both principal (m−1)×(n−1) submatrices of the Toeplitz matrix are identical. An efficient parallel implementation of the algorithm is possible.
    Type of Medium: Electronic Resource
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 42
    Electronic Resource
    Electronic Resource
    Springer
    Numerische Mathematik 48 (1986), S. 499-523 
    ISSN: 0945-3245
    Keywords: AMS(MOS): 65F10 ; CR: G1.3
    Source: Springer Online Journal Archives 1860-2000
    Topics: Mathematics
    Notes: Summary We derive new estimates for the rate of convergence of the conjugate gradient method by utilizing isolated eigenvalues of parts of the spectrum. We present a new generalized version of an incomplete factorization method and compare the derived estimates of the number of iterations with the number actually found for some elliptic difference equations and for a similar problem with a model empirical distribution function.
    Type of Medium: Electronic Resource
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 43
    Electronic Resource
    Electronic Resource
    Springer
    Journal of mathematical biology 23 (1985), S. 75-101 
    ISSN: 1432-1416
    Keywords: Population dynamics ; Optimal harvesting ; Optimal control ; Hyperbolic systems ; Pontryagin's principle
    Source: Springer Online Journal Archives 1860-2000
    Topics: Biology , Mathematics
    Notes: Abstract In this paper, Pontryagin's principle is proved for a fairly general problem of optimal control of populations with continuous time and age variable. As a consequence, maximum principles are developed for an optimal harvesting problem and a problem of optimal birth control.
    Type of Medium: Electronic Resource
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 44
    Electronic Resource
    Electronic Resource
    Springer
    Journal of optimization theory and applications 45 (1985), S. 265-293 
    ISSN: 1573-2878
    Keywords: Optimal control ; control differential systems involving impulses ; jumps in the state variables
    Source: Springer Online Journal Archives 1860-2000
    Topics: Mathematics
    Notes: Abstract The existence of an optimal control is proved in a problem where the criterion functional and the equation of motion contain a control that is a Lebesgue-Stieltjes measure. Due to nonlinearities in the problem, it is necessary to postulate a condition implying that large atoms of the control measures are sparse and that their derivatives, wherever they exist, have a uniform bound.
    Type of Medium: Electronic Resource
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 45
    Electronic Resource
    Electronic Resource
    Springer
    Journal of optimization theory and applications 56 (1988), S. 227-243 
    ISSN: 1573-2878
    Keywords: Optimal control ; nondifferentiable control problems ; optimality conditions ; Fritz John optimality conditions ; duality ; converse duality ; continuous programming ; mathematical programming
    Source: Springer Online Journal Archives 1860-2000
    Topics: Mathematics
    Notes: Abstract Optimality conditions and duality results are obtained for a class of control problems having a nondifferentiable term in the integrand of the objective functional. These results generalize many well-known results in optimal control theory involving differentiable functions, and also provide a relationship with certain nondifferentiable mathematical programming problems. Some extensions concerning the unified treatment of optimal control theory and continuous programming are also mentioned. Finally, a control problem containing an arbitrary norm, along with its appropriate norm, is given.
    Type of Medium: Electronic Resource
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 46
    Electronic Resource
    Electronic Resource
    Springer
    Journal of optimization theory and applications 53 (1987), S. 429-449 
    ISSN: 1573-2878
    Keywords: Optimal control ; discrete systems ; convex problems ; state and control constraints ; sensitivity analysis ; right-differentiability with respect to a parameter
    Source: Springer Online Journal Archives 1860-2000
    Topics: Mathematics
    Notes: Abstract A family of optimal control problems for discrete systems that depend on a real parameter is considered. The problems are strongly convex and subject to state and control constraints. Some regularity conditions are imposed on the constraints. The control problems are reformulated as mathematical programming problems. It is shown that both the primal and dual optimal variables for these problems are right-differentiable functions of a parameter. The right-derivatives are characterized as solutions to auxiliary quadratic control problems. Conditions of continuous differentiability are discussed, and some estimates of the rate of convergence of the difference quotients to the respective derivatives are given.
    Type of Medium: Electronic Resource
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 47
    Electronic Resource
    Electronic Resource
    Springer
    Journal of optimization theory and applications 45 (1985), S. 517-531 
    ISSN: 1573-2878
    Keywords: Optimal control ; parameter optimization ; power systems ; load frequency control
    Source: Springer Online Journal Archives 1860-2000
    Topics: Mathematics
    Notes: Abstract A new approach for designing a linear regulator for the problem of load frequency control (LFC) of interconnected power systems is developed. The control is specified to be of proportional-plus-integral (P-I) form and is only a function of the measurable states. The LFC problem is formulated as a parameter optimization problem.
    Type of Medium: Electronic Resource
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 48
    Electronic Resource
    Electronic Resource
    Springer
    Journal of optimization theory and applications 57 (1988), S. 513-517 
    ISSN: 1573-2878
    Keywords: Optimal control ; population growth ; logistic equation
    Source: Springer Online Journal Archives 1860-2000
    Topics: Mathematics
    Notes: Abstract The rate at which a population should grow is determined by finding the best trade-off between the loss due to the deviation from a target population size and the loss associated to the growing effort. It is also shown that, in the case of infinite-time horizon and quadratic loss functions, the optimal growth is logistic.
    Type of Medium: Electronic Resource
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 49
    Electronic Resource
    Electronic Resource
    Springer
    Journal of optimization theory and applications 46 (1985), S. 493-504 
    ISSN: 1573-2878
    Keywords: Optimal control ; envelope theorem ; differential inclusions ; nondifferentiable functions ; production-employment model
    Source: Springer Online Journal Archives 1860-2000
    Topics: Mathematics
    Notes: Abstract In optimal control problems involving nondifferentiable functions of the state variable, the adjoint differential inclusion can be formulated by either use of the Hamiltonian or the maximized Hamiltonian. In this paper, we solve a production-employment model in which the latter approach must be utilized, since the former does not enable one to determine the optimal policy.
    Type of Medium: Electronic Resource
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 50
    Electronic Resource
    Electronic Resource
    Springer
    Journal of optimization theory and applications 59 (1988), S. 445-465 
    ISSN: 1573-2878
    Keywords: Optimal control ; linear-quadratic systems ; unknown targets ; points of information
    Source: Springer Online Journal Archives 1860-2000
    Topics: Mathematics
    Notes: Abstract The situation considered is of optimally controlling a deterministic system from a given state to an initially unknown targety in a fixed time interval [T 0,T]. The target will be certainly known at a random time τ in [T 0,T]. The controller knows the distributions ofy and τ. We derive the Bellman equation for the problem, prove a verification theorem for it, and demonstrate how the distribution τ influences the optimal control. We show that, in the linear-quadratic case, the optimal control is given by a feedback law that does not depend on the distribution of τ.
    Type of Medium: Electronic Resource
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 51
    Electronic Resource
    Electronic Resource
    Springer
    Journal of optimization theory and applications 45 (1985), S. 73-88 
    ISSN: 1573-2878
    Keywords: Optimal control ; nonconvexity ; relaxed controls ; parabolic systems ; relaxed minimum principle ; approximations ; descent methods
    Source: Springer Online Journal Archives 1860-2000
    Topics: Mathematics
    Notes: Abstract In this paper, we consider an optimal control problem for distributed systems governed by parabolic equations. The state equations are nonlinear in the control variable; the constraints and the cost functional are generally nonconvex. Relaxed controls are used to prove existence and derive necessary conditions for optimality. To compute optimal controls, a descent method is applied to the resulting relaxed problem. A numerical method is also given for approximating a special class of relaxed controls, notably those obtained by the descent method. Convergence proofs are given for both methods, and a numerical example is provided.
    Type of Medium: Electronic Resource
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 52
    Electronic Resource
    Electronic Resource
    Springer
    Journal of optimization theory and applications 45 (1985), S. 505-516 
    ISSN: 1573-2878
    Keywords: Optimal control ; parameter optimization ; power systems ; load frequency control
    Source: Springer Online Journal Archives 1860-2000
    Topics: Mathematics
    Notes: Abstract A new technique is described for designing an optimal controller for a system whose dynamical equations contain a backlash element. The approach is applied to the problem of load frequency control (LFC) of a single area steam power system.
    Type of Medium: Electronic Resource
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 53
    Electronic Resource
    Electronic Resource
    Springer
    Journal of optimization theory and applications 46 (1985), S. 441-453 
    ISSN: 1573-2878
    Keywords: Optimal control ; singular perturbations ; autonomous systems ; feedback control ; Hamilton-Jacobi-Bellmann equation
    Source: Springer Online Journal Archives 1860-2000
    Topics: Mathematics
    Notes: Abstract This paper shows how to construct a feedback control law for a class of singularly perturbed autonomous optimization problems. The control law is expressed as a single power series in the small parameter ∈ representing the ratio of the two effective time scales of the problem. The present approach avoids the need of expansion matching. The method is applied to a constant-speed interception problem. Comparison of numerical results with the exact solution shows an excellent agreement.
    Type of Medium: Electronic Resource
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 54
    Electronic Resource
    Electronic Resource
    Springer
    Journal of optimization theory and applications 47 (1985), S. 465-481 
    ISSN: 1573-2878
    Keywords: Optimal control ; tracking systems ; computational methods
    Source: Springer Online Journal Archives 1860-2000
    Topics: Mathematics
    Notes: Abstract Closed-from solutions are derived for a class of tracking problems including a linear optimal regulator and a prefilter for a time-invariant plant. The solutions for the prefilter equation and state trajectory, coupled by the Riccati equation, are exponentially related to the stability matrix of the plant. A computational procedure is presented in recursive form when the desired output state dynamics is assumed linear and time-invariant. Several examples are given for illustration.
    Type of Medium: Electronic Resource
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 55
    Electronic Resource
    Electronic Resource
    Springer
    Journal of optimization theory and applications 48 (1986), S. 265-287 
    ISSN: 1573-2878
    Keywords: Optimal control ; Hamilton-Jacobi theory ; infinite horizon ; equivalent variational problems
    Source: Springer Online Journal Archives 1860-2000
    Topics: Mathematics
    Notes: Abstract In this paper, we extend Carathéodory's concept of equivalent variational problems to infinite-horizon optimal control problems. In such a setting, the usual concept of a minimum need not exist, and we therefore consider a weaker definition of optimality, known as catching up optimality. The extension presented here leads us to a Hamilton-Jacobi theory for infinite-horizon optimal control problems that closely parallels the classical work of Carathéodory as well as providing sufficient conditions for optimality. Finally, we show that the results given here subsume several previously known results as a special case.
    Type of Medium: Electronic Resource
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 56
    Electronic Resource
    Electronic Resource
    Springer
    Journal of optimization theory and applications 48 (1986), S. 289-302 
    ISSN: 1573-2878
    Keywords: Optimal control ; Pontryagin's maximum principle ; optimization of liver kinetics
    Source: Springer Online Journal Archives 1860-2000
    Topics: Mathematics
    Notes: Abstract The enzyme activities in the liver are described by a system of ordinary differential equations. A certain substance in the blood is transformed twice by different kinds of enzymes. The mathematical problem is to determine the distribution of the enzymes along the blood flow so that the outflow concentration of the once-transformed form of the substance is as small as possible. This problem has been considered before and solved for particular types of enzyme kinetics. In this paper, we solve the problem for more general types of kinetics (including substrate- inhibition kinetics). The methods used are also different, in that the problem is considered as a problem of optimal control and Pontryagin's maximum principle is applied to derive necessary conditions.
    Type of Medium: Electronic Resource
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 57
    Electronic Resource
    Electronic Resource
    Springer
    Journal of optimization theory and applications 49 (1986), S. 207-225 
    ISSN: 1573-2878
    Keywords: Optimal control ; existence theory ; catching-up optimal solutions ; infinite-horizon problems
    Source: Springer Online Journal Archives 1860-2000
    Topics: Mathematics
    Notes: Abstract In this paper, we are concerned with the question of the existence of optimal solutions for infinite-horizon optimal control problems of Lagrange type. In such problems, the objective or cost functional is described by an improper integral. As dictated by applications arising in mathematical economics, we do nota priori assume that this improper integral converges. This leads us to consider a weaker type of optimality, known as catching-up optimality. The results presented here utilize the classical convexity and seminormality conditions typically imposed in the existence theory for the case of finite intervals. These conditions are significantly weaker than those imposed by other authors; as a consequence, their existence results are contained as special cases of the results presented here. The method of proof utilizes the Carathéodory-Hamilton-Jacobi theory previously developed by the author for infinite-horizon optimal control problems.
    Type of Medium: Electronic Resource
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 58
    Electronic Resource
    Electronic Resource
    Springer
    Journal of optimization theory and applications 49 (1986), S. 449-461 
    ISSN: 1573-2878
    Keywords: Optimal control ; multireservoir hydroelectric power systems ; optimal release policies ; series reservoirs
    Source: Springer Online Journal Archives 1860-2000
    Topics: Mathematics
    Notes: Abstract The optimal monthly operating policy of a multireservoir hydroelectric power system is a stochastic nonlinear problem. This paper presents a new method for determining the optimal monthly operating policy of a power system ofn reservoirs in series on a river taking into account the stochasticity of the river flows. Functional optimization techniques and minimum-norm formulation have been used to find the optimal release policy of the system. Results for a numerical example composed of four reservoirs are presented.
    Type of Medium: Electronic Resource
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 59
    Electronic Resource
    Electronic Resource
    Springer
    Journal of optimization theory and applications 50 (1986), S. 383-395 
    ISSN: 1573-2878
    Keywords: Optimal control ; multireservoir hydroelectric power systems ; optimal release policy ; parallel reservoirs
    Source: Springer Online Journal Archives 1860-2000
    Topics: Mathematics
    Notes: Abstract Long-term optimal operation of a multireservoir system is complex because it is a dynamic problem (present decisions for one reservoir depend on future decisions for all reservoirs); the optimal operating policy for one reservoir depends not only on its own energy content, but also on the corresponding content of each one of the other reservoirs; it is a highly stochastic problem with respect to the reservoir inflows and it is a nonlinear problem. This paper presents a new method for determining the optimal monthly operating policy of a power system consisting of multireservoirs on a multiriver system taking into account the stochasticity of the river flows. Functional optimization techniques and minimum norm formulation have been used. Results for a numerical example composed of three rivers with four reservoirs, three reservoirs, and two reservoirs on each river, respectively, are presented.
    Type of Medium: Electronic Resource
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 60
    Electronic Resource
    Electronic Resource
    Springer
    Journal of optimization theory and applications 50 (1986), S. 463-477 
    ISSN: 1573-2878
    Keywords: Optimal control ; multireservoir hydroelectric power system ; optimal release policy ; variable head ; functional analysis
    Source: Springer Online Journal Archives 1860-2000
    Topics: Mathematics
    Notes: Abstract In this paper, functional analysis and minimum norm formulation are applied to maximize the total benefits from two hydro reservoirs. The hydroelectric power generation is treated as a nonlinear function; water head variation and stochasticity of the river flows are included. The resulting problem has a nonlinear objective function and linear constraints. The proposed method is computationally efficient, compared to previous techniques. Numerical results are presented for widely different water conditions for an actual system in operation.
    Type of Medium: Electronic Resource
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 61
    Electronic Resource
    Electronic Resource
    Springer
    Journal of optimization theory and applications 51 (1986), S. 41-62 
    ISSN: 1573-2878
    Keywords: Optimal control ; existence theory ; finite optimality ; infinite horizons ; continuous dependence
    Source: Springer Online Journal Archives 1860-2000
    Topics: Mathematics
    Notes: Abstract In this paper, we investigate the existence of finitely optimal solutions for the Lagrange problem of optimal control defined on [0, ∞) under weaker convexity and seminormality hypotheses than those of previous authors. The notion of finite optimality has been introduced into the literature as the weakest of a hierarchy of types of optimality that have been defined to permit the study of Lagrange problems, arising in mathematical economics, whose cost functions either diverge or are not bounded below. Our method of proof requires us to analyze the continuous dependence of finite-interval Lagrange problems with respect to a prescribed terminal condition. Once this is done, we show that a finitely optimal solution can be obtained as the limit of a sequence of solutions to a sequence of corresponding finite-horizon optimal control problems. Our results utilize the convexity and seminormality hypotheses which are now classical in the existence theory of optimal control.
    Type of Medium: Electronic Resource
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 62
    Electronic Resource
    Electronic Resource
    Springer
    Journal of optimization theory and applications 51 (1986), S. 307-320 
    ISSN: 1573-2878
    Keywords: Optimal control ; tracking systems ; computational methods
    Source: Springer Online Journal Archives 1860-2000
    Topics: Mathematics
    Notes: Abstract Feedback control laws are derived for a class of optimal finite time tracking problems with terminal constraints. Analytical solutions are obtained for the feedback gain and the closed-loop response trajectory. Such formulations are expressed in recursive forms so that a real-time computer implementation becomes feasible. An example involving the feedback slewing of a flexible spacecraft is given to illustrate the validity and usefulness of the formulations.
    Type of Medium: Electronic Resource
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 63
    Electronic Resource
    Electronic Resource
    Springer
    Journal of optimization theory and applications 51 (1986), S. 355-362 
    ISSN: 1573-2878
    Keywords: Optimal control ; piecewise periodic feedback gains ; two-point boundary-value problems ; continuous-time Riccati equation ; discrete-time Riccati equation ; continuous-time Lyapunov equation
    Source: Springer Online Journal Archives 1860-2000
    Topics: Mathematics
    Notes: Abstract A noniterative algebraic method is presented for solving differential Riccati equations which satisfy two-point boundary-value problems. This class of numerical problems arises in quadratic optimization problems where the cost functionals are composed of both continuous and discrete state penalties, leading to piecewise periodic feedback gains. The necessary condition defining the solution for the two-point boundary value problem is cast in the form of a discrete-time algebraic Riccati equation, by using a formal representation for the solution of the differential Riccati equation. A numerical example is presented which demonstrates the validity of the approach.
    Type of Medium: Electronic Resource
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 64
    Electronic Resource
    Electronic Resource
    Springer
    Journal of optimization theory and applications 60 (1989), S. 173-190 
    ISSN: 1573-2878
    Keywords: Optimal control ; singular control ; singular control order
    Source: Springer Online Journal Archives 1860-2000
    Topics: Mathematics
    Notes: Abstract Pontryagin's maximum principle gives no information about a singular optimal control if the problem is linear. This survey shows how candidate singular optimal controls may be found for linear and nonlinear problems. A theorem is given on the maximum order of a linear singular problem.
    Type of Medium: Electronic Resource
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 65
    Electronic Resource
    Electronic Resource
    Springer
    Journal of optimization theory and applications 61 (1989), S. 123-136 
    ISSN: 1573-2878
    Keywords: Optimal control ; linear estimation ; state estimation ; least absolute value ; optimal estimation
    Source: Springer Online Journal Archives 1860-2000
    Topics: Mathematics
    Notes: Abstract This paper presents a new technique for solving the problem of linear static state estimation, based on weighted least absolute value (WLAV). A set ofm optimality equations is obtained, wherem=number of measurements, based on minimizing a WLAV performance index involvingn unknown state variables,m〉n. These equations are solved using the left pseudo-inverse transformation, least-square sense, to obtain approximately the residual of each measurement. Ifk is the rank of the matrixH,k=n, we choose among the optimality equations a number of equations equal to the rankk and having the smallest residuals. The solution of thesen equations inn unknowns yields the best WLAV estimation. A numerical example is reported; the results for this example are obtained by using both WLS and WLAV techniques. It is shown that the best WLAV approximation is superior to the best WLS approximation when estimating the true form of data containing some inaccurate observations.
    Type of Medium: Electronic Resource
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 66
    Electronic Resource
    Electronic Resource
    Springer
    Journal of optimization theory and applications 61 (1989), S. 221-245 
    ISSN: 1573-2878
    Keywords: Optimal control ; descriptor systems ; singular systems ; discrete systems
    Source: Springer Online Journal Archives 1860-2000
    Topics: Mathematics
    Notes: Abstract The present paper deals with the investigation of the linear quadratic optimal control problem for the discrete-time descriptor systemsEx k+1=Ax k +Bu k , whereE is in general a singular matrix and the system structure is in general noncausal. The problem is considered in its general form, having singular cost matrices and cross-weighting term in the cost functional. The key idea for the solution approach is the use of the Weierstrass theorem for regular pencils, combined with a suitable permutation transformation, to form a base for the image ofE. The optimization problem is solved by forcing causality to the Hamiltonian equations, which are produced by considering the entireN-stage process as a large system of linear equations. The feedback gain matrix is obtained as a manifold which is generated by the intersection of two other manifolds.
    Type of Medium: Electronic Resource
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 67
    Electronic Resource
    Electronic Resource
    Springer
    Journal of optimization theory and applications 61 (1989), S. 359-376 
    ISSN: 1573-2878
    Keywords: Optimal control ; nonconvex differential inclusions ; endpoint constraints ; maximum principle
    Source: Springer Online Journal Archives 1860-2000
    Topics: Mathematics
    Notes: Abstract We consider an optimization problem with endpoint constraints associated with a nonconvex differential inclusion. We give a necessary condition of the maximum principle type for a solution of the problem. Following the approach from Ref. 1, the condition is stated in terms of single-valued selections of the convexified right-hand side of the inclusion.
    Type of Medium: Electronic Resource
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 68
    Electronic Resource
    Electronic Resource
    Springer
    Journal of optimization theory and applications 62 (1989), S. 289-310 
    ISSN: 1573-2878
    Keywords: Optimal control ; sufficient conditions ; state constraints ; nonconvex optimization
    Source: Springer Online Journal Archives 1860-2000
    Topics: Mathematics
    Notes: Abstract The sufficient optimality conditions of Zeidan for optimal control problems (Refs. 1 and 2) are generalized such that they are applicable to problems with pure state-variable inequality constraints. We derive conditions which neither assume the concavity of the Hamiltonian nor the quasiconcavity of the constraints. Global as well as local optimality conditions are presented.
    Type of Medium: Electronic Resource
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 69
    Electronic Resource
    Electronic Resource
    Springer
    Journal of optimization theory and applications 62 (1989), S. 489-513 
    ISSN: 1573-2878
    Keywords: Optimal control ; regular problems ; state constraints ; method of region analysis
    Source: Springer Online Journal Archives 1860-2000
    Topics: Mathematics
    Notes: Abstract A method of region analysis is developed for solving a class of optimal control problems with one state and one control variable, including state and control constraints. The performance index is strictly convex with respect to the control variable, while this variable appears only linearly in the state equation. The convexity or linearity assumption of the performance index or the state equation with respect to the state variable is not required.
    Type of Medium: Electronic Resource
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 70
    Electronic Resource
    Electronic Resource
    Springer
    Journal of optimization theory and applications 63 (1989), S. 1-22 
    ISSN: 1573-2878
    Keywords: Optimal control ; continuous inequality state constraints ; computational algorithms ; control parametrization
    Source: Springer Online Journal Archives 1860-2000
    Topics: Mathematics
    Notes: Abstract In this paper, we consider a class of optimal control problems involving continuous inequality state constraints. This class of optimal control problems can be solved using the technique developed in a paper by Goh and Teo, where a simple constraint transcription is used to convert continuous inequality state constraints into equivalent equality terminal state constraints. However, that constraint transcription has the disadvantage that the equality terminal state constraints so obtained do not satisfy the usual constraint qualification. Thus, convergence is not guaranteed and some oscillation may exist in numerical computation. The aim of this paper is to use a new constraint transcription together with the concept of control parametrization to devise a new computational algorithm for solving this general class of constrained optimal control problems. This new algorithm is much more stable numerically, as we have successfully overcome the above-mentioned disadvantages.
    Type of Medium: Electronic Resource
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 71
    Electronic Resource
    Electronic Resource
    Springer
    Journal of optimization theory and applications 53 (1987), S. 219-235 
    ISSN: 1573-2878
    Keywords: Optimal control ; existence theory ; sporadically catching-up optimality ; infinite-horizon optimal control problems
    Source: Springer Online Journal Archives 1860-2000
    Topics: Mathematics
    Notes: Abstract In this paper, we extend the existence theory of Brock and Haurie concerning the existence of sporadically catching-up optimal solutions for autonomous, infinite-horizon optimal control problems. This notion of optimality is one of a hierarchy of types of optimality that have appeared in the literature to deal with optimal control problems whose cost functionals, described by an improper integral, either diverge or are unbounded below. Our results rely on the now classical convexity and seminormality hypotheses due to Cesari and are weaker than those assumed in the work of Brock and Haurie. An example is presented where our results are applicable, but those of the above-mentioned authors do not.
    Type of Medium: Electronic Resource
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 72
    Electronic Resource
    Electronic Resource
    Springer
    Journal of optimization theory and applications 54 (1987), S. 403-411 
    ISSN: 1573-2878
    Keywords: Optimal control ; Green's theorem method ; infinite-horizon problems ; most rapid approach
    Source: Springer Online Journal Archives 1860-2000
    Topics: Mathematics
    Notes: Abstract In this paper, we show that in some previous literature the most rapid approach theorems are incomplete. This is proved by providing a counterexample. We also introduce an additional condition that guarantees the optimality of the most rapid approach path.
    Type of Medium: Electronic Resource
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 73
    Electronic Resource
    Electronic Resource
    Springer
    Journal of optimization theory and applications 46 (1985), S. 265-293 
    ISSN: 1573-2878
    Keywords: Optimal control ; terminal equality constraints ; multiplier methods ; augmented Lagrangians ; linearly convergent update formulas
    Source: Springer Online Journal Archives 1860-2000
    Topics: Mathematics
    Notes: Abstract In this paper, a method is proposed for the numerical solution of optimal control problems with terminal equality constraints. The multiplier method is employed to deal with the terminal equality constraints. It is shown that a sequence of control functions, which converges to the optimal control, is obtained by the alternate update of control functions and multipliers.
    Type of Medium: Electronic Resource
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
Close ⊗
This website uses cookies and the analysis tool Matomo. More information can be found here...