A New Methodology for Solving Multiobjective Chance-Constrained Problems: An Application on IoT Systems

)is 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 aij. Twomethods are proposed to obtain the deterministic equivalent of chance-constrained problem. )e 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. )e 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.


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.
is 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 socalled "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. erefore, 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]. e classical linear programming problem, which is a specific class of mathematical programming problems, is formulated as follows: max z(x) � n j�1 c j x j , n j�1 a ij x j ≤ b i , i � 1, ..., m, x j ≥ 0, j � 1, ..., n. (1) Note that all coefficients (i.e., technologic coefficients a ij , right-hand side values b i , and objective function coefficients c j for j � 1, .., n and i � 1, ..., m) 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.
is 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 chanceconstrained problem, the original constraints n j�1 a ij x j ≤ b i , i � 1, ..., m, are replaced with the chance-constrained programming formulation P( n j�1 a ij x j ≤ b i ) ≥ 1 − δ i , i � 1, ..., m. e term "chance-constrained" is used because of the constraints of P(.), which designate the probability of the event in parentheses. e parameters δ i (0 < δ i ≤ 1) are assumed to be deterministically fixed. Chance-constrained programming is one of the major approaches for addressing random coefficients in mathematical programming problems [1,[4][5][6][7][8][9].
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. Agpak 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. ey 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 chanceconstrained goal programming. Zhang et al. [29] investigated the fixed-charge solid transportation problem under uncertainty.
ree models under different criteria were presented, which include expected value model, chanceconstrained 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. e 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. e 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 chanceconstraint problems with gamma random variables. e 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. e 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. ey 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 chanceconstraint 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 a kj are gamma random variables, n j�1 a kj x j 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 F(.) 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. erefore, probabilistic chanceconstraint corresponding to the direct distribution function is itself deterministic equivalent.
In this paper, we consider a multiobjective chanceconstrained stochastic programming problem, where the technologic coefficients a ij are independent gamma random variables with known means and variances. e 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 α 2 β parameters converge to normal distribution with parameters μ � αβ and σ 2 � α 2 β and β ⟶ ∞ and X ∼ N(μ, σ 2 ) [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. is 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. e 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. ese devices have limitations on energy usage, data transfer rate, and continuous usage time. ese 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.

Chance-Constrained Programming
An ordinary linear programming model is said to be chanceconstrained 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. us, 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: (2)

Mathematical Problems in Engineering
Note that δ i are specified probabilities and that c lj , l � 1, 2, ..., K, a ij , and b i 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 a ij parameters are gamma random variables. Here, it is assumed that the decision variables x i are deterministic.

The Deterministic Equivalent of Chance-Constraint Models
In this section, we derive a deterministic equivalent of model (2), where a ij , j � 1, ..., n, are independent random variables with known means and variances and c lj and b i are constants. e k th chance constraint is given in model (2).
Note that each a kj 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 n j�1 a kj x j , which is the linear combination of a kj . 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. e 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, n ≥ 5), it becomes very difficult to obtain the probability P( n j�1 a kj x j ≤ b k ). erefore, in this case Lyapunov's central limit theorem can be used to obtain an approximate solution.

Method I: Transformation of a Chance-Constrained
Programming Model into a Deterministic Model Using the Distribution Directly-Multivariate Transformation. Let X 1 , X 2 , . . . , X n be a random vector with a probability density which is defined as follows: Suppose that A 0 , A 1 , . . . , A n is from a partition of A. e set A 0 , which may be an empty set, satisfies en, for each j, the inverse function from B to A j can be obtained. Denote the j th inverse as follows: . , x n )). Let J j denote the Jacobian computed from the j th inverse. Assuming that these Jacobians do not vanish identically on B, we obtain the following representation of the joint probability density To obtain the deterministic equivalent of chance constraint (3), transformation methods can be used. Let a k1 , a k2 , ..., a kn be independent random variables. Consider the random variables (U 1 , U 2 , . . . , U n ), which are defined as follows: By deriving the probability density functions of the random variables U 1 , U 2 , ..., U n , the deterministic equivalents of chance constraints can be obtained. Hence, chance constraint (3) can be converted to a deterministic equivalent as follows:

Method II: Transformation of a Chance-Constrained
Programming Model into a Deterministic Model Using Lyapunov's Central Limit eorem. Let X 1 , X 2 , ..., X n be independent random variables. Suppose that each X j has a finite expected value and finite variance. Suppose also that the third central moments satisfy .., n) and the following Lyapunov condition [35]: en, the normalized sum of the random variables X j converges to a standard normal random variable as n ⟶ ∞.

4
Mathematical Problems in Engineering e 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 a ij are independent gamma variables with different parameters.

The Deterministic Model for the Gamma Distribution
Definition 1. Let a k1 , a k2 , ..., a kn be independent random variables with parameters (α kj , β kj ). e probability density function for each random variable a kj can be stated as follows: Method I: we start by considering a mathematical model with only two random variables. For j � 1, 2, chance constraint (3) can be expressed as follows: Note that a k1 and a k2 are independent gamma random variables with Γ(α k1 , β k1 ) and Γ(α k2 , β k2 ), respectively. To calculate the deterministic equivalent of chance constraint (13), we must first find the density of the random variable a k1 x 1 + a k2 x 2 . erefore, we define random variables U 1 and U 2 as follows: Note that the random variables Y j have a gamma distribution, with Γ(α kj , β kj x j ) � Γ(α kj , β * kj ). Let V 2 be an auxiliary variable used to find the probability density function of U 2 .
Note that if y 1 � u 2 − v 2 and y 2 � v 2 , then |J| � 1. e joint probability density function of U 2 and V 2 is as follows: Mathematical Problems in Engineering e marginal probability density function of U 2 can be computed as follows: After integrating the above integral, we obtain the following expression: Note that the term hypergeom is the generalized hypergeometric function defined by with By substituting equation (19) into (18), the probability density function of U 2 can be rewritten as follows: e distribution function of U 2 , which is the left side of the chance constraint given in (13) can be rewritten as follows: us, the deterministic equivalent of chance constraint (13) is as follows: We now consider a chance-constraint with three random variables. e k th chance constraint involving these random variables can be presented as follows: Note that a k1 , a k2 , and a k3 are gamma random variables with Γ(α k1 , β k1 ), Γ(α k2 , β k2 ), and Γ(α k3 , β k3 ), respectively. We define the random variable U 3 and auxiliary variable V 3 as follows: Note that if u 2 � u 3 − v 3 and y 3 � v 3 , then |J| � 1. e joint probability density function of U 3 and V 3 is as follows: e marginal probability density function of U 3 can be computed as follows: e distribution function of U 3 , 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 n random variables. To derive the marginal probability density function of U n , we define the random variable U n and auxiliary variable V n as follows: U n � a k1 x 1 + a k2 x 2 + · · · + a kn x n � Y 1 + Y 2 + · · · + Y n � U n−1 + Y n , V n � Y n � a kn x n . (29) Note that if u n−1 � u n − v n and y n � v n , then |J| � 1. e joint probability density function of U n and V n can be stated as follows: e marginal probability density function of U n can be computed as follows: Hence, the deterministic equivalent of the left side of chance constraint (3) is us, we have reached a generalized result for n random variables. e deterministic equivalent of chance constraint (3) can be easily calculated with equation (32) for known parameters of gamma distributions.

Mathematical Problems in Engineering
Method II: in this method, we apply Lyapunov's central limit theorem to the gamma distributions with large values of n. As such, we are able to transform chance constraint (3) to deterministic equivalent (11).

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 f l (x * ), l � 1, ..., K, which is the feasible solution of each lth objective function with respect to f l (x). 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 p � 1 or p � 2. While the value of p is increasing, total distance between f l (x * ) and f l (x) is decreasing.

