• Title/Summary/Keyword: inverse matrices

Search Result 119, Processing Time 0.026 seconds

Numerical Verification of Hybrid Optimization Technique for Finite Element Model Updating (유한요소모델개선을 위한 하이브리드 최적화기법의 수치해석 검증)

  • Jung, Dae-Sung;Kim, Chul-Young
    • Journal of the Earthquake Engineering Society of Korea
    • /
    • v.10 no.6 s.52
    • /
    • pp.19-28
    • /
    • 2006
  • Most conventional model updating methods must use mathematical objective function with experimental modal matrices and analytical system matrices or must use information about the gradient or higher derivatives of modal properties with respect to each updating parameter. Therefore, most conventional methods are not appropriate for complex structural system such as bridge structures due to stability problem in inverse analysis with ill-conditions. Sometimes, moreover, the updated model may have no physical meaning. In this paper, a new FE model updating method based on a hybrid optimization technique using genetic algorithm (GA) and Holder-Mead simplex method (NMS) is proposed. The performance of hybrid optimization technique on the nonlinear problem is demonstrated by the Goldstein-Price function with three local minima and one global minimum. The influence of the objective function is evaluated by the case study of a simulated 10-dof spring-mass model. Through simulated case studies, finally, the objective function is proposed to update mass as well as stiffness at the same time. And so, the proposed hybrid optimization technique is proved to be an efficient method for FE model updating.

Study on Class Separability Measure for Radar Signals (레이다 신호의 클래스 분리도 측정을 위한 연구)

  • Jeong, Seong-Jae;Lee, Seung-Jae;Kim, Kyung-Tae
    • The Journal of Korean Institute of Electromagnetic Engineering and Science
    • /
    • v.29 no.2
    • /
    • pp.128-137
    • /
    • 2018
  • In this paper, we propose a novel class separability measure for radar signals. To reduce the sensitivity of the relative aspect angle between a target and radar, to evaluate the discriminatory power of radar signals, the proposed method first calculates the correlation coefficients between two radar cross sections (RCSs) or linearly shifts one-dimensional (1D) radar signals (i.e., high-resolution range profiles (HRRPs)), or rotates two 2D radar signals (i.e., inverse synthetic aperture radar (ISAR) images). Then, it uses the maximum correlation coefficient when two radar signals are best aligned. Next, the proposed method obtains new correlation-based discriminant matrices (CDM) using maximum correlation coefficients. Finally, the cumulative distribution function (CDF) in the CDM and the value corresponding to the specific probability in the CDF are obtained, and this value represents the discriminatory power of the radar signal. Experimental results show that the proposed method can accurately measure the target separability.

A PRECONDITIONER FOR THE LSQR ALGORITHM

  • Karimi, Saeed;Salkuyeh, Davod Khojasteh;Toutounian, Faezeh
    • Journal of applied mathematics & informatics
    • /
    • v.26 no.1_2
    • /
    • pp.213-222
    • /
    • 2008
  • Iterative methods are often suitable for solving least squares problems min$||Ax-b||_2$, where A $\epsilon\;\mathbb{R}^{m{\times}n}$ is large and sparse. The well known LSQR algorithm is among the iterative methods for solving these problems. A good preconditioner is often needed to speedup the LSQR convergence. In this paper we present the numerical experiments of applying a well known preconditioner for the LSQR algorithm. The preconditioner is based on the $A^T$ A-orthogonalization process which furnishes an incomplete upper-lower factorization of the inverse of the normal matrix $A^T$ A. The main advantage of this preconditioner is that we apply only one of the factors as a right preconditioner for the LSQR algorithm applied to the least squares problem min$||Ax-b||_2$. The preconditioner needs only the sparse matrix-vector product operations and significantly reduces the solution time compared to the unpreconditioned iteration. Finally, some numerical experiments on test matrices from Harwell-Boeing collection are presented to show the robustness and efficiency of this preconditioner.

  • PDF

A Robust Preconditioner on the CRAY-T3E for Large Nonsymmetric Sparse Linear Systems

  • Ma, Sangback;Cho, Jaeyoung
    • Journal of the Korean Society for Industrial and Applied Mathematics
    • /
    • v.5 no.1
    • /
    • pp.85-100
    • /
    • 2001
  • In this paper we propose a block-type parallel preconditioner for solving large sparse nonsymmetric linear systems, which we expect to be scalable. It is Multi-Color Block SOR preconditioner, combined with direct sparse matrix solver. For the Laplacian matrix the SOR method is known to have a nondeteriorating rate of convergence when used with Multi-Color ordering. Since most of the time is spent on the diagonal inversion, which is done on each processor, we expect it to be a good scalable preconditioner. Finally, due to the blocking effect, it will be effective for ill-conditioned problems. We compared it with four other preconditioners, which are ILU(0)-wavefront ordering, ILU(0)-Multi-Color ordering, SPAI(SParse Approximate Inverse), and SSOR preconditioner. Experiments were conducted for the Finite Difference discretizations of two problems with various meshsizes varying up to 1024 x 1024, and for an ill-conditioned matrix from the shell problem from the Harwell-Boeing collection. CRAY-T3E with 128 nodes was used. MPI library was used for interprocess communications. The results show that Multi-Color Block SOR and ILU(0) with Multi-Color ordering give the best performances for the finite difference matrices and for the shell problem only the Multi-Color Block SOR converges.

  • PDF

ON POSITIVE DEFINITE SOLUTIONS OF A CLASS OF NONLINEAR MATRIX EQUATION

  • Fang, Liang;Liu, San-Yang;Yin, Xiao-Yan
    • Bulletin of the Korean Mathematical Society
    • /
    • v.55 no.2
    • /
    • pp.431-448
    • /
    • 2018
  • This paper is concerned with the positive definite solutions of the nonlinear matrix equation $X-A^*{\bar{X}}^{-1}A=Q$, where A, Q are given complex matrices with Q positive definite. We show that such a matrix equation always has a unique positive definite solution and if A is nonsingular, it also has a unique negative definite solution. Moreover, based on Sherman-Morrison-Woodbury formula, we derive elegant relationships between solutions of $X-A^*{\bar{X}}^{-1}A=I$ and the well-studied standard nonlinear matrix equation $Y+B^*Y^{-1}B=Q$, where B, Q are uniquely determined by A. Then several effective numerical algorithms for the unique positive definite solution of $X-A^*{\bar{X}}^{-1}A=Q$ with linear or quadratic convergence rate such as inverse-free fixed-point iteration, structure-preserving doubling algorithm, Newton algorithm are proposed. Numerical examples are presented to illustrate the effectiveness of all the theoretical results and the behavior of the considered algorithms.

Energy efficient joint iterative SIC-MMSE MIMO detection (에너지 효율적 반복 SIC-MMSE MIMO 검출)

  • Ngayahala, F.C. Kamaha;Ahmed, Saleem;Kim, Sooyoung
    • Journal of Satellite, Information and Communications
    • /
    • v.10 no.1
    • /
    • pp.22-28
    • /
    • 2015
  • In this paper, we propose a new computationally efficient joint iterative multi-input multi-output (MIMO) detection scheme using a soft interference cancellation and minimum mean squared-error (SIC-MMSE) method. The critical computational burden of the SIC-MMSE scheme lies in the multiple inverse operations of the complex matrices. We find a new way which requires only a single matrix inversion by utilizing the Taylor series expansion of the matrix, and thus the computational complexity can be reduced. The computational complexity reduction increases as the number of antennas is increased. The simulation results show that our method produces almost the same performances as the conventional SIC-MMSE with reduced computational complexity.

