• Title/Summary/Keyword: Least squares solution

Search Result 206, Processing Time 0.022 seconds

Inversion of Acoustical Properties of Sedimentary Layers from Chirp Sonar Signals (Chirp 신호를 이용한 해저퇴적층의 음향학적 특성 역산)

  • 박철수;성우제
    • The Journal of the Acoustical Society of Korea
    • /
    • v.18 no.8
    • /
    • pp.32-41
    • /
    • 1999
  • In this paper, an inversion method using chirp signals and two near field receivers is proposed. Inversion problems can be formulated into the probabilistic models composed of signals, a forward model and noise. Forward model to simulate chirp signals is chosen to be the source-wavelet-convolution planewave modeling method. The solution of the inversion problem is defined by a posteriori pdf. The wavelet matching technique, using weighted least-squares fitting, estimates the sediment sound-speed and thickness on which determination of the ranges for a priori uniform distribution is based. The genetic algorithm can be applied to a global optimization problem to find a maximum a posteriori solution for determined a priori search space. Here the object function is defined by an L₂norm of the difference between measured and modeled signals. The observed signals can be separated into a set of two signals reflected from the upper and lower boundaries of a sediment. The separation of signals and successive applications of the genetic algorithm optimization process reduce the search space, therefore improving the inversion results. Not only the marginal pdf but also the statistics are calculated by numerical evaluation of integrals using the samples selected during importance sampling process of the genetic algorithm. The examples applied here show that, for synthetic data with noise, it is possible to carry out an inversion for sedimentary layers using the proposed inversion method.

  • PDF

Statistical review and explanation for Lanchester model (란체스터 모형에 대한 통계적 고찰과 해석)

  • Yoo, Byung Joo
    • The Korean Journal of Applied Statistics
    • /
    • v.33 no.3
    • /
    • pp.335-345
    • /
    • 2020
  • This paper deals with the problem of estimating the log-transformed linear regression model to fit actual battle data from the Ardennes Campaign of World War II into the Lanchester model. The problem of determining a global solution for parameters and multicollinearity problems are identified and modified by examining the results of previous studies on data. The least squares method requires attention because a local solution can be found rather than a global solution if considering a specific constraint or a limited candidate group. The method of exploring this multicollinearity problem can be confirmed by a statistic known as a variance inflation factor. Therefore, the Lanchester model is simplified to avoid these problems, and the combat power attrition rate model was proposed which is statistically significant and easy to explain. When fitting the model, the dependence problem between the data has occurred due to autocorrelation. Matters that might be underestimated or overestimated were resolved by the Cochrane-Orcutt method as well as guaranteeing independence and normality.

An Instantaneous Integer Ambiguity Resolution for GPS Real-Time Structure Monitoring (GPS 실시간 구조물 모니터링을 위한 반송파 관측데이터 순간미지정수 결정)

  • Lee, Hungkyu
    • KSCE Journal of Civil and Environmental Engineering Research
    • /
    • v.34 no.1
    • /
    • pp.341-353
    • /
    • 2014
  • In order to deliver a centimeter-level kinematic positioning solution with GPS carrier-phase measurements, it is prerequisite to use correctly resolved integer ambiguities. Based on the mathematical modeling of GPS network with application of its geometrical constraints, this research has investigated an instantaneous ambiguity resolution procedure for the so-called 'integer constrained least-squares' technique which can be effectively implemented in real-time structure monitoring. In this process, algorithms of quality control for the float solutions and hypothesis tests using the constrained baseline for the ambiguity validation are included to enhance reliability of the solutions. The proposed procedure has been implemented by MATLAB, the language of technical computing, and processed field trial data obtained at a cable-stayed bridge to access its real-world applicability. The results are summarized in terms of ambiguity successful rates, impact of the stochastical models, and computation time to demonstrate performance of the instantaneous ambiguity resolution proposed.

