• Title/Summary/Keyword: Newton-Raphson algorithm

Search Result 173, Processing Time 0.027 seconds

A Variable Latency Newton-Raphson's Floating Point Number Reciprocal Square Root Computation (가변 시간 뉴톤-랍손 부동소수점 역수 제곱근 계산기)

  • Kim Sung-Gi;Cho Gyeong-Yeon
    • The KIPS Transactions:PartA
    • /
    • v.12A no.5 s.95
    • /
    • pp.413-420
    • /
    • 2005
  • The Newton-Raphson iterative algorithm for finding a floating point reciprocal square mot calculates it by performing a fixed number of multiplications. In this paper, a variable latency Newton-Raphson's reciprocal square root algorithm is proposed that performs multiplications a variable number of times until the error becomes smaller than a given value. To find the rediprocal square root of a floating point number F, the algorithm repeats the following operations: '$X_{i+1}=\frac{{X_i}(3-e_r-{FX_i}^2)}{2}$, $i\in{0,1,2,{\ldots}n-1}$' with the initial value is '$X_0=\frac{1}{\sqrt{F}}{\pm}e_0$'. The bits to the right of p fractional bits in intermediate multiplication results are truncated and this truncation error is less than '$e_r=2^{-p}$'. The value of p is 28 for the single precision floating point, and 58 for the double precision floating point. Let '$X_i=\frac{1}{\sqrt{F}}{\pm}e_i$, there is '$X_{i+1}=\frac{1}{\sqrt{F}}-e_{i+1}$, where '$e_{i+1}{<}\frac{3{\sqrt{F}}{{e_i}^2}}{2}{\mp}\frac{{Fe_i}^3}{2}+2e_r$'. If '$|\frac{\sqrt{3-e_r-{FX_i}^2}}{2}-1|<2^{\frac{\sqrt{-p}{2}}}$' is true, '$e_{i+1}<8e_r$' is less than the smallest number which is representable by floating point number. So, $X_{i+1}$ is approximate to '$\frac{1}{\sqrt{F}}$. Since the number of multiplications performed by the proposed algorithm is dependent on the input values, the average number of multiplications Per an operation is derived from many reciprocal square root tables ($X_0=\frac{1}{\sqrt{F}}{\pm}e_0$) with varying sizes. The superiority of this algorithm is proved by comparing this average number with the fixed number of multiplications of the conventional algorithm. Since the proposed algorithm only performs the multiplications until the error gets smaller than a given value, it can be used to improve the performance of a reciprocal square root unit. Also, it can be used to construct optimized approximate reciprocal square root tables. The results of this paper can be applied to many areas that utilize floating point numbers, such as digital signal processing, computer graphics, multimedia, scientific computing, etc.

Computer simulation system of robot manipulator motion (로보트 매니퓰레이터 운동의 컴퓨터 시뮬레이션 시스템)

  • 김창부;윤장로
    • 제어로봇시스템학회:학술대회논문집
    • /
    • 1991.10a
    • /
    • pp.539-544
    • /
    • 1991
  • In order to verify robot motions for a desired work, it is necessary to visualize it on a computer screen. This paper presents a simulation algorithm for robot manipulator motion. Kinematic description is based on the Denavit- Hartenberg link representation. In order to be applied to various types of the robot manipulator, inverse kinematics make use of the Newton-Raphson iterative method with the least squares method. Joint variables are interpolated by the lowest polynomial segment satisfying acceleration continuity. The robot motions are generated and then animated on a computer screen in the form of skeleton type.

  • PDF

Structural Damage Identification by Using the Structurally Damped Spectral Element Model (구조감쇠가 고려된 스펙트럴요소 모델을 이용한 구조손상규명)

  • Kim, Jung-Soo;Cho, Joo-Yong;Lee, U-Sik
    • Proceedings of the KSR Conference
    • /
    • 2004.10a
    • /
    • pp.121-126
    • /
    • 2004
  • In this paper, a nonlinear structural damage identification algorithm is derived by taking into account the structurally damped spectral element model thinking over a real situation. The structural damage identification analyses are conducted by using the Newton-Raphson method. It is found that, in general Structural Damage Identification by using the Structurally Damped Spectral Element Model provides the same exact damage identification results when compared with the results obtained by the structurally undamped spectral model.

  • PDF

Inhomogeneous Poisson Intensity Estimation via Information Projections onto Wavelet Subspaces

  • Kim, Woo-Chul;Koo, Ja-Yong
    • Journal of the Korean Statistical Society
    • /
    • v.31 no.3
    • /
    • pp.343-357
    • /
    • 2002
  • This paper proposes a method for producing smooth and positive estimates of the intensity function of an inhomogeneous Poisson process based on the shrinkage of wavelet coefficients of the observed counts. The information projection is used in conjunction with the level-dependent thresholds to yield smooth and positive estimates. This work is motivated by and demonstrated within the context of a problem involving gamma-ray burst data in astronomy. Simulation results are also presented in order to show the performance of the information projection estimators.

FINDING A ZERO OF THE SUM OF TWO MAXIMAL MONOTONE OPERATORS WITH MINIMIZATION PROBLEM

  • Abdallah, Beddani
    • Nonlinear Functional Analysis and Applications
    • /
    • v.27 no.4
    • /
    • pp.895-902
    • /
    • 2022
  • The aim of this paper is to construct a new method for finding the zeros of the sum of two maximally monotone mappings in Hilbert spaces. We will define a simple function such that its set of zeros coincide with that of the sum of two maximal monotone operators. Moreover, we will use the Newton-Raphson algorithm to get an approximate zero. In addition, some illustrative examples are given at the end of this paper.

Classification of Ultrasonic NDE Signals Using the Expectation Maximization (EM) and Least Mean Square (LMS) Algorithms (최대 추정 기법과 최소 평균 자승 알고리즘을 이용한 초음파 비파괴검사 신호 분류법)

  • Kim, Dae-Won
    • Journal of the Korean Society for Nondestructive Testing
    • /
    • v.25 no.1
    • /
    • pp.27-35
    • /
    • 2005
  • Ultrasonic inspection methods are widely used for detecting flaws in materials. The signal analysis step plays a crucial part in the data interpretation process. A number of signal processing methods have been proposed to classify ultrasonic flaw signals. One of the more popular methods involves the extraction of an appropriate set of features followed by the use of a neural network for the classification of the signals in the feature spare. This paper describes an alternative approach which uses the least mean square (LMS) method and exportation maximization (EM) algorithm with the model based deconvolution which is employed for classifying nondestructive evaluation (NDE) signals from steam generator tubes in a nuclear power plant. The signals due to cracks and deposits are not significantly different. These signals must be discriminated to prevent from happening a huge disaster such as contamination of water or explosion. A model based deconvolution has been described to facilitate comparison of classification results. The method uses the space alternating generalized expectation maximiBation (SAGE) algorithm ill conjunction with the Newton-Raphson method which uses the Hessian parameter resulting in fast convergence to estimate the time of flight and the distance between the tube wall and the ultrasonic sensor. Results using these schemes for the classification of ultrasonic signals from cracks and deposits within steam generator tubes are presented and showed a reasonable performances.

An Analysis Algorithm to Overcome the Singularity of Time Integrations for Dynamics Problems (동역학 문제의 시간적분 특이성을 극복하기 위한 해석 알고리듬)

  • 엄기상;윤성호
    • Proceedings of the Computational Structural Engineering Institute Conference
    • /
    • 2004.04a
    • /
    • pp.1-8
    • /
    • 2004
  • For the linearized differential algebraic equation of the nonlinear constrained system, exact initial values of the acceleration are needed to solve itself. It may be very troublesome to perform the inverse operation for obtaining the incremental quantities since the mass matrix contains the zero element in the diagonal. This fact makes the mass matrix impossible to be positive definite. To overcome this singularity phenomenon the mass matrix needs to be modified to allow the feasible application of predictor and corrector in the iterative computation. In this paper the proposed numerical algorithm based on the modified mass matrix combines the conventional implicit algorithm, Newton-Raphson method and Newmark method. The numerical example presents reliabilities for the proposed algorithm via comparisons of the 4th order Runge-kutta method. The proposed algorithm seems to be satisfactory even though the acceleration, Lagrange multiplier, and energy show unstable behaviour. Correspondingly, it provides one important clue to another algorithm for the enhancement of the numerical results.

  • PDF

New Power Flow Calculation Using Improved Genetic Algorithm (개선된 유전 알고리즘을 이용한 새로운 전력조류계산)

  • Chae, Myung-Suck;Lee, Tae-Hyung;Shin, Joong-Rin;Im, Han-Seok
    • Proceedings of the KIEE Conference
    • /
    • 1999.11a
    • /
    • pp.43-51
    • /
    • 1999
  • The power flow calculations(PFc) are the most important and powerful tools in power systems engineering. The conventional power flow problem is solved generally with numerical methods such as Newton-Raphson(NR). The conventional numerical method generally have some convergency problem, which is sensitive to initial value, and numerical stability problem concerned with jacobian matrix inversion. This paper presents a new PFc algorithm based on the improved genetic algorithm (IGA) which can overcome the disadvantages mentioned above. The parameters of GA, with dynamical hierarchy of the coding system, are improved to make GA a practical algorithm in the problem of real system. Some case studies with test bus system also present to show the performance of proposed algorithm. The results of proposed algorithm are compared with the results of PFc obtained using a conventional NR method.

  • PDF

Intelligent Optimization Algorithm Approach to Image Reconstruction in Electrical Impedance Tomography (지능 최적 알고리즘을 이용한 전기임피던스 단층촬영법의 영상복원)

  • Kim, Ho-Chan;Boo, Chang-Jin;Lee, Yoon-Joon
    • Proceedings of the KIEE Conference
    • /
    • 2002.11c
    • /
    • pp.513-516
    • /
    • 2002
  • In electrical impedance tomography(EIT), various image reconstruction algorithms have been used in order to compute the internal resistivity distribution of the unknown object with its electric potential data at the boundary. Mathematically the EIT image reconstruction algorithm is a nonlinear ill-posed inverse problem. This paper presents two intelligent optimization algorithm techniques such as genetic algorithm and simulated annealing for the solution of the static EIT inverse problem. We summarize the simulation results for the three algorithm forms: modified Newton-Raphson, genetic algorithm, and simulated annealing.

  • PDF

Partitioning likelihood method in the analysis of non-monotone missing data

  • Kim Jae-Kwang
    • Proceedings of the Korean Statistical Society Conference
    • /
    • 2004.11a
    • /
    • pp.1-8
    • /
    • 2004
  • We address the problem of parameter estimation in multivariate distributions under ignorable non-monotone missing data. The factoring likelihood method for monotone missing data, termed by Robin (1974), is extended to a more general case of non-monotone missing data. The proposed method is algebraically equivalent to the Newton-Raphson method for the observed likelihood, but avoids the burden of computing the first and the second partial derivatives of the observed likelihood Instead, the maximum likelihood estimates and their information matrices for each partition of the data set are computed separately and combined naturally using the generalized least squares method. A numerical example is also presented to illustrate the method.

  • PDF