
The crude oil market, as a complex evolutionary nonlinear driving system, is by nature a highly noisy, nonlinear and deterministic chaotic series of price series. In this paper, a computational intelligence-based portfolio model is constructed to forecast crude oil prices using weekly price data of West Texas intermediate crude oil (WTI) crude oil futures from 2011 to 2021. First, the WTI crude oil price series are decomposed using the ensemble empirical modal decomposition method (EEMD) and the set of component series is reconstructed using the cluster analysis method. Second, the reconstructed series are modeled and predicted using neural network models such as time-delay neural network (TDNN), extreme learning machine (ELM), multilayer perceptron (MLP) and the GM (1, 1) gray prediction algorithm and the output of the model with the best prediction effect for each component is integrated. Finally, the EGARCH model is used to further optimize the predictive power of the combined model and output the final predicted values. The results show that the combined model based on computational intelligence has higher forecasting accuracy than single models such as GM (1, 1), ARIMA, MLP and the combined EEMD-ELM model for forecasting crude oil futures prices.
Citation: Ming Li, Ying Li. Research on crude oil price forecasting based on computational intelligence[J]. Data Science in Finance and Economics, 2023, 3(3): 251-266. doi: 10.3934/DSFE.2023015
[1] | Yefu Zheng, Jun Xu, Hongzhang Chen . TOPSIS-based entropy measure for intuitionistic trapezoidal fuzzy sets and application to multi-attribute decision making. Mathematical Biosciences and Engineering, 2020, 17(5): 5604-5617. doi: 10.3934/mbe.2020301 |
[2] | Qiwen Wang, Guibao Song, Xiuxia Yang . Mixed-attitude three-way decision model for aerial targets: Threat assessment based on IF-VIKOR-GRA method. Mathematical Biosciences and Engineering, 2023, 20(12): 21514-21536. doi: 10.3934/mbe.2023952 |
[3] | Yongfeng Yin, Routing Zhang, Qingran Su . Threat assessment of aerial targets based on improved GRA-TOPSIS method and three-way decisions. Mathematical Biosciences and Engineering, 2023, 20(7): 13250-13266. doi: 10.3934/mbe.2023591 |
[4] | Fankang Bu, Jun He, Haorun Li, Qiang Fu . Interval-valued intuitionistic fuzzy MADM method based on TOPSIS and grey correlation analysis. Mathematical Biosciences and Engineering, 2020, 17(5): 5584-5603. doi: 10.3934/mbe.2020300 |
[5] | Nian Zhang, Yifan Zhou, Qiang Pan, Guiwu Wei . Multi-attribute decision-making method with triangular fuzzy numbers based on regret theory and the catastrophe progression method. Mathematical Biosciences and Engineering, 2022, 19(12): 12013-12030. doi: 10.3934/mbe.2022559 |
[6] | Bo Sun, Ming Wei, Wei Wu, Binbin Jing . A novel group decision making method for airport operational risk management. Mathematical Biosciences and Engineering, 2020, 17(3): 2402-2417. doi: 10.3934/mbe.2020130 |
[7] | Muhammad Akram, G. Muhiuddin, Gustavo Santos-García . An enhanced VIKOR method for multi-criteria group decision-making with complex Fermatean fuzzy sets. Mathematical Biosciences and Engineering, 2022, 19(7): 7201-7231. doi: 10.3934/mbe.2022340 |
[8] | Nian Zhang, Xue Yuan, Jin Liu, Guiwu Wei . Stochastic multiple attribute decision making with Pythagorean hesitant fuzzy set based on regret theory. Mathematical Biosciences and Engineering, 2023, 20(7): 12562-12578. doi: 10.3934/mbe.2023559 |
[9] | Yuting Zhu, Wenyu Zhang, Junjie Hou, Hainan Wang, Tingting Wang, Haining Wang . The large-scale group consensus multi-attribute decision-making method based on probabilistic dual hesitant fuzzy sets. Mathematical Biosciences and Engineering, 2024, 21(3): 3944-3966. doi: 10.3934/mbe.2024175 |
[10] | Sumera Naz, Muhammad Akram, Mohammed M. Ali Al-Shamiri, Mohammed M. Khalaf, Gohar Yousaf . A new MAGDM method with 2-tuple linguistic bipolar fuzzy Heronian mean operators. Mathematical Biosciences and Engineering, 2022, 19(4): 3843-3878. doi: 10.3934/mbe.2022177 |
The crude oil market, as a complex evolutionary nonlinear driving system, is by nature a highly noisy, nonlinear and deterministic chaotic series of price series. In this paper, a computational intelligence-based portfolio model is constructed to forecast crude oil prices using weekly price data of West Texas intermediate crude oil (WTI) crude oil futures from 2011 to 2021. First, the WTI crude oil price series are decomposed using the ensemble empirical modal decomposition method (EEMD) and the set of component series is reconstructed using the cluster analysis method. Second, the reconstructed series are modeled and predicted using neural network models such as time-delay neural network (TDNN), extreme learning machine (ELM), multilayer perceptron (MLP) and the GM (1, 1) gray prediction algorithm and the output of the model with the best prediction effect for each component is integrated. Finally, the EGARCH model is used to further optimize the predictive power of the combined model and output the final predicted values. The results show that the combined model based on computational intelligence has higher forecasting accuracy than single models such as GM (1, 1), ARIMA, MLP and the combined EEMD-ELM model for forecasting crude oil futures prices.
Because the social and economic environment are more complex than ever, many practical decision making problems are difficult to solve, and it is a need for us to develop multiple criteria/ attribute decision making methods to deal with various problems. Recent years, multi-criteria decision making (MCDM) and multi-attribute decision making (MADM) methods have been received great attention in solving various engineering problems [1,2,3,4,5]. In actual decision process, many qualitative attributes (criteria, indicators) in practical decision-making problems are difficult to quantify [6]. Some of these attributes are inherently ambiguous, such as the evaluation index of "problem-solving ability" in human resources management, which is itself ambiguous. The fuzziness of evaluation attributes is due to the influence of subjective and objective factors such as decision maker's time, energy and incomplete understanding of objective things, such as "safety and reliability" as an evaluation index in aircraft type selection. Zadeh [7] proposed the concept of fuzzy sets to solve these fuzzy phenomena in real-life problems. Because Zadeh's fuzzy sets can better depict the fuzziness of the objective things themselves and the fuzziness of human thinking, the fuzzy multi-attribute decision-making methods based on fuzzy sets have been greatly developed and applied in various fields [8,9].
Zadeh's fuzzy set has only one degree of membership and can only reflect two aspects of information, which results in some limitations in many real decision making problems. For example, with the increasing complexity of social and economic activities, decision makers often hesitate to some extent or show a certain degree of lack of knowledge in decision-making, which make the evaluation results show three aspects: support, against and hesitation. Because Zadeh's fuzzy set has only one degree of membership and can only reflect two aspects of information, it is difficult to describe the above situation better. So Atanassov [10,11] extended Zadeh's theory of fuzzy set, and put forward the concept of intuitionistic fuzzy (IF) set and interval-valued intuitionistic fuzzy (IVIF) set. Compared with Zadeh's fuzzy set, IF set and IVIF set can describe the concept of "either one or the other" by adding a non-membership parameter and considering the information of membership degree, non-membership degree and hesitation degree. Therefore, they can better describe the hesitation and uncertainty degree of decision makers' judgments, and then they can describe the model of objective things more delicately. Because IF sets and IVIF sets are more flexible than fuzzy sets in dealing with fuzziness and uncertainty, they are widely used in decision-making field [12,13,14,15,16,17]. Some other extensions of Zadeh's fuzzy set are also proposed and attached great attention to many authors. For example, picture fuzzy set [18], single valued neutrosophic set [19], intuitionistic fuzzy soft set [20,21,22], Pythagorean fuzzy soft set [23], interval type-2 fuzzy set [24] and q-rung orthopair fuzzy set [25,26] are all very famous successful extensions of Zadeh's fuzzy set.
In order to use the concept of intuitionistic fuzzy sets to define an uncertain quantity or a quantity difficult to quantify, some scholars put forward the concept of continuous intuitionistic fuzzy numbers. For example, Shu et al. [27] first introduced the concept of intuitionistic triangular fuzzy number and applied it to fault tree analysis. On this basis, Wang and Zhang [28] developed the concept of intuitionistic trapezoidal fuzzy number (ITFN). In recent years, the ITFNs have attracted many scholars' attention [29,30,31].
Distance measure is an important information measure in the study of intuitionistic fuzzys. But there are few distance formulas for ITFNs proposed in documents. Wang and Zhang [28] proposed a distance formula for two ITFNs. But we find that Wang and Zhang’s distance formula exists shortcomings, for example, $\tilde A$ = < [0.2, 0.3, 0.4, 0.5], 0.1, 0.7 > and $\tilde B$ = < [0.4, 0.6, 0.8, 1.0], 0.1, 0.9 > are two ITFNs. When we use Wang and Zhang’s distance formula, we obtain the distance of $\tilde A$ and $\tilde B$ is 0, which is counter-intuitive. Then this paper will proposed a new distance to overcome this shortcomings.
Most of the existing MADM models under intuitionistic trapezoidal fuzzy environment are based on the expected utility theory assuming that decision makers are completely rational. However, in the real MADM process, decision makers often have subjective risk preferences for alternatives, so it is important to consider the real-life attitude of decision makers in the MADM model. The prospect theory, first proposed by Kahneman and Tversky [32], is a challenge to expected utility theory, as it explains the phenomena that expected utility theory can not explain. It is more in line with the actual decision-making mode of human beings to introduce the cumulative prospect theory into the fuzzy decision model [33,34].
There are few documents discussed the application of prospect theory in intuitionistic trapezoidal fuzzy MADM problems. Chen et al. [35] developed a decision making method combining with parameter estimation and score functions for solving the MADM problem under intuitionistic trapezoidal fuzzy environment. But they defined a prospect value function of ITFNs based on the difference of score functions between ITFN. In prospect theory, a prospect value function is mainly defined by distance, which can get a well discrimination [36]. Thus this paper will first propose a new distance for ITFNs, then use it to define prospect values, and furthermore develop a new decision-making method for solving MADM problem under intuitionistic trapezoidal fuzzy environment with incomplete information based on cumulative prospect theory.
Some basic concepts and properties of ITFN and cumulative prospect theory are reviewed. To overcome the limitation of Zadeh’s fuzzy set, Atanassov proposed the concept of intuitionistic fuzzy set, which is illustrated in Definition 1.
Definition 1 [10]. Let $X$ be a given universal set, a set $A = \{ < {x_i}, {\mu _A}({x_i}), {\upsilon _A}({x_i}) > |{x_i} \in X\} $is called an intuitionistic fuzzy set, where ${\mu _A}({x_i})$ and ${\upsilon _A}({x_i})$ are the membership degree and non-membership degree of ${x_i}$, respectively. Here, $0 \leqslant {\mu _A}({x_i}) + {\upsilon _A}({x_i}) \leqslant 1$ for any${x_i} \in X$ and an index ${\pi _A}({x_i}) = 1 - {\mu _A}({x_i}) - {\upsilon _A}({x_i})$ is called the hesitancy degree of ${x_i}$. For convenience, if there is only one element in $X$, we call $A$ intuitionistic fuzzy number.
In order to use the concept of intuitionistic fuzzy sets to define an uncertain quantity or a quantity difficult to quantify, Wang and Zhang [28] introduced the concept of ITFN, which is illustrated in Definition 2.
Definition 2. Let $a, b, c, d \in R$, $0 \leqslant {\mu _{\tilde A}} \leqslant 1$, $0 \leqslant {\upsilon _{\tilde A}} \leqslant 1$ and${\mu _{\tilde A}} + {\upsilon _{\tilde A}} \leqslant 1$. A fuzzy number $\tilde A = < [a, b, c, d]; {\mu _{\tilde A}}, {\upsilon _{\tilde A}} > $ is called an ITFN, if its membership degree and non-membership degree functions are ([28]).
${\lambda _{\tilde A}}(x) = \left\{ \begin{gathered} {\mu _{\tilde A}}(x - a)/(b - a), \;\;a \leqslant x < b \\ \quad \quad \;{\mu _{\tilde A}}, \quad \quad \quad b \leqslant x \leqslant c \\ {\mu _{\tilde A}}(d - x)/(d - c), \;\;c < x \leqslant d \\ \quad \quad \;0, \quad \quad \quad \;other \\ \end{gathered} \right., $ | (1) |
and
${\gamma _{\tilde A}}(x) = \left\{ \begin{gathered} \frac{{b - x + {\upsilon _{\tilde A}}(x - a)}}{{b - a}}, \;\;a \leqslant x < b \\ \quad \quad {\upsilon _{\tilde A}}, \quad \quad \;\;\;b \leqslant x \leqslant c \\ \frac{{x - c + {\upsilon _{\tilde A}}(d - x)}}{{d - c}}, c < x \leqslant d \\ \quad \quad 1, \quad \quad \quad \;\quad other \\ \end{gathered} \right., $ | (2) |
respectively.
To rank the order of ITFNs, Wang and Zhang [29] introduced the concepts of the score and accuracy functions of an ITFN$\tilde A = < [a, b, c, d]; {\mu _{\tilde A}}, {\upsilon _{\tilde A}} > $, which is similar to that concepts of intuitionistic fuzzy numbers. Based on score function $S(\tilde A) = \frac{1}{4}(a + b + c + d) \times ({\mu _{\tilde A}} - {\upsilon _{\tilde A}})$ and accuracy function $H(\tilde A) = \frac{1}{4}(a + b + c + d) \times ({\mu _{\tilde A}} + {\upsilon _{\tilde A}})$. Score function $S(\tilde A)$ can represent a support degree if we choose a vote case as example, while accuracy function$H(\tilde A) = \frac{1}{4}(a + b + c + d) \times (1 - {\pi _{\tilde A}})$can reflect the deterministic information of $\tilde A$ in some degree. Thus Wang and Zhang [29] proposed a ranking rule as in Definition 3.
Definition 3 [29]. Let $\tilde A = < [a, b, c, d]; {\mu _{\tilde A}}, {\upsilon _{\tilde A}} > $ be an ITFN. Then the score and accuracy functions of $\tilde A$ are $S(\tilde A)$ and $H(\tilde A)$, respectively. Then for two ITFN$\tilde A$and $\tilde B$,
(1) If $S(\tilde A) > S(\tilde B)$, then $\tilde A > \tilde B$
(2) If $S(\tilde A) = S(\tilde B)$ and
(ⅰ) $H(\tilde A) = H(\tilde B)$, then $\tilde A = \tilde B$;
(ⅱ) $H(\tilde A) > H(\tilde B)$,then $\tilde A > \tilde B$.
As a very useful tool for describing the number of fuzzy numbers, $\alpha $-cut set is very important for studing the properties of fuzzy numbers. Let $\tilde A = < [{a_1}, {a_2}, {a_3}, {a_4}]; {\mu _{\tilde A}}, {\upsilon _{\tilde A}} > $ be an ITFN, and the corresponding $\alpha $-cut is ${\tilde a_\alpha } = \{ x|{\mu _{\tilde a}}(x) \geqslant \alpha \} $, $\beta $-cut is ${\tilde a_\beta } = \{ x|{\upsilon _{\tilde a}}(x) \leqslant \beta \} $. Here $0 \leqslant \alpha \leqslant {\mu _A}$, ${\upsilon _A} \leqslant \beta \leqslant 1$ and $0 \leqslant \alpha + \beta \leqslant 1$. By straightforward calculation, we can get
$\begin{array}{l} {{\tilde A}_\alpha } = [\tilde A_\alpha ^L, \tilde A_\alpha ^U] = [{a_1} + \frac{{({a_2} - {a_1})\alpha }}{{{\mu _{\tilde A}}}}, {a_4} - \frac{{({a_4} - {a_3})\alpha }}{{{\mu _{\tilde A}}}}], \\ {{\tilde A}_\beta } = [\tilde A_\beta ^L, \tilde A_\beta ^U] = [{a_1} + \frac{{({a_2} - {a_1})(1 - \beta )}}{{1 - {\upsilon _{\tilde A}}}}, {a_4} - \frac{{({a_4} - {a_3})(1 - \beta )}}{{1 - {\upsilon _{\tilde A}}}}]. \\ \end{array} $ | (3) |
Prospect theory, firstly proposed by Kahneman and Tversky [32], can objectively reflect decision makers’ behavior and the theory has been developed into solving MADM problems [33,34]. This paper will use the prospect theory to determine weights of indices with the help of a new distance measure of ITFNs. Some key elements of prospect theory are the prospect value and value function, which are given in the following definitions.
Definition 4 [32]. Prospect value, determined by value functions and decision weight, is defined as follows:
$V = \sum\limits_{j = 1}^n {\pi ({p_j}} )\nu ({x_j})$ | (4) |
where $\pi ({p_j})$ is the probability weight function which is monotone increasing with respect to ${p_j}$, and $\nu (x)$ is the value function.
Definition 5. Value function, first proposed by Kahneman and Tversky [32], has the following formula
$\nu (x) = \left\{ \begin{gathered} \quad \;{x^\alpha }\quad , x > 0 \\ - \theta {( - x)^\beta }, x < 0 \\ \end{gathered} \right.$ | (5) |
where $\nu (x)$ comes from the subjective feeling of the decision maker and the variable $x$ is the gains ($x > 0$) or the losses ($x < 0$), parameters $\alpha, \beta \in (0, 1)$ are the concave-convex degree of the region value power function of the gains and the losses, respectively. $\theta $ shows that the region value power function is more steeper for the losses than for the gains, and $\theta > 1$ means the losses aversion.
Definition 6. The probability weight, proposed by Kahneman and Tversky [32] has the following form:
$\pi (p) = \frac{{{p^\gamma }}}{{{{({p^\gamma } + {{(1 - p)}^\gamma })}^{1/\gamma }}}}$ | (6) |
where $\gamma > 0$ is a parameter measuring degree of curvature.
Let $A$ and $B$ be two fuzzy numbers. Scholars have constructed many distance formulas of $A$ and $B$, such as Euclidian distance, Hamming distance [37]. Grzegorzewski [38] proposed a new distance formula as follows:
Let $A$ and $B$ be two fuzzy number, the corresponding $\alpha $-cuts are $[{({A_L})_\alpha }, {({A_U})_\alpha }]$ and $[{({B_L})_\alpha }, {({B_U})_\alpha }]$, respectively. The distance measure of $A$ and $B$ is defined as
$d(A, B) = {\left( {\int_0^1 {{{({{({A_L})}_\alpha } - {{({B_L})}_\alpha })}^2}d\alpha } + \int_0^1 {{{({{({A_U})}_\alpha } - {{({B_U})}_\alpha })}^2}d\alpha } } \right)^{1/2}}$ | (7) |
The distance formula (7) has received great attention, for the reason that many authors used it to obtain nearest interval, triangular and trapezoidal approximation of a fuzzy number preserving ambiguity [39,40]. Motivated by the distance formula (7), this section will construct a new distance formula for ITFNs.
For two any ITFNs $\tilde A$ and $\tilde B$, the most existing distance formulas $d(\tilde A, \tilde B)$ could not consider the influence of $\alpha - {\rm{cut}}$ and $\beta - {\rm{cut}}$, which are two important concepts of fuzzy sets. There exist limitations of most distance formulas when we use them to discriminate some intuitionistic fuzzy numbers [41]. This section will consider the $\alpha - {\rm{cut}}$ and $\beta - {\rm{cut}}$ into the construction of the new distance formula. The new distance measure of $\tilde A$ and $\tilde B$ is constructed as follows:
$\begin{array}{l} d(\tilde A, \tilde B) = {\left( {\int_0^1 {{{(\tilde A_\alpha ^L - \tilde B_\alpha ^L)}^2}d\alpha } + \int_0^1 {{{(\tilde A_\alpha ^U - \tilde B_\alpha ^U)}^2}d\alpha } } \right)^{1/2}} \\ \quad \quad \quad \quad + {\left( {\int_0^1 {{{(\tilde A_\beta ^L - \tilde B_\beta ^L)}^2}d\beta } + \int_0^1 {{{(\tilde A_\beta ^U - \tilde B_\beta ^U)}^2}d\beta } } \right)^{1/2}} \\ \quad \quad \quad \quad + \frac{1}{2}(|{\mu _{\tilde A}} - {\mu _{\tilde B}}| + |{v_{\tilde A}} - {v_{\tilde B}}| + |{\pi _{\tilde A}} - {\pi _{\tilde B}}|) \\ \end{array} $ | (8) |
Lemma 1. (Cauchy-Schwarz Inequality) Let $f(x)$ and $g(x)$ be two integral functions on [$a, b$], then
${\left[ {\int_a^b {f(x)g(x)dx} } \right]^2} \leqslant \int_a^b {{f^2}(x)dx} \cdot \int_a^b {{g^2}(x)dx} $ | (9) |
By straightforward calculation, we have the following result seen in Lemma 2.
Lemma 2. For four positive real numbers:$x_1, x_2, y_1, y_2,$ then
$\sqrt {{x_1}{y_1}} + \sqrt {{x_2}{y_2}} \leqslant \sqrt {{x_1} + {x_2}} \cdot \sqrt {{y_1} + {y_2}}. $ | (10) |
Lemma 3. Let $f_1,f_2,g_1,g_2$ and $g$ are real functions of variable $x$, and they are integral on [$a,b$], set
${d_1} = \sqrt {\int_a^b {{{({f_1} - {f_2})}^2}} dx + \int_a^b {{{({g_1} - {g_2})}^2}} dx} , $ |
${d_2} = \sqrt {\int_a^b {{{({f_1} - f)}^2}} dx + \int_a^b {{{({g_1} - g)}^2}} dx} , $ |
${d_3} = \sqrt {\int_a^b {{{(f - {f_2})}^2}} dx + \int_a^b {{{(g - {g_2})}^2}} dx} .$ |
Then
${d_1} \leqslant {d_2} + {d_3}$ |
Proof.
$\begin{array}{l} d_1^2 = \int_a^b {{{({f_1} - {f_2})}^2}} dx + \int_a^b {{{({g_1} - {g_2})}^2}} dx \\ \;\;\, {\kern 1pt} = \int_a^b {{{({f_1} - f + f - {f_2})}^2}} dx + \int_a^b {{{({g_1} - g + g - {g_2})}^2}} dx \\ \;\;{\kern 1pt} {\kern 1pt} {\kern 1pt} = \int_a^b {{{({f_1} - f)}^2}} dx + \int_a^b {{{(f - {f_2})}^2}} dx + \int_a^b {{{({g_1} - g)}^2}} dx + \int_a^b {{{(g - {g_2})}^2}} dx \\ \quad \; + 2\int_a^b {({f_1} - f)(f - {f_2})} dx + 2\int_a^b {({g_1} - g)(g - {g_2})} dx \\ \end{array} $ |
According to Lemma 1 and Lemma 2, we have
$\begin{array}{l} d_1^2 \leqslant d_2^2 + d_3^2{\rm{ + 2}}\sqrt {\int_a^b {{{({f_1} - f)}^2}} dx{\rm{ + }}\int_a^b {{{({g_1} - g)}^2}} dx} {\rm{ + 2}}\sqrt {\int_a^b {{{(f - {f_2})}^2}} dx{\rm{ + }}\int_a^b {{{(g - {g_2})}^2}} dx} \\ \quad {\rm{ = (}}{d_2} + {d_2}{)^2} \\ \end{array} $ |
Then we can easily complete the proof.
Theorem 1. Let $\tilde A$, $\tilde B$ and $\tilde C$ be three ITFNs, then $d(\tilde A, \tilde B)$ is a distance measure of $\tilde A$ and $\tilde B$. That is, the following three conditions are satisfied: (ⅰ)$d(\tilde A, \tilde B) \geqslant 0$; (ⅱ)$d(\tilde A, \tilde B) = d(\tilde B, \tilde A)$ and (ⅲ)$d(\tilde A, \tilde B) \leqslant d(\tilde A, \tilde C) + d(\tilde B, \tilde C)$.
Proof. Obviously, (ⅰ) and (ⅱ) are true. Now we prove (ⅲ). Set,
$d(\tilde A, \tilde B) = {d_1}(\tilde A, \tilde B) + {d_2}(\tilde A, \tilde B) + {d_3}(\tilde A, \tilde B)$ |
where
${d_1}(\tilde A, \tilde B) = {\left( {\int_0^1 {{{(\tilde A_\alpha ^L - \tilde B_\alpha ^L)}^2}d\alpha } + \int_0^1 {{{(\tilde A_\alpha ^U - \tilde B_\alpha ^U)}^2}d\alpha } } \right)^{1/2}}, $ |
${d_2}(\tilde A, \tilde B) = {\left( {\int_0^1 {{{(\tilde A_\beta ^L - \tilde B_\beta ^L)}^2}d\beta } + \int_0^1 {{{(\tilde A_\beta ^U - \tilde B_\beta ^U)}^2}d\beta } } \right)^{1/2}}, $ |
${d_3}(\tilde A, \tilde B) = \frac{1}{2}(|{\mu _{\tilde A}} - {\mu _{\tilde B}}| + |{v_{\tilde A}} - {v_{\tilde B}}| + |{\pi _{\tilde A}} - {\pi _{\tilde B}}|).$ |
According to Lemma 3 and ${d_3}(\tilde A, \tilde B)$ is a formula of Minkowski distance, we can immediately get (ⅲ).
Therefore $d(\tilde A, \tilde B)$ is a distance measure of $\tilde A$ and $\tilde B$.
Let $\tilde A$ and $\tilde B$ be two ITFNs, where $\tilde A = < [{a_1}, {a_2}, {a_3}, {a_4}]; {\mu _{\tilde A}}, {\upsilon _{\tilde A}} > $ and $\tilde B = < [{b_1}, {b_2}, {b_3}, {b_4}]; {\mu _{\tilde B}}, {\upsilon _{\tilde B}} > $. Let ${S_1} = 1 + {\mu _{\tilde A}} - {\upsilon _{\tilde A}}$ and ${S_2} = 1 + {\mu _{\tilde B}} - {\upsilon _{\tilde B}}$. Then Wang and Zhang’s distance measure [28] of $\tilde A$ and $\tilde B$ is
${d_{WZ}}(\tilde A, \tilde B) = \frac{1}{8}(|{S_1}{a_1} - {S_2}{b_1}| + |{S_1}{a_2} - {S_2}{b_2}| + |{S_1}{a_3} - {S_2}{b_3}| + |{S_1}{a_4} - {S_2}{b_4}|)$ | (11) |
Example 1. Let us consider two ITFNs = < [0.2, 0.3, 0.4, 0.5], 0.1, 0.7 > and = < [0.4, 0.6, 0.8, 1.0], 0.1, 0.9 > . Then, using Wang and Zhang’s distance measure (11), we obtain the distance of and is 0, which is counter-intuitive. Here we use the proposed distance formula (8) to calculate the distance between them. According to Eq (3), we have
$\begin{array}{l} {{\tilde A}_\alpha } = [\tilde A_\alpha ^L, \tilde A_\alpha ^U] = [0.2 + \alpha , 0.5 - \alpha ], \\ {{\tilde A}_\beta } = [\tilde A_\beta ^L, \tilde A_\beta ^U] = [0.2 + \frac{1}{3}(1 - \beta ), 0.5 - \frac{1}{3}(1 - \beta )], \\ {{\tilde B}_\alpha } = [\tilde B_\alpha ^L, \tilde B_\alpha ^U] = [0.4 + 2\alpha , 1 - 2\alpha ], \\ {{\tilde B}_\beta } = [\tilde B_\beta ^L, \tilde B_\beta ^U] = [0.4 + 2(1 - \beta ), 1 - 2(1 - \beta )]. \\ \end{array} $ |
Substituting them into formula (8), we can easily get the calculation result $d(\tilde A, \tilde B) = 2.2917$. This result shows that the proposed distance formula has a better discrimination ability than Wang and Zhang’s distance measure (11).
Definition 7 [34]. Let two $\tilde A$ and $\tilde B$ be two then we can define the value function of $\tilde B$ by using $\tilde A$ as a reference point as follows:
$\nu (\tilde B) = \left\{ \begin{gathered} \;\;\;{(d(\tilde A, \tilde B))^\alpha }\;\;, \tilde B \geqslant \tilde A \\ - \theta {(d(\tilde A, \tilde B))^\beta }, \;\tilde B < \tilde A \\ \end{gathered} \right., $ | (12) |
Where $d(\tilde A, \tilde B)$ is the distance of $\tilde A$ and $\tilde B$.
In real-life decision-making problems, different attitudes are often accompanied by different decision results due to different conditions and external environments of decision-makers. Some scholars have considered the influence of different attitudes of decision-makers for solving MADM problems under interval number and triangular fuzzy number environment. But few literature considered the behavior factors of decision maker in solving MADM problem in which attribute values are expressed by ITFNs.
In this section, we will propose a new decision method of intuitionistic trapezoidal fuzzy model, which takes the decision maker's attitude into consideration based on prospect theory combining with the new distance measure of ITFNs. The detail decision process is shown in the Figure 2.
Now, we first introduce the intuitionistic trapezoidal fuzzy multi-attribute decision making problem as follows:
(ⅰ) Alternative set $\{ {A_1}, {A_2}, ..., {A_m}\} $ and attribute set $O = \{ {o_1}, {o_2}, ..., {o_n}\} $;
(ⅱ) Because the importance of different attributes is often different, and ${w_j}$ is the importance degree of attribute ${o_j}$, and they satisfy $\sum\limits_{j = 1}^n {{w_j}} = 1$ and ${w_j} \geqslant 0, $ $j = 1, 2, ..., n$.
(ⅲ) For alternative ${A_i}$, the evaluation value of ${A_i}$ with respect to ${o_j}$ is an ITFN ${\tilde a_{ij}} = < [{a_{ij1}}, {a_{ij2}}, {a_{ij3}}, {a_{ij4}}]; {\mu _{{{\tilde a}_{ij}}}}, {v_{{{\tilde a}_{ij}}}} > $, then one can get a fuzzy decision making matrix $\tilde A = {({\tilde a_{ij}})_{m \times n}}$.
(ⅳ) The weight information is partially known, and the set of incomplete information is defined as $H$.
The task of decision maker is to give a rational ranking result for given alternatives. Now we develop a new decision making method based on prospect theory combining with the new distance measure. The specific decision-making steps are as follows:
Step 1. Normalize $\tilde A = {({\tilde a_{ij}})_{m \times n}}$ to $\tilde R = {({\tilde r_{ij}})_{m \times n}}$.
In order to eliminate the influence of different physical dimensions on decision results, it is necessary to standardize the decision matrix. Set $M = \{ 1, 2, ..., m\}, N = \{ 1, 2, ..., n\} $. ${I_1}$ and ${I_2}$ represent the subscript sets of benefit type and cost type respectively. The matrix $\tilde A = {({\tilde a_{ij}})_{m \times n}}$ is transformed into a normalized matrix $\tilde R = {({\tilde r_{ij}})_{m \times n}}$ by using the following formulas, where ${\tilde r_{ij}} = < [{r_{ij1}}, {r_{ij2}}, {r_{ij3}}, {r_{ij4}}]; {\mu _{{{\tilde r}_{ij}}}}, {v_{{{\tilde r}_{ij}}}} > $,
${r_{ijk}} = \frac{{{a_{ijk}} - \mathop {\min }\limits_i \{ {a_{ij1\;}}\} }}{{\mathop {\max }\limits_i \{ {a_{ij4\;}}\} - \mathop {\min }\limits_i \{ {a_{ij1\;}}\} }}, i \in M, j \in {I_1}, k = 1, 2, 3, 4, $ | (13) |
${r_{ijk}} = \frac{{\mathop {\max }\limits_i \{ {a_{ij4\;}}\} - {a_{ij(4 - k)}}}}{{\mathop {\max }\limits_i \{ {a_{ij4\;}}\} - \mathop {\min }\limits_i \{ {a_{ij1\;}}\} }}, i \in M, j \in {I_2}, k = 1, 2, 3, 4$ | (14) |
and
${\mu _{{{\tilde r}_{ij}}}} = {\mu _{{{\tilde a}_{ij}}}}, {v_{{{\tilde r}_{ij}}}} = {v_{{{\tilde a}_{ij}}}}.$ | (15) |
Step 2. Define the positive ideal solution (PIS) and negative ideal solution (PIS).
PIS is defined as ${\tilde r^ + } = \{ \tilde r_1^ +, \tilde r_2^ +, ..., \tilde r_n^ + \} $and NIS is defined as ${\tilde r^ - } = \{ \tilde r_1^ -, \tilde r_2^ -, ..., \tilde r_n^ - \} $, where $\tilde r_j^ + = < (1, 1, 1, 1), 1, 0 > $ and $\tilde r_j^ - = < (0, 0, 0, 0), 0, 1 > $.
Step 3. Determine the gain and loss values of the attributes of each alternative.
When the NIS is taken as the reference point, the alternative ${A_i}$ is better than the NIS. For the decision maker, he/she is faced with benefits. According to the prospect theory, he/she is averse in risk at this time. Then the positive prospect value is defined by ${v^ + }({r_{ij}}) = {(d({\tilde r_{ij}}, r_j^ -))^\alpha }$, where
$\begin{array}{l} d({{\tilde r}_{ij}}, \tilde r_j^ - ) = {\left( {\int_0^1 {{{(({{\tilde r}_{ij}})_\alpha ^L - 0)}^2}d\alpha } + \int_0^1 {{{(({{\tilde r}_{ij}})_\alpha ^U - 0)}^2}d\alpha } } \right)^{1/2}} \\ \quad \quad \quad \quad + {\left( {\int_0^1 {{{(({{\tilde r}_{ij}})_\beta ^L - 0)}^2}d\beta } + \int_0^1 {{{(({{\tilde r}_{ij}})_\beta ^U - 0)}^2}d\beta {\kern 1pt} } } \right)^{1/2}} \\ \quad \quad \quad \quad + \frac{1}{2}(|{\mu _{{{\tilde r}_{ij}}}}| + |{\upsilon _{{{\tilde r}_{ij}}}} - 1| + |{\mu _{{{\tilde r}_{ij}}}} + {\upsilon _{{{\tilde r}_{ij}}}}|) \\ \end{array}.$ | (16) |
The distance decision matrix from the alternative ${x_i}$ to the negative ideal point is obtained as follows:
${D^ - } = {(d({\tilde r_{ij}}, \tilde r_j^ - ))_{m \times n}} = \left[ {\begin{array}{*{20}{c}} {d({{\tilde r}_{11}}, \tilde r_1^ - )}&{d({{\tilde r}_{12}}, \tilde r_2^ - )}& \cdots &{d({{\tilde r}_{1n}}, \tilde r_n^ - )} \\ {d({{\tilde r}_{21}}, \tilde r_1^ - )}&{d({{\tilde r}_{22}}, \tilde r_2^ - )}& \cdots &{d({{\tilde r}_{2n}}, \tilde r_n^ - )} \\ \vdots & \vdots & \ddots & \vdots \\ {d({{\tilde r}_{m1}}, \tilde r_1^ - )}&{d({{\tilde r}_{m2}}, \tilde r_2^ - )}& \cdots &{d({{\tilde r}_{mn}}, \tilde r_n^ - )} \end{array}} \right]$ | (17) |
When the PIS is taken as the reference point, the alternative ${A_i}$ is inferior to the PIS. For the decision maker, he/she is faced with loss. According to the prospect theory, he/she is pursuing risk averse at this time. At this time, the negative prospect value is defined by ${v^ - }({\tilde r_{ij}}) = - \theta {(d({\tilde r_{ij}}, r_j^ +))^\beta }$, where
$\begin{array}{l} d({{\tilde r}_{ij}}, \tilde r_j^ + ) = {\left( {\int_0^1 {{{(({{\tilde r}_{ij}})_\alpha ^L - 1)}^2}d\alpha } + \int_0^1 {{{(({{\tilde r}_{ij}})_\alpha ^U - 1)}^2}d\alpha } } \right)^{1/2}} \\ \quad \quad \quad + {\left( {\int_0^1 {{{(({{\tilde r}_{ij}})_\beta ^L - 1)}^2}d\beta } + \int_0^1 {{{(({{\tilde r}_{ij}})_\beta ^U - 1)}^2}d\beta } } \right)^{1/2}} \\ \quad \quad \quad + \frac{1}{2}(|{\mu _{{{\tilde r}_{ij}}}} - 1| + |{\upsilon _{{{\tilde r}_{ij}}}}| + |{\mu _{{{\tilde r}_{ij}}}} + {\upsilon _{{{\tilde r}_{ij}}}}|) \\ \end{array} .$ | (18) |
Then the distance decision matrix from the alternative ${x_i}$ to the positive ideal point is obtained as follows:
${D^ + } = {(d({\tilde r_{ij}}, \tilde r_j^ + ))_{m \times n}} = \left[ {\begin{array}{*{20}{c}} {d({{\tilde r}_{11}}, \tilde r_1^ + )}&{d({{\tilde r}_{12}}, \tilde r_2^ + )}& \cdots &{d({{\tilde r}_{1n}}, \tilde r_n^ + )} \\ {d({{\tilde r}_{21}}, \tilde r_1^ + )}&{d({{\tilde r}_{22}}, \tilde r_2^ + )}& \cdots &{d({{\tilde r}_{2n}}, \tilde r_n^ + )} \\ \vdots & \vdots & \ddots & \vdots \\ {d({{\tilde r}_{m1}}, \tilde r_1^ + )}&{d({{\tilde r}_{m2}}, \tilde r_2^ + )}& \cdots &{d({{\tilde r}_{mn}}, \tilde r_n^ + )} \end{array}} \right]$ | (19) |
Step 4. Determine the positive prospect and negative prospect values.
When the NIS is taken as the reference point, the alternative ${A_i}$ is better than the NIS. For the decision maker, he is faced with benefits. According to the prospect theory, the decision maker is risk averse at this time. At this time, the positive prospect value is defined by ${v^ + }({r_{ij}}) = {(d({\tilde r_{ij}}, r_j^ -))^\alpha }$ and abbreviated by $v_{ij}^ + $.
When the PIS is taken as the reference point, the alternative ${A_i}$ is inferior to the PIS. For the decision maker, he is faced with loss. According to the prospect theory, the decision maker is pursuing risk at this time. At this time, the negative prospect value is defined by ${v^ - }({\tilde r_{ij}}) = - \theta {(d({\tilde r_{ij}}, r_j^ +))^\beta }$ and abbreviated by $v_{ij}^ - $.
Step 5. Compute the comprehensive prospect value of ${x_i}$.
Assume that ${\pi ^ + }({w_j})$ and ${\pi ^ - }({w_j})$ are the prospect weight functions, which represent the weights of decision maker facing to the gains and losses, respectively. Then the comprehensive prospect value, noted by ${V_i}$, of alternative ${A_i}$ is the sum of positive prospect values and negative prospect values. Then we can get
${V_i} = \sum\limits_{j = 1}^n {{\nu ^ + }({{\tilde r}_{ij}})} {\pi ^ + }({w_j}) + \sum\limits_{j = 1}^n {{\nu ^ - }({{\tilde r}_{ij}})} {\pi ^ - }({w_j}), $ | (20) |
Where
${\pi ^\delta }({w_j}) = \frac{{w_j^{\gamma \delta }}}{{{{[w_j^{\gamma \delta } + {{(1 - {w_j})}^{\gamma \delta }}]}^{1/\gamma \delta }}}}$ | (21) |
and $v_{ij}^ + $ and $v_{ij}^ - $are the value functions. According to [34], in this paper we will adopt the parameter values as follows:
$\alpha = \beta = 0.88, \theta = 2.25, {\gamma ^ + } = 0.61, {\gamma ^ - } = 0.69.$ |
Step 6. Determine the optimal attribute weights.
Let $w = ({w_1}, {w_2}, ..., {w_n})$ and$w \in H$. For each alternative${A_i}$, its comprehensive prospect value is always the larger the better, so we need maximum $V = ({V_1}, {V_2}, ..., {V_m})$ under weights information. Then we can construct the following optimization model
$\begin{array}{l} \max \;V = ({V_1}, {V_2}, ..., {V_m}) \\ s.t.\left\{\begin{array}{l} w \in H \\ \sum\limits_{j = 1}^n {{w_j}} = 1 \\ {w_j} \geqslant 0, j = 1, 2, ..., n \\ \end{array} \right. \\ \end{array} $ | (22) |
Because of the fair competition among the alternatives, we establish the following optimization model
$\begin{array}{l} \max \;\tilde V = \sum\limits_{i = 1}^m {{V_i}} = \sum\limits_{i = 1}^m {\sum\limits_{j = 1}^n {v_{ij}^ + } {\pi ^ + }({w_j})} + \sum\limits_{i = 1}^m {\sum\limits_{j = 1}^n {v_{ij}^ - } {\pi ^ - }({w_j})} \\ \quad s.t.\left\{ \begin{array}{l} w \in H \\ \sum\limits_{j = 1}^n {{w_j}} = 1 \\ {w_j} \geqslant 0, j = 1, 2, ..., n \\ \end{array} \right. \\ \end{array} $ | (23) |
The optimal weight vector ${w^ * } = (w_1^ *, w_2^ *, ..., w_n^ *)$ can be obtained by MATLAB or Lingo software.
Step 7. Calculate the optimal comprehensive prospect value of each alternative $A_i(i=1,2,3,…,m)$ as
$V_i^ * = \sum\limits_{j = 1}^n {v_{ij}^ + } {\pi ^ + }(w_j^ * ) + \sum\limits_{j = 1}^n {v_{ij}^ - } {\pi ^ - }(w_j^ * )$ | (24) |
By sorting $V_i^ * $ $(i = 1, 2, ..., m)$ from large to small, the sorting results and the optimal scheme of the whole alternative set can be obtained.
Taking the example [29] of supplier selection for an engine parts manufacturing company to illustrate the performance of the proposed method. There are five alternative suppliers${A_1}, {A_2}, {A_3}, {A_4}, {A_5}$, and the five evaluation attributes: supply capability (${o_1}$), delivery capability (${o_2}$), service quality (${o_3}$), influence (${o_4}$), scientific research strength (${o_5}$). Decision makers give the intuitionistic trapezoidal decision matrix provided in Table 1. The attributes’ weight information of given by decision makers is partially known, $H = {\rm{\{ (}}{w_1}, {w_2}, \ldots, {w_5})|0.10 \leqslant {w_1} \leqslant 0.15, 0.10 \leqslant {w_2}{\rm{, }}{w_5} \leqslant 0.30, $ $0.05 \leqslant {w_3} \leqslant 0.20, $ $0.08 \leqslant {w_4} \leqslant 0.18, {w_2} \geqslant {w_5} \geqslant {w_1} \geqslant {w_4} \geqslant {w_3}\} $. Our task is to determine the best supplier.
${o_1}$ | ${o_2}$ | ${o_3}$ | ${o_4}$ | ${o_5}$ | |
${A_1}$ | <[1,2,3,4]; 0.7, 0.3> | <[5,6,7,8]; 0.7, 0.3> | <[3,4,5,6]; 0.6, 0.2> | <[4,5,7,8]; 0.6, 0.3> | <[4,5,7,8]; 0.8, 0.0> |
${A_2}$ | <[2,3,4,5]; 0.6, 0.3> | <[6,7,8,9]; 0.8, 0.1> | <[4,5,6,7]; 0.8, 0.1> | <[3,4,5,6]; 0.7, 0.3> | <[6,7,8,9]; 0.6, 0.3> |
${A_3}$ | <[1,2,3,5]; 0.6, 0.4> | <[4,6,7,8]; 0.6, 0.3> | <[3,4,5,6]; 0.5, 0.5> | <[4,5,6,7]; 0.8, 0.1> | <[5,6,7,8]; 0.8, 0.2> |
${A_4}$ | <[2,3,4,6]; 0.6, 0.2> | <[5,6,7,8]; 0.8, 0.2> | <[2,3,5,6]; 0.6, 0.4> | <[3,4,5,7]; 0.6, 0.3> | <[4,6,7,8]; 0.6, 0.3> |
${A_5}$ | <[2,3,4,5]; 0.8, 0.2> | <[4,5,6,7]; 0.9, 0.0> | <[3,4,5,6]; 0.8, 0.2> | <[3,5,7,8]; 0.7, 0.1> | <[4,5,6,7]; 0.8, 0.0> |
Now we give the detail calculation steps of this example.
Step 1. Normalize $\tilde A$ by the normalization formulas (13)–(15), we can get $\tilde R$ shown in Table 2.
${o_1}$ | ${o_2}$ | ${o_3}$ | ${o_4}$ | ${o_5}$ | |
${A_1}$ | <[0, 0.2, 0.4, 0.6];0.7, 0.3> | <[0.2, 0.4, 0.6, 0.8;0.7, 0.3> | <[0.2, 0.4, 0.6, 0.8;0.6, 0.2> | <[0.2, 0.4, 0.8, 1.0];0.6, 0.3> | <[0, 0.2, 0.6, 0.8];0.8, 0.0> |
${A_2}$ | <[0.2, 0.4, 0.6, 0.8];0.6, 0.3> | <[0.4, 0.6, 0.8, 1.0];0.8, 0.1> | <[0.4, 0.6, 0.8, 1.0];0.8, 0.1> | <[0, 0.2, 0.4, 0.6];0.7, 0.3> | <[0.4, 0.6, 0.8, 1.0];0.6, 0.3> |
${A_3}$ | <[0, 0.2, 0.4, 0.8];0.6, 0.4> | <[0, 0.4, 0.6, 0.8];0.6, 0.3> | <[0.2, 0.4, 0.6, 0.8];0.5, 0.5> | <[0.2, 0.4, 0.6, 0.8];0.8, 0.1> | <[0.2, 0.4, 0.6, 0.8];0.8, 0.2> |
${A_4}$ | <[0.2, 0.4, 0.6, 1.0];0.6, 0.2> | <[0.2, 0.4, 0.6, 0.8];0.8, 0.2> | <[0, 0.2, 0.6, 0.8];0.6, 0.4> | <[0, 0.2, 0.4, 0.8];0.6, 0.3> | <[0, 0.4, 0.6, 0.8];0.6, 0.3> |
${A_5}$ | <[0.2, 0.4, 0.6, 0.8];0.8, 0.2> | <[0, 0.2, 0.4, 0.6];0.9, 0.0> | <[0.2, 0.4, 0.6, 0.8];0.8, 0.2> | <[0, 0.4, 0.8, 1.0];0.7, 0.1> | <[0, 0.2, 0.4, 0.6];0.8, 0.0> |
Step2. Determine the PIS ${\tilde r^ + } = (\tilde r_1^ +, \tilde r_2^ +, ..., \tilde r_5^ +)$, where $\tilde r_j^ + = < (1, 1, 1, 1), 1, 0 > $ and NIS${\tilde r^ - } = (\tilde r_1^ -, \tilde r_2^ -, ..., \tilde r_5^ -)$, where $\tilde r_j^ - = < (0, 0, 0, 0), 0, 1 > $ for all $j = 1, 2, ..., 5$.
Step3. Calculate the distance measures between each attribute of alternative ${A_i}$ with that of positive and negative ideal points, and these values are shown in Tables 3 and 4.
$\tilde r_1^ + $ | $\tilde r_2^ + $ | $\tilde r_3^ + $ | $\tilde r_4^ + $ | $\tilde r_5^ + $ | |
${A_1}$ | 2.3425 | 1.8007 | 1.9004 | 1.7514 | 2.0912 |
${A_2}$ | 1.8947 | 1.2112 | 1.2112 | 2.3425 | 1.3768 |
${A_3}$ | 2.4250 | 2.0042 | 1.9787 | 1.7174 | 1.7122 |
${A_4}$ | 1.8850 | 1.9004 | 2.3298 | 2.1994 | 1.8467 |
${A_5}$ | 1.9004 | 2.2583 | 1.9004 | 1.9187 | 2.4497 |
$\tilde r_1^ - $ | $\tilde r_2^ - $ | $\tilde r_3^ - $ | $\tilde r_4^ - $ | $\tilde r_5^ - $ | |
${A_1}$ | 1.6859 | 2.2007 | 2.3004 | 2.5511 | 2.4058 |
${A_2}$ | 2.1947 | 2.9549 | 2.9549 | 1.6859 | 2.7382 |
${A_3}$ | 1.6805 | 2.1873 | 1.9787 | 2.4174 | 2.3122 |
${A_4}$ | 2.4239 | 2.3004 | 1.9221 | 1.8514 | 2.1776 |
${A_5}$ | 2.3004 | 2.0181 | 2.3004 | 2.7504 | 2.0002 |
Step 4. Calculate ${v^ + }({r_{ij}})$ and${v^ - }({\tilde r_{ij}})$, and we can get the positive prospect matrix and negative prospect matrix as shown in Tables 5 and 6.
${o_1}$ | ${o_2}$ | ${o_3}$ | ${o_4}$ | ${o_5}$ | |
${A_1}$ | 1.5835 | 2.0019 | 2.0816 | 2.2799 | 2.1653 |
${A_2}$ | 1.9971 | 2.5946 | 2.5946 | 1.5835 | 2.4264 |
${A_3}$ | 1.5790 | 1.9912 | 1.8231 | 2.1744 | 2.0909 |
${A_4}$ | 2.1795 | 2.0816 | 1.7771 | 1.7195 | 1.9834 |
${A_5}$ | 2.0816 | 1.8550 | 2.0816 | 2.4359 | 1.8405 |
${o_1}$ | ${o_2}$ | ${o_3}$ | ${o_4}$ | ${o_5}$ | |
${A_1}$ | –4.7589 | –3.7754 | –3.9589 | –3.6844 | –4.3065 |
${A_2}$ | –3.9484 | –2.6632 | –2.6632 | –4.7589 | –2.9811 |
${A_3}$ | –4.9061 | –4.1485 | –4.1021 | –3.6213 | –3.6116 |
${A_4}$ | –3.9306 | –3.9589 | –4.7360 | –4.5020 | –3.8601 |
${A_5}$ | –3.9589 | –4.6080 | –3.9589 | –3.9924 | –4.9499 |
Step 5. Set ${\gamma ^ + } = 0.61, {\gamma ^ - } = 0.69$. Because of solving the nonlinear optimization model (23), we first choose the initial point ${w_0}$ = (0.15;0.25;0.10;0.10;0.20). With the help of matlab software, the optimal weight vector is obtained as ${w^ * } = (0.15, 0.30, 0.10, 0.15, 0.30)$.
Step 6. Calculate the optimal comprehensive prospect value of each alternative ${A_i}$, and we get $V_1^ * = $ –1.7631, $V_2^ * = $ –0.9588, $V_3^ * = $ –1.7610, $V_4^ * = $ –1.8069, $V_5^ * = $ –2.0599.
Step 7. According to the order of the optimal comprehensive prospect value, the ranking result is: ${A_2} > {A_3} > {A_1} > {A_4} > {A_5}$. Therefore, the optimal alternative is ${A_2}$.
In Example 1, the optimal alternative is in consistent with the ranking results of Wang and Zhang [29], but other alternatives’ order is different. The reason is that the they didn’t consider psychological factors of decision makers in the decision process.
Most of the existing MADM models under intuitionistic trapezoidal fuzzy environment are based on the expected utility theory assuming that decision makers are completely rational. However, in the actual MADM process, decision makers often have subjective risk preferences, such as psychological and behavioral factors for alternatives. So it is important to consider the risk attitude of decision makers in the decision process. Prospect theory can describe human's psychological behavior under uncertain conditions more truthfully, and it can explain the phenomena that expected utility theory can not explain. Our decision making method is considering the psychological factors of decision makers based on prospect theory. The main contribution of this article is as follows:
This article firstly constructs a new distance formula of ITFNs, and the distance measure can well recognize the difference between two TIFNs. Then aiming at the problem of MADM with incomplete attribute weight information and attribute evaluation values expressed by ITFNs, a new MADM method based on prospect theory combining with the new proposed distance formula is developed. According to the distance formula, the value function of ITFNs is given, and the positive and negative ideal solutions are taken as reference, based on these, we adopt different ways to deal with "gain" and "loss" so that decision makers can get more accurate decision weights. The intuitionistic trapezoidal fuzzy MADM method based on prospect theory proposed in this paper has the characteristics of clear concept, simple calculation process and easy software implementation, and can be easily applied in practice.
The proposed decision making method can be used to solve other MADM problems, such as the investment decision, project evaluation, etc. Also, the proposed method has its limitation: How to determine the most suitable values of parameters in prospect theory. In the future, we will further develop similar constructing method proposed in this paper into the study of new distance formulas for trapezoidal fuzzy neutrosophic number, which is a useful generalized intuitionistic fuzzy number.
This research was financially supported by National Natural Science Foundation of China (No.71661012), Key Projects of Science Research of Jiangxi Educational Committee (Grant No. GJJ170496 and GJJ160629). The authors are grateful to the reviewers for their careful reading of the manuscript and the suggestions which are helpful to the improvement of the paper.
The authors declared that they have no conflicts of interest to this work.
[1] |
Abdoos AA (2016) A new intelligent method based on combination of VMD and ELM for short term wind power forecasting. Neurocomputing 203: 111–120. https://doi.org/10.1016/j.neucom.2016.03.054 doi: 10.1016/j.neucom.2016.03.054
![]() |
[2] |
Akram QF (2020) Oil price drivers, geopolitical uncertainty and oil exporters' currencies. Energ Econo 89: 104801. https://doi.org/10.1016/j.eneco.2020.104801 doi: 10.1016/j.eneco.2020.104801
![]() |
[3] |
Baumeister C, Kilian L (2015) Forecasting the real price of oil in a changing world: a forecast combination approach. J Bu Econ Stat 33: 338–351. https://doi.org/10.1080/07350015.2014.949342 doi: 10.1080/07350015.2014.949342
![]() |
[4] |
Baumeister C, Kilian L (2016) Understanding the Decline in the Price of Oil since June 2014. J Assoc Enviro Resour Econ 3: 131–158. https://doi.org/10.1086/684160 doi: 10.1086/684160
![]() |
[5] |
Bouoiyour J, Selmi R, Hammoudeh S, et al. (2019) What are the categories of geopolitical risks that could drive oil prices higher? Acts or threats? Energ Econ 84: 104523. https://doi.org/10.1016/j.eneco.2019.104523 doi: 10.1016/j.eneco.2019.104523
![]() |
[6] |
Caldara D, Iacoviello M (2022) Measuring geopolitical risk. Am Econ Rev 112: 1194–1225. https://doi.org/10.1257/aer.20191823 doi: 10.1257/aer.20191823
![]() |
[7] |
Chen HH, Chen M, Chiu C C (2016) The integration of artificial neural networks and text mining to forecast gold futures prices. Commu Stat-Simul C 45: 1213–1225. https://doi.org/10.1080/03610918.2013.786780 doi: 10.1080/03610918.2013.786780
![]() |
[8] |
He M, Zhang Y, Wen D, et al. (2021) Forecasting crude oil prices: A scaled PCA approach. Energ Econ 97: 105189. https://doi.org/10.1016/j.eneco.2021.105189 doi: 10.1016/j.eneco.2021.105189
![]() |
[9] | Huifeng L (2017) Price forecasting of stock index futures based on a new hybrid EMD-RBF neural network model. Agro Food Ind Hi-Tech 28: 1744–1747. |
[10] |
Morana C (2001) A semiparametric approach to short-term oil price forecasting. Energ Econ 23: 325–338. https://doi.org/10.1016/S0140-9883(00)00075-X doi: 10.1016/S0140-9883(00)00075-X
![]() |
[11] |
Phan DHB, Narayan PK, Gong Q (2021) Terrorist attacks and oil prices: Hypothesis and empirical evidence. Int Re Financ Anal 74: 101669. https://doi.org/10.1016/j.irfa.2021.101669 doi: 10.1016/j.irfa.2021.101669
![]() |
[12] |
Plourde A, Watkins GC (1998) Crude oil prices between 1985 and 1994: how volatile in relation to other commodities? Resour Energy Econ 20: 245–262. https://doi.org/10.1016/S0928-7655(97)00027-4 doi: 10.1016/S0928-7655(97)00027-4
![]() |
[13] | Tian LH, Tan DK (2015) Analysis of the factors influencing crude oil prices: financial speculation or Chinese demand? Economics (Quarterly) 14: 961–982. |
[14] |
Wang C (2016) Forecast on price of agricultural futures in China based on ARIMA model. Asian Agric Res 8: 9. https://doi.org/10.22004/ag.econ.253255 doi: 10.22004/ag.econ.253255
![]() |
[15] |
Wu Z, Huang NE (2009) Ensemble empirical mode decomposition: a noise-assisted data analysis method. Adv Adapt Data Analysis 1: 1–41. https://doi.org/10.1142/S1793536909000047 doi: 10.1142/S1793536909000047
![]() |
[16] |
Yin L, Yang Q (2016) Predicting the oil prices: Do technical indicators help? Energ Econ 56: 338–350. https://doi.org/10.1016/j.eneco.2016.03.017 doi: 10.1016/j.eneco.2016.03.017
![]() |
1. | Jawad Ali, Zia Bashir, Tabasam Rashid, Wali Khan Mashwani, A q-rung orthopair hesitant fuzzy stochastic method based on regret theory with unknown weight information, 2022, 1868-5137, 10.1007/s12652-022-03746-8 | |
2. | Xiaoli Tian, Jiangshui Ma, Liu Li, Zeshui Xu, Ming Tang, Development of prospect theory in decision making with different types of fuzzy sets: A state-of-the-art literature review, 2022, 615, 00200255, 504, 10.1016/j.ins.2022.10.016 | |
3. | Junli Cao, Fangfang Xu, Qiangyi Li, Entropy-Based Fuzzy TOPSIS Method for Investment Decision Optimization of Large-Scale Projects, 2022, 2022, 1687-5273, 1, 10.1155/2022/4381293 |
${o_1}$ | ${o_2}$ | ${o_3}$ | ${o_4}$ | ${o_5}$ | |
${A_1}$ | <[1,2,3,4]; 0.7, 0.3> | <[5,6,7,8]; 0.7, 0.3> | <[3,4,5,6]; 0.6, 0.2> | <[4,5,7,8]; 0.6, 0.3> | <[4,5,7,8]; 0.8, 0.0> |
${A_2}$ | <[2,3,4,5]; 0.6, 0.3> | <[6,7,8,9]; 0.8, 0.1> | <[4,5,6,7]; 0.8, 0.1> | <[3,4,5,6]; 0.7, 0.3> | <[6,7,8,9]; 0.6, 0.3> |
${A_3}$ | <[1,2,3,5]; 0.6, 0.4> | <[4,6,7,8]; 0.6, 0.3> | <[3,4,5,6]; 0.5, 0.5> | <[4,5,6,7]; 0.8, 0.1> | <[5,6,7,8]; 0.8, 0.2> |
${A_4}$ | <[2,3,4,6]; 0.6, 0.2> | <[5,6,7,8]; 0.8, 0.2> | <[2,3,5,6]; 0.6, 0.4> | <[3,4,5,7]; 0.6, 0.3> | <[4,6,7,8]; 0.6, 0.3> |
${A_5}$ | <[2,3,4,5]; 0.8, 0.2> | <[4,5,6,7]; 0.9, 0.0> | <[3,4,5,6]; 0.8, 0.2> | <[3,5,7,8]; 0.7, 0.1> | <[4,5,6,7]; 0.8, 0.0> |
${o_1}$ | ${o_2}$ | ${o_3}$ | ${o_4}$ | ${o_5}$ | |
${A_1}$ | <[0, 0.2, 0.4, 0.6];0.7, 0.3> | <[0.2, 0.4, 0.6, 0.8;0.7, 0.3> | <[0.2, 0.4, 0.6, 0.8;0.6, 0.2> | <[0.2, 0.4, 0.8, 1.0];0.6, 0.3> | <[0, 0.2, 0.6, 0.8];0.8, 0.0> |
${A_2}$ | <[0.2, 0.4, 0.6, 0.8];0.6, 0.3> | <[0.4, 0.6, 0.8, 1.0];0.8, 0.1> | <[0.4, 0.6, 0.8, 1.0];0.8, 0.1> | <[0, 0.2, 0.4, 0.6];0.7, 0.3> | <[0.4, 0.6, 0.8, 1.0];0.6, 0.3> |
${A_3}$ | <[0, 0.2, 0.4, 0.8];0.6, 0.4> | <[0, 0.4, 0.6, 0.8];0.6, 0.3> | <[0.2, 0.4, 0.6, 0.8];0.5, 0.5> | <[0.2, 0.4, 0.6, 0.8];0.8, 0.1> | <[0.2, 0.4, 0.6, 0.8];0.8, 0.2> |
${A_4}$ | <[0.2, 0.4, 0.6, 1.0];0.6, 0.2> | <[0.2, 0.4, 0.6, 0.8];0.8, 0.2> | <[0, 0.2, 0.6, 0.8];0.6, 0.4> | <[0, 0.2, 0.4, 0.8];0.6, 0.3> | <[0, 0.4, 0.6, 0.8];0.6, 0.3> |
${A_5}$ | <[0.2, 0.4, 0.6, 0.8];0.8, 0.2> | <[0, 0.2, 0.4, 0.6];0.9, 0.0> | <[0.2, 0.4, 0.6, 0.8];0.8, 0.2> | <[0, 0.4, 0.8, 1.0];0.7, 0.1> | <[0, 0.2, 0.4, 0.6];0.8, 0.0> |
$\tilde r_1^ + $ | $\tilde r_2^ + $ | $\tilde r_3^ + $ | $\tilde r_4^ + $ | $\tilde r_5^ + $ | |
${A_1}$ | 2.3425 | 1.8007 | 1.9004 | 1.7514 | 2.0912 |
${A_2}$ | 1.8947 | 1.2112 | 1.2112 | 2.3425 | 1.3768 |
${A_3}$ | 2.4250 | 2.0042 | 1.9787 | 1.7174 | 1.7122 |
${A_4}$ | 1.8850 | 1.9004 | 2.3298 | 2.1994 | 1.8467 |
${A_5}$ | 1.9004 | 2.2583 | 1.9004 | 1.9187 | 2.4497 |
$\tilde r_1^ - $ | $\tilde r_2^ - $ | $\tilde r_3^ - $ | $\tilde r_4^ - $ | $\tilde r_5^ - $ | |
${A_1}$ | 1.6859 | 2.2007 | 2.3004 | 2.5511 | 2.4058 |
${A_2}$ | 2.1947 | 2.9549 | 2.9549 | 1.6859 | 2.7382 |
${A_3}$ | 1.6805 | 2.1873 | 1.9787 | 2.4174 | 2.3122 |
${A_4}$ | 2.4239 | 2.3004 | 1.9221 | 1.8514 | 2.1776 |
${A_5}$ | 2.3004 | 2.0181 | 2.3004 | 2.7504 | 2.0002 |
${o_1}$ | ${o_2}$ | ${o_3}$ | ${o_4}$ | ${o_5}$ | |
${A_1}$ | 1.5835 | 2.0019 | 2.0816 | 2.2799 | 2.1653 |
${A_2}$ | 1.9971 | 2.5946 | 2.5946 | 1.5835 | 2.4264 |
${A_3}$ | 1.5790 | 1.9912 | 1.8231 | 2.1744 | 2.0909 |
${A_4}$ | 2.1795 | 2.0816 | 1.7771 | 1.7195 | 1.9834 |
${A_5}$ | 2.0816 | 1.8550 | 2.0816 | 2.4359 | 1.8405 |
${o_1}$ | ${o_2}$ | ${o_3}$ | ${o_4}$ | ${o_5}$ | |
${A_1}$ | –4.7589 | –3.7754 | –3.9589 | –3.6844 | –4.3065 |
${A_2}$ | –3.9484 | –2.6632 | –2.6632 | –4.7589 | –2.9811 |
${A_3}$ | –4.9061 | –4.1485 | –4.1021 | –3.6213 | –3.6116 |
${A_4}$ | –3.9306 | –3.9589 | –4.7360 | –4.5020 | –3.8601 |
${A_5}$ | –3.9589 | –4.6080 | –3.9589 | –3.9924 | –4.9499 |
${o_1}$ | ${o_2}$ | ${o_3}$ | ${o_4}$ | ${o_5}$ | |
${A_1}$ | <[1,2,3,4]; 0.7, 0.3> | <[5,6,7,8]; 0.7, 0.3> | <[3,4,5,6]; 0.6, 0.2> | <[4,5,7,8]; 0.6, 0.3> | <[4,5,7,8]; 0.8, 0.0> |
${A_2}$ | <[2,3,4,5]; 0.6, 0.3> | <[6,7,8,9]; 0.8, 0.1> | <[4,5,6,7]; 0.8, 0.1> | <[3,4,5,6]; 0.7, 0.3> | <[6,7,8,9]; 0.6, 0.3> |
${A_3}$ | <[1,2,3,5]; 0.6, 0.4> | <[4,6,7,8]; 0.6, 0.3> | <[3,4,5,6]; 0.5, 0.5> | <[4,5,6,7]; 0.8, 0.1> | <[5,6,7,8]; 0.8, 0.2> |
${A_4}$ | <[2,3,4,6]; 0.6, 0.2> | <[5,6,7,8]; 0.8, 0.2> | <[2,3,5,6]; 0.6, 0.4> | <[3,4,5,7]; 0.6, 0.3> | <[4,6,7,8]; 0.6, 0.3> |
${A_5}$ | <[2,3,4,5]; 0.8, 0.2> | <[4,5,6,7]; 0.9, 0.0> | <[3,4,5,6]; 0.8, 0.2> | <[3,5,7,8]; 0.7, 0.1> | <[4,5,6,7]; 0.8, 0.0> |
${o_1}$ | ${o_2}$ | ${o_3}$ | ${o_4}$ | ${o_5}$ | |
${A_1}$ | <[0, 0.2, 0.4, 0.6];0.7, 0.3> | <[0.2, 0.4, 0.6, 0.8;0.7, 0.3> | <[0.2, 0.4, 0.6, 0.8;0.6, 0.2> | <[0.2, 0.4, 0.8, 1.0];0.6, 0.3> | <[0, 0.2, 0.6, 0.8];0.8, 0.0> |
${A_2}$ | <[0.2, 0.4, 0.6, 0.8];0.6, 0.3> | <[0.4, 0.6, 0.8, 1.0];0.8, 0.1> | <[0.4, 0.6, 0.8, 1.0];0.8, 0.1> | <[0, 0.2, 0.4, 0.6];0.7, 0.3> | <[0.4, 0.6, 0.8, 1.0];0.6, 0.3> |
${A_3}$ | <[0, 0.2, 0.4, 0.8];0.6, 0.4> | <[0, 0.4, 0.6, 0.8];0.6, 0.3> | <[0.2, 0.4, 0.6, 0.8];0.5, 0.5> | <[0.2, 0.4, 0.6, 0.8];0.8, 0.1> | <[0.2, 0.4, 0.6, 0.8];0.8, 0.2> |
${A_4}$ | <[0.2, 0.4, 0.6, 1.0];0.6, 0.2> | <[0.2, 0.4, 0.6, 0.8];0.8, 0.2> | <[0, 0.2, 0.6, 0.8];0.6, 0.4> | <[0, 0.2, 0.4, 0.8];0.6, 0.3> | <[0, 0.4, 0.6, 0.8];0.6, 0.3> |
${A_5}$ | <[0.2, 0.4, 0.6, 0.8];0.8, 0.2> | <[0, 0.2, 0.4, 0.6];0.9, 0.0> | <[0.2, 0.4, 0.6, 0.8];0.8, 0.2> | <[0, 0.4, 0.8, 1.0];0.7, 0.1> | <[0, 0.2, 0.4, 0.6];0.8, 0.0> |
$\tilde r_1^ + $ | $\tilde r_2^ + $ | $\tilde r_3^ + $ | $\tilde r_4^ + $ | $\tilde r_5^ + $ | |
${A_1}$ | 2.3425 | 1.8007 | 1.9004 | 1.7514 | 2.0912 |
${A_2}$ | 1.8947 | 1.2112 | 1.2112 | 2.3425 | 1.3768 |
${A_3}$ | 2.4250 | 2.0042 | 1.9787 | 1.7174 | 1.7122 |
${A_4}$ | 1.8850 | 1.9004 | 2.3298 | 2.1994 | 1.8467 |
${A_5}$ | 1.9004 | 2.2583 | 1.9004 | 1.9187 | 2.4497 |
$\tilde r_1^ - $ | $\tilde r_2^ - $ | $\tilde r_3^ - $ | $\tilde r_4^ - $ | $\tilde r_5^ - $ | |
${A_1}$ | 1.6859 | 2.2007 | 2.3004 | 2.5511 | 2.4058 |
${A_2}$ | 2.1947 | 2.9549 | 2.9549 | 1.6859 | 2.7382 |
${A_3}$ | 1.6805 | 2.1873 | 1.9787 | 2.4174 | 2.3122 |
${A_4}$ | 2.4239 | 2.3004 | 1.9221 | 1.8514 | 2.1776 |
${A_5}$ | 2.3004 | 2.0181 | 2.3004 | 2.7504 | 2.0002 |
${o_1}$ | ${o_2}$ | ${o_3}$ | ${o_4}$ | ${o_5}$ | |
${A_1}$ | 1.5835 | 2.0019 | 2.0816 | 2.2799 | 2.1653 |
${A_2}$ | 1.9971 | 2.5946 | 2.5946 | 1.5835 | 2.4264 |
${A_3}$ | 1.5790 | 1.9912 | 1.8231 | 2.1744 | 2.0909 |
${A_4}$ | 2.1795 | 2.0816 | 1.7771 | 1.7195 | 1.9834 |
${A_5}$ | 2.0816 | 1.8550 | 2.0816 | 2.4359 | 1.8405 |
${o_1}$ | ${o_2}$ | ${o_3}$ | ${o_4}$ | ${o_5}$ | |
${A_1}$ | –4.7589 | –3.7754 | –3.9589 | –3.6844 | –4.3065 |
${A_2}$ | –3.9484 | –2.6632 | –2.6632 | –4.7589 | –2.9811 |
${A_3}$ | –4.9061 | –4.1485 | –4.1021 | –3.6213 | –3.6116 |
${A_4}$ | –3.9306 | –3.9589 | –4.7360 | –4.5020 | –3.8601 |
${A_5}$ | –3.9589 | –4.6080 | –3.9589 | –3.9924 | –4.9499 |