Next Article in Journal
An Algorithm for the Numerical Integration of Perturbed and Damped Second-Order ODE Systems
Previous Article in Journal
A Risk-Aversion Approach for the Multiobjective Stochastic Programming Problem
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Stability of Solutions for Parametric Inverse Nonlinear Cost Transportation Problem

1
Department of Mathematics and Statistics, College of Science, Taif University, P.O. Box 11099, Taif 21944, Saudi Arabia
2
Department of Basic Engineering Science, Faculty of Engineering, Menofia University, Shebin El-Kom 32511, Egypt
3
Department of Basic Science, Higher Technological Institute, Tenth of Ramadan City 44629, Egypt
*
Author to whom correspondence should be addressed.
Mathematics 2020, 8(11), 2027; https://doi.org/10.3390/math8112027
Submission received: 30 September 2020 / Revised: 7 November 2020 / Accepted: 10 November 2020 / Published: 14 November 2020

Abstract

:
This paper investigates the solution for an inverse of a parametric nonlinear transportation problem, in which, for a certain values of the parameters, the cost of the unit transportation in the basic problem are adapted as little as possible so that the specific feasible alternative become an optimal solution. In addition, a solution stability set of these parameters was investigated to keep the new optimal solution (feasible one) is unchanged. The idea of this study based on using a tuning parameters λ R m in the function of the objective and input parameters υ R l in the set of constraint. The inverse parametric nonlinear cost transportation problem P ( λ , υ ) , where the tuning parameters λ R m in the objective function are tuned (adapted) as less as possible so that the specific feasible solution x has been became the optimal ones for a certain values of υ R l , then, a solution stability set of the parameters was investigated to keep the new optimal solution x unchanged. The proposed method consists of three phases. Firstly, based on the optimality conditions, the parameter λ R m are tuned as less as possible so that the initial feasible solution x has been became new optimal solution. Secondly, using input parameters υ R l resulting problem is reformulated in parametric form P ( υ ) . Finally, based on the stability notions, the availability domain of the input parameters was detected to keep its optimal solution unchanged. Finally, to clarify the effectiveness of the proposed algorithm not only for the inverse transportation problems but also, for the nonlinear programming problems; numerical examples treating the inverse nonlinear programming problem and the inverse transportation problem of minimizing the nonlinear cost functions are presented.

1. Introduction

In the last 20 years, the community of optimization has shown a significant interest in the field of inverse optimization problems. Examples of implementation of inverse optimization in real life have been investigated in various fields such as: traffic, Geophysics, monotonic regression, portfolio and so forth. In any optimization scenario, there are some parameters in the function of the objective and in the functions of the set constraints. When dealing with this problem, generally it is assumed that all these parameters are known but in real situations, there are many cases where the values of the parameters are not precisely known but we may have some fuzzy estimates of the values of these parameters and also have an optimal solution from the past experience or past practice. In these situations, inverse optimization problem can be implemented to adapt the values of the parameters values as less as possible so that the specific solution becomes the optimal one.
Recently, inverse optimization problem is an important new area of investigation involving study and research for the mathematics community [1,2,3,4,5]. In optimization sciences, the coefficients of the cost function are not precisely exactly known; it is truly acceptable to assign a feasible solution x for minimizing objective function nearly optimal, if there exist some nearby cost vector, such that the feasible solution x is an optimal solution for minimizing the objective function. Suppose that the formulation of an optimization problem model consists of a set of optimization parameters (e.g., cost, time,…, etc.), denoted by λ R m , so we can call these parameters as a tuning parameter, which needs to be tuned so as to x become an optimal solution. Some researcher investigated some applications of inverse optimization problems such as inverse of minimum spanning tree problems and inverse of shortest path problem [5,6,7,8]. Several papers have appeared recently studying some inverse combinatorial problems, where a desired given solution x which is feasible and not an optimal for the model and it is needed to adapt the cost coefficients as little as possible, so that the given solution x becomes a new optimal solution of the model. Some application of the inverse optimization problems is discussed in References [5,6]. Ahuja and Orlin [9] present some studies in the field of inverse optimization problems; also they investigated many applications in network analysis. Zhang and Liu [10] present various inverse linear programming models and they also further investigated it in Reference [11].
Huang and Liu [12] investigated some applications of the inverse of linear programming. Some various applications of the inverse problem have been studied by Amin and Emrouznejad [13]. Inverse linear programming and inverse quadratic programming using perturbation methods was investigated by Zhang et al. [14,15]. Scheafer [16] and Wang [17] discussed the inverse of the integer programming problem. Some basic notions in the convex parametric programming with its qualitative analysis was presented in References [18,19,20], also they discussed the solvability set and stability set.
Transportation problems are considered as special kinds of optimization problems. They present real-world activities that are managed with logistics. it includes transportation with manufacturing products in several supplies to several destinations. The aim is to achieve the minimum total transportation cost that will satisfy the demands at various destinations [21]. Furthermore, a few researchers have been studied variety of inverse transportation problem due to their practical motivation. Implementation of inverse optimization have been applied in various fields [2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20,21,22,23,24], where it is used to measure operational variance in transit operators, to detect shifts in travel/traffic objectives in system security and risk management, to learn mechanism in autonomous vehicles and so forth. Andrew et al. [24] present a systematic method to derive obscure costs from observations with noisy data of the optimal transportation plans. They implement a formulation of the problem based on graph theory, where nodes represent countries of graphs and assign nonzero weight on the edges between adjacent countries which have a common border. Thai [25] investigate the implementation of inverse optimization to define two types of problems in transportation; he described how inverse optimization and robust optimization can be implemented to find actual time of travel, with noisy information data on travel times. Sanjay et al. [23] investigate the capacitated transportation problem and its inverse; in this problem the transportation cost unit of some products in the basic problem are adapted as less as possible so that the specific feasible solution becomes an optimal solution. Sanjay [26] presents the transportation problem and its inverse, where he investigates the optimizing the ratio of linear objective subject to the set of linear equality constraints and non-negative constraints. Dequan et al. [27] give a brief overview of the inverse optimization problems of the general linear programming (LP). Xu et al. [22] presents a new inverse optimization models and supporting algorithms to learn the parameters of heterogeneous travelers’ route behavior to infer shared network state parameters he proves that the method can obtain unique dual prices for a network shared by these agents in polynomial time.
The main goal of this proposal is to study the stability of solutions for parametric inverse convex nonlinear programming. This work is formal, an advanced extension of our work [28], where the inverse model for nonlinear programming problem NLP are investigated, also different norms L1, L2 and L∞ were implemented in the solution process. The proposed method consists of three phases. Firstly, based on the optimality conditions, tuning parameter λ R m is adjusted as less as possible so that the specific initial solution x becomes the new optimal solution. Secondly, with the help of input parameters, υ R l the resulting problem is reformulated in parametric form P ( υ ) . Finally, based on the stability notions, the availability domain of the input parameters was detected to keep its optimal solution unchanged. The proposal is structured as follows. Section 2 introduces problem formulation of the nonlinear transportation. Section 3 discusses the inverse of nonlinear programming problem (NLP). Section 4 investigated the stability notions. Section 5 summarizes the solution procedure. Numerical examples of the inverse optimization problems are described in Section 6. Finally, the results are concluded in the last section.

