Abstract
Interval type-2 fuzzy set (IT2FS) offers great ability to depict high order information in reality while dealing with both extrinsic and intrinsic facets of uncertainty. In this paper, we try to develop a general framework of information measures and a flexible multiple attribute decision-making (MADM) approach in the interval type-2 fuzzy information environment. First, we propose the fuzzy factor, hesitant factor and interval factor to quantify the fuzziness, hesitancy and interval information of one IT2FS, respectively. An interval type-2 fuzzy cross-entropy has been initiated based on these three factors to measure the discrimination degree of uncertain information between two IT2FSs. Meanwhile, we exploit the axiomatic principles of interval type-2 fuzzy entropy and study the inherent relationship between cross-entropy and entropy measures. Then, some parameterized information measures are naturally investigated, and the decomposition formula suggests that the interval type-2 fuzzy entropy could be expressed as the weighted average of the fuzzy entropy, hesitant entropy and interval entropy. Finally, we construct two programming models based on the maximizing cross-entropy principle to determine attribute weights, and a novel MADM procedure is proposed and applied to a case study on the banks’ liquidity risk evaluation.
Keywords
Introduction
Fuzzy sciences originated almost 50 years ago since Zadeh [33] introduced the concept of fuzzy set (FS) in 1965, and Merigó et al. [21] presented a general overview of research in the fuzzy sciences using bibliometric analysis method. In fuzzy sciences, type-2 fuzzy set (T2FS) theory, initiated by Zadeh [35] in 1975, has been viewed as an effective extension of FS due to its great ability to cope with uncertain nature of real world. In order to clearly characterize the complex structure of a T2FS, Mendel & John [19] proposed a new Representation Theorem, which provides us a convenient way to analysis the uncertain information of T2FSs. Mendel [20] further developed the fundamental theory and extended related operational laws of T2FSs. Nevertheless, we usually confront with high computational overhead when T2FSs are used to solve some practical problems. To resolve these difficulties, an interval type-2 fuzzy set (IT2FS) [18] is expressed as a special case of a general type-2 fuzzy set, which can be just considered the primary membership grade, as the secondary membership degree equals to 1. Because of their easiness and reduced computational effort in comparison with general T2FSs, IT2FSs have received increasing attention both in theory and practice fields.
How to quantify the uncertain information of a T2FS and its extensions is of great importance. Wu & Mendel [30] defined some uncertainty measures, such as centroid, cardinality, fuzziness, variance and skewness, for IT2FSs. Zhai & Mendel [36] extended these measures to general type-2 fuzzy sets. Greenfield and John [12, 13] used vertical slice and invariance under vertical translation to measure uncertainty of T2FSs and IT2FSs, respectively. Some other information measures are also studied in T2FSs, such as similarity [15], inclusion and subsethood measure [27]. Fuzzy entropy, initiated by Zadeh [34] in 1968, has been viewed as one of important tools to provide the foundation for analogical reasoning and could directly describe uncertainty of FSs. Then, fuzzy entropy has been extended to other extensions of FSs. Such as Szmidt & Kacprzyk [26] redefined the constructive axioms of intuitionistic fuzzy entropy, but they cannot express the equilibrium state of supportability and opposability if neutral evidences are indicated in the hesitancy degree. Pal et al. [22] pointed out that existing uncertainty measures cannot capture all fuzziness and lack of knowledge and proposed a generating family of measures. Xu & Xia [31], Zhao et al. [37] introduced some hesitant fuzzy entropies in hesitant fuzzy sets (HFSs). Hwang et al. [15] presented similarity, inclusion and entropy of T2FSs by use of Sugeno integral. However, they only consider the fuzziness of secondary membership degree while neglecting the fuzziness of primary membership grade, and they also have neglected the dispersity of primary membership degree. Up to now, few scholars have extended the entropy concept to IT2FSs and proposed the axiomatic principle of interval type-2 fuzzy entropy. By analyzing the structural characteristics of IT2FSs and on the basis of previous works, we need to point out that the uncertainty of an IT2FS may be divided into three parts: fuzziness, hesitancy and interval information. Based on this thought, we propose a constructive principle of interval type-2 fuzzy entropy.
The classical cross-entropy is used to depict the discrimination information between FSs or their extensions. For instance, Shang & Jiang [25], Vlachos & Sergiadis [28] discussed fuzzy cross-entropy and intuitionistic fuzzy cross-entropy and studied the relationship between these two measures. Xu & Xia [31] introduced a hesitant fuzzy cross-entropy in HFSs. It is worth noting that Mao et al. [17] and Yao et al. [32] proposed two novel kinds of cross-entropies to measure the discrimination degree of uncertain information between two IFSs and T2FSs, respectively. In reality, we need to discriminate and quantify different amount of uncertain information since they usually represent some cost-type facets, such as risk, cost, etc. Inspired by Mao et al. [17] and Yao et al. [32], we try to introduce an interval type-2 fuzzy cross-entropy in IT2FSs and study the inherent relation between cross-entropy and entropy of IT2FSs.
An increased number of multiple attribute or criteria decision-making methods (MADM or MCDM) have been developed in interval type-2 fuzzy environment. Wang et al. [29] studied the multiple attribute group decision-making (MAGDM) models by use of weighted arithmetic averaging operator and ranking-value measure. Some other scholars further extended ranking value methods and applied them in MADM problems, such as Chen & Wang [5], Chen et al. [6], Chen [7], Qin & Liu [23]. Meanwhile, Hu et al. [14] proposed a new approach based on possibility degree to solve multi-criteria decision-making problems in IT2FSs. Qin et al. [24] extended the VIKOR method based on prospect theory to accommodate interval type-2 fuzzy circumstances. Celik et al. [1] gave a comprehensive review of MCDM approaches based on IT2FSs. In addition, some researches have focus on other interval type-2 fuzzy MADM (or MAGDM or MCDM) approaches based on a variety of “traditional” decision-making techniques, such as AHP [4, 16], TOPSIS [3], PROMETHEE [8], ELECTRE [9], Preference relationship [2], and some other methods [10, 11]. However, the decision process in most of these works depends on aggregating operators, and different operators may lead to a contradict results. Meanwhile, the attribute information of these studies is often down to completely known or partly known, few scholars consider the case that the attribute weights are completely unknown. In this paper, we divide the attribute information into two cases: completely unknown and partly known, and construct programming models based on the maximizing cross-entropy principle to determine attribute weights under different cases. We also propose a new expected score function, which considers information comes from both lower membership function (LMF) and upper membership function (UMF), to avoid using aggregating operators. The step-by-step MADM procedure is developed and applied to a case study on banks’ liquidity risk evaluation. Some comparisons with other existing methods are made to verify the validity of our approach.
This paper is organized as follows. In Section 2, we briefly introduce some basic concepts. In Section 3, the interval type-2 fuzzy cross-entropy and entropy are proposed, and the relationship between these two measures is also investigated. We exploit some parameterized information measures for IT2FSs in Section 4. A novel MADM approach and expected score function are presented and applied to a case study concerning banks’ liquidity risk evaluation in Section 5. Section 6 concludes the whole paper.
Preliminaries
In this section, we first give a brief review of some concepts related to T2FSs and IT2FSs, and their basic operations, which will be used in the following sections.
In addition, an IT2FS is completely determined by the union of all primary memberships which is called the footprint of uncertainty (FOU), the FOU of A denoted by FOU (A) may be expressed as follows:
For an IT2FS A, let A
U
(x) and A
L
(x) be the upper membership function (UMF) and the lower membership function (LMF) of A, respectively. That is to say, for any x ∈ X, we have
Therefore, the FOU of an IT2FS A can be also expressed in the following form:
For the sake of simplicity, we denote by A = (A
L
, A
U
) a trapezoidal interval type-2 fuzzy set (TIT2FS) in X, as shown in Fig. 1, where

A trapezoidal interval type-2 fuzzy set A.
Since the trapezoidal interval type-2 fuzzy set is the typical case of a general IT2FS, many scholars have used the trapezoidal IT2FSs rather than the general IT2FSs to analysis the practical problems, such as Qin and Liu [23], Qin et al. [24], Wang et al. [29], Wu and Mendel [30]. Therefore, the following sections are also discussed based on the trapezoidal IT2FSs. For an IT2FS A, whose UMF A
U
(x) and LMF A
L
(x) are defined by
In general, all IT2FSs on X could be denoted by IT2FS (X), and the complement set (denoted by A
U
) of an IT2FS A could be defined as [20]:
Up to now, many kinds of information measures have been constructed to describe uncertainty of fuzzy set and its extensions, such as centroid, cardinality, variance, skewness, similarity, distance and correlation degree, etc. However, among these measures, entropy is considered as one of direct and effective tools to quantify uncertainty. For examples, Shang & Jiang [25] for FSs, Vlachos & Sergiadis [28], Mao et al. [17] for IFSs, Xu & Xia [31], Zhao et al. [37] for HFSs, Takáč [27], Yao et al. [32] for T2FSs. Unfortunately, few scholars have developed these entropy measures in IT2FSs. Inspired by Mao et al. [17] and Yao et al. [32], we try to construct an interval type-2 fuzzy cross-entropy to measure the discrimination degree of uncertainty between two IT2FSs and induce a new interval type-2 fuzzy entropy to quantify uncertainty of one IT2FS.
Some preparations
As we know, different extensions of FS have different types of uncertainty, for example, a FS’s uncertainty is its fuzziness, an IFS’s uncertainty includes its fuzziness and intuitionism, and a T2FS’s uncertainty is consist of fuzziness and hesitancy.
By analyzing characteristic of an IT2FS, we may point out that the uncertain information of an IT2FS could be divided into three facets: fuzziness, hesitancy and interval information. Specifically, fuzziness is reflected by the average closeness between primary membership degree and 1/2 in UMF and LMF, the hesitancy is displayed with dispersity degree of all primary membership degrees in UMF and LMF, and the interval information is determined by area of region between UMF and LMF. In order to quantify these three types of uncertainty in an IT2FS, we propose three new concepts (namely, the fuzzy factor, hesitant factor and interval factor) as follows:
For a given IT2FS A ∈ IT2FS (X), we denote the fuzzy factor, hesitant factor and interval factor by Δ
A
, σ
A
and S
A
, respectively, where
It’s easy to verify that Δ A U , Δ A L , σ A U , σ A L ∈ [0, 1/ - 2] ⇒ Δ A , σ A ∈ [0, 1/ - 2] and S A ∈ [0, 1]. Moreover, we also have Δ A = Δ A C , σ A = σ A C , S A = S A C . It’s worth noting that the fuzziness of A becomes stronger while Δ A decreases, its hesitancy increases with σ A and interval information increases with S A . So the whole uncertainty of an IT2FS could be adequately characterized by a three-tuple (Δ*, σ*, S*).
In reality, it’s of great importance to accurately measure uncertainty of IT2FSs since it often represents cost, risk, etc. However, the traditional cross-entropy is used to measure the discrimination information, which includes determination information and uncertain information. Mao et al. [17] and Yao et al. [32] proposed a novel cross-entropy of IFSs and T2FSs to describe their discrimination of uncertain information, respectively. Inspired by above references, we try to develop an analytical framework to measure the discrimination degree of uncertain information between two IT2FSs based on three-tuple (Δ*, σ*, S*).
However, one can observe that CE (A, B) is not symmetric with respect to its arguments. Then a symmetric form is given as
DE (A, B) = DE (B, A); DE (A, B) = DE (A
C
, B) = DE (A, B
C
) = DE (A
C
, B
C
); 0 ≤ DE (A, B) ≤3 ln 2.
(3) Since Δ A , σ A ∈ [0, 1/ - 2] and S A ∈ [0, 1], by use of Shannon’s inequality, we have CE (A, B) ≥0 ⇒ DE (A, B) ≥0. Meanwhile, the relation CE (A, B) ≤ (Δ A + σ A + S A ) ln 2 ≤ 3 ln 2/ - 2 is true, so we have DE (A, B) ≤3 ln 2. □
Since Zadeh [34] proposed fuzzy entropy in 1968 to describe FS’s fuzziness, this concept has gradually become one of powerful tools to measure uncertainty while its has been extended to other extensions of FSs both in its form and connotation, such as Szmidt & Kacprzk [26], Pal et al. [22], Mao et al. [17] developed intuitionistic fuzzy entropy, Xu & Xia [31], Zhao et al. [37] defined hesitant fuzzy entropy, and Hwang et al. [15], Takáč [27] have given some basic axioms of entropy for T2FSs. Based on these previous works, we now propose some axiomatic principles of an interval type-2 fuzzy entropy as follows.
E (A) =0 iff (if and only if) A is a crisp set, i.e., Δ
A
= 1/ - 2, σ
A
= S
A
= 0; E (A) ≤1 for any A ∈ IT2FS (X); E (A) = E (A
C
); f (Δ
A
, σ
A
, S
A
) is a real-valued continuous function, which decreases with the first variable Δ
A
and increases with second variable σ
A
and third variable S
A
.
The classical relationship between entropy and cross-entropy have been studied by Vlaschos & Sergiadis [28], Xu & Xia [31], that is, E (A) = θ × DE (A, A C ), where θ is a normalized parameter. However, this relation may be false in our models due to DE (A, A C ) =0 for any A ∈ IT2FS (X). Motivated by Mao et al. [17], we denote a crisp set by C*, i.e., Δ C * = 1/ - 2, σ C * = 0, S C * = 0 then DE (A, C*) could be viewed as a method to measure the amount of uncertainty of an IT2FS A, which agrees with the intuitive concept of interval type-2 fuzzy entropy.
(2) According to Theorem 3.1, we know E (A) ≤1 and E (A) = E (A C ).
(3) In order to prove E (A) satisfies the principle (4) in Definition 3.2, it suffices to prove that the following function:
Let x = Δ A , y = σ A and z = S A , then E (A) = g (Δ A , σ A , S A )/ - 3 ln 2 decreases with Δ A and increases as σ A and S A . □
Vlachos & Sergiadis [28] had proposed an unified formulation to analysis the relation of entropy, fuzziness and intuitionism for IFSs. Mao et al. [17] also gave a decomposition formula of intuitionistic fuzzy entropy. Based on these works, we try to construct a similar decomposition formula of interval type-2 fuzzy entropy to study inherent connections of fuzziness, hesitancy and interval information. In order to clearly illustrate this relationship, we introduce three kinds of entropies to discuss different effects of each factor (the fuzzy factor, hesitant factor and interval factor).
E
F
(A) =0 iff Δ
A
= 1/ - 2; E
F
(A) =1 iff Δ
A
= 0, namely, A
U
(x) = A
L
(x) = 1/ - 2, ∀ x ∈ X; E
F
(A) = E
F
(A
C
); f (Δ
A
) decreases with the fuzzy factor Δ
A
.
E
H
(A) =0 iff σ
A
= 0, namely A
U
(x) and A
L
(x) are constant; E
H
(A) =1 iff σ
A
= 1/ - 2; E
H
(A) = E
H
(A
C
); f (σ
A
) increases with the hesitant factor σ
A
.
E
I
(A) =0 iff S
A
= 0, namely, A
U
(x) = A
L
(x) , ∀ x ∈ X; E
I
(A) =1 iff A
U
(x) =1, A
L
(x) =0; E
I
(A) = E
I
(A
C
); f (S
A
) increases with the hesitant factor S
A
.
According to the above definitions, fuzzy entropy E F (A) measures uncertain information caused by fuzziness of A, hesitant entropy E H (A) describes uncertainty induced by its hesitancy, and interval entropy E I (A) depicts the uncertainty arise from its interval information. If we rewrite Equation (6) as
By Definitions 3.3, 3.4 and 3.5, it’s easy to verify that E F (A), E H (A) and E I (A) are fuzzy entropy, hesitant entropy and interval entropy of A, respectively. From Equation (7), we may conclude that an interval type-2 fuzzy entropy E (A) could be viewed as the arithmetic average of fuzzy E F (A), hesitant entropy E H (A) and interval entropy E I (A).
It’s evident to find that the coefficients of E F (A), E H (A) and E I (A) are all 1/3 in Equation (7). This indicates that E F (A), E H (A) and E I (A) have same effect on E (A) in the view of Mathematics, which implies that fuzziness, hesitancy and interval information have the same position to the global uncertainty of one IT2FS. Naturally, if we consider these three facets have different effects on the whole uncertainty, we need introduce some parameters to adjust and control them.
Obviously, CEp,q,r (A, B) is also not symmetric with its arguments, so a symmetric parameterized type-2 fuzzy cross-entropy could be obtained in similar way:
DEp,q,r (A, B) = DEp,q,r (B, A); DEp,q,r (A, B) = DEp,q,r (A
C
, B) = DEp,q,r (A, B
C
) = DEp,q,r (A
C
, B
C
); 0 ≤ DEp,q,r (A, B) ≤ (p + q + r) ln 2.
In similar way, the parameterized interval type-2 fuzzy cross-entropy may induce a parameterized interval type-2 fuzzy entropy as following theorem.
Obviously, Ep,q,r (A) fulfils the principle (4) if we are able to show that the function
Let x = Δ A , y = σ A and z = S A , Ep,q,r (A) = g (Δ A , σ A , S A )/ - [(p + q + r) ln 2] decreases with Δ A and increases with σ A and S A . □
We can also divide Equation (10) into three parts:
We could verify that
In general, if p > q > r, it shows that fuzziness of A has a stronger effect on the whole uncertainty than its hesitancy and interval information, and other situations may also be explained in the similar way. However, we should choose the consistent parameters p, q and r in the same practical problem, and make comparison for different events under a fixed entropy distribution.
In this section, we solve multiple attribute decision-making (MADM) problems under interval type-2 fuzzy environment by use of above information measures. We suppose there are m alternatives x
i
(1 ≤ i ≤ m) and n attributes e
j
(1 ≤ j ≤ n) with the attribute weight vector
A new MADM approach based on information measures
Based on above notations, a new MADM method by use of the proposed information measures could be described in the following steps, and the procedure of the MADM approach can be shown in Fig. 2.

