scholarly journals On the Evaluation of Determinants, the Formation of their Adjugates, and the Practical Solution of Simultaneous Linear Equations

1933 ◽  
Vol 3 (3) ◽  
pp. 207-219 ◽  
Author(s):  
A. C. Aitken

There are various methods in existence for the practical solution of a set of simultaneous equationsSome of these methods are appropriate to special systems, as for example to the axisymmetric “normal equations” of Least Squares. In many applications, however, as in problems of statistical correlation of many variables, it may be desired not merely to solve a given set of equations but to obtain as much knowledge as possible about the system or matrix of coefficients; perhaps to evaluate its determinant and various minors, such as principal minors, possibly also to determine the elements of the adjugate matrix, or the reciprocal matrix. The examination of the sign of successive principal minors of an axisymmetric determinant, in order to find the signature of the corresponding quadratic form, is a case in point; and there are many such applications.

1935 ◽  
Vol 54 ◽  
pp. 12-16 ◽  
Author(s):  
A. C. Aitken

This paper concludes the study of fitting polynomials by Least Squares, treated in two previous papers. The problem being concerned with the minimum of a positive definite quadratic form, it makes for conciseness to use matrix notation. We shall therefore adopt the following conventions :—The n values of the variable x, of the data u0, u1, …, un−1, of certain polynomials qr(x) entering into the solution, and so on, will be regarded compositely as vectors. They will be imagined as having their components or elements disposed in column array, but when written in full will be written horizontally, to save space, enclosed by curled brackets. Row vectors, when written out in full, will be enclosed by square brackets. In the shorter notation we shall write, for example, u, x for column vectors, u′, x′ for the row vectors obtained by transposition. The vectors occurring in the problem will be the following:—


1945 ◽  
Vol 49 (411) ◽  
pp. 104-111
Author(s):  
F. J. Turton

The application of strain energy or slope-deflection methods in the analysis of redundant structures leads to a number of simultaneous linear equations with numerical coefficients; the equations may be obtained in such order that each successive equation contains one new unknown, until all the unknowns are so included. This is the only condition essential for the method to be described in the present paper, but the labour is much reduced in slope-deflection and strain energy applications by the fact that most (or all) of the equations contain very few of the unknowns. The method to be given reduces the solving of these equations to a column of successive evaluations, followed by the solution, by algebraic methods, of a small number of simultaneous equations; and a final column of evaluations. In the remaining paragraphs a number of problems are examined to show how the equations may be obtained in suitable sequence for the method to apply. Following an application to the determination of secondary stresses, the operations involved in the moment-distribution method and in this method are compared. A numerical example is worked out in the simple case of §2, and it is shown how any order of mathematical accuracy in the roots may be ensured, provided that sufficient figures have been retained to permit that accuracy.


A physical interpretation is made of various complicated formulae which have been given for the thermal conductivity A and viscosity n of mixtures of gases. The interpretation is based on the recognition of two principal effects operating in the transport of heat or momentum through gaseous mixtures. The first (and larger) effect is that molecules of one species impede transport of heat or momentum by other species. The second effect is a transfer of the transport of heat (or momentum) from one species to another. When the transfer of transport is neglected, equations of the form proposed by Sutherland (1895) and Wassiljewa (1904) follow immediately. For the thermal conductivity (symbols used are defined in the main text): A = E{ n i /( n i A i -1 +Ea ij n j )] There is an identical expression for the viscosity, though the values of a ij are different from those for the thermal conductivity. The significance to be given to each term of the sum over i is that of a quotient of (i) a force driving conduction, proportional to n i , and (ii) a resistance due to species i(n i A i -1 ) and to all other species (En j a ij ). Here A i -1 is the resistance offered by species i to its own transport, a ij and the resistance offered by the species j to transport by i . When the transfer of transport is taken into account, two simultaneous linear equations have to be solved for a binary mixture, and the solution has the form of the quotient of quadratics familiar in the theoretical analysis of mixtures of monatomic gases. For a mixture of N constituents, the resulting expression for A appears as a quotient of determinants. Application of the same principles also throws light on the transport of internal energy in mixtures of polyatomic gases. The total thermal conductivity may be divided into contributions from each species, and each such contribution may be further subdivided into internal and translational parts: thus for a binary mixture it is necessary to replace the pair of simultaneous equations by four. Such linear equations represent a direct generalization of the equations of Mason & Monchick (1962) for a simple gas. In an appendix, the physically significant parameters of the generalized approach employed here are compared explicitly with the predictions of rigorous analysis for mixtures of monatomic gases.


2009 ◽  
Vol 24 (S1) ◽  
pp. S16-S21 ◽  
Author(s):  
Balder Ortner

The sin2ψ method can be formulated as a single system of simultaneous linear equations. Using this it is easy to show that the sin2ψ method is not a least-squares method. It further helps to compare the accuracies of the stress tensors obtained by the sin2ψ method and the method of least squares. Quantitative comparisons have been made for different fictitious measurements. It is shown that the unnecessary loss in accuracy by using the sin2ψ method is quite significant and by no means negligible. The same course of action has been applied to compare the so-called Dölle-Hauk method with a least-squares method; the result is similar. Some other methods for X-ray stress determination, most often similar to the sin2ψ method, and their shortcomings are also discussed briefly, together with the corresponding, more effective and more versatile least-squares method.


1976 ◽  
Vol 157 (2) ◽  
pp. 489-492 ◽  
Author(s):  
I A Nimmo ◽  
G L Atkins

