Research Article: Extension of Wolfe Method For Solving Quadratic Programming With Interval Coefficients
Research Article: Extension of Wolfe Method For Solving Quadratic Programming With Interval Coefficients
Research Article
Extension of Wolfe Method for Solving Quadratic
Programming with Interval Coefficients
Copyright © 2017 Syaripuddin et al. This is an open access article distributed under the Creative Commons Attribution License,
which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
Quadratic programming with interval coefficients developed to overcome cases in classic quadratic programming where the
coefficient value is unknown and must be estimated. This paper discusses the extension of Wolfe method. The extended Wolfe
method can be used to solve quadratic programming with interval coefficients. The extension process of Wolfe method involves the
transformation of the quadratic programming with interval coefficients model into linear programming with interval coefficients
model. The next step is transforming linear programming with interval coefficients model into two classic linear programming
models with special characteristics, namely, the optimum best and the worst optimum problem.
programming with interval coefficients. Furthermore, the The model in (2a)–(2c) is solved by means of transform-
linear programming with interval coefficients which has been ing the linear programming with interval coefficients into two
obtained from the transformation will be solved by using the classic linear programming models with the special charac-
method in [8]. teristics, namely, the best optimum and the worst optimum
This paper is organized as follows. Section 2 discusses problems. The best optimum problem has properties of best
interval arithmetic operations. In Section 3, a general form version on the objective function and maximum feasible area
of linear programming with interval coefficients is stated. In on the constraint function. On the other hand, the worst
Section 4, a general form of quadratic programming with optimum problem has a characteristic that it is the worst
interval coefficients is stated. Extension of Wolfe method version of the objective function and the minimum feasible
as one method of solving the quadratic programming with area of the constraint function.
interval coefficients is discussed in Section 5, whereas Sec- Chinneck and Ramadan [8] provide a rule to determine
tion 6 discusses numerical examples, and Section 7 provides the best and worst optimum problem in a linear program-
some concluding remarks. ming problem with interval coefficients. The constraints of
linear programming with interval coefficients which have
an inequality sign (≤) in (2b) have characteristics of the
2. Interval Arithmetic
maximum feasible area and the minimum feasible area which
The basic definition and properties of interval number and is given by the following theorem.
interval arithmetic can be seen at Moore [3], Alefeld and
Herzberg [11], and Hansen [12]. Theorem 4 (Chinneck and Ramadan [8]). Suppose that we
have an interval inequality given ∑𝑛𝑗=1 [𝑎𝑗𝐼 , 𝑎𝑗𝑆 ]𝑥𝑗 ≤ [𝑏𝐼 , 𝑏𝑆 ],
Definition 1. A closed real interval 𝑥 = [𝑥𝐼 , 𝑥𝑆 ] denoted by 𝑥 where 𝑥𝑗 ≥ 0. Then, ∑𝑛𝑗=1 𝑎𝑗𝐼 𝑥𝑗 ≤ 𝑏𝑆 is maximum feasible area
is a real interval number which can be defined completely by and ∑𝑛𝑗=1 𝑎𝑗𝑆 𝑥𝑗 ≤ 𝑏𝐼 is minimum feasible area.
𝑥 = [𝑥𝐼 , 𝑥𝑆 ] = {𝑥 ∈ R | 𝑥𝐼 ≤ 𝑥 ≤ 𝑥𝑆 ; 𝑥𝐼 , 𝑥𝑆 ∈ R} , (1) The objective function of linear programming with
interval coefficients for the case of maximizing (2a) has
where 𝑥𝐼 and 𝑥𝑆 are called infimum and supremum, respec- characteristics of the best version and the worst version of the
tively. objective function is expressed in the following theorem.
Definition 2. A real interval number 𝑥 = [𝑥𝐼 , 𝑥𝑆 ] is called Theorem 5 (Chinneck and Ramadan [8]). If 𝑧 = ∑𝑛𝑗=1 [𝑐𝑗𝐼 ,
degenerate, if 𝑥𝐼 = 𝑥𝑆 .
𝑐𝑗𝑆 ]𝑥𝑗 is the objective function for 𝑥𝑗 ≥ 0, then ∑𝑛𝑗=1 𝑐𝑗𝑆 𝑥𝑗 ≥
Definition 3. Let 𝑥 = [𝑥𝐼 , 𝑥𝑆 ] and 𝑦 = [𝑦𝐼 , 𝑦𝑆 ]; then ∑𝑛𝑗=1 𝑐𝑗𝐼 𝑥𝑗 , where ∑𝑛𝑗=1 𝑐𝑗𝑆 𝑥𝑗 is the best version of the objective
function and ∑𝑛𝑗=1 𝑐𝑗𝐼 𝑥𝑗 is the worst version of the objective
1. 𝑥 + 𝑦 = [𝑥𝐼 + 𝑦𝐼 , 𝑥𝑆 + 𝑦𝑆 ] (addition), function.
2. 𝑥 − 𝑦 = [𝑥𝐼 , 𝑥𝑆 ] − [𝑦𝐼 , 𝑦𝑆 ] = [𝑥𝐼 , 𝑥𝑆 ] + [−𝑦𝑆 , −𝑦𝐼 ] =
[𝑥𝐼 − 𝑦𝑆 , 𝑥𝑆 − 𝑦𝐼 ] (subtraction),
4. Quadratic Programming with
3. 𝑥 ⋅ 𝑦 = [min{𝑥𝐼 𝑦𝐼 , 𝑥𝐼 𝑦𝑆 , 𝑥𝑆 𝑦𝐼 , 𝑥𝑆 𝑦𝑆 }, max{𝑥𝐼 𝑦𝐼 , 𝑥𝐼 𝑦𝑆 , Interval Coefficients
𝑥𝑆 𝑦𝐼 , 𝑥𝑆 𝑦𝑆 }] (multiplication),
The general form of quadratic programming with interval
4. 𝑥/𝑦 = 𝑥(1/𝑦) = [𝑥𝐼 , 𝑥𝑆 ][1/𝑦𝑆 , 1/𝑦𝐼 ], 0 ∉ 𝑦 (division).
coefficients introduced by Li and Tian [5] is defined as
follows:
3. Linear Programming with
Interval Coefficients Maximize 𝑧
The general form of linear programming with interval coeffi- 𝑛
where 𝑥𝑗 ∈ R, [𝑐𝑗𝐼 , 𝑐𝑗𝑆 ], [𝑏𝑖𝐼 , 𝑏𝑖𝑆 ], and [𝑎𝑖𝑗𝐼 , 𝑎𝑖𝑗𝑆 ] ∈ 𝐼(R). 𝑥𝑗 ≥ 0, 𝑗 = 1, 2, . . . , 𝑛, (3c)
Journal of Applied Mathematics 3
where 𝑥𝑗 ∈ R, [𝑐𝑗𝐼 , 𝑐𝑗𝑆 ], [𝑞𝑖𝑗𝐼 , 𝑞𝑖𝑗𝑆 ], [𝑏𝑖𝐼 , 𝑏𝑖𝑆 ], [𝑎𝑖𝑗𝐼 , 𝑎𝑖𝑗𝑆 ] ∈ 𝐼(R), to the variables and equating to zero (KKT necessary condi-
∑𝑛𝑗=1 ∑𝑛𝑘=1 [𝑞𝑗𝑘𝐼 , 𝑞𝑗𝑘𝑆 ]𝑥𝑗 𝑥𝑘 is negative semidefinite, and 𝐼(R) tions) (see [13, 14]).
are the set of all interval numbers in R. 𝑛
The model as shown in (3a)–(3c) is a generalization of 𝜕𝐿
= [𝑐𝑗𝐼 , 𝑐𝑗𝑆 ] + ∑ [𝑞𝑗𝑘𝐼 , 𝑞𝑗𝑘𝑆 ] 𝑥𝑘
the model in [4]. The coefficient of the objective function 𝜕𝑥𝑖 𝑘=1
and constraints of the quadratic programming with interval
𝑚
coefficients model in (3a)–(3c) have an interval form. The (5a)
− ∑ [𝑎𝑖𝑗𝐼 , 𝑎𝑖𝑗𝑆 ] 𝜆 𝑖 + 𝜇𝑗 = 0,
idea to solve the model is the extension of Wolfe method. 𝑖=1
This method focuses on how to transform the quadratic
programming with interval coefficients in (3a)–(3c) into 𝑗 = 1, 2, . . . , 𝑛,
linear programming with interval coefficients in (2a)–(2c).
Furthermore, linear programming with interval coefficients 𝜕𝐿
= −2𝜆 𝑖 𝑦𝑖 = 0, 𝑖 = 1, 2, . . . , 𝑚, (5b)
obtained from the transformation is done using the method 𝜕𝑦𝑖
in [8].
𝜕𝐿
Extension of Wolfe method is the main result in this = −2𝜇𝑗 𝑟𝑗 = 0, 𝑗 = 1, 2, . . . , 𝑛, (5c)
𝜕𝑟𝑖
paper. The fundamental difference between the extensions of
Wolfe method and the method in [5] is, on the extension of 𝑛
𝜕𝐿
Wolfe method, quadratic programming model with interval = ∑ [𝑎𝑖𝑗𝐼 , 𝑎𝑖𝑗𝑠 ] 𝑥𝑗 + 𝑦𝑖 2 − [𝑏𝑖𝐼 , 𝑏𝑖𝑆 ]
𝜕𝜆 𝑖 𝑗=1 (5d)
coefficients is transformed into linear programming with
interval coefficients, while, in [5], the model of quadratic
programming with interval coefficients is maintained. = 0, 𝑖 = 1, 2, . . . , 𝑚,
𝜕𝐿
= 𝑥𝑗 − 𝑟𝑗2 = 0, 𝑗 = 1, 2, . . . , 𝑛, (5e)
5. Extension of Wolfe Method 𝜕𝜇𝑖
(1) Given a quadratic programming problem with inter- (3) The optimum value of the quadratic programming
val coefficients in (3a)–(3c), Extension of Wolfe with interval coefficients is obtained by combining the
method is based on (3a)–(3c) equivalent to the linear optimum value from the worst and the best optimum
programming with interval coefficients in (8a)–(8e). problem; that is, 𝑧 = [𝑧𝐼 , 𝑧𝑆 ].
(2) Use Theorems 4 and 5 for transforming the linear
programming with interval coefficients in (8a)–(8e) Algorithm 6 shows that the best and the worst optimum
into two classic linear programming models with problem are linear programming models added by comple-
special characteristics; namely, mentary conditions. Thus, both problems can be solved by
simplex method.
(a) the best optimum problem is
Minimize 𝑧𝑆 = V1 + V2 + ⋅ ⋅ ⋅ + V𝑛 (9a) 6. Numerical Example
According to Li and Tian [5], for the solution of the model in [−1, 1] 𝑥1 + [20, 40] 𝑥2 + [3, 3] 𝜆 1 + [4, 6] 𝜆 2 − 𝜇2
(11a)–(11c), the best optimum problem is 𝑧𝑆 = −0.9, 𝑥1 = 0.3, (13c)
+ V2 = [−3, −2]
and 𝑥2 = 0, the worst optimum problem is 𝑧𝐼 = −6.25, 𝑥1 =
1, 25, and 𝑥2 = 0, and the optimum value is 𝑧 = [𝑧𝐼 , 𝑧𝑆 ] = [1, 2] 𝑥1 + 3𝑥2 ≤ [1, 10] (13d)
[−6.25, −0.9].
This paper presents only the maximization problem [−2, 8] 𝑥1 + [4, 6] 𝑥2 ≤ [4, 6] (13e)
so that any minimization problem will be converted into
maximization problem, the simple procedure to convert a 𝑥1 , 𝑥2 ≥ 0. (13f)
minimization problem to a maximization problem and vice
versa. Simply multiply the objective function of a minimiza- We apply Algorithm 6 for transforming linear programming
tion problem by −1 converting it into a maximization problem with interval coefficients model in ((13a)–(13f)) into two clas-
and vice versa. sic linear programming models with special characteristics,
namely, the best optimum and the worst optimum problem.
Maximize 𝑧 The result of the transformation is shown in Table 1.
So, the optimum value of the quadratic programming
= [6, 10] 𝑥1 + [−3, −2] 𝑥2 with interval coefficients is obtained by combining the opti-
(12a) mum value from the worst and the best optimum problem;
+ [−1, 1] 𝑥1 𝑥2 + [−10, −4] 𝑥12 that is, 𝑧 = [𝑧𝐼 , 𝑧𝑆 ] = [0.9, 6.25]. This solution gives the same
value as obtained by Li and Tian [5].
+ [−20, −10] 𝑥22
7. Conclusion
subject to
This paper presents an extension of Wolfe method. The
[1, 2] 𝑥1 + 3𝑥2 ≤ [1, 10] (12b) extension of Wolfe method performed by transforming the
quadratic programming with interval coefficients model
[−2, 8] 𝑥1 + [4, 6] 𝑥2 ≤ [4, 6] (12c) into linear programming with interval coefficients model.
Furthermore, linear programming with interval coefficients
𝑥1 , 𝑥2 ≥ 0. model is transformed into two classic linear programming
(12d)
models using Algorithm 6. The extension of Wolfe method
has a particular benefit: the final model is linear program-
We apply the extension of Wolfe method for transforming
ming. Hence, it can be solved by the simplex method.
quadratic programming with interval coefficients model in
((12a)–(12d)) into linear programming with interval coeffi-
cients model. We have Conflicts of Interest
Minimize 𝑧 = V1 + V2 (13a) The authors declare that there are no conflicts of interest
regarding the publication of this paper.
subject to
References
[8, 20] 𝑥1 + [−1, 1] 𝑥2 + [1, 2] 𝜆 1 + [−2, 8] 𝜆 2 − 𝜇1
(13b) [1] F. S. Hillier and G. J. Lieberman, Introduction to operations
+ V1 = [6, 10] research, Holden-Day, Inc., Oakland, Calif., Third edition, 1980.
6 Journal of Applied Mathematics
International
Journal of Journal of
Mathematics and
Mathematical
#HRBQDSDĮ,@SGDL@SHBR
Sciences