Threat evaluation is extremely important to decision makers in many situations, such as military application and physical protection systems. In this paper, a new threat assessment model based on interval number to deal with the intrinsic uncertainty and imprecision in combat environment is proposed. Both objective and subjective factors are taken into consideration in the proposed model. For the objective factors, the genetic algorithm (GA) is used to search out an optimal interval number representing all the attribute values of each object. In addition, for the subjective factors, the interval Analytic Hierarchy Process (AHP) is adopted to determine each object’s threat weight according to the experience of commanders/experts. Then a discounting method is proposed to integrate the objective and subjective factors. At last, the ideal of Technique for Order Preference by Similarity to Ideal Solution (TOPSIS) is applied to obtain the threat ranking of all the objects. A real application is used to illustrate the effectiveness of the proposed model.

1. Introduction

In the real world, it is always complex and difficult to make a decision and to deal with threat assessment in most situations, because the problems we are faced with are always intrinsic imprecision and uncertainty which is typically represented by the threat assessment for multiobject and multiattribute with uncertain information. Hence, it is meaningful to develop a reasonable model to solve this problem. The threat assessment of multiobject and multiattribute has been researched by many researchers. Generally speaking, threat assessment can be seen as a decision-making problem under uncertain environment. We first introduce some related works about decision making and threat assessment.

1.1. Decision Making

Decision making has been the popular topic of research in many fields since the concept was proposed. Saaty [1, 2] introduces the AHP which is a multicriteria decision-making approach in which factors are arranged in a hierarchic structure. Hall and Davis [3] introduce a value-based decision-making model which suggests that multiple perspectives may be achieved by considering a foundation of individual values. Zhang et al. [4] present a method for solving the stochastic MCDM problem and propose a new concept of stochastic dominance degree (SDD). Chiadamrong [5] introduces the concept of fuzzy set theory to overcoming the precision-based evaluation for manufacturing strategies selection. Hon et al. [6] make use of a procedure to solve the multiple attributes and multiple hierarchical system under fuzzy environment. This procedure is based on gradient eigenvector method or rating the fuzzy criteria weighting and the max-min paired elimination method for aggregation. Forgionne and Kohli [7] present a methodology that objectively rates decision technology system journals across multiple and relevant dimensions. Through this process, the AHP-based rating model and the DSS-delivery mechanism are the major contribution. Fan et al. [8] propose a method based on pairwise comparisons of alternatives with random evaluations to solve stochastic MCDM problem. By computing the probabilities of all possible results of pairwise comparisons and using the identification rule, superior, indifferent, and inferior probabilities on pairwise comparisons of alternatives are estimated. Al-Najjar and Alsyouf [9] make use of the advantage of both fuzzy logic and multiple criteria decision making to introduce an evaluation methodology. Wang [10] proposes a new multiple criteria decision-making method with incomplete certain information based on ternary AHP; in this process, the ternary comparison matrix of the alternatives is constructed to normalize priority vector of the alternatives. Similar researches are conducted by Chan et al. [11]. Due to the efficient modelling and handling of uncertain information, fuzzy set theory [1215] and evidence theory are also widely used in decision making [1623]. Some typical evidence theories with AHP are published in [24, 25].

1.2. Threat Assessment

