# Neutrosophic Number Optimization Models and Their Application in the Practical Production Process

Authors: Angyan Tu; Jun Ye; Bing Wang
Journal: Journal of Mathematics
 (2021)

---

## Abstract

In order to simplify the complex calculation and solve the difficult solution problems of neutrosophic number optimization models (NNOMs) in the practical production process, this paper presents two methods to solve NNOMs, where Matlab built-in function “fmincon()” and neutrosophic number operations (NNOs) are used in indeterminate environments. Next, the two methods are applied to linear and nonlinear programming problems with neutrosophic number information to obtain the optimal solution of the maximum/minimum objective function under the constrained conditions of practical productions by neutrosophic number optimization programming (NNOP) examples. Finally, under indeterminate environments, the fit optimal solutions of the examples can also be achieved by using some specified indeterminate scales to fulfill some specified actual requirements. The NNOP methods can obtain the feasible and flexible optimal solutions and indicate the advantage of simple calculations in practical applications.


---

## Body

## 1. Introduction

Traditional inventory models [1–4] and production planning models [5–7] involve deterministic constrained functions and/or objective functions in deterministic environments. Nevertheless, uncertainty is nearly universal in real world. Therefore, many uncertain optimization methods were proposed for optimization problems with uncertain variables, interval numbers, stochastic, and fuzzy logics [8–15]. In many applied fields, such as management, engineering, and design problems, uncertain programming has been broadly carried out so far. In order to obtain the optimal crisp values of the objective function and the optimal feasible crisp solutions of the decision variables, the constrained functions and/or objective functions are usually changed into some crisp or deterministic programming problems in existing uncertain programming approaches. So, the aforementioned transformed methods are not really meaningful indeterminate approaches because the real indeterminate optimization problems can only indicate indeterminate solutions rather than optimal crisp solutions in indeterminate environments. Nevertheless, indeterminate programming problems imply the corresponding indeterminate optimal values of the objective function and indeterminate optimal solutions for the decision variables under indeterminate environments. So, it is necessary to find some fit optimization approaches for dealing with indeterminate programming problems with indeterminate solutions.Smarandache [16–18] is a pioneer of indeterminacy theories which provide the new minds to solve indeterminacy problems. He adopted the imaginary value denoted by I and then introduced a neutrosophic number (NN) z = x + yI for x, y ∈ R (R: the set of all real numbers) composed of the determinate part x and indeterminate part yI. As for describing indeterminate and incomplete information, obviously, NNs in the indeterminacy theories are a useful mathematical tool. With the development of indeterminacy theories, NNs were also applied to fault diagnosis [19, 20] and decision making [21, 22] under indeterminate environments.Further, thick function or interval function named neutrosophic function, neutrosophic precalculus, and neutrosophic calculus were provided by Smarandache [23] in 2015, where thick function e: S ⟶ E(S) (E(S) is the set of all interval functions) as the form of an interval function e(x) = [e1(x), e2(x)]. The indeterminate function was applied in engineering problems successfully. For example, Ye et al. [24, 25] and Chen et al. [26, 27] proposed expressions of neutrosophic function and applied NNs in analyzing the joint roughness coefficient. Later, Ye [28] used neutrosophic linear equations of NNs to solve traffic flow problems.At present, neutrosophic linguistic numbers, hesitant neutrosophic linguistic numbers, and their aggregation operators were applied to multiattribute decision making [29–31].But in real situations, affected by each kind subjective and objective reasons, such as absences of precise information judged by decision makers or experts, loss of data, and measurement errors, there exist some indeterminate problems. As for the concepts of NNs, NN functions containing indeterminacyI can represent the indeterminate problems with partial certainty and partial uncertainty under indeterminate environments. Ye [32] and Jiang and Ye [33] introduced NN nonlinear and linear programming models and their preliminary solution methods. However, existing methods for solving complex NN optimization problems imply some difficulty and calculational complexity in their solution process. Inspired by the previous solution methods, this paper first selects the models of practical applications in production process, such as inventory models and production planning models. Then, NN nonlinear and linear mathematical models and their solution methods (Matlab built-in function “fmincon()” and operations of NNs) are built with indeterminacy I as our preliminary application study. Finally, real examples of NN linear programming (NN-LP) and NN nonlinear programming (NN-NP) problems illustrate the feasibility of the proposed methods. The advantage of the proposed methods is that the optimization calculations are simple and effective in practical applications.The remainder of this paper is organized as follows. Section2 depicts some concepts and their operations of NNs. Section 3 first introduces NN-NP problems with an inventory mathematical model and model formation and then uses two methods (Matlab built-in function “fmincon()” and operations of NNs) to solve the NN-NP problems in indeterminate setting. Section 4 presents NN-LP problems with the production planning mathematical model and model formation and then applies two methods regarding the Matlab built-in function “fmincon()” and operations of NNs to solve the solutions in the NN-NP problems and to show the simplicity and effectiveness of the proposed NN-LP methods. Conclusions and future research are provided in Section 5.


## 2. Mathematical Preliminaries

### 2.1. Some Concepts and Their Operations of Neutrosophic Numbers (NNs)

The concept of NN was first proposed by Smarandache [34, 35], which consists of two parts (a determinate part and an indeterminate part). He defined the mathematical expression form z = x + yI for x, y ∈ R, where R represents all real numbers and I is indeterminacy. So, it is conveniently used in indeterminate environments.For example, consider that a NN isz = 13 + 5I. Then, its determinate part value is 13 and its indeterminate part value is 5I. When I ∈ [0, 0.5], it is equivalent to z ∈ [13, 15.5] for sure z ≥ 13.Letz1 = x1 + y1I and z2 = x2 + y2I be two NNs. Then, Smarandache [34, 35] gave their operations of NNs in the following:(1)
z1 + z2 = x1 + x2 + (y1 + y2)I.(2)
z1– z2 = x1– x2 + (y1– y2)I.(3)
z1 × z2 = x1x2 + (x1y2 + x2y1 + y1y2)I, in particular, when z1 = 0 and z2 = I, we get the equation with 0 × I = 0.(4)
z12 = (x1 + y1I)2 = x12 + (2x1y1 + y12)I, in particular, when z1 = I, we get the equation with I^2 = I.(5)
z1/z2=x1+y1I/x2+y2I=x1/x2+x2y1−x1y2/x2x2+y2I for x2 ≠ 0 and x2 ≠ –y2.(6)
z1=x1+y1I=x1−x1+x1+y1Ix1−x1−x1+y1I−x1+x1+x1+y1I−x1+x1−x1+y1I.


### 2.2. Example