Finite Element Analysis of Shape Rolling Process using Destributive Parallel Algorithms on Cray T3E (병렬 컴퓨터를 이용한 형상 압연공정 유한요소 해석의 분산병렬처리에 관한 연구)

  • Gwon, Gi-Chan;Yun, Seong-Gi
    • Transactions of the Korean Society of Mechanical Engineers A
    • /
    • v.24 no.5 s.176
    • /
    • pp.1215-1230
    • /
    • 2000
  • Parallel Approaches using Cray T3E which is NIPP (Massively Parallel Processors) machine are presented for the efficient computation of the finite element analysis of 3-D shape rolling processes. D omain decomposition method coupled with parallel linear equation solver is used. Domain decomposition is applied for obtaining element tangent stifffiess matrices and residual vectors. Direct and iterative parallel algorithms are used for solving the linear equations. Direct algorithm is_parallel version of direct banded matrix solver. For iterative algorithms, the well-known preconditioned conjugate gradient solver with Jacobi preconditioner is also employed. Moreover a new effective iterative scheme with block inverse matrix preconditioner, which is named by present authors, is presented and its results are compared with the one using Jacobi preconditioner. PVM and MPI are used for message passing and synchronization between processors. The performance and efficiency of each algorithm is discussed and comparisons are made among different algorithms.

Log-density Ratio with Two Predictors in a Logistic Regression Model (로지스틱 회귀모형에서 이변량 정규분포에 근거한 로그-밀도비)

  • Kahng, Myung Wook;Yoon, Jae Eun
    • The Korean Journal of Applied Statistics
    • /
    • v.26 no.1
    • /
    • pp.141-149
    • /
    • 2013
  • We present methods for studying the log-density ratio that enables the selection of the predictors and the form to be included in the logistic regression model. Under bivariate normal distributional assumptions, we investigate the form of the log-density ratio as a function of two predictors. If two covariance matrices are equal, then the crossproduct and quadratic terms are not needed. If the variables are uncorrelated, we do not need the crossproduct terms, but we still need the linear and quadratic terms. We also explore other conditions in which the crossproduct and quadratic terms are not needed in the logistic regression model.

Structural Damage Detection through System Identification (시스템 동정을 통한 구조물의 결함 탐지)

  • Koh, Bong-Hwan;Nagarajaiah, S.;Phan, M.Q.
    • Proceedings of the Korean Society for Noise and Vibration Engineering Conference
    • /
    • 2006.05a
    • /
    • pp.1223-1228
    • /
    • 2006
  • This paper presents an experimental investigation of a recently developed Kronecker Product (KP) method to determine the type, location, and intensity of structural damage from an identified state-space model of the system. Although this inverse problem appears to be highly nonlinear, the system mass, stiffness, and damping matrices are identified through a series of transformations, and with the aid of the Kronecker product, only linear operations are involved in the process. Since a state-space model can be identified directly from input-output data, an initial finite element model and/or model updating are not required. The test structure is a two-degree-of-freedom torsional system in which mass and stiffness are arbitrarily adjustable to simulate various conditions of structural damage. This simple apparatus demonstrates the capability of the damage detection method by not only identifying the location and the extent of the damage, but also differentiating the nature of the damage. The potential applicability of the KP method for structural damage identification is confirmed by laboratory test.

  • PDF

Spectral encapsulation of OFDM systems based on orthogonalization for short packet transmission

  • Kim, Myungsup;Kwak, Do Young;Kim, Ki-Man;Kim, Wan-Jin
    • ETRI Journal
    • /
    • v.42 no.6
    • /
    • pp.859-871
    • /
    • 2020
  • A spectrally encapsulated (SE) orthogonal frequency-division multiplexing (OFDM) precoding scheme for wireless short packet transmission, which can suppress the out-of-band emission (OoBE) while maintaining the advantage of the cyclic prefix (CP)-OFDM, is proposed. The SE-OFDM symbol consists of a prefix, an inverse fast Fourier transform (IFFT) symbol, and a suffix generated by the head, center, and tail matrices, respectively. The prefix and suffix play the roles of a guard interval and suppress the OoBE, and the IFFT symbol has the same size as the discrete Fourier transform symbol in the CP-OFDM symbol and serves as an information field. Specifically, as the center matrix generating the IFFT symbol is orthogonal, data and pilot symbols can be allocated to any subcarrier without distinction. Even if the proposed precoder is required to generate OFDM symbols with spectral efficiency in the transmitter, a corresponding decoder is not required in the receiver. The proposed scheme is compared with CP-OFDM in terms of spectrum, OoBE, and bit-error rate.