Table of Contents Author Guidelines Submit a Manuscript
Journal of Applied Mathematics
Volume 2014 (2014), Article ID 478241, 6 pages
Research Article

P-Fuzzy Diffusion Equation Using Rules Base

1Federal University of Piauí, 64063-040 Teresina, PI, Brazil
2Federal University of ABC, 09210-170 Santo André, SP, Brazil
3Federal Institute of Maranhão, 65635-468 Timon, MA, Brazil
4Federal University of Mato Grosso, 78075-202 Cuiabá, MT, Brazil

Received 12 February 2014; Accepted 23 April 2014; Published 6 July 2014

Academic Editor: Laurent Gosse

Copyright © 2014 Jefferson Leite 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.


We propose a fuzzy system that simulates dispersion of individuals whose movements are described by diffusion. We will use only the position of the population as an input variable for describing the process. We emphasize that the classical diffusion equation along with its analytical solution in no time was used for obtaining our solution.

1. Introduction

The differential equations and deterministic differences are a powerful tool for modeling phenomena whose state variables are subject to changes over time. However, for the deterministic modeling be efficient is necessary to have a deep well of the relationship between the variables and their variations knowledge. It is the knowledge of the phenomenon that makes it possible to choose the functions that determine the variations with respect to the state (value) of the variable. In many cases, however, this relationship between variables and variations is only partially known, which make the deterministic model is less applicable [15].

On the other hand, models of fuzzy variational equations though behaving subjectivities are also not applicable modeling phenomena with partially known relationships. This comes from the fact that these models are derived from deterministic models. Subjectivity supported by fuzzy equations refers to uncertainties as the initial states of fuzziness demographic variables and parameters of fuzziness environmental. In general, both types of fuzziness are present in equations of population dynamics [69].

The p-fuzzy systems incorporate subjective information in both variables as the variations and their relationships with the variables and are therefore a very useful tool for modeling phenomena whose behavior is partially known.

The fuzzy systems are generally the result of a generalization of the classical systems; that is, in this approach the uncertain concepts are incorporated into these systems. A central feature of fuzzy systems is that they are based on the concept of fuzzy partition information. The use of fuzzy sets allows a generalization of information that is associated with the introduction of imprecision ignoring the phenomena. In essence, the representation of information in fuzzy systems tries to imitate the process of human reasoning, considering heuristic knowledge and information across the disconnected principle [10].

In this work, we describe a diffusive process without the use of their analytical solution, using dynamical systems p-fuzzy and given a rule base. It is worth noting that the results obtained in terms of solution are very similar to the deterministic.

2. Preliminaries

A subset (classic) of the universe can be represented by its characteristic function given by

Thus, the membership function describes completely the set , since it indicates which elements of the set are elements of .

Allowing a kind of relaxation in the image set of the membership function of a set, Zadeh mathematically formulated the concept of fuzzy set [11]. We define a fuzzy subset of , or simply fuzzy set, through the membership function , where represents the degree of membership of the element to fuzzy subset .

A fuzzy rule-based system has basically four components: an input processor (or fuzzificator), a set of linguistic rules, a model of fuzzy inference, and a processor output (or defuzzificator), generating a real number as output. Figure 1 illustrates the controller fuzzy.

Figure 1: Structure of the fuzzy controller.

The fuzzification is the process by which the input values of the system are converted to fuzzy sets with respective ranges of values where they are defined. It is a mapping of the field of real numbers to fuzzy field.

We can define the fuzzy roles by structures of the form If Then used to describe specific situations that can be subjected to analysis a panel of experts, whose inference leads to a desired result. The background to define a fuzzy region in the space of input variables describes a system and condition, whereas the result define a region in the space of output variables of the system and describe a conclusion or an action that can be drawn when the premises occur. A fuzzy rule is a unit able to capture some specific knowledge.

A set of rules (or rule base) can describe a system in its various possibilities, fulfilling the role of translating mathematically the information of basis of knowledge of the fuzzy system. The rule base systems fuzzy (RBSF), in this case called fuzzy controllers, has four modules: the fuzzification module, the module based on linguistic rules, fuzzy inference, and defuzzification module. These modules are connected as shown in Figure 2. We call system p-fuzzy iterative system where is almost linear; that is, , and is obtained by a system based on fuzzy rules.

Figure 2: Structure of operation of a p-fuzzy system.

The rule base is a set consisting of fuzzy rules that relates the linguistic terms of the input variables and output variables. The rule base is considered as an element of a member of the fuzzy controller core. Each rule base satisfies the following structure: where and are fuzzy sets that represent linguistics terms of input variables and output variables, respectively. The expression in means . Both the fuzzy set and can be a cartesian product of fuzzy sets; that is, and . In this case, each fuzzy set, and , is a linguistic term for the th input variable and -th output variable, and the expression in means that

It is the definition of the rule base the information of the phenomenon under study are used. For each state defined by the linguistic terms of the input variable is a rule base. Thus, the more linguistic terms more details are incorporated in the model.

The relationship between the linguistic variable is characterized MIN by the operator; that is, each rule is considered a fuzzy relation where the degree of membership for each is

The relationship between each rule is characterized by the maximum operator, that is, fuzzy relation , which is determined by the rule base model is obtained by joining each individual rule, so that for each pair we have where represents the operator MIN.

Now we want to find each entry a corresponding action; that is, for a set of input data, we determine a set output data. Mamdani method, the membership function of , is given by

If the entry is a classic unitary set, then and . Therefore, the previous expression results in and therefore we have the fuzzy set representing the action of each input (Figure 3).

Figure 3: Mamdani inference engine with two linguistic input variables and one output.

The role of the defuzzificador is to convert each completion obtained by the method of inference in a real number which is the best action to take. For systems p-fuzzy, the actual number is obtained by defuzzification added to the value assumed by input variable in instant , feeding the iterative system.

One of the main methods of defuzzification is the center of mass, for continuous variables, which is given by

This defuzzification method will be used throughout this paper. Note that the fuzzy controller can be seen as a function , since given an input value, there is a single output value.

3. The Model

In this section, we are interested in developing a base of rules that enable us to find a solution to a graphical problem involving diffusion. For this, we use only the position population, initial condition, and population growth. By these ways we may estimate the population density in an instant without necessarily using the analytical solution problem.

Classical models of population dynamics and/or epidemiology, in overall, are given by a system of differential equations. In this case, the parameters of the models are often taken as mean values obtained from one set of data such that the model is to be deterministically known. However, admitting uncertainty due to partial knowledge, which is common in biological phenomena, an alternative is to model such knowledge from a set of rules of the form if-then.

It is common to adopt an equation to represent the dynamic system, where the field is variation, from which the evolution of the system is studied. Actually we can ask the following question: how can we analyze system (10) if it is partially known? The response is to adopt a linguistic model capable of capable of capturing the information availables from the model with a specialist. Propose a methodology to estimate solutions to differential equations using fuzzy controllers in which the state variables are the inputs and outputs are changes of state.

Thus, consider how linguistic variables to position of the population (distance to origin): low positive (Bp), mean positive (Mp), mean high positive (MAp), positive high (Ap), low negative (Bn), mean negative (Mn), mean high negative (MAn) and negative high (An). Where the positive or negative terms means distance from the origin to the right or left respectively. Thus, Figure 4 represents graphically these linguistic variables.

Figure 4: Linguistic variables for population position.

Likewise, as output variable, consider the variation population and linguistic variables: low positive , positive average , average high positive , high positive , low negative , average negative , average high negative , and high negative . In Figure 5 these plots represent the linguistic variables.

Figure 5: Linguistic variables for population growth.

Considering the known results about diffusion process, consider the following of fuzzy rules:(a)if the position is of the individuals low positive then the variation of population is low positive ;(b)if the position is of the individuals positive average then the variation of population is positive average ;(c)if the position is of the individuals average high positive then the variation of population is average high positive ;(d)if the position is of the individuals high positive then the variation of population is high positive ;(e)if the position is of the individuals low negative then the variation of population is low negative ;(f)if the position is of the individuals average negative then the variation of population is average negative ;(g)if the position is of the individuals average high negative then the variation of population is average high negative ;(h)if the position is of the individuals high negative then the variation of population is high negative .

Figure 6 shows the graph of all variation of dependent variable (change in population) as a function of independent variable (position of the population) in the Mamdani model.

Figure 6: Curve generated from the rules of the previous figure using Mamdani.

So with this rule base, the Mamdani controller and defuzzification given by the center of mass, the p-fuzzy system in leads to the path illustrated by the sequence of Figures 7, 8, 9, 10, 11, 12, 13, and 14. To this end, we consider as a condition initial .

Figure 7: Solution to .
Figure 8: Solution to .
Figure 9: Solution to .
Figure 10: Solution to .
Figure 11: Solution to .
Figure 12: Solution to .
Figure 13: Solution to .
Figure 14: Solution to .

Another important observation is that time in our problem means the number of iterations in simulated MATLAB. This information is of great importance because it reduces the amount of data needed to describe the problem. We can thus write p-fuzzy dynamical system in the following form:

As for each iteration , we have a curve; then, we can plot a surface that will be our graph p-fuzzy solution of the diffusion equation. This graph is shown in Figure 15.

Figure 15: Union of p-fuzzy solutions.

4. Conclusion

The most interesting in this process is not possible to know what the best model, it is the deterministic model or the p-fuzzy model, since the results are very similar as we can see in Figure 15. The modeling assumes that not alwayshas real data, the intuition or common sense can guide the formulation of models.

From the educational point of view the best model is secondary because you can always do better than the previous one and you can always imagine different situations for the same phenomenon.

We show that it is possible to use a FRBS to model the behavior of the population density of a species when you want to take into account the diffusion of individuals.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.


  1. P. Diamond, “Stability and periodicity in fuzzy differential equations,” IEEE Transactions on Fuzzy Systems, vol. 8, no. 5, pp. 583–590, 2000. View at Publisher · View at Google Scholar · View at Scopus
  2. P. Diamond, “Brief note on the variation of constants formula for fuzzy differential equations,” Fuzzy Sets and Systems, vol. 129, no. 1, pp. 65–71, 2002. View at Publisher · View at Google Scholar · View at MathSciNet · View at Scopus
  3. A. Kandel and W. J. Byatt, “Fuzzy differential equations,” in Proceedings of the International Conference on Cybernetics and Society, Tokyo, Japan, 1981.
  4. M. T. Malinowski, “On random fuzzy differential equations,” Fuzzy Sets and Systems, vol. 160, no. 21, pp. 3152–3165, 2009. View at Publisher · View at Google Scholar · View at MathSciNet · View at Scopus
  5. J. D. Murray, Mathematical Biology I. An Introduction, vol. 1, Springer, New York, NY, USA, 2002. View at MathSciNet
  6. M. S. Cecconello, R. C. Bassanezi, A. V. Brandão, and J. Leite, “Periodic orbits for fuzzy flows,” Fuzzy Sets and Systems, vol. 230, pp. 21–38, 2013. View at Publisher · View at Google Scholar · View at MathSciNet · View at Scopus
  7. M. S. Cecconello, J. Leite, R. C. Bassanezi, J. de Deus, and M. Silva, “About projections of solutions for fuzzy differential equations,” Journal of Applied Mathematics, vol. 2013, Article ID 184950, 9 pages, 2013. View at Publisher · View at Google Scholar · View at Scopus
  8. J. D. Silva, J. Leite, R. C. Bassanezi, and M. S. Cecconello, “Stationary points—I: one-dimensional p-fuzzy dynamical systems,” Journal of Applied Mathematics, vol. 2013, Article ID 495864, 11 pages, 2013. View at Publisher · View at Google Scholar
  9. M. T. Mizukoshi, L. C. Barros, and R. C. Bassanezi, “Stability of fuzzy dynamic systems,” International Journal of Uncertainty, Fuzziness and Knowledge-Based Systems, vol. 17, no. 1, pp. 69–83, 2009. View at Publisher · View at Google Scholar · View at MathSciNet · View at Scopus
  10. J. Leite, Sistemas dinâmicos fuzzy aplicado a processos difusivos [Ph.D. thesis], Universidade Estadual de Campinas, Campinas, Brazil.
  11. L. A. Zadeh, “Fuzzy sets,” Information and Control, vol. 8, pp. 338–353, 1965. View at Google Scholar · View at MathSciNet · View at Scopus