• Geophysics and Geophysical Exploration
• /
• v.10 no.2
• /
• pp.124-130
• /
• 2007

Least squares ($l^2$ norm) solutions of seismic inversion tend to be very sensitive to data points with large errors. The $l^1$ norm minimization gives more robust solutions, but usually with higher computational cost. Iteratively reweighted least squares (IRLS) method gives efficient approximate solutions of these $l^1$ norm problems. I propose an efficient implementation of the IRLS method for a hybrid $l^1/l^2$ minimization problem that behaves as $l^2$ norm fit for small residual and $l^1$ norm fit for large residuals. The proposed algorithm shows more robust characteristics to the decision of the threshold value than the l1 norm IRLS inversion does with respect to the threshold value to avoid singularity.

• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• v.14 no.2
• /
• pp.125-140
• /
• 2010

This paper develops a least-squares approach to the solution of the optimal control problem for the Navier-Stokes equations. We recast the optimality system as a first-order system by introducing velocity-flux variables and associated curl and trace equations. We show that a least-squares principle based on $L^2$ norms applied to this system yields optimal discretization error estimates in the $H^1$ norm in each variable.

• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• v.11 no.4
• /
• pp.55-68
• /
• 2007

First-order least-squares method of a distributed optimal control problem for the incompressible Navier-Stokes equations is considered. An optimality system for the optimal solution are reformulated to the equivalent first-order system by introducing velocity-flux variables and then the least-squares functional corresponding to the system is defined in terms of the sum of the squared $L^2$ norm of the residual equations of the system. The optimal error estimates for least-squares finite element approximations are obtained.

• The Journal of Korean Institute of Communications and Information Sciences
• /
• v.29 no.2C
• /
• pp.272-282
• /
• 2004

This paper introduces a regularized mixed norm multi-channel image restoration algorithm using both within-and between- channel deterministic information. For each channel a functional which combines the least mean squares (LMS), the least mean fourth (LMF), and a smoothing functional is proposed. We introduce a mixed norm parameter that controls the relative contribution between the LMS and the LMF, and a regularization parameter defining the degree of smoothness of the solution, where both parameters are updated at each iteration according to the noise characteristics of each channel. The novelty of the proposed algorithm is that no knowledge of the noise distribution for each channel is required and that the parameters mentioned above are adjusted based on the partially restored image.

• Proceedings of the IEEK Conference
• /
• 2000.09a
• /
• pp.399-402
• /
• 2000

In this paper, we present a regularized mixed norm multichannel image restoration algorithm. The problem of multichannel restoration using both within- and between- channel deterministic information is considered. For each channel a functional which combines the least mean squares (LMS), the least mean fourth(LMF), and a smoothing functional is proposed, We introduce a mixed norm parameter that controls the relative contribution between the LMS and the LMF, and a regularization parameter that defines the degree of smoothness of the solution, both updated at each iteration according to the noise characteristics of each channel. The novelty of the proposed algorithm is that no knowledge of the noise distribution for each channel is required, and the parameters mentioned above are adjusted based on the partially restored image.

• The Journal of Korea Robotics Society
• /
• v.10 no.4
• /
• pp.230-244
• /
• 2015

Continuous-path motion control such as resolved motion rate control requires online solving of the inverse differential kinematics for a robot. However, the solution space of the inverse differential kinematics related to Jacobian J is not well-established. In this paper, the solution space of inverse differential kinematics is analyzed through categorization of mapping conditions between joint velocities and end-effector velocity of a robot. If end-effector velocity is within the column space of J, the solution or the minimum norm solution is obtained. If it is not within the column space of J, an approximate solution by least-squares is obtained. Moreover, this paper introduces an improved mapping diagram showing orthogonality and mapping clearly between subspaces, and concrete examples numerically showing the concept of several subspaces. Finally, a solver and graphics user interface (GUI) for inverse differential kinematics are developed using MATLAB, and the solution of inverse differential kinematics using the GUI is demonstrated for a vertically articulated robot.

• Communications of the Korean Mathematical Society
• /
• v.10 no.1
• /
• pp.235-249
• /
• 1995

Let K be a bounded linear operator from Hilbert space $H_1$ into Hilbert space $H_2$. When numerically solving the first kind equation Kf = g, one usually picks n orthonormal functions $\phi_1, \phi_2,...,\phi_n$ in $H_1$ which depend on the numerical method and on the problem, see Varah [12] for more details. Then one findes the unique minimum norm element $f_M \in M$ that satisfies $\Vert K f_M - g \Vert = inf {\Vert K f - g \Vert : f \in M}$, where M is the linear span of $\phi_1, \phi_2,...,\phi_n$. Such a solution $f_M \in M$ is called the M-solution of K f = g. Some methods for finding the M-solution of K f = g were proposed by Banks [2] and Marti [9,10]. See [5,6,8] for convergence results comparing the M-solution of K f = g with $f_0$, the least squares solution of minimum norm (LSSMN) of K f = g.

• The Journal of Korean Institute of Communications and Information Sciences
• /
• v.33 no.8C
• /
• pp.583-590
• /
• 2008

In this paper, we present a spatially adaptive post-processing algorithm for H.264 video coding standard to remove blocking artifacts. The loop filter of H.264 increases computational complexity of the encoder. Furthermore it doesn't clearly remove the blocking artifacts, resulting in over-blurring. For overcoming them, we combine the projection method with the Constraint Least Squares(CLS) method to restore the high quality image. To reflect the Human Visual System, we adopt the weight norm CLS method. Particularly pixel location-based local variance and laplacian operator are newly defined for the CLS method. In addition, the fact that correlation among adjoining pixels is high is utilized to constrain the solution space when the projection method is applied. Quantization Index(QP) of H.264 is also used to control the degree of smoothness. The simulation results show that the proposed post-processing filter works better than the loop filter of H.264 and converges more quickly than the CLS method.

• 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.