2. Problem Formulation

Nonlinear transportation problem (NTP) is a special case of nonlinear programming problem (NLP), the formulation of the NTP is more specific, especially in terms of decision variables and the set of constraints. The goal of the nonlinear transportation problem is to optimize (minimize) the vector of nonlinear transportation cost function, in addition to meeting demand, supply and transporting constraints. The standard formulation of the optimization model for the parametric nonlinear transportation problem is given as follows:
M i n i m i z e   C T = r = 1 m j = 1 n f r j ( λ r j x r j ) , S u b j e c t   t o :   j = 1 n x r j = v i , r = 1 , 2 , , m   r = 1 m x r j = b j , j = 1 , 2 , , n   x r j 0 r , j
where, objective function variable CT exemplifies the total parametric transportation cost for the transportation problem of a single commodity from sources m to destinations n , f r j ( λ r j x r j ) denotes the cost function of transporting flow x r j which represent the flow from source i to destination j ,   v r ,   b j are parametric capacities of each source i and each destination j respectively, λ i j R m + n are set of parameters at the cost function m is the number of sources and n is the number of destination. Problem (1) could be written as a general parametric nonlinear programming problem, which reformulated as follows:
P ( λ , υ ) : min   f ( x , λ ) s u b j e c t   t o :   x M = { x R n | g i ( x , υ ) 0 , i = 1 , 2 , k }
where, the objective function f ( x , λ ) = r = 1 m j = 1 n f r j ( λ r j x r j ) , is a convex on M ,   λ R m , is a vector of tuning parameters and j = 1 n x r j = a i , r = 1 , 2 , , m ,   r = 1 m x r j = b j , j = 1 , 2 , , n   , x r j 0   r , j is g i ( x , υ ) , i = 1 , 2 , k such that g i ( x , υ ) , i = 1 , 2 , k is a finite-dimensional convex function of class c 1 depending on the input parameter υ R l . For a predetermined υ * (assigned by the decision maker) the problem P ( λ , υ * ) is transformed to P ( λ ) . Also, it is assumed that M is a non-empty set, such that, the feasibility is guaranteed. Suppose that f * the optimum value for P ( λ ) corresponding to the decision variables x * and tuning parameter λ * . To give the precise formulation of the problem, we present the following definitions.
Definition 1.
Solvability set S ( λ ) of P ( λ ) is defined as:
S ( λ ) = { λ R m | min x M   f ( x , λ )   e x i s t s } .
Let x is the initial solution (feasible) for the problem P ( λ ) , knowing that x is not an optimal solution of P ( λ ) problem. In inverse programming, we need to adjust as less as possible the tuning parameter λ so that the predetermined initial feasible solution x (which is not an optimal solution for the original model) becomes a new optimal solution of P ( λ ) , where λ is the new value of tuning parameters. First, let us define the expected domain of λ which is denoted by D ( λ ) , as follows:
Definition 2.
Suppose that λ S ( λ ) with the corresponding optimal solution x ο , then we have D ( λ ) which defined by
D ( λ ο ) = { λ ο S ( λ ο ) |     f ( x ο , λ ο ) = min x M   f ( x , λ ο ) } .
Theorem 1.
D ( λ ) is nonempty for P ( λ ) the problem.
Proof. 
Since it is assumed that the feasibility is guaranteed for P ( λ ) the problem, that is, { λ R m | min x M f ( x , λ )   e x i s t } then, D ( λ ) is not empty. □
Theorem 2
([29]). Let the solution ( x , λ ) be an initial solution (feasible) of the problem P ( λ ) , knowing that I = { i : g i ( x ) = 0 } , Also, suppose that the two functions f ( . ) and g i ( . ) for i I are both differentiable at ( x   , λ   ) and that the function g i is continuous for i I at ( x , λ ) , then if ( x , λ ) is a local solution for the problem P ( λ ) , then there exist scalars u 0 , u i such that,
u ο f ( x ο , λ ) + i I u i g i ( x ο , λ ) = 0 ,     u ο , u i 0 ,         f o r             i I ( u ο , u l ) ( 0 , 0 )
where u i for i I are the components of the vector u ¯ l , furthermore, if g i ( . ) for i I are also differentiable at ( x , λ ) for, i I , then the forgoing conditions can be formulated in the following form:
u ο f ( x ο , λ ) + i I u i g i ( x ο ) = 0 , u i g i ( x 0 ) = o         f o r         i = 1 , 2 , , k u ο , u i 0 ,           f o r         i = 1 , 2 , , k . (   u ο , u l ) ( 0 , 0 )
where, u i   f o r   i = 1 , 2 , , k are the components of vector u .
Proof. 
Since ( x , λ ) solve P ( λ ) locally, then there exists no vector d ¯ . Such that f ( x , λ ) t d ¯ < 0   a n d   g i ( x ) t < 0 for each i I . Now, let A is a matrix whose rows are f ( x , λ ) t d ¯   a n d   g i ( x ) t i I . By Gordon’s theorem [22], there exists a non-zero vector p ¯ 0 such that A t p ¯ = 0 . Denoting the component of p ¯ , u 0 ,   a n d   u i for i I , the first part of the result follows, the equivalent form of the necessary conditions is readily obtained by letting u i = 0 for i I . □
To determine D ( λ ) of P ( λ ) problem, the following theorems are presented.
Theorem 3
([22]). For a certain υ , let ( x , λ ) be an initial feasible solution of P ( λ ) , knowing that I = { i : g i ( x ) = 0 } ., Suppose that f ( . )   a n d   g i ( . ) for i I are differentiable at ( x , λ ) and that g i ( . ) for i I is continuous at ( x , λ ) . Furthermore, suppose that g i ( x ) for i I are linearly independent. If ( x , λ ) solves P ( λ ) locally, then there exist scalars u i for i I such that
f ( x ο , λ ) + i I u i g i f ( x ο , λ ) = 0 , u i 0 ,           i I
If g i ( . ) , i I are also differentiable at ( x , λ ) , then the forgoing conditions can be reformulated in the following equivalent form:
f ( x ο , λ ) + i I u i g i ( x ο , λ ) = 0 , u i g i ( x 0 ) = o f o r         , i = 1 , 2 , , k u i 0 f o r         , i = 1 , 2 , , k .
Proof. 
By Theorem 2, there exist scalars u ο , u i for i I , not all equal to zero, such that
u ο f ( x ο , λ ) + i I u i g i ( x ο , λ ) = 0 , u ο , u i 0 ,       f o r     i I
Knowing that, u ο > 0 , because the system (8) would be contradict with the assumption of linear independence of g i ( x , λ ) for i I if   u ο = 0 . The first part of the theorem then follows by letting u i = u i u o for each i I . The equivalent form of system (6) follows by letting u i = 0 for i I . □
Now, let the system as follows:
f ( x , λ ο ) x α + i = 1 k u i g r ( x ) x α = 0 , α = 1 , 2 , n ,   λ ο R m u i g i ( x ) = o , i = 1 , 2 , , k u i 0 , i = 1 , 2 , , k .
From system (9), we can determine D ( λ ) . Since, this system (9) represents independence n + k equations in k + m unknowns and λ R m , which are linear in u i , i = 1 , 2 , , k and nonlinear in λ , then we can obtain u i , i = 1 , 2 , , k and λ explicitly after substituting x by x ο .