There are two NNsz1 = 5 + 3I and z2 = 2 + 5I. Then, we can obtain the following results according to the above operations:(1)
z1 + z2 = x1 + x2 + (y1 + y2)I = 5 + 2 + (3 + 5)I = 7 + 8I.(2)
z1−z2 = x1−x2 + (y1−y2)I = 5−2 + (3−5)I = 3−2I.(3)
z1 × z2 = x1x2 + (x1y2 + x2y1 + y1y2)I = 5 × 2 + (5 × 5 + 3 × 2 + 3 × 5)I = 10 + 46I.(4)
z12 = (x1 + y1I)2 = x12 + (2x1y1 + y12)I = 52 + (2 ×  5 × 3 + 32)I = 25 + 39I, z22 = (x2 + y2I)2 = x22 + (2x2y2 + y22)I = 22 + (2 × 2 × 5 + 52)I = 4 + 45I.(5)
z1/z2=x1+y1I/x2+y2I=x1/x2+x2y1–x1y2/x2x2+y2I=5/2+2×3−5×5/22+5I=2.5−1.3571I.(6)
z1=x1+y1I=x1−x1+x1+y1Ix1−x1−x1+y1I−x1+x1+x1+y1I−x1+x1−x1+y1I=5−5+5+3I=2.2361−5.0645I5−5−5+3I=2.2361+0.5924I−5+5+5+3I=−2.2361+5.0645I−5+5−5+3I=−2.2361−0.5924Iz2=x2+y2I=x2−x2+x2+y2Ix2−x2−x2+y2I−x2+x2+x2+y2I−x2+x2−x2+y2I=2−2+2+5I=1.4142−4.0600I2−2−2+5I=1.4142+1.2315I−2+2+2+5I=−1.4142+4.0600I−2+2−2+5I=−1.4142−1.2315I.


## 2.1. Some Concepts and Their Operations of Neutrosophic Numbers (NNs)

The concept of NN was first proposed by Smarandache [34, 35], which consists of two parts (a determinate part and an indeterminate part). He defined the mathematical expression form z = x + yI for x, y ∈ R, where R represents all real numbers and I is indeterminacy. So, it is conveniently used in indeterminate environments.For example, consider that a NN isz = 13 + 5I. Then, its determinate part value is 13 and its indeterminate part value is 5I. When I ∈ [0, 0.5], it is equivalent to z ∈ [13, 15.5] for sure z ≥ 13.Letz1 = x1 + y1I and z2 = x2 + y2I be two NNs. Then, Smarandache [34, 35] gave their operations of NNs in the following:(1)
z1 + z2 = x1 + x2 + (y1 + y2)I.(2)
z1– z2 = x1– x2 + (y1– y2)I.(3)
z1 × z2 = x1x2 + (x1y2 + x2y1 + y1y2)I, in particular, when z1 = 0 and z2 = I, we get the equation with 0 × I = 0.(4)
z12 = (x1 + y1I)2 = x12 + (2x1y1 + y12)I, in particular, when z1 = I, we get the equation with I^2 = I.(5)
z1/z2=x1+y1I/x2+y2I=x1/x2+x2y1−x1y2/x2x2+y2I for x2 ≠ 0 and x2 ≠ –y2.(6)
z1=x1+y1I=x1−x1+x1+y1Ix1−x1−x1+y1I−x1+x1+x1+y1I−x1+x1−x1+y1I.


## 2.2. Example

There are two NNsz1 = 5 + 3I and z2 = 2 + 5I. Then, we can obtain the following results according to the above operations:(1)
z1 + z2 = x1 + x2 + (y1 + y2)I = 5 + 2 + (3 + 5)I = 7 + 8I.(2)
z1−z2 = x1−x2 + (y1−y2)I = 5−2 + (3−5)I = 3−2I.(3)
z1 × z2 = x1x2 + (x1y2 + x2y1 + y1y2)I = 5 × 2 + (5 × 5 + 3 × 2 + 3 × 5)I = 10 + 46I.(4)
z12 = (x1 + y1I)2 = x12 + (2x1y1 + y12)I = 52 + (2 ×  5 × 3 + 32)I = 25 + 39I, z22 = (x2 + y2I)2 = x22 + (2x2y2 + y22)I = 22 + (2 × 2 × 5 + 52)I = 4 + 45I.(5)
z1/z2=x1+y1I/x2+y2I=x1/x2+x2y1–x1y2/x2x2+y2I=5/2+2×3−5×5/22+5I=2.5−1.3571I.(6)
z1=x1+y1I=x1−x1+x1+y1Ix1−x1−x1+y1I−x1+x1+x1+y1I−x1+x1−x1+y1I=5−5+5+3I=2.2361−5.0645I5−5−5+3I=2.2361+0.5924I−5+5+5+3I=−2.2361+5.0645I−5+5−5+3I=−2.2361−0.5924Iz2=x2+y2I=x2−x2+x2+y2Ix2−x2−x2+y2I−x2+x2+x2+y2I−x2+x2−x2+y2I=2−2+2+5I=1.4142−4.0600I2−2−2+5I=1.4142+1.2315I−2+2+2+5I=−1.4142+4.0600I−2+2−2+5I=−1.4142−1.2315I.


## 3. Neutrosophic Number Nonlinear Programming (NN-NP)

### 3.1. NN-NP Mathematical Model

The usual mathematical model of NN-NP is represented in the following form [36, 37]:(1)MinFx,Is.t.gix,I≤ci,i=1,2,...,p,hjx,I=0,j=1,2,...,q,x∈Zn,where g1x,I,g2x,I,...,gpx,I,h1x,I,h2x,I: Zn⟶Z(Z is the set of all NNs), and I ∈ [IL, IU](the interval range of I).


### 3.2. Inventory Mathematical Model [38]

#### 3.2.1. Notations

The following notations are used in the inventory model.Three decision variables:(i)
D: demand/unit/time(ii)
Qp: production quantity/batch(iii)
Cs: setup cost/unit/timeExcept the above cost variable CS, three other cost variables are(i)
Cta: total average cost/unit/time(ii)
Ctp: total production cost/cycle(iii)
Ch: time depending on holding cost/unit/itemOther time and space variables:(i)
T: every cycle of length(ii)
Q(t): inventory level at time t (t ≥ 0)(iii)
S: total storage space area(iv)
s0: space area/unit/quantity.


#### 3.2.2. Assumptions

The inventory model is developed by considering the following assumptions:(i)
Only one item is involved in the inventory system.(ii)
The replenishment occurs with the near instantaneous response.(iii)
The startup time can be ignored.(iv)
The demand rate at any time is constant.(v)
The total production costCtp is related to the setup cost CS and production quantity QP.(vi)
Holding cost is the time depended function.


