作者:
MARAZZI, AUNIV LAUSANNE
INST MED SOCIALE & PREVENT DIV STAT & INFORMAT CH-1005 LAUSANNE SWITZERLAND
A preliminary step for the computation of bounded influence regression is the determination of weights for the points x i in the factor space. This can be achieved through the computation of a robust affinely invarian...
详细信息
A preliminary step for the computation of bounded influence regression is the determination of weights for the points x i in the factor space. This can be achieved through the computation of a robust affinely invariant covariance matrix (with fixed location) and the computation of robust distances of the x i to the origin, in the metric defined by the covariance matrix. For a well known fixed point algorithm convergence has been proved. However, its speed is not very satisfactory in practice. Newton-like and conjugategradient approaches have been suggested by P.J. Huber. This paper reviews the various proposals, discusses their implementations, some of their convergence properties and compares empirically their performances.
We develop a new effective method of determining the rotation R and the stretches U and V in the polar decomposition F RU = VR of the deformation gradient. The method is based on a minimum property of R to have the sm...
详细信息
We develop a new effective method of determining the rotation R and the stretches U and V in the polar decomposition F RU = VR of the deformation gradient. The method is based on a minimum property of R to have the smallest "distance" from F in the Euclidean norm. The proposed method does not require to perform any square root and/or inverse operations. With each F having nine independent components we associate a 4 x 4 symmetric traceless matrix Q. The rotation is described by quaternion parameters from which a quadrivector X is formed. It is shown that X corresponding to R maximizes X(T)QX over all X satisfying (XX)-X-T = 1. We prefer to equivalently maximize the Rayleigh quotient Y(T)QY/(YY)-Y-T over all non-vanishing Y and to deduce X by subsequent normalization X = Y/root(YY)-Y-T. The maximization of the Rayleigh quotient is performed by a conjugate gradient algorithm with all iterative steps carried out by explicit closed formulae. Efficiency and accuracy of the method is illustrated by a numerical example.
A time-domain hyperbolic Radon transform based method for separating multicomponent seismic data into P-P and P-SV wavefields is presented. This wavefield separation method isolates P-P and P-SV wavefields in the Rado...
详细信息
A time-domain hyperbolic Radon transform based method for separating multicomponent seismic data into P-P and P-SV wavefields is presented. This wavefield separation method isolates P-P and P-SV wavefields in the Radon panel due to their differences in slowness, and an inverse transform of only part of the data leads to separated wavefields. A problem of hyperbolic Radon transform is that it works in the time domain entailing the inversion of large operators which is prohibitively time-consuming. By applying the conjugate gradient algorithm during the inversion of hyperbolic Radon transform, the computational cost can be kept reasonably low for practical application. Synthetic data examples prove that P-P and P-SV wavefield separation by hyperbolic Radon transform produces more accurate separated wavefields compared with separation by high-resolution parabolic Radon transform, and the feasibility of the proposed separation scheme is also verified by a real field data example.
The existing structural rating system in the maintenance steering group-3 is facing difficulties with developing new aircraft. There is always a lack of available data to help make proper inspection intervals. By taki...
详细信息
The existing structural rating system in the maintenance steering group-3 is facing difficulties with developing new aircraft. There is always a lack of available data to help make proper inspection intervals. By taking advantage of the powerful learning ability and data fusion capability of the backpropagation network, this paper aims to develop a new decision support system for the determination of scheduled maintenance inspection intervals for aircraft structures susceptible to accidental damage. As long as various data are available, regardless of engineering experience, experimental values, and in-service reliability data feedback, the backpropagation network has been proven able to make acceptable predictions after training, making accidental damage rating a dynamic and flexible procedure.
The article addresses the problem of efficient analysis of the statistics of initial and corrected shape distortions in space structures. Two approaches for improving efficiency are considered. One is an adjoint techn...
详细信息
Today the adjustment of structural models is an essential step in the modeling of complex structures. In this paper, we are interested in the improvement of finite element models. Our approach is a parametric updating...
详细信息
Today the adjustment of structural models is an essential step in the modeling of complex structures. In this paper, we are interested in the improvement of finite element models. Our approach is a parametric updating using modal test results, which supply eigenvalues and associated eigenvectors. It is based on the computation of the error measure on the constitutive relation and allows us to correct both the stiffness and the mass matrices. In particular, this paper shows how this tuning strategy can improve a given finite element model when the measures are noisy. Several simulation examples illustrate the behavior of this method.
A new method of coherent summation of laser beams is proposed. Two-dimension numerical simulation based on the conjugate gradient algorithm shows that the combination efficiency is about 94.8%, and an error less than ...
详细信息
A new method of coherent summation of laser beams is proposed. Two-dimension numerical simulation based on the conjugate gradient algorithm shows that the combination efficiency is about 94.8%, and an error less than 2.4% with a 16-level phase quantized diffraction optical element. Error analysis has also been performed. (C) 2003 Elsevier Science Ltd. All rights reserved.
Tsunamis are gravitational, i.e. gravity-controlled waves generated by a given motion of the bottom. There are different natural phenomena, such as submarine slumps, slides, volcanic explosions, earthquakes, etc. that...
详细信息
Tsunamis are gravitational, i.e. gravity-controlled waves generated by a given motion of the bottom. There are different natural phenomena, such as submarine slumps, slides, volcanic explosions, earthquakes, etc. that can lead to a tsunami. This paper deals with the case where the tsunami source is an earthquake. The mathematical model studied here is based on shallow water theory, which is used extensively in tsunami modeling. The inverse problem consists of determining an unknown initial tsunami source q(x, y) by using measurements f(m) (t) of the height of a passing tsunami wave at a finite number of given points (x(m), y(m)), m = 1,2,..., M, of the coastal area. The proposed approach is based on the weak solution theory for hyperbolic PDEs and adjoint problem method for minimization of the corresponding cost functional. The adjoint problem is defined to obtain an explicit gradient formula for the cost functional J(q) = parallel to Aq - F parallel to(2), F = (f(1),..., f(M)). Numerical algorithms are proposed for the direct as well as adjoint problems. conjugate gradient algorithm based on explicit gradient formula is used for numerical solution of the inverse problem. Results of computational experiments presented for the synthetic noise free and random noisy data in real scale illustrate bounds of applicability of the proposed approach, also its efficiency and accuracy. (C) 2014 IMACS. Published by Elsevier B.V. All rights reserved.
Through the use of equivalent modeling and simple hand calculations, an optimization of a truss structure was achieved. In addition, valuable insight was gained into the relative performance of a tubular and solid mem...
详细信息
Through the use of equivalent modeling and simple hand calculations, an optimization of a truss structure was achieved. In addition, valuable insight was gained into the relative performance of a tubular and solid member truss structure. This approach could be extended to more complicated structures, resulting in inexpensive preliminary optimization. However, the analyst must be careful in selecting his equivalent models in order to avoid erroneous results.
作者:
Liu, YangUniv Maryland
Dept Human Dev & Quantitat Methodol 12308 Benjamin Bldg3942 Campus Dr College Pk MD 20742 USA
In exploratory factor analysis, latent factors and factor loadings are seldom interpretable until analytic rotation is performed. Typically, the rotation problem is solved by numerically searching for an element in th...
详细信息
In exploratory factor analysis, latent factors and factor loadings are seldom interpretable until analytic rotation is performed. Typically, the rotation problem is solved by numerically searching for an element in the manifold of orthogonal or oblique rotation matrices such that the rotated factor loadings minimize a pre-specified complexity function. The widely used gradient projection (GP) algorithm, although simple to program and able to deal with both orthogonal and oblique rotation, is found to suffer from slow convergence when the number of manifest variables and/or the number of latent factors is large. The present work examines the effectiveness of two Riemannian second-order algorithms, which respectively generalize the well-established truncated Newton and trust-region strategies for unconstrained optimization in Euclidean spaces, in solving the rotation problem. When approaching a local minimum, the second-order algorithms usually converge superlinearly or even quadratically, better than first-order algorithms that only converge linearly. It is further observed in Monte Carlo studies that, compared to the GP algorithm, the Riemannian truncated Newton and trust-region algorithms require not only much fewer iterations but also much less processing time to meet the same convergence criterion, especially in the case of oblique rotation.
暂无评论