3. Inverse Nonlinear Programming Problem

The purpose of the proposed inverse optimization approach is to adapt the tuning parameter value from λ *   to λ so that the specific feasible (given) solution x becomes an optimal solution of P ( λ ) by dealing with the problem
I P ( λ ) : min   λ λ * s u b j e c t   t o :   D ( λ ) = { λ S :     f ( x ο , λ ) = min x M   f ( x , λ ) }
Let denote a vector norm, such as L1, L2 or L∞ norm and N is a domain of λ - parameters that satisfy { λ S ( λ ) :     f ( x ο , λ ) = min x M   f ( x , λ ) } , it is clear that the domain of the inverse nonlinear programming problem is the same domain of D ( λ ) , so that problem (11) can be reformulated as follows:
Q ( λ ) : min   λ ο λ * s u b j e c t   t o :   λ ο D ( λ ο )
Definition 3.
The set of optimality solution of problem (12) which represented by N o p t ( λ , x ο ) is defined as follows: N o p t ( λ , x ο ) = { λ R m | min λ D ( λ ) λ λ * e x i s t s } .
Note that: N o p t N
Theorem 4.
The set of optimality solutions N o p t ( λ , x ο ) is closed and convex
Proof. 
  • First: Convexity:
Let, λ 1 N o p t ( λ , x ο ) then
( 1 w ) min λ D ( λ ) λ λ * = ( 1 w ) λ 1 λ *   , w 1 ,
and let λ 2 N o p t ( λ , x ο ) , then
w min λ D ( λ ) λ λ * = w λ 2 λ *   , w 0
Adding (12), (13) we get,
w λ 2 λ * + ( 1 w ) λ 1 λ * = min λ D ( λ ) λ λ * , 0 w 1
Then form (14), we have,
( 1 w ) λ 1 o + w λ 2 ο λ *   = min λ D ( λ ο ) λ ο λ * , 0 w 1
Then from (15), we can say that, ( 1 w ) λ 1 + w λ 2 N o p t ( λ , x ο )
Hence, N o p t ( λ , x ο ) is convex.
  • Second: Closeness:
Let { λ ¯ n } N o p t ( λ , x ο ) be a sequence of points λ n which converges to λ ¯ , that is, lim n λ ¯ n = λ ¯
Sin ce ,   min λ D ( λ ) λ λ * = λ n λ * , n = 1 , 2 , ,   then ,   lim n min λ D ( λ ) λ λ * = lim n λ n λ *
From the continuity of the norm, it follows that
lim n λ n λ * = λ ¯ λ * ,   then   λ ¯ N o p t ( λ , x ο ) = min λ D ( λ ) λ λ *
Then, N o p t ( λ , x ο ) is closed. □

4. The Stability of the Optimal Solution in the Decision Space

Studying the stability to the problem involves finding the value range of the input parameter to keep the optimal solution unchanged (the optimal solution in the decision space remains unchanged), which denoted as the stability of the first type [27].
For any, λ * R m the problem P ( λ * , υ ) is transformed to P ( υ ) , which have only input parameter.
P ( υ ) : M i n m i z e   f ( x ) s . t .   x M = { x R n / g i ( x , υ ) 0 , i = 1 , 2 , , k } .
Definition 4.
The set of parameters υ R l (input parameters) for the problem P ( υ ) is defined by
U = { υ R l | M ϕ }  
Theorem 5.
The set U is convex.
Proof. 
If U is the whole space R l , then the proof is clear. Otherwise, suppose that υ 1 , υ 2 U , then there exist x 1 , x 2 R n such that g i ( x 1 , υ 1 ) 0     ,   g i ( x 2 , υ 2 ) 0 and for 0 w 1 , ( 1 w ) g i ( x 1 , υ 1 ) 0   ,   w g i ( x 2 , υ 2 ) 0 . Since   g r ( x , υ ) are convex in x , υ , then, we have,
g i [ ( 1 w ) x 1 + w x 2 , ( 1 w ) υ 1 + w υ 2 ]   ( 1 w ) g i ( x 1 , υ 1 ) + w g i ( x 2 , υ 2 )
and hence, ( 1 w ) υ 1 + w υ 2 U . Then U is convex. □
Definition 5.
The set V is denoted as the solvability set of P ( υ ) and it is defined as
V = { υ R l | m o p t ( υ ) ϕ } .
where, m o p t ( υ ) = { x * R n | f ( x * ) = min x M f ( x ) } . Note that m o p t M .
Definition 6.
Let υ ¯ V with a corresponding optimal solution x ο , then the stability set of the first type of P ( υ ) corresponding to the solution x ο which is defined as follows:
G ( x o ) = { υ R m | f ( x o ) = min x M ( υ ) f ( x ) }
If f ( x ) and g i ( x , υ )   ,   i = 1 , 2 , , k are set of functions belong to class C 1 on R n . Let υ V with a corresponding optimal point x ο ; then from the stability of problem Q ( λ ) , there exists u i , i = 1 , 2 , , k such that ( x ο , υ ) solves the Kuhn- Tucker condition problem [22], which can be described as follows:
f ( x o ) x j + r = 1 k g i ( x o , υ ¯ ) x j = 0   ,   i = 1 , 2 , , k ,   j = 1 , 2 , , n u i ¯ g i ( x o , υ ¯ ) = 0 , , i = 1 , 2 , , k g i ( x o , υ ¯ ) 0 , i = 1 , 2 , , k u i ¯ 0 , i = 1 , 2 , , k
To determine the set G ( x ο ) , let us consider the system:
f ( x o ) x j + r = 1 k u i g i ( x o , υ ) x j = 0           ,     i = 1 , 2 , , k ,     j = 1 , 2 , , n
Which represents n equations in k + m unknowns, which are linear in u and nonlinear in υ , we denote to this system by φ ( x ο , u , υ ) = 0 .
If m = n ,   φ ( x ο , u , υ ) = 0 has a solution, υ φ ( x ο , u , υ ) exists and is a continuous function of ( u , υ ) and [ φ ( x ο , u , υ ) ] 1 exists. The solution of φ ( x ο , u , υ ) = 0 can be expressed explicitly as υ = ψ ( u ) where ψ is n -dimensional vector function. The value of u i ,   i = 1 , 2 , , k in such a way that ( u r , υ ) solves Kuhn-Tucker problem where ( u , υ ) solves the system
f ( x o ) x j + r = 1 k u i g i ( x o , υ ) x j = 0         ,     i = 1 , 2 , , k ,     j = 1 , 2 , , n
The following cases were considered [20]
(i)  u i > 0 , i J ( 1 , 2 , , k ) , u i = 0 ,   i ( 1 , 2 , , k ) J .
We define the set,
G J ( x ο ) = { υ R l / g i ( x ο , υ ) = 0 , g i ( x ο , υ ) 0 , i j }
And we define,
G 1 ( x ο ) = J q G J ( x ο )
where q is a proper subset of J
(ii)  u i = 0     , i = 1 , 2 , k
We define the set,
G 2 ( x o ) = { v R l / g i ( x o , υ ) 0 , i = 1 , 2 , , k }
(iii)  u i > 0 , i = 1 , 2 , k
We define the set
G 3 ( x ο ) = { v D R m / g i ( x ο , υ ) = 0 , i = 1 , 2 , , k }
From Kuhn-Tucker sufficient optimality [22], it follows that the sets G i ( x ο ) ,   i = 1 , 2 , 3 or the union of some or all of them depending on the values of u .

