• KSII Transactions on Internet and Information Systems (TIIS)
• /
• 제12권4호
• /
• pp.1504-1526
• /
• 2018

For uplink multi-user massive MIMO systems, conventional minimum mean square error (MMSE) linear detection method achieves near-optimal performance when the number of antennas at base station is much larger than that of the single-antenna users. However, MMSE detection involves complicated matrix inversion, thus making it cumbersome to be implemented cost-effectively and rapidly. In this paper, we first summarize in detail the state-of-the-art simplified MMSE detection algorithms that circumvent the complicated matrix inversion and hence reduce the computation complexity from ${\mathcal{O}}(K^3)$ to ${\mathcal{O}}(K^2)$ or ${\mathcal{O}}(NK)$ with some certain performance sacrifice. Meanwhile, we divide the simplified algorithms into two categories, namely the matrix inversion approximation and the classical iterative linear equation solving methods, and make comparisons between them in terms of detection performance and computation complexity. In order to further optimize the detection performance of the existing detection algorithms, we propose more proper solutions to set the initial values and relaxation parameters, and present a new way of reconstructing the exact effective noise variance to accelerate the convergence speed. Analysis and simulation results verify that with the help of proper initial values and parameters, the simplified matrix inversion based detection algorithms can achieve detection performance quite close to that of the ideal matrix inversion based MMSE algorithm with only a small number of series expansions or iterations.

• 지구물리와물리탐사
• /
• 제8권2호
• /
• pp.129-136
• /
• 2005

감쇠최소자승법은 각종 물리탐사 자료에 가장 널리 사용되는 역산법이다. 일반적으로 최소자승법에서 최소화되는 목적함수는 자료오차(data misfit)와 모델제한자의 합으로 주어진다. 따라서 역산에서 자료오차와 모델제한자는 함께 중요한 역할을 담당한다. 하지만 역산에 관한 대부분의 연구는 주로 모델제한자의 설정방법과 적절한 라그랑지 곱수의 선정방법에 치중되어 왔다. 일반적으로 자료획득시 자료가 갖는 표준편차를 자료가중값의 계산에 사용하는 것이 추천되고 있지만, 실제 현장조사에서는 자료의 표준편차는 좀처럼 측정되지 않으며, 대부분의 역산에서 자료가중행렬은 어쩔 수 없이 단위행렬로 간주된다. 본 논문에서는 자료분해능행렬과 그 분산함수를 분석하여 자동적으로 계산된 자료가중행렬을 사용하는 역산법을 개발하였다. EACB법이라 명명한 이 역산법에서는 분해능이 높은 자료에는 높은 가중값을, 작은 자료에는 작은 가중값을 부여한다. 개발된 EACB 역산법을 전기비저항 토모그피법에 적용한 결과, 보다 안정적이고 분해능이 향상된 결과를 얻을 수 있었다.

• 한국지구물리탐사학회:학술대회논문집
• /
• 한국지구물리탐사학회 2007년도 공동학술대회 논문집
• /
• pp.329-336
• /
• 2007

For scaling of the gradient of misfit function, we develop a new pseudo-Hessian matrix constructed by combining amplitude field and pseudo-Hessian matrix. Since pseudo- Hessian matrix neglects the calculation of the zero-lag auto-correlation of impulse responses in the approximate Hessian matrix, the pseudo-Hessian matrix has a limitation to scale the gradient of misfit function compared to the approximate Hessian matrix. To validate the new pseudo- Hessian matrix, we perform frequency-domain elastic full waveform inversion using this Hessian matrix. By synthetic experiments, we show that the new pseudo-Hessian matrix can give better convergence to the true model than the old one does. Furthermore, since the amplitude fields are intrinsically obtained in forward modeling procedure, we do not have to pay any extra cost to compute the new pseudo-Hessian. We think that the new pseudo-Hessian matrix can be used as an alternative of the approximate Hessian matrix of the Gauss-Newton method.

• 지구물리와물리탐사
• /
• 제17권2호
• /
• pp.95-103
• /
• 2014

MT 자료에서 3차원 전기비저항 구조와 정적효과를 동시에 구하기 위한 실용적인 역산법을 소개한다. 이 방법은 감도행렬이 필요한 Gauss-Newton법을 기반으로 하고 반복과정에서 Broyden의 방식으로 감도를 수정하는 것을 기본으로 하고 있다. 이 논문에서는 합성 MT 자료에 대한 역산실험을 통해 근사역산법의 성능과 정적효과에 대한 가중치에 대해 검토하였다. 해석적으로 구해지는 초기감도를 Broyden의 방식으로 수정하는 역산법은 초기감도를 끝까지 쓰는 역산법보다 자료오차를 줄이는데 효과적이었다. 그리고 완전한 감도행렬을 반복 중간에서 단 한번만 사용하는 근사역산법으로서는 반복 전반부에서 완전한 감도를 사용할 때 자료오차를 가장 많이 줄이는 것으로 나타났다. 정적효과에 대한 가중치는 어느 특정 한계값 이하로 선택하면 최종 자료오차에는 결정적인 영향을 주지 않는다. 합성 MT 자료에 대한 실험 결과 이 역산법은 정적효과가 포함된 MT 자료로부터 3차원 전기비저항 구조를 재현하는데 효과적임을 확인하였다.

• 한국지진공학회:학술대회논문집
• /
• 한국지진공학회 2002년도 춘계 학술발표회 논문집
• /
• pp.20-27
• /
• 2002

Conventional Levenberg-Marquardt's nonlinear inversion method is simply modified by taking into account the second derivatives of the Hessian matrix so as to give robust inversion results. The weight of the second derivative terms is determined by the value of so-called λ in Levenberg-Marquardt's method. The new inversion method is applied to observed data from small-to-moderate earthquakes to simultaneously evaluate the modes parameters of the stochastic point-source model in and around the Korean Peninsula. Best estimates of the stochastic model parameters are obtained along with their statistics and compared with the previous results. Overall characteristics of the model parameters are found to be more of those of interplate than intraplate tectonic region.

• 한국소음진동공학회:학술대회논문집
• /
• 한국소음진동공학회 2000년도 추계학술대회논문집
• /
• pp.201-206
• /
• 2000

In this work transfer function synthesis method based on FRF data of each substructure is investigated for a complex structure composed of many substructures. Though the transfer function synthesis method has superiority to analyze the characteristics of interfaces among substructures effectively, many problems arise in the computation process, especially matrix inversion process. Due to computational problems, the error between the data obtained by test and the predictions through computations is inevitable. So in this paper, computational aspects in the transfer function synthesis method are examined through a steering system problem of passenger car. For the FBS method, frequency response functions of 3 substructures are measured experimentally. Effects of several parameters such as matrix inversion method, connection conditions between substructures and off-diagonal terms on system response are studied numerically.

• 한국지구과학회지
• /
• 제37권5호
• /
• pp.277-285
• /
• 2016

We introduce a depth scaling strategy to improve the accuracy of frequency-domain elastic full waveform inversion (FWI) using the new pseudo-Hessian matrix for seismic data without low-frequency components. The depth scaling strategy is based on the fact that the damping factor in the Levenberg-Marquardt method controls the energy concentration in the gradient. In other words, a large damping factor makes the Levenberg-Marquardt method similar to the steepest-descent method, by which shallow structures are mainly recovered. With a small damping factor, the Levenberg-Marquardt method becomes similar to the Gauss-Newton methods by which we can resolve deep structures as well as shallow structures. In our depth scaling strategy, a large damping factor is used in the early stage and then decreases automatically with the trend of error as the iteration goes on. With the depth scaling strategy, we can gradually move the parameter-searching region from shallow to deep parts. This flexible damping factor plays a role in retarding the model parameter update for shallow parts and mainly inverting deeper parts in the later stage of inversion. By doing so, we can improve deep parts in inversion results. The depth scaling strategy is applied to synthetic data without lowfrequency components for a modified version of the SEG/EAGE overthrust model. Numerical examples show that the flexible damping factor yields better results than the constant damping factor when reliable low-frequency components are missing.

• 전기의세계
• /
• 제23권3호
• /
• pp.43-52
• /
• 1974

The matrix inversion is very inefficient for computing direct solutions of the large spare systems of linear equations that arise in many network problems as a large electrical power system. Optimally ordered triangular factorization of sparse matrices is more efficient and offers the other important computational advantages in some applications with this method. The direct solutions are computed from sparse matrix factors instead of a full inverse matrix, thereby gaining a significant advantage is speed and computer memory requirements. In this paper, it is shown that the sparse matrix method is superior to the inverse matrix method to solve the linear equations of large sparse networks. In addition, it is shown that the sparse matrix method is superior to the inverse matrix method to solve the linear equations of large sparse networks. In addition, it is shown that the solutions may be applied directly to sove the load flow in an electrical power system. The result of this study should lead to many aplications including short circuit, transient stability, network reduction, reactive optimization and others.

• 지구물리와물리탐사
• /
• 제6권2호
• /
• pp.95-100
• /
• 2003

본 연구에서는 복잡한 지질구조에 대해서도 신속하고 효율적으로 주시를 계산할 수 있는 Straight Ray Technique(SRT)을 이용한 반사주시 토모그래피 역산 알고리듬을 개발하였다. 역산을 위한 초기 속도모델은 지층경계면에 임피던스 변화를 갖는 상속도 모델을 사용하였다. 실제 속도모델의 반사주시와 초기 속도모델의 반사주시 차이를 계산하여 각각의 요소마다 주시의 오차를 줄이는 방법인 가우스-뉴튼 알고리듬을 이용하여 역산온 수행하였다. 자코비안의 요소는 파선이 지나가는 거리함수로 구성되며, 이를 최소자승형태의 근사 헤시안 행렬로 구성하여 역산을 수행하였다. 역산시 해가 수렴할 수 있도록 근사 헤시안 행렬의 대각성분에 일정한 감쇠인자를 더하였다. 역산된 속도모델을 이용하여 Kirchhoff구조보정을 실시한 결과 실제 속도모델구조에 근사한 단면영상을 얹을 수 있었다.

• 대한자원환경지질학회:학술대회논문집
• /
• 대한자원환경지질학회 2002년도 춘계 공동학술발표회
• /
• pp.167-169
• /
• 2002

The extended Born, or localized nonlinear approximation of integral equation (IE) solution has been applied to inverting single-hole electromagnetic (EM) data using a cylindrically symmetric model. The extended Born approximation is less accurate than a full solution but much superior to the simple Born approximation. When applied to the cylindrically symmetric model with a vertical magnetic dipole source, however, the accuracy of the extended Born approximation is greatly improved because the electric field is scalar and continuous everywhere. One of the most important steps in the inversion is the selection of a proper regularization parameter for stability. Occam's inversion (Constable et al., 1987) is an excellent method for obtaining a stable inverse solution. It is extremely slow when combined with a differential equation method because many forward simulations are needed but suitable for the extended Born solution because the Green's functions, the most time consuming part in IE methods, are repeatedly re-usable throughout the inversion. In addition, the If formulation also readily contains a sensitivity matrix, which can be revised at each iteration at little expense. The inversion algorithm developed in this study is quite stable and fast even if the optimum regularization parameter Is sought at each iteration step. Tn this paper we show inversion results using synthetic data obtained from a finite-element method and field data as well.