# A Genetic-Algorithms-Based Approach for Programming Linear and Quadratic Optimization Problems with Uncertainty

Authors: Weihua Jin; Zhiying Hu; Christine W. Chan
Journal: Mathematical Problems in Engineering
 (2013)

---

## Abstract

This paper proposes a genetic-algorithms-based approach as an all-purpose problem-solving method for operation programming problems under uncertainty. The proposed method was applied for management of a municipal solid waste treatment system. Compared to the traditional interactive binary analysis, this approach has fewer limitations and is able to reduce the complexity in solving the inexact linear programming problems and inexact quadratic programming problems. The implementation of this approach was performed using the Genetic Algorithm Solver of MATLAB (trademark of MathWorks). The paper explains the genetic-algorithms-based method and presents details on the computation procedures for each type of inexact operation programming problems. A comparison of the results generated by the proposed method based on genetic algorithms with those produced by the traditional interactive binary analysis method is also presented.


---

## Body

## 1. Introduction

Economic optimization in operation planning of municipal solid waste management was first proposed in the 1960s [1]. Since then, different models of waste management planning have been proposed, which include linear programming [2–5], mixed integer linear programming, dynamic programming, multiobjective programming [6–8], and hybrids of these methods combined with probability, fuzzy, and inexact analyses [9–12]. The objectives of these waste management models include reduction of total cost, protection of the environment, and reuse of waste material and energy.In real-life engineering problems, the available information often cannot be represented as deterministic numbers or distribution functions. Instead, it is often possible to represent the available information as inexact numbers, which can be readily used in the inexact programming models. This may be due to the fact that decision makers often prefer  to have an inexact representation of uncertainty than provide a specification for distributions of fuzzy sets [13–17]. For operational planning, the inexact analysis approach typically treats the uncertain parameters as intervals with known upper and lower bounds but unclear distributions. A major advantage of inexact analysis in operation planning is that variation of system performance and decision variables can be investigated by solving relatively simple submodels.Research work on different kinds of inexact programming, such as inexact linear programming (ILP), inexact quadratic programming (IQP), inexact integer programming (IIP), inexact dynamic programming (IDP), and inexact multiobjective programming (IMOP) [6, 14, 15, 18–26], has been conducted. This paper presents an alternative heuristic-based method, which involves generic linear and quadratic programming with inexact information; the approach adopted involves the use of genetic algorithms (GA).This paper is organized as follows. Section2 presents the background of this research, which includes an introduction of the Genetic Algorithm Solver of MATLAB used for implementing the proposed method and the concepts of ILP, IQP, and their interactive binary analysis solution method [18]. Section 3 discusses the methodology of the proposed genetic-algorithms-based methods for solving inexact liner problems and inexact quadratic problems. Section 4 presents the solution of the IQP problem of solid waste disposal planning as a case study.


## 2. Background

Linear programming and nonlinear programming are considered powerful optimization tools suitable for modeling and solving complex optimization problems in engineering. To handle uncertainty in real world data, inexact parameters and constraints are combined with various kinds of optimization techniques. Huang et al. [18–21] proposed two inexact nonlinear programming methods by introducing internal and fuzzy numbers into the quadratic programming (QP) frameworks. The methods of inexact quadratic programming (IQP) and inexact-fuzzy quadratic programming are applicable for operation planning of solid waste management systems. Often a detailed solution of IQP involves a large number of direct comparisons to interactively identify the uncertain relationships among the objective function and decision variables, whether the problems are medium-sized or larger-scaled. When these methods are applied to complicated and nonlinear problems, the number of direct comparisons can become exponential. In such a situation, we suggest that GA are a feasible problem-solving method.GA have been applied as the optimization techniques for solving complex and nonlinear problems in operations research, industrial engineering, and management science. The GA method is a suitable optimization tool especially for solving problems, which involve nonsmooth and multi-modal search spaces. An engineering problem that has traditionally been solved as an IQP problem often involves a large and uneven search space, for which a global optimal solution is often not required. Instead, we suggest that the GA-based method is a more effective problem-solving approach than some traditional inexact programming methods.


### 2.1. Genetic Algorithm Solver for MATLAB

For implementation of genetic algorithms, the Genetic Algorithm Solver of Global Optimization Toolbox (GASGOT), developed by MATLAB (Trademark of MathWorks), has been adopted. GASGOT implements simulated evolution in the MATLAB environment using both binary and floating point representations and ordered base representation. This enables flexible implementation of the genetic operators, selection functions, termination functions, and evaluation functions. GASGOT was developed by the Department of Industrial Engineering of North Carolina State University as a toolbox of MATLAB. Hence, it runs in a MATLAB workspace and can be easily invoked by other programs. GASGOT supports implementation of binary chromosomes, binary mutation, and simple crossover. For floating point representation, the operators of uniform mutation, nonuniform mutation, multinonuniform mutation, boundary mutation, simple crossover, arithmetic crossover, and heuristic crossover are defined.The GASGOT is adopted as the problem-solving engine of both the GA linear program and GA nonlinear program; all the applications and numeric examples were calculated using this solver in MATLAB.


### 2.2. Inexact Linear Programming and Its Problem-Solving Approach

To support decision making involving uncertainties, Huang et al. [20, 22] proposed an interactive binary analysis to solve the inexact linear programming (ILP) problem.A typical ILP problem can be expressed as follows:(1)Maxf±=∑j=1n[cj±xj±],s.t.∑j=1naij±xj±≤bi±,i=1,2,…m,xj±≥0,j=1,2,…,n,
where aij±, bi±, and cj± are inexact parameters and xj± is an inexact variable. It is assumed that an optimal solution exists. For an inexact number g±∈[g-,g+], g+ and g- are the upper and lower bounds, respectively.The traditional binary solution procedure is specified as follows.Step 1.
Group symbols for inexact coefficientscj±; let former k1 coefficients be positive and latter k2 be negative; k1+k2=n.Step 2.
Define the upper and lower bounds of the objective function asf+ andf-.Step 3.
Define absolute values and signs for the coefficients of the constraintsaij±.Step 4.
Define the relationships between the decision variablesxj± and the absolute value of the coefficients of the constraints |aij±|.Step 5.
Formulate constraints corresponding to the upper and lower bounds of the objective functionf+ and f-.Step 6.
When the right-hand side of the constraintsbi are also inexact numbers, define the relationships between f± and bi±.Step 7.
Specify the two submodels.For a detailed description of the procedure, see Huang et al. [20, 22, 23].


### 2.3. Inexact Quadratic Programming and Its Problem-Solving Approach

A typical IQP problem is formulated as follows:(2)Maxf±=∑j=1n[cj±xj±+dj±(xj±)2],s.t.∑j=1naij±xj±≤bi±,i=1,2,…m,xj±≥0,j=1,2,…,n,
where aij±, bij±, cij±, and dij± are inexact parameters, xj± is an inexact variable, and it is assumed that an optimal solution exists.The solution procedure is similar to that of ILP but involves more complexity and computation; the main steps of the solution procedure are listed as follows; for a detailed description, see Huang et al. [18, 24, 25].Step 1.
Group symbols for inexact coefficientscj± and dij±; when cj± and dij± have the same signs, similar to the ILP, let former k1 coefficients be positive and latter k2 be negative; k1+k2=n. When cj± and dij± have different signs, 2n combinations of the upper and lower bounds of xj± have to be formulated for the objective function, which will require a large number of computations.Step 2.
Define the upper and lower bounds of the objective function asf+ and f-.Step 3.
Define the absolute values and signs for the coefficients of the constraintsaij± and the relationships between the decision variables xj± and the absolute value of the coefficients of the constraints |aij±|.
When somexj- corresponds to f+ and some xj+ corresponds to f-, the specification of the constraints aij±xj±≤bi± requires a comparison of the contribution of xj+ and xj- groups to the sum ∑j=1naij±xj±, when f+ is desired. When the problem is complex, a direct comparison of the dominance of xj+ and xj- becomes impossible; then some simplification and assumption need to be considered, which will affect the quality of the result.Step 4.
Formulate the constraints corresponding to the upper and lower bounds of the objective functionf+ and f-.Step 5.
When the right-hand side of the constraintsbi are also inexact numbers, define the relationships between f± and bi±.Step 6.
Specify the two submodels.The genetic-algorithms-based methods to solve the above inexact linear problem and inexact quadratic problem will be presented in the next section, and the results from the GA-based methods will be compared to those generated using the traditional approach in [18, 20, 22–25].


## 2.1. Genetic Algorithm Solver for MATLAB

For implementation of genetic algorithms, the Genetic Algorithm Solver of Global Optimization Toolbox (GASGOT), developed by MATLAB (Trademark of MathWorks), has been adopted. GASGOT implements simulated evolution in the MATLAB environment using both binary and floating point representations and ordered base representation. This enables flexible implementation of the genetic operators, selection functions, termination functions, and evaluation functions. GASGOT was developed by the Department of Industrial Engineering of North Carolina State University as a toolbox of MATLAB. Hence, it runs in a MATLAB workspace and can be easily invoked by other programs. GASGOT supports implementation of binary chromosomes, binary mutation, and simple crossover. For floating point representation, the operators of uniform mutation, nonuniform mutation, multinonuniform mutation, boundary mutation, simple crossover, arithmetic crossover, and heuristic crossover are defined.The GASGOT is adopted as the problem-solving engine of both the GA linear program and GA nonlinear program; all the applications and numeric examples were calculated using this solver in MATLAB.


## 2.2. Inexact Linear Programming and Its Problem-Solving Approach

To support decision making involving uncertainties, Huang et al. [20, 22] proposed an interactive binary analysis to solve the inexact linear programming (ILP) problem.A typical ILP problem can be expressed as follows:(1)Maxf±=∑j=1n[cj±xj±],s.t.∑j=1naij±xj±≤bi±,i=1,2,…m,xj±≥0,j=1,2,…,n,
where aij±, bi±, and cj± are inexact parameters and xj± is an inexact variable. It is assumed that an optimal solution exists. For an inexact number g±∈[g-,g+], g+ and g- are the upper and lower bounds, respectively.The traditional binary solution procedure is specified as follows.Step 1.
Group symbols for inexact coefficientscj±; let former k1 coefficients be positive and latter k2 be negative; k1+k2=n.Step 2.
Define the upper and lower bounds of the objective function asf+ andf-.Step 3.
Define absolute values and signs for the coefficients of the constraintsaij±.Step 4.
Define the relationships between the decision variablesxj± and the absolute value of the coefficients of the constraints |aij±|.Step 5.
Formulate constraints corresponding to the upper and lower bounds of the objective functionf+ and f-.Step 6.
When the right-hand side of the constraintsbi are also inexact numbers, define the relationships between f± and bi±.Step 7.
Specify the two submodels.For a detailed description of the procedure, see Huang et al. [20, 22, 23].


## 2.3. Inexact Quadratic Programming and Its Problem-Solving Approach

A typical IQP problem is formulated as follows:(2)Maxf±=∑j=1n[cj±xj±+dj±(xj±)2],s.t.∑j=1naij±xj±≤bi±,i=1,2,…m,xj±≥0,j=1,2,…,n,
where aij±, bij±, cij±, and dij± are inexact parameters, xj± is an inexact variable, and it is assumed that an optimal solution exists.The solution procedure is similar to that of ILP but involves more complexity and computation; the main steps of the solution procedure are listed as follows; for a detailed description, see Huang et al. [18, 24, 25].Step 1.
Group symbols for inexact coefficientscj± and dij±; when cj± and dij± have the same signs, similar to the ILP, let former k1 coefficients be positive and latter k2 be negative; k1+k2=n. When cj± and dij± have different signs, 2n combinations of the upper and lower bounds of xj± have to be formulated for the objective function, which will require a large number of computations.Step 2.
Define the upper and lower bounds of the objective function asf+ and f-.Step 3.
Define the absolute values and signs for the coefficients of the constraintsaij± and the relationships between the decision variables xj± and the absolute value of the coefficients of the constraints |aij±|.
When somexj- corresponds to f+ and some xj+ corresponds to f-, the specification of the constraints aij±xj±≤bi± requires a comparison of the contribution of xj+ and xj- groups to the sum ∑j=1naij±xj±, when f+ is desired. When the problem is complex, a direct comparison of the dominance of xj+ and xj- becomes impossible; then some simplification and assumption need to be considered, which will affect the quality of the result.Step 4.
Formulate the constraints corresponding to the upper and lower bounds of the objective functionf+ and f-.Step 5.
When the right-hand side of the constraintsbi are also inexact numbers, define the relationships between f± and bi±.Step 6.
Specify the two submodels.The genetic-algorithms-based methods to solve the above inexact linear problem and inexact quadratic problem will be presented in the next section, and the results from the GA-based methods will be compared to those generated using the traditional approach in [18, 20, 22–25].


## 3. Methodology

### 3.1. Genetic-Algorithms-Based Method for Solution of ILP Problems (GAILP)

