Inverse Solution to the Coupled Ordinary Differential Equation

Author(s):  
P. Venkataraman

The identification of the actual form of the constant coefficient coupled differential equations and their boundary conditions, from two sets of discrete data points, is possible through a unique two-step approach developed in this paper. In the first step, the best Bezier function is fitted to the data. This allows an effective approximation of the data and the required number of derivatives for the entire range of the independent variable. In the second step, the known derivatives are introduced in a generic model of the coupled differential equation. This generic form includes two types of unknowns, real numbers and integers. The real numbers are the coefficients of the various terms in the differential equations, while the integers are exponents of the derivatives. The unknown exponents and coefficients are identified using an error formulation. Two examples are solved. The given data is exact, smooth and they represent solutions to coupled linear differential equations. The solution is obtained through discrete programming. Three methods are presented. The first is limited enumeration, which is useful if the coefficients belong to a limited set of discrete values. The second is global search using the genetic algorithm for a larger choice of coefficient values. The third uses a state space integrator driven by the genetic algorithm, to minimize the error between known data and that obtained from numerical integration.

Author(s):  
P. Venkataraman

A challenging inverse problem is to identify the smooth function and the differential equation it represents from uncertain data. This paper extends the procedure previously developed for smooth data. The approach involves two steps. In the first step the data is smoothed using a recursive Bezier filter. For smooth data a single application of the filter is sufficient. The final set of data points provides a smooth estimate of the solution. More importantly, it will also identify smooth derivatives of the function away from the edges of the domain. In the second step the values of the function and its derivatives are used to establish a specific form of the differential equation from a particular class of the same. Since the function and its derivatives are known, the only unknowns are parameters describing the structure of the differential equations. These parameters are of two kinds: the exponents of the derivatives and the coefficients of the terms in the differential equations. These parameters can be determined by defining an optimization problem based on the residuals in a reduced domain. To avoid the trivial solution a discrete global search is used to identify these parameters. An example involving a third order constant coefficient linear differential equation is presented. A basic simulated annealing algorithm is used for the global search. Once the differential form is established, the unknown initial and boundary conditions can be obtained by backward and forward numerical integration from the reduced region.


1955 ◽  
Vol 51 (4) ◽  
pp. 604-613
Author(s):  
Chike Obi

1·1. A general problem in the theory of non-linear differential equations of the second order is: Given a non-linear differential equation of the second order uniformly almost periodic (u.a.p.) in the independent variable and with certain disposable constants (parameters), to find: (i) the non-trivial relations between these parameters such that the given differential equation has a non-periodic u.a.p. solution; (ii) the number of periodic and non-periodic u.a.p. solutions which correspond to each such relation; and (iii) explicit analytical expressions for the u.a.p. solutions when they exist.


Author(s):  
P. Venkataraman

The simplest inverse boundary value problem is to identify the differential equation and the boundary conditions from a given set of discrete data points. For an ordinary differential equation, it would involve finding a function, which when expressed through some function of itself and its derivatives, and integrated using particular boundary conditions would generate the given data. Parametric Bezier functions are excellent candidates for these functions. They allow efficient approximation of data and its derivative content. The Bezier function is smooth and continuous to a high degree. In this paper the best Bezier function to fit the data represents this function which is being sought. This Bezier approximation also determines the boundary conditions. Next, a generic form of the differential equation is assumed. The Bezier function and its derivatives are then used in this generic form to establish the exponents and coefficients of the various terms in the actual differential equation. The paper looks at homogeneous ordinary differential equations and shows it can recover the exact form of both linear and nonlinear differential equations. Two examples are presented. The first example uses data from the Bessel equation, representing a linear equation. The second example uses the data from the Blassius equation which is nonlinear. In both cases the exact form of the equation is identified from the given discrete data.


2020 ◽  
Vol 0 (0) ◽  
Author(s):  
Kusano Takaŝi ◽  
Jelena V. Manojlović

AbstractWe study the asymptotic behavior of eventually positive solutions of the second-order half-linear differential equation(p(t)\lvert x^{\prime}\rvert^{\alpha}\operatorname{sgn}x^{\prime})^{\prime}+q(% t)\lvert x\rvert^{\alpha}\operatorname{sgn}x=0,where q is a continuous function which may take both positive and negative values in any neighborhood of infinity and p is a positive continuous function satisfying one of the conditions\int_{a}^{\infty}\frac{ds}{p(s)^{1/\alpha}}=\infty\quad\text{or}\quad\int_{a}^% {\infty}\frac{ds}{p(s)^{1/\alpha}}<\infty.The asymptotic formulas for generalized regularly varying solutions are established using the Karamata theory of regular variation.


1931 ◽  
Vol 27 (4) ◽  
pp. 546-552 ◽  
Author(s):  
E. C. Bullard ◽  
P. B. Moon

A mechanical method of integrating a second-order differential equation, with any boundary conditions, is described and its applications are discussed.


Filomat ◽  
2019 ◽  
Vol 33 (13) ◽  
pp. 4013-4020
Author(s):  
Jianren Long ◽  
Sangui Zeng

We investigate the [p,q]-order of growth of solutions of the following complex linear differential equation f(k)+Ak-1(z) f(k-1) + ...+ A1(z) f? + A0(z) f = 0, where Aj(z) are analytic in C? - {z0}, z0 ? C. Some estimations of [p,q]-order of growth of solutions of the equation are obtained, which is generalization of previous results from Fettouch-Hamouda.


2013 ◽  
Vol 21 (2) ◽  
pp. 35-52
Author(s):  
Benharrat Belaïdi ◽  
Habib Habib

Abstract In this paper, we investigate the order and the hyper-order of growth of solutions of the linear differential equation where n≥2 is an integer, Aj (z) (≢0) (j = 1,2) are entire functions with max {σ A(j) : (j = 1,2} < 1, Q (z) = qmzm + ... + q1z + q0 is a nonoonstant polynomial and a1, a2 are complex numbers. Under some conditions, we prove that every solution f (z) ≢ 0 of the above equation is of infinite order and hyper-order 1.


2018 ◽  
Vol 20 (04) ◽  
pp. 1750038
Author(s):  
Andrei Minchenko ◽  
Alexey Ovchinnikov

Motivated by developing algorithms that decide hypertranscendence of solutions of extensions of the Bessel differential equation, algorithms computing the unipotent radical of a parameterized differential Galois group have been recently developed. Extensions of Bessel’s equation, such as the Lommel equation, can be viewed as homogeneous parameterized linear differential equations of the third order. In this paper, we give the first known algorithm that calculates the differential Galois group of a third-order parameterized linear differential equation.


1975 ◽  
Vol 27 (3) ◽  
pp. 508-512
Author(s):  
G. B. Gustafson ◽  
S. Sedziwy

Consider the wth order scalar ordinary differential equationwith pr ∈ C([0, ∞) → R ) . The purpose of this paper is to establish the following:DECOMPOSITION THEOREM. The solution space X of (1.1) has a direct sum Decompositionwhere M1 and M2 are subspaces of X such that(1) each solution in M1\﹛0﹜ is nonzero for sufficiently large t ﹛nono sdilatory) ;(2) each solution in M2 has infinitely many zeros ﹛oscillatory).


2017 ◽  
Vol 6 (2) ◽  
pp. 18-37 ◽  
Author(s):  
Vijaya Lakshmi V. Nadimpalli ◽  
Rajeev Wankar ◽  
Raghavendra Rao Chillarige

In this article, an innovative Genetic Algorithm is proposed to find potential patches enclosing roots of real valued function f:R→R. As roots of f can be real as well as complex, the function is reframed on to complex plane by writing it as f(z). Thus, the problem now is transformed to finding potential patches (rectangles in C) enclosing z such that f(z)=0, which is resolved into two components as real and imaginary parts. The proposed GA generates two random populations of real numbers for the real and imaginary parts in the given regions of interest and no other initial guesses are needed. This is the prominent advantage of the method in contrast to various other methods. Additionally, the proposed ‘Refinement technique' aids in the exhaustive coverage of potential patches enclosing roots and reinforces the selected potential rectangles to be narrow, resulting in significant search space reduction. The method works efficiently even when the roots are closely packed. A set of benchmark functions are presented and the results show the effectiveness and robustness of the new method.


Sign in / Sign up

Export Citation Format

Share Document