### 3.3. Model Formation

As shown in Figure1, in every time period T, the value of the production quantity Q(t) decreases from Qp to zero. The slope of the line is constant negative D and denoted by dQt/dt=−D0≤t≤T.Figure 1
Crisp inventory model.The total average cost of the cycleT (denoted by Cta) consists of three sections: setup cost (denoted by C1), holding cost (denoted by C2), and production cost (denoted by C3).(2)Cta=C1+C2+C3.Because we have the equationQ(t) = Qp−Dt, we obtain the cycle T, T=Qp/D.(3)C1=CsT=CsDQp,C2=∫0TChQtdtT=∫0TetQtdtT=eQp3/6D2TC3=CtpT=fCs−xQp−yT=fDCsxQp1+yeQp26D.Based on equations (2) and (3), we obtain the following equation:(4)Cta=CsDQp+eQp26D+fDCsxQp1+y.So, the inventory model is constructed as follows:(5)MinCtaD,Cs,Qp=CsDQp+eQp26D+fDCsxQp1+y,s.t.s0Qp≤S,D,Cs,Qp>0.


### 3.4. Solution Corresponding to Matlab Built-In Function “fmincon()”

In order to conveniently calculate the solutions, we simplify some parameters and set some constants with history records, wheree = 18, f = 5, x = 1, y = 3, s0 = 200, and S = 1100. When we assume D = x1, Cs = x2, and Qp = x3, we can obtain the following mathematical model:(6)MinCtax1,x2,x3=x1x2x3+18x326x1+5x1x2x34,s.t.200x3≤1100with maximum allowable tolerance400x1,x2,x3≥0.Assumea1 = x1−40.496I, a2 = x2 + 0.058I, and a3 = x3−2I; then, equation (6) can be expressed in the following form:(7)s.t.200x3−2I≥1100with maximum allowable tolerance400,x1−40.496I>0,x2+0.058I>0,x3−2I>0.According to the de-neutrosophication technique proposed by Ye [39] and considering I = 0 or 0.5 or 1 as the minimum or moderate or maximum indeterminacy, we can obtain three optimal solutions as follows:(1)
x1∗ = 80.615, x2∗ = 0.097, x3∗ = 7.500, and fx∗,I  = 4.187 for I = 0.(2)
x1∗ = 60.367, x2∗ = 0.126, x3∗ = 6.5, and fx∗,I  = 4.343 for I = 0.5.(3)
x1∗ = 40.119, x2∗ = 0.155, x3∗ = 5.5, and fx∗,I  = 4.525 for I = 1.Clearly, using the indeterminacyI ∈ [0, 1], different optimal results are revealed. The optimal solutions of the optimization problem are x1∗ = [40.119, 80.615], x2∗ = [0.097, 0.155], and x3∗ = [5.5, 7.5] for fx∗,I  = [4.187, 4.525], which show the interval optimal ranges.


### 3.5. Solution Corresponding to Operations of NNs

According to the front optimal solutions, we assumea1 = x1 + y1I = 80.615−40.496I, a2 = x2 + y2I = 0.097 + 0.058I, and a3 = x3−y3I = 7.500−2I, and then we give the results by equation (9):(8)MinCtaa1,a2,a3=a1a2a3+18a326a1+5a1a2a34=x1+y1Ix2+y2Ix3+y3I+18x3+y3I26x1+y1I+5x1+y1Ix2+y2Ix3+y3I4=x1x2+x1y2+x2y1+y1y2Ix3+y3I+3x32+2x3y3+y32Ix1+y1I+5x1+y1Ix2+y2Ix32+2x3y3+y32I2=x1x2+x1y2+x2y1+y1y2Ix3+y3I+3x32+2x3y3+y32Ix1+y1I+5x1+y1Ix2+y2Ix34+22x3y3+y32x32+2x3y3+y322I5x1x2x34+5x2x34y1−x1x222x3y3+y32x32+2x3y3+y322+y2x34+y222x3y3+y32x32+2x3y3+y322x2x34x2+y2x34+y32x3+y32x32+2x3y3+y32I.Becausea1 = 80.615−40.496I, a2 = 0.097 + 0.058 I, and a3 = 7.500−2I, we can get x1 = 80.615, y1 = −40.496, x2 = 0.097, y2 = 0.058, x3 = 7.5, and y3 = −2. Then, we calculate the three costs, respectively, as follows:Setup cost:(9)C1=x1+y1Ix2+y2Ix3+y3I=1.407+0.087I.Holding cost:(10)C2=18x3+y3I26x1+y1I=2.093+0.169IProduction cost:(11)C3=5x1+y1Ix2+y2Ix3+y3I4=1.047+0.082IThen, we add the three costs and obtain the total costCta with equation (2) as follows:(12)Cta=C1+C2+C3=x1+y1Ix2+y2Ix3+y3I+18x3+y3I26x1+y1I+5x1+y1Ix2+y2Ix3+y3I4=1.047+0.087I+2.093+0.169I+1.047+0.082I=4.187+0.338I.So, the calculational results validate that the same solution is obtained by using the two methods of both the Matlab built-in function “fmincon” and the operations of NNs, which arex1∗ = [40.119, 80.615], x2∗ = [0.097, 0.155], and x3∗ = [5.5, 7.5] for fx∗,I  = [4.187, 4.525]. We also obtain every cost C1 = [1.047, 1.134], C2 = [ 2.093, 2.262], and C3 = [1.047, 1.129], which are the interval optimal ranges.


## 3.1. NN-NP Mathematical Model

The usual mathematical model of NN-NP is represented in the following form [36, 37]:(1)MinFx,Is.t.gix,I≤ci,i=1,2,...,p,hjx,I=0,j=1,2,...,q,x∈Zn,where g1x,I,g2x,I,...,gpx,I,h1x,I,h2x,I: Zn⟶Z(Z is the set of all NNs), and I ∈ [IL, IU](the interval range of I).


## 3.2. Inventory Mathematical Model [38]

### 3.2.1. Notations

The following notations are used in the inventory model.Three decision variables:(i)
D: demand/unit/time(ii)
Qp: production quantity/batch(iii)
Cs: setup cost/unit/timeExcept the above cost variable CS, three other cost variables are(i)
Cta: total average cost/unit/time(ii)
Ctp: total production cost/cycle(iii)
Ch: time depending on holding cost/unit/itemOther time and space variables:(i)
T: every cycle of length(ii)
Q(t): inventory level at time t (t ≥ 0)(iii)
S: total storage space area(iv)
s0: space area/unit/quantity.


### 3.2.2. Assumptions

The inventory model is developed by considering the following assumptions:(i)
Only one item is involved in the inventory system.(ii)
The replenishment occurs with the near instantaneous response.(iii)
The startup time can be ignored.(iv)
The demand rate at any time is constant.(v)
The total production costCtp is related to the setup cost CS and production quantity QP.(vi)
Holding cost is the time depended function.


## 3.2.1. Notations

The following notations are used in the inventory model.Three decision variables:(i)
D: demand/unit/time(ii)
Qp: production quantity/batch(iii)
Cs: setup cost/unit/timeExcept the above cost variable CS, three other cost variables are(i)
Cta: total average cost/unit/time(ii)
Ctp: total production cost/cycle(iii)
Ch: time depending on holding cost/unit/itemOther time and space variables:(i)
T: every cycle of length(ii)
Q(t): inventory level at time t (t ≥ 0)(iii)
S: total storage space area(iv)
s0: space area/unit/quantity.


## 3.2.2. Assumptions

The inventory model is developed by considering the following assumptions:(i)
Only one item is involved in the inventory system.(ii)
The replenishment occurs with the near instantaneous response.(iii)
The startup time can be ignored.(iv)
The demand rate at any time is constant.(v)
The total production costCtp is related to the setup cost CS and production quantity QP.(vi)
Holding cost is the time depended function.


## 3.3. Model Formation

As shown in Figure1, in every time period T, the value of the production quantity Q(t) decreases from Qp to zero. The slope of the line is constant negative D and denoted by dQt/dt=−D0≤t≤T.Figure 1
Crisp inventory model.The total average cost of the cycleT (denoted by Cta) consists of three sections: setup cost (denoted by C1), holding cost (denoted by C2), and production cost (denoted by C3).(2)Cta=C1+C2+C3.Because we have the equationQ(t) = Qp−Dt, we obtain the cycle T, T=Qp/D.(3)C1=CsT=CsDQp,C2=∫0TChQtdtT=∫0TetQtdtT=eQp3/6D2TC3=CtpT=fCs−xQp−yT=fDCsxQp1+yeQp26D.Based on equations (2) and (3), we obtain the following equation:(4)Cta=CsDQp+eQp26D+fDCsxQp1+y.So, the inventory model is constructed as follows:(5)MinCtaD,Cs,Qp=CsDQp+eQp26D+fDCsxQp1+y,s.t.s0Qp≤S,D,Cs,Qp>0.


## 3.4. Solution Corresponding to Matlab Built-In Function “fmincon()”

In order to conveniently calculate the solutions, we simplify some parameters and set some constants with history records, wheree = 18, f = 5, x = 1, y = 3, s0 = 200, and S = 1100. When we assume D = x1, Cs = x2, and Qp = x3, we can obtain the following mathematical model:(6)MinCtax1,x2,x3=x1x2x3+18x326x1+5x1x2x34,s.t.200x3≤1100with maximum allowable tolerance400x1,x2,x3≥0.Assumea1 = x1−40.496I, a2 = x2 + 0.058I, and a3 = x3−2I; then, equation (6) can be expressed in the following form:(7)s.t.200x3−2I≥1100with maximum allowable tolerance400,x1−40.496I>0,x2+0.058I>0,x3−2I>0.According to the de-neutrosophication technique proposed by Ye [39] and considering I = 0 or 0.5 or 1 as the minimum or moderate or maximum indeterminacy, we can obtain three optimal solutions as follows:(1)
x1∗ = 80.615, x2∗ = 0.097, x3∗ = 7.500, and fx∗,I  = 4.187 for I = 0.(2)
x1∗ = 60.367, x2∗ = 0.126, x3∗ = 6.5, and fx∗,I  = 4.343 for I = 0.5.(3)
x1∗ = 40.119, x2∗ = 0.155, x3∗ = 5.5, and fx∗,I  = 4.525 for I = 1.Clearly, using the indeterminacyI ∈ [0, 1], different optimal results are revealed. The optimal solutions of the optimization problem are x1∗ = [40.119, 80.615], x2∗ = [0.097, 0.155], and x3∗ = [5.5, 7.5] for fx∗,I  = [4.187, 4.525], which show the interval optimal ranges.


## 3.5. Solution Corresponding to Operations of NNs

According to the front optimal solutions, we assumea1 = x1 + y1I = 80.615−40.496I, a2 = x2 + y2I = 0.097 + 0.058I, and a3 = x3−y3I = 7.500−2I, and then we give the results by equation (9):(8)MinCtaa1,a2,a3=a1a2a3+18a326a1+5a1a2a34=x1+y1Ix2+y2Ix3+y3I+18x3+y3I26x1+y1I+5x1+y1Ix2+y2Ix3+y3I4=x1x2+x1y2+x2y1+y1y2Ix3+y3I+3x32+2x3y3+y32Ix1+y1I+5x1+y1Ix2+y2Ix32+2x3y3+y32I2=x1x2+x1y2+x2y1+y1y2Ix3+y3I+3x32+2x3y3+y32Ix1+y1I+5x1+y1Ix2+y2Ix34+22x3y3+y32x32+2x3y3+y322I5x1x2x34+5x2x34y1−x1x222x3y3+y32x32+2x3y3+y322+y2x34+y222x3y3+y32x32+2x3y3+y322x2x34x2+y2x34+y32x3+y32x32+2x3y3+y32I.Becausea1 = 80.615−40.496I, a2 = 0.097 + 0.058 I, and a3 = 7.500−2I, we can get x1 = 80.615, y1 = −40.496, x2 = 0.097, y2 = 0.058, x3 = 7.5, and y3 = −2. Then, we calculate the three costs, respectively, as follows:Setup cost:(9)C1=x1+y1Ix2+y2Ix3+y3I=1.407+0.087I.Holding cost:(10)C2=18x3+y3I26x1+y1I=2.093+0.169IProduction cost:(11)C3=5x1+y1Ix2+y2Ix3+y3I4=1.047+0.082IThen, we add the three costs and obtain the total costCta with equation (2) as follows:(12)Cta=C1+C2+C3=x1+y1Ix2+y2Ix3+y3I+18x3+y3I26x1+y1I+5x1+y1Ix2+y2Ix3+y3I4=1.047+0.087I+2.093+0.169I+1.047+0.082I=4.187+0.338I.So, the calculational results validate that the same solution is obtained by using the two methods of both the Matlab built-in function “fmincon” and the operations of NNs, which arex1∗ = [40.119, 80.615], x2∗ = [0.097, 0.155], and x3∗ = [5.5, 7.5] for fx∗,I  = [4.187, 4.525]. We also obtain every cost C1 = [1.047, 1.134], C2 = [ 2.093, 2.262], and C3 = [1.047, 1.129], which are the interval optimal ranges.


## 4. Production Planning Mathematical Model

### 4.1. NN-LP Mathematical Model

The usual mathematical model of NN-LP is similar to mathematical model (1), so we omit it.


### 4.2. Production Planning Mathematical Model

#### 4.2.1. Notations

The following notations are used in the production planning model.Nine decision variables:(i)
a1 to a6: product quantities of six plans of type I(ii)
a7 to a8: product quantities of two plans of type II(iii)
a9: product quantities of two plans of type IIIObjective function:(i)
b: maximum profit


#### 4.2.2. Assumptions

The production planning model is developed by considering the following assumptions:(i)
Every product must pass two working procedures: A and B.(ii)
The startup time of two working procedures can be ignored.(iii)
Product quantities are only affected by validity time of machines.(iv)
The demand rate at any time is constant.


### 4.3. Model Formation

As shown in Table1, we consider an application in production planning studied by Hu [40]. A company manufactures three types of products: Types I, II, and III. All types must pass two working procedures: A and B. We consider that procedure A can be operated on machine A1 or A2 and procedure B can be operated on the machines B1, B2, and B3. Type I can be operated on all machines of procedure A and procedure B; Type II can be operated on all machines of procedure A and only machine B1 of procedure B; Type III can be operated on only machine A2 of procedure A and machine B2 of procedure B. Our aim is to schedule the optimal production planning, which can pursue for the maximum profits. All used data are listed in Table 1, including required procedure time of every working procedure, processing fees, material cost, and selling price per unit. So, Type I has six plans to produce products, along with (A1, B1) or (A1, B2) or (A1, B3) or (A2, B1) or (A2, B2) or (A2, B3), respectively. Similarly, we consider the product quantities of the six plans a1, a2, a3, a4, a5, and a6, respectively. Type II has two plans to produce products, along with (A1, B1) or (A2, B1), and Type III has one plan to produce products, along with (A2, B2). We consider the product quantities of the remaining three plans a7, a8, and a9, respectively. So, we can get the following objective function.Table 1
All data of three types of products.
MachineProductValidity time (hours/machine)Processing fees ($/hour/machine)IIIIIIA14.5 + 1.7I8 +I—5600 + 700I0.04 + 0.02IA26.7 + 1.8I8.6 + 1.4I11-I9000 + 700I0.02 + 0.01IB15.6–0.1I7.8 + 1.2I—3700 + 1100I0.05 + 0.02IB23.5 + 2.5I—10+2I6000 + 1000I0.10 + 0.02IB36.7 + 1.3I——3500 + 1500I0.04 + 0.02IMaterial cost ($/piece)0.23 + 0.03I0.30 + 0.07I0.30 + 0.05I——Selling price ($/piece)1.20 + 0.03I1.60 + 0.5I2.30 + 0.3I(13)b=1.20+0.03I−0.23+0.03I×a1+a2+a3+a4+a5+a6+1.60+0.5I−0.30+0.07I×a7+a8+2.30+0.3I−0.30+0.05I×a9−0.04+0.02I×4.5+1.7I×a1+a2+a3+8+I×a7−0.02+0.01I×0.67+0.8I×a4+a5+a6+8.6+1.4I×a8+11−I×a9−0.05+0.02I×5.6−0.1I×a1+a4+7.8+1.2I×a7+a8−0.10+0.02I×3.5+2.5I×a2+a5+10+2I×a9−0.04+0.02I×6.7+1.3I×a3+a6.So, we get the followed production planning mathematical model:(14)maxbs.t.4.5+1.71×x1+x2+x3+8+I×x7≤5600+700I,6.7+1.8I×x4+x5+x6+8.6+1.4I×x8+11−I×x9≤0.02+0.01I,5.6−0.1I×x1+x4+7.8+1.2I×x7+x8≤0.05+0.02I,3.5+2.5I×x2+x5+10+2I×x9≤0.10+0.02I,6.7+1.3I×x3+x6≤0.04+0.02I.


### 4.4. Solution regarding Matlab Built-In Function “fmincon()”

According to the de-neutrosophication technique proposed by Ye [37] and considering I = 0 or 0.5 or 1 as the minimum or moderate or maximum indeterminacy, we can obtain three optimal solutions as follows:(1)
x1∗ = 0, x2∗ = 778.508, x3∗ = 465.936, x4∗ = 0, x5∗ = 677.953, x6∗ = 56.452, x7∗ = 0, x8∗ = 474.359, x9∗ = 0, and fx∗,I  = 1297.389 for I = 0.(2)
x1∗ = 0, x2∗ = 0, x3∗ = 578.231, x4∗ = 0, x5∗ = 0, x6∗ = 0, x7∗ = 167.732, x8∗ = 338.221, x9∗ = 590.909, and fx∗,I  = 940.871 for I = 0.5.(3)
x1∗ = 0, x2∗ = 0, x3∗ = 625, x4∗ = 0, x5∗ = 0, x6∗ = 0, x7∗ = 146.667, x8∗ = 386.667, x9∗ = 583.333, and fx∗,I  = 762.717 for I = 1.Clearly, using the indeterminacyI ∈ [0, 1], different optimal results are revealed. The optimal solutions of the optimization problem are x1∗ = [0, 0], x2∗ = [0, 778.508], x3∗ = [465.936, 625], x4∗ = [0, 0], x5∗ = [0, 677.953], x6∗ = [0, 56.452], x7∗ = [0, 146.667], x8∗ = [386.667, 474.359], andx9∗ = [0, 583.333] for fx∗,I  = [762.717, 1297.389], which shows the interval optimal ranges.


### 4.5. Solution regarding Operations of NNs

According the front optimal solutions, we next calculate the nine relation formulas of the indeterminacy I and variablesa1, a2, a3, a4, a5, a6 , a7, a8, and a9. For example, let us calculate a3 = 465.936 + 159.064I. Firstly, according to three points (0, 465.936), (0.5, 578.231), and (1, 625), we obtain the linear equation (a3 = 159.06I + 476.86). Next we amend the intercept of trend curve on the vertical coordinate. The other linear equations are obtained in the same way. So, a1 = x1 + y1I = 0 + 0I = 0, a2 = x2 + y2I = 778.508−778.508I, a3 = x3 + y3I = 465.936 + 159.064I, a4 = x4 + y4I = 0 + 0I = 0, a5 = x5 + y5I = 677.953−677.953I, a6 = x6 + y6I = 56.452−56.452I, a7 = x7−y7I = 0 + 146.667I, a8 = x8 + y8I = 474.359−87.692I, and a9 = x9 + y9I = 0 + 583.333I; then, we calculate the results of equation (13) as follows:b = [(1.2 + 0.03I) − (0.23 + 0.03I)] × (a1 + a2 + a3 + a4 + a5 + a6)  +  [(1.60 + 0.5I)  −  (0.30 + 0.07I)]  ×  (a7 + a8)  + [(2.30 + 0.3I)  −  (0.30 + 0.05I)]  ×  a9 − (0.04 + 0.02I)  ×  [(4.5 + 1.7I) × (a1 + a2 + a3)  +  (8 + I) × a7]  −  (0.02 + 0.01I) × [(6.7 + 1.8I) × (a4 + a5 + a6)  +  (8.6 + 1.4I)  ×  a8  +  (11  −  I)  ×  a9]  −  (0.05 + 0.02I)  ×  [(5.6  –  0.1I) × (a1+a4) + (7.8 + 1.2I) × (a7 + a8)] − (0.10 + 0.02I)  ×  [(3.5 + 2.5I) × (a2+a5)  +  (10 + 2I) × a9] − (0.04 + 0.02I)  ×  [(6.7 + 1.3I)  ×  (a3 + a6)]  =  0.97 × (a1 + a2 +a3 + a4 + a5 + a6) + (1.30 + 0.43I)  ×  (a7 + a8) + (2.0 + 0.25I) × a9 − (0.04 + 0.02I) × [(4.5 + 1.7I) × (a1 + a2 + a3) + (8 + I) × a7] − (0.02 + 0.01I)  ×  [(6.7 + 1.8I)  ×  (a4 +  a5 +  a6)  +  (8.6  + 1.4I) × a8 + (11 – I) × a9] − (0.05 + 0.02I) × [(5.6  –  0.1I)  ×  (a1+ a4) + (7.8 + 1.2I) × (a7 + a8)] − (0.10 + 0.02I)  ×  [(3.5 + 2.5I)× (a2 + a5)+(10 + 2I) × a9]  –  (0.04  +  0.02I)  ×  [(6.7 + 1.3I)  × (a3 + a6)] = 0.97 × (a2 + a3 + a5 + a6)  +  (1.30 + 0.43I)  ×  (a7 +a8) + (2.0 + 0.25I) × a9 − (0.04  +  0.02I) × [(4.5 + 1.7I) × (a2 + a3) + (8 + I) × a7] − (0.02 + 0.01I) × [(6.7 + 1.8I) × (a5 + a6) + (8.6 + 1.4I) × a8 + (11 – I) × a9]  −  (0.05 + 0.02I) × [(7.8 + 1.2I)× (a7 + a8)]  −  (0.10 + 0.02I) × [(3.5 + 2.5I) × (a2 + a5) + (10 + 2I) × a9] − (0.04 + 0.02I) × [(6.7 + 1.3I)  ×  (a3 + a6)]  =  0.97 × (778.508  −  778.508I  +  465.936  +  159.064I  +  677.953  − 677.953I + 56.452 – 56.452I) + (1.30 + 0.43I) × (0 + 146.667I+ 474.359 − 87.692I) + (2.0 + 0.25I) × (0 + 583.333I) − (0.04+ 0.02I) × [(4.5 + 1.7I) × (778.508 − 778.508I + 465.936 + 159.064I) + (8 + I) × (0 + 146.667I)] − (0.02 + 0.01I) × [(6.7 + 1.8I) × (0 + 677.953 − 677.953I + 56.452 − 56.452I) + (8.6 + 1.4I)  ×  (474.359  −  87.692I)  +  (11  −  I)  ×  (0  +  583.333I)]  − (0.05  +  0.02I)  ×  [(7.8  +  1.2I)  ×  (0  +  146.667I  +  474.359  − 87.692I)] − (0.10 + 0.02I) × [(3.5 + 2.5I)  × (778.508 – 778.508I + 677.953 − 677.953I) + (10 + 2I) × (0 + 583.333I)]−  (0.04  +  0.02I)  ×  [(6.7  +  1.3I)  ×  (465.936  +  159.064I + 56.452 − 56.452I)]  =  0.97  ×  (1978.849  −  1353.849I) + (1.3 + 0.43I)  ×  (474.359  +  58.975I) + (2 + 0.25I) × (0 + 583.333I) − (0.04 + 0.02I) × (5599.998 − 404.995I) − (0.02 + 0.01I) × (9000.001  +  699.9991I)  − (0.05 + 0.02I) × (3700 + 1100.006I)− (0.10 + 0.02I) × (5097.614  +  1902.383I)  −  (0.04  +  0.02I) × (3500 + 1500I) = 1919.484  −  1313.234I + 616.6667 + 306.001I + 1312.499I − 224 − 87.7I  −  180  −  111.000I  −  185 − 509.761 − 330.238I − 140 − 160.000I = 1297.389 − 534.672I.So, these calculational results validate that the same solution is obtained by using the two methods of both the Matlab built-in function “fmincon()” and the operations of NNs, which arex1∗ = [0, 0], x2∗ = [0, 778.508], x3∗ = [465.936, 625], x4∗ = [0, 0], x5∗ = [0, 677.953], x6∗ = [0, 56.452], x7∗ = [0, 146.667], x8∗ = [386.667, 474.359], andx9∗ = [0, 583.333] for fx∗,I  = [762.717, 1297.389] and show the interval optimal ranges.


## 4.1. NN-LP Mathematical Model

The usual mathematical model of NN-LP is similar to mathematical model (1), so we omit it.


## 4.2. Production Planning Mathematical Model

### 4.2.1. Notations

The following notations are used in the production planning model.Nine decision variables:(i)
a1 to a6: product quantities of six plans of type I(ii)
a7 to a8: product quantities of two plans of type II(iii)
a9: product quantities of two plans of type IIIObjective function:(i)
b: maximum profit


### 4.2.2. Assumptions

The production planning model is developed by considering the following assumptions:(i)
Every product must pass two working procedures: A and B.(ii)
The startup time of two working procedures can be ignored.(iii)
Product quantities are only affected by validity time of machines.(iv)
The demand rate at any time is constant.


## 4.2.1. Notations

The following notations are used in the production planning model.Nine decision variables:(i)
a1 to a6: product quantities of six plans of type I(ii)
a7 to a8: product quantities of two plans of type II(iii)
a9: product quantities of two plans of type IIIObjective function:(i)
b: maximum profit


## 4.2.2. Assumptions

The production planning model is developed by considering the following assumptions:(i)
Every product must pass two working procedures: A and B.(ii)
The startup time of two working procedures can be ignored.(iii)
Product quantities are only affected by validity time of machines.(iv)
The demand rate at any time is constant.


## 4.3. Model Formation

As shown in Table1, we consider an application in production planning studied by Hu [40]. A company manufactures three types of products: Types I, II, and III. All types must pass two working procedures: A and B. We consider that procedure A can be operated on machine A1 or A2 and procedure B can be operated on the machines B1, B2, and B3. Type I can be operated on all machines of procedure A and procedure B; Type II can be operated on all machines of procedure A and only machine B1 of procedure B; Type III can be operated on only machine A2 of procedure A and machine B2 of procedure B. Our aim is to schedule the optimal production planning, which can pursue for the maximum profits. All used data are listed in Table 1, including required procedure time of every working procedure, processing fees, material cost, and selling price per unit. So, Type I has six plans to produce products, along with (A1, B1) or (A1, B2) or (A1, B3) or (A2, B1) or (A2, B2) or (A2, B3), respectively. Similarly, we consider the product quantities of the six plans a1, a2, a3, a4, a5, and a6, respectively. Type II has two plans to produce products, along with (A1, B1) or (A2, B1), and Type III has one plan to produce products, along with (A2, B2). We consider the product quantities of the remaining three plans a7, a8, and a9, respectively. So, we can get the following objective function.Table 1
All data of three types of products.
MachineProductValidity time (hours/machine)Processing fees ($/hour/machine)IIIIIIA14.5 + 1.7I8 +I—5600 + 700I0.04 + 0.02IA26.7 + 1.8I8.6 + 1.4I11-I9000 + 700I0.02 + 0.01IB15.6–0.1I7.8 + 1.2I—3700 + 1100I0.05 + 0.02IB23.5 + 2.5I—10+2I6000 + 1000I0.10 + 0.02IB36.7 + 1.3I——3500 + 1500I0.04 + 0.02IMaterial cost ($/piece)0.23 + 0.03I0.30 + 0.07I0.30 + 0.05I——Selling price ($/piece)1.20 + 0.03I1.60 + 0.5I2.30 + 0.3I(13)b=1.20+0.03I−0.23+0.03I×a1+a2+a3+a4+a5+a6+1.60+0.5I−0.30+0.07I×a7+a8+2.30+0.3I−0.30+0.05I×a9−0.04+0.02I×4.5+1.7I×a1+a2+a3+8+I×a7−0.02+0.01I×0.67+0.8I×a4+a5+a6+8.6+1.4I×a8+11−I×a9−0.05+0.02I×5.6−0.1I×a1+a4+7.8+1.2I×a7+a8−0.10+0.02I×3.5+2.5I×a2+a5+10+2I×a9−0.04+0.02I×6.7+1.3I×a3+a6.So, we get the followed production planning mathematical model:(14)maxbs.t.4.5+1.71×x1+x2+x3+8+I×x7≤5600+700I,6.7+1.8I×x4+x5+x6+8.6+1.4I×x8+11−I×x9≤0.02+0.01I,5.6−0.1I×x1+x4+7.8+1.2I×x7+x8≤0.05+0.02I,3.5+2.5I×x2+x5+10+2I×x9≤0.10+0.02I,6.7+1.3I×x3+x6≤0.04+0.02I.


## 4.4. Solution regarding Matlab Built-In Function “fmincon()”

According to the de-neutrosophication technique proposed by Ye [37] and considering I = 0 or 0.5 or 1 as the minimum or moderate or maximum indeterminacy, we can obtain three optimal solutions as follows:(1)
x1∗ = 0, x2∗ = 778.508, x3∗ = 465.936, x4∗ = 0, x5∗ = 677.953, x6∗ = 56.452, x7∗ = 0, x8∗ = 474.359, x9∗ = 0, and fx∗,I  = 1297.389 for I = 0.(2)
x1∗ = 0, x2∗ = 0, x3∗ = 578.231, x4∗ = 0, x5∗ = 0, x6∗ = 0, x7∗ = 167.732, x8∗ = 338.221, x9∗ = 590.909, and fx∗,I  = 940.871 for I = 0.5.(3)
x1∗ = 0, x2∗ = 0, x3∗ = 625, x4∗ = 0, x5∗ = 0, x6∗ = 0, x7∗ = 146.667, x8∗ = 386.667, x9∗ = 583.333, and fx∗,I  = 762.717 for I = 1.Clearly, using the indeterminacyI ∈ [0, 1], different optimal results are revealed. The optimal solutions of the optimization problem are x1∗ = [0, 0], x2∗ = [0, 778.508], x3∗ = [465.936, 625], x4∗ = [0, 0], x5∗ = [0, 677.953], x6∗ = [0, 56.452], x7∗ = [0, 146.667], x8∗ = [386.667, 474.359], andx9∗ = [0, 583.333] for fx∗,I  = [762.717, 1297.389], which shows the interval optimal ranges.


## 4.5. Solution regarding Operations of NNs

According the front optimal solutions, we next calculate the nine relation formulas of the indeterminacy I and variablesa1, a2, a3, a4, a5, a6 , a7, a8, and a9. For example, let us calculate a3 = 465.936 + 159.064I. Firstly, according to three points (0, 465.936), (0.5, 578.231), and (1, 625), we obtain the linear equation (a3 = 159.06I + 476.86). Next we amend the intercept of trend curve on the vertical coordinate. The other linear equations are obtained in the same way. So, a1 = x1 + y1I = 0 + 0I = 0, a2 = x2 + y2I = 778.508−778.508I, a3 = x3 + y3I = 465.936 + 159.064I, a4 = x4 + y4I = 0 + 0I = 0, a5 = x5 + y5I = 677.953−677.953I, a6 = x6 + y6I = 56.452−56.452I, a7 = x7−y7I = 0 + 146.667I, a8 = x8 + y8I = 474.359−87.692I, and a9 = x9 + y9I = 0 + 583.333I; then, we calculate the results of equation (13) as follows:b = [(1.2 + 0.03I) − (0.23 + 0.03I)] × (a1 + a2 + a3 + a4 + a5 + a6)  +  [(1.60 + 0.5I)  −  (0.30 + 0.07I)]  ×  (a7 + a8)  + [(2.30 + 0.3I)  −  (0.30 + 0.05I)]  ×  a9 − (0.04 + 0.02I)  ×  [(4.5 + 1.7I) × (a1 + a2 + a3)  +  (8 + I) × a7]  −  (0.02 + 0.01I) × [(6.7 + 1.8I) × (a4 + a5 + a6)  +  (8.6 + 1.4I)  ×  a8  +  (11  −  I)  ×  a9]  −  (0.05 + 0.02I)  ×  [(5.6  –  0.1I) × (a1+a4) + (7.8 + 1.2I) × (a7 + a8)] − (0.10 + 0.02I)  ×  [(3.5 + 2.5I) × (a2+a5)  +  (10 + 2I) × a9] − (0.04 + 0.02I)  ×  [(6.7 + 1.3I)  ×  (a3 + a6)]  =  0.97 × (a1 + a2 +a3 + a4 + a5 + a6) + (1.30 + 0.43I)  ×  (a7 + a8) + (2.0 + 0.25I) × a9 − (0.04 + 0.02I) × [(4.5 + 1.7I) × (a1 + a2 + a3) + (8 + I) × a7] − (0.02 + 0.01I)  ×  [(6.7 + 1.8I)  ×  (a4 +  a5 +  a6)  +  (8.6  + 1.4I) × a8 + (11 – I) × a9] − (0.05 + 0.02I) × [(5.6  –  0.1I)  ×  (a1+ a4) + (7.8 + 1.2I) × (a7 + a8)] − (0.10 + 0.02I)  ×  [(3.5 + 2.5I)× (a2 + a5)+(10 + 2I) × a9]  –  (0.04  +  0.02I)  ×  [(6.7 + 1.3I)  × (a3 + a6)] = 0.97 × (a2 + a3 + a5 + a6)  +  (1.30 + 0.43I)  ×  (a7 +a8) + (2.0 + 0.25I) × a9 − (0.04  +  0.02I) × [(4.5 + 1.7I) × (a2 + a3) + (8 + I) × a7] − (0.02 + 0.01I) × [(6.7 + 1.8I) × (a5 + a6) + (8.6 + 1.4I) × a8 + (11 – I) × a9]  −  (0.05 + 0.02I) × [(7.8 + 1.2I)× (a7 + a8)]  −  (0.10 + 0.02I) × [(3.5 + 2.5I) × (a2 + a5) + (10 + 2I) × a9] − (0.04 + 0.02I) × [(6.7 + 1.3I)  ×  (a3 + a6)]  =  0.97 × (778.508  −  778.508I  +  465.936  +  159.064I  +  677.953  − 677.953I + 56.452 – 56.452I) + (1.30 + 0.43I) × (0 + 146.667I+ 474.359 − 87.692I) + (2.0 + 0.25I) × (0 + 583.333I) − (0.04+ 0.02I) × [(4.5 + 1.7I) × (778.508 − 778.508I + 465.936 + 159.064I) + (8 + I) × (0 + 146.667I)] − (0.02 + 0.01I) × [(6.7 + 1.8I) × (0 + 677.953 − 677.953I + 56.452 − 56.452I) + (8.6 + 1.4I)  ×  (474.359  −  87.692I)  +  (11  −  I)  ×  (0  +  583.333I)]  − (0.05  +  0.02I)  ×  [(7.8  +  1.2I)  ×  (0  +  146.667I  +  474.359  − 87.692I)] − (0.10 + 0.02I) × [(3.5 + 2.5I)  × (778.508 – 778.508I + 677.953 − 677.953I) + (10 + 2I) × (0 + 583.333I)]−  (0.04  +  0.02I)  ×  [(6.7  +  1.3I)  ×  (465.936  +  159.064I + 56.452 − 56.452I)]  =  0.97  ×  (1978.849  −  1353.849I) + (1.3 + 0.43I)  ×  (474.359  +  58.975I) + (2 + 0.25I) × (0 + 583.333I) − (0.04 + 0.02I) × (5599.998 − 404.995I) − (0.02 + 0.01I) × (9000.001  +  699.9991I)  − (0.05 + 0.02I) × (3700 + 1100.006I)− (0.10 + 0.02I) × (5097.614  +  1902.383I)  −  (0.04  +  0.02I) × (3500 + 1500I) = 1919.484  −  1313.234I + 616.6667 + 306.001I + 1312.499I − 224 − 87.7I  −  180  −  111.000I  −  185 − 509.761 − 330.238I − 140 − 160.000I = 1297.389 − 534.672I.So, these calculational results validate that the same solution is obtained by using the two methods of both the Matlab built-in function “fmincon()” and the operations of NNs, which arex1∗ = [0, 0], x2∗ = [0, 778.508], x3∗ = [465.936, 625], x4∗ = [0, 0], x5∗ = [0, 677.953], x6∗ = [0, 56.452], x7∗ = [0, 146.667], x8∗ = [386.667, 474.359], andx9∗ = [0, 583.333] for fx∗,I  = [762.717, 1297.389] and show the interval optimal ranges.


## 5. Conclusion

This paper first introduced some concepts and their operations of NNs with indeterminacyI. Next, we built a mathematical model with constrained conditions and then constructed the corresponding inventory model and production planning model. Finally, we obtained the optimal solutions by using the two methods of the Matlab built-in function “fmincon()” and the operations of NNs to solve the NN-NP and NN-LP problems with constrained conditions as preliminary application study in indeterminate setting. The final results show that the two methods obtained the same effective solutions, but the former needs the Matlab built-in function along with the simple calculational process, while the latter needs a lot of operations of NNs along with the complex calculational process. Some contributions in this study are that (1) different methods can obtain the same optimal results, (2) the NN-NP and NN-LP methods provided the new application ways for engineering management, (3) the NN-NP and NN-LP methods are more suitable than other ones under uncertain environments as the generalization of traditional programming methods, and (4) the two approaches can obtain the interval solutions for avoiding determinate solutions of traditional programming methods.Obviously, the proposed NN-LP and NN-NP methods can handle indeterminate and/or determinate mathematical programming problems, which are the generalization of existing uncertain or certain linear and nonlinear programming methods. As the preliminary application study in this paper, however, there exist a lot of mathematical solution methods and proof problems along with some complexity/difficulty in the nonlinear programming problems which need to be studied further. Hence, as our future works, one is to further analyze the two presented methods of this paper from the mathematical problems, such as the convexity problem in the nonlinear programming, the stability and solution range problems regarding the changeability of NNs, and the sensitivities of NNs on the solution results, and then NN-LP and NN-NP approaches will be extended to other ﬁelds, such as engineering design and management science.



---

*Source: 6668711-2021-04-19.xml*