Application of the Newly Proposed Method to an IoT End-System
e Internet of ings [39] is an Internet-based information infrastructure which enables several smart "things" to communicate with each other and servers. ese "things" are in general small devices with computing and communicating skills where most of the time measures some value. ey need a power supply to operate, and since there are distributed, they often rely on battery power. is energy limitation requires some kind of intelligent power consumption and thus a good way of battery usage optimization. e battery usage optimization process has a number of constraints, and the result should maximize the total amount of data sending time. e 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. e 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. us, maximizing the data sending time directly improves the quality of server-side jobs. ere 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.
e 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. e total of this time cannot exceed 3 seconds a day because the battery has a 3-second limit daily. e 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. e total of this time cannot exceed 2 seconds without changing the measurement equipment because of its limitations. e 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. e total of this time cannot exceed 1 second on consecutive sending operations because the device cannot handle that much memory. is application requires 3 objective functions. e 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. e 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. e third objective function includes sending the data with data overheads which include the data types such as measurement or location. e final one has the most transmitted data at each time, and the objective is maximizing these transmitted data overall. e 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: 8 Mathematical Problems in Engineering Note that the objective function coefficients, the righthand side values, and the specified probabilities are deter-  Table 1.
From Method I, we obtain the following deterministic problem: Note that nonnegative constraint x j ≥ 0 is replaced with x j > 0 to avoid dividing by zero in equation (36). From Method II, we obtain the following deterministic problem: Mathematical Problems in Engineering en, we solve the deterministic models (36) and (37) for each objective function and derive the solution shown in Table 2.
en, 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 p � 1 and p � 2 in the multiobjective model. e 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. is method uses an approach which converts a multiobjective problem into a single-objective problem. us, it is possible to achieve a single solution for the multiobjective problem [40]. e 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 d − l and d + l , l � 1, 2, 3 , 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. e result of the goal programming solution is given in Table 4.
Membership function method is used for Pareto optimal solution. e procedure is as follows: first, structure membership function μ l (f l (x)) is generated for every objective function f l (x). Second, μ l (f l (x)) is used as the new objective function for the structure of the new multiobjective programming problem. en, 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: , l � 1, 2, . . . , K.
Pareto optimal solutions for model (36) is given in Table 5.
If model (36) is analyzed for Pareto optimal solutions, optimal solution for f 1 is obtained when w 1 � 1 and w 2 � w 3 � 0. According to this, functions are resulted as f 1 � 0.408286, f 2 � 2.24557, and f 3 � 2.857998. Optimal solution for f 2 is obtained when w 2 � 1 and w 1 � w 3 � 0. According to this, functions are resulted as f 1 � 0.405686, f 2 � 2.247886, and f 3 � 2.873029. Optimal solution for f 3 is obtained when w 3 � 1 and w 2 � w 1 � 0. According to this, functions are resulted as f 1 � 0.396843, f 2 � 2.236849, and    (36) are given in Figure 1.
If model (37) is analyzed for Pareto optimal solutions, optimal solution for f 1 is obtained when w 1 � 1 and w 2 � w 3 � 0. According to this, functions are resulted as f 1 � 0.3598632, f 2 � 1.919270, and f 3 � 2.399088. Optimal solution for f 2 is obtained when w 2 � 1 and w 1 � w 3 � 0. According to this, functions are resulted as f 1 � 0.3546770, f 2 � 1.946318, and f 3 � 2.473929. Optimal solution for f 3 is obtained when w 3 � 1 and w 2 � w 1 � 0. According to this, functions are resulted as f 1 � 0.3458346, f 2 � 1.935281, and f 3 � 2.487224. Graphical representations of Pareto optimal solution for model (37) are given in Figure 1.

Results and Discussions
In this paper, we have considered a chance-constrained programming problem under the assumption that the random variables a ij are independent with gamma distributions. e 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 n independent random variables. en, using the probability density function, the chance constraints are transformed into deterministic constraints. For large values of n, 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 n � 2, 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 n.
For the multiobjective case, solutions are presented in Table 5. In the case of p � 2, 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. e 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: 1 st objective function: Method I (p � 2) is the optimal solution, where f 1 � 0.3802155�0.38. e result demonstrates that, with sending minimum data as only integer value, the maximum total data sent can reach 0.38 MB of data.  (36)). e 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
e datasets generated and/or analyzed during the current study are available from the corresponding author on reasonable request.

Conflicts of Interest
e authors declare that they have no conflicts of interest.   (37)).