Abstract

This study presents a newly developed methodology to transform the chance-constrained problem into a deterministic problem and then solving this multiobjective deterministic problem with the proposed method. Chance-constrained problem contains independent gamma random variables that are denoted as . Two methods are proposed to obtain the deterministic equivalent of chance-constrained problem. The first of the methods is directly based on using the distribution, and the second consists of normalizing probabilistic constraints using Lyapunov’s central limit theorem. An algorithm which uses the Global Criterion Method is developed to solve the multiobjective deterministic equivalent of chance-constrained problem. The methodology is applied to a real-life engineering problem that consists of an IoT device and its data sending process. Using Lyapunov’s central limit theorem for large numbers of random variables is found to be more appropriate.

1. Introduction

Most decision problems, whether real world or hypothetical, can be defined with a mathematical programming model. Unfortunately, in many practical systems, it is very difficult to gather all necessary information to derive such models. This difficulty arises from the uncertainty of many of the factors in modeling equations. One possible reason for this uncertainty is the randomness of these factors, and stochastic programming is a branch of mathematical programming that theoretically and empirically focuses on “conditional extremum problems under incomplete information about random coefficients” [1].

In terms of the investigation and optimization of complicated systems, the deterministic approach broadly dominates projection, planning, and management processes in economic and military sciences. However, by assuming that all initial data exhibit the same basic characteristics, the deterministic approach ignores the probabilistic properties of complicated systems. On the one hand, these deterministic systems emphasize objective laws. On the other hand, the deterministic approach is constantly undermined by random factors in real-life data. As such, a stochastic approach should be more useful for most real-life problems [2].

Assume that a decision-maker wants to find the so-called “absolutely best” decision that corresponds to the minimum (or maximum) of a suitable objective function, while satisfying a given set of feasibility constraints. If the decision models are very complex, standard optimization strategies cannot be used to solve them directly. Therefore, a global optimization technique should be used. Global optimization applications are popular in all disciplines, including science, engineering, economy, and applied science. Global optimization is aimed at finding the best global solution of a model, in the presence of multiple local optima [3].

The classical linear programming problem, which is a specific class of mathematical programming problems, is formulated as follows:

Note that all coefficients (i.e., technologic coefficients , right-hand side values , and objective function coefficients for ) are deterministic. However, when at least one coefficient is a random variable, the problem becomes a stochastic programming problem.

Stochastic programming problems can be classified into the following three categories: probability distribution problems, chance-constrained problems, and recourse problems.

This paper attempts to analyze a class of chance-constrained stochastic programming problems. Chance-constrained programming assumes the presence of random data variations and permits constraint violations up to specified probability limits. To transform equation (1) into a chance-constrained problem, the original constraints , , are replaced with the chance-constrained programming formulation , . The term “chance-constrained” is used because of the constraints of , which designate the probability of the event in parentheses. The parameters are assumed to be deterministically fixed. Chance-constrained programming is one of the major approaches for addressing random coefficients in mathematical programming problems [1, 49].

Chance-constrained programming was first developed by Charnes and Cooper [5]. A variant of the Charnes and Cooper model, which utilized joint chance constraints, is given by Miller and Wagner [10]. Prekopa [11] introduced a more general case in which the right-hand side is allowed to have stochastically dependent components. Contini [12] developed an algorithm for stochastic goal programming which assumes that the random variable coefficients are normally distributed with known means and variances. In doing so, he transformed the stochastic problem into an equivalent deterministic quadratic programming problem. Rather than assuming that the coefficients are normally distributed, Sengupta [8] assumed a chi-squared distribution, which has a nonnegative range and is more applicable to economic models of production and investment. Jagannathan [7] obtained equivalent equations for chance-constrained programming models based on coefficient matrices with normally distributed elements and dependent random right-hand side elements. Regarding multiobjective stochastic linear programming problems, Stancu-Minasian [13] suggested a minimum risk approach, while both Leclercq [14] and Teghem et al. [15] have proposed interactive methods. Extensive discussion on the convex properties of stochastic problems can be found in [16]. Meanwhile, a fuzzy approach to stochastic programming has been presented by Mohan and Nguyen [17]. An interactive fuzzy satisfying method was also developed by Sakawa et al. [18]. Yang and Wen [19] presented a chance-constrained programming model for transmission system planning in the competitive electricity market environment. Huang [20] provided two types of credibility-based chance-constrained models for portfolio selection with fuzzy returns. Ağpak and Gökçen [21] developed new mathematical models for stochastic traditional and U-type assembly lines with a chance-constrained 0-1 integer programming technique. Henrion and Strugarek [22] investigated the convexity of chance constraints with independent random variables. Parpas et al. [23] proposed a stochastic algorithm for the global optimization of chance-constrained problems. They assumed that the probability measure used to evaluate the constraints is known only through its moments. Xu et al. [24] developed a robust hybrid stochastic chance-constraint programming model to support municipal solid waste management under uncertainty. Ben Abdelaziz and Masri [25] proposed a chance-constrained approach and a compromise programming approach to transform the multiobjective stochastic linear program with partial linear information on the probability distribution into its equivalent uniobjective problem. Goyal and Ravi [26] presented a polynomial time approximation scheme for the chance-constrained knapsack problem when item sizes are normally distributed and independent of other items. Zhang and Peng [27] proposed uncertain optimal assignment problem in which the uncertain factor is described by uncertainty theory. In order to construct an uncertain programming model, the concept of α-optimal assignment is proposed. After that, -optimal model is given. Chen et al. [28] investigated an uncertain bicriteria solid transportation problem. Based on these equivalence relations, the optimal transportation plans were obtained by expected value goal programming and chance-constrained goal programming. Zhang et al. [29] investigated the fixed-charge solid transportation problem under uncertainty. Three models under different criteria were presented, which include expected value model, chance-constrained programming model, and measure-chance programming model. Majumder et al. [30] proposed an uncertain multiobjective Chinese postman problem under the framework of uncertainty theory, for an undirected connected network. The model is then solved using two classical multiobjective techniques, the Global Criterion Method and fuzzy programming method, and two multiobjective genetic algorithms. Majumder et al. [31] proposed a biobjective rough‐fuzzy quadratic minimum spanning tree problem which is modeled using rough‐fuzzy chance‐constrained programming technique. The linear and quadratic weights of the proposed model are considered as two different objectives which are optimized simultaneously. Atalay and Apaydın [4] proposed a new method using Esseen inequality to obtain deterministic equivalent of chance-constraint problems with gamma random variables. The method in this paper uses estimation of the distance between distribution of sum of these independent random variables having gamma distribution and normal distribution; probabilistic constraint obtained via Esseen inequality has been made deterministic using the approach suggested by Polya. The main purpose of the article is to transform the chance-constrained model into a deterministic model based on the Esseen inequality. As a result, the upper bounds of the chance constraints are derived by the Esseen inequality and approximate deterministic equivalent of the model is developed. Sarkar et al. [32] have reduced the chance constraints with independent gamma parameters to deterministic constraints. They have used geometric inequality and some other inequalities for obtained deterministic equivalent of chance constraints where the objective function may be considered as linear or nonlinear with deterministic cost coefficient. In contrast to the two papers discussed above, in our paper, distribution function is directly used and multiobjective mathematical model solution is obtained with the use of Global Criterion Method. Distribution function is on the left-hand side of the chance constraint. Besides, Pareto optimal solutions are also provided. Furthermore, deterministic equivalent of chance constraint which is derived by direct distribution function is compared to the approximate model which is derived using Lyapunov’s central limit theorem.

To the best of our knowledge, there is not any study that proposes an algorithm for solving multiobjective chance-constraint problems with gamma random variables by directly using the distribution and then using Lyapunov’s central limit theorem and Global Criterion Method. Instead of using quantile rules, we directly used the distribution function. In this case, since are gamma random variables, also becomes gamma random variables. Instead of using quantile function, distribution function of sum of independent random variables is directly used. According to this, distribution function is obtained and it is assigned as the left-hand side value of the chance constraint. In the mathematical model, the parameters of the chance constraint which are random variables have a distribution function and this distribution function provides the direct calculation of probability. Therefore, probabilistic chance-constraint corresponding to the direct distribution function is itself deterministic equivalent.

In this paper, we consider a multiobjective chance-constrained stochastic programming problem, where the technologic coefficients are independent gamma random variables with known means and variances. The decision why we use the gamma distribution is motivated by the fact that it has relationships with gamma and normal distributions. If univariate distribution relationship is examined, Gamma random variables consisting and parameters converge to normal distribution with parameters and and and [33]. Bivariate gamma distributions, such as by marginal transformation, are mainly derived from the bivariate normal in some fashion. A univariate chi-squared distribution can be derived from one or more independent and identically distributed normal variables. Moreover, a chi-squared random variable is a special case of gamma. As a result of this, a bivariate gamma model is related to the bivariate normal one [34].

In addition, other distributions such as the exponential and chi-square distributions can be closely approximated by a gamma distribution. For our chance-constrained model, we developed a deterministic equivalent model without random variables. A transformation was implemented to convert the chance-constrained programming problem into a deterministic equivalent problem. This deterministic equivalent problem is obtained by two methods, that is, by deriving the direct distribution and then normalizing probabilistic constraints using Lyapunov’s central limit theorem [35]. After converting the problem into a deterministic model, we used global optimization to propose an algorithm that utilizes the Global Criterion Method to solve the multiobjective problems. The solution was applied to a real-life engineering problem. IoT devices and distributed computing are recently a hot topic, and there are fairly new problems arising rapidly following the emerging new hardware and software running on these IoT devices. These devices have limitations on energy usage, data transfer rate, and continuous usage time. These limitations and the need for maximum data transfer amount make the application a suitable problem for the model. We have applied the model for such an IoT application in this study. Using numerical examples, we compared the objective functions of models obtained from these exact and approximate methods.

2. Chance-Constrained Programming

An ordinary linear programming model is said to be chance-constrained if its linear constraints are associated with a set of probability measures that indicate the degree to which constraints can be violated [36]. As a stochastic programming method, chance-constrained programming includes fixing certain appropriate levels for random constraints. Thus, it is generally used for modeling technical or economic systems, particularly economic planning, input control, structural design, inventory, and air and water quality management problems. Under chance-constrained programming, each constraint can be realized with a certain probability.

When some or all of the coefficients of the model are independent random variables, the chance-constrained programming formulation requires only that each constraint must hold for most of these combinations [6, 8, 37].

A multiobjective chance-constrained programming problem can be stated as follows:

Note that are specified probabilities and that , , , and are independent random variables. A general CCP problem is defined with model (2). In the model, one of the parameters, two of the parameters, or three of the parameters can be random variables. However, the main focus of this work is on the case where parameters are gamma random variables. Here, it is assumed that the decision variables are deterministic.

3. The Deterministic Equivalent of Chance-Constraint Models

In this section, we derive a deterministic equivalent of model (2), where , are independent random variables with known means and variances and and are constants. The kth chance constraint is given in model (2).

Note that each has a finite third central moment. If the deterministic equivalent of constraint (3) can be derived, the problem can be converted by determining the distribution of , which is the linear combination of . To derive the deterministic equivalent of the chance constraint (3), we propose two methods, namely, deriving the direct distribution and using Lyapunov’s central limit theorem. Both methods are first defined, and then they are adapted to the gamma distribution.

The most suitable method for deriving the deterministic equivalent of a chance-constrained programming model is to find the direct distribution. However, as the number of decision variables increases to a large value of n (more precisely, ), it becomes very difficult to obtain the probability . Therefore, in this case Lyapunov’s central limit theorem can be used to obtain an approximate solution.

3.1. Method I: Transformation of a Chance-Constrained Programming Model into a Deterministic Model Using the Distribution Directly—Multivariate Transformation

Let be a random vector with a probability density function . Let . Consider a new random vector, , which is defined as follows:

Suppose that is from a partition of . The set , which may be an empty set, satisfies . The transformation is a one-to-one transformation from onto for each . Then, for each , the inverse function from to can be obtained. Denote the jth inverse as follows:

For , the jth inverse yields a unique such that . Let denote the Jacobian computed from the jth inverse. Assuming that these Jacobians do not vanish identically on , we obtain the following representation of the joint probability density function for [38]:

To obtain the deterministic equivalent of chance constraint (3), transformation methods can be used. Let be independent random variables. Consider the random variables which are defined as follows:

By deriving the probability density functions of the random variables , the deterministic equivalents of chance constraints can be obtained. Hence, chance constraint (3) can be converted to a deterministic equivalent as follows:

3.2. Method II: Transformation of a Chance-Constrained Programming Model into a Deterministic Model Using Lyapunov’s Central Limit Theorem

Let be independent random variables. Suppose that each has a finite expected value and finite variance. Suppose also that the third central moments satisfy and the following Lyapunov condition [35]:

Then, the normalized sum of the random variables converges to a standard normal random variable as .

The deterministic equivalent of chance constraint (3) can be obtained using normalization as follows:

In the following section, we provide a deterministic equivalent of constraint (3) for the two methods mentioned above when coefficients are independent gamma variables with different parameters.

4. The Deterministic Model for the Gamma Distribution

Definition 1. Let be independent random variables with parameters . The probability density function for each random variable can be stated as follows:Method I: we start by considering a mathematical model with only two random variables. For , chance constraint (3) can be expressed as follows:Note that and are independent gamma random variables with and , respectively. To calculate the deterministic equivalent of chance constraint (13), we must first find the density of the random variable . Therefore, we define random variables and as follows:Note that the random variables have a gamma distribution, with  = . Let be an auxiliary variable used to find the probability density function of .Note that if and , then . The joint probability density function of and is as follows:The marginal probability density function of can be computed as follows:After integrating the above integral, we obtain the following expression:Note that the term is the generalized hypergeometric function defined bywithBy substituting equation (19) into (18), the probability density function of can be rewritten as follows:The distribution function of , which is the left side of the chance constraint given in (13) can be rewritten as follows:Thus, the deterministic equivalent of chance constraint (13) is as follows:We now consider a chance-constraint with three random variables. The kth chance constraint involving these random variables can be presented as follows:Note that , , and are gamma random variables with , , and , respectively. We define the random variable and auxiliary variable as follows:Note that if and , then . The joint probability density function of and is as follows:The marginal probability density function of can be computed as follows:The distribution function of , which is on the left side of the chance constraints given in (24), can be written as follows:Now, we consider a chance constraint involving random variables. To derive the marginal probability density function of , we define the random variable and auxiliary variable as follows:Note that if and , then . The joint probability density function of and can be stated as follows:The marginal probability density function of can be computed as follows:Hence, the deterministic equivalent of the left side of chance constraint (3) isThus, we have reached a generalized result for random variables. The deterministic equivalent of chance constraint (3) can be easily calculated with equation (32) for known parameters of gamma distributions.
Method II: in this method, we apply Lyapunov’s central limit theorem to the gamma distributions with large values of . As such, we are able to transform chance constraint (3) to deterministic equivalent (11).

5. A New Algorithm for Multiobjective Chance-Constrained Problem Procedure

We now present a methodology for solving a multiobjective chance-constrained stochastic programming model. To derive a deterministic model, we use the proposed Methods I and II, as presented in Section 4.Step 1: convert the given chance-constraint programming problem into an equivalent deterministic programming problem as discussed in Section 4.Step 2: solve the deterministic problem obtained in Step 1, using only one objective at a time in model (2) and ignoring the others.Step 3: calculate , which is the feasible solution of each th objective function with respect to . Find the corresponding values of the objective functions at each solution for comparison.Step 4: using the solutions obtained from Step 3 and the deterministic constraint obtained from Step 1, solve the following global criterion problem:Set or . While the value of is increasing, total distance between and is decreasing.

6. Application of the Newly Proposed Method to an IoT End-System

The Internet of Things [39] is an Internet-based information infrastructure which enables several smart “things” to communicate with each other and servers. These “things” are in general small devices with computing and communicating skills where most of the time measures some value. They need a power supply to operate, and since there are distributed, they often rely on battery power. This energy limitation requires some kind of intelligent power consumption and thus a good way of battery usage optimization.

The battery usage optimization process has a number of constraints, and the result should maximize the total amount of data sending time. The data sending time is directly related to the accuracy of updating the current state of the nodes on the server. According to this, at each sending process, the node sends either the current measurement or its current location to the server. The measurement data contain a measured height in centimeters, and location data contain world coordinates of the node. Once the data are sent to the server, the value on the server is updated and the number of data sending processes provides a better result on the server side for further calculations. Thus, maximizing the data sending time directly improves the quality of server-side jobs. There are certain constraints while maximizing the data sending time. Height change can occur with a random frequency through a day. Similarly, device location can change with a random frequency. With these in mind, constraints are as follows.

The first constraint is as follows: height change occurs with a random frequency every day and the assumption is that the height change is a gamma random variable with the parameters provided in Table 1. Every time a static duration of time is required for sending process. Similarly, location change occurs with a random frequency every day, and every time a static duration of time is required for sending process. The total of this time cannot exceed 3 seconds a day because the battery has a 3-second limit daily.

The second constraint is as follows: height change occurs with a random frequency every day and the assumption is that the height change is a gamma random variable with the parameters provided in Table 1. Every time a static duration of time is required for sending process. Similarly, location change occurs with a random frequency every day, and every time a static duration of time is required for sending process. The total of this time cannot exceed 2 seconds without changing the measurement equipment because of its limitations.

The third constraint is as follows: height change occurs with a random frequency every day and the assumption is that the height change is a gamma random variable with the parameters provided in Table 1. Every time a static duration of time is required for sending process. Similarly, location change occurs with a random frequency every day, and every time a static duration of time is required for sending process. The total of this time cannot exceed 1 second on consecutive sending operations because the device cannot handle that much memory.

This application requires 3 objective functions. The first objective function includes sending the data only as integer values where the total transmitted data at each time are the least and the objective is maximizing these transmitted data overall. The second objective function includes sending the data as float values where the total transmitted data are more than just sending integer values at each time and the objective is maximizing these transmitted data overall. The third objective function includes sending the data with data overheads which include the data types such as measurement or location. The final one has the most transmitted data at each time, and the objective is maximizing these transmitted data overall.

The mathematical model of the above engineering problem with multiobjective, chance-constrained programming problem by using the methodology presented in Section 5 is given as follows:

Note that the objective function coefficients, the right-hand side values, and the specified probabilities are determined as , , and , respectively. In this problem, is a gamma random variable with the characteristics shown in Table 1.

From Method I, we obtain the following deterministic problem:

Note that nonnegative constraint is replaced with to avoid dividing by zero in equation (36). From Method II, we obtain the following deterministic problem:

Then, we solve the deterministic models (36) and (37) for each objective function and derive the solution shown in Table 2.

Then, we apply the Global Criterion Method as the solution technique to consider the objectives together. Deterministic constraints that are obtained by these two methods were examined one by one for both and in the multiobjective model. The solutions obtained by substituting the values given in Table 2 in model (34) in Step 4 are shown in Table 3.

Goal programming (GP), proposed by Charnes and Cooper, is a classical method to solve multiobjective problems, which is widely used in the literature. This method uses an approach which converts a multiobjective problem into a single-objective problem. Thus, it is possible to achieve a single solution for the multiobjective problem [40]. The main idea of GP is to reach predefined goals for every single objective. Deviations of goals are tried to be minimized in GP because it is impossible to obtain a solution for every objective since the goals conflict with each other. Models (36) and (37) with and , , are the positive deviations and negative deviations that are modeled and solved as goal programming which is proposed as an alternative to the Global Criterion Method which is given in Table 3. The result of the goal programming solution is given in Table 4.

Membership function method is used for Pareto optimal solution. The procedure is as follows: first, structure membership function is generated for every objective function . Second, is used as the new objective function for the structure of the new multiobjective programming problem. Then, this converts the new multiobjective programming problem into a single-objective programming problem by using appropriate methods. Finally, single-objective programming problem is solved to get the optimal solution, which is also the Pareto optimal solution of the original multiobjective programming problem [41]. Structure of the membership function for every objective function is as follows:

Pareto optimal solutions for model (36) is given in Table 5.

If model (36) is analyzed for Pareto optimal solutions, optimal solution for is obtained when  = 1 and  =  = 0. According to this, functions are resulted as  = 0.408286,  = 2.24557, and  = 2.857998. Optimal solution for is obtained when  = 1 and  =  = 0. According to this, functions are resulted as  = 0.405686,  = 2.247886, and  = 2.873029. Optimal solution for is obtained when  = 1 and  =  = 0. According to this, functions are resulted as  = 0.396843,  = 2.236849, and  = 2.886325. Graphical representations of Pareto optimal solution for model (36) are given in Figure 1.

If model (37) is analyzed for Pareto optimal solutions, optimal solution for is obtained when  = 1 and  =  = 0. According to this, functions are resulted as  = 0.3598632,  = 1.919270, and  = 2.399088. Optimal solution for is obtained when  = 1 and  =  = 0. According to this, functions are resulted as  = 0.3546770,  = 1.946318, and  = 2.473929. Optimal solution for is obtained when  = 1 and  =  = 0. According to this, functions are resulted as  = 0.3458346,  = 1.935281, and  = 2.487224. Graphical representations of Pareto optimal solution for model (37) are given in Figure 1.

7. Results and Discussions

In this paper, we have considered a chance-constrained programming problem under the assumption that the random variables are independent with gamma distributions. The main purpose of this paper is to transform the chance-constrained model into a deterministic model based on two proposed methods. Method I provides exact solutions for the conversion of a chance constraint into a deterministic constraint. Using Method I, we found the probability density function of a linear combination of independent random variables. Then, using the probability density function, the chance constraints are transformed into deterministic constraints. For large values of , we proposed Method II, which is an approximation method. As can be seen in the example in Table 1, the numerical results for the exact and approximate models are rather close to each other. In this case, we first used Method I; however, because of the complex and difficult calculations in Method I, Method II can also be used.

As can be seen in Table 2, for multiobjective problems, we calculated the objective value by considering all objective functions and then summarized the solutions of the deterministic model under the two methods. We compared the solutions of the deterministic model according to the decision variables and objective functions. Even for the special case with , the objective values obtained from Method I and Method II are very similar. According to Lyapunov’s central limit theorem, it is obvious that we can obtain more approximate solutions for larger values of .

For the multiobjective case, solutions are presented in Table 5. In the case of , the decrease in deviations from the ideal solutions can be observed in objective values. As such, the numerical results from these two methods are very close to one another. For multiobjective problems, the Global Criterion Method provides more effective solutions after converting chance constraints to deterministic constraints.

Considering the results obtained from the real-life problem in Section 6, Table 3 is created as a summary. The multiobjective problems that are solved by Global Criterion Method are also solved by goal programming, and the results are given in Table 4. Also, Pareto optimal solutions for these problems are obtained, and the results are given in Tables 5 and 6, respectively. Graphical representations of these Pareto optimal solutions are also provided. According to the results provided in the table, optimal solutions for the tree objective functions can be summarized as follows:1st objective function: Method I () is the optimal solution, where  = 0.38021550.38. The result demonstrates that, with sending minimum data as only integer value, the maximum total data sent can reach 0.38 MB of data.2nd objective function: Method I () is the optimal solution, where  = 2.1917252.19. The result demonstrates that, with sending medium size data as float values, the maximum total data sent can reach 2.19 MB of data.3rd objective function: Method I () is the optimal solution, where  = 2.8651352.87. The result demonstrates that, with sending minimum data as float values and their headers, the maximum total data sent can reach 2.87 MB of data.

The obtained results show that the methods provided can be applicable to such problems and the solutions prove that an optimal design can be achieved by our study. By comparing the values of objective functions, the designers may benefit the most applicable solution to their problems.

Data Availability

The datasets generated and/or analyzed during the current study are available from the corresponding author on reasonable request.

Conflicts of Interest

The authors declare that they have no conflicts of interest.