Mathematical Problems in Engineering

Volume 2012 (2012), Article ID 259494, 15 pages

doi:10.1155/2012/259494

Research Article

Approximate Analytic and Numerical Solutions to Lane-Emden Equation via Fuzzy Modeling Method

De-Gang Wang,¹ Wen-Yan Song,² Peng Shi,3,4 and Hamid Reza Karimi⁵

¹School of Control Science and Engineering, Dalian University of Technology, Dalian 116024, China

²Department of Quantitative Economics, Dongbei University of Finance and Economics, Dalian 116025, China

³School of Electrical and Electronic Engineering, The University of Adelaide, Adelaide, SA 5005, Australia

⁴School of Engineering and Science, Victoria University, Melbourne, VIC 8001, Australia

⁵Department of Engineering, Faculty of Technology and Science, University of Agder, 4898 Grimstad, Norway

Received 19 August 2012; Accepted 21 October 2012

Academic Editor: Mohammed Chadli

Copyright © 2012 De-Gang Wang et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

A novel algorithm, called variable weight fuzzy marginal linearization (VWFML) method, is proposed. This method can supply approximate analytic and numerical solutions to Lane-Emden equations. And it is easy to be implemented and extended for solving other nonlinear differential equations. Numerical examples are included to demonstrate the validity and applicability of the developed technique.

1. Introduction

Lane-Emden equations are used to describe singular initial value problems (IVPs) relating to second-order ordinary differential equations which have been used to model several phenomena in mathematical physics, thermodynamics, fluid mechanics, and astrophysics such as the theory of stellar structure, the thermal behavior of a spherical cloud of gas, isothermal gas spheres, and theory of thermionic currents.

Lane-Emden equations, first introduced by Jonathan Homer Lane in 1870 and further explored in detail by Emden, have the following form:

subject to the conditions

where and are constants and is a continuous real valued function.

Since Lane-Emden type equations have significant applications in many scientific fields, various forms of have been investigated in many research works. Among them, many attentions have been carried on the generalized Lane-Emden type equations,

subject to condition (1.2), where is a continuous real valued function, and .

Many different methods have been used to obtain solutions for the generalized Lane-Emden equations. Wazwaz [1, 2] got approximate solutions by using the Adomian decomposition method (ADM) and obtained the analytic solutions of some equations. But it may be an intricate problem to calculate the so-called Adomian polynomials involved in ADM sometimes. He [3, 4] developed a more convenient analytical technique, called the homotopy perturbation method (HPM). Chowdhury and Hashim [5] and Yildirim and Öziş [6] gave the solutions for a class of singular second-order IVPs of Lane-Emden type by using HPM. Sajid et al. [7] pointed out that HPM is a special case of the homotopy analysis method (HAM) and that it is valid only for weakly nonlinear problems. Liao [8] and Van Gorder and Vajravelu [9] used HAM to increase the radius of convergence of series solutions for Lane-Emden equations. Recently, Yiğider et al. studied a numerical method for solving Lane-Emden type equations by Padé approximation in [10]. Generally, when all the above cited analytical approaches are used to solve Lane-Emden equation, a truncated power series solution of the true solution is obtained. By the methods such as HPM and HAM, a series of newly achievements on the analytical solving for some nonlinear differential equations have been proposed recently. By HAM, Ziabakhsh et al. [11, 12] studied the natural convection of a non-Newtonian fluid between two infinite parallel vertical flat plates and the effects of the non-Newtonian nature of fluid on the heat transfer. Jalaal et al. [13, 14] investigated the settling behavior of solid particles using HPM, which show the capability and effectiveness of the method and exhibit new application of it further.

Besides, many soft computing technologies are developed to deal with all types of models for dynamic systems [15–18]. It is important to note that fuzzy modeling technology can transfer data information into a mathematical model which can approximate the original system with high accuracy [19–24]. Li et al. [25] used fuzzy modeling method to approximate the solutions of a class of autonomous differential equation. In order to obtain the analytical solution of the fuzzy system, Li et al. [26] introduced fuzzy marginal linearization method. Further, Li et al. [27] proposed fuzzy inference modeling (FIM) method to approximate the time-variant system. Wang et al. [28] proposed a dynamic fuzzy inference modeling (DFIM) method and proved that fuzzy system generalized by this method was universal approximators to the solutions of some nonautonomous differential equations.

However, above fuzzy modeling technology could not be used to solve the differential equation (1.3). And when and in (1.3) are unknown and only input-output data of them are obtained, how to obtain the corresponding solution is an interesting question. Motivated by this fact, the aim of this paper is to propose a novel fuzzy modeling method to solve the Lane-Emden equation. This paper is organized as follows. In Section 2, we introduce some preliminary knowledge. In Section 3, we propose a novel fuzzy modeling technology and use it to obtain the approximate analytical solutions of the Lane-Emden equation. Some examples are used to illustrate the validity of the proposed method in Section 4. Finally, conclusions are presented in Section 5.

2. Preliminaries and Basic Ideas of FIM

Firstly, we introduce some basic concepts which will be used in sequel.

Definition 2.1. A fuzzy set of is a function from the reference set to the unit interval , that is, , .

Definition 2.2. Let be a group of normal fuzzy sets of , where is the peak point of , that is, . If satisfies the condition: (for all ) and for all , then is called a fuzzy partition of .

Definition 2.3. The mappings from to are variable weights if the following conditions hold

(a)for any , ;(b), .

Example 2.4. Let . The following mappings are variable weights on : , .

In the following, we will introduce FIM method. Consider a second order ordinary differential equation

where is the independent variable, is the unknown function, , , , and are the universes of , , , and , respectively, and is a real-valued continuous function, which explicitly contains the independent variable . Obviously (1.3) is the special case of (2.1). In this paper, we assume that , , , and are real number intervals, that is, , , , and .

Let    be a group of known data from (2.1), which satisfy the following two conditions:

(1),  ,  ;(2),  .

Then, we use above data information to construct fuzzy rule base. For each index , we, respectively, take , , and as the peak points of fuzzy sets , , and , such that is a fuzzy partition of , is a fuzzy partition of , and is a fuzzy partition of . Similarly, we take as peak point of fuzzy set , such that is a fuzzy partition of . In this way, fuzzy rules based on data information are represented as follows:

For fuzzy rules (2.2), fuzzy system generalized by FIM method can be expressed by

where the characteristic function is defined as

Remark 2.5. It can be seen that fuzzy system (2.3) is a nonlinear differential equation with variable coefficients. When and are, respectively, chosen as triangular membership functions, in each local region , (2.3) is changed into an autonomous differential equation with constant coefficients. This fact means that (2.3) is a two-order differential equation with piecewise constant coefficients. In [27], it is proved that solutions of (2.3) can approximate the numerical solutions of some nonautonomous differential equations with high accuracy.

3. Approximate Analytic Solutions of the Lane-Emden Equation Based on Variable Weight Fuzzy Marginal Linearization Method

From (2.3), we find that in each local region it is still a nonlinear differential equation. Hence, it is difficult for us to get the corresponding analytical solution for it. In this section, we propose a novel fuzzy modeling method, called variable weight fuzzy marginal linearization (VMFML) method, and utilize this technology to obtain the approximate analytical solution for Lane-Emden equation. For simplicity, we introduce some denotations.

Let be a partition of , where , and . Similarly, is another partition of , where . Further, we, respectively, divide and into pieces and let

where ; ; , , ; . The characteristic functions on and are, respectively, denoted as and , that is,

In the following, we will introduce the basic idea of VMFML method.

For any , without loss of generality, we assume that . For any and , by fuzzy marginal linearization technology, we take as triangular membership function and as rectangle-shaped membership function, then fuzzy system (2.3) can be changed into

Similarly, when is chosen as rectangle-shaped membership function and is chosen as triangular membership function, fuzzy system (2.3) can be changed into

Furthermore, we take sum of the right side of expressions (3.3) and (3.4) and subtract a constant, the corresponding fuzzy system in the local region is represented as

where

By characteristic function, we can obtain fuzzy system on as follows:

Remark 3.1. It is easy to see that (3.9) is a piecewise linear differential equation. In each local region , (3.9) transfers into a linear differential equation and the corresponding coefficients can be computed by input-output data of the original system. By expressions (3.5)–(3.8), it is easy to prove that the right-hand side is an interpolation function of , that is, , ; ; . By numerical analysis theory, we know that when we obtain enough data information on (2.1), can approximate with the specified accuracy. This fact means that (3.6) can be used to describe the nonlinear differential equation (2.1).

Then, we will solve the initial value problem of (3.9).

Given an initial value problem as

subject to the conditions and .

First, we can determine the local region which the initial vector locates in. Suppose that . By (3.9) the corresponding piecewise equation can be written as

and the corresponding coefficients can be computed by expressions (3.6)–(3.8). Since it is a linear ode, we can get the corresponding analytical solution and denote it as . Then, we take as the initial vector and search the next region which it moves in. And we can solve the corresponding piecewise linear equation with the initial value . In this way, by transferring initial value piece by piece and solving the corresponding piecewise linear equation, we can obtain analytical solution of (3.9).

On the other hand, if we take    as the local region of universe , using fuzzy marginal linearization technology, we can obtain another fuzzy system as follows:

where

Obviously, both (3.9) and (3.12) can describe (2.1) and the analytical solutions of them can be obtained. The differences between them are that the coefficients of (3.9) in the local region are computed by data , , ,.

On the other hand, the corresponding coefficients of (3.12) in the local region are deduced by data , , , and . For a given initial value and , by transferring initial value technology, we can solve the corresponding analytical solutions of (3.9) and (3.12) denote them as and , respectively. In order to describe (2.1) better, we take variable weighted sum of and , and denote

as the approximation analytical solution of (2.1), where and .

In this way, we use variable weight fuzzy marginal linearization (VWFML) technology to obtain the approximation analytical solution for (2.1).

Next, we also take (2.1) as example to summarize basic processes of VWFLM method.

Step 1. Determine the universes of , , and and denote them as , , and .

Step 2. Divide the universes and let , , be the partition points of , , and , where .

Step 3. By (2.1), compute the corresponding output data .

Step 4. Using expressions (3.6)–(3.8) and data information   , deduce the coefficients of piecewise equation and construct the corresponding fuzzy system (3.9).

Step 5. For given initial value, by transferring initial value technology, solve above fuzzy system and denote the solution as .

Step 6. By expression (3.13) and data information   , we can also deduce the coefficients of piecewise equation and construct another fuzzy system (3.12).

Step 7. Similarly, by transferring initial value technology, we can solve this fuzzy system (3.12) and denote the solution as .

Step 8. We take weighted sum of and , that is,

where and .

Remark 3.2. In many physical problems when differential equation models are unknown and only some data information of the investigated systems are known, by VWFML method we can still set up dynamic models and obtain corresponding approximation analytical solutions for the problems.

4. Numerical Results

In this section, we will use VWFML method to solve four Lane-Emden equations. The former two equations are of original type Lane-Emden equations and the latter two are of generalized type.

Example 4.1. Consider the linear singular initial value problem:

subject to the initial conditions and . The exact solution for this equation is .

The solution, which is generated from this Lane-Emden equation of classical astrophysics, has a proven physical foundation and can be used to calculate the radius of electron in an electromagnetic mass model.

Let , , and . Figure 1, respectively, shows curves of exact solution and approximate solutions obtained by VWFML method, DFIM method, and FIM method. Table 1 shows the coefficients of piecewise (3.9) and (3.12). In Figure 2, the absolute error curves among them are given, where “crossed line” denotes the absolute error curve between the approximate solution obtained by VWFML method and the exact solution, “” denotes the absolute error curve between the approximate solution obtained by FIM method and the exact solution, and “–*–” denotes the absolute error curve between the approximate solution obtained by DFIM method and the exact solution. Table 2 presents the errors among the exact solution and the approximate solutions obtained by them.

Table 1: Coefficients of (3.9) and (3.12).

Example 4.1 ()

1
(0, 0.25]
0
0
0
−1
−8
9.0781

2
(0.25, 0.5]
−1.0000
−8.0000
9.0781
−1.0000
−4.0000
12.3750

3
(0.5, 0.75]
−1.0000
−4.0000
12.3750
−1.0000
−2.6667
15.9844

4
(0.75, 1]
−1.0000
−2.6667
15.9844
−1.0000
−2.0000
20.0000

Table 2: Errors of different approximate solutions for Example 4.1.

Example 4.1 ()

Error of VWFML
Error of DFIM
Error of FIM

0.1
0.0067
0.0052
0.0125

0.2
0.0087
0.0059
0.0176

0.3
−0.0027
0.0028
0.0035

0.4
0.0029
−0.0035
0.0041

0.5
0.0010
−0.0345
−0.0027

0.6
0.0008
−0.0689
−0.0081

0.7
0.0072
−0.0991
−0.0070

0.8
0.0027
−0.1306
−0.0169

0.9
0.0107
−0.1678
−0.0146

1
0.0117
−0.2096
−0.0193

Figure 1: Solution curves of Example 4.1.

Figure 2: Error curves of Example 4.1.

In this simulation, the approximation results are satisfying. In fact, the number of divided local regions is not large, which means that we can get all the coefficients of (3.9) and (3.12) without a great deal of computation. Accordingly, the approximate analytic solution can be found for Example 4.1. From Figure 2 and Table 2, it is clear that for solving this equation VWFML technology is more effective than DFIM and FIM, and the corresponding amount of error data is smaller.

Example 4.2. Consider isothermal gas spheres equation:

subject to the initial conditions , .

This type of equation has been used to model the thermal behavior of a spherical cloud of gas acting under the mutual attraction of its molecules. Isothermal gaseous sphere, in which the temperature remains constant, subjects to the classical laws of thermodynamics when one seeks to determine the density and electric force of an electron gas in the neighborhood of a hot body in thermal equilibrium. It is worthy to notice that this equation is nonlinear and has no analytic solution.

From [1, 5], we know that the 6-term series solutions of this equation obtained by ADM and HPM are the same. Obviously, in Figure 3 we find that solutions obtained by VWFML method are very close to solutions of ADM method, where “” denotes the solution obtained by ADM method and “crossed line” denotes the solution obtained by VWFML method.

Figure 3: Solution curves of Example 4.2.

For the consistency of the comparison and without loss of generality, chosen the same initial value and the same independent variable values, Table 3 shows the approximations of solutions for Example 4.2, respectively, obtained by ADM method and VWFML method. Further, in Figure 4, the curves obtained by the above two methods are given. Since Lane-Emden equations are singular initial value problems, the accuracy of solutions near zero point is important, whatever the solution is gotten from approximate analytical method or from numerical method. From simulation results, we can see that if the number of local regions increased in VWFML method, the accuracy of the approximation solution can be improved evidently, especially when the independent variable is near the initial value . Besides, the solutions of Lane-Emden equations converge rapidly in a very small region , and in Figure 3 the approximate curve follows the trend appropriately. This implies that VWFML method can describe the isothermal gas spheres equation with high accuracy degree.

Table 3: Solutions of Example 4.2.

6-term ADM

VWFML solution
VWFML solution

solution

()
()

0.1
−0.0016658
0.0019471
0.0014253

0.2
−0.0066534
−0.0062487
−0.0062547

0.3
−0.014933
−0.012383
−0.014411

0.4
0.026456
−0.023321
−0.025868

0.5
−0.041154
−0.037287
−0.04051

0.6
−0.058945
−0.054402
−0.058256

0.7
−0.079728
−0.074936
−0.079011

0.8
−0.10339
−0.098012
−0.10264

0.9
−0.12981
−0.12426
−0.12904

1
−0.15886
−0.15302
−0.15806

Figure 4: Error curves of Example 4.2.

Example 4.3. Consider the linear initial value problem:

subject to the initial conditions , . The exact solution is .

Let , , and . The curves of exact solution of Example 4.3 and the corresponding approximate solution obtained by VWFML method are shown in Figure 5, where “” denotes the exact solution and “crossed line” denotes the approximate solution obtained by VWFML method.

Figure 5: Solution curves of Example 4.3.

Example 4.4. We consider the nonlinear initial value problem:

subject to the initial conditions , . The exact solution is .

Let , , and . The interval of independent variable is chosen as . The curves of exact solution of Example 4.4 and the corresponding approximate solution obtained by VWFML method are shown in Figure 6, where “” denotes the exact solution and “crossed line” denotes the approximate solution obtained by VWFML method. Furthermore, the interval of independent variable is extended as and let , , and . The simulation result is shown in Figure 7.

Figure 6: Solution curves of Example 4.4.

Figure 7: Solution curves of Example 4.4  .

Remark 4.5. From Examples 4.3 and 4.4, it can be seen that VWFML method is effective for solving generalized Lane-Emden type differential equations. Furthermore, if the solving interval is extended, VWFML method can also provide approximate solutions in the larger domains with high accuracy degree.

Remark 4.6. Comparison with some analytical solutions, solutions of VWFML methods are totally dependent on the dividing points of independent variable interval, which means that the proposed technique can be presented in a general way. In particular, when the objective equation is unknown and only some data information can be obtained, equation determined by VWFML method can be solved analytically and the corresponding solution can approximate the solution of objective equation with high accuracy.

5. Conclusion

In this paper, we apply VWFML method to obtain the approximate analytical and numerical solution for Lane-Emden type differential equation. Some numerical examples show that by relatively minor data information, solutions obtained by VWFML method can approximate the corresponding solutions of Lane-Emden type equations with high accuracy. This means that VWFML method can be utilized to solve and analyze complex nonlinear differential equations in practical application.

Acknowledgments

This work is partially supported by the National Natural Science Foundation of China (61104038, 71201019, 61074044, 61174058), the National 973 Basic Research Program of China (2009CB320602, 2012CB215202), and the 111 Project (B12018).

References

A.-M. Wazwaz, “A new algorithm for solving differential equations of Lane-Emden type,” Applied Mathematics and Computation, vol. 118, no. 2-3, pp. 287–310, 2001.

A.-M. Wazwaz, “A new method for solving singular initial value problems in the second-order ordinary differential equations,” Applied Mathematics and Computation, vol. 128, no. 1, pp. 45–57, 2002.

J. H. He, “Homotopy perturbation method: a new nonlinear analytical technique,” Applied Mathematics and Computation, vol. 135, no. 1, pp. 73–79, 2003.

J. H. He, Non-Perturbative Methods for Strongly Nonlinear Problems, Die Deutsche Bibliothek, Germany, 2006.

M. S. H. Chowdhury and I. Hashim, “Solutions of a class of singular second-order IVPs by homotopy-perturbation method,” Physics Letters A, vol. 365, no. 5-6, pp. 439–447, 2007.

A. Yildirim and T. Öziş, “Solutions of singular IVPs of Lane-Emden type by homotopy perturbation method,” Physics Letters A, vol. 369, no. 1-2, pp. 70–76, 2007.

M. Sajid, T. Hayat, and S. Asghar, “Comparison between the HAM and HPM solutions of thin film flows of non-Newtonian fluids on a moving belt,” Nonlinear Dynamics, vol. 50, no. 1-2, pp. 27–35, 2007.

S. Liao, “A new analytic algorithm of Lane-Emden type equations,” Applied Mathematics and Computation, vol. 142, no. 1, pp. 1–16, 2003.

R. A. Van Gorder and K. Vajravelu, “Analytic and numerical solutions to the Lane-Emden equation,” Physics Letters A, vol. 372, no. 39, pp. 6060–6065, 2008.

M. Yiğider, K. Tabatabaei, and E. C. Çelik, “The numerical method for solving differential equations of Lane-Emden type by Pade approximation,” Discrete Dynamics in Nature and Society, vol. 2011, Article ID 479396, 9 pages, 2011.

Z. Ziabakhsh and G. Domairry, “Analytic solution of natural convection flow of a non-Newtonian fluid between two vertical flat plates using homotopy analysis method,” Communications in Nonlinear Science and Numerical Simulation, vol. 14, no. 5, pp. 1868–1880, 2009.

Z. Ziabakhsh, G. Domairry, and H. Bararnia, “Analytical solution of non-Newtonian micropolar fluid flow with uniform suction/blowing and heat generation,” Journal of the Taiwan Institute of Chemical Engineers, vol. 40, no. 4, pp. 443–451, 2009.

M. Jalaal, D. D. Ganji, and G. Ahmadi, “Analytical investigation on acceleration motion of a vertically falling spherical particle in incompressible Newtonian media,” Advanced Powder Technology, vol. 21, no. 3, pp. 298–304, 2010.

M. Jalaal and D. D. Ganji, “On unsteady rolling motion of spheres in inclined tubes filled with incompressible Newtonian fluids,” Advanced Powder Technology, vol. 22, no. 1, pp. 58–67, 2011.

S. K. Nguang, P. Shi, and S. Ding, “Fault detection for uncertain fuzzy systems: an LMI approach,” IEEE Transactions on Fuzzy Systems, vol. 15, no. 6, pp. 1251–1262, 2007.

L. Wu, X. Su, P. Shi, and J. Qiu, “A new approach to stability analysis and stabilization of discrete-time T-S fuzzy time-varying delay systems,” IEEE Transactions on Systems, Man, and Cybernetics, Part B, vol. 41, no. 1, pp. 273–286, 2011.

M. A. Z. Raja, I. M. Qureshi, and J. A. Khan, “Swarm intelligence optimized neural networks for solving fractional differential equations,” International Journal of Innovative Computing, Information and Control, vol. 7, no. 11, pp. 6301–6318, 2011.

L. Wu, X. Su, P. Shi, and J. Qiu, “Model approximation for discrete-time state-delay systems in the TS fuzzy framework,” IEEE Transactions on Fuzzy Systems, vol. 19, no. 2, pp. 366–378, 2011.

J. Zhang, P. Shi, and Y. Xia, “Robust adaptive sliding-mode control for fuzzy systems with mismatched uncertainties,” IEEE Transactions on Fuzzy Systems, vol. 18, no. 4, pp. 700–711, 2010.

Q. Zhou, P. Shi, J. Lu, and S. Xu, “Adaptive output feedback fuzzy tracking control for a class of nonlinear systems,” IEEE Trans on Fuzzy Systems, vol. 19, no. 5, pp. 972–982, 2011.

S. Goto, T. Yamauchi, S. Katafuchi et al., “Coal type selection for power plants through fuzzy inference using coal and fly ash property,” International Journal of Innovative Computing, Information and Control, vol. 7, no. 1, pp. 29–38, 2011.

Y. Sato, T. Miki, and K. Honda, “A multi-core processor based real-time multi-modal emotion extraction system employing fuzzy inference,” International Journal of Innovative Computing, Information and Control, vol. 7, no. 8, pp. 4603–4620, 2011.

C. Y. Chang, C. H. Hsu, W. J. Wang, et al., “Operation characteristic control of direct methanol fuel cell liquid feed fuel system using adaptive Neuro-Fuzzy inference system model,” International Journal of Innovative Computing, Information and Control, vol. 8, no. 6, pp. 4177–4187, 2012.

A. N. Ahmed, A. E. Shae, O. A. Karim, and A. E. Shae, “An augmented wavelet De-Noising technique with Neuro-Fuzzy inference system for water quality prediction,” International Journal of Innovative Computing, Information and Control, vol. 8, no. 10 B, pp. 7055–7082, 2012.

H. Li, J. Wang, and Z. Miao, “Modelling on fuzzy control systems,” Science in China. Series A, vol. 45, no. 12, pp. 1506–1517, 2002.

H. Li, J. Wang, and Z. Miao, “Marginal linearization method in modeling on fuzzy control systems,” Progress in Natural Science, vol. 13, no. 7, pp. 489–496, 2003.

H. X. Li, W. Y. Song, X. H. Yuan, and Y. C. Li, “Time-varying system modeling method based on fuzzy inference,” Journal of Systems Science and Mathematical Sciences, vol. 29, no. 8, pp. 1109–1128, 2009.

D.-G. Wang, W.-Y. Song, and H.-X. Li, “Analysis and design of time-variant fuzzy systems based on dynamic fuzzy inference,” Computers & Mathematics with Applications, vol. 60, no. 3, pp. 464–489, 2010.