5. Solution Procedure

The following are the main steps of our method that are used to find the inverse nonlinear programming problem and to investigate the stability of the solution in the decision space. The main steps can be stated as follows:
Phase 1:
Obtain the value of the tuning parameters λ R m , so that the given (determined) feasible solution x ο becomes the optimal ones.
Step 0.
For certain input parameters υ R l , the problem P ( λ , υ ) is transformed to P ( λ )
Step 1.
Obtain the optimal solution x * R n and the corresponding optimum value f * for the problem P ( λ ) for a certain λ * parameter.
Step 2.
Choose the desired feasible decision variables x ο M which determined by the decision maker (DM).
Step 3.
Obtain D ( λ ) condition of P ( λ ) .
Step 4.
Formulate the problem Q ( λ ) .
Step 5.
Solve the problem Q ( λ ) to obtain the vector λ 0 with three different main definitions of the norm (as L1, L2 or L∞ norm).
Step 6.
Formulate the inverse nonlinear programming problem { P ( λ ) |   λ = λ 0 } .
Phase 2:
Formulate parametric inverse nonlinear programming problem
With the help of input parameters, υ R l the resulting problem is reformulated in parametric form P ( υ ) = { P ( λ , υ ) | λ = λ ο } as in Equation (16).
Phase 3:
Stability analysis
Based on the stability notions, the availability domain of the input parameters υ R l was found to keep its optimal solution unchanged.
Step 1.
Formulate parametric problem P ( υ ) = { P ( λ , υ ) | λ = λ } .
Step 2.
Construct the KKT as in Equation (17).
Step 3.
Determine the values of Lagrange multipliers.
Step 4.
Determination of the availability domain of the input parameters υ R l , according to Equations (18)–(20).

6. Numerical Simulation

To validate our method, three inverse parametric nonlinear programming examples are given, having tuning parameters λ R m at the objective functions and input parameters υ R l in the constraint and a transportation application are presented.

6.1. Classical Benchmark Examples