A GA, as a heuristic search algorithm, has been adopted for solving the aforementioned ILP problem. In the GA approach, the upper and lower bounds of the inexact numbers of coefficientsaij±, bi±, and cj± can be determined by substituting the initial suboptimal decision variables into the objective function. f+and f-can be calculated directly without any uncertainty in the coefficients. This approach is called the genetic-algorithms-based method for solving ILP problems or the GAILP method.GAILP has been designed to include three stages, which are discussed as follows.The objective of the first stage is to get an initial suboptimalxjs for the following problem, which is transformed from the ILP problem defined in (1):
(3)Maxf±=∑j=1n[cjrxjs],s.t.∑j=1naijrxjs≤bir,i=1,2,…,m,xj≥0,j=1,2,…,n,
where aijr, bir, and cjr are random numbers that satisfy the continuous uniform distribution in the intervals of [aij-,aij+], [bi-,bi+], and [cj-,cj+], respectively. Then, the problem is solved by the GA linear program solving engine of GASGOT, which uses the objective function in (3) as the positive term of the fitness function and the constraints of (1) as the negative punishment terms. Thus, a suboptimal solution fs can be identified and the corresponding decision variables of xjs are also obtained.In the second stage, the inexact coefficients ofaij±,  bi±, and  cj± will be determined. Let the determined coefficients corresponding tof+ be aij±+, bi±+, and cj±+ and those corresponding to f- be aij±-,  bi±-, and  cj±-. These two sets of coefficients can be obtained using the following method.Substitutingxjs into (1) will convert (1) into the following equation:
(4)Maxf±=∑j=1n[cj±xjs],s.t.∑j=1naij±xjs≤bi±,i=1,2,…,m,xjs≥0,j=1,2,…,n.
To identify the coefficients aij±,  bi±, and cj± corresponding to f±, a set of objective functions needs to be constructed and solved. Since xjs are suboptimal variables, which tend to make the objective function closer to f+, consider aij±,  bi±, and  cj± as variables; then the objective function of (5) can be constructed so as to find cj±+:
(5)Maxf±=∑j=1n[cj±xjs],s.t.∑j=1naij±xjs≤bi±,i=1,2,…,m.
The coefficients cj±+ are considered as corresponding to f+.At the same time, the objective function presented in (6) can be constructed so as to find cj±-:
(6)Minf±=∑j=1n[cj±xjs],s.t.∑j=1naij±xjs≤bi±,i=1,2,…,m.
There are two kinds of decision schemes for inexact programming problems, which are the conservative schemes and optimistic schemes [26]. The former assumes less risk than the latter, so that, for a maximization objective function, planning for the lower bound of an objective value represents the conservative scheme, and planning for the upper bound of an objective value represents the optimistic scheme [26]. In terms of constraints, the conservative scheme involves more rigorous or stringent constraints, and the optimistic scheme adopts more tolerant ones.Thus, the problem of searching foraij±+ and  bi±+ of the optimistic scheme and corresponding to the upper bound of the objective value of f+ can be represented as follows:
(7)Max∑j=1nabs(aij±xjs-bi±),s.t.∑j=1naij±xjs≤bi±,i=1,2,…,m.
The problem
(8)Min∑j=1nabs(aij±xjs-bi±),s.t.∑j=1naij±xjs≤bi±,i=1,2,…,m,
will give aij±- and bi±- of the conservative scheme, corresponding to the lower bound of the objective value of f-.Hence, the values ofaij±+, bi±+, and cj±+and aij±-, bi±-, and cj±- can be calculated.In the third stage, the problem represented in (1) is converted into the following two subproblems.Forf+,
(9)Maxf+=∑j=1n[cj±+xj±],s.t.∑j=1naij±+xj±≤bi±+,i=1,2,…,m,xj±≥0,j=1,2,…,n.Forf-,
(10)Maxf-=∑j=1n[cj±-xj±],s.t.∑j=1naij±-xj±≤bi±-,i=1,2,…,m,xj±≥0,j=1,2,…,n.
This step eliminates the inexact parameters in (1) and generates instead (9) and (10) as typical linear programming (LP) problems, which can be solved easily using the traditional methods. Generally speaking, the interactive binary algorithm of Huang et al.  [20, 22] can be used for solving inexact linear problems reliably and relatively quickly for many real-life decision-making scenarios in the engineering field. However, this binary algorithm has some limitations. One of them, for example, is the limitation that the upper and lower bounds of an inexact coefficient cannot have different signs. By contrast, the GAILP does not have this kind of limitation because the GA method does not depend on any assumed distribution of the inexact parameter.The following inexact linear problem demonstrates how GAILP is able to handle this situation:(11)Maxf±=[2,6]x1±+[7.8,10.2]x2±,s.t.x1±+[-0.8,1.2]x2±≤4.6,[-6,-8.5]x1±+12x2±≤[-2,3],x1±≥0,x2±≥0.
In stage one, the suboptimal fs and the corresponding x1s and x2s are calculated:
(12)fs=112,x1s=10.1,x2s=5.8.
In stage two, substituting x1s and x2s into (11), the two sets of problems for determining aij±+,  bi±+, and  cj±+ and aij±-,  bi±-, and  cj±- are constructed forf+ and f-, respectively.Forf+, cj±+can be determined by solving the following problem:
(13)Maxf+=c1±+10.1+c2±+5.8,s.t.10.1+[-0.8,1.2]5.8≤4.6,[-6,-8.5]10.1+125.8≤[-2,3],c1±+≤[2,6],c2±+≤[7.8,10.2].
And aij±+ and bi±+ can be determined by solving the following problem:
(14)Max4.6-(x1s+a12x2s),Maxb2-(a21x1s+12x2s),a12=[-0.8,1.2],a21=[-6,-8.5],b2=[-2,3].
For f-, cj±-can be determined by solving the following problem:
(15)Minf-=c1±-10.1+c2±-5.8,s.t.10.1+[-0.8,1.2]5.8≤4.6,[-6,-8.5]10.1+125.8≤[-2,3],c1±-≤[2,6],c2±-≤[7.8,10.2].
And aij±- and bi±- can be determined by solving the following problem:
(16)Min4.6-(x1s+a12x2s),Minb2-(a21x1s+12x2s),a12=[-0.8,1.2],a21=[-8.5,-6],b2=[-2,3].
By solving (13), (14), (15), and (16), the coefficients are calculated as follows:
(17)a12±+=-0.8,a21±+=-8.5,b2±+=3,c1±+=6,c2±+=10.2,a12±-=1.2,a21±-=-6,b2±-=-2,c1±-=2,c2±-=7.8.
In stage three, the two submodels for the optimistic scheme and the conservative scheme are constructed, respectively, as follows.For the optimistic scheme,f+,
(18)Maxf+=6x1++10.2x2+,s.t.x1+-0.8x2+≤4.6,-8.5x1++12x2+≤3,x1+≥0,x2+≥0.
For the conservative scheme, f-,
(19)Maxf-=2x1-+7.8x2-,s.t.x1-+1.2x2-≤4.6,-6x1-+12x2-≤2,x1-≥0,x2-≥0.
The final result can be found by solving the previous two submodels:
(20)f+=149.04,x1+=11.08,x2+=8.1,f-=16.4,x1+=3,x2+=1.33.
From the previous calculation, it can be seen that the GAILP method can be used without any assumption of the upper and lower bounds of the inexact coefficients. In fact, this method effectively extends the scope of problems solvable using the methods of the inexact linear programming problem. Therefore, the GAILP method is more adaptable for real world applications of optimization problems with uncertainty. In the next section, this method will be extended to solve the inexact quadratic problems.A sample inexact linear programming problem in [22] is as follows:
(21)Maxf±=c1x1±+c2x2±,s.t.a11x1±+a12x2±≤b1,a21x1±+a22x2±≤b2,
where c1=[26,30],c2=[-6,-5.5],
(22)a11=[8,10],a12=[-14,-12],b1=[3.8,4.2],a21=[2.4,2.8],a22=[3.4,4],b2=6.5.
By using the traditional interactive binary algorithm, two submodels are obtained:
(23)Maxf+=30x1+-5.5x2-,s.t.8x1+-14x2-≤4.2,2.4x1++4x2-≤6.5,x1+≥0,x2-≥0,Maxf-=26x1--6.0x2+,s.t.10x1--12x2+≤3.8,2.8x1-+3.4x2+≤6.5,x1-≥0,x2+≥0.
The result was f+=45.78, x1=1.64, and x2=0.64; f-=30.77, x1=1.37, and x2=0.79.For a detailed description of the problem, see [22].By using the GAILP method as stated in (21), the result can be calculated with the following objective functions:
(24)Maxf+=30x1+-5.5x2+,s.t.8x1+-14x2+≤4.2,2.4x1++3.4x2+≤6.5,x1+≥0,x2+≥0,Maxf-=26x1--6.0x2-,s.t.10x1--12x2-≤3.8,2.8x1-+4x2-≤6.5,x1-≥0,x2-≥0.
The result was f+=48.15, x1=1.73, and x2=0.69; f-=29.15, x1=1.29, and x2=0.72.The GAILP method generates a solution, which is different from that obtained using the interactive binary analysis proposed in [22]. A comparison of the results will be discussed as follows.For thef+ optimistic scheme, the GAILP method can generate a result that is guaranteed to be as close as possible to the upper bound of the constraints. Hence, the maximized value of the objective function is greater than that produced by the interactive binary analysis. For the f- conservative scheme, the GAILP method has a higher probability of satisfying the requirements of the constraints as close as possible to the lowest limit. Hence, the maximized objective value is smaller.In Figures1, 2, 3, and 4 the bold lines denote the boundaries of the constraints, which limit the possible values for x1 and x2 to the left lower area. The constraint a11x1±+a12x2±≤b1 is shown in these figures as the grey bold solid lines, which is the same for both the interactive binary analysis and the GAILP method. The dark bold dotted lines represent the constraint of a21x1±+a22x2±≤b2 given by the interactive binary analysis, and the dark bold solid lines represent the same constraint given by the proposed GAILP method.Figure 1
Optimistic scheme,f+.Figure 2
Zoom-in of the optimistic scheme,f+.Figure 3
Conservative scheme,f-.Figure 4
Zoom-in of the conservative scheme,f-.The boundaries, together with thex1 and x2 axes, enclose the entire area defined by the constraints. The objective functions f+=30x1+-5.5x2+ or f-=26x1--6x2- are groups of parallel lines, as shown in the figures by the thin solid and dotted lines. According to different values of x1 and x2, these objective function lines would have different intercepts on both axes. These constraints restrict the objective function lines to cross with the constraints area, so that, at some vertex, the objective function would reach its extreme (i.e., maximized or minimized) values.In Figures1, 2, 3, and 4, the thin dotted lines are given by the interactive binary analysis, and the thin solid lines represent the objective functions given by the proposed GAILP method. The legend for Figure 1 to Figure 4 is shown in Table 1.Table 1
Legend for Figures1 to 4.
The constrainta21x1±+a22x2±≤b2 given by interactive binary analysis
The constrainta21x1±+a22x2±≤b2 given by GAILP
The constrainta11x1±+a12x2±≤b1
Objective function line given by interactive binary analysis
Objective function line given by GAILP