A threat is defined as any act, entity, event, or phenomenon with the potential to harm a person or thing. In other words, a threat is a source of potential harm. Sometimes the word hazard or risk is used as a synonym for threat. The term “threat assessment” can be broadly interpreted as evaluation of impending danger or harm by a person, group, circumstance, or set of conditions. It is a pattern of activities involving detection and analysis of threat stimuli and the situations in which the threat is encountered. Threat or risk assessment serving as the foundation of regulatory decision making whether to take actions to reduce the loss has been researched by many researchers. Pedrycz et al. [26] introduce a scheme of risk assessment on a basis of classification results produced for some experimental data capturing the history of previous threat cases. Through this procedure, information granulation (fuzzy clustering) was used to reveal the structural relationships in the experimental data. van de Walle and Rutkowski [27] designed and developed FURIA (fuzzy relational incident analysis), a prototype fuzzy decision support system allowing individual group members to compare their individual assessment of a decision alternative to the assessment of the other group members. Looney and Liang [28] propose a multiphase approach to middle and high level data fusion with an application to situation and threat assessment, threat assessment with fuzzy belief network. Khadam and Kaluarachachi [29] propose a decision analysis framework which integrates probabilistic health risk assessment into comprehensive, yet simple, cost-based multicriteria decision analysis framework. Assmuth and Hildén [30] present an analysis of frameworks for information on integrated risk assessment and management, which focuses on issues related to environmental and health risks of chemicals and their regulation at the EU level. Chen and Ma [31] propose a method combining the cost of reducing uncertainty with the selection of risk assessment models for remediation decision of site contamination; in this study, they explore how to choose among equally plausible models for risk-related decisions, which is closely related to model uncertainty. Jeong et al. [32] introduce a risk assessment method which is developed by using risk matrix and fuzzy inference logic, on the basis of the radiological and nonradiological hazard for a decommissioning safety of a nuclear facility. Bonano et al. [33] developed and demonstrated a single framework integrating risk assessment and decision analysis methods for evaluating, ranking, and selecting preferred remediation alternatives at a contaminated site; in the decision analysis framework, the Analytic Hierarchy Process was selected as a technique to obtain stakeholder working group input. van Duijne et al. [34] outline the preconditions for risk assessment method and aim to clarify the difficulties and pitfalls in the subsequent stages of risk assessment. Kalantarnia et al. [35] demonstrate the use of Bayesian theory in quantitative risk assessment (QRA) and its application as a useful tool in dynamic risk assessment to prevent accidents, which overcomes the major disadvantages of QRA: its inability to update risk during the life of a process. In the past few years, the field of group decision and threat evaluation has been researched hotly, like in Deng et al. [36].

Evaluation for aerial objects is an important part in the aerial defensive operation, which can be described as the assessment of multiobject and multiattribute under uncertain environment. Due to the limitation of detection means or secrecy, the effective information we can obtain from the threat objects is little. Generally, the interval number can be used to represent the uncertain or imprecise information. In this paper, a new model is proposed to rank the threat objects based on interval numbers, which takes the subjective and objective factors into consideration. Two examples are used to illustrate the effectiveness of the proposed threat evaluation model.

This paper is organized as follows. Section 2 mainly introduces several concepts of interval numbers, the classical genetic algorithm, interval AHP, and the TOPSIS. Section 3 introduces the model of threat assessment in detail. Section 4 introduces a simple numerical example and a real application to illustrate the effectiveness of the proposed assessment model. A conclusion is made in Section 5.

2. Preliminaries

As mentioned above, there are many math tools to handle uncertain information, such as fuzzy sets theory [3740] and evidence theory [41]. In addition, interval number is also a simple and efficient method. In this section, the basic concepts of interval number are briefly introduced.

2.1. Interval Number

The interval approach is originally developed by Moore [42, 43] in order to keep track of the lower and upper bounds to the exact result when carrying out numerical calculation on digital computers with a finite number of significant digits. Recently, it has been suggested to use interval numbers to represent the uncertainty or probability resulting from instinct of the incomplete knowledge or cognition [4446]. In this paper, interval numbers are used to represent the uncertain value of threat attributes.

Definition 1. Interval number is defined as an ordered pair ; particularly, if , then is obviously equal to a real number, and the interval numbers we discussed in this paper are all positive interval numbers which can be described as , where is positive real number set.