GPS Baseline Estimation of the $2^{nd}$ Order Geodetic Control Network (2등 측지기준점 GPS 관측데이터의 기선벡터 추정)

  • Lee, Young-Jin;Lee, Hung-Kyu;Kwon, Chan-Oh;Cha, Sang-Heon
    • Journal of the Korean Society of Surveying, Geodesy, Photogrammetry and Cartography
    • /
    • v.26 no.2
    • /
    • pp.157-164
    • /
    • 2008
  • GPS baseline analysis is a mathematical procedure which estimates a baseline vector from carrier-phase double-differenced observations. Least squares technique is generally applied for the processing and integer ambiguities in the observations should be resolved to obtain maximum accuracy of the solution. In GPS control surveying, after assembling the baseline solutions into a network, adjustment is performed to derive final coordinate sets of unknown points. This paper deals with details of GPS baseline analysis for the $2^{nd}$ order national geodetic network adjustment. After reviewing GPS campaigns carried out by National Geographic Information Institute (NGII) and their observations. technical issues and considerations for the GPS baseline analysis are presented with emphasis of selecting the processing strategies and software. Finally, the analyzed results will be evaluated by examining the close of figures formed by joining the processed baseline vectors.

The Kinematic Comparison of Energy Walking and Normal Walking (에너지보행과 일반보행의 운동학적 비교)

  • Shin, Je-Min;Jin, Young-Wan
    • Korean Journal of Applied Biomechanics
    • /
    • v.16 no.4
    • /
    • pp.61-71
    • /
    • 2006
  • The purpose of this study was to compare kinematic characteristics on the limbs at 3 different walking speed during the energy and the normal walking. Eight subjects performed energy walking and normal walking at the slow speed(65 beats/min), the normal speed(115beats/min), the fast speed(160 beats/min). The 3-d angle was calculated by vector projected with least squares solution with three-dimensional cinematography(Motion Analysis corporation). The range of motion was calculated on the trunk, shoulder, elbow, hip, knee joint. The results showed that stride length was no difference of the two walking pattern. The duration of support phase was also no difference of the two walking pattern. The range of motion of shoulder joint significantly increased in the sagittal and frontal planes, and the range of motion of elbow joint significantly increased as the energy walking. The range of motion of hip joint had no significant difference in the any planes in changing of walking speed. But the most remarkable difference of the two walking patterns revealed at the trunk. The range of flexion/extension angle had significant increasing $2.36^{\circ}$ at normal speed, and the range of the right/left flexion angle had significant increasing below $4^{\circ}$ at the 3 walking speed, and The range of rotation angle had significant increasing $7.35^{\circ}$, $9.22^{\circ}$, respectively at the normal and slow speed. But there was no significant difference of range of motion at the hip and knee joints between energy walking and normal walking.

Direction-of-Arrival Estimation in Broadband Signal Processing : Rotation of Signal Subspace Approach (광대역 신호 처리에서의 도래각 추정 : Rotation of Signal Subspaces 방법)

  • Kim, Young-Soo
    • Journal of the Korean Institute of Telematics and Electronics
    • /
    • v.26 no.7
    • /
    • pp.166-175
    • /
    • 1989
  • In this paper, we present a method which is based on the concept of the rotation of subspaces. This method is highly related to the angle (or distance) between subspaces arising in many applications. An effective procedures is first derived for finding the optimal transformation matrix which rotates one subspace into another as closely as possible in the least squares sense , and then this algorithm is applied to the solution to general direction-of-arrival estimation problem of multiple broadband plane waves which may be a mixture of incoherent, partially coherent or coherent. In this typical application, the rotation of signal subspaces (ROSS) algorithm is effectively developed to achieve the high performance in the active systems for the case in which the noise field remains invariant with the measurement of the array spectral density matrix (or data matrix). It is not uncommon to observe this situation in sonar systems. The advantage of this techniques is not to require the preliminary processing and spatial prefiltering which is used in Wang-Kaveh's CSS focusing method. Furthermore, the array's geometry is not restricted. Simulation results are presented to illustrate the high performance achieved with this new approach relative to that obtained with Wang-Kaveh's CSS focusing method for incoherent sources and forward-backward spatial smoothed MUSIC for coherent sources including the signal eigenvector method (SEM).

  • PDF

An Efficient Adaptive Digital Filtering Algorithm for Identification of Second Order Volterra Systems (이차 볼테라 시스템 인식을 위한 효율적인 적응 디지탈 필터링 알고리즘)

  • Hwang, Y.S.;Mathews, V.J.;Cha, I.W.;Youn, D.H.
    • The Journal of the Acoustical Society of Korea
    • /
    • v.7 no.4
    • /
    • pp.98-109
    • /
    • 1988
  • This paper introduces an adaptive nonlinear filtering algorithm that uses the sequential regression(SER) method to update the second order Volterra filter coefficients in a recursive way. Conventionally, the SER method has been used to invert large matrices which result from direct application of Wiener filter theory to the Volterra filter. However, the algorithm proposed in this paper uses the SER approach to update the least squares solution which is derived for Gaussian input signals. In such an algorithm, the size of the matrix to be inverted is smaller than that of conventional approaches, and hence the proposed method is computationally simpler than conventional nonlinear system identification techniques. Simulation results are presented to demonstrate the performance of the proposed algorithm.

  • PDF

