Building Mathematical Models for Multicriteria and Multiobjective Applications 2019
View this Special IssueResearch Article  Open Access
A Novel TOPSIS Method Based on Improved Grey Relational Analysis for Multiattribute DecisionMaking Problem
Abstract
Multiattribute decisionmaking (MADM) problem is difficult to assess because of the large number of attribute indices and the diversity of data distribution. Based on the understanding of data dispersion degree, a new grey TOPSIS method for MADM is studied. The main idea of this paper is to redefine the grey relational analysis through the dispersion of data distribution and redesign the TOPSIS by using the improved grey relational analysis. As a classical multiattribute decision analysis method, traditional TOPSIS does not consider the data distribution of the degree of dispersion and aggregation when it is compared with the optimal and worst alternative solutions. In view of the limitations of traditional TOPSIS, this paper has made two major improvements to TOPSIS. Firstly, the new grey relational analysis is applied to evaluate the grey positive relational degree between each alternative and the optimal solution and compute the grey negative relational degree between each alternative and the worst solution. Secondly, the weights of every attribute index about the optimal and worst solutions are put forward based upon the distance standard deviation and the average distance. Finally, the comprehensive grey TOPSIS is utilized to analyze the ranking of weapon selection problem. The numerical results verify the feasibility of the improved grey relational analysis and also highlight the practicability of the grey comprehensive TOPSIS.
1. Introduction
In recent years, with the development of society, the research on multiattribute decisionmaking (MADM) is becoming more and more concerned about its simplicity, validity, and accuracy [1–9]. Generally speaking, multiattribute decisionmaking problem often involves complicated external environment and many different attributes. The TOPSIS (techniques for order preference by similarity to ideal solution) is a useful and powerful method for dealing with multiattribute decisionmaking (MADM) problems which is presented by Hwang and Yoon [10]. The classical TOPSIS approach is based on the suggestion of the chosen alternative which requires consideration of the distance between the positive ideal solution and the negative ideal solution. The most suitable alternative should be closest to the positive ideal solution and the farthest from the negative ideal solution. At present, with the deepening realization about TOPSIS, some extended TOPSIS techniques have been used widely for MADM problems to overcome the vagueness in the judgments made by the assessors [11–21]. In order to overcome the problem of index weight calculation, AHP is the main auxiliary method of TOPSIS, and fuzzy sets are introduced together to facilitate the quantification of the two pairs of indicators [11–13]. Some researchers have improved TOPSIS by using the intuitionistic fuzzy number, triangular fuzzy number, vague set, intuitionistic fuzzy entropy, and other uncertainty tools, which can solve MADM problems with intuitionistic fuzzy environment or linguistic fuzziness [7, 14–18].
In fact, in the process of multiattribute decisionmaking with TOPSIS, there are two main problems to solve. On the one hand, the weights of different indexes need to be determined when the alternatives are compared with the optimal or worst schemes. As a relatively simple method, analytic hierarchy process (AHP) has been successfully applied to the weight calculation of TOPSIS [11–13, 19], but subjective interference also exists due to expert scoring and pairwise comparison. The technique of analytic network process (ANP) is used to obtain the relative weights to reduce the number of pairwise comparison on the basis of AHP [20]. So the ANP is the development and supplement of AHP. In addition, the method of comprehensive evaluation, principal component analysis, and other methods are all important methods to calculate weights, which are applied in TOPSIS [21–23]. On the other hand, the optimal or worst value for every attribute should be selected, and a quantitative measure method should be constructed to compare with each other. Fuzzy mathematics provides a quantitative calculation method for describing uncertain things [11–15, 18, 19], but there is a lot of uncertainty and subjectivity in the designing of fuzzy sets. Professor Deng proposed the grey system theory which is one of powerful methods for data analysis with partly known and partly unknown information in uncertain environment [24, 25]. Distance standard deviation is an important numerical characteristic which reflects the distribution of data. Very recently, it is worth noting that distance standard deviation is introduced into the grey confidence interval estimation for small samples [26]. This study shows that distance standard deviation contributes to describe the degree of dispersion between different samples for grey system. However, in the grey relational analysis, the value of distance standard deviation has not been explored. In view of this, this paper will introduce the distance standard deviation factor in the grey relational analysis, so as to take full account of the numerical characteristics of different series of data.
Based on the deepening realization about the generality of multiattribute decisionmaking problems and the core of TOPSIS method, part of the work has been completed in the following. In Section 2, the sample distance standard deviation is defined between different system behavior sequences, and then the improved grey relational analysis based on sample distance standard deviation is constructed, which paved the way to further discussion. In Section 3, the innovative new TOPSIS method is proposed. The degree of aggregation between different data is served as an important decision factor for TOPSIS method. The grey relational analysis is improved on the basis of the distance standard deviation and applied to quantify the correlation between different sequence data. The distance standard deviation is also used to calculate the weights between the sequence data and the optimal scheme or the worst case. Section 4 gives three applications of the developed approach, and the numerical steps and results are clearly shown. The last section summarizes the whole paper.
2. The Improved Grey Relational Analysis Based on Dispersion Degree
At present, the grey relational analysis has been applied more commonly to many scientific decisionmaking problems, such as computer science, engineering and other related fields. It provides a tool to evaluate grey relational degree based on distances between reference sequence and comparative sequence. In fact, only considering the maximum and minimum distance between data is not enough. Generally speaking, the sample aggregation characteristics reflect the inherent character of data, which should be considered in the calculation of grey relational degree. The numerical characteristics of the attribute data distribution can reflect the degree of dispersion of the data distribution. If the degree of dispersion of the data is small, the distribution is relatively concentrated, and this attribute is relatively important. And the standard deviation of sample distance of data coincides with this characteristic. Therefore, the sample distance standard deviation as an important concept that embodies the degree of data distribution density will be brought into the calculation of grey relational degree.
Definition 1. For the given system behavior sequence is as follows:The distance between and is defined as follows:For , and for every , letwhere = denotes the sample distance standard deviation between and , and is defined as the average distance between and .
Then is called grey relational degree with sample distance standard deviation, in which and are called distinguishing coefficient, , and ’s value satisfies the following requirements:
Theorem 2. The improved grey relational degree with sample distance standard deviation satisfies the following three axioms of grey correlation:(1)The property of normality: .(2)The property of closeness: the greater the , the closer of to .(3)The property of pair symmetry: .
Proof. It is obvious that for system behavior sequence . So the theorem can be proved in two cases.
When or , from (3), we can receive the result that for every , then in these circumstances, we can obtain the result that by (4). , if and only if .
For system behavior sequences and , is a constant. Because is monotonically decreasing and is monotonically decreasing, the property of closeness is received.
At last, because of , the property of pair symmetry is obvious.
If and , then . So, by (3). Because , we can receive the result that for every , then in this case.
Because and are constant which determined by system behavior sequence . It is obvious that the closeness is established.
Because , the property of pair symmetry can be received.
Remark 3. The new grey relational analysis takes full advantage of the inherent regularity of data distribution and introduces the dispersion of data distribution into them. The distance standard deviation is more likely to reflect the overall dispersion of the data distribution than the distance.
Remark 4. In the case of or , that is, there are two system behavior sequence and that are parallel or overlapped, the distance between and is fully considered by (3). Expressly if two system behavior sequence data are parallel, then the smaller the distance is, the greater the similarity is, and the greater the grey relational degree is.
3. The Novel TOPSIS with Improved Grey Relational Analysis to Solve Multiattribute DecisionMaking Problem
This section tries to improve traditional TOPSIS with new grey relational analysis mentioned above. The degree of dispersion and aggregation between different data will serve as an important decision factor for TOPSIS. The abovementioned grey relational analysis will be used to calculate the grey relational degree between the alternative scheme and the optimal or worst case scheme. At the same time, the standard deviation of distance will be taken as the decisive factor of weight calculation in order to determine the weight of each index for MADM problems.
Step 1. Standardization of the value of original decision matrix : this will help eliminating the influence of dimension and magnitude in data processing. The original decision matrix can be expressed as follows:where is the th attribute value of the th evaluation object. For multiattribute decisionmaking problem, the normalized decision matrix is after standardization of .
For the benefit attributeFor the cost attribute
Step 2. To pick out the optimal and worst alternatives, the optimal alternative is denoted as , whereThe worst alternative denoted as , whereFor the normalized decision matrix , the normalized optimal alternative can be denoted as , and the normalized worst alternative can be denoted as . It is not difficult to calculate that , and by (7) and (8), respectively.
Step 3. To calculate the sample distance standard deviation for different attribute, the symbol used to denote the sample distance standard deviation between the th attribute vector and the th attribute optimal value, whereThe symbol is used to denote the sample distance standard deviation between the th attribute vector and the th attribute worst value, where
Step 4. To compute the grey relational degree and . There are some different cases based on (3).
When and for every , , the symbol is the grey positive relational degree between and , whereAt the same case, the symbol is the grey negative relational degree between and , where When or , for every , , the grey positive relational degree can be computed as follows:And the grey negative relational degree between and can be received as follows:
Step 5. To determine weights and for the grey positive relational degree and grey negative relational degree, similar to the above, there are two cases below.
When , , , and for every , , and can be obtained as follows:where the sample distance standard deviation and can be determined by (11) and (13).
Otherwise, the weights and can be calculated as follows:where the average distance and can be determined by (12) and (14).
The smaller the distance standard deviation of attribute data or the average distance is, the stronger the regularity of attribute data is, and the bigger its weight is. On the contrary, the larger the distance standard deviation of attribute data or the average distance is, the larger the distribution range of attribute data is, the weaker the regularity is, and the smaller the weight is. It is obvious to find that the greater the is, the closer alternative attribute value and the optimal alternative attribute value are. Similarly, the weight has the same nature. The greater the is, the closer alternative attribute value and the worst value are.
Step 6. To calculate the grey TOPSIS grade, the grey positive relational grade between the ith evaluation alternative and the optimal solution can be received as follows:The grey negative relational grade between the ith evaluation alternative and the worst solution can be received as follows:And the comprehensive grey TOPSIS evaluation result is defined as follows:It is obviously that the comprehensive grey TOPSIS evaluation result reflects the closeness between alternative and positive ideal solution and negative ideal solution. The larger the is, the closer is to positive ideal solution, and the smaller the is, the farther is to negative ideal solution. The reason behind this is that after the first dimensionless standardized processing of the original data, the improved TOPSIS method is used for evaluation, and the obtained grey relational degree and weights have the function of adaptive adjustment with the change of data. On the one hand, the grey positive relational degree can reflect the closeness between the data and the same attribute data of the positive ideal solution, while the grey negative relational degree can reflect the closeness between the data and the same attribute data of the negative ideal solution. On the other hand, the weights of different attributes can feedback the importance of the attributes based on the numerical characteristics of the attribute data.
4. Numerical Results
In this section, three examples will be used to demonstrate the effectiveness of the methods proposed above. The first example shows that the improved grey relational analysis is more comprehensive than the traditional grey relational analysis, which can avoid some unreasonable results. The second one is used to illustrate the rationality of the improved grey relational analysis. The novel TOPSIS with improved grey relational analysis will be applied for assessment of weapon system in the second example. In addition, Example 3 will illustrate that the method of this paper can avoid the occurrence of the rank reversal phenomenon, which is superior to the traditional TOPSIS method.
Example 1. The following example will try to solve the grey relational degree with improved method in this paper. Take four system behavior sequences for example, , , , and . Figure 1 illustrates the changing trend curves for different system behavior sequences.
Solution
Step 1. Because the data of four system behavior sequences are in a relatively small range , the data preprocessing can be ignored.
According to (2), compute the distance , respectively, (). The symbol is used to represent the distance vector between and , .Step 2. Compute the average distance and the sample distance standard deviation .Step 3. According to (3) and (4), the improved grey relational degree based on discrete degree of data in statistics can be solved as follows:The results show that , which indicates is more relevant to compared with others. The two sets of data and are parallel and nooverlap, so they can be considered different, and the grey relational degree cannot be valued at one. Obviously, the valuing and sorting of the improved grey relational analysis results are superior to the classical method shown in Table 1. In this example, some parameters are taken as follows: , , , and .

Example 2. In this example, seven surfacetoair missile weapon system alternatives will be evaluated by grey TOPSIS method in the following. The multiattribute weapon system is evaluated from the combination of traditional grey relational analysis and TOPSIS in reference paper [27], but no convincing calculation index weight is given. In order to obtain the best selection of weapon systems, [28] developed fuzzy analytic hierarchy process by entropy weight to evaluate the weight. Metin, Serkan, and Nevzat [29] utilized AHP and TOPSIS to assess weapon systems. AHP is more difficult because of the comparison between the two factors, then this paper from the distribution of different indicators to consider the dispersion degree and facilitate the given weight. A good surfacetoair missile weapon system depends on a lot of attributes, such as the maximum speed of missile (), the maximum speed of target (), the maximum overload of target (), the highest boundary of killing range (), the farthest boundary of killing range (), the number of targets that can simultaneity be shot by one weapon system (), the single shot kill probability of missiles (), the reaction time of missile weapon system (), the lowest boundary of killing range (), the launching weight of missiles (), and the nearest boundary of killing range () as listed in Table 2 [27]. Based on the above facts, the optimal selection of multiple weapons is very important and complicated for evaluation and assessment of weapon system.
Solution
Step 1. The normalized decision matrix can be calculated by (7) and (8) as follows: Step 2. Use (9) and (10) to determine the optimal alternative and worst alternative .Step 3. The average distance and distance standard deviation are calculated, respectively, by (11)(14).where , , , .
Step 4. To determine the indicator weights,where , .
Step 5. The grey positive relational degree and grey negative relational degree can be obtained by (15) and (16), . Then, we can get the grey positive relational grade and grey negative relational grade .Step 6. At last, the comprehensive grey TOPSIS evaluation results can be derived as follows:For the sake of convenience, the data in Table 3 retains three significant digits. Obviously, the presented method in this paper is consistent with the sorting results of the reference method [27] and TOPSIS. According to the comprehensive grey TOPSIS evaluation results, the order is . In general, the ranking results of the methods are consistent, but locally, the change of evaluation values is different by Figure 2. The attribute value of the single shot kill probability of missiles (P) of X_{5}, X_{6}, and X_{7} is relatively close in Table 2. And the weights and are the biggest in the same kind of attribute index weights. And the weights of attribute indexes with large difference are relatively small. These lead to a close ranking of X_{5}, X_{6}, and X_{7}.


Example 3. For the multiattribute missile selection problem, this example will delete or add alternative to evaluate whether there is a rank reversal phenomenon based on the original data of Example 2. The traditional TOPSIS method usually uses Euclidean distance to measure, and the weights of different attributes are fixed, which easily causes the phenomenon of rank reversal [30]. In the decisionmaking process of the proposed method, when the alternatives are compared with the positive ideal solution and the negative ideal solution, the attributes have two sets of weights, and the weights can be dynamically adjusted with the change of the density of data distribution. Now we drop out of alternative set based on the second example. Then we add a new alternative to the original set of Example 2. The main results are shown in Table 4.
In the case of dropping out of alternative set, the indicator weights are updated to the following results:In the second case, the weights and for the grey positive relational degree and grey negative relational degree can be determined as follows:The above weights are two sets of weights given after considering both positive and negative ideal solutions, which can be adjusted automatically according to the characteristics of data distribution.
Based on Example 2, the ranking results of the original alternative set areSimilarly, we can obtain the following results with dropping out of alternative set:Then for the last case, we can receive the following ranking:It is not difficult to find that the ranking of the original alternatives is unchanged under the two different circumstances, so there is no rank reversal. We believe that no rank reversal is a normal phenomenon, so the novel TOPSIS method in the paper is more in line with people’s thinking logic.
Since the improved TOPSIS method in this paper has two sets of attribute weights during the evaluation process, the weights can be dynamically adjusted with the change of data distribution density no matter in reducing the evaluation objects or increasing the evaluation objects. The final result of the evaluation presented to us is that the original data will not have the phenomenon of rank reversal. Correspondingly, the traditional TOPSIS method is prone to rank reversal because the attribute weights cannot be adjusted dynamically. The specific information under the two different circumstances is shown in Table 4.

5. Conclusions
The traditional TOPSIS method is limited because it only takes the distance between the data into account, without introducing the degree of data dispersion. In addition, the classical grey relational analysis method does not study the aggregation and dispersion degree of different series of data. Based on this, a novel grey TOPSIS method developed from an improved grey relational analysis and a new indicator weight determination method, which well evaluated different objects. The improved grey relational analysis method makes full use of the discrete degree of the different sequence data and considers the overall distribution of the data, so as to avoid the shortage of the maximum and minimum distance of data. At the same time, the concept of sample distance standard deviation is proposed to solve the problem of different indicator weights determination. And there examples have successfully tested the effectiveness of the proposed method. Especially, the ranking of multiple weapons with MADM typical features verifies the feasibility for practical MADM problem. In the process of solving MADM, the present method in this paper overcomes the rank reversal problem of classical TOPSIS. It is worth noting that the proposed TOPSIS method depends on quantitative test data. As a future work, qualitative data and semantic fuzziness will also be introduced into the improved TOPSIS method to solve more complex multiattribute decision problems. In addition, we will reveal the necessity of introducing grey relational analysis into MADM problems further.
Data Availability
Data used to support the findings of this study are included within the article.
Conflicts of Interest
The authors declare that they have no conflicts of interest.
Acknowledgments
This research is supported by the National Natural Science Foundation of China (11801173, 11702094, and 61472137), the Scientific Research Project of Higher Education in Hebei Education Department (Z2017044), and the Provincial Science and Technology Foundation of Hebei (162176438). The study is also supported by Key Discipline of Probability Theory and Mathematical Statistics of North China Institute of Science and Technology (06DV09) and Open Project of Key Laboratory of Data Science and Application in Hebei Province (HBSJQ0708).
References
 G. T. Temur, “A novel multi attribute decision making approach for location decision under high uncertainty,” Applied Soft Computing, vol. 40, pp. 674–682, 2016. View at: Publisher Site  Google Scholar
 H. Bayram and R. Şahin, “A simulation based multiattribute group decision making technique with decision constraints,” Applied Soft Computing, vol. 49, pp. 629–640, 2016. View at: Publisher Site  Google Scholar
 Y. Zhang, W. Wang, and A. Bernard, “Embedding multiattribute decision making into evolutionary optimization to solve the manyobjective combinatorial optimization problems,” The Journal of Grey System, vol. 28, no. 3, pp. 124–143, 2016. View at: Google Scholar
 Q. Pang, H. Wang, and Z. Xu, “Probabilistic linguistic term sets in multiattribute group decision making,” Information Sciences, vol. 369, pp. 128–143, 2016. View at: Publisher Site  Google Scholar
 A. Baykasoğu, K. Subulan, and F. S. Karaslan, “A new fuzzy linear assignment method for multiattribute decision making with an application to spare parts inventory classification,” Applied Soft Computing, vol. 42, pp. 1–17, 2016. View at: Publisher Site  Google Scholar
 D. Xie, Y. Lu, C. Gu, and F. Li, “Scheme design considering network topology and multiattribute decisionmaking for under frequency load shedding,” International Journal of Electrical Power & Energy Systems, vol. 88, pp. 42–54, 2017. View at: Publisher Site  Google Scholar
 Y. He, H. Chen, Z. He, and L. Zhou, “Multiattribute decision making based on neutral averaging operators for intuitionistic fuzzy information,” Applied Soft Computing, vol. 27, pp. 64–76, 2015. View at: Publisher Site  Google Scholar
 S. Singh and M. S. Dasgupta, “Evaluation of research on CO2 transcritical work recovery expander using multi attribute decision making methods,” Renewable & Sustainable Energy Reviews, vol. 59, pp. 119–129, 2016. View at: Publisher Site  Google Scholar
 J. Qin, X. Liu, and W. Pedrycz, “Frank aggregation operators and their application to hesitant fuzzy multiple attribute decision making,” Applied Soft Computing, vol. 41, pp. 428–452, 2016. View at: Publisher Site  Google Scholar
 C. L. Hwang and K. Yoon, Multiple Attribute Decision Making: Methods and Applications, vol. 186, Springer, Heidelberg, Germany, 1981. View at: MathSciNet
 M. P. Amiri, “Project selection for oilfields development by using the AHP and fuzzy TOPSIS methods,” Expert Systems with Applications, vol. 37, no. 9, pp. 6218–6224, 2010. View at: Publisher Site  Google Scholar
 C.C. Sun, “A performance evaluation model by integrating fuzzy AHP and fuzzy TOPSIS methods,” Expert Systems with Applications, vol. 37, no. 12, pp. 7745–7754, 2010. View at: Publisher Site  Google Scholar
 T. Paksoy, N. Y. Pehlivan, and C. Kahraman, “Organizational strategy development in distribution channel management using fuzzy AHP and hierarchical fuzzy TOPSIS,” Expert Systems with Applications, vol. 39, no. 3, pp. 2822–2841, 2012. View at: Publisher Site  Google Scholar
 Z.L. Yue, “TOPSISbased group decisionmaking methodology in intuitionistic fuzzy setting,” Information Sciences, vol. 277, pp. 141–153, 2014. View at: Publisher Site  Google Scholar  MathSciNet
 T. Wang, J. Liu, J. Li, and C. Niu, “An integrating OWA–TOPSIS framework in intuitionistic fuzzy settings for multiple attribute decision making,” Computers & Industrial Engineering, vol. 98, pp. 185–194, 2016. View at: Publisher Site  Google Scholar
 Z. Pei, “A note on the TOPSIS method in MADM problems with linguistic evaluations,” Applied Soft Computing, vol. 36, article no. 3048, pp. 24–35, 2015. View at: Publisher Site  Google Scholar
 S. Zhou, W. Liu, and W. Chang, “An improved TOPSIS with weighted hesitant vague information,” Chaos, Solitons & Fractals, vol. 89, pp. 47–53, 2016. View at: Publisher Site  Google Scholar
 D. Joshi and S. Kumar, “Intuitionistic fuzzy entropy and distance measure based TOPSIS method for multicriteria decision making,” Egyptian Informatics Journal, vol. 15, no. 2, pp. 97–104, 2014. View at: Publisher Site  Google Scholar
 F. Torfi, R. Z. Farahani, and S. Rezapour, “Fuzzy AHP to determine the relative weights of evaluation criteria and Fuzzy TOPSIS to rank the alternatives,” Applied Soft Computing, vol. 10, no. 2, pp. 520–528, 2010. View at: Publisher Site  Google Scholar
 H. J. Shyur, “COTS evaluation using modified TOPSIS and ANP,” Applied Mathematics and Computation, vol. 177, no. 1, pp. 251–259, 2006. View at: Publisher Site  Google Scholar
 J. Li, J. Zhou, and Y. Zhang, “Comprehensive evaluation of energy efficiency based on TOPSIS for protocols of collaborative wireless sensor networks,” International Journal of Distributed Sensor Networks, vol. 2015, Article ID 640736, 13 pages, 2015. View at: Publisher Site  Google Scholar
 L.I. Tong, C.H. Wang, and H.C. Chen, “Optimization of multiple responses using principal component analysis and technique for order preference by similarity to ideal solution,” The International Journal of Advanced Manufacturing Technology, vol. 27, no. 34, pp. 407–414, 2005. View at: Publisher Site  Google Scholar
 Z. Zhang, Y. Wang, and Z. Wang, “A grey TOPSIS method based on weighted relational coefficient,” The Journal of Grey System, vol. 26, no. 2, pp. 112–123, 2014. View at: Google Scholar
 D. JuLong, “Control problems of grey systems,” Systems & Control Letters, vol. 1, no. 5, pp. 288–294, 1982. View at: Publisher Site  Google Scholar
 L. Sifeng, Y. Yingjie, W. Lifeng et al., Grey System Theory and Its Application, Science Press, Beijing, China, 7th edition, 2014.
 Y. Wenguang, W. Yunjie, and W. Jianmin, “Grey confidence interval estimation of small samples based on numerical characteristics,” The Journal of Grey System, vol. 29, no. 2, pp. 113–124, 2017. View at: Google Scholar
 H. Gu and B. Song, “Study on effectiveness evaluation of weapon systems based on grey relational analysis and TOPSIS,” Journal of Systems Engineering and Electronics, vol. 20, no. 1, pp. 106–111, 2009. View at: Google Scholar
 M. Dağdeviren, S. Yavuz, and N. Kilinç, “Weapon selection using the AHP and TOPSIS methods under fuzzy environment,” Expert Systems with Applications, vol. 36, no. 4, pp. 8143–8151, 2009. View at: Publisher Site  Google Scholar
 D.L. Mon, C.H. Cheng, and J.C. Lin, “Evaluating weapon system using fuzzy analytic hierarchy process based on entropy weight,” Fuzzy Sets and Systems, vol. 62, no. 2, pp. 127–134, 1994. View at: Publisher Site  Google Scholar
 Y.M. Wang and Y. Luo, “On rank reversal in decision analysis,” Mathematical and Computer Modelling, vol. 49, no. 56, pp. 1221–1229, 2009. View at: Publisher Site  Google Scholar  MathSciNet
Copyright
Copyright © 2019 Wenguang Yang and Yunjie Wu. 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.