The procedure of the new MADM approach.
By means of Lagrange multiplier method and normalization, we have
Similarly, if we use the parameterized cross-entropy do construct the model
The model
In order to illustrate the feasibility and reasonability of the proposed MADM method, we shall give a practical example about liquidity risk evaluation for banks. We suppose the decision maker expect to form linguistic terms (Table 1) to give the linguistic value to evaluate liquidity risk with interval type-2 fuzzy information. Table 1 shows the linguistic terms “Very Low” (VL), “Low” (L), “Medium Low” (ML), “Medium” (M), “Medium High” (MH), “High”’ (H), “Very High” (VH) and their corresponding interval type-2 fuzzy sets, respectively.
Linguistic terms and their corresponding IT2FSs
Linguistic terms and their corresponding IT2FSs
The financial crisis of 2007 that triggered the turmoil in the financial markets has demonstrated the central role that effective processes for managing liquidity risk play in maintaining both the stability of individual banks and the soundness of the entire banking system. In the narrow sense, liquidity risk is the potential inability of a bank to meet punctually and in a cost-effective way its envisaged contractual payment obligations when they fall due. Thus the liquidity for the bank takes on an essentially protective meaning that is able to continuously ensure a state of technical solvency, and liquidity risk evaluation is of great importance for banks all over the world.
Here, we suppose there are five banks (x1, x2, ⋯ , x5) will be evaluated and four attributes will be considered in reality: the ability to ensure at all times an adequate corresponding balance between cash inflows and cash outflows (e1), the ability to coordinate the issuing by the bank of short, medium and long term financing instruments (e2), the ability to optimize the costs of refinancing, striking a trade-off balance between liquidity and profitability (e3), and the ability to optimize, for banks structured as banking groups, the intra-group management of cash flows, with the aim of reducing dependence on external financial requirements, by means of cash pooling techniques or other optimization instruments (e4). The evaluation information on the five banks with respect to each attribute is characterized by IT2FSs, which are contained in the interval type-2 fuzzy original matrix, as shown in Table 2.
The evaluation information on five banks
The realization of decision-making steps is shown as follows:
Attribute weights under different situation
If the attribute weights are partly known (Case 2), and assume the information set of attribute weightsis ω ∈ Φ = {0.2 ≤ ω1 ≤ 0.3, 0.1 ≤ ω2 ≤ 0.2, 0.2 ≤ ω3 ≤ 0.4 0.2 ≤ ω4 ≤ 0.4}, we may obtain the weights by solving the model
The collective scores under different entropies
It is worth noting that the bigger the score of a certain bank, the lower the corresponding liquidity risk. Therefore, we find that the ranking results are consistent while using Equations (14 and 15), and the bank 2 has the lowest liquidity risk and the bank 4 has the highest liquidity risk, that is to say bank 2 has the strongest solvency ability and stability. If we impose some information on attribute weights, we find the overall conclusion is changeless while the position between bank 1 and bank 5 is reversed, which implies that the weights information will play a role in ranking results.
In order to verify the validity of the proposed MADM approach, we make some comparisons with other previous methods, such as ranking value method [6, 23], possibility degree method [14], by calculating the above example. In order to keep the consistency of weights information, we let weight vector ω = {0.3, 0.1, 0.4, 0.2} in these previous methods.
(1) According to the ranking value formula proposed by Chen et al. [6], and we may establish the ranking matrix and calculate the average agreement (AD), then the comprehensive ranking value of each bank is shown as follows:
Thus the ranking order of all banks is x2 ≻ x1 ≻ x5 ≻ x3 ≻ x4.
(2) By use of combined ranking value formula proposed by Qin & Liu [23] and IT2OWA operator, we have
So all banks are ranked as x2 ≻ x5 ≻ x1 ≻ x3 ≻ x4.
(3) Based on the possibility degree proposed by Hu et al. [14] and TIT2-WAA operator, we can obtain the possibility degree matrix P as follows:
Then we obtain the ranking vector r = (0.1008, 0.1037, 0.0985, 0.0948, 0.1022), so the ranking result is x2 ≻ x5 ≻ x1 ≻ x3 ≻ x4.
All ranking results including both he proposed method and previous approaches could be shown in Table 5.
Comparisons of four MADM methods
Comparisons of four MADM methods
From Table 5, it is clear that most of four methods have the same results, this verifies the method we proposed is reasonable and validity in this paper. We can also find some advantages while comparing with other previous methods as follows: The method of determining attribute weights is more reasonable. We construct two programming models based on the maximizing cross-entropy principle, which divides the weights information into two cases: completely unknown and partly known. However, all weights in Chen et al. [6] are assumed equal and not consider the importance of attribute information. Qin & Liu [23] and Hu et al. [14] only consider the situation that weights are partly known. Meanwhile, the maximizing cross-entropy method considers simultaneously at least three aspect of information of both LMF and UMF: fuzziness, hesitancy and interval information. Comparing with these three methods, the flexible of our approach is much higher than previous methods, and our method has a solid theoretical foundation. The proposed MADM method does not depend on the aggregating operators. In fact, different aggregating operators may loss different levels of information, which might lead to some contradict results. Such as Chen et al. [6] adopted arithmetic operations, Qin & Liu [23] used IT2OWA operator and Hu et al. [14] considered TIT2-WAA operator. However, our method transforms the decision matrix to a score matrix and avoids using aggregating operators, which changes the complex interval type-2 fuzzy information into the simple real numbers. Comparing with existing methods, the computational complexity of our approach is much lower and the decision-making process is more straightforward.
According to the above stated comparisons and analysis, our proposed MADM approach is better than other three methods.
Interval type-2 fuzzy sets (IT2FSs) have received increasing attention due to its great ability to deal with imprecise and ambiguous information in real word. The purpose of this paper is to develop a general framework of interval type-2 fuzzy information measures and propose a new MADM approach in interval type-2 fuzzy information environment. The main results could be concluded as follows:
Firstly, we point out, for the first time, that uncertainty of an IT2FS is composed of its fuzziness, hesitancy and interval information, which could be quantified by fuzzy factor, hesitant factor and interval factor. Thus, we initiate an interval type-2 fuzzy cross-entropy and entropy measures to describe the discrimination degree of uncertainty between two IT2FSs and uncertainty level of one IT2FS, respectively. We also study the inherent connection between two kinds of entropies. Meanwhile, we extend these information measures to parameterized case, and find interval type-2 fuzzy entropy could be viewed as the weighted average of fuzzy entropy, hesitant entropy and interval entropy.
Secondly, we construct two optimal models based on the maximizing cross-entropy principle to determine the attribute weights, and develop a novel MADM method in interval type-2 fuzzy information environment by using the proposed information measures. Meanwhile, a case study concerning banks’ liquidity risk evaluation is provided to illustrate the applicability and flexibility of the above approaches.
In future, we shall extend the information measures to general type-2 fuzzy environment, with expectations that could be made applicable to other similar decision-making problems, such as supplier selection, public transportation evaluation and social network analysis and so on.
Footnotes
Acknowledgments
The authors are highly grateful to any anonymous for their careful reading and insightful comments, and the views and opinions expressed are those of the authors. The work is supported by the National Nature Science Foundation of China (No. 71473036), the Talent Introduction Project of Anhui University and the Natural Science Key Project of Anhui Sanlian University (No. kjzd 2016 001).