1. Descriptions are given of two ways for fitting non-linear equations by least-squares criteria to experimental data. One depends on solving a set of non-linear simultaneous equations, and the other on Taylor's theorem. 2. It is shown that better parameter estimates result when an equation with two or more non-linear parameters is fitted to all the sets of data simultaneously than when it is fitted to each set in turn.


Author(s):  
J. M. Taylor

SynopsisIt has been known for some time that certain least-squares problems are “ill-conditioned”, and that it is therefore difficult to compute an accurate solution. The degree of ill-conditioning depends on the basis chosen for the subspace in which it is desired to find an approximation. This paper characterizes the degree of ill-conditioning, for a general inner-product space, in terms of the basis.The results are applied to least-squares polynomial approximation. It is shown, for example, that the powers {1, z, z2,…} are a universally bad choice of basis. In this case, the condition numbers of the associated matrices of the normal equations grow at least as fast as 4n, where n is the degree of the approximating polynomial.Analogous results are given for the problem of finite interpolation, which is closely related to the least-squares problem.Applications of the results are given to two algorithms—the Method of Moments for solving linear equations and Krylov's Method for computing the characteristic polynomial of a matrix.


2019 ◽  
Vol 149 (5) ◽  
pp. 1241-1290
Author(s):  
Pierre-Yves Bienvenu

AbstractWe prove asymptotic formulae for sums of the form $$\sum\limits_{n\in {\open z}^d\cap K} {\prod\limits_{i = 1}^t {F_i} } (\psi _i(n)),$$where K is a convex body, each Fi is either the von Mangoldt function or the representation function of a quadratic form, and Ψ = (ψ1, …, ψt) is a system of linear forms of finite complexity. When all the functions Fi are equal to the von Mangoldt function, we recover a result of Green and Tao, while when they are all representation functions of quadratic forms, we recover a result of Matthiesen. Our formulae imply asymptotics for some polynomial patterns in the primes. For instance, they describe the asymptotic behaviour of the number of k-term arithmetic progressions of primes whose common difference is a sum of two squares.The paper combines ingredients from the work of Green and Tao on linear equations in primes and that of Matthiesen on linear correlations amongst integers represented by a quadratic form. To make the von Mangoldt function compatible with the representation function of a quadratic form, we provide a new pseudorandom majorant for both – an average of the known majorants for each of the functions – and prove that it has the required pseudorandomness properties.


1935 ◽  
Vol 39 (292) ◽  
pp. 349-351 ◽  
Author(s):  
J. Morris

Simultaneous equations of three or more variables are notoriously trouble some to solve numerically and, moreover, are frequently sensitive for certain relative values of the variables. This usually precludes the use of the slide rule and in consequence resort has to be had to accurate working to an increasing number of decimal places involving much labour even with a calculating machine.A method is here given which is applicable to certain classes of equations of frequent occurrence in mathematical physics and engineering science, which method is a combination of what is knovwi as the Iteration process and the means of applying this principle to end moment distribution devised by Professor Hardy Cross.


Author(s):  
Andrii Sohor ◽  
◽  
Markiian Sohor ◽  

The most reliable method for calculating linear equations of the least squares principle, which can be used to solve incorrect geodetic problems, is based on matrix factorization, which is called a singular expansion. There are other methods that require less machine time and memory. But they are less effective in taking into account the errors of the source information, rounding errors and linear dependence. The methodology of such research is that for any matrix A and any two orthogonal matrices U and V there is a matrix Σ, which is determined from the ratio. The idea of a singular decomposition is that by choosing the right matrices U and V, you can convert most elements of the matrix to zero and make it diagonal with non-negative elements. The novelty and relevance of scientific solutions lies in the feasibility of using a singular decomposition of the matrix to obtain linear equations of the least squares method, which can be used to solve incorrect geodetic problems. The purpose of scientific research is to obtain a stable solution of parametric equations of corrections to the results of measurements in incorrect geodetic problems. Based on the performed research on the application of the singular decomposition method in solving incorrect geodetic problems, we can summarize the following. A singular expansion of a real matrix is any factorization of a matrix with orthogonal columns , an orthogonal matrix and a diagonal matrix , the elements of which are called singular numbers of the matrix , and the columns of matrices and left and right singular vectors. If the matrix has a full rank, then its solution will be unique and stable, which can be obtained by different methods. But the method of singular decomposition, in contrast to other methods, makes it possible to solve problems with incomplete rank. Research shows that the method of solving normal equations by sequential exclusion of unknowns (Gaussian method), which is quite common in geodesy, does not provide stable solutions for poorly conditioned or incorrect geodetic problems. Therefore, in the case of unstable systems of equations, it is proposed to use the method of singular matrix decomposition, which in computational mathematics is called SVD. The SVD singular decomposition method makes it possible to obtain stable solutions of both stable and unstable problems by nature. This possibility to solve incorrect geodetic problems is associated with the application of some limit τ, the choice of which can be made by the relative errors of the matrix of coefficients of parametric equations of corrections and the vector of results of geodetic measurements . Moreover, the solution of the system of normal equations obtained by the SVD method will have the shortest length. Thus, applying the apparatus of the singular decomposition of the matrix of coefficients of parametric equations of corrections to the results of geodetic measurements, we obtained new formulas for estimating the accuracy of the least squares method in solving incorrect geodetic problems. The derived formulas have a compact form and make it possible to easily calculate the elements and estimates of accuracy, almost ignoring the complex procedure of rotation of the matrix of coefficients of normal equations.


Sign in / Sign up

Export Citation Format

Share Document