Abstract
The current paper presents a soft-sensor method based on belief-rule-base (BRB) system for solving the problem of tipping paper permeability measurement in the tobacco industry. Firstly, BRB is utilized to establish a model between the feature variables in the tipping paper image and the corresponding paper permeability obtained by the traditional measuring device. Unlike the traditional case of BRB, this paper adds the output attribute as the optimization parameters. In this way, the feasible solution space can be enlarged to obtain an effective BRB model. Second, in order to find the reasonable parameters of BRB in a complex nonconvex solution space, an enhanced differential evolutionary (DE) algorithm is developed to train BRB, which not only embeds a simplex method to stress the balance between the global and local search but also designs a perturbation operation and an adaptively selected mutation strategy to maintain the diversity of search direction. The test results and comparisons based on the data collected from a cigarette factory in China show that the presented method is effective and robust.
1. Introduction
Cigarette smoke contains tar, nicotine, and other harmful substances. The research on the technology in reducing the tar and harmful components of cigarette smoke has been the important research and application field in the tobacco industry since the 1970s. Now many countries have made strict restrictions on the content of tar, nicotine, and CO in commercial cigarettes. There are some existing tar reduction techniques, such as the use of tobacco with low tar, the improvement of tobacco formula and cigarette structure, the addition of expanded tobacco, and the application of tipping paper. Since the perforations on tipping paper are helpful for introducing air into cigarette smoke to dilute tar, the use of perforated tipping paper to package cigarette butt has become one of the most common methods for reducing tar. Considering the air permeability of tipping paper is directly related to the dilution of smoke and have greatly influenced on the cigarette taste, the tobacco industry must strictly control the tipping paper permeability [1, 2]. As people’s requisition for health and cigarette quality is increasing constantly, the permeability measuring technology is becoming more and more important and valuable.
The tipping paper permeability is a measurement of the amount of air per unit of time that flows through the vertical given area of tipping paper under a certain pressure condition, whose unit is (). A representative of existing detection devices is BASTAN Company’s online air permeability tester. It utilizes the conversion of the photoelectric signal to obtain a measurement value of permeability. Nevertheless, the existing online testing equipments are very expensive and complicated to repair, which have become a burden on enterprise. Therefore, the development of a more economical, simple, and effective detection system is necessary for enterprise.
Recently, with the technological advances in CCD sensors and integrated semiconductors, industrial machine vision cameras can deliver great features and functionality at outstanding price-performance ratios. In fact, the sizes of the perforations on tipping paper directly affect the permeability. In addition, the surface roughness and the thickness of tipping paper also have some influence on the permeability. These factors can be utilized to build a model for measuring it. Figure 1 shows a cigarette with a wrapping of perforated tipping paper on its filter, where the perforations have small diameters (30–90 μm). The gray image of the corresponding tipping paper is given in Figure 2. The light points in this image are the perforations produced by a laser-beam drilling machine. The area of perforations can be calculated by using image processing technique, and the surface roughness as well as the thickness of tipping paper can be reflected by the gray value of tipping paper’s image. Thus, we are motivated to design a permeability measuring device based on machine vision system. That is, an industrial camera with microlens can be utilized to capture the image of tipping paper with microperforations in a certain small region, and a soft-sensor model can be built between the intermediate variables and the tipping paper permeability. Since a high-quality camera and other hardwares are easy to buy directly from manufacturers, the key issue for this new device is how to develop an effective soft-sensor method.


BRB is a kind of expert system in nature, which can effectively use various types of information to establish the nonlinear model between input and output. The belief-rule-base concept and its inference methodology were proposed by Yang et al. [3], which are based on the evidential reasoning approach (ER) [4]. Compared with traditional rule-based systems, a BRB system can not only be used to analyze the decision problems by handling both quantitative data and qualitative information with uncertainties but also can capture more complicated nonlinear causal relationships between antecedent attributes and consequents. By now, it has been successfully used in the safety analysis of offshore systems [5, 6], leakage detection of oil pipeline [7, 8], product life assessment [4], evaluation of gastric cancer [9–11], hidden behavior prediction [12], safety assessment for complex system [13, 14], and information fusion [15]. Due to BRB’s wide applications in a variety of fields, the current paper presents a novel soft-sensor method based on BRB for measuring the permeability of tipping paper.
In the traditional methods, BRB systems are built by the decision maker according to the experience and knowledge or other original model. However, for complex systems, it is often unable to use the traditional methods to establish an accurate analytical model. A BRB represents functional mappings between inputs and outputs. Compared with the traditional IF-THEN rule, it provides a more informative and realistic knowledge expression. In BRB, the input data, the belief structure, the attribute weight, the rule weight, the activation weight, and all the activation belief rules are combined to generate the corresponding outputs. Thus, it can be seen that these BRB parameters are important factors to affect the performance of the BRB system. In the traditional methods, the experts set the parameters of BRB according to their experience and knowledge. However, for the complexity and uncertainty of many decision-making problems, in reality, it is often unable to use the traditional method to establish an accurate model, which limits the ability of BRB to simulate the actual system. Therefore, in order to improve the modeling precision of a complex system, it is necessary to train and optimize the parameters of BRB.
The parameters’ optimization or training for BRB is a nonlinear constrained optimization problem (NCOP), whose objective function is a complex nonlinear formula. From equations (4)–(11), it can be seen that the objective function is a multivariate compound one with strong nonconvexity. Since there is still no general method for obtaining global optimum of NCOP, the NCOP of BRB is usually solved by using an effective gradient-based solver (i.e., ) in MATLAB. Yang et al. [16] first proposed the generic BRB learning model and a parameter optimization method by using . An essential feature of is sensitively dependent on initial conditions, which determines that is applicable to small-scale optimization problem but is not suitable for solving complex problems. After that, Zhou et al. [17] proposed an online update algorithm based on expectation maximization estimation algorithm, but the effectiveness of the algorithm is related to the completeness of the data and the set of referenced values. Subsequently, Chen et al. [18] added the premise attribute as the new training parameters and used to optimize the system. Chang and Zhang [19] developed a self-learning parameter training algorithm based on the optimization step and gradient method. The accuracy of BRB is improved, but the algorithm has the defects of large implementation complexity and low convergence. Chang et al. [20] first used standard differential evolution algorithm to train the parameter learning model, which can effectively improve the efficiency and precision of BRB. However, all the above papers focus on the traditional case of BRB with predetermined referenced values of the input and output attributes. From an optimization point of view, the more adjustable parameters or decision variables the BRB model has the larger area of feasible solutions it owns. That is, the BRB model with more decision variables is helpful in finding a more reasonable model by using optimization methods. Moreover, in many real-life applications, not all the expert knowledge is explicit, which means the overuse of it may decrease the accuracy of BRB. Thus, how to model an effective BRB with less expert knowledge or predetermined parameters is an important issue. Therefore, this paper adds the referenced values of input and output attributes as the optimization parameters of the BRB system. Obviously, this new BRB model is more flexible than the traditional one and the latter is just a subset of the former.
Traditional mathematical approaches are usually applied to optimize the parameters of BRB, which are guaranteed to converge only under certain convexity assumptions and are sensitive to the initial solution. Unfortunately, the BRB model is nonconvex and has more than one local optimal solution. That is to say, mathematical approaches are limited to solving its simple and small-scale cases. When the referenced values of both input and output attributes are added as the optimization parameters, the NCOP of BRB becomes more complex and its feasible solution region is enlarged. Under such circumstances, it is difficult for the traditional approaches to obtain high-quality solutions within a reasonable time. Recently, some intelligent algorithms have been proposed to deal with different complex optimization problems, such as differential evolutionary (DE) algorithm and estimation of distribution algorithm (EDA) for flow shop scheduling problems [21–24], particle swarm optimization (PSO) algorithm for multiobjective optimization problem [25], DE and PSO for parametric identification problem [26], and evolutionary algorithm (EC) and genetic algorithm (GA) for interval multiobjective optimization problems [27, 28]. These algorithms utilize their own evolutionary mechanisms to perform an effective search in complex solution space and obtain very satisfactory results. Among them, differential evolutionary (DE) algorithm is one of the most successful evolutionary optimization techniques, which uses a simple operator to create new candidate solutions and a one-to-one competition scheme to select the new candidate [29]. Although having some advantages in finding the global or near-global optimum, the local search ability of DE is relatively weak and sometimes the searching process of DE may be trapped in local minima [30]. NM simplex method (NMS) was first introduced by Nelder and Mead [31] and is one of the most widely used local direct search methods for nonlinear unconstrained optimization. NMS is a derivative-free line search method that was specially designed for solving traditional unconstrained problems of minimization type, like nonlinear least-squares problem, non-linear simultaneous equations, and function minimization. It has been successfully applied to the low dimensional unconstrained optimization problem [32]. Thus, in this paper, an enhanced DE is proposed for the NCOP of extended BRB. The contribution of this paper can be outlined in the following aspects: (1)The presented BRB model adds the referenced values of input and output attributes as the optimization parameters to construct a more flexible model, which is beneficial to obtaining a more accurate model in an extended solution space(2)To avoid the drawback of traditional mathematical methods, a hybrid differential evolutionary algorithm, namely MNMSDE, is proposed to efficiently determine the extended model’s parameters, which combines DE, NMS, perturbation operation, and an adaptive Meta-Lamarckian learning strategy for improving the global and local search ability(3)A soft-sensor method based on extended BRB model and MNMSDE is presented, whose effectiveness is verified by the data collected from a Chinese cigarette factory
The structure of this paper is as follows. In the second part, BRB is described. The third part introduces the optimization problem and the proposed MNMSDE algorithm. The fourth part is the experimental analysis. Finally, the work in this paper is summarized in the last part.
2. Problem Description
2.1. Belief Rule Base
A BRB consists of a series of belief rules, and the expression of the kth rule is as follows: with a rule weight and attribute weights , where denotes the kth rule of BRB. Each rule is provided with a rule weight , which indicates the importance of it relative to other rules. is the input of BRB. represents the referenced values of the prerequisite attribute of the kth rule. represents the number of prerequisite attributes. represents the number of belief rules in BRB. represents the weight of , which reflects the importance of relative to the other premises. is the referenced value of the output in the kth rule. is the belief degree assessed to .
2.2. The Reasoning Mechanism of BRB
Belief-rule-base inference methodology is based on ER. ER can be directly applied to combine activated belief rules and generate final conclusions as follows. First, transform the inputs into the probability mass corresponding to the referenced values and generate a series of . Then, aggregate all activated rules to generate the combined degree of belief in each possible consequent by ER. Finally, the output is generated by aggregating all and , which can be represented as follows:
Detail reasoning process is as follows. (1)Transform Input. Using the information transformation technology, the matching degrees of the inputs relative to each reference level can be described as where is the belief degree assessed to input . The calculation method is as follows: (2)Calculate the Activation Weight. is the activation weight of , which measures the degree to which the kth rule is weighted and activated. The calculation formulas of the activation weight of the kth rule can be calculated by the following formula: where is the relative weight of the kth rule, and represents the ith antecedent attribute weight of the kth rule. If >0, the kth rule is activated. is generally set to 1.(3)Aggregate the Belief Rules Using the Evidential Reasoning Approach. Using ER to combine all the rules and the following output can be obtained as follows: where is calculated by equation (7), and is the belief degrees relative to the evaluation result .(4)Calculate the Expected Utility. Assuming the utility of evaluation result is , then the expected utility of can be calculated by the following formula:
3. Optimization Model and Algorithm
3.1. Optimization Model of BRB
are the adjustable parameters in BRB. The objective of the training is to minimize by adjusting the parameters . is 1 in this paper.
The processes of BRB system parameter optimization are shown in Figure 3. The optimization model for training BRB parameters is defined as follows:

In order to obtain a more stable and effective BRB, some related optimization methods have been put forward [5, 8, 16], mainly for the optimization of parameters and structures. Predecessors had taken the optimization of the belief degree, the prerequisite attribute weight, and the rule weight into account, but the referenced values are usually given by experts according to their experience. Each belief rule is obtained by adding the belief degrees, the prerequisite attribute weights and the rule weights to a traditional IF-THEN rule, whose parameters are largely influenced by their corresponding referenced values. This means the relatively reasonable referenced values are conductive to build a more accurate BRB model. Therefore, in this paper, we consider adding all referenced values as the parameters to be optimized.
Since the hole area and the average gray value in an image largely determine the perforated tipping paper permeability, we chose these two factors as inputs and the air permeability as output to establish a BRB. Denote as the reference level of the hole area, the reference level of gray value, and the reference level of air permeability.
Here is an example of the evaluation level of the output. and represent the lower and upper bounds of each referenced value of the output, respectively. When the referenced values are trained, the new referenced values are generated in the referenced value interval as follows: where and the referenced values are arranged in sequence according to the order of small to large.
The column vector of the parameters to be optimized is as follows: where and the dimension DM of this problem is .
In (11), , and represents the expected output utility of the BRB system. The smaller the is, the more accurate the BRB system can simulate the actual system. According to [4], the definitions of constraints are given as follows: (1)The rule weight is changed between 1 and 0, i.e., (2)The belief degree shall not be greater than 1 or less than 0, i.e., (3)If the rule is complete, the rule outputs part of confidence is equal to 1; otherwise, it is less than 1, i.e., (4)represents the jth referenced value of , and it should be in the range of input data, i.e., where represents the lower limit of the input range and represents the upper limit of the input range.(5) represents the jth referenced value of , and it should be in the range of output data, i.e., where represents the lower limit of the output range and represents the upper limit of the output range.
3.2. MNMSDE for BRB Model
3.2.1. Differential Evolution Algorithm
Differential evolution algorithm is an intelligent algorithm that is based on population evolution, and it solves the optimization problem by the cooperation and competition among the individuals in the population. In this paper, the objective function is .
Denote as the generation, the individual at generation , and the mutant vector relative to . The initial population is randomly generated by the following formula:
The three core operations of DE are as follows.
Mutation operation: DE creates a mutant vector by employing the following mutation operation: where is randomly selected, and . is the scaling factor that is usually less than 1.
Crossover operation: after mutation, the crossover operation is applied to replace partial variables of with the corresponding variables of to generate a trial vector by the following formula: where is the jth evaluation of a random number uniformly distributed in the range of [0, 1], is a randomly chosen index from the set , and is the crossover probability.
Selection actions: the greedy selection method is used to select the better individual between and for creating the next generation, which is shown as follows:
3.2.2. Adaptive Meta-Lamarckian Learning Strategy
In the recent years, methods for automated tuning of mutation strategies attract increasing attention [33–35], it has been verified that the application of the multimutation strategies can enhance the performance of DE. Ong and Keane [36] proposed an adaptive Meta-Lamarckian learning strategy to dynamically decide which neighborhood to be chosen for local search, which is conducive to improve search performance and reduce the probability of utilizing an inappropriate local search method. This learning strategy widely is used in local search [37, 38]. Inspired by the work of Ong and Keane, we dynamically determine the selected probability of each mutation strategy, which can gradually evolve the most suitable learning strategy at different learning stages. That is, we select three mutation strategies as candidates: “DE/best/1”, “DE/rand/1” and “DE/rand to best/1”. In the initial stage of evolution, the selected probability of each mutation strategy is set to the same value. As evolution continues, is updated at every generations. In each generation, the mutation strategy adopted by each trial individual is determined by the roulette wheel rule. Denote the population size, the total times of be selected at generation , the total success times of at generation , and the success probability of . Summing up the above symbols, we present the procedure of calculating as follows:
Step 1. Set , , , and for ;
Step 2. Set and ;
Step 3. Set and apply a roulette wheel rule based on to select one strategy for the individual ;
Step 4. If , then
if and the trial individual is better than , then
if , then ;
if and the trial individual is better than , then
if , then ;
if and the trial individual is better than , then
Step 5. If , then go to Step 3;
Step 6. If , then
for ,
for ,
and , ;
Step 7. If , then ;
Step 8. If , then go to Step 2.
Step 7 is used to prevent a mutation strategy from being completely eliminated in the process of evolution.
3.2.3. NM Simplex Method
The NM simplex method forms a simplex by setting individuals as vertices, where is the size of the vertices. Through a series of geometric operations (i.e., reflection, expansion, contraction, and reduction), a new vertex can be created. After each geometric operation, the current worst vertex is replaced by the new vertex if the latter is better than the former. With the help of these geometric operations, the simplex can improve itself and come closer to the optimum. The steps of the simplex method are concluded as below:
Step 1. Set and ;
Step 2. Sort and renumber all the individuals by ;
Step 3. Reflection: calculate the reflection point as follows:
where the geometric center and the reflection coefficient is greater than 0;
if , then
and go to Step 7;
if , then go to Step 4;
if , then go to Step 5.
Step 4. Expansion: generate the expansion point as follows:
where the expansion coefficient is greater than 1;
if , then
otherwise,
Go to Step 7;
Step 5. Outside contraction: if , then produce the external contraction point as follows:
where is the contraction coefficient between 0~1;
if , then
and go to Step 7;
otherwise,
go to Step 6;
Inside contraction: if , then produce the internal contraction point by the following formula:
if , then
and go to Step 7;
otherwise,
go to Step 6;
Step 6. Shrink: calculate new vertices as follows: where is the collapse coefficient between 0~1;
Step 7. If , then go to Step 2;
Step 8. Output .
3.2.4. MNMSDE
For the n-dimensional optimization problem, the initial population is generated randomly. At each generation, the new individuals are generated by using the standard DE’s operations with the adaptive Meta-Lamarckian learning strategy in Subsection 3.2.3, and then the NM simplex method is applied to the new population for stressing local search. The details are given by the following steps:
Step 1. Set ;
Step 2. Reindex the individuals in the new population in the increasing order of their fitness values;
Step 3. Select the worst individuals as the vertexes to perform the NM simplex method in Subsection 3.2.3. In this stage, we use the best individuals in the population to find the geometric center to guide the simplex search. is calculated by the following formula:
Step 4. Replace (i.e., the worst individuals) with .
For the purpose of increasing the diversity of the population and preventing the algorithm from falling into local optimum, we add a perturbation operation in the algorithm. That is, at the beginning of each generation, the standard deviation of fitness function values is calculated as a criterion for judging the degree of individual differences. If the standard deviation is less than the set value , randomly select half individuals of the population and update them by the following formula: where is a randomly chosen index from the set , represents the index of the selected individuals, is a randomly selected individual, and is the best individual at generation .
In addition, when any individual violates the constraints of BRB’s corresponding NCOP, its corresponding elements (see Equation (13)) should be repaired by the following equations:
On the one hand, the local search and global search ability get the balance due to the combination of the geometric search based on NMS and the evolutionary search based on DE. On the other hand, as NMS uses the center of the good individuals in the population to implement geometric search, the local search ability will be enhanced.
The MNMSDE algorithm’s visual diagram is shown in Figure 4.

4. Experimental Study
4.1. BRB Model Construction
The number of referential points of each input and output determines the size of the rule base, a reasonable size of the BRB is conducive to the accuracy of the results. In this paper, for each input and output attributes, 5 referenced values are selected as follows:
Suppose we choose the following referenced values:
According to the above-referenced values, the prerequisite attribute of input of each rule is calculated by formulas (4)–(6), and the belief degree of its corresponding input and output is given by the expert. Then the tipping paper permeability BRB detection model can be constructed. The rule base is shown in Table 1. As there are 25 belief rules in BRB and 5 referenced values for the two inputs and the corresponding output, the total number of optimization parameters is 159. Precisely speaking, the optimization parameters include the rule weight , the belief degree , the referenced values of input attribute , and the output attribute .
4.2. Experimental Results and Comparisons
In order to verify the effectiveness of the proposed model and method, we used 1000 data samples, in which 800 samples are randomly selected for training and the rest of the samples are used for testing. The training and testing samples are collected from a Chinese cigarette factory. Each sample includes the hole area, the average gray value, and the corresponding permeability on a unit area. The former two are the input values, and they are calculated by using the image captured with an industrial camera of our designed device. The last one is the output value, which is obtained via traditional instruments. All compared algorithms are coded in MATLAB, and experiments are executed on an Intel 2.6GHz PC with 8GB RAM. To make a fair comparison, each algorithm is independently run 30 times under the same evaluation times.
For the purpose of demonstrating the effectiveness of MNMSDE, we carry out some comparisons with [16] and standard DE [20]. The mutation strategy used in standard DE is “DE/rand/1”. Moreover, two variants of MNMSDE are also compared, whose abbreviations are as follows: (1)MNMSDE_V1. Standard DE with Meta-Lamarckian learning strategy(2)MNMSDE_V2. MNMSDE_V1 with NM simplex method
The traditional BRB model optimized by , standard DE, MNMSDE, MNMSDE_V1, and MNMSDE_V2 are denoted as F-TBRB, DE-TBRB, MNMSDE-TBRB, MNMSDE_V1-TBRB, and MNMSDE_V2-TBRB, respectively. Similarly, the new BRB model optimized by , standard DE, MNMSDE, MNMSDE_V1, and MNMSDE_V2 are denoted as F-BRB, DE-BRB, MNMSDE-BRB, MNMSDE_V1-BRB, and MNMSDE_V2-BRB, respectively.
Based on the previous experiments, parameters of MNMSDE and its variants are suggested as follows: (1)The population size can be set to [50, 80], which is useful in ensuring the diversity of the population(2)To keep the width of global search, the scaling factor and the crossover probability can be set to [0.5, 0.7] and [0.8, 0.9], respectively(3)To enhance the depth of local search, the reflection coefficient , the expansion coefficient , and the contraction coefficient can be set to [0.8, 1.1], [1.6, 2.2], and [0.4, 0.6], respectively(4)The evaluation times can be set to [42,000, 45,000], which is helpful in getting a reasonable model
In the experiment, DE, MNMSDE_V1, MNMSDE_V2, and MNMSDE use the same parameters as follows: , , , , , and the evaluation times are 44,000.
Let MaxRE denote the maximum relative error, MinRE the minimum relative error, ARE the average relative error, RMSE the root mean square error, RSD the relative standard deviation, and ASD the absolute standard deviation. Test results are shown in Figures 5–17 and Tables 2 and 3, the optimized model parameters in Table 1 and (41)~(43) are displayed in the attachment.













Figures 5 and 6 show the permeability detection or measurement results of original BRB and MNMSDE-BRB, respectively. From Figures 5 and 6, it is shown that MNMSDE-BRB has better fitting and predicted precision than original BRB. That is to say, the special designed DE (i.e., MNMSDE) with the new BRB model has the ability to build a more reasonable model for the considered problem.
The relative errors are listed in Table 2. From Table 2, it can be seen that the new BRBs with referenced values optimized perform better than the traditional BRBs (i.e., TBRBs). This manifests the validity and rationality of adding the referenced values of both input and output attributes as the optimization parameters of BRB. As for DE-BRB, MNMSDE-BRB, MNMSDE_V1-BRB, and MNMSDE_V2-BRB, the sequence (DE-BRB, MNMSDE_V1-BRB, MNMSDE_V2-BRB, MNMSDE-BRB) can be obtained in the decreasing order of their corresponding ARE/MaxRE/MinRE values, which demonstrates the effectiveness by incorporating the Meta-Lamarckian learning strategy, NM simplex method, and the perturbation operation into standard DE. Moreover, Figures 7–16 show the absolute and relative errors of testing samples for the BRBs and TBRBs trained by different optimization algorithms, and Figure 17 shows the evolutionary curves of DE, MNMSDE_V1, MNMSDE_V2, and MNMSDE for training the new BRB model. It can be seen from Figures 7–17 that MNMSDE-BRB provides the best fitting precision among the compared methods. So, it is concluded that our proposed MNMSDE-BRB are more effective than the other compared BRB-based methods.
In order to reflect the dispersion degree of the testing samples and the corresponding estimated permeability values, the RSD, ASD, and RMSE values of the compared methods are given in Table 3. The smaller the RSD/ASD/RMSE is, the more stable the soft-sensor method is. So, it is concluded from Table 3 that MNMSDE-BRB is a robust method.
To sum up, MNMSDE-BRB is an effective and robust method for tipping paper permeability measurement.
5. Conclusion
With a view to modelling an effective BRB with less expert knowledge for the tipping paper permeability measurement problem, this paper proposes a novel soft-sensor method based on BRB and DE. Simulations and comparisons based on real data of a China tobacco factory demonstrate the effectiveness and robustness of the proposed method.
There are two main features of the proposed method. First, in order to improve the reasoning or fitting ability of the traditional BRB model and reduce the dependence on the expert knowledge, the referenced values of both input and output attributes are added as the optimization parameters of the traditional one, which can construct a more flexible model for the considered problem. Second, a hybrid DE algorithm, namely, MNMSDE, is designed to find the reasonable parameters of the BRB model in a more complex solution space, which incorporates the NM simplex method, the perturbation operation, and the Meta-Lamarckian learning strategy into a standard DE to enhance its search ability. That is, we consider the improvements of both the accuracy of the BRB model and the effectiveness of the optimization algorithm.
The proposed MNMSDE-BRB has been used in our designed permeability measuring device. This new device has so far been applied to a famous tobacco factory in China and proved to be quite effective and economical. In our future work, we will generalize MNMSDE-BRB to the diagnosis problem of respiratory disease, which needs to dynamically change the size of the belief rule in accordance with the actual situation and remove the noise in the data.
Data Availability
The partial data used to support the findings of this study are available from the corresponding author upon request.
Conflicts of Interest
The authors declare that they have no conflicts of interest.
Acknowledgments
This paper was partially supported by the National Science Foundation of China under Grants 51665025 and 71601180, and the Applied Basic Research Foundation of Yunnan Province under Grant 2015FB136.