Definition 2. Let , be any two interval numbers; then some of their arithmetic operations are defined as follows:(1) if and only if and ;(2).
The problem of a multiobjective threat assessment can be described with the mathematical model. The object set can be represented with and the attribute set is . Decision matrix is composed of attribute values (interval number , ; ). The threat goals are measured by the attribute . And the attribute is supplied by specific researches according to the different situation.

Definition 3. The attribute of each object is composed of two types, namely, “cost type” and “benefit type.” Let represent their corresponding subscript set, and the “cost type” is a kind of attribute which becomes better with its value increasing while “benefit type” is a kind of attribute which becomes better with its value decreasing. To eliminate the influence of decision results in different physical dimension, is normalized to with the procedure defined as follows:

Definition 4. Let be the set of interval numbers in , and let the distance between two interval numbers and be defined by Tran and Duckstein [47] asIt can be proved that is a distance on . First, . Symmetry is transparent. If , then . The triangle inequality follows from the fact that the function to be integrated in (2) is the square of Euclidean distance.
Although only the lower and upper bound values of the two interval numbers appear in (3), which is derived from (2) for operational purpose, the integral in (2) shows that this distance takes into account every point in both intervals when computing the distance between those two interval numbers. It is different from most existing distance measures for interval numbers which often use only the lower and upper bound values [4850]. Paper [51] proposes a distance measure for intervals which also considers every point of both intervals. Its general form, however, is too complicated and the authors later restricted the measure to a particular case with a finite number of considered values for operational purpose.

2.2. Genetic Algorithms

The simple genetic algorithm (Algorithm 1) over populations defined as multisets consists of individual binary -tuples with fitness values . For the creation of offspring individual in each generation random genetic operators like crossover and mutation are applied to parental individuals which are selected according to their fitness values as follows. The population is initialised appropriately, for example, by randomly choosing individuals in . GA is also used to cope with the issue of decision making by lots of researchers like Chan et al. [52].


(1) Selection. For fitness-proportional selection each individual in the population is selected with a probability which is directly proportional to its fitness . The individual selection probability is then given by

(2) Crossover. The crossover operator takes two selected individual -tuples and and randomly generates an offspring -tuple according to .

For -pointcrossover a position is chosen with uniform probability and the two -tuples and are generated, one of which is randomly chosen as offspring . The 1-point crossover operator is applied with crossover probability . In case of uniform crossover each bit of the offspring is chosen with probability from the parental -tuple and with probability from the parental -tuple or vice versa.

(3) Mutation. In the simple genetic algorithm with binary -tuples the bitwise mutation operator is defined by randomly flipping each bit of the individual -tuple with small mutation probability . A typical value is .

2.3. Interval AHP

Assume that a decision maker provides interval judgements instead of precise judgements for a pairwise comparison matrix. For example, it could be judged that th criterion is between and times as important as th criterion with and being nonnegative real numbers and . Then, an interval comparison matrix can be expressed as where and for all . The above interval comparison matrix can be split into two crisp nonnegative matrices: where .

Let be the principal eigenvalue of , which is an interval number. The principal right eigenvector corresponding to is denoted by , where and . Therefore, the eigenvalue problem can be written as namely, According to interval arithmetic, (8) can be expressed as which is equivalent to the following two equations: The vector of weight is , where


Assume that a MCDM problem has alternatives, , and decision criteria/attributes, . Each alternative is evaluated with respect to the criteria/attributes. All the values/ratings assigned to the alternatives with respect to each criterion form a decision matrix denoted by . Let be the relative weight vector about the criteria, satisfying . Then the TOPSIS method can be summarized as follows.

(1) Normalize the decision matrix using the following equation: where is the normalized criteria/attribute value/rating.

(2) Calculate the weighted normalized decision matrix : where is the relative weight of the th criterion or attribute, and .

(3) Determine the ideal and negative-ideal solutions: where and are the sets of benefit criteria/attributes and cost criteria/attributes, respectively.

(4) Calculate the Euclidean distances of each alternative from the ideal solution and the negative-ideal solution, respectively:

(5) Calculate the relative closeness of each alternative to the idea solution. The relative closeness of the alternative with respect to is defined as

(6) Rank the alternatives according to the relative closeness to the ideal solution. The bigger the , the better the alternative . The best alternative is the one with the greatest relative closeness to the ideal solution.

3. Mathematic Model of Threat Assessment

The mathematic model of threat assessment is detailed in this section. As shown in Figure 1, the input of our model includes two parts, namely, subjective factors and objective factors, and the output of this model is the ranking of all the threat objects.

3.1. Integration of Objective Factors

Due to the imprecision of the sensors in the complex situation, sometimes, we should make use of interval numbers to represent the attributive values of the threat targets. In order to measure the different threat degree, the problem arises regarding how best we can aggregate this interval number into a general interval number. The main idea of our proposed method can be shown as follows: this paper integrates the fuzzy numbers through the distance among interval numbers to find a special interval number whose distance is the optimal value from all other normalized interval numbers with the searching method of genetic algorithms. The binary function is constructed with the summation of all the distance between and each normalized interval number. Hence, the problem is mapped into finding the optimal solution of under the constraint with the method of genetic algorithms. This process can be mathematically described by where is the interval value of attribute of the th object. is the interval distance between and . and are the boundary range of and .

Now the procedure of aggregating interval numbers (objective factors) is made by GA in detail as follows.

(1) Determination of the Encoding Method. The first step is to determine the encoding strategy, in other words, how to represent the data between and with the computer language. In this paper, the binary encoding strategy is adopted to represent the chromosome. The length of chromosome depends on the accuracy of encoding. Assume that the domain of variable is and the accuracy of encoding is the digit after decimal point. The length of a binary string variable can be calculated as follows: where is to represent the precision of the result and means the maximum integer which is not more than .

(2) Determination of the Decoding Method. The decoding of chromosome is to translate the chromosome from binary data to decimal data in the interval number accordingly. Assume that the binary encoding of a chromosome is represented with . Firstly, the binary string should be translated to decimal number , which can be denoted as follows: Then, the final decoding data in the interval number can be calculated as where is the encoding length of the chromosome.

(3) Construction of the Initial Population. According to (18), the encoding length of chromosome for each variable can be obtained. Hence, the total length of chromosome for variable weight can be accumulated to be composed of the length of a single long chromosome. For the point position of each chromosome, the method of generating population can be denoted as follows:

(4) Determination of the Adaptive Function and Adaptive Value. Generally, the adaptive function is designed according to the objective function , and the adaptive function is denoted as . In order to lay the foundation for calculating the selected probability of each individual behind, the optimalizing direction of adaptive function should adapt to the incremental direction of adaptive value. Due to the value of objective function , the adaptive function is defined as In this paper, the objective function is determined by according to (17).

(5) Determination of the Selection Criteria. In this paper, the proportional selection strategy of adaptive value is introduced, and the proportion of every individual is defined as selected probability . Assume that the population whose scale is n is given as and the adaptive of is given as ; then the selected probability is denoted as Then the accumulative probability of every chromosome is denoted as follows: After random data , is generated and the selection of chromosome for a new population can be made if.

(6) Determination of the Genetic Operators. Assume that we have a population including four individuals described as follows:

After conducting several Roulette Wheel tests, assume that chromosome occupies the most area of the whole circle and chromosome occupies the least area of the whole circle. According to the selecting criteria, chromosome is selected to make a reproduction, while chromosome fell into disuse:

The crossover operator in this paper adopts the strategy of a single cutting crossover. This method considers the two flanks of the cutting into two substrings; then the right substrings should be exchanged with each other to get two new individuals. If the crossover probability , it means that 25% of the chromosomes on average exchange each other:

Mutation operator is to change some gene of chromosome with a tiny possibility. If the mutation probability , it means that 1% of all the genes are expected to mutate.