To examine the proposed inverse optimization method, three examples were chosen from the literature.
Example 1.
Given the nonlinear programming problem P ( λ , υ ) having tuning parameters λ R 2 in the function in the objective and input parameters υ R 2 in the function of the constraint,
min f ( x , λ ) = x 1 2 + x 2 2 λ 1 x 1 + λ 2 x 2 s . t . :   υ 1 x 1 + x 2 υ 2 2 0 , υ 1 2 x 1 2 + υ 2 x 2 2   υ 1 0 , x 1 , x 2 0 .
Step 0.
For certain input parameters υ = ( 1 , 1 ) , P ( λ , υ ) is transformed to P ( λ )
min   f ( x , λ ) = x 1 2 + x 2 2 λ 1 x 1 + λ 2 x 2 , s . t . :   x 1 + x 2 1 , x 1 2 + x 2 2 1 , x 1 , x 2 0 .
Step 1.
The optimal solution x * of the problem P ( λ ) is found x * = ( 1 , 0 ) at f * = 1 and λ * = ( 2 , 0 ) ,
Step 2.
The desired feasible decision variables are x = ( 0 , 1 )
Step 3.
Obtain D ( λ ) condition of NLP as follows:
2 x 1 + λ 1 ο + u 1 + 2 x 1 u 1 = 0   , 2 x 2 + λ 2 ο + u 1 + 2 x 2 u 2 = 0   , u 1 ( x 1 + x 2 1 ) = 0   , u 2 ( x 1 2 + x 2 2 1 ) = 0 , x 1 , x 2 , u 1 , u 2 0 .
Substituting with the desired feasible decision variables x ο = ( 0 , 1 ) we get,
λ 1 ο + u 1 = 0 , 2 + λ 2 ο + u 1 + 2 u 2 = 0 , u 1 ( 1 1 ) = 0   , u 2 ( 1 1 ) = 0 , u 1 , u 2 0 .
From system (21), we get D ( λ ) as follows
D ( λ ) = { ( λ 1 , λ 2 ) R 2 | λ 1 0 , λ 2 λ 1 2 }
Step 4.
Formulate the Q ( λ ) problem
min λ λ * , s u b j e c t   t o :   λ D ( λ ) .
Substituting λ * = ( 2 , 0 ) , f * = 1 , D ( λ ) = { ( λ 1 , λ 2 ) R 2 | λ 1 0 , λ 2 λ 1 2 }
Step 5.
Using L 2 -norm to solve the problem Q ( λ ) ) to obtain the vector λ 0 as follows:
min   ( λ 1 + 2 ) 2 ( λ 2 ) 2 , s u b j e c t   t o :       λ 1 0 ,       2 + λ 2 λ 1 0 ,       1 + λ 2 = 1 .
Solving this NLP using “Lindo” software we get λ 1 = 0 ,   λ 2 = 2 then we get λ ο = ( 0 , 2 )
Phase 2:
Formulate parametric inverse nonlinear programming problem
Substituting by λ 1 = 0 ,   λ 2 = 2 the P ( υ ) takes the form
min   f ( x , λ ) = x 1 2 + x 2 2 2 x 2 , s . t . :   υ 1 x 1 + x 2 υ 2 2 0 , υ 1 2 x 1 2 + υ 2 x 2 2   υ 1 0 , x 1 , x 2 0 .
Phase 3:
Stability analysis
Construct the KKT. conditions for P ( υ ) with the optimal solution x ο = ( 0 , 1 )
u 1 υ 1 = 0 , 2 2 + u 1 + 2 u 2 υ 2 = 0   ,   1 υ 2 2 0 ,   υ 2 υ 1 0 , u 1 ( + 1 υ 2 2 ) = 0 , u 2 ( υ 2 υ 1 ) = 0 , u 1 , u 2 0 .
Then we get the stability set of the first type as follows:
G ( 0 , 1 ) = { υ R 2 / υ 1 υ 2 0 , 1 υ 2 2 = 0 } { υ R 2 / 1 υ 2 2 0 , υ 1 υ 2 = 0 }
Example 2.
Consider the nonlinear programming problem P ( λ , ν ) having tuning parameters λ R 2 in the function of the objective and input parameters υ R 2 in the functions of the constraint,
min f ( x , λ ) = 4 x 1 + λ 1 x 2 + λ 2 x 1 2 2 x 1 x 2 2 x 2 2 , s . t . :   υ 1 x 1 + υ 2 x 2 1 0 , x 1 , x 2 0 .
Step 0.
For certain input parameters υ = ( 1 , 2 ) the problem P ( λ , υ ) is transformed to P ( λ )
min   f ( x , λ ) = 4 x 1 + λ 1 x 2 + λ 2 x 1 2 2 x 1 x 2 2 x 2 2 , s . t . :   x 1 + 2 x 2 2 , x 1 , x 2 0 .
Step 1.
The optimal solution x * of P ( λ ) is found x * = ( 1 3 , 5 6 ) with z * = 4.16 , for λ * = ( 6 , 2 ) ,
Step 2.
The desired feasible decision variables are x ο = ( 3 2 , 1 4 )
Step 3.
Obtain D ( λ ) condition of P ( λ ) as follows:
4 2 λ 0 2 x 1 + 2 x 2 + u = 0 λ 0 1 + 2 x 1 + 4 x 2 + 2 u = 0 u ( x 1 + 2 x 2 2 ) = 0 u 0 .
Substituting x 0 = ( x 1 , x 2 ) = ( 3 2 , 1 4 ) we get,
4 3 λ 2 + 1 2 + u = 0 , λ 1 + 3 + 1 + 2 u = 0 , u ( 3 2 + 1 2 2 ) = 0 , u 0 .
From system (22) the D ( λ ) is as follows:
D ( λ ) = { ( λ 1 , λ 2 ) R 2 / λ 1 4 , λ 2 7 6 }
Step 4.
The desired feasible decision variables are x 0 = ( 3 2 , 1 4 )
Step 5.
Formulate the Q ( λ ) problem as follows
min λ λ * , s u b j e c t   t o : λ D ( λ ) .
Substituting λ * = ( 6 , 2 ) f * = 4.16 , D ( λ ) = { ( λ 1 , λ 2 ) R 2 / λ 1 4 , λ 2 7 6 } .
Step 6.
Using L 2 - norm to solve the problem Q ( λ ) to obtain the vector λ 0 as follows:
min   ( λ 1 6 ) 2 + ( λ 2 + 2 ) , s u b j e c t   t o :   λ 1 4 , λ 2 7 6 .
Solving this NLP using “Lindo” software we get λ 1 = 5.066 ,   λ 2 = 0.989 so we get λ 0 = ( 5.066 , 0.989 )
Phase 2:
Formulate parametric inverse nonlinear programming problem
Substituting by λ 1 0 = 5.066 ,   λ 2 0 = 0.989 the P ( υ ) takes the form:
min   f ( x , λ ) = 4 x 1 + 5.006 x 2 0989 x 1 2 2 x 2 2 2 x 1 x 2 , s . t . :     υ 1 x 1 + υ 2 x 2 1 0 , x 1 , x 2 0 .
Phase 3:
Stability analysis
Construct the KKT. Conditions for P ( υ ) with the optimal solution x 0 = ( 3 2 , 1 4 )
4 + 2.967 + 1 2 + u υ 1 = 0 , 5.006 + 3 + 1 + 2 u υ 2 = 0 , 3 2 υ 1 + 1 4 υ 2 2 0 , u ( 3 2 υ 1 + 1 2 υ 2 2 ) = 0 , u 0 .
Then we get the stability set of the first type as follows:
G ( 3 2 , 1 4 ) = { υ R 1 / 3 2 υ 1 + 1 2 υ 2 2 = 0 }
Example 3.
Given the nonlinear programming problem P ( λ , υ ) having tuning parameters λ R 2 in the objective functions and input parameters υ R 2 in the constraint,
min   f ( x , λ ) = λ 1 x 1 2 + λ 2 x 2 , s u b j e c t   t o : υ 1 x 1 + x 2 2 υ 2 2 o , x 1 + υ 1 2 x 2 + υ 2 0 x 1 , x 2 0 .
Step 0.
For certain input parameters υ = ( 1 , 0 ) the problem P ( λ , υ ) is transformed to P ( λ )
Step 1.
The optimal solution x * of P ( λ ) the problem is   x * = ( 0 , 0 ) with f * = 0 , for λ * = ( 2 , 0 ) ,
Step 2.
The desired feasible decision variables are x * = ( 0.5 , 0.5 )
Step 3.
Obtain D ( λ ) condition of NLP as follows:
2 x 1 λ 1 ο + u 1 u 2 = 0 , λ 2 ο + 2 u 1 x 2 + u 2 = 0 , u 1 ( x 1 + x 2 2 ) = 0 , u 2 ( x + x 2 ) = 0 , x 1 , x 2 , u 1 , u 2 0 ,
Substituting with the desired feasible decision variables x ο = ( 0.5 , 0.5 ) we get,
λ 1 ο + u 1 u 2 = 0 , λ 2 ο + u 1 + u 2 = 0 , u 1 ( 1 2 + 1 4 ) = 0 , u 2 ( 1 2 1 2 ) = 0 , u 1 , u 2 0 ,
From system (23), we get D ( λ ) is as follows
D ( λ ) = { ( λ 1 , λ 2 ) R 2 / λ 1 0 , λ 2 0 , λ 1 + λ 2 = 0 }
Step 4.
Formulate the Q ( λ ) problem
min λ λ * s u b j e c t   t o :   λ D ( λ )
Substituting λ * = ( 2 , 0 ) , f * = 0 , D ( λ ) = { ( λ 1 , λ 1 ) = R 2 / λ 1 0 , λ 2 0 , λ 1 + λ 2 = 0 }
Step 5
Using L 2 -norm to solve the problem I P ( λ ) ) to obtain the vector λ 0 as follows:
min   ( λ 1 + 2 ) 2 + λ 2 , s u b j e c t   t o :           λ 1 + λ 2 = 0 , λ 1 0 , λ 2 0 .
Solving this NLP using “LINDO” software we get λ 1 = 1 ,   λ 2 = 1 so we get λ 1 ο = ( 1 , 1 )
Phase 2:
Formulate parametric inverse nonlinear programming problem
Substituting by ( λ 1 , λ 2 ) = ( 1 , 1 ) , letting us form the parametric P ( υ ) problem: parameters in the constraints as follows:
min   f ( x , λ ) = x 1 2 x 2 , s u b j e c t   t o : υ 1 x 1 + x 2 2 υ 2 2 o , x 1 + υ 1 2 x 2 + υ 2 0     x 1 , x 2 0 .
Phase 3:
Stability analysis
Construct the KKT. Conditions for P ( υ ) with the optimal solution   x ο = ( 1 2 , 1 2 )
1 + u 1 υ 1 u 2 = 0 , 1 + u 1 + u 2 υ 1 2 = 0 , u 1 ( 1 2 υ 1 + 1 4 υ 2 2 ) = 0 , u 2 ( 1 2 + 1 2 υ 1 2 + υ 2 ) = 0 , u 1 , u 2 0 ,
Then we get the stability set of the first type as follows:
G ( 1 2 , 1 2 ) = { υ R 2 / 1 2 υ 1 υ 2 2   1 4 ,   1 2 υ 1 2 + υ 2 = 1 2 } { υ R 2 /   1 2 υ 1 2 + υ 2 1 2 , 1 2 υ 1 υ 2 2 =   1 4 }
On solving the previous examples by the given approach, we stress the following:
At the first example, when the input parameters υ = ( 1 , 1 ) and the tuning parameters λ * = ( 2 , 0 ) , the optimal solution was x * = ( 1 , 0 ) but the decision making need to have the point x = ( 0 , 1 ) as an optimal one, this methodology not only achieve the goal but also detect the stability set of the first type, that is used to define available range of these input parameters G ( 0 , 1 ) = { υ R 2 / υ 1 υ 2 0 ,   1 υ 2 2 = 0 } { υ R 2 /   1 υ 2 2 0 , υ 1 υ 2 = 0 } that keep the predetermined point x * = ( 1 , 0 ) as optimal solution.
At the second one, when the input parameters υ = ( 1 , 2 ) and the tuning parameters λ * = ( 6 , 2 ) , then the corresponding optimal solution was x * = ( 1 3 , 5 6 ) but the decision making wish that the point x ο = ( 3 2 , 1 4 ) to be an optimal one, by using this method the goal was achieved in addition, the stability set of the first type was detected to define available range of these input parameters G ( 3 2 , 1 4 ) = { υ R 1 / 3 2 υ 1 + 1 2 υ 2 2 = 0 } that keep the point x ο = ( 3 2 , 1 4 ) as an optimal one.