### 3.2. Genetic-Algorithms-Based Method for Solving IQP Problems (GAIQP)

The GAILP method can be extended to solve the inexact quadratic programming (IQP) problems or other more complicated inexact nonlinear programming problems. The typical IQP problem was presented in Section2 as  (2).In stage one, to obtain an initial suboptimalxjs from a problem transformed from the IQP problem we use the following:
(25)Maxf=∑j=1n[cjrxj+djr(xj)2],s.t.∑j=1naijrxjr≤bjr,i=1,2,…,m,xj≥0,j=1,2,…,n,
where aijr, bir, cjr, and djr are random numbers that satisfy the continuous uniform distribution in the intervals [aij-,aij+], [bi-,bi+], [cj-,cj+], and [dj-,dj+]. Then, a suboptimal solution fscan be identified, and the corresponding decision variables xjs are also obtained.In the second stage, substitutingxjs into the formula in (2) converts (2) into the following formula:
(26)Maxf±=∑j=1n[cj±xjs+dj±(xjs)2],s.t.∑j=1naij±xjs≤bi±,i=1,2.…,m,xjs≥0,j=1,2,…,n.
To determine the coefficients aij±, bi±, cj±, and dj± corresponding to f± we use the following:
(27)Maxf±=∑j=1n[cj±xjs+dj±(xjs)2],s.t.∑j=1naij±xjs≤bi±,i=1,2,…,m,(28)Minf±=∑j=1n[cj±xjs+dj±(xjs)2],s.t.∑j=1naij±xjs≤bi±,i=1,2,…,m.
To determine aj±+ and bi±+ of the optimistic scheme corresponding to the upper limit of the objective value of f+, we have the following:
(29)Max∑j=1nabs(aij±xjs-bi±),s.t.∑j=1naij±xjs≤bi±,i=1,2,…,m.
To obtain aj±- and  bi±-, we use the following:
(30)Min∑j=1nabs(aij±xjs-bi±),s.t.∑j=1naij±xjs≤bi±,i=1,2,…,m.
In the third stage, the problem expressed in (2) has been converted into the following two subproblems.Forf+,
(31)Maxf+=∑j=1n[cj±+xj±+dj±+(xj±)2],s.t.∑j=1naij±+xj±≤bi±+,i=1,2,…,m,xj±≥0,j=1,2,…,n.
For f-,
(32)Maxf-=∑j=1n[cj±-xj±+dj±-(xj±)2],s.t.∑j=1naij±-xj±≤bi±-,i=1,2,…,m,xj±≥0,j=1,2,…,n.
The inexact information has been incorporated in these two subproblems. These two subproblems, as typical nonlinear programming problems, can be solved by the GA nonlinear program solver engine of GASGOT.This method is applied to an IQP problem that was originally proposed by [27]. This IQP problem can be expressed as follows:
(33)Maxf±=[16,18]x1±-[12,14](x1±)2-[4,5]x2±+[14,15](x2±)2,s.t.[4.5,5.5]x1±+[1.8,2.2]x2±≤[1.8,2.1],x1±+[1.8,2.2]x2±≤[0.9,1.1],x1±≥0,x2±≥0.
In stage one, suboptimal variables can be calculated using the GA nonlinear program solver engine of GASGOT:
(34)x1s=0.26497,x2s=0.37772,fs=4.397.
In stage two, x1s and x2s are used to construct the objective functions expressed in (27), (28), (29), and (30) in order to determine the coefficients aij±+, bi±+, cj±+, and dj±+ and aij±-, bi±-, cj±-, and dj±-. By solving (27), (28), (29), and (30), the coefficients are identified as follows:
(35)a11±+=4.5,a12±+=1.8,a21±+=1.8,b1±+=2.1,b2±+=1.1,c1±+=18,c2±+=4,d1±+=12,d2±+=15,a11±-=5.5,a12±-=2.2,a21±-=2.2,b1±-=1.8,b2±-=0.9,c1±-=16,c2±-=5,d1±-=14,d2±-=14.
In stage three, the problems of (31) and (32) are generated as follows:
(36)Maxf+=18x1±-12(x1±)2-4x2±+15(x2±)2,s.t.4.5x1±+1.8x2±≤2.1,x1±+1.8x2±≤1.1,x1±≥0,x2±≥0,Maxf-=16x1±-14(x1±)2-5x2±+14(x2±)2,s.t.5.5x1±+2.2x2±≤1.8,x1±+2.2x2±≤0.9,x1±≥0,x2±≥0.
Solving the above two problems, the solution of this sample problem is f±=[2.6371,5.4234], x1±=[0.2441,0.2857], and x2±=[0.20792,0.45239]. As a comparison, the solution given by [27] is f±=[2.59,5.42],x1±=[0.238,0.286], and x2±=[0.224,0.452].To enhance the solution, the GAIQP engine can be reconfigured in stage three. For example, the maximum genetics generations number can be increased. However, this may not be necessary as the above generated solution is sufficiently satisfactory for many practical engineering problems.


## 3.1. Genetic-Algorithms-Based Method for Solution of ILP Problems (GAILP)

A GA, as a heuristic search algorithm, has been adopted for solving the aforementioned ILP problem. In the GA approach, the upper and lower bounds of the inexact numbers of coefficientsaij±, bi±, and cj± can be determined by substituting the initial suboptimal decision variables into the objective function. f+and f-can be calculated directly without any uncertainty in the coefficients. This approach is called the genetic-algorithms-based method for solving ILP problems or the GAILP method.GAILP has been designed to include three stages, which are discussed as follows.The objective of the first stage is to get an initial suboptimalxjs for the following problem, which is transformed from the ILP problem defined in (1):
(3)Maxf±=∑j=1n[cjrxjs],s.t.∑j=1naijrxjs≤bir,i=1,2,…,m,xj≥0,j=1,2,…,n,
where aijr, bir, and cjr are random numbers that satisfy the continuous uniform distribution in the intervals of [aij-,aij+], [bi-,bi+], and [cj-,cj+], respectively. Then, the problem is solved by the GA linear program solving engine of GASGOT, which uses the objective function in (3) as the positive term of the fitness function and the constraints of (1) as the negative punishment terms. Thus, a suboptimal solution fs can be identified and the corresponding decision variables of xjs are also obtained.In the second stage, the inexact coefficients ofaij±,  bi±, and  cj± will be determined. Let the determined coefficients corresponding tof+ be aij±+, bi±+, and cj±+ and those corresponding to f- be aij±-,  bi±-, and  cj±-. These two sets of coefficients can be obtained using the following method.Substitutingxjs into (1) will convert (1) into the following equation:
(4)Maxf±=∑j=1n[cj±xjs],s.t.∑j=1naij±xjs≤bi±,i=1,2,…,m,xjs≥0,j=1,2,…,n.
To identify the coefficients aij±,  bi±, and cj± corresponding to f±, a set of objective functions needs to be constructed and solved. Since xjs are suboptimal variables, which tend to make the objective function closer to f+, consider aij±,  bi±, and  cj± as variables; then the objective function of (5) can be constructed so as to find cj±+:
(5)Maxf±=∑j=1n[cj±xjs],s.t.∑j=1naij±xjs≤bi±,i=1,2,…,m.
The coefficients cj±+ are considered as corresponding to f+.At the same time, the objective function presented in (6) can be constructed so as to find cj±-:
(6)Minf±=∑j=1n[cj±xjs],s.t.∑j=1naij±xjs≤bi±,i=1,2,…,m.
There are two kinds of decision schemes for inexact programming problems, which are the conservative schemes and optimistic schemes [26]. The former assumes less risk than the latter, so that, for a maximization objective function, planning for the lower bound of an objective value represents the conservative scheme, and planning for the upper bound of an objective value represents the optimistic scheme [26]. In terms of constraints, the conservative scheme involves more rigorous or stringent constraints, and the optimistic scheme adopts more tolerant ones.Thus, the problem of searching foraij±+ and  bi±+ of the optimistic scheme and corresponding to the upper bound of the objective value of f+ can be represented as follows:
(7)Max∑j=1nabs(aij±xjs-bi±),s.t.∑j=1naij±xjs≤bi±,i=1,2,…,m.
The problem
(8)Min∑j=1nabs(aij±xjs-bi±),s.t.∑j=1naij±xjs≤bi±,i=1,2,…,m,
will give aij±- and bi±- of the conservative scheme, corresponding to the lower bound of the objective value of f-.Hence, the values ofaij±+, bi±+, and cj±+and aij±-, bi±-, and cj±- can be calculated.In the third stage, the problem represented in (1) is converted into the following two subproblems.Forf+,
(9)Maxf+=∑j=1n[cj±+xj±],s.t.∑j=1naij±+xj±≤bi±+,i=1,2,…,m,xj±≥0,j=1,2,…,n.Forf-,
(10)Maxf-=∑j=1n[cj±-xj±],s.t.∑j=1naij±-xj±≤bi±-,i=1,2,…,m,xj±≥0,j=1,2,…,n.
This step eliminates the inexact parameters in (1) and generates instead (9) and (10) as typical linear programming (LP) problems, which can be solved easily using the traditional methods. Generally speaking, the interactive binary algorithm of Huang et al.  [20, 22] can be used for solving inexact linear problems reliably and relatively quickly for many real-life decision-making scenarios in the engineering field. However, this binary algorithm has some limitations. One of them, for example, is the limitation that the upper and lower bounds of an inexact coefficient cannot have different signs. By contrast, the GAILP does not have this kind of limitation because the GA method does not depend on any assumed distribution of the inexact parameter.The following inexact linear problem demonstrates how GAILP is able to handle this situation:(11)Maxf±=[2,6]x1±+[7.8,10.2]x2±,s.t.x1±+[-0.8,1.2]x2±≤4.6,[-6,-8.5]x1±+12x2±≤[-2,3],x1±≥0,x2±≥0.
In stage one, the suboptimal fs and the corresponding x1s and x2s are calculated:
(12)fs=112,x1s=10.1,x2s=5.8.
In stage two, substituting x1s and x2s into (11), the two sets of problems for determining aij±+,  bi±+, and  cj±+ and aij±-,  bi±-, and  cj±- are constructed forf+ and f-, respectively.Forf+, cj±+can be determined by solving the following problem:
(13)Maxf+=c1±+10.1+c2±+5.8,s.t.10.1+[-0.8,1.2]5.8≤4.6,[-6,-8.5]10.1+125.8≤[-2,3],c1±+≤[2,6],c2±+≤[7.8,10.2].
And aij±+ and bi±+ can be determined by solving the following problem:
(14)Max4.6-(x1s+a12x2s),Maxb2-(a21x1s+12x2s),a12=[-0.8,1.2],a21=[-6,-8.5],b2=[-2,3].
For f-, cj±-can be determined by solving the following problem:
(15)Minf-=c1±-10.1+c2±-5.8,s.t.10.1+[-0.8,1.2]5.8≤4.6,[-6,-8.5]10.1+125.8≤[-2,3],c1±-≤[2,6],c2±-≤[7.8,10.2].
And aij±- and bi±- can be determined by solving the following problem:
(16)Min4.6-(x1s+a12x2s),Minb2-(a21x1s+12x2s),a12=[-0.8,1.2],a21=[-8.5,-6],b2=[-2,3].
By solving (13), (14), (15), and (16), the coefficients are calculated as follows:
(17)a12±+=-0.8,a21±+=-8.5,b2±+=3,c1±+=6,c2±+=10.2,a12±-=1.2,a21±-=-6,b2±-=-2,c1±-=2,c2±-=7.8.
In stage three, the two submodels for the optimistic scheme and the conservative scheme are constructed, respectively, as follows.For the optimistic scheme,f+,
(18)Maxf+=6x1++10.2x2+,s.t.x1+-0.8x2+≤4.6,-8.5x1++12x2+≤3,x1+≥0,x2+≥0.
For the conservative scheme, f-,
(19)Maxf-=2x1-+7.8x2-,s.t.x1-+1.2x2-≤4.6,-6x1-+12x2-≤2,x1-≥0,x2-≥0.
The final result can be found by solving the previous two submodels:
(20)f+=149.04,x1+=11.08,x2+=8.1,f-=16.4,x1+=3,x2+=1.33.
From the previous calculation, it can be seen that the GAILP method can be used without any assumption of the upper and lower bounds of the inexact coefficients. In fact, this method effectively extends the scope of problems solvable using the methods of the inexact linear programming problem. Therefore, the GAILP method is more adaptable for real world applications of optimization problems with uncertainty. In the next section, this method will be extended to solve the inexact quadratic problems.A sample inexact linear programming problem in [22] is as follows:
(21)Maxf±=c1x1±+c2x2±,s.t.a11x1±+a12x2±≤b1,a21x1±+a22x2±≤b2,
where c1=[26,30],c2=[-6,-5.5],
(22)a11=[8,10],a12=[-14,-12],b1=[3.8,4.2],a21=[2.4,2.8],a22=[3.4,4],b2=6.5.
By using the traditional interactive binary algorithm, two submodels are obtained:
(23)Maxf+=30x1+-5.5x2-,s.t.8x1+-14x2-≤4.2,2.4x1++4x2-≤6.5,x1+≥0,x2-≥0,Maxf-=26x1--6.0x2+,s.t.10x1--12x2+≤3.8,2.8x1-+3.4x2+≤6.5,x1-≥0,x2+≥0.
The result was f+=45.78, x1=1.64, and x2=0.64; f-=30.77, x1=1.37, and x2=0.79.For a detailed description of the problem, see [22].By using the GAILP method as stated in (21), the result can be calculated with the following objective functions:
(24)Maxf+=30x1+-5.5x2+,s.t.8x1+-14x2+≤4.2,2.4x1++3.4x2+≤6.5,x1+≥0,x2+≥0,Maxf-=26x1--6.0x2-,s.t.10x1--12x2-≤3.8,2.8x1-+4x2-≤6.5,x1-≥0,x2-≥0.
The result was f+=48.15, x1=1.73, and x2=0.69; f-=29.15, x1=1.29, and x2=0.72.The GAILP method generates a solution, which is different from that obtained using the interactive binary analysis proposed in [22]. A comparison of the results will be discussed as follows.For thef+ optimistic scheme, the GAILP method can generate a result that is guaranteed to be as close as possible to the upper bound of the constraints. Hence, the maximized value of the objective function is greater than that produced by the interactive binary analysis. For the f- conservative scheme, the GAILP method has a higher probability of satisfying the requirements of the constraints as close as possible to the lowest limit. Hence, the maximized objective value is smaller.In Figures1, 2, 3, and 4 the bold lines denote the boundaries of the constraints, which limit the possible values for x1 and x2 to the left lower area. The constraint a11x1±+a12x2±≤b1 is shown in these figures as the grey bold solid lines, which is the same for both the interactive binary analysis and the GAILP method. The dark bold dotted lines represent the constraint of a21x1±+a22x2±≤b2 given by the interactive binary analysis, and the dark bold solid lines represent the same constraint given by the proposed GAILP method.Figure 1
Optimistic scheme,f+.Figure 2
Zoom-in of the optimistic scheme,f+.Figure 3
Conservative scheme,f-.Figure 4
Zoom-in of the conservative scheme,f-.The boundaries, together with thex1 and x2 axes, enclose the entire area defined by the constraints. The objective functions f+=30x1+-5.5x2+ or f-=26x1--6x2- are groups of parallel lines, as shown in the figures by the thin solid and dotted lines. According to different values of x1 and x2, these objective function lines would have different intercepts on both axes. These constraints restrict the objective function lines to cross with the constraints area, so that, at some vertex, the objective function would reach its extreme (i.e., maximized or minimized) values.In Figures1, 2, 3, and 4, the thin dotted lines are given by the interactive binary analysis, and the thin solid lines represent the objective functions given by the proposed GAILP method. The legend for Figure 1 to Figure 4 is shown in Table 1.Table 1
Legend for Figures1 to 4.
The constrainta21x1±+a22x2±≤b2 given by interactive binary analysis
The constrainta21x1±+a22x2±≤b2 given by GAILP
The constrainta11x1±+a12x2±≤b1
Objective function line given by interactive binary analysis
Objective function line given by GAILP


## 3.2. Genetic-Algorithms-Based Method for Solving IQP Problems (GAIQP)

The GAILP method can be extended to solve the inexact quadratic programming (IQP) problems or other more complicated inexact nonlinear programming problems. The typical IQP problem was presented in Section2 as  (2).In stage one, to obtain an initial suboptimalxjs from a problem transformed from the IQP problem we use the following:
(25)Maxf=∑j=1n[cjrxj+djr(xj)2],s.t.∑j=1naijrxjr≤bjr,i=1,2,…,m,xj≥0,j=1,2,…,n,
where aijr, bir, cjr, and djr are random numbers that satisfy the continuous uniform distribution in the intervals [aij-,aij+], [bi-,bi+], [cj-,cj+], and [dj-,dj+]. Then, a suboptimal solution fscan be identified, and the corresponding decision variables xjs are also obtained.In the second stage, substitutingxjs into the formula in (2) converts (2) into the following formula:
(26)Maxf±=∑j=1n[cj±xjs+dj±(xjs)2],s.t.∑j=1naij±xjs≤bi±,i=1,2.…,m,xjs≥0,j=1,2,…,n.
To determine the coefficients aij±, bi±, cj±, and dj± corresponding to f± we use the following:
(27)Maxf±=∑j=1n[cj±xjs+dj±(xjs)2],s.t.∑j=1naij±xjs≤bi±,i=1,2,…,m,(28)Minf±=∑j=1n[cj±xjs+dj±(xjs)2],s.t.∑j=1naij±xjs≤bi±,i=1,2,…,m.
To determine aj±+ and bi±+ of the optimistic scheme corresponding to the upper limit of the objective value of f+, we have the following:
(29)Max∑j=1nabs(aij±xjs-bi±),s.t.∑j=1naij±xjs≤bi±,i=1,2,…,m.
To obtain aj±- and  bi±-, we use the following:
(30)Min∑j=1nabs(aij±xjs-bi±),s.t.∑j=1naij±xjs≤bi±,i=1,2,…,m.
In the third stage, the problem expressed in (2) has been converted into the following two subproblems.Forf+,
(31)Maxf+=∑j=1n[cj±+xj±+dj±+(xj±)2],s.t.∑j=1naij±+xj±≤bi±+,i=1,2,…,m,xj±≥0,j=1,2,…,n.
For f-,
(32)Maxf-=∑j=1n[cj±-xj±+dj±-(xj±)2],s.t.∑j=1naij±-xj±≤bi±-,i=1,2,…,m,xj±≥0,j=1,2,…,n.
The inexact information has been incorporated in these two subproblems. These two subproblems, as typical nonlinear programming problems, can be solved by the GA nonlinear program solver engine of GASGOT.This method is applied to an IQP problem that was originally proposed by [27]. This IQP problem can be expressed as follows:
(33)Maxf±=[16,18]x1±-[12,14](x1±)2-[4,5]x2±+[14,15](x2±)2,s.t.[4.5,5.5]x1±+[1.8,2.2]x2±≤[1.8,2.1],x1±+[1.8,2.2]x2±≤[0.9,1.1],x1±≥0,x2±≥0.
In stage one, suboptimal variables can be calculated using the GA nonlinear program solver engine of GASGOT:
(34)x1s=0.26497,x2s=0.37772,fs=4.397.
In stage two, x1s and x2s are used to construct the objective functions expressed in (27), (28), (29), and (30) in order to determine the coefficients aij±+, bi±+, cj±+, and dj±+ and aij±-, bi±-, cj±-, and dj±-. By solving (27), (28), (29), and (30), the coefficients are identified as follows:
(35)a11±+=4.5,a12±+=1.8,a21±+=1.8,b1±+=2.1,b2±+=1.1,c1±+=18,c2±+=4,d1±+=12,d2±+=15,a11±-=5.5,a12±-=2.2,a21±-=2.2,b1±-=1.8,b2±-=0.9,c1±-=16,c2±-=5,d1±-=14,d2±-=14.
In stage three, the problems of (31) and (32) are generated as follows:
(36)Maxf+=18x1±-12(x1±)2-4x2±+15(x2±)2,s.t.4.5x1±+1.8x2±≤2.1,x1±+1.8x2±≤1.1,x1±≥0,x2±≥0,Maxf-=16x1±-14(x1±)2-5x2±+14(x2±)2,s.t.5.5x1±+2.2x2±≤1.8,x1±+2.2x2±≤0.9,x1±≥0,x2±≥0.
Solving the above two problems, the solution of this sample problem is f±=[2.6371,5.4234], x1±=[0.2441,0.2857], and x2±=[0.20792,0.45239]. As a comparison, the solution given by [27] is f±=[2.59,5.42],x1±=[0.238,0.286], and x2±=[0.224,0.452].To enhance the solution, the GAIQP engine can be reconfigured in stage three. For example, the maximum genetics generations number can be increased. However, this may not be necessary as the above generated solution is sufficiently satisfactory for many practical engineering problems.


## 4. Case Study

To illustrate the proposed method, the problem of solid waste disposal planning presented in [25] has been recalculated using the GAIQP method. In this case study, the system involves three cities. As shown in Figure 5, the planning horizon is 15 years, which is divided equally into three periods. A landfill and an incinerator are available for the disposal of the municipal solid waste. The landfill has an existing capacity of [2.05,2.30]×106 t, and the incinerator has a capacity of [500,600] t/d. The incinerator generates residues of approximately 30% of the incoming waste streams, and its revenue from energy sale is $[15,25] per ton combusted. The waste generation rate of each city, operating costs of the facilities, and the waste transportation costs are summarized and shown in Table 2.Table 2
Waste generation rates and operating costs.
Period 1 (k=1)
Period 2 (k=2)
Period 3 (k=3)
Waste generation (t/d)
W
G
1
k
±
[260, 340]
[310, 390]
[360, 440]
W
G
2
k
±
[160, 240]
[185, 265]
[210, 290]
W
G
3
k
±
[260, 340]
[260, 340]
[310, 390]
Operation cost ($/t)
O
P
1
k
±
[35, 45]
[40, 60]
[50, 80]
O
P
2
k
±
[55, 75]
[60, 85]
[65, 95]
City-to-landfill transportation cost ($/t)
T
R
11
k
-
-
0.0123
x
+
14.58
-
0.0135
x
+
16.04
-
0.0148
x
+
17.64
T
R
11
k
+
-
0.0163
x
+
19.40
-
0.0179
x
+
21.34
-
0.0197
x
+
23.48
T
R
12
k
-
-
0.0106
x
+
12.65
-
0.0117
x
+
13.92
-
0.0129
x
+
15.31
T
R
12
k
+
-
0.0142
x
+
16.87
-
0.0156
x
+
18.56
-
0.0172
x
+
20.41
T
R
13
k
-
-
0.0129
x
+
15.30
-
0.0141
x
+
16.83
-
0.0156
x
+
18.52
T
R
13
k
+
-
0.0172
x
+
20.49
-
0.0189
x
+
22.53
-
0.0208
x
+
24.79
City-to-incinerator transportation cost ($/t)
T
R
21
k
-
-
0.0097
x
+
11.57
-
0.0107
x
+
12.73
-
0.0118
x
+
14.00
T
R
21
k
+
-
0.0130
x
+
15.42
-
0.0143
x
+
16
/
97
-
0.0157
x
+
18.66
T
R
22
k
-
-
0.0102
x
+
12.17
-
0.0113
x
+
13.39
-
0.0124
x
+
14.73
T
R
22
k
+
-
0.0136
x
+
16.15
-
0.0149
x
+
17.76
-
0.0164
x
+
19.54
T
R
23
k
-
-
0.0089
x
+
10.60
-
0.0098
x
+
11.67
-
0.0108
x
+
12.83
T
R
23
k
+
-
0.0118
x
+
14.10
-
0.0130
x
+
15.51
-
0.0143
x
+
17.06
Incinerator-to-landfill transportation cost ($/t)
F
T
k
-
-
0.0048
x
+
5.71
-
0.0053
x
+
6.28
-
0.0058
x
+
6.91
F
T
k
+
-
0.0064
x
+
7.62
-
0.0070
x
+
8.38
-
0.0077
x
+
9.33Figure 5
Diagram of municipalities and waste management facilities in the case study.The objective of the optimization problem in this case study is to minimize the total costs by allocating waste flow between cities and facilities. In [25], this IQP model was formulated as follows:
(37)Minf±=∑i=12∑j=13∑k=13Lk(TRijk±+OPik±)xijk±+∑j=13∑k=13Lk[FE(FTk±+OP2k±)-REk±]x2jk±(whereTRijk±=αijk±xijk±+βijk±,FTjk±=σk±x2jk±FE+δk±),s.t.∑j=13∑k=13Lk(x1jk±+FEx2jk±)≤TL±(landfill  capacity  constraint),∑j=13x2jk±≤TE±∀k(incineratorcapacity  constraint),∑i=12xijk±≥WGjk±∀j,k(waste  disposaldemand  constraint),xijk±≥0∀i,j,k(nonnegativity  constraint),
whereFE is residue flow rate from incinerator to landfill (it is 0.3 in this case),FTk± is transportation cost for residue from incinerator to landfill during period k ($/t),i is type of waste management facility (i=1,2, where i=1 for landfill, 2 for incinerator),j is city, j=1,2,3,k is time period, k=1,2,3,Lk is length of period k, L1=L2=L3=365*5 (day),OPik± is operating cost of facility i during period k ($/t),REk± is revenue from incinerator during period k ($/t),RE1±=RE2±=RE3±=[15,25],TE± is capacity of incinerator (t/d),TL± is capacity of landfill (t),TLijk± is transportation cost for waste from city j to facility i during period k ($/t),WGjk± is waste generation rate in city j during period k (t/d),xijk± is waste flow from city j to facility i during period k (t/d),αijk± is slope of transportation cost curve for waste from city j to facility i, during period k,βijk± is Y-intercept of transportation cost curve for waste from city j to facility i during period k,δk± is Y-intercept of transportation cost curve for residue from incinerator to landfill during period k,σk± is Slope of transportation cost curve for residue from incinerator to landfill during period k.This interval quadratic programming problem has 18 variables xijk and 13 constraints. The GAIQP method is applied to solve this problem according to the three stages presented in Section 3.2.Stage  1. Based on the data shown in Table 2, the model formulated at stage 1 is as follows:
(38)minf=365×5×{x2332[44.6,64.4]x111-[0.0123,0.0163]x1112+[56.04,81.34]x112-[0.0135,0.0179]x1122+[67.64,103.48]x113-[0.0148,0.0197]x1132+[42.65,61.9]x121-[0.0106,0.0142]x1212+[53.9,78.56]x122-[0.0117,0.0156]x1222+[65.3,100.4]x123-[0.0129,0.0172]x1232+[45.3,65.5]x131-[0.0129,0.0172]x1312+[56.83,82.53]x132-[0.0141,0.0189]x1322+[58.52,104.79]x133-[0.0156,0.0208]x1332+[60.6,100.4]x211-[0.0111,0.015]x2112+[64.2,115.3]x212-[0.0122,0.0164]x2122+[74.4,126.73]x213-[0.0134,0.018]x2132+[54.3,103.44]x221-[0.0117,0.0157]x2212+[66.2,102.76]x222-[0.013,0.0172]x2222+[72.13,114.54]x223-[0.0143,0.0186]x2232+[60.3,99.1]x231-[0.0102,0.0137]x2312+[65.41,108.7]x232-[0.0104,0.0148]x2322+[71.38,133.25]x233-[0.012,0.0143]x2332},s.t.365×5×(x111+0.3x211+x112+0.3x212+x113+0.3x213+x121+0.3x221+x122+0.3x222+x123+0.3x223+x131+0.3x231+x132+0.3x232+x123+0.3x233)≤[2.05×106,2.30×106],x211+x221+x231≤[500,600],x212+x222+x232≤[500,600],x213+x223+x233≤[500,600],x111+x211≥[260,340],x112+x212≥[310,390],x113+x213≥[360,440],x121+x221≥[160,240],x122+x222≥[185,265],x131+x231≥[260,340],x132+x232≥[260,340],x133+x233≥[310,390].
Considering interval numbers as random numbers, which can be determined between the interval numbers’ lower and upper endpoints, the GAIQP method is applied to find a suboptimal solution fs. With an initial population size of 500, after 2000 generations, one suboptimal solution has been found: fs=288×106, and the concomitant variables are
(39)x111s=250,x112s=315,x113s=376,x121s=3,x122s=220,x123s=57,x131s=8,x132s=10,x133s=5,x211s=45,x212s=40,x213s=10,x221s=188,x222s=12,x223s=201,x231s=278,x232s=293,x233s=322.Stage  2. Substitute xijks into (38) to identify those coefficients corresponding to f+ and f-.Stage  3. In this stage, the problem expressed in (37) has been converted into two subproblems: f+ and f-.The solution found by this method isf±=[2.405×108,5.133×108]. This is close to the result given by [25], which is, for comparison purposes, f±=[2.39×108,5.14×108]. This case study indicates that the GAIQP method can be configured to deal with large-scaled and complex engineering problems and give a satisfactory solution. If the parameters of the genetic algorithms of the nonlinear program solver engine of GASGOT are further tuned, a better solution can be generated.


## 5. Conclusions

Two genetic-algorithms-based methods have been proposed and applied for identifying an all-purpose optimization solution for the inexact linear programming and Inexact quadratic programming problems. The two methods are called GAILP and GAIQP. The Genetic Algorithm Solver of MATLAB was the implementation environment of the proposed methods. Compared to the GAILP and GAIQP methods, the traditional problem-solving method has limitations due to the complexity involved in selecting the upper or lower bounds of variables and parameters when the subobjective functions are being constructed. The complexity arises due to the extensive computation and necessary assumptions and simplification. The solution procedures of the proposed GA-based optimization methods do not involve any such assumption or simplification, and the quality of the result is guaranteed. The GAIQP method has been applied to a case study that deals with municipal solid waste management taken from [25]. A comparison of the results shows that the proposed GA-based heuristic optimization approach is able to handle more complicated quadratic relationships involving uncertainty and provide better results.The GA-based heuristic optimization approach is a flexible approach, which can be extended to find solutions for various types of operation programming scenarios. It can also be used as an all-purpose algorithm for economic optimizations. In the future, the methods of GAILP and GAIOP will be further developed for handling more complex inexact nonlinear problems.



---

*Source: 272491-2013-12-12.xml*