(7) Determine the Process Parameter and Terminal Condition. Through the prior test, the crossover probability is between 0.4 and 0.99, the mutation probability is between 0.0001 and 0.01, the scale of population is between 20 and 100, and the terminated condition may be determined by the precise iterated generation or determined by the minimum bias , which satisfies where is the maximum fitness value and is the objective fitness value. If the judgement condition is not satisfied with the terminal condition, then go to step (4).

Example 1. Three interval numbers are as follows:Assume that the importance of each interval number is the same. First, we suppose that the integrated interval number of the three fuzzy numbers is represented with , . So the distance between and each number can be shown as follows according to the distance of the interval numbers:Hence, the integration of three interval numbers is mapped into solving the minimum of the formula under the constraint of .
Assume that the accuracy of encoding is 0.0001, the crossover probability is 0.2, and the mutation probability is 0.05. First, the length of chromosome can be calculated as . Then the encoding length of and is 32. The adaptive function is . After 300 times of repeat, the convergent result is (5.4949, 7.5023) and the minimal value of the objective function is 4.4315.

3.2. Integration of Subjective Factors

For the subjective factors, the interval AHP (interval eigenvector) is applied to deal with the object’s threat weight according to experience of the commanders/experts. A numerical example is used to describe the process as follows.

Example 2. Suppose that we get a pairwise comparison matrix through interval judgements which is denoted as Then the interval comparison matrix can be split into two crisp nonnegative matrices: Then the principal eigenvalue of can be calculated: , where and . The principal right eigenvector corresponding to is , where So, the vector of weight is where and .

3.3. Integration of Objective and Subjective Factors

In this part, the method of discounting will be used to integrate the objective and subjective factors. Through the two procedures above, for each threat object, the integrative interval number describing the objective factor can be integrated by the convergent process of genetic algorithm. Hence, assume that there are threat objects; optimal interval number will be used to represent the different threat degree; they can be formally denoted as

Sometimes, the commanders need to evaluate the threat degree of each target according to their previous experience in order to control the situation. Hence, through constructing the interval decision matrix and calculating the interval eigenvector, the subjective interval numbers representing the evaluation of experts can be described as

Then the fusion weight of all the threat targets can be calculated according to the method of discounting the objective and subjective factors, which is denoted as where . In the process of decision, the discounting factors and may be modified to adapt the importance of objective factors comparing with subjective factors. If the influence from experts’ subjective evaluation is dominant over the influence from the sensors’ detection in some complex situation, then the discounting factor should satisfy , .

3.4. Ranking the Threat Objects

In this part, the method of TOPSIS will be made use of to deal with the ranking of the threat objects according to their different threat degree. The classical TOPSIS method is a technique for order preference by similarity to ideal solution. The ideal solution (also called positive ideal solution) is a solution that maximizes the benefit criteria/attributes and minimizes the cost criteria/attributes, whereas the negative ideal solution (also called anti-ideal solution) maximizes the cost criteria/attributes and minimizes the benefit criteria/attributes. The so-called benefit criteria/attributes are those for maximization, while the cost criteria/attributes are those for minimization. The best alternative is the one which is closest to the ideal solution and farthest from the negative ideal solution.

Because the normalized interval numbers all distribute in , the two extreme bounds can be considered as the anti-ideal solution and the ideal solution. That is to say, the interval number can be the anti-ideal solution and the interval number can be the ideal solution.

Then the distance between the ideal solution and each target’s weight from the fusion process of integration with the method of discounting can be calculated. It is not difficult to know that the target nearest to the ideal solution is the most dangerous. Hence the ranking of threat objects can be made easily according to their distance to the ideal solution or the anti-ideal solution.

Example 3. Assume that an interval number is (data from integrating objective and subjective factors); then the distance (ideal distance) can be calculated with the interval distance formula ; similarly, the distance (anti-ideal distance) also can be calculated with interval distance ; at last, the relative closeness can be calculated as

4. Numerical Examples

4.1. The First Numerical Example

In this part, a numerical example [53] is used to illustrate the proposed threat model. Suppose that assessments of the threat objects are all determined with the objective factors (). The interval decision matrix can be denoted as in Table 1; there are three alternatives (, , ) and two criteria (, ). The decision maker wants to choose an alternative that has minimum and maximum . Now we order the three alternatives with the proposed model above.

First, the interval decision matrix can be normalized through (1). The normalized decision matrix can be denoted as in Table 2.

Second, through the proposed objective aggregating method with GA, the aggregated results can be converged as in Table 3.

At last, the priority of the three alternatives can be determined by TOPSIS, and the result can be shown as in Table 4. The ranking of three alternatives is , which is the same as the result of the paper [53].

This numerical example illustrates the effectiveness of the proposed model without considering the subjective factors. The aggregation of the objective factors is always optimum with GA. In the following part, another numerical example is presented to illustrate the effectiveness of the proposed threat model with considering both objective factors and subjective factors.

4.2. The Second Numerical Example

In this part, a real application in aerial object threat will be introduced to describe the operational procedure of threat assessment. The threat objects and their attributes can be shown as in Table 5.

In Table 5, , and mean “type of target,” “speed of target,” “time of arrival,” and “shortcut of airway.” At the same time are benefit attributes, are cost attributes, and the numerical value is described in Table 6.

In order to eliminate the influence resulting from different measures, the normalization can be made through (1). The normalized result of Table 5 can be described as in Table 7.

Through the objective fusion method by GA, the integrated results can be converged as in Table 8.

In the following part, we will take the commanders’ experience into consideration. The interval decision by experts for the five threat objects is described as

As for the subjective factors, the method of interval eigenvector is used to generate the interval threat weights of all targets. The fusion process of subjective factors is like Example 2, and the fusion weight of threat objects by experts’ decision is

Hence, at present, the integration between objective factors and subjective factors can be made by the method of discounting. It can be denoted as in Table 9.

According to the idea of TOPSIS, we can get the ideal distance and anti-ideal distance of each target’s fusion weight, and the ideal interval number is and the anti-ideal interval number is , so we can get the relative closeness of each target described as in Table 10. According to the relative closeness , the ranking of all threat targets is (threat degree from large to small). This result is shown by Figure 2 explicitly.

As presented above, the proposed threat model can not only be applied to the objective issue, but can also be easily extended to considering the subjective issue from knowledge of the domain experts. When the objective discounting factor and the subjective factor , the proposed threat model can be absolutely applied to the objective issues. When the objective discounting factor and the subjective factor , the proposed threat model can be absolutely applied to the subjective issues. Otherwise, the proposed threat model is applied in both subjective and objective issues according to the optimum discounting factors ( and ).

5. Conclusion

Threat assessment is very important in many fields. This paper proposes a new model of threat assessment based on interval number. Both subjective and objective factors are taken into consideration in this threat assessment model. GA is used to determine the objective factors. In addition, interval AHP is applied to determine subjective factors. Then a method of discounting is proposed to integrate the subjective weight and objective weight. At last, the idea of TOPSIS is adopted to rank the objects according to their threat degree. A numerical example and a real application are used to illustrate the effectiveness of the proposed threat model.

Conflict of Interests

The author declares that there is no conflict of interests regarding the publication of this paper.


The author greatly appreciates the reviewers’ valuable suggestions. This work is partially supported by National High Technology Research and Development Program of China (863 Program) (Grant no. 2013AA013801), National Natural Science Foundation of China (Grant no. 61174022), Specialized Research Fund for the Doctoral Program of Higher Education (Grant no. 20131102130002), R&D Program of China (2012BAH07B01), and the Open Funding Project of State Key Laboratory of Virtual Reality Technology and Systems, Beihang University (Grant no. BUAA-VR-14KF-02).