Variational Homotopy Perturbation Method for Solving Riccati Type Differential Problems

Abstract

In this paper, a Variational homotopy perturbation method is proposed to solve nonlinear Riccati differential equation. By combining the Variational Iteration Method and the Homotopy Perturbation Method, this technique possesses a fast convergence rate with high accuracy. The results reveal that the proposed method is very effective and simple.

Share and Cite:

Kashkari, B. and Saleh, S. (2017) Variational Homotopy Perturbation Method for Solving Riccati Type Differential Problems. Applied Mathematics, 8, 893-900. doi: 10.4236/am.2017.87070.

1. Introduction

The Riccati equation plays a great role in blueprint and analysis the linear and nonlinear optimal control problems. Numerical Solution of this equation has been acquired by applying Adomian’s decomposition method  , homotopy Analysis method HAM  , variational iteration method VIM  and homotopy perturbation method HPM  . HPM introduced by He  , it can solve a large class of nonlinear problems activity, accurately and easily.

The application of HPM on nonlinear problems has been implemented by scientists and engineers, because this method is to continuously deform a difficult problem under study into a simple problem easy to solve. VIM proved by Ji- Huan He  . It is simple and powerful method for solving a broad type of nonlinear Problem. It was shown that this method is operative and reliable analytic and numerical purposes. The method gives rapidly convergent successive approximation of the exact solution if such solution existed.

The nonlinear Riccati differential equation  has following form

$\begin{array}{l}L\left(x\right)=A\left(x\right){u}^{2}\left(x\right)+B\left(x\right)u\left(x\right)+C\left(x\right),\text{\hspace{0.17em}}\text{\hspace{0.17em}}0\le x\le X\\ u\left(0\right)=\alpha \end{array}$ (1)

where $L=\frac{\text{d}}{\text{d}x}$ or $\frac{{\text{d}}^{2}}{\text{d}{x}^{2}}$ , $A\left(x\right),B\left(x\right)$ and $C\left(x\right)$ are continuous functions

and $\alpha$ is an arbitrary constant.

We organize the following paper as follows. In Section 2, we present the VIM, while in Section 3, we present the HPM. In Section 4, we apply the VHPM to solve quadratic Riccati equation. Moreover, we find solutions of some examples by VHPM in Section 5.

The results reveal that the proposed method is very effective and simple. We end this paper by conclusion that reveal that these methods are very effective and convenient for solving nonlinear Riccati equations.

2. Variational Iteration Method

To illustrate the basic concepts of the VIM we consider the following differential equation  

$L\left(u\right)+N\left(u\right)=f\left(x\right)$ (2)

where L is a linear operator, N is a nonlinear operator, and $f\left(x\right)$ is an inhomogeneous term. Then, we can construct the correct functional as follows

${u}_{n+1}={u}_{n}+{\int }_{0}^{x}\lambda \left(\xi \right)\left[L\left({u}_{n}\right)+N\left({\stackrel{˜}{u}}_{n}\right)-f\left(\xi \right)\right]\text{d}\xi$ (3)

where $\lambda$ is a general Lagrangian multiplier defined as 

$\lambda \left(x,t\right)=\frac{{\left(-1\right)}^{m}}{\left(m-1\right)!}{\left(x-t\right)}^{m-1},\text{\hspace{0.17em}}\text{\hspace{0.17em}}m\ge 1$ (4)

And ${\stackrel{˜}{u}}_{n}$ are restricted variation which means $\delta {\stackrel{˜}{u}}_{n}=0$ . Consequently, the solution $u=\underset{n\to \infty }{\mathrm{lim}}{u}_{n}$ .

3. Homotopy Perturbation Method

To explain this method, we construct the following function 

$A\left(u\right)-f\left(x\right)=0,\text{\hspace{0.17em}}\text{\hspace{0.17em}}x\in \Omega$ (5)

With boundary condition

$B\left(u,\partial u/\partial n\right)=0,\text{\hspace{0.17em}}\text{\hspace{0.17em}}x\in \Gamma$ (6)

where $A$ is a general differential operator, $B$ is a boundary operator, $f\left(x\right)$ is a known analytical function. The operator $A$ can be decomposed into two operators $L$ and $N$ , where $L$ is a linear operator and $N$ is a nonlinear operator.

By using the homotopy technique, we construct a homotopy $u\left(x,p\right):\Omega ×\left[0,1\right]\to ℝ$ which are satisfies

$H\left(u,p\right)=\left(1-p\right)\left[L\left(u\right)-L\left({u}_{0}\right)\right]+p\left[L\left(u\right)+N\left(u\right)-f\left(x\right)\right]=0$ (7)

or

$H\left(u,p\right)=L\left(u\right)-L\left({u}_{0}\right)+p\left[L\left({u}_{0}\right)+N\left(u\right)-f\left(x\right)\right]=0$ (8)

where $p\in \left[0,1\right]$ is an embedding parameter. ${u}_{0}$ is an initial approximation of solution of equation.

$L\left(u\right)+N\left(u\right)-f\left(x\right)=0$ (9)

We have

$H\left(u,0\right)=L\left(u\right)+L\left({u}_{0}\right)=0,\text{\hspace{0.17em}}H\left(u,1\right)=A\left(u\right)-f=0$ (10)

The solution can be written as a power series in $p$ , $u\left(x\right)=\underset{p\to 1}{\mathrm{lim}}u={u}_{0}+p{u}_{1}+{p}^{2}{u}_{2}+\cdots$ .

4. Variational Homotopy Perturbation Method

In this section, we apply the VHPM to Riccati Equation (1), we start this method by applying HPM in Equation (8) on Equation (1), we get

$L{u}_{n}\left(x\right)-L{u}_{0}\left(x\right)+p\left[L{u}_{0}\left(x\right)-A\left(x\right){u}_{n}^{2}-B\left(x\right){u}_{n}-C\left(x\right)\right]=0$ (11)

Now we using correction functional in Equation (5) to get

${u}_{n+1}={u}_{n}+{\int }_{0}^{x}\lambda \left(\xi \right)\left[L{u}_{n}-L{u}_{0}+p\left(L{u}_{0}-A\left(\xi \right){u}_{n}^{2}-B\left(\xi \right){u}_{n}-C\left(\xi \right)\right)\right]\text{d}\xi$ (12)

We can obtain

${u}_{n+1}={u}_{0}+p{\int }_{0}^{x}\lambda \left(\xi \right)\left[L{u}_{0}-A\left(\xi \right){u}_{{}^{n}}^{2}-B\left(\xi \right){u}_{n}-C\left(\xi \right)\right]\text{d}\xi$ (13)

Now we can rewrite Equation (13) in the form

$\underset{n=0}{\overset{\infty }{\sum }}{p}^{n}{u}_{n}={u}_{0}+p\underset{0}{\overset{x}{\int }}\left[\lambda \left(\xi \right)\left(L{u}_{0}-N\underset{n=0}{\overset{\infty }{\sum }}{p}^{n}{\stackrel{˜}{u}}_{n}-C\left(\xi \right)\right)\right]\text{d}\xi$ (14)

As we see, the procedure is formulated by the coupling of VIM and HPM    . A comparison of like powers of $p$ give solutions of various orders.

5. Numerical Examples

5.1. Example

Consider the following classical Riccati differential equation

${u}^{\prime }\left(x\right)=-{u}^{2}+2u+1$ (15)

With initial condition $u\left(0\right)=0$ .

For the above differential equation, the exact solution  is previously known to be

$u\left(x\right)=1+\sqrt{2}\mathrm{tanh}\left(\sqrt{2}x+\frac{1}{2}\mathrm{log}\left(\frac{\sqrt{2}-1}{\sqrt{2}+1}\right)\right)$ (16)

The Taylor expansion of $u\left(x\right)$ about $x=0$ gives

$u\left(x\right)=x+{x}^{2}+\frac{1}{3}{x}^{3}-\frac{1}{3}{x}^{4}-\frac{7}{15}{x}^{5}-\frac{7}{45}{x}^{6}+\frac{53}{315}{x}^{7}+\cdots$ (17)

Suppose that the initial approximation is ${u}_{0}=x$ .

To solve Equation (15), by the VHPM we substitution it in Equation (14)

$\underset{n=0}{\overset{\infty }{\sum }}{p}^{n}{u}_{n}={u}_{0}-p\underset{0}{\overset{x}{\int }}\left[-2\underset{n=0}{\overset{\infty }{\sum }}{p}^{n}{u}_{n}+{\left(\underset{n=0}{\overset{\infty }{\sum }}{p}^{n}{u}_{n}\right)}^{2}\right]\text{d}\xi$ (18)

Here $\lambda =-1$ .

By comparing the coefficient of like powers of $p$ , we have

$\begin{array}{l}{p}^{\left(0\right)}:{u}_{0}=x\\ {p}^{\left(1\right)}:{u}_{1}=-{\int }_{0}^{x}\left(-2{u}_{0}+{u}_{0}^{2}\right)\text{d}\xi ={x}^{2}-\frac{1}{3}{x}^{3}\\ {p}^{\left(2\right)}:{u}_{2}=-{\int }_{0}^{x}\left(-2{u}_{1}+2{u}_{0}{u}_{1}\right)\text{d}\xi =\frac{2}{3}{x}^{3}-\frac{2}{3}{x}^{4}+\frac{2}{15}{x}^{5}\\ {p}^{\left(3\right)}:{u}_{3}=-{\int }_{0}^{x}\left(-2{u}_{2}+2{u}_{0}{u}_{2}+{u}_{1}^{2}\right)\text{d}\xi =\frac{1}{3}{x}^{4}-\frac{11}{15}{x}^{5}+\frac{17}{45}{x}^{6}-\frac{17}{315}{x}^{7}\\ \text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}⋮\end{array}$ (19)

The other components of the VHPM can be determined in similar way. Finally, the approximate solution of Equation (15) is $u={u}_{0}+{u}_{1}+{u}_{2}+{u}_{3}+\cdots$ . Which converge to the exact solution in Equation (16).

5.2. Example

Consider the following quadratic Riccati differential equation

${u}^{\prime }\left(x\right)=-2{\text{e}}^{x}{u}^{2}+2{\text{e}}^{2x}u+{\text{e}}^{x}-{\text{e}}^{3x}$ (20)

With initial condition $u\left(0\right)=1$ .

For the above differential equation, the exact solution  is previously known to be

$u\left(x\right)={\text{e}}^{x}$ (21)

The Taylor expansion of $u\left(x\right)$ about $x=0$ gives

$u\left(x\right)=1+x+\frac{1}{2}{x}^{2}+\frac{1}{6}{x}^{3}+\frac{1}{24}{x}^{4}+\frac{1}{120}{x}^{5}+\frac{1}{720}{x}^{6}+\frac{1}{5040}{x}^{7}+\cdots$ (22)

Suppose that the initial approximation is ${u}_{0}=x+1$ .

To solve Equation (20), by the VHPM we substitution it in Equation (14), then we get

$\underset{n=0}{\overset{\infty }{\sum }}{p}^{n}{u}_{n}={u}_{0}-p\underset{0}{\overset{x}{\int }}\left[1+\underset{n=0}{\overset{\infty }{\sum }}\frac{{\left(2\xi \right)}^{n}}{n!}\left(\underset{n=0}{\overset{\infty }{\sum }}{p}^{n}{u}_{n}\right)-2\underset{n=0}{\overset{\infty }{\sum }}\frac{{\left(2\xi \right)}^{n}}{n!}{\left(\underset{n=0}{\overset{\infty }{\sum }}{p}^{n}{u}_{n}\right)}^{2}+\underset{n=0}{\overset{\infty }{\sum }}\frac{{\left(3\xi \right)}^{n}}{n!}-\underset{n=0}{\overset{\infty }{\sum }}\frac{{\xi }^{n}}{n!}\right]\text{d}\xi$ (23)

Here $\lambda =-1$ .

By comparing the coefficient of like powers of $p$ , we have

$\begin{array}{l}{p}^{\left(0\right)}:{u}_{0}=x+1\\ {p}^{\left(1\right)}:{u}_{1}=\frac{1}{2}{x}^{2}+\frac{1}{6}{x}^{3}+\frac{1}{24}{x}^{4}-\frac{1}{24}{x}^{5}-\frac{49}{270}{x}^{6}-\frac{37}{270}{x}^{7}-\frac{1091}{40320}{x}^{8}\\ \text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}+\frac{247}{45360}{x}^{9}-\frac{1}{50400}{x}^{10}\\ {p}^{\left(2\right)}:{u}_{2}=\frac{1}{10}{x}^{5}+\frac{5}{36}{x}^{6}+\frac{13}{126}{x}^{7}+\frac{23}{480}{x}^{8}+\frac{19}{2592}{x}^{9}-\frac{3389}{302400}{x}^{10}-\cdots \\ {p}^{\left(3\right)}:{u}_{3}=-\frac{1}{20}{x}^{5}-\frac{5}{72}{x}^{6}-\frac{13}{252}{x}^{7}-\frac{1}{120}{x}^{8}+\frac{59}{1728}{x}^{9}+\frac{32119}{604800}{x}^{10}+\cdots \\ \text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}⋮\end{array}$ (24)

The other components of the VHPM can be determined in similar way. Finally, the approximate solution of Equation (20) is $u={u}_{0}+{u}_{1}+{u}_{2}+{u}_{3}+\cdots$ . Which converge to the exact solution in Equation (21).

5.3. Example

Consider the Riccati Type Painleve’s First Transcendent equation 

${u}^{″}\left(x\right)=6{u}^{2}+\mu x,\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\mu =1$ (25)

With initial conditions $u\left(0\right)=1,{u}^{\prime }\left(0\right)=0$ .

The above differential equation without known exactly solutions and we suppose that the initial approximation is ${u}_{0}=1+{x}^{2}$ .

To solve Equation (25) by the VHPM, substitution it in Equation (14), then we get

$\underset{n=0}{\overset{\infty }{\sum }}{p}^{n}{u}_{n}={u}_{0}+p\underset{0}{\overset{x}{\int }}\lambda \left(\xi \right)\left[2-6{\left(\underset{n=0}{\overset{\infty }{\sum }}{p}^{n}{u}_{n}\right)}^{2}-\xi \right]\text{d}\xi$ (26)

In this example $\lambda =\left(\xi -x\right)$ .

By comparing the coefficient of like powers of $p$ , we have

$\begin{array}{l}{p}^{\left(0\right)}:{u}_{0}=1+{x}^{2}\\ {p}^{\left(1\right)}:{u}_{1}=2{x}^{2}+\frac{1}{6}{x}^{3}+{x}^{4}+\frac{1}{5}{x}^{6}\\ {p}^{\left(2\right)}:{u}_{2}=2{x}^{4}+\frac{1}{10}{x}^{5}+\frac{6}{5}{x}^{6}+\frac{1}{21}{x}^{7}+\frac{9}{35}{x}^{8}+\frac{2}{75}{x}^{10}\\ {p}^{\left(3\right)}:{u}_{3}=\frac{8}{5}{x}^{6}+\frac{13}{105}{x}^{7}+\frac{1877}{1680}{x}^{8}+\frac{11}{210}{x}^{9}+\frac{11}{35}{x}^{10}+\frac{17}{1925}{x}^{11}\\ \text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}+\frac{254}{5775}{x}^{12}+\frac{{x}^{14}}{325}\\ \text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}⋮\end{array}$ (27)

The other components of the VHPM can be determined in similar way. Finally, the approximate solution of Equation (25) is $u\approx {u}_{0}+{u}_{1}+{u}_{2}+{u}_{3}+\cdots$ . In Table 1 we present the comparison between the approximate solution founded by VHPM with Truncated Taylor series(TTS)  , and Rational approximation (RA)  .

5.4. Example

Consider the Riccati Type Painleve’s Second Transcendent equation 

${u}^{″}\left(x\right)=2{u}^{3}+xu+\mu ,\text{\hspace{0.17em}}\text{\hspace{0.17em}}\mu =1$ (28)

With initial conditions $u\left(0\right)=1,{u}^{\prime }\left(0\right)=0$ .

The above differential equation without known exactly solutions and we suppose that the initial approximation is ${u}_{0}=1+{x}^{2}$ .

To solve Equation (28) by the VHPM, substitution it in Equation (14), then we get

$\underset{n=0}{\overset{\infty }{\sum }}{p}^{n}{u}_{n}={u}_{0}+p\underset{0}{\overset{x}{\int }}\lambda \left(\xi \right)\left[1-2{\left(\underset{n=0}{\overset{\infty }{\sum }}{p}^{n}{u}_{n}\right)}^{3}-\xi \underset{n=0}{\overset{\infty }{\sum }}{p}^{n}{u}_{n}\right]\text{d}\xi$ (29)

Table 1. Comparison between the approximate solution $u=\underset{i=0}{\overset{20}{\sum }}{u}_{i}$ with TTS and RA.

And we have $\lambda =\left(\xi -x\right)$ .

By comparing the coefficient of like powers of $p$ , we have

$\begin{array}{l}{p}^{\left(0\right)}:{u}_{0}=1+{x}^{2}\\ {p}^{\left(1\right)}:{u}_{1}=\frac{1}{2}{x}^{2}+\frac{1}{6}{x}^{3}+\frac{1}{2}{x}^{4}+\frac{1}{20}{x}^{5}+\frac{1}{5}{x}^{6}+\frac{1}{28}{x}^{8}\\ {p}^{\left(2\right)}:{u}_{2}=\frac{1}{4}{x}^{4}+\frac{3}{40}{x}^{5}+\frac{11}{36}{x}^{6}+\frac{1}{51}{x}^{7}+\frac{41}{224}{x}^{8}+\frac{1}{224}{x}^{9}\\ \text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}+\frac{131}{2100}{x}^{10}+\frac{47}{15400}{x}^{11}+\frac{19}{1540}{x}^{12}+\frac{3}{2548}{x}^{14}\\ {p}^{\left(3\right)}:{u}_{3}=\frac{1}{10}{x}^{6}+\frac{17}{420}{x}^{7}+\frac{383}{2240}{x}^{8}+\frac{703}{12960}{x}^{9}+\frac{6887}{50400}{x}^{10}\\ \text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}+\frac{4013}{123200}{x}^{11}+\frac{1039}{15840}{x}^{12}+\frac{2759}{257400}{x}^{13}+\frac{117269}{5605600}{x}^{14}\\ \text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}+\frac{151}{77000}{x}^{15}+\frac{74933}{16816800}{x}^{16}+\frac{28671}{190590400}{x}^{17}+\frac{5653}{9529520}{x}^{18}\\ \text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}+\frac{15}{387296}{x}^{20}\\ \text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}⋮\end{array}$ (30)

The other components of the VHPM can be determined in similar way. Finally, the approximate solution of Equation (28) is $u\approx {u}_{0}+{u}_{1}+{u}_{2}+{u}_{3}+\cdots$ . In Table 2 we present the comparison between the approximate solution founded by VHPM with TTS  , and RA  .

6. Conclusion

In this paper, we studied the solution of nonlinear Riccati differential equation. We have applied a recently introduced technique called the VHPM to solve this nonlinear differential equation. This method is more efficient and simpler. The results of the method exhibit excellent agreement with the exact solution. The

Table 2. Comparison between the approximate solution $u=\underset{i=0}{\overset{20}{\sum }}{u}_{i}$ with TTS and RA.

comparison between the numerical results with TTS and RA in Problems 3 and 4 of validates the accuracy of the VHPM method for problems without known exactly solutions that have been advanced for solving Riccati equation shows that the new technique is reliable and powerful.

Conflicts of Interest

The authors declare no conflicts of interest.     customer@scirp.org +86 18163351462(WhatsApp) 1655362766  Paper Publishing WeChat 