The Bezier Control Points Method for Solving Delay Differential Equation

In this paper, Bezier surface form is used to find the approximate solution of delay differential equations (DDE’s). By using a recurrence relation and the traditional least square minimization method, the best control points of residual function can be found where those control points determine the approximate solution of DDE. Some examples are given to show efficiency of the proposed method.


Introduction
Delay differential equations are type of differential equations where the time derivatives at the current time depend on the solution, and possibly its derivatives, at previous times.A class of such equations, which involve derivatives with delays as well as the solution itself has been called neutral DDEs over the past century (see [1,2]).
The basic theory concerning the stable factors and works on fundamental theory, e.g., existence and uniqueness of solutions, was presented in [1,2].Since then, DDE have been extensively studied in recent decades and a great number of monographs have been published including significant works on dynamics of DDEs by Hale and Lunel [3], on stability by Niculescu [4], and so on.The interest in study of DDEs is caused by the fact that many processes have time-delays and have been models for better representations by systems of DDEs in science, engineering, economics, etc.Such systems, however, are still not feasible to actively analyze and con-trol precisely, thus, the study of systems of DDEs has actively been conducted over the recent decades (see [1,2]).
In this paper, we show a novel strategy by using the Bezier curves to find the approximate solution for delay differential equations by Bezier curves.Other numerical methods for DDEs are available in (see [5][6][7][8]).In Section 2 delay differential equations will be introduced.Example of Time-Delay System will be stated in Section 3. In Section 4 delay differential equations with proportional delay will be introduced.Bezier curves and degree eleva-tion will be stated in Sections 5 and 6 respectively.In Section 7 solution of delay differential equation using Bezier control points presented and aforementioned method will be implemented on it.In Section 8, solved numerical examples, showed the efficiency and reliability of the method.Finally, Section 9 will give a conclusion briefly.

Delay Differential Equations
Most delay differential equations that arise in population dynamics and epidemiology model intrinsically nonnegative quantities.Therefore it is important to establish that nonnegative initial data give rise to nonnegative solutions.Consider the following with a single delay h > 0. Assume that f t u y and

 
, , u t u y are continuous on R 3 .Let s f R  be given and let and satisfying (2.1) on systems is quite natural since there must be finite period of time following a decision for its effects to appear.In one model [9] of aggregate economy, we let   Y t be the income which can split into consumption   C t

 
, investment I t and autonomous expenditure. Thus where c is a consumption coefficient.From (3.1) we get It is assumed that there is finite interval of time between ordering and delivery of capital equipment following a decision to invest  .

 
In terms of the stock of capital assets we have Economic rationale implies that is determined by the rate of saving (proportional to ) and by the capital stock .This means that where   0 ,   and ε is a trend factor.Combining (3.4) and (3.5), we obtain: By (3.3) and (3.7), we arrive at Finally, it follows from (3.5), (3.6) and (3.8) that which expresses the formation of the rate of delivery of the new equipment.This is a typical functional differential equation (FDE) of retarded type.

Delay Differential Equations with Proportional Delay
In this paper, approximate analytical solutions with high accuracy can be obtained by carrying out in the Bezier control points method.Consider the following neutral functional-differential equation with proportional delays (see [10][11][12]), The existence and the uniqueness of the analytic solution of the multi-pantograph equation are proved in [13], the Dirichlet series solution is constructed, and the sufficient condition of the asymptotic stability for the analytic solution is obtained.It is proved that the θ-methods with a variable stepsize are asymptotically stable if Some numerical examples are given to show the properties of the θ-methods.
In order to apply the Bezier control points method, we rewrite Equation (4.1) as Neutral functional-differential equations with proportional delays represent a particular class of delay differential equation.Such functional-differential equations play an important role in the mathematical modeling of real world phenomena [14].Obviously, most of these equations cannot be solved exactly.It is therefore necessary to design efficient numerical methods to approximate their solutions.Ishiwata et al. used the rational approximation method [15] and the collocation method [16] to compute numerical solutions of delay differential equations with proportional delays.Hu et al. [17] applied linear multistep methods to compute numerical solutions for neutral delay differential equations.Wang et al. obtained approximate solutions for neutral delay differential equations by continuous Runge-Kutta methods [18] and oneleg θ-methods [13,19].

Bezier Curves
A Bezier curve of degree n can be defined as follows (see [19]): where , are the Bern- . The Bezier coefficient is called the control point (see Figure 1).In particular , , P P an explicit Bezier curve by (see [20,21]).
   , t C t 

Degree Elevation
Suppose we were designing with Bezier curve as described, and use a Bezier polygon of degree n to approximate the desired given shape.Suppose the degree polygon dose not feat neatly the desired shape.
One way to proceed in such a situation is to increase the flexibility of the polygon by adding another vertex (control point) to it.As a first step, one might want to add another vertex, yet leave the desired curve of the shape unchanged, this corresponds to raising the degree of the Bezier curve by one (see Figure 2).Therefore, we are looking for a curve with control vertices 0 1 n  , , P P  that describes the same curve of the shape as the original polygon (see [21][22][23][24][25] for more details).
0 n We rewrite our given Bezier curve as The upper index of the first sum may be extended to n + 1, since the corresponding term is zero.The summation indices of the second sum may be shifted to index 1 and n + 1, but one may choose the lower index zero since only a zero term is added.Thus we have Combining both sums and computing coefficients  yields: where i is the control point of the Bezier curve  C t when it is elevated to degree n + 1.Now, the new control polygon consists of n + 2 control points.

Solution of Delay Differential Equation Using Bezier Control Points
Consider the following boundary value problem where L is differential operator with proportional delay,  is also a polynomial in t, and (k = 0, 1, ••• , m) [26].
We propose to represent the approximate solution of (7.1)   u t in Bezier form.The choice of the Bezier form rather than the B-spline form is due to the fact that Copyright © 2012 SciRes.

ICA
the Bezier form is easier to symbolically carry out the operations of multiplication, comparison and degree elevation than B-spline form.We choose the sum of squares of the Bezier control points of the residual to be the measure quantity.Minimizing this quantity gives the approximate solution.So, the obvious spotlight is in the following, if the minimizing of the quantity is zero, so the residual function is zero, which implies that the solution is the exact solution.We call this approach the control-point-based method.The detailed steps of the method are as follows (see [24]):  Step 1. Choose a degree n and symbolically express the solution   u t in the degree   where the control points are to be determined. Step 2. Substituting the approximate solution      u u t into the differential Equation (7.1), we gain the residual function So the residual function can be expressed in Bezier form as well, where the control points 0 1 k are linear functions in the unknowns i .These functions are derived using the operations of multiplication, degree elevation and differentiation for Bezier form.
Then F is also a function of , , , .
by some optimization techniques, such as Lagrange multipliers method, we can be used to solve (7.4). Step 5. Substituting the minimum solution back into (7.2) arrives at the approximate solution to the differential equation.

Numerical Examples
In this part, we used the mentioned control-point-based method on Bezier control points to solve DDE's and system of DDE's.Example 8.1.As a practical example, we consider Evens and Raslan [6] the following pantograph delay equation:


The exact solution is . Now we try to find a degree two approximate solution.Let Substituting it into the above delay differential equation Thus the approximate solution is In Figure 3 compare approximated and exact value of   u t .Figure 4 shows the residual function.Example 8.2.Consider the previous example with degree raising in Bezier control points.
Let In Figure 5 compare approximated and exact value of   u t .Figure 6 shows the residual function.Example 8.3.Consider the following second order linear DDE (see [5]):   .Figure 8 shows the residual function.Example 8.4.Consider the following second order linear DDE (see [10]):  where the exact solution is   4 .

Conclusion
In this paper, we use the control-point-based method to solve delay differential equations.In this method, firstly, Copyright © 2012 SciRes.ICA the rough solution is expressed in Bezier form, then the residual function is minimized to find the best approximate solution.Some examples are given to verify the reliability and efficiency of the proposed method.
Bezier curve.The control polygon of a Bezier curve comprise of the line segments

Figure 1 .
Figure 1.A degree three Bezier curve and its control polygon.

Figure 5 .
Figure 5. Approximate and exact solution of u(t) for Example 8.2.
the exact solution is .

Figure 7 .
Figure 7. Approximate and exact solution u(t) for Example 8.3.

Figure 9 , 8 . 5 .
Figure 9, compare the exact and approximate solution of.Figure10shows the residual function.Example 8.5.In this example the following first order linear DDE's is considered (see[13]):     u t u t u t   

Figure 9 .
Figure 9. Approximate and exact solution u(t) for Example 8.4.

Figure 11
shows the approximate value of .