Development of MLS Difference Method for Material Nonlinear Problem (MLS차분법을 이용한 재료비선형 문제 해석)

  • Yoon, Young-Cheol
    • Journal of the Computational Structural Engineering Institute of Korea
    • /
    • v.29 no.3
    • /
    • pp.237-244
    • /
    • 2016
  • This paper presents a nonlinear Moving Least Squares(MLS) difference method for material nonlinearity problem. The MLS difference method, which employs strong formulation involving the fast derivative approximation, discretizes governing partial differential equation based on a node model. However, the conventional MLS difference method cannot explicitly handle constitutive equation since it solves solid mechanics problems by using the Navier's equation that unifies unknowns into one variable, displacement. In this study, a double derivative approximation is devised to treat the constitutive equation of inelastic material in the framework of strong formulation; in fact, it manipulates the first order derivative approximation two times. The equilibrium equation described by the divergence of stress tensor is directly discretized and is linearized by the Newton method; as a result, an iterative procedure is developed to find convergent solution. Stresses and internal variables are calculated and updated by the return mapping algorithm. Effectiveness and stability of the iterative procedure is improved by using algorithmic tangent modulus. The consistency of the double derivative approximation was shown by the reproducing property test. Also, accuracy and stability of the procedure were verified by analyzing inelastic beam under incremental tensile loading.

Performance Evaluation of a Time-domain Gauss-Newton Full-waveform Inversion Method (시간영역 Gauss-Newton 전체파형 역해석 기법의 성능평가)

  • Kang, Jun Won;Pakravan, Alireza
    • Journal of the Computational Structural Engineering Institute of Korea
    • /
    • v.26 no.4
    • /
    • pp.223-231
    • /
    • 2013
  • This paper presents a time-domain Gauss-Newton full-waveform inversion method for the material profile reconstruction in heterogeneous semi-infinite solid media. To implement the inverse problem in a finite computational domain, perfectly-matchedlayers( PMLs) are introduced as wave-absorbing boundaries within which the domain's wave velocity profile is to be reconstructed. The inverse problem is formulated in a partial-differential-equations(PDE)-constrained optimization framework, where a least-squares misfit between measured and calculated surface responses is minimized under the constraint of PML-endowed wave equations. A Gauss-Newton-Krylov optimization algorithm is utilized to iteratively update the unknown wave velocity profile with the aid of a specialized regularization scheme. Through a series of one-dimensional examples, the solution of the Gauss-Newton inversion was close enough to the target profile, and showed superior convergence behavior with reduced wall-clock time of implementation compared to a conventional inversion using Fletcher-Reeves optimization algorithm.

Derivation of An Empirical Formula for Determining Water Content of Mixed Uranyl Nitrate-Thorium Nitrate Solutions (질산(窒酸)우라늄-질산(窒酸)토륨 혼합용액중(混合溶液中)의 물함량(含量) 결정식(決定式) 유도(誘導))

  • Min, Duck-Kee;Choi, Byung-Il;Ro, Seung-Gy;Eom, Tae-Yoon;Kim, Zong-Goo
    • Journal of Radiation Protection and Research
    • /
    • v.11 no.1
    • /
    • pp.22-28
    • /
    • 1986
  • Densities of a large number of mixed uranyl nitrate-thorium nitrate solutions were measured with pycnometer. By the least squares analysis of the experimental result, an empirical formula for determining water content of mixed uranyl nitrate-thorium nitrate solutions as functions of uranium concentration, thorium concentration and nitric acid normality is derived; $W=1.0-0.358\;C_u-0.4538\;C_{Th}-0.0307\;H^+$ where $W,\;C_u,\;C_{Th},\;and\;H^+$ stand for water content(g/cc), uranium concentration (g/cc), thorium concentration (g/cc), and nitric acid normality, respectively. Water contents of the mixed uranyl nitrate-thorium nitrate solutions are calculated by using the empirical formula, and compared with the values calculated by Bouly's equation in which an additional data, solution density, is required. The two results show good agreements within 2.7%.

  • PDF