6.2. An Application: Transportation Problem Application

To examine the proposed inverse optimization method, nonlinear transportation problem was chosen. When the unit of the transportation cost on a specific road is nonlinear depending on the number of the transported units, then the transportation problem is called nonlinear transportation problem (NTP). Investigating for the optimal solution of NTP has been one of the important fields of intensive study on supply chain process. This section proposed an algorithm for inverse transportation problem of minimizing maximizing the nonlinear cost functions. The approach can be used to obtain the modified values of objective coefficients such that the specific (given) feasible solution becomes an optimal solution. A transportation network model shown in Figure 1 has two factories, factory 1 and factory 2 which represent the source nodes; on the other hand, the destination nodes represent warehouse 1, warehouse 2 and warehouse 3, any existing by an arc. The cost of each unit product unit, through specific path is represented by the numbers on that arc. Cost, supply and demand values are presented in Table 1.
P ( λ , υ ) : M i n   F = λ 1 x 11 2 + 2 x 12 2 + 7 x 13 2 + λ 2 x 21 2 + 9 x 22 2 + 3 x 23 2 s . t . :           x 11 + x 12 + x 13 υ 1 , x 21 + x 22 + x 23 υ 2 , x 11 + x 21 = 11 , x 12 + x 22 = 3 , x 13 + x 23 = 14 , x i j 0
Step (0).
For certain υ = ( 13 , 15 ) the problem P ( λ , υ ) is transformed to P ( λ )
P ( λ ) M i n   F = λ 1 x 11 2 + 2 x 12 2 + 7 x 13 2 + λ 2 x 21 2 + 9 x 22 2 + 3 x 23 2 s . t . :           x 11 + x 12 + x 13 13 , x 21 + x 22 + x 23 15 , x 11 + x 21 = 11 , x 12 + x 22 = 3 , x 13 + x 23 = 14 , x i j 0
Step (1).
the optimal solution is
x 11 = 11 ,   x 12 = 0 ,   x 21 = 0 ,   x 22 = 3 ,   x 13 = 2 ,   x 23 = 12 ,   F * = 622   With   λ * = ( 1 , 1 )
Step (2).
the desired feasible solution is
x = ( x 11 ,   x 12 , x 21 , x 22 ,   x 13 , x 23 ) x 11 = 11 ,   x 12 = 2 ,   x 21 = 0 ,   x 22 = 1 ,   x 13 = 0 ,   x 23 = 14
Step (3).
Formulate the K.K.T conditions to get the domain of λ such that x is optimal one
2 λ 1 x 11 + u 1 + v 1 = 0 4 x 11 + u 1 + v 2 = 0 14 x 13 + u 1 + v 3 = 0 2 λ 2 x 21 + u 2 + v 1 = 0 18 x 22 + u 2 + v 2 = 0 9 x 23 + u 2 + v 3 = 0 u 1 ( x 11 + x 12 + x 13 13 ) = 0 u 2 ( x 21 + x 22 + x 23 15 ) = 0 u 1 , u 2 0 ,   v 1 , v 2 , v 3 0
By Substituting with x we get
D ( λ ) = { ( λ 1 , λ 2 ) R 2 | λ 1 10 0 , λ 2 14 0 }
Step (4).
Formulate D ( λ ) to find the of values λ 1 , λ 2
min λ λ * s . t . λ D ( λ )
For L 2 norm
min λ λ * = min ( λ 1 1 ) 2 + ( λ 2 1 ) 2 s . t . λ D ( λ )
Solving this NLP using “LINDO” software we get λ 1 = 10 ,   λ 2 = 14 , then λ = ( 10 , 14 )
Phase 2:
Formulate parametric inverse nonlinear transportation problem
Substituting by λ = ( 10 , 14 ) , letting us form the parametric P(u) problem: parameters in the constraints as follows:
P ( υ ) : M i n   F = 10 x 11 2 + 2 x 12 2 + 7 x 13 2 + 14 x 21 2 + 9 x 22 2 + 3 x 23 2 s . t . :           x 11 + x 12 + x 13 υ 1 , x 21 + x 22 + x 23 υ 2 , x 11 + x 21 = 11 , x 12 + x 22 = 3 , x 13 + x 23 = 14 , x i j 0
Phase 3:
Stability analysis
Construct the K.K.T of P ( υ ) with optimal solution
x = ( x 11 ,   x 12 , x 21 , x 22 ,   x 13 , x 23 ) x 11 = 11 ,   x 12 = 2 ,   x 21 = 0 ,   x 22 = 1 ,   x 13 = 0 ,   x 23 = 14
and λ = ( 10 , 14 ) to be as follows:
24 x 11 + u 1 + v 1 = 0 4 x 12 + u 1 + v 2 = 0 14 x 13 + u 1 + v 3 = 0 28 x 21 + u 2 + v 1 = 0 18 x 22 + u 2 + v 2 = 0 9 x 23 + u 2 + v 3 = 0 u 1 ( x 11 + x 12 + x 13 v 1 ) = 0 u 2 ( x 21 + x 22 + x 23 v 2 ) = 0 u 1 , u 2 0 v 1 , v 2 , v 3 0
At x ,
264 + u 1 + v 1 = 0 8 + u 1 + v 2 = 0 u 1 + v 3 = 0 u 2 + v 1 = 0 36 + u 2 + v 2 = 0 126 + u 2 + v 3 = 0 u 1 ( 11 + 2 + 0 v 1 ) = 0 u 2 ( 0 + 1 + 14 v 2 ) = 0 u 1 , u 2 0 v 1 , v 2 , v 3 0
It is clear that u 1 , u 2 0 then we get the stability set of the first type at x as follows
G ( 12 , 14 ) = { v R 2 | 13 v 1 = 0 , 12 v 2 = 0 }
It is clear that, for this application, we not only justify the cost function parameters as little as possible so that the specific feasible solution ( x 11 = 11 ,   x 12 = 2 , x 21 = 0 , x 22 = 1 ,   x 13 = 0 , x 23 = 14 ) becomes an optimal one but also, a solution stability set of parameters was investigated to keep the new optimal solution is unchanged.

7. Conclusions

The inverse optimization problem is an interesting field for both academic scientist and real-life applications. Implementation the inverse optimization and adapting the cost function parameters as little as possible so that the specific feasible solution becomes an optimal spatially in nonlinear domain is not easy, so keeping that solution with different sets of certain parameters is valuable. Nonlinear transportation problem (NTP) is a special case of nonlinear programming problem (NLP), the formulation of the NTP is more specific, especially in terms of decision variables and the set of constraints.
This manuscript proposed a methodology for finding the inverse problem of convex nonlinear programming problem P ( λ , υ ) having tuning parameters λ R m in the function of the objective and input parameters υ R l in the functions in the set of constraint. The proposed method consists of three phases. Firstly, based on the optimality conditions, tuning parameters λ R m are tuned as less as possible so that the given initial feasible solution x ο becomes the optimal ones. Secondly, using input parameters, υ R l the resulting problem is reformulated in parametric form P ( υ ) . Then, based on the stability notions, the availability domain of the input parameters was detected to keep its optimal solution unchanged.
Finally, to validate and demonstrate the advantage of the new approach, three nonlinear programming examples and nonlinear transportation problem application are provided for the sake of illustration. On solving the transportation problem by the given approach, we summarize the result as follows:
When the input parameters υ = ( 13 , 15 ) and the tuning parameters λ * = ( 1 , 1 ) the optimal solution was x 11 = 11 ,   x 12 = 0 ,   x 21 = 0 ,   x 22 = 3 ,   x 13 = 2 ,   x 23 = 12 , but the decision making need that the solution x 11 = 11 ,   x 12 = 2 ,   x 21 = 0 ,   x 22 = 1 ,   x 13 = 0 ,   x 23 = 14 to be an optimal one, this methodology not only achieve the goal but also detect the stability set of the first type, which is used to define a available range of these input parameters G ( 12 , 14 ) = { v R 2 | 13 v 1 = 0 , 12 v 2 = 0 } that keep the solution x 11 = 11 ,   x 12 = 2 ,   x 21 = 0 ,   x 22 = 1 ,   x 13 = 0 ,   x 23 = 14 as an optimal one. From the above study, the following may be concluded:
  • A solution of a parametric inverse transportation problem is introduced.
  • The paper deals with Parametric nonlinear programming having tuning parameters λ R m in the objective and input parameters υ R l in the constraints.
  • An inverse model for the proposed problem was investigated.
  • Solution stability of the problem was investigated to retain its optimal solution.
  • Numerical examples are provided for the sake of illustration.
For the future work, this method can be extended to nonlinear Large-Scale Inverse transportation Problems and its applications in IoT.

Author Contributions

Conceptualization, A.A.A.M. and Y.A.-E.; Methodology, A.A.A.M. and Y.A.-E.; Investigation, A.A.A.M. and Y.A.-E.; Resources, A.A.A.M. and Y.A.-E.; Writing-Editing, A.A.A.M. and Y.A.-E. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Acknowledgments

The authors would like to express their sincere gratitude to the Deanship of Scientific Research, in Taif University for funding Taif University Researcher Supporting Project number (TURSP-2020/48), Taif University, Taif, Saudi Arabia.

Conflicts of Interest

The authors declare no conflict of interests.

Abbreviations

λ Tuning parameters
υ Input parameters
x The new optimal solution
P ( λ , υ ) The inverse parametric nonlinear programming problem
P ( υ ) Parametric nonlinear programming problem
CTThe total parametric transportation cost
f i j ( λ i j x i j ) The cost transportation function
x i j Transportation flow
ν i , b j Parametric capacities of each source i and each destination j
S ( λ ) The Solvability set
D ( λ ) The expected domain of λ
λ The new value of tuning parameters
λ * The initial value of tuning parameters
x * The decision variables at tuning parameters λ *
υ * The initial value of the input parameters
V The solvability set of P ( υ )
I P ( λ ) The inverse optimization problem
Q ( λ ) The equivalent inverse optimization problem
N o p t ( λ , x ο ) The set of optimality solution of Q ( λ ) problem
U The set of the input parameters
G i ( x ο ) The set of parameters that kept x ο optimal solution

References

  1. Ghobadi, K.; Lee, T.; Mahmoudzadeh, H.; Terekhov, D. Robust inverse optimization. Oper. Res. Lett. 2018, 46, 339–344. [Google Scholar] [CrossRef] [Green Version]
  2. Babier, A.; Chan, T.C.; Lee, T.; Mahmood, R.; Terekhov, D. A unified framework for model fitting and evaluation in inverse linear optimization. arXiv 2019, arXiv:1804.04576. [Google Scholar]
  3. Chan, T.C.; Kaw, N. Inverse optimization for the recovery of constraint parameters. Eur. J. Oper. Res. 2020, 282, 415–427. [Google Scholar] [CrossRef] [Green Version]
  4. Chow, J.Y.; Ritchie, S.G.; Jeong, K. Nonlinear inverse optimization for parameter estimation of commodity-vehicle-decoupled freight assignment. Transp. Res. Part E Logist. Transp. Rev. 2014, 67, 71–91. [Google Scholar] [CrossRef]
  5. Burton, D.; Toint, P.L. On an instance of the inverse shortest paths problem. Math. Program. 1992, 53, 45–61. [Google Scholar] [CrossRef] [Green Version]
  6. Xu, S.; Zhang, J. An inverse problem of the weighted shortest path problem. Jpn. J. Ind. Appl. Math. 1995, 12, 47–59. [Google Scholar] [CrossRef]
  7. Zhang, J.; Liu, Z.; Ma, Z. On the inverse problem of minimum spanning tree with partition constraints. Math. Methods Oper. Res. 1996, 44, 171–187. [Google Scholar] [CrossRef]
  8. Zhang, J.; Liu, Z. A general model of some inverse optimization problems and its solution method under l8 norm. J. Comb. Optim. 2002, 6, 207–227. [Google Scholar] [CrossRef]
  9. Ahuja, R.K.; Orlin, J.B. Inverse Optimization. Oper. Res. 2001, 49, 771–783. [Google Scholar] [CrossRef]
  10. Zhang, J.; Liu, Z. Calculating some inverse linear programming problems. J. Comput. Appl. Math. 1996, 72, 261–273. [Google Scholar] [CrossRef] [Green Version]
  11. Zhang, J.; Liu, Z. A further study on inverse linear programming problems. J. Comput. Appl. Math. 1999, 106, 345–359. [Google Scholar] [CrossRef]
  12. Huang, S.; Liu, Z. On the inverse problem of linear programming and its application to minimum weight perfect k-matching. Eur. J. Oper. Res. 1999, 112, 421–426. [Google Scholar] [CrossRef]
  13. Amin, G.R.; Emrouznejad, A. Inverse forecasting: A new approach for predictive modeling. Comput. Ind. Eng. 2007, 53, 491–498. [Google Scholar] [CrossRef]
  14. Zhang, J.; Zhang, L.-W.; Xiao, X. A Perturbation approach for an inverse quadratic programming problem. Math. Methods Oper. Res. 2010, 72, 379–404. [Google Scholar] [CrossRef]
  15. Jiang, Y.; Xiao, X.; Zhang, L.; Zhang, J. A perturbation approach for a type of inverse linear programming problems. Int. J. Comput. Math. 2011, 88, 508–516. [Google Scholar] [CrossRef]
  16. Scheafer, A.J. Inverse integer programming. Opt. Lett. 2009, 3, 483–489. [Google Scholar] [CrossRef]
  17. Wang, L. Cutting plane algorithms for the inverse mixed integer linear programming problem. Oper. Res. Lett. 2009, 37, 114–116. [Google Scholar] [CrossRef]
  18. Osman, M.S.A. Qualitative Analysis of Basic Notions in Parametric Convex Programming II: Parameters in the Objective Function. Aplikace Matematiky 1997, 22, 333–348. [Google Scholar] [CrossRef]
  19. Osman, M.S.A.; EL-Hefni, M.R. Determination of the Stability Sets in Parametric Convex Programming Problems. In Proceedings of the Second Conference on Operations Research and its Military Applications, Military Technical College (M.T.C.), Cairo, Egypt, 17–19 November 1987. [Google Scholar]
  20. Osman, M.S.A.; El-Hefny, M.R.; El-Ariny, A.K.H. Implementing stability results in solving large scale convex programming problems. In Analysis and Optimization of Systems. Lecture Notes in Control and Information Sciences; Bensoussan, A., Lions, J.L., Eds.; Springer: Berlin/Heidelberg, Germany, 1988; Volume 111, pp. 543–556. [Google Scholar]
  21. Mathur, N.; Srivastava, P.K.; Paul, A. Algorithms for solving fuzzy transportation problem. Int. J. Math. Oper. Res. 2018, 12, 190. [Google Scholar] [CrossRef]
  22. Xu, S.J.; Nourinejad, M.; Lai, X.; Chow, J.Y. Network Learning via Multiagent Inverse Transportation Problems. Transp. Sci. 2018, 52, 1347–1364. [Google Scholar] [CrossRef] [Green Version]
  23. Jain, S.J.S. An Inverse Capacitated Transportation Problem. IOSR J. Math. 2013, 5, 24–27. [Google Scholar] [CrossRef]
  24. Stuart, A.M.; Wolfram, M.-T. Inverse Optimal Transport. SIAM J. Appl. Math. 2020, 80, 599–619. [Google Scholar] [CrossRef] [Green Version]
  25. Nguyan, T.D. Applictions of Robust and Inverse Optimization in Transportation. Master’s Thesis, Massachusetts Institute of Technology, Cambridge, MA, USA, September 2010. [Google Scholar]
  26. Sanjay, J. An Inverse Transportation Problem with the Linear Fractional Objective Function. AMO Adv. Model. Optim. 2013, 15, 677–678. [Google Scholar]
  27. Yang, D.; Gao, D. Study on Some Kinds of the Inverse Transportation Problems. In Proceedings of the Eighth International Conference of Chinese Logistics and Transportation Professionals (ICCLTP), Chengdu, China, 8–10 October 2008. [Google Scholar]
  28. Abo-Elnaga, Y.; Mousa, A.A. Inverse optimization for nonlinear programming problems. Wulfenia J. 2014, 21, 307–317. [Google Scholar]
  29. Bazaraa, M.S.; Sherali, H.D.; Shetty, C.M. Nonlinear Programming: Theory and Algorithms, 3rd ed.; Wiley: Hoboken, NJ, USA, 2006. [Google Scholar]
Figure 1. The transportation problem.
Figure 1. The transportation problem.
Mathematics 08 02027 g001
Table 1. Transportation problem tableau.
Table 1. Transportation problem tableau.
DemandWHS 1WHS 2WHS 3Supply
Supply
FAC 1 λ 1 x 11 2 x 12 7 x 13 7 x 13 13
x 11 x 12 x 13
FAC 2 λ 2 x 21 9 x 22 3 x 23 15
x 21 x 22 x 23
Demand1131428
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Mousa, A.A.A.; Abo-Elnaga, Y. Stability of Solutions for Parametric Inverse Nonlinear Cost Transportation Problem. Mathematics 2020, 8, 2027. https://doi.org/10.3390/math8112027

AMA Style

Mousa AAA, Abo-Elnaga Y. Stability of Solutions for Parametric Inverse Nonlinear Cost Transportation Problem. Mathematics. 2020; 8(11):2027. https://doi.org/10.3390/math8112027

Chicago/Turabian Style

Mousa, Abd Allah A., and Yousria Abo-Elnaga. 2020. "Stability of Solutions for Parametric Inverse Nonlinear Cost Transportation Problem" Mathematics 8, no. 11: 2027. https://doi.org/10.3390/math8112027

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop