• International Journal of Control, Automation, and Systems
• /
• v.6 no.2
• /
• pp.160-170
• /
• 2008

Temperature control of an enclosed thermal system which has many applications including Rapid Thermal Processing (RTP) of semiconductor wafers showed an input-constraint violation for nonlinear controllers due to inherent strong coupling between the elements [1]. In this paper, a constrained nonlinear optimal control design is developed, which accommodates input constraints using the linear algebraic equivalence of the nonlinear controllers, for the temperature control of an enclosed thermal process. First, it will be shown that design of nonlinear controllers is equivalent to solving a set of linear algebraic equations-the linear algebraic equivalence of nonlinear controllers (LAENC). Then an input-constrained nonlinear optimal controller is designed based on that LAENC using the constrained linear least squares method. Through numerical simulations, it is demonstrated that the proposed controller achieves the equivalent performances to the classical nonlinear controllers with less total energy consumption. Moreover, it generates the practical control solution, in other words, control solutions do not violate the input-constraints.

• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• v.24 no.3
• /
• pp.243-291
• /
• 2020

In this paper, multi-block generalized backward differentiation methods for numerical solutions of ordinary differential and differential algebraic equations are introduced. This class of linear multi-block methods is implemented as multi-block boundary value methods (MB₂ VMs). The root distribution of the stability polynomial of the new class of methods are determined using the Wiener-Hopf factorization of a matrix polynomial for the purpose of their correct implementation. Numerical tests, showing the potential of such methods for output of multi-block of solutions of the ordinary differential equations in the new approach are also reported herein. The methods which output multi-block of solutions of the ordinary differential equations on application, are unlike the conventional linear multistep methods which output a solution at a point or the conventional boundary value methods and multi-block methods which output only a block of solutions per step. The MB₂ VMs introduced herein is a novel approach at developing very large scale integration methods (VLSIM) in the numerical solution of differential equations.

• 제어로봇시스템학회:학술대회논문집
• /
• 1995.10a
• /
• pp.432-435
• /
• 1995

In this paper a new approach to obtain the solution of the linear-quadratic Gaussian control problem for singularly perturbed discrete-time stochastic systems is proposed. The alogorithm proposed is based on exploring the previous results that the exact solution of the global discrete algebraic Riccati equations is found in terms of the reduced-order pure-slow and pure-fast nonsymmetric continuous-time algebraic Riccati equations and, in addition, the optimal global Kalman filter is decomposed into pure-slow and pure-fast local optimal filters both driven by the system measurements and the system optimal control input. It is shown that the optimal linear-quadratic Gaussian control problem for singularly perturbed linear discrete systems takes the complete decomposition and parallelism between pure-slow and pure-fast filters and controllers.

• Kyungpook Mathematical Journal
• /
• v.28 no.2
• /
• pp.185-196
• /
• 1988

We study existence of polynomial integrating factors and solutions F(x, y)=c of first order nonlinear differential equations. We characterize the homogeneous case, and give algorithms for finding existence of and a basis for polynomial solutions of linear difference and differential equations and rational solutions or linear differential equations with polynomial coefficients. We relate singularities to nature of the solution. Solution of differential equations in closed form to some degree might be called more an art than a science: The investigator can try a number of methods and for a number of classes of equations these methods always work. In particular integrating factors are tricky to find. An analogous but simpler situation exists for integrating inclosed form, where for instance there exists a criterion for when an exponential integral can be found in closed form. In this paper we make a beginning in several directions on these problems, for 2 variable ordinary differential equations. The case of exact differentials reduces immediately to quadrature. The next step is perhaps that of a polynomial integrating factor, our main study. Here we are able to provide necessary conditions based on related homogeneous equations which probably suffice to decide existence in most cases. As part of our investigations we provide complete algorithms for existence of and finding a basis for polynomial solutions of linear differential and difference equations with polynomial coefficients, also rational solutions for such differential equations. Our goal would be a method for decidability of whether any differential equation Mdx+Mdy=0 with polynomial M, N has algebraic solutions(or an undecidability proof). We reduce the question of all solutions algebraic to singularities but have not yet found a definite procedure to find their type. We begin with general results on the set of all polynomial solutions and integrating factors. Consider a differential equation Mdx+Ndy where M, N are nonreal polynomials in x, y with no common factor. When does there exist an integrating factor u which is (i) polynomial (ii) rational? In case (i) the solution F(x, y)=c will be a polynomial. We assume all functions here are complex analytic polynomial in some open set.

• International Journal of Automotive Technology
• /
• v.6 no.4
• /
• pp.375-381
• /
• 2005

This paper presents linear algebraic equations in the form of recursive formula to compute elastokinematic characteristics of a suspension system. Conventional methods of elastokinematic analysis are based on nonlinear kinematic constrant equations and force equilibrium equations for constrained mechanical systems, which require complicated and time-consuming implicit computing methods to obtain the solution. The proposed linearized elastokinematic equations in the form of recursive formula are derived based on the assumption that the displacements of elastokinematic behavior of a constrained mechanical system under external forces are very small. The equations can be easily computerized in codes, and have the advantage of sharing the input data of existing general multi body dynamic analysis codes. The equations can be applied to any form of suspension once the type of kinematic joints and elastic components are identified. The validity of the method has been proved through the comparison of the results from established elastokinematic analysis software. Error estimation and analysis due to piecewise linear assumption are also discussed.

• Coupled systems mechanics
• /
• v.1 no.2
• /
• pp.155-163
• /
• 2012

This paper presents analytical approximate solutions for the initial post-buckling deformation of the pipes in oil and gas wells. The governing differential equation with sinusoidal nonlinearity can be reduced to form a third-order-polynomial nonlinear equation, by coupling of the well-known Maclaurin series expansion and orthogonal Chebyshev polynomials. Analytical approximations to the resulting boundary condition problem are established by combining the Newton's method with the method of harmonic balance. The linearization is performed prior to proceeding with harmonic balancing thus resulting in a set of linear algebraic equations instead of one of non-linear algebraic equations, unlike the classical method of harmonic balance. We are hence able to establish analytical approximate solutions. The approximate formulae for load along axis, and periodic solution are established for derivative of the helix angle at the end of the pipe. Illustrative examples are selected and compared to "reference" solution obtained by the shooting method to substantiate the accuracy and correctness of the approximate analytical approach.

• Journal of Electrochemical Science and Technology
• /
• v.7 no.4
• /
• pp.293-297
• /
• 2016

Here we propose a new equation by which we can estimate the baseline for measuring the peak current of the reverse curve in a cyclic voltammogram. A similar equation already exists, but it is a linear algebraic equation that over-simplifies the voltammetric curve and may cause unpredictable errors when calculating the baseline. In our study, we find a quadratic algebraic equation that acceptably reflects the complexity included in a voltammetric curve. The equation is obtained from a laborious numerical analysis of cyclic voltammetry simulations using the finite element method, and not from the closed form of the mathematical equation. This equation is utilized to provide a virtual baseline current for the reverse peak current. We compare the results obtained using the old linear and new quadratic equations with the theoretical values in terms of errors to ascertain the degree to which accuracy is improved by the new equation. Finally, the equations are applied to practical cyclic voltammograms of ferricyanide in order to confirm the improved accuracy.

• Journal of the Korean Society for Precision Engineering
• /
• v.21 no.7
• /
• pp.76-82
• /
• 2004

The solutions of neighboring optimal control are typically obtained using the sweep method or transition matrices. Due to the numerical integration, however, the gain matrix can become infinite as time go to final one in the transition matrices, and the Riccati solution can become infinite when the final time free. To overcome these disadvantages, this paper proposes the pseudospectral Legendre method which is to first discreteize the linear boundary value problem using the global orthogonal polynomial, then transforms into an algebraic equations. Because this method is not necessary to take any integration of transition matrix or Riccati equation, it can be usefully used in real-time operation. Finally, its performance is verified by the numerical example for the space vehicle's orbit transfer.

• Kyungpook Mathematical Journal
• /
• v.56 no.4
• /
• pp.1141-1160
• /
• 2016

In this paper, a symbolic algorithm for solving a regular initial value problem (IVP) for a system of linear differential-algebraic equations (DAEs) with constant coeffcients has been presented. Algebra of integro-differential operators is employed to express the given system of DAEs. We compute a canonical form of the given system which produces another simple equivalent system. Algorithm includes computing the matrix Green's operator and the vector Green's function of a given IVP. Implementation of the proposed algorithm in Maple is also presented with sample computations.

• Journal of Korean Society of Forest Science
• /
• v.90 no.2
• /
• pp.210-216
• /
• 2001

Pinus densiflora S. et Z. has widely been distributed, and is one of the important main foret resources in Korea. Diameter and height growth patterns were estimated using non-linear algebraic difference equation, which requires two-measurement times $T_1$ and $T_2$. To maximize data use, all possible measurement interval data were derived using Lag and Put statements in the SAS. In results, of the algebraic difference equations applied, the Schumacher and the Gompertz polymorphic equations for diameter and height, respectively showed the higher precision of the fitting. In order to allow more precise estimation of growth than those of the basic Schumacher and the Gompertz, further refinement that combine biological realism as input into the equation would be necessary.