Next Article in Journal
A Functional Data Analysis Approach for the Detection of Air Pollution Episodes and Outliers: A Case Study in Dublin, Ireland
Previous Article in Journal
Multi-Attribute Group Decision Making Based on Multigranulation Probabilistic Models with Interval-Valued Neutrosophic Information
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

New Framework for Quality Function Deployment Using Linguistic Z-Numbers

1
School of Business, Central South University, Changsha 410012, China
2
Business School, Guilin University of Technology, Guilin 541004, China
*
Author to whom correspondence should be addressed.
Mathematics 2020, 8(2), 224; https://doi.org/10.3390/math8020224
Submission received: 24 December 2019 / Revised: 4 February 2020 / Accepted: 6 February 2020 / Published: 10 February 2020
(This article belongs to the Section Fuzzy Sets, Systems and Decision Making)

Abstract

:
Quality function deployment (QFD) is a useful design quality control tool in service enterprises and manufacturing enterprises. However, there are several issues in extant QFD frameworks, that is, in three aspects: description of evaluation information, weight determination of expert team members (TMs), and weight identification of customer requirements (CRs). In order to address these issues, a novel QFD framework is first proposed utilizing linguistic Z-numbers (LZNs) with integrated subjective and objective weights of TMs and CRs. The LZNs can represent uncertain information and the reliability of information in a specific way while the fuzzy numbers cannot. Moreover, the order relation analysis (G1) method and improved maximum consensus (MC) method are developed to get the subjective and objective weights of TMs, respectively. Further, the step-wise weight assessment ratio analysis (SWARA) method and statistical distance (SD) method are studied to acquire combined weights of CRs. Next, the proposed QFD framework is applied to a case of logistics service provider, which illustrates the availability and utility of the framework. Then, a sensitivity analysis is conducted to prove the reliability of the framework. Finally, two comparative analyses are performed to declare the advantages of the framework. Results indicate the proposed QFD framework is better than existing models.

1. Introduction

As one of the most famous design quality control technologies, quality function deployment (QFD) [1] is an effective and efficient way to convert customer requirements into production demand [2]. This facilitates research and development (R&D) companies to design and manufacture customer-driven products. Generally speaking, there are four phases in classical QFD models, including product overall design phase, part design phase, process design phase, and production system design phase [3]. In these phases, customer requirements are translated into engineering characteristics, part characteristics, process plans, and production requirements, respectively [4]. Further, a matrix is required in each phase to construct these relationships which is called as house of quality (HOQ).
The HOQ in the first phase is the most fundamental and central one [5] which is illustrated in Figure 1. In particular, the seven components of HOQ are customer requirements (CRs), engineering characteristics (ECs), relationship between CRs and ECs, importance of CRs, importance of ECs, correlation among ECs, and benchmarking analysis [4]. As far as we know, numerous researchers have integrated QFD models with multi-criteria decision-making (MCDM) theories to rank the ECs [6,7]. However, although the usefulness of QFD models has improved a little, three deficiencies still exist in these studies which are formulated below: (1) the evaluation information of relationships between CRs and ECs was given with crisp values or various kinds of fuzzy numbers [8]. Nevertheless, in most actual circumstances, the information provided by experts may be unreliable or a bit reliable because of their lacks of professional knowledge and limited time to make decisions; (2) as a rule, the ranking process of ECs requires several QFD experts to participate in the consultations. However, the weights of expert team members (TMs) were provided directly without any discussions [2], which was not scientific and rational. What’s more, it was inconsistent with actuality to determine the experts’ weights by an either subjective or objective method [9]; (3) the ranking results may be inaccurate in real-life decision-making processes when only subjective or objective weighting method employed [10,11].
To address the above mentioned deficiencies, a new framework is essential for QFD models. At that point, the motivations of this paper are itemized as follows:
(1)
Considering the fuzzy information and reliability of this information at the same time, Z-numbers are suitable tools to describe the evaluation values [12,13]. Nonetheless, the operations of Z-numbers are complex [14,15]. Meanwhile, the processing of Z-numbers into various fuzzy numbers may lose the original information. Therefore, in this paper, linguistic Z-numbers (LZNs) [16,17] are utilized to represent the evaluation information.
(2)
To dispose the improper ways to determine the weights of TMs with single preference, the proposed QFD model takes the subjective and objective weights into considerations. Since order relation analysis (G1) method [18] and maximum consensus (MC) method [19] are practical means to identify the subjective and objective weights, separately. Further, MC method is improved due to the issue that it merely thinks over the positive ideal point. Thus, the G1 method and improved MC method are developed to obtain combined weights of TMs.
(3)
In order to fill the research gap about the combined weights of TMs, SWARA method [20] and extended statistical distance (SD) method [21] are mixed to acquire the hybrid weights of CRs.
This paper aims to propose a novel framework for QFD using LZNs with integrated weights of TMs and CRs. This framework contains three stages. The first stage is to construct the decision matrix, the second is to determine the combined weights of TMs and CRs, and the third is to rank ECs in terms of their importance. In sum, the contributions of this study are listed below:
(1)
The novel framework for QFD is set up. The sensitivity analysis and comparative analyses indicate the reliabilities, effectiveness, and practicality of the proposed QFD model. Namely, the proposed QFD framework is superior to existing QFD models.
(2)
LZNs are absorbed to express the evaluation information, which involves the fuzzy information and the reliability of fuzzy information simultaneously.
(3)
The G1 method and improved MC method are utilized to get the combined weights of experts. And SWARA method and extended SD method are used to gain the combined weights of CRs. These methods are appropriate for QFD to be applied in real life.
(4)
The MC approach is improved by considering both the positive and negative ideal points. The weights derived by improved MC method may be more precise than traditional MC method.
(5)
The technique for order performance by similarity to ideal solution (TOPSIS) method is extended with LZNs which enriches the theories of LZNs.
The rest of this study is established in the following order. In Section 2, the literatures about QFD models are reviews. Section 3 introduces the definitions of Z-numbers and linguistic terms sets. Moreover, the relevant theories of LZNs are presented. In Section 4, the concrete QFD framework is proposed. Section 5 implements the proposed QFD framework in a case to verify the availability and validity of it. Finally, the conclusions are drawn in Section 6.

2. Literature Review

With the development of the social economy and enlarged volume of information, the application of multiple criteria decision-making (MCDM) has become increasingly wide. Many popular and topical MCDM methods have been proposed by researchers, such as: SAW (simple additive weighting method), AHP (analytic hierarchy process) [22], TOPSIS (technique for order of preference by similarity to ideal solution) [23], ELECTRE (elimination and choice expressing reality) [24], VIKOR (VlseKriterijumska Optimizacija I Kompromisno Resenje) [25], PROMETHEE (preference ranking organization method for enrichment evaluations) [26], GDM (grey decision-making) [27], SWARA [20], and BWM [28].
The QFD model, integrating with fuzzy theories and MCDM methods, has been widely investigated by many scholars [29]. For instance, analytical network process (ANP) and fuzzy logic were utilized to structure the intensive QFD framework by Zaim et al. [30]. Abdel-Basset, Manogaran [31] presented a QFD method where analytical hierarchy process (AHP) was adopted to determine the weights of CRs. Hsu et al. [29] developed a new QFD approach with fuzzy Delphi method. Wu et al. [32] proposed a two-phase robust fuzzy QFD procedure to assess the service station for nail catering service. Wu et al. [3] provided a hybrid QFD model with decision-making trial and evaluation laboratory (DEMATEL) technique and VIKOR method within hesitant fuzzy environment. In addition, the fuzzy QFD-MCDM frameworks were applied to various fields, such as supply chain management strategies formulation problems [2,33], risk management of hazardous material transportation process [34], and green supplier selection problems [6].
However, as indicated above, the ways to determine weights of QFD experts and customer requirements were inaccurate in extant studies. That is, the weights of experts and CRs were either supposed to be known completely or calculated by only one type of approach: subjective approach or objective approach. To be specific, in subjective weight estimation methods, the weights information is obtained based on humans’ subjective judgments or experience [35,36]. In contrast with the subjective methods, the objective weights are derived merely by the data given in the decision matrices [37,38]. Therefore, it may cause the intense subjectivity or the absence of subjective preference if only the subjective method or objective method is used in weights determination.
There was an issue in the information representation of extant fuzzy QFD-MCDM models, i.e., the models cannot reserve the irresponsible evaluation integrally. Hence, Z-numbers are the pretty helpful instruments to denote this information [12]. Two ways were investigated to handle Z-numbers. One is to dispose Z-numbers straightforward with the operations of them [15]. However, Wang et al. [16] pointed out that the operations are too intricate to be employed in common life. The other is to transform Z-numbers into other forms of fuzzy numbers. Nevertheless, Wang et al. [16] indicated that this way may lead to information distortion. As a result, Wang et al. [16] defined the linguistic Z-numbers (LZNs) to cover these problems.
In this paper, the evaluation information of relationships between CRs and ECs is depicted by LZNs so that the information can be kept perfectly. What’s more, both subjective and objective methods are proposed to derive the combined weights of QFD experts and CRs. The outcomes show that the proposed QFD framework outperforms state-of -the-art QFD methods.

3. Preliminaries

In this section, we review some basic notions related to linguistic Z-numbers and the operational law, distance, and comparison method of linguistic Z-numbers are introduced. Since linguistic Z-numbers can be an integrated formation of Z-numbers and linguistic terms sets [16], we briefly present the definitions of Z-numbers and linguistic terms sets.
Definition 1
([12]).Let X be a real-valued uncertain variable. A Z-number ( A , B ) is an ordered pair of fuzzy numbers which is related to X , where A is a fuzzy restriction on the values that X is allowed to take, and B is a reliability measure of A .
In general, the two components A and B   are described in nature languages. For example, the Z-number ( A , B ) = (about 45 million, quite sure) can be used to evaluate the population of Spain ( X ).
Definition 2
([39]).Let S = { s a | a = 0 , 1 , , 2 t } be a finite ordered discrete term set, where S a represents a possible value for a linguistic variable. We require that s a and s b satisfy the following conditions:
(1)
The set S is ordered: s a > s b if and only if a > b ;
(2)
Negation operator: neg ( s a ) = s 2 t a .
Definition 3
([16]).Let X be a universe of discourse, S 1 = { s 0 , s 1 , , s 2 p } and S 2 = { s 0 , s 1 , , s 2 q } be two linguistic term sets. A linguistic Z-number set (LZNS) in X is composed of the fuzzy restriction A ( x ) and the reliability measure B ( x ) , where A ( x ) S 1 and B ( x ) S 2 . Then, a LZNS can be defined as follows:
Z = { x , A ( x ) , B ( x ) | x X }
It should be noted that two LTSs, S 1 and S 2 , denote different linguistic preference information, separately. In particular, a LZNS is regarded as a linguistic Z-number (LZN) when X has only one element which is illustrated by Z = x , A ( x ) , B ( x ) . In order to facilitate the representation, a LZN can be also denoted by Z = A , B , where A and B are linguistic terms, respectively.
Definition 4
([40]).Let S = { s i | i = 0 , 1 , , 2 t } be a linguistic term set, and θ i is a numeric value ranging from 0 to 1. The linguistic scale function (LSF) conducts the mapping from s i to θ i which is defined below:
f : s i θ i ( i = 0 , 1 , , 2 t )
where 0 θ 1 θ 2 θ 2 t 1 . A LSF is a strictly monotonically increasing function regarding to label i . Simultaneously, symbol θ i expresses the preference information of decision-makers (DMs) when linguistic term s i is used. Thus, LSF represents the semantics of s i . There are three types of LSFs are shown as follows:
(1)
The first LSF is defined in terms of subscript function s u b ( s i ) = i . Further, the appraisement scale for linguistic term is even.
f 1 ( s i ) = θ i = i 2 t ( i = 0 , 1 , , 2 t )
(2)
The second LSF is developed according to the exponential scale. Moreover, from the middle of given linguistic term to both ends, the absolute deviation between adjacent linguistic terms increases.
f 2 ( s i ) = θ i = { γ t γ t i 2 γ t 2 ( i = 0 , 1 , , t ) γ t + γ i t 2 2 γ t 2 ( i = t + 1 , t + 2 , , 2 t )
As we can see in Equation (4), there is a parameter γ which can be identified by experts’ advice or subjective approach. On the basis of lots of experiment researches, Bao, Lian [41] explained that the value γ usually is derived in the interval [ 1.36 , 1.4 ] .
(3)
The third LSF is proposed based on prospect theory’s value function. Prospect theory [42], is a very successful descriptive model about decision making under risk. This theory acknowledges value gains and losses differently, placing more weight on perceived gains versus perceived losses. This theory is more consistent with the practical decision, and it can explain a lot of things that the expected utility theory can’t. From the middle of given linguistic term to both ends, the absolute deviation between adjacent linguistic terms decreases.
f 3 ( s i ) = θ i = { t α ( t i ) α 2 t α ( i = 0 , 1 , , t ) t β + ( i t ) β 2 t β ( i = t + 1 , t + 2 , , 2 t )
where α , β 0 , 1 are the concavo-convex degree of the subjective value function in the area of the gain and loss, respectively. On the basis of empirical data, Kahneman and Tversky [42] gave α = β = 0.88 .
For information integrity and calculation convenience in the process of decision-making, these three aforementioned functions are supposed to be continuous, that is, f * : s i θ i ( 0 i 2 t ) . Further, the inverse function of f * is marked as f * 1 .
There are also some LSFs based on other distances could be considered, such as fuzzy divergence [43] and fuzzy similarity [44].
In this study, we assume γ = 1.4 ,   α = β = 0.88 , and the characteristics of three LSFs are displayed in Figure 2.
Definition 5
([16]).Let z i = ( A i , B i ) and z j = ( A j , B j ) be two LZNs, and f * and g * be one of three LSFs denoted by Equations (1) and (2), separately. The operational laws of LZNs are defined as follows:
(1)
n e g ( z i ) = ( f * 1 ( f * ( A 2 p ) f * ( A i ) ) , g * 1 ( g * ( B 2 q ) g * ( B i ) ) ) ;
(2)
z i z j = ( f * 1 ( f * ( A i ) + f * ( A j ) ) , g * 1 ( f * ( A i ) × g * ( B i ) + f * ( A j ) × g * ( B j ) f * ( A i ) + f * ( A j ) ) ) ;
(3)
z i z j = ( f * 1 ( f * ( A i ) f * ( A j ) ) , g * 1 ( g * ( B i ) g * ( B j ) ) ) ;
(4)
λ z i = ( f * 1 ( λ f * ( A i ) ) , B i ) , where λ 0 ;
(5)
z i λ = ( f * 1 ( f * ( A i ) λ ) , g * 1 ( g * ( B i ) λ ) ) , where λ 0 .
Definition 6
([16]).Let z i = ( A i , B i ) and z j = ( A j , B j ) be two LZNs, the score function of z i is defined below:
S F ( z i ) = f * ( A i ) × g * ( B i )
The accuracy function of z i is defined as follows:
A F ( z i ) = f * ( A i ) × ( 1 g * ( B i ) )
Therefore, the comparison approach can be proposed by the following rules:
(1)
If A i > A j and B i > B j , then z i is strictly superior than z j , denoted by z i > z j ;
(2)
If S F ( z i ) > S F ( z j ) or S F ( z i ) = S F ( z j ) and A F ( z i ) > A F ( z j ) , then z i is superior than z j , denoted by z i > z j .
(3)
If S F ( z i ) = S F ( z j ) and A F ( z i ) = A F ( z j ) , then z i is equal to z j , denoted by z i z j .
(4)
If S F ( z i ) = S F ( z j ) and A F ( z i ) < A F ( z j ) or S F ( z i ) < S F ( z j ) , then z i is inferior than z j , denoted by z i < z j .
Definition 7
([16]).Let z i = ( A i , B i ) and z j = ( A j , B j ) be two LZNs, and f * and g * be two LSFs. The distance between z i and z j is defined as follows:
d ( z i , z j ) = 1 2 ( | f * ( A i ) × g * ( B i ) f * ( A j ) × g * ( B j ) | + max { | f * ( A i ) f * ( A j ) | , | g * ( B i ) g * ( B j ) | } ) .

4. The Proposed Model for QFD

This section develops a hybrid decision-making model to obtain the importance degrees of engineering characteristics (ECs). In this model, both subjective and objective weights are considered for QFD team members (TMs) and customer requirements (CRs) at the same time. There are three stages included in the proposed model. The first stage is to determine the subjective weights of QFD experts TMs by order relation analysis (G1) method and objective weights of TMs by improved maximum consensus (MC) method, and the second stage is to estimate the subjective weights of CRs by SWARA method and objective weights of CRs by statistical distance (SD) method. In addition, the last stage is to prioritize the ECs according to the importance ratings of them. For convenience, Figure 3 is provided to visualize the decision-making procedures in the developed model.
In order to facilitate the description of the proposed model, we assume that there are m ECs ( E C i , i = 1 , 2 , , m ) with respect to n CRs ( C R j , j = 1 , 2 , , n ) to be evaluated by h QFD experts TMs ( T M k , k = 1 , 2 , , h ) . The subjective and objective weight vectors of TMs are λ = ( λ 1 , λ 2 , , λ h ) and λ = ( λ 1 , λ 2 , , λ h ) , respectively. Further, the subjective and objective weight vectors of CRs are w = ( w 1 , w 2 , , w n ) and w = ( w 1 , w 2 , , w n ) , separately. Moreover, the equalities k = 1 h λ = k = 1 h λ = 1 and j = 1 n w = j = 1 n w = 1 . What’s more, the combined weights of TMs and CRs are calculated by the following equations:
λ k = τ λ k + ( 1 τ ) λ k ,   for   k = 1 , 2 , , h ,
w j = ξ w j + ( 1 ξ ) w j ,   for   j = 1 , 2 , , n ,
where 0 < τ , ξ < 1 , λ k and w j are the integrated weights of T M k and C R j , respectively.
Furthermore, let U k = [ u i j k ] m × n ( k = 1 , 2 , , h ) be the evaluation matrix given by T M k where u i j k is assessed value of E C i with regard to C R j . In addition, assume V = [ v i j ] m × n be the collective decision matrix in which v i j = u i j k λ k .

4.1. Identifying the Integrated Weights of TMs

In this part, the order relation analysis (G1) method and improved maximum consensus (MC) method are utilized to capture the integrated subjective and objective weights of QFD team members.

4.1.1. Determining the Subjective Weights of TMs

A deal of subjective weights determining methods has conducted by many scholars. For example, AHP method was employed to derive the criteria weights in the case of tourism environmental impact assessment [35], What’s more, Wang, Zhang [45] used multi-criteria group decision-making method to hotel building energy efficiency retrofit project selection. However, all of these subjective weights determining methods have lots of complexity and inconvenience. Thus, the G1 method was proposed by Ruan and Yang [46] as a succinct method to get the subjective weights. Due to G1 method’s characteristics, it was applied by Xu, Liu [18] to identify criteria weight vector. In this paper, we utilize G1 method to figure out the weights of TMs subjectively. The detailed steps are shown below.
Step O1. Rank the TMs on the basis of their importance.
The ranking order of TMs, which is denoted by a set T M = { T M 1 , , T M k , , T M h } , should be given decreasingly by the leader of them according to their importance.
Step O2. Identify the relative importance degree between adjacent TMs.
Based on a description of the relative importance degree as shown in Table 1, the relative importance degree θ k = λ k 1 λ k of the QFD expert T M k 1 and T M k can be determined for k = 2 , 3 , , h .
Step O3. Compute the subjective weights of TMs.
The subjective weights of TMs can be derived by the following equations:
λ h = ( 1 + l = 2 h k = l h θ k ) 1
λ k = l = k + 1 h θ l λ h

4.1.2. Determining the Objective Weights of TMs

There is a mass of subjectivity when obtaining the weights of QFD team members with only subjective method. Therefore, to avoid much subjectivity, it is quite necessary to derive the weights of TMs as a supplement by an objective method. Inspired by the maximum consensus (MC) method in the current literature [19,47], this study proposes an improved MC method which takes both positive and negative ideal assessment into account simultaneously, by contrast, the traditional MC method in related literatures only considers the positive ideal point.
Step M1. Determine the positive and negative ideal assessment matrices.
The average of all QFD team members’ evaluation can be deemed as the positive ideal assessment, and the positive ideal assessment matrix is denoted by the following formula:
U + = [ u i j + ] m × n = [ 1 h k = 1 h u i j k ] m × n
There are two negative ideal assessment matrices. One is composed with the maximum evaluation value of an EC with respect to a certain CR, and the other consists of the minimum assessment value of an EC regarding a certain CR. Moreover, these two negative ideal assessment matrices can be represented as follows:
U m a x = [ u i j m a x ] m × n = [ m a x k { u i j k } ] m × n
U m i n = [ u i j m i n ] m × n = [ m i n k { u i j k } ] m × n
Step M2. Calculate the distances between every TM’s evaluation and ideal assessments.
Based on the distance defined in Definition 7, the differences between personal evaluation and ideal assessments are shown below:
d i f f ( U k , U + ) = 1 n × m i = 1 m j = 1 n d ( u i j k , u i j + )
d i f f ( U k , U m a x ) = 1 n × m i = 1 m j = 1 n d ( u i j k , u i j m a x )
d i f f ( U k , U m i n ) = 1 n × m i = 1 m j = 1 n d ( u i j k , u i j m i n )
Step M3. Compute the objective weights of QFD team members.
Obviously, the bigger weights can be assigned if the TM’s evaluation is closer to positive ideal assessment and farther from any negative ideal assessments. Thus, the objective weights of QFD experts TMs can be identified by the expression below:
λ k = 1 d i f f ( U k , U + ) + min { d i f f ( U k , U m a x ) , d i f f ( U k , U m i n ) } k = 1 h ( 1 d i f f ( U k , U + ) ) + k = 1 h min { d i f f ( U k , U m a x ) , d i f f ( U k , U m i n ) }
where k = 1 , 2 , , h .

4.2. Identifying the Integrated Weights of CRs

This subsection develops SWARA method and statistical distance (SD) method to obtain subjective and objective weights of CRs, separately. And the integrated weights of CRs are acquired ultimately.

4.2.1. Determining the Subjective Weights of CRs

To the best of our knowledge, the SWARA method is adopted widely to derive the weight vector [20,48]. For example, the SWARA method was employed to prospect for polymetallic mineralization [49]. Further, Nakhaei, Bitarafan [50] utilized the SWARA method to assess the weights of indexes and sub-indexes. Due to SWARA method’s practical applicability, this study constructs subjective weights of CRs with this utility method.
Step S1. Rank the CRs according to their importance.
Through QFD experts team members’ discussions, the order of CRs concerning importance degree can be provided and the most important CR is listed at the first place. The order is represented by C R = { C R 1 , , C R j , , C R n } .
Step S2. Calculate the weighted average comparative importance value of adjacent CRs.
Each QFD expert T M k shall give a comparative importance value of adjacent CRs s j k for k = 1 , 2 , , h and j = 2 , 3 , , n . Moreover, the symbol s j = λ k s j k denotes weighted average comparative importance value in which λ k is the combined weights of TMs.
Step S3. Determine characteristics of the comparative importance and recalculated weights.
The characteristic of the comparative importance x j and recalculated weight y j for j = 1 , 2 , , n can be computed by Equations (17) and (18)
x j = s j + 1 ,
y j = y j 1 x j .
Step S4. Determine the final weight.
On the basis of the recalculated weight, the final subjective weight w j for j = 1 , 2 , , n can be calculated by the following formula:
w j = y j j = 1 n y j .

4.2.2. Determining the Objective Weights of CRs

It is not sufficient to only consider the QFD experts’ subjective preference on the weights of CRs. There are many objective weights determination methods, such as maximizing deviation method [51] and entropy weight method [52]. Compared with these methods, statistical distance (SD) method is simple during the process of weight determination. Motivated by relevant researches about SD method, this paper utilizes extended SD method to identify the objective weights of CRs under linguistic Z-number environment.
According to the distance of LZNs introduced in Definition 7, the weights of CRs can be determined by the following formulae:
v j = 1 m i = 1 m v i j ,
σ j 2 = 1 m i = 1 m d ( v i j , v j ) 2 ,
w j = σ j 2 j = 1 n σ j 2 .

4.3. Prioritizing All the ECs

In order to rank engineering characteristics (ECs), this subsection extends TOPSIS method into LZNs environment. Next, the algorithm of the extended TOPSIS method is listed below.
Step T1. Identify the positive ideal solution (PIS) and negative ideal solution (NIS).
In terms of the comparison method in Definition 6, PIS and NIS can be defined as:
v j + = m a x i v i j , i = 1 , 2 , , m   for   j = 1 , 2 , , n ,
v j = m i n i v i j , i = 1 , 2 , , m   for   j = 1 , 2 , , n .
Step T2. Calculate the distances between CRs and both PIS and NIS.
Using the distance of LZNs in Definition 7, the distances between CRs and both PIS and NIS can be computed by:
d i + = j = 1 n w j d ( v i j , v j + ) ,   for   i = 1 , 2 , , m ,
d i = j = 1 n w j d ( v i j , v j ) ,   for   i = 1 , 2 , , m ,
where w j is the combined weights of C R j .
Step T3. Compute the closeness degrees.
The closeness degree can be calculated by Equation (27)
C i = d i d i + d i + ,   for   i = 1 , 2 , , m ,
It easily can be seen that the bigger is C i , the prior is E C i at the ranking order.

4.4. The Decision-Making Procedures

The brief steps of the proposed approach are summarized below:
Step 1. Establish the decision matrices.
The evaluation values u i j k should be provided for the relationships between engineering characteristics E C i ( i = 1 , 2 , , m ) and customer requirements C R j ( j = 1 , 2 , , n ) by the QFD team members T M k ( k = 1 , 2 , , h ) , in which u i j k is in the form of LZNs.
Step 2. Determine the combined weights of TMs.
This paper employs Equations (8) and (9) to compute the subjective weights of TMs. Then, Equations (10)–(16) are utilized to calculate the objective weights of TMs. Finally, Equation (6) is used to determine the combined weights of TMs λ = ( λ 1 , λ 2 , , λ h ) .
Step 3. Determine the combined weights of CRs.
Equations (17)–(19) are used to calculate the subjective weights of CRs. Next, Equations (20)–(22) are applied to compute the objective weights of CRs. Finally, Equation (7) is utilized to identify the integrated weights of CRs w = ( w 1 , w 2 , , w n ) .
Step 4. Calculate the closeness degree.
This paper utilizes Equations (23)–(27) to obtain the closeness degree to positive ideal solution C i ( i = 1 , 2 , , m ) .
Step 5. Obtain the ranking orders of ECs.
The ECs are ranked according to the closeness degree C i .

5. Illustration Example

In this section, a logistics service provider (LSP) [53] is absorbed as a case example to illustrate the effectiveness and efficiency of the proposed QFD model.

5.1. Problem Description

As the market changes rapidly and the customer demand shifts unpredictably, the LSP faces more and more serious challenges. For instance, many clients request better service with shorter delivery time, greater order accuracy and lower service cost. That is, in order to maintain core competitiveness, LSP should transform the traditional operations to just-in-time service. In this sense, it is essential to have research on the operations by taking customer requirements into accounts. Owing to the features of QFD model, this paper uses the proposed MCDM approach to determine the importance ratings of operation characteristics according to customer service demands.
In the case, three experts are selected as the QFD team members, including a senior manager of quality control T M 1 , a representative of the client T M 2 , and a professor of industry engineering T M 3 . After a mass of investigations and discussions, six customer requirements ( C R 1 , C R 2 , C R 3 , C R 4 , C R 5 , C R 6 ) and eight engineering characteristics ( E C 1 , E C 2 , E C 3 , E C 4 , E C 5 , E C 6 , E C 7 , E C 8 ) are defined for LSP as shown in Table 2.

5.2. Implementation of the Proposed Model

In this part, the proposed QFD model considering combined subjective and objective weights is applied for the case of LSP.
Step 1. Establish the decision matrices.
Since LZNs can express fuzzy and uncertain information perfectly, it is employed to evaluate the relationships between CRs and ECs by QFD expert team members. For a LZN Z = A , B , the fuzzy restriction A is a linguistic term in set S 1 :
S 1 = { s 0 = Very   low ,   s 1 = low ,   s 2 = Medium   low , s 3 = Medium , s 4 = Medium   high ,   s 5 = High , s 6 = Very   High } ;
Meanwhile, the reliability measure B is a linguistic term in set S 2 :
S 2 = { s 0 = Uncertain ,   s 1 = Slightly   Uncertain , s 2 = Medium , s 3 = Slightly   Sure ,   s 4 = Sure }
The relationship evaluation information between CRs and ECs is represented in Table 3 given by three QFD experts.
Step 2. Determine the combined weights of TMs.
(1)
The CEO of LSP thinks that the senior manager of quality control T M 1 accounts for the maximum weights and the professor of industry engineering T M 3 is designate as the minimum weights. Further, the relative importance degrees provided by CEO are θ 2 = 1.2 and θ 3 = 1.8 according to Table 1. By the means of Equations (8) and (9), the subjective weight vector of TMs are λ = ( 0.4355 , 0.3629 , 0.2016 ) T .
(2)
On the basis of the operational method and comparison method of LZNs in Definition 5 and
(3)
Definition 6, the ideal assessments can be obtained as shown in Table 4, Table 5 and Table 6. The differences between personal evaluation of TMs and ideal assessments can be computed by Equations (13)–(15). The calculated results are listed below:
d i f f ( U 1 , U + ) = 0.1758 ;   d i f f ( U 2 , U + ) = 0.1166 ;   d i f f ( U 3 , U + ) = 0.1626 ;
d i f f ( U 1 , U m a x ) = 0.1997 ;   d i f f ( U 2 , U m a x ) = 0.1814 ;   d i f f ( U 3 , U m a x ) = 0.1749 ;
d i f f ( U 1 , U m i n ) = 0.1623 ;   d i f f ( U 2 , U m i n ) = 0.1823 ;   d i f f ( U 3 , U m i n ) = 0.1984 .
By Equation (16), the objective weights vector of TMs can be determined a λ = ( 0.3220 , 0.3476 , 0.3304 ) T .
(4)
From the above analyses, we assume that the preference coefficient τ is 0.5, the integrated weights of TMs can be identified by Equation (6) which is λ = ( 0.3788 , 0.3552 , 0.2660 ) T .
Step 3. Determine the combined weights of CRs.
(1)
With QFD experts’ consultations, CRs are sorted according to their importance and the descending order is Lead-time C R 1 , Accuracy C R 2 , Reliability C R 5 , Fill rate C R 3 , Flexibility C R 4 , and Frequency C R 6 . Next, the comparative importance values between adjacent CRs are provided which are declared in Table 7. Finally, the subjective weights of CRs can be derived by Equations (17) and (19), the results are shown in Table 8. That is, CRs’ subjective weight vector is w = ( 0.32 , 0.24 , 0.12 , 0.10 , 0.15 , 0.07 ) T .
(2)
The comprehensive decision matrix, as shown in Table 9, can be obtained by integrating the weights of QFD experts and their evaluation values. According to the Equations (20)–(22), the objective weights vector of CRs is w = ( 0.1641 , 0.2669 , 0.1353 , 0.1546 , 0.1519 , 0.1272 ) T .
(3)
We suppose that the preference coefficient ξ is 0.5, then the combined weights vector of CRs is calculated by Equation (7) which is w = ( 0.2421 , 0.2534 , 0.1277 , 0.1273 , 0.1509 , 0.0986 ) T .
Step 4. Calculate the closeness degree.
Drawing on the comparison method in Definition 6, the PIS and NIS of engineering characteristics (EC) are acquired as follows:
v + = { s 5.11 , s 2.75 , s 5.64 , s 3.37 , s 3.89 , s 3.05 , s 4.51 , s 2.15 , s 4.98 , s 2.61 , s 4.89 , s 2.64 } ;
v = { s 1.89 , s 1.02 , s 0.36 , s 4.00 , s 0.36 , s 0.00 , s 1.77 , s 0.40 , s 0.98 , s 0.27 , s 0.36 , s 3.00 } .
The closeness degrees C i can be computed for i = 1 , 2 , , m by Equations (25)–(27) which are illustrated in Table 10.
Step 5. Obtain the ranking orders of ECs.
In terms of the closeness degrees of ECs, the ranking order of ECs is E C 1 E C 2 E C 8 E C 4 E C 3 E C 5 E C 7 E C 6 . Note that the most important three ECs are Just-in time E C 1 , Accuracy E C 2 , and Frequency E C 8 .

5.3. Sensitivity Analysis

This subsection studies a sensitivity analysis of different LSFs is.
In the proposed QFD model, three LSFs are absorbed to deal with linguistic terms of the utilized LZNs. If different LSFs are employed, different ranking results of ECs may be obtained. Thus, we should alter the LSFs to see how the ranking orders are varied. What’s more, for a LZN Z = A , B , the fuzzy restriction A is a linguistic term which expresses the correlation degree between CRs and ECs in this paper while the reliability measure B is the uncertainty estimation of A provided by TMs. Namely, it is worth noting that two parts of LZNs describe disparate evaluation preferences. Therefore, the two components of LZNs may choose different LSFs.
On the one hand, we use the first LSF f 1 to handle the fuzzy restriction A and change the LSFs to dispose reliability measure B . On the other hand, different LSFs are utilized to transform linguistic information A to a numeric value and the LSF of linguistic term B is always f 1 . The results are offered in Table 11.
As the results shown in Table 11, no matter how the utilized LSFs of the fuzzy restriction A and reliability measure B change, the Just-in time E C 1 , Accuracy E C 2 , and Frequency E C 8 are always the three engineering characteristics needed to pay close attention to mostly. As a consequence, the robustness of proposed QFD model can be concluded.

5.4. Comparative Analysis

In this subsection, two cases of comparative analysis are conducted to demonstrate the validity and practicality of the developed model in this study. One is to compare the proposed method with state-of -the-art MCDM method using LZNs. The other is to make a comparison between the proposed QFD model and other QFD models.
Case 1: The proposed method is compared with existing methods under LZNs environment.
At present, only one MCDM method is investigated with LZNs [16]. In this method, Wang et al. defined the distance measure of LZNs and extended the TODIM method based on Choquet Integral and distance measure. Hence, Wang et al.’s method is applied in the LSP example to find the most significant engineering characteristics. However, since the criteria weights were confirmed by a subjective method by Wang et al., the weights are assigned as the combined weights of CRs in this study. And the ranking orders of ECs are presented in Table 12.
From Table 12, the most vital engineering characteristic is E C 1 and the least crucial engineering characteristic is E C 6 for the both methods with LZNs. As a result, the feasibility of proposed method is proved. Meanwhile, there are two differences between these two methods. The first is that the computations of two methods are diverse. Wang et al.’s method needed to calculate the dominance degree of one engineering characteristic over another engineering characteristic with respect to a certain customer requirement. In other words, it should make a comparison between each pair of engineering characteristic regarding every customer requirement. However, if a lot of ECs are adopted in the case of LSP, it will cause the calculated amounts of Wang et al.’s method increase tremendously. By contrast, the proposed method only requires comparing engineering characteristics with PIS and NIS which haves less calculated quantities. The second is that the weights determining methods of CRs are distinctive in these two methods using LZNs. Wang et al. employed a subjective method to identify the weights. Nevertheless, in the practical processes, it is not enough to derive weights only by a subjective method. The objective weights determining methods are also wanted just like the proposed method does in this paper. In conclusion, the proposed method is more likely to be applied in real-life cases than Wang et al.’s method. Thus, the utility of proposed method is validated.
Case 2: The proposed QFD model is compared with extant QFD models.
For all we know, many MCDM approaches have been integrated into conventional QFD models to address the problems of them [3,6,54]. In order to demonstrate the strength of proposed QFD model, these QFD models are selected to do a comparative analysis. However, Yazdani et al. [6], Liao and Kao [54], Wu et al. [3], and this paper use crisp numbers (CNs), triangular fuzzy numbers (TFNs), hesitant fuzzy numbers (HFNs), and LZNs for evaluation information representations, separately. Namely, it is hard to compare these QFD models in the examples of this study due to different information expressions. Thus, we only make a comparison of the frameworks among these integrated QFD models. Furthermore, the frameworks of QFD models [3,6,54] are declared in Table 13.
As we can see in Table 13, Yazdani et al. [6] utilized crisp numbers to describe the relationship between CRs and ECs. Meanwhile, Liao and Kao [54] and Wu et al. [3] employed fuzzy numbers for information description. By contrast, this paper uses LZNs to denote the correlation between CRs and ECs. As for the weights of TMs, Wu et al. [3] only took objective weighting method into accounts while this paper considers both subjective and objective weighting method. What’s more, as for the weights of CRs, Yazdani et al. [6], Liao and Kao [54], and Wu et al. [3] determined the weights of CRs by merely subjective weighting method while this paper derives a combined subjective and objective weights of CRs as well.
Through the aforementioned comparisons, the advantages of proposed QFD models can be summarized as follows:
(1)
Compared with crisp numbers and fuzzy numbers, LZNs can take fuzzy information and the reliability of information into considerations at the same time. Owing to powerful ability of LZNs’ expressions, the proposed framework using LZNs can depict evaluation information perfectly.
(2)
It is inadvisable that only adopting objective method to obtain weights of TMs. Because some data is be hard to obtain and quantify. The QFD experts have discrepant professional backgrounds and areas, the subjective preferences to the importance degree of expert TMs should also be thought of. With regard to the weights of CRs, it may cause large deviation when only subjective weighting method is employed. Therefore, the QFD model with integrated subjective and objective weights of CRs may be more reasonable.
(3)
Compared with other MCDM method using LZNs, the proposed method only requires comparing engineering characteristics with PIS and NIS and avoids the problem: with the dimension increased, the calculated quantity rapidly does.

6. Conclusions

Since LZNs can deal with the uncertain and unreliable information well which exists in actual situations commonly, it is imperative to investigate the QFD framework utilizing LZNs. In this study, at first place, the weights of team members are identified by G1 method and improved MC method, respectively. Then, SWARA and SD methods are put forward to derive the subjective and objective weights of customer requirements. At this point, the combined weights of experts and customer requirements are determined. Next, the extended TOPSIS is proposed to aggregate the correlation between customer requirements and engineering characteristics. Finally, a case of logistics service provider is studied by the proposed QFD model to verify the efficiency and applicability of it. Moreover, a sensitivity analysis of LSFs is conducted to prove the robustness of our model. A comparison with the methods under LZN environment, and comparative analysis are conducted, the results indict the proposed QFD model has reliability. The way to determine the weights of CRs is a combined subjective and objective method, and the proposed method is more likely to be applied in real-life cases.
There are three main contributions in this paper. First, this paper develops an improved MC method which takes both positive and negative ideal assessment into account. Second, the novel QFD model is proposed by using LZNs and considering the combined weights of experts and customer requirements. Because LZNs can take fuzzy information and the reliability of information into considerations simultaneously, the proposed method can depict evaluation information more flexible and precisely. Third, the proposed QFD framework is utilized to discover the most significant engineering characteristics, which helps the logistics service providers to enhance their services. Compared with other MCDM method using LZNs, it only requires comparing engineering characteristics with PIS and NIS, which have fewer calculated quantities.
In the future, researches can be made on these aspects. First of all, the interactive relationships among customer requirements could be paid attention to in the stage of weights determination and information fusion. Subsequently, both qualitative and quantitative data can be considered in the meantime to strengthen the completeness of information expression. Ultimately, the proposed QFD can be applied to improve operational services or product design, such as medical services [55], industrial wastewater management [56], and energy selection [57].

Author Contributions

C.S. and J.-Q.W. developed methodology and provided the original draft preparation. J.-Q.W. and J.-B.L. analyzed the data. All authors have read and agreed to the published version of the manuscript.

Funding

This work is supported by the National Natural Science Foundation of China (Grant No. 71871228).

Conflicts of Interest

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

References

  1. Akao, Y. New product development and quality assurance–quality deployment system. Stand. Qual. Control 1972, 25, 7–14. [Google Scholar]
  2. Osiro, L.; Lima-Junior, F.R.; Carpinetti, L.C.R. A group decision model based on quality function deployment and hesitant fuzzy for selecting supply chain sustainability metrics. J. Clean. Prod. 2018, 183, 964–978. [Google Scholar] [CrossRef]
  3. Wu, S.-M.; Liu, H.-C.; Wang, L.-E. Hesitant fuzzy integrated MCDM approach for quality function deployment: A case study in electric vehicle. Int. J. Prod. Res. 2017, 55, 4436–4449. [Google Scholar] [CrossRef]
  4. Dincer, H.; Yüksel, S.; Martinez, L. Balanced scorecard-based Analysis about European Energy Investment Policies: A hybrid hesitant fuzzy decision-making approach with Quality Function Deployment. Expert Syst. Appl. 2019, 115, 152–171. [Google Scholar] [CrossRef]
  5. Onar, S.Ç.; Büyüközkan, G.; Öztayşi, B.; Kahraman, C. A new hesitant fuzzy QFD approach: An application to computer workstation selection. Appl. Soft Comput. 2016, 46, 1–16. [Google Scholar] [CrossRef]
  6. Yazdani, M.; Chatterjee, P.; Zavadskas, E.K.; Zolfani, S.H. Integrated QFD-MCDM framework for green supplier selection. J. Clean. Prod. 2017, 142, 3728–3740. [Google Scholar] [CrossRef]
  7. Dolgun, L.E.; Köksal, G. Effective use of quality function deployment and Kansei engineering for product planning with sensory customer requirements: A plain yogurt case. Qual. Eng. 2018, 30, 569–582. [Google Scholar] [CrossRef]
  8. Yazdani, M.; Kahraman, C.; Zarate, P.; Onar, S.C. A fuzzy multi attribute decision framework with integration of QFD and grey relational analysis. Expert Syst. Appl. 2019, 115, 474–485. [Google Scholar] [CrossRef] [Green Version]
  9. Chowdhury, M.M.H.; Quaddus, M.A. A multiple objective optimization based QFD approach for efficient resilient strategies to mitigate supply chain vulnerabilities: The case of garment industry of Bangladesh. Omega 2015, 57, 5–21. [Google Scholar] [CrossRef]
  10. Jia, W.; Liu, Z.; Lin, Z.; Qiu, C.; Tan, J. Quantification for the importance degree of engineering characteristics with a multi-level hierarchical structure in QFD. Int. J. Prod. Res. 2016, 54, 1627–1649. [Google Scholar] [CrossRef]
  11. Song, W.; Ming, X.; Han, Y. Prioritising technical attributes in QFD under vague environment: A rough-grey relational analysis approach. Int. J. Prod. Res. 2014, 52, 5528–5545. [Google Scholar] [CrossRef]
  12. Zadeh, L.A. A note on Z-numbers. Inf. Sci. 2011, 181, 2923–2932. [Google Scholar] [CrossRef]
  13. Qiao, D.; Shen, K.-W.; Wang, J.-Q.; Wang, T.-L. Multi-criteria PROMETHEE method based on possibility degree with Z-numbers under uncertain linguistic environment. J. Ambient Intell. Humaniz. Comput. 2019, 1–15. [Google Scholar] [CrossRef]
  14. Shen, K.-W.; Wang, X.-K.; Qiao, D.; Wang, J.-Q. Extended Z-MABAC method based on regret theory and directed distance for regional circular economy development program selection with Z-information. IEEE Trans. Fuzzy Syst. 2019. [Google Scholar] [CrossRef]
  15. Qiao, D.; Wang, X.-K.; Wang, J.-Q.; Chen, K. Cross Entropy for Discrete Z-numbers and Its Application in Multi-Criteria Decision-Making. Int. J. Fuzzy Syst. 2019, 21, 1786–1800. [Google Scholar] [CrossRef]
  16. Wang, J.-Q.; Cao, Y.-X.; Zhang, H.-Y. Multi-criteria decision-making method based on distance measure and Choquet integral for linguistic Z-numbers. Cogn. Comput. 2017, 9, 827–842. [Google Scholar] [CrossRef]
  17. Huang, J.; Xu, D.-H.; Liu, H.-C.; Song, M.-S. A New Model for Failure Mode and Effect Analysis Integrating Linguistic Z-numbers and Projection Method. IEEE Trans. Fuzzy Syst. 2019. [Google Scholar] [CrossRef]
  18. Xu, F.; Liu, J.; Lin, S.; Yuan, J. A VIKOR-based approach for assessing the service performance of electric vehicle sharing programs: A case study in Beijing. J. Clean. Prod. 2017, 148, 254–267. [Google Scholar] [CrossRef]
  19. Wu, X.; Liao, H. A consensus-based probabilistic linguistic gained and lost dominance score method. Eur. J. Oper. Res. 2019, 272, 1017–1027. [Google Scholar] [CrossRef]
  20. Zavadskas, E.K.; Bausys, R.; Juodagalviene, B.; Garnyte-Sapranaviciene, I. Model for residential house element and material selection by neutrosophic MULTIMOORA method. Eng. Appl. Artif. Intell. 2017, 64, 315–324. [Google Scholar] [CrossRef] [Green Version]
  21. Liu, H.-C.; You, J.-X.; Chen, S.; Chen, Y.-Z. An integrated failure mode and effect analysis approach for accurate risk assessment under uncertainty. IIE Trans. 2016, 48, 1027–1042. [Google Scholar] [CrossRef]
  22. Saaty, T.L.; Kearns, K.P. The Analytic Hierarchy Process; RWS: Cologne, Germany, 1985. [Google Scholar]
  23. Hwang, C.-L.; Yoon, K. Multiple Attribute Decision Making: Methods and Applications a State-of-the-Art Survey; Springer Science & Business Media: Berlin, Germany, 1981; Volume 186. [Google Scholar]
  24. Roy, B. Classement et choix en présence de points de vue multiples. Rev. Fr. D’informatique Rech. Oper. 1968, 2, 57–75. [Google Scholar] [CrossRef]
  25. Opricovic, S.; Tzeng, G.-H. Compromise solution by MCDM methods: A comparative analysis of VIKOR and TOPSIS. Eur. J. Oper. Res. 2004, 156, 445–455. [Google Scholar] [CrossRef]
  26. Mareschal, B.; Brans, J.P.; Vincke, P. PROMETHEE: A New Family of Outranking Methods in Multicriteria Analysis; ULB—Universite Libre de Bruxelles: Bruxelles, Belgium, 1984. [Google Scholar]
  27. Julong, D. Introduction to grey system theory. J. Grey Syst. 1989, 1, 1–24. [Google Scholar]
  28. Rezaei, J. Best-worst multi-criteria decision-making method. Omega 2015, 53, 49–57. [Google Scholar] [CrossRef]
  29. Hsu, C.-H.; Chang, A.-Y.; Luo, W. Identifying key performance factors for sustainability development of SMEs–integrating QFD and fuzzy MADM methods. J. Clean. Prod. 2017, 161, 629–645. [Google Scholar] [CrossRef]
  30. Zaim, S.; Sevkli, M.; Camgöz-Akdağ, H.; Demirel, O.F.; Yayla, A.Y.; Delen, D. Use of ANP weighted crisp and fuzzy QFD for product development. Expert Syst. Appl. 2014, 41, 4464–4474. [Google Scholar] [CrossRef]
  31. Abdel-Basset, M.; Manogaran, G.; Mohamed, M.; Chilamkurti, N. Three-way decisions based on neutrosophic sets and AHP-QFD framework for supplier selection problem. Future Gener. Comput. Syst. 2018, 89, 19–30. [Google Scholar] [CrossRef]
  32. Wu, X.; Nie, L.; Xu, M. Robust fuzzy quality function deployment based on the mean-end-chain concept: Service station evaluation problem for rail catering services. Eur. J. Oper. Res. 2017, 263, 974–995. [Google Scholar] [CrossRef]
  33. Abdel-Basset, M.; Mohamed, R.; Zaied, A.E.-N.H.; Smarandache, F. A hybrid plithogenic decision-making approach with quality function deployment for selecting supply chain sustainability metrics. Symmetry 2019, 11, 903. [Google Scholar] [CrossRef] [Green Version]
  34. Yang, Q.; Chin, K.-S.; Li, Y.-L. A quality function deployment-based framework for the risk management of hazardous material transportation process. J. Loss Prev. Process Ind. 2018, 52, 81–92. [Google Scholar] [CrossRef]
  35. Tian, C.; Peng, J.-J.; Zhang, W.-Y.; Zhang, S.; Wang, J.-Q. Tourism environmental impact assessment based on improved AHP and picture fuzzy PROMETHEE II methods. Technol. Econ. Dev. Econ. 2019, 1–24. [Google Scholar] [CrossRef]
  36. Yu, S.M.; Wang, J.; Wang, J.Q. An extended TODIM approach with intuitionistic linguistic numbers. Int. Trans. Oper. Res. 2018, 25, 781–805. [Google Scholar] [CrossRef]
  37. Wang, J.Q.; Han, Z.Q.; Zhang, H.Y. Multi-criteria Group Decision-Making Method Based on Intuitionistic Interval Fuzzy Information. Group Decis. Negot. 2014, 23, 715–733. [Google Scholar] [CrossRef]
  38. Wang, L.; Wang, X.-K.; Peng, J.-J.; Wang, J.-Q. The differences in hotel selection among various types of travellers: A comparative analysis with a useful bounded rationality behavioural decision support model. Tour. Manag. 2020, 76, 103961. [Google Scholar] [CrossRef]
  39. Herrera, F.; Herrera-Viedma, E.; Verdegay, J.L. A model of consensus in group decision making under linguistic assessments. Fuzzy Sets Syst. 1996, 78, 73–87. [Google Scholar] [CrossRef]
  40. Zhang, H.-Y.; Peng, H.-G.; Wang, J.; Wang, J.-Q. An extended outranking approach for multi-criteria decision-making problems with linguistic intuitionistic fuzzy numbers. Appl. Soft Comput. 2017, 59, 462–474. [Google Scholar] [CrossRef]
  41. Bao, G.-Y.; Lian, X.-L.; He, M.; Wang, L.-L. Improved two-tuple linguistic representation model based on new linguistic evaluation scale. Control Decis. 2010, 25, 780–784. [Google Scholar]
  42. Kahneman, D.; Tversky, A. Prospect theory: An analysis of decision under risk. Econom. J. Econom. Soc. 1979, 47, 263–291. [Google Scholar] [CrossRef] [Green Version]
  43. Versaci, M.; La Foresta, F.; Morabito, F.C.; Angiulli, G. A fuzzy divergence approach for solving electrostatic identification problems for NDT applications. Int. J. Appl. Electromagn. Mech. 2018, 57, 133–146. [Google Scholar] [CrossRef]
  44. Versaci, M. Fuzzy approach and Eddy currents NDT/NDE devices in industrial applications. Electron. Lett. 2016, 52, 943–945. [Google Scholar] [CrossRef]
  45. Wang, L.; Zhang, H.; Wang, J.; Wu, G. Picture fuzzy multi-criteria group decision making method to hotel building energy efficiency retrofit project selection. Rairo Oper. Res. 2018. [Google Scholar] [CrossRef] [Green Version]
  46. Ruan, C.; Yang, J. Software quality evaluation model based on weighted mutation rate correction incompletion G1 combination weights. Math. Probl. Eng. 2014, 2014. [Google Scholar] [CrossRef]
  47. Zhang, X.; Xu, Z. Soft computing based on maximizing consensus and fuzzy TOPSIS approach to interval-valued intuitionistic fuzzy group decision making. Appl. Soft Comput. 2015, 26, 42–56. [Google Scholar] [CrossRef]
  48. Zolfani, S.H.; Chatterjee, P. Comparative evaluation of sustainable design based on Step-wise Weight Assessment Ratio Analysis (SWARA) and best worst method (BWM) methods: A perspective on household furnishing materials. Symmetry 2019, 11, 74. [Google Scholar] [CrossRef] [Green Version]
  49. Shahsavar, S.; Rad, A.J.; Afzal, P.; Nezafati, N.; Aghdam, M.A. Prospecting for polymetallic mineralization using step-wise weight assessment ratio analysis (SWARA) and fractal modeling in Aghkand Area, NW Iran. Arab. J. Geosci. 2019, 12, 248. [Google Scholar] [CrossRef]
  50. Nakhaei, J.; Bitarafan, M.; Lale Arefi, S.; Kapliński, O. Model for rapid assessment of vulnerability of office buildings to blast using SWARA and SMART methods (a case study of swiss re tower). J. Civ. Eng. Manag. 2016, 22, 831–843. [Google Scholar] [CrossRef] [Green Version]
  51. Selvachandran, G.; Quek, S.; Smarandache, F.; Broumi, S. An extended technique for order preference by similarity to an ideal solution (TOPSIS) with maximizing deviation method based on integrated weight measure for single-valued neutrosophic sets. Symmetry 2018, 10, 236. [Google Scholar] [CrossRef] [Green Version]
  52. Hu, J.; Yang, Y.; Zhang, X.; Chen, X. Similarity and entropy measures for hesitant fuzzy sets. Int. Trans. Oper. Res. 2018, 25, 857–886. [Google Scholar] [CrossRef] [Green Version]
  53. Oláh, J.; Karmazin, G.; Pető, K.; Popp, J. Information technology developments of logistics service providers in Hungary. Int. J. Logist. Res. Appl. 2018, 21, 332–344. [Google Scholar] [CrossRef]
  54. Liao, C.-N.; Kao, H.-P. An evaluation approach to logistics service using fuzzy theory, quality function development and goal programming. Comput. Ind. Eng. 2014, 68, 54–64. [Google Scholar] [CrossRef]
  55. Hu, J.; Zhang, X.; Yang, Y.; Liu, Y.; Chen, X. New doctors ranking system based on VIKOR method. Int. Trans. Oper. Res. 2020, 27, 1236–1261. [Google Scholar] [CrossRef]
  56. Shen, K.-W.; Li, L.; Wang, J.-Q. Circular economy model for recycling waste resources under government participation: A case study in industrial waste water circulation in China. Technol. Econ. Dev. Econ. 2020, 26, 21–47. [Google Scholar] [CrossRef]
  57. Peng, H.; Wang, J.; Zhang, H. Multi-criteria outranking method based on probability distribution with probabilistic linguistic information. Comput. Ind. Eng. 2020, 106318. [Google Scholar] [CrossRef]
Figure 1. The house of quality (HOQ).
Figure 1. The house of quality (HOQ).
Mathematics 08 00224 g001
Figure 2. Graphical representation of three LSFs.
Figure 2. Graphical representation of three LSFs.
Mathematics 08 00224 g002
Figure 3. The conceptual framework of the proposed model for QFD.
Figure 3. The conceptual framework of the proposed model for QFD.
Mathematics 08 00224 g003
Table 1. The description of the relative importance degree.
Table 1. The description of the relative importance degree.
θ k Description
1.0 T M k 1 is equally important as T M k
1.2 T M k 1 is slightly more important than T M k
1.4 T M k 1 is more important than T M k
1.6 T M k 1 is strongly more important than T M k
1.8 T M k 1 is extremely more important than T M k
Table 2. CRs and ECs in the case study.
Table 2. CRs and ECs in the case study.
CRsCustomer RequirementsECsEngineering Characteristics
C R 1 Lead-time E C 1 Just-in time
C R 2 Accuracy E C 2 Information technology
C R 3 Fill rate E C 3 Order picking optimization
C R 4 Flexibility E C 4 Demand forecasting methods
C R 5 Reliability E C 5 Quality of services
C R 6 Frequency E C 6 Customer relationship management
E C 7 Warehouses lay-out optimization
E C 8 Information sharing and mutual trust
Table 3. The relationship evaluation information between CRs and ECs.
Table 3. The relationship evaluation information between CRs and ECs.
E C 1 E C 2 E C 3 E C 4 E C 5 E C 6 E C 7 E C 8
T M 1 C R 1 s 5 , s 3 s 6 , s 4 s 2 , s 4 s 0 , s 2 s 1 , s 0 s 1 , s 3 s 2 , s 3 s 1 , s 2
C R 2 s 4 , s 3 s 3 , s 4 s 0 , s 1 s 0 , s 2 s 2 , s 2 s 1 , s 1 s 0 , s 2 s 6 , s 4
C R 3 s 1 , s 3 s 4 , s 3 s 0 , s 1 s 3 , s 2 s 1 , s 0 s 2 , s 2 s 0 , s 4 s 1 , s 3
C R 4 s 5 , s 3 s 4 , s 1 s 0 , s 1 s 1 , s 2 s 0 , s 2 s 3 , s 4 s 0 , s 3 s 1 , s 4
C R 5 s 4 , s 4 s 4 , s 1 s 1 , s 2 s 2 , s 3 s 1 , s 2 s 0 , s 1 s 1 , s 2 s 3 , s 1
C R 6 s 4 , s 1 s 2 , s 2 s 4 , s 3 s 1 , s 1 s 5 , s 2 s 1 , s 0 s 4 , s 2 s 0 , s 4
T M 2 C R 1 s 3 , s 4 s 5 , s 1 s 3 , s 2 s 2 , s 3 s 2 , s 2 s 2 , s 0 s 1 , s 2 s 1 , s 3
C R 2 s 6 , s 2 s 2 , s 3 s 2 , s 1 s 1 , s 3 s 1 , s 2 s 2 , s 2 s 1 , s 4 s 5 , s 2
C R 3 s 2 , s 2 s 2 , s 4 s 1 , s 0 s 4 , s 3 s 2 , s 1 s 3 , s 1 s 2 , s 4 s 3 , s 2
C R 4 s 4 , s 2 s 3 , s 2 s 1 , s 2 s 3 , s 4 s 2 , s 1 s 5 , s 2 s 1 , s 2 s 2 , s 4
C R 5 s 6 , s 2 s 5 , s 2 s 0 , s 3 s 1 , s 2 s 2 , s 2 s 2 , s 0 s 1 , s 3 s 4 , s 2
C R 6 s 3 , s 2 s 1 , s 2 s 5 , s 2 s 3 , s 2 s 4 , s 3 s 2 , s 1 s 3 , s 0 s 1 , s 3
T M 3 C R 1 s 3 , s 3 s 4 , s 3 s 2 , s 1 s 0 , s 3 s 1 , s 2 s 3 , s 1 s 1 , s 3 s 2 , s 2
C R 2 s 5 , s 3 s 1 , s 2 s 3 , s 2 s 2 , s 4 s 1 , s 1 s 0 , s 4 s 0 , s 2 s 6 , s 4
C R 3 s 4 , s 1 s 1 , s 3 s 0 , s 2 s 5 , s 4 s 3 , s 3 s 2 , s 1 s 1 , s 2 s 5 , s 1
C R 4 s 3 , s 1 s 4 , s 1 s 2 , s 3 s 2 , s 2 s 4 , s 0 s 6 , s 1 s 3 , s 2 s 4 , s 4
C R 5 s 5 , s 2 s 4 , s 1 s 1 , s 2 s 2 , s 3 s 1 , s 3 s 1 , s 1 s 2 , s 2 s 3 , s 4
C R 6 s 2 , s 4 s 3 , s 0 s 6 , s 3 s 4 , s 1 s 5 , s 1 s 1 , s 3 s 2 , s 2 s 0 , s 2
Table 4. The positive ideal assessment.
Table 4. The positive ideal assessment.
E C 1 E C 2 E C 3 E C 4 E C 5 E C 6 E C 7 E C 8
C R 1 s 3.67 , s 3.27 s 5.00 , s 2.73 s 2.33 , s 2.29 s 0.67 , s 3.00 s 1.33 , s 1.50 s 2.00 , s 1.00 s 1.33 , s 2.75 s 1.33 , s 2.25
C R 2 s 5.00 , s 2.60 s 2.00 , s 3.33 s 1.67 , s 1.60 s 1.00 , s 3.67 s 1.33 , s 1.75 s 1.00 , s 1.67 s 0.33 , s 4.00 s 5.67 , s 3.41
C R 3 s 2.33 , s 1.57 s 2.33 , s 3.29 s 0.33 , s 0.00 s 4.00 , s 3.17 s 2.00 , s 1.83 s 2.33 , s 1.29 s 1.00 , s 3.33 s 3.00 , s 1.56
C R 4 s 4.00 , s 2.17 s 3.67 , s 1.27 s 1.00 , s 2.67 s 2.00 , s 3.00 s 2.00 , s 0.33 s 4.67 , s 2.00 s 1.33 , s 2.00 s 2.33 , s 4.00
C R 5 s 5.00 , s 2.53 s 4.33 , s 1.38 s 0.67 , s 2.00 s 1.67 , s 2.80 s 1.33 , s 2.25 s 1.00 , s 0.33 s 1.33 , s 2.25 s 3.33 , s 2.30
C R 6 s 3.00 , s 2.00 s 2.00 , s 1.00 s 5.00 , s 2.67 s 2.67 , s 1.38 s 4.67 , s 1.93 s 1.33 , s 1.25 s 3.00 , s 1.33 s 0.33 , s 3.00
Table 5. The negative ideal assessment—maximum value.
Table 5. The negative ideal assessment—maximum value.
E C 1 E C 2 E C 3 E C 4 E C 5 E C 6 E C 7 E C 8
C R 1 s 5 , s 3 s 6 , s 4 s 2 , s 4 s 2 , s 3 s 2 , s 2 s 3 , s 1 s 2 , s 3 s 2 , s 2
C R 2 s 5 , s 3 s 3 , s 4 s 3 , s 2 s 2 , s 4 s 2 , s 2 s 2 , s 2 s 1 , s 4 s 6 , s 4
C R 3 s 4 , s 1 s 4 , s 3 s 1 , s 0 s 5 , s 4 s 3 , s 3 s 2 , s 2 s 2 , s 4 s 3 , s 2
C R 4 s 5 , s 3 s 3 , s 2 s 2 , s 3 s 3 , s 4 s 2 , s 1 s 3 , s 4 s 3 , s 2 s 4 , s 4
C R 5 s 4 , s 4 s 5 , s 2 s 1 , s 2 s 2 , s 3 s 2 , s 2 s 1 , s 1 s 2 , s 2 s 3 , s 4
C R 6 s 2 , s 4 s 2 , s 2 s 6 , s 3 s 3 , s 2 s 4 , s 3 s 1 , s 3 s 4 , s 2 s 1 , s 3
Table 6. The negative ideal assessment—minimum value.
Table 6. The negative ideal assessment—minimum value.
E C 1 E C 2 E C 3 E C 4 E C 5 E C 6 E C 7 E C 8
C R 1 s 3 , s 3 s 5 , s 1 s 2 , s 1 s 0 , s 3 s 1 , s 0 s 2 , s 0 s 1 , s 2 s 1 , s 2
C R 2 s 4 , s 3 s 1 , s 2 s 0 , s 1 s 0 , s 2 s 1 , s 1 s 0 , s 4 s 0 , s 2 s 5 , s 2
C R 3 s 1 , s 3 s 1 , s 3 s 0 , s 2 s 3 , s 2 s 1 , s 0 s 2 , s 1 s 0 , s 4 s 1 , s 3
C R 4 s 3 , s 1 s 4 , s 1 s 0 , s 1 s 1 , s 2 s 0 , s 2 s 6 , s 1 s 0 , s 3 s 1 , s 4
C R 5 s 5 , s 2 s 4 , s 1 s 0 , s 3 s 1 , s 2 s 1 , s 2 s 0 , s 1 s 1 , s 2 s 3 , s 1
C R 6 s 4 , s 1 s 3 , s 0 s 5 , s 2 s 1 , s 1 s 5 , s 1 s 1 , s 0 s 3 , s 0 s 0 , s 2
Table 7. The comparative importance values between adjacent CRs.
Table 7. The comparative importance values between adjacent CRs.
C R 1 2 C R 2 5 C R 5 3 C R 3 4 C R 4 6
T M 1 0.200.600.100.000.50
T M 2 0.400.500.300.400.70
T M 3 0.500.600.400.300.40
Table 8. The final subjective weights of CRs.
Table 8. The final subjective weights of CRs.
Weighted   Average   Comparative   Importance   Value s j The   Characteristic   of   the   Comparative   Importance   x j The   Recalculated   Weight   y j The   Final   Subjective   Weight   w j
C R 1 -1.001.000.32
C R 2 0.351.350.740.24
C R 5 0.561.560.470.15
C R 3 0.251.250.380.12
C R 4 0.221.220.310.10
C R 6 0.541.540.200.07
Table 9. The comprehensive decision matrix.
Table 9. The comprehensive decision matrix.
E C 1 E C 2 E C 3 E C 4 E C 5 E C 6 E C 7 E C 8
C R 1 s 3.76 , s 3.28 s 5.11 , s 2.75 s 2.36 , s 2.42 s 0.71 , s 3.00 s 1.36 , s 1.44 s 1.89 , s 1.02 s 1.38 , s 2.74 s 1.27 , s 2.28
C R 2 s 4.98 , s 2.57 s 2.11 , s 3.41 s 1.51 , s 1.53 s 0.89 , s 3.60 s 1.38 , s 1.81 s 1.09 , s 1.65 s 0.36 , s 4.00 s 5.64 , s 3.37
C R 3 s 2.15 , s 1.68 s 2.49 , s 3.29 s 0.35 , s 0.00 s 3.89 , s 3.05 s 1.89 , s 1.65 s 2.36 , s 1.32 s 0.98 , s 3.46 s 2.77 , s 1.66
C R 4 s 4.11 , s 2.27 s 3.64 , s 1.29 s 0.89 , s 2.60 s 1.98 , s 3.08 s 1.97 , s 0.40 s 4.51 , s 2.15 s 1.15 , s 2.00 s 2.15 , s 4.00
C R 5 s 4.98 , s 2.61 s 4.36 , s 1.41 s 0.64 , s 2.00 s 1.64 , s 2.78 s 1.36 , s 2.20 s 0.98 , s 0.27 s 1.27 , s 2.28 s 3.36 , s 2.14
C R 6 s 3.11 , s 1.86 s 1.91 , s 1.16 s 4.89 , s 2.64 s 2.51 , s 1.42 s 4.64 , s 2.02 s 1.36 , s 1.11 s 3.11 , s 1.32 s 0.36 , s 3.00
Table 10. The closeness degrees of ECs.
Table 10. The closeness degrees of ECs.
E C 1 E C 2 E C 3 E C 4 E C 5 E C 6 E C 7 E C 8
d i + 0.17230.27900.50290.47650.50460.50690.55330.3080
d i 0.49010.38450.28060.30540.21950.18290.23010.4213
C i 0.73990.57940.35810.39060.30310.26520.29370.5777
Table 11. Ranking results with different LSFs.
Table 11. Ranking results with different LSFs.
f * ( A ) g * ( B ) Ranking Results
f * ( A ) = f 1 ( A ) g * ( B ) = f 1 ( B ) E C 1 E C 2 E C 8 E C 4 E C 3 E C 5 E C 7 E C 6
f * ( A ) = f 1 ( A ) g * ( B ) = f 2 ( B ) E C 1 E C 2 E C 8 E C 4 E C 3 E C 5 E C 6 E C 7
f * ( A ) = f 1 ( A ) g * ( B ) = f 3 ( B ) E C 1 E C 2 E C 8 E C 4 E C 3 E C 5 E C 7 E C 6
f * ( A ) = f 1 ( A ) g * ( B ) = f 1 ( B ) E C 1 E C 2 E C 8 E C 4 E C 3 E C 5 E C 7 E C 6
f * ( A ) = f 2 ( A ) g * ( B ) = f 1 ( B ) E C 1 E C 8 E C 2 E C 4 E C 3 E C 7 E C 5 E C 6
f * ( A ) = f 3 ( A ) g * ( B ) = f 1 ( B ) E C 1 E C 2 E C 8 E C 4 E C 3 E C 5 E C 7 E C 6
Table 12. The comparison between the proposed method and existing methods utilizing LZNs.
Table 12. The comparison between the proposed method and existing methods utilizing LZNs.
MethodRank Orders
Wang et al.’s method ( θ = 1 ) [16] E C 1 E C 2 E C 8 E C 4 E C 3 E C 7 E C 5 E C 6
The proposed method ( f * ( A ) = f 1 ( A ) , g * ( B ) = f 1 ( B ) ) E C 1 E C 2 E C 8 E C 4 E C 3 E C 5 E C 7 E C 6
Table 13. The comparative analysis of integrated QFD models.
Table 13. The comparative analysis of integrated QFD models.
Yazdani et al. [6]Liao and Kao [54]Wu et al. [3]This Paper
Information representationCNsTFNsHFNsLZNs
Weights of TMsSubjective weighting method×××
Objective weighting method××
Weights of CRsSubjective weighting method
Objective weighting method×××

Share and Cite

MDPI and ACS Style

Song, C.; Wang, J.-Q.; Li, J.-B. New Framework for Quality Function Deployment Using Linguistic Z-Numbers. Mathematics 2020, 8, 224. https://doi.org/10.3390/math8020224

AMA Style

Song C, Wang J-Q, Li J-B. New Framework for Quality Function Deployment Using Linguistic Z-Numbers. Mathematics. 2020; 8(2):224. https://doi.org/10.3390/math8020224

Chicago/Turabian Style

Song, Chao, Jian-Qiang Wang, and Jun-Bo Li. 2020. "New Framework for Quality Function Deployment Using Linguistic Z-Numbers" Mathematics 8, no. 2: 224. https://doi.org/10.3390/math8020224

APA Style

Song, C., Wang, J. -Q., & Li, J. -B. (2020). New Framework for Quality Function Deployment Using Linguistic Z-Numbers. Mathematics, 8(2), 224. https://doi.org/10.3390/math8020224

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop