Research Article  Open Access
Ianire Taboada, Fidel Liberal, "A Novel Scheduling Index Rule Proposal for QoE Maximization in Wireless Networks", Abstract and Applied Analysis, vol. 2014, Article ID 647157, 14 pages, 2014. https://doi.org/10.1155/2014/647157
A Novel Scheduling Index Rule Proposal for QoE Maximization in Wireless Networks
Abstract
This paper deals with the resource allocation problem aimed at maximizing users’ perception of quality in wireless channels with timevarying capacity. First of all, we model the subjective qualityaware scheduling problem in the framework of Markovian decision processes. Then, given that the obtaining of the optimal solution of this model is unachievable, we propose a simple scheduling index rule with closedform expression by using a methodology based on Whittle approach. Finally, we analyze the performance of the achieved scheduling proposal in several relevant scenarios, concluding that it outperforms the most popular existing resource allocation strategies.
1. Introduction
Undoubtedly, the use of mobile Internet applications has notably increased over the last years, which has led to the growth of the demand for wireless bandwidth. Hence, one of the fundamental challenges that network providers nowadays face is how to efficiently share radio resources among users’ traffic flows. In wireless links channel capacity evolves over time due to the intrinsic degradations of this medium, which has motivated the investigation of scheduling problems in timevarying channels.
Traditional scheduling strategies for resource allocation are oriented to objective quality parameters such as delay. Nevertheless, considering the importance and the necessity of network resource allocation for maximizing users’ subjective quality of service or quality of experience (QoE) [1], QoEdriven scheduling becomes essential for network providers.
Thus, motivated by the necessity of obtaining an implementable QoEaware scheduler in wireless channels with timevarying capacity, in this paper we aim at characterizing in closedform a novel channelaware scheduler for the problem of maximizing users’ perceived quality. We focus on a scenario where traffic flows arrive and depart upon service completion.
1.1. Background
Channelaware or opportunistic schedulers give priority to users in good channel conditions. Although several channelaware strategies exist in the literature, Max Rate and Proportional Fair [2] among the most popular, the achievement of the optimal solution for timevarying scheduling optimization problems is rather difficult and unknown. Without any doubt, most analyzed resource allocation problems aim at minimizing the average delay of user traffic flows [3–7].
The best contribution found in the field of optimizing delay in a timevarying context is found in [6], which proposes a simple indexbased heuristic scheduler that performs satisfactorily. In this work, the authors consider exponentially distributed flow sizes and formulate the optimal scheduling problem as a Markov decision process (MDP) [8]. Since this model is PSPACEhard [9], it does not allow obtaining a tractable solution, and the authors use approximations for its resolution. This way, using the theory of restless bandits [10] and Whittle approach [11], the authors develop an index rule called potential improvement (PI), which has the property of being asymptotically optimal and fluidoptimality [12]. Nonetheless, analyzing the problem of minimizing mean delay for generally distributed flow sizes in randomly timevarying channels remains still open.
Furthermore, even though a few QoE and channelaware scheduling proposals exist [13–15], no significant analytical framework is provided for these heuristics. In fact, there is no solution based on relevant mathematical methods such as MDP for modelling and Whittle method for its resolution in order to obtain a simple scheduling index rule.
On the other hand, it is worth mentioning the work carried out in [16], which is pioneer in introducing QoEawareness in a MDP. Moreover, since the proposed model is analytically unsolvable in general, they achieve an easy, tractable, and wellperforming QoEaware index rule based on Gittins approach [4] for constant channel capacity, which they denominate Gittins MOS (GM). It is known that Gittins proposed a method to obtain the socalled Gittins index which originally minimized mean delay for constant channel capacity. However, as concluded in [17], Gittins approach is unsuitable for solving problems with timevarying channel.
1.2. Scope and Contributions
Since the stochastic and dynamic resource allocation problem of subjective quality maximization in channels with timevarying capacity is analytically and computationally unfeasible for finding an optimal solution, the main objective of this work is to design a simple and tractable heuristic priority scheduling rule using analytical tools that have had a great contribution in the optimization area.
Thus, the main contribution of this paper is threefold.(i)Firstly, we propose a QoEaware MDP model in a timevarying channel context.(ii)Secondly, in order to achieve a tractable solution for the aforementioned QoE and channelaware model, we focus on designing a simple heuristic index rule using Whittle approach.(iii)Thirdly, although for many years exponential flow size distributions have been considered for traffic modelling in order to simplify the resolution of scheduling optimization problems, as a step forward, in this work we take into account size distributions that better capture the real world patterns.
The remainder of the paper is organized as follows. First of all, we characterize the QoEaware scheduling problem for the timevarying channel context in Section 2. Then, Section 3 describes the extension of the basic QoEaware MDP model to timevarying channels. In Section 4 we achieve the Whittle indexbased solution for the proposed MDP, and its performance is evaluated in Section 5. Finally, Section 6 gathers the main conclusions of this paper.
2. Problem Description
We analyze the problem of maximizing average QoE in timevarying channels. Even though this study is applicable to any timevarying channel context, we focus on wireless networks. In particular, we centre on a wireless downlink data channel in a single cell system. In this way, at the beginning of each transmission time interval (TTI), the scheduler located in the base station makes decisions in order to choose a user traffic flow to transmit.
2.1. User
Each user in the system is uniquely associated with the flow related to its requested content, with its QoE characteristic and with its wireless channel. We use a generic concept of flow, which refers to an amount of data that can be displayed to (received by) the end user as a standalone element.
We consider a system with a fixed number of users, without arrivals of new users. This assumption simplifies the mathematical model. However, being aware of the impact of arrivals on the performance of scheduling, we will analyze the performance of scheduling strategies in the presence of arrivals.
2.1.1. Flow Size
Flows are characterized by their random size , the total amount of bits to be transferred. Sizes are assumed to be independently distributed with , where , , and are the corresponding cumulative distribution function, density function, and hazard rate, respectively.
We focus on the important class of size distributions with a decreasing hazard rate, particularly we assume Pareto distributed flow sizes. It is known that Internet traffic flow sizes are properly modelled by means of Pareto distributions [18], which better fit Internet flow sizes compared with the simplistic exponential approach. We define the Pareto distribution with shape parameter and scale parameter whose density function for all is
Note that, even though we use Pareto distributed flow sizes, the results provided are valid for any size distribution with a decreasing hazard rate.
2.1.2. QoE Characteristic
The time needed to transmit to channel all the bits from a user flow has a direct impact on user’s perceived quality. We consider that this delay is the main cause of subjective quality distortion. Nowadays, most of the Internet traffic is transported over TCP, which implements the packet recovery mechanisms to avoid application level losses. In TCPbased services, the predominant source of losses at the application is packets arriving later than their playout time. Hence, it is generally useful to provide a delaydriven QoEawareness.
This way, we quantify QoE by a delay dependent Mean Opinion Score (MOS) [19] utility function. User satisfaction is evaluated in MOS scale of 1 (very poor quality) to 5 (excellent quality).
As shown in QoE studies [20], perceived quality is maximum until a low delay value, . Then, from this delay threshold, QoE degradation happens until a high delay value, . Above this delay threshold service is unaffordable, resulting in a MOS value of 1. The values of these delay thresholds depend on several factors [21], such as the flow size, the service type, the network type and the user expertise.
We illustrate an example of the explained behavior in Figure 1.
2.1.3. Channel Model
The channel quality associated to a user evolves randomly and independently of its evolution history. We assume that the channel can be in conditions. Therefore, user is in channel condition with probability , having . In channel state , bits are transmitted to channel, and transmission rates are .
2.2. Server
The objective of the scheduler is finding a policy that solves the following optimization problem: According to the already described aspects about the scheduling problem, we consider a set of admissible policies, , which fulfils the next properties.(i)Delay dependency: the time that a flow remains in the system, , is known.(ii)QoEawareness: the delaydependent MOS function that characterizes user subjective quality is known.(iii)Channelawareness: channel quality indicator (CQI) information sent from mobile users is used as user’s instantaneous channel condition information.(iv)Nonanticipation: similarly to most current IP systems, flow size is unknown. However, existent nonanticipating sizebased disciplines make use of flow attained service [4, 7], the bits that have been transferred of a flow. Due to the memoryless property of the exponential distributions, for these distributions, the hazard rate is independent of attained service. Nevertheless, the considered Pareto distributions have a decreasing hazard rate with attained service [22]; this is, completion probability, , is higher for lower attained service values. This way, attained service, , is used as size information.(v)Preemption: we assume that the server is preemptive, that is, at every decision epoch, it is permitted to suspend the service of a flow whose transmission is unfinished.(vi)Single service: only the transmission of a flow is allowed in each TTI.
3. MDP Formulation
In this section we formulate the problem described in Section 2 as a MDP. First, we provide the QoEaware MDP model of each user . Then, we formulate the optimization problem for the joint MDP model, which takes into account all the users in the system.
3.1. MDP Model of a Job
The action space, , of user is binary; action 0 means not serving, and if action 1 is chosen bits from this user being in channel state are transmitted to channel. Besides, each user is defined by tuple . These elements are defined as follows.(i) is the set of all states for a user , which is classified in three groups.(a)Unfinished states: all the states in which flow transmission is uncompleted belong to this group. These states have three dimensions, denoted as , with attained service level , delay level , and channel condition as components. The attained services that correspond to these attained service levels are multiples of the first not null transmission rate, whereas the delays that correspond to delay levels are multiples of TTI.(b)Reward states: these states are reached once flow transmission has finished, and in these delaydependent states the QoEaware reward is given. This way, aimed at maximizing average MOS, a delay dependant QoEaware reward giving only happens once the whole flow has been transmitted, since from subjective quality perspective the MOS function only makes sense once flow transmission is completed (see [16] for more details).(c)Final state: the chain of states ends in the absorbing state .(ii), where is the expected oneslot reward received for user at state if action is decided at the beginning of a TTI, where is normalized MOS in the range .(iii), where is the expected oneslot work done for user at state if action is decided at the beginning of a TTI, (iv), where is the probability of user of moving from state to state if action is decided at the beginning of a TTI, where . Observe that the probability of changing from channel state to channel state corresponds to the steadystate probability due to the memoryless property of channel evolution.
Thus, the dynamics of user is captured by state process and action process . Moreover, we highlight that the proposed MDP model has restless property, which admits evolution and rewards even if not chosen for transmission due to the delay increase and the channel stochastic evolution.
Figure 2 shows part of the peruser state diagram of the new QoEaware MDP model, which summarizes and relates the aforedescribed components.
3.2. Optimization Problem
Besides, the optimization problem (2) associated with such MDP can be written as Note that we consider the undiscounted reward case, where the discounted factor is 1, being reward deterioration in time implicit in the delaydependent MOS function.
Nonetheless, optimally solving problem (6) is unfeasible in general. On the one hand, when channel capacity is constant the impossibility of obtaining the exact QoEaware solution is revealed in [16], and hence, since the timevarying channel model is more complex, the achievement of the optimal solution is unviable in this case. On the other hand, it is known that the allocation constraint causes intractability in similar problems [11], and therefore, in our case having threedimensional states finding the optimal solution is extremely difficult or even unachievable. In this way, once seeing the impossibility of obtaining the optimal solution for the considered problem in this work, we will focus on achieving an approximate solution. For that purpose, in the next section we will propose a heuristic based on Whittle method.
4. QoEAware Index Rule Proposal Based on Whittle Approach
In this part we provide a Whittle index rule type approximate solution to problem (6) based on Whittle method. The idea of Whittle [11] consists in achieving a function that measures the dynamic service priority in order to obtain a simple scheduling index rule. Whittle indices are computed in isolated way for each user, and they measure the expected efficiency of serving a user in each state.
Note that we are only interested in the indices of states as once flow transmission is completed scheduling makes no sense. Thereby, the optimization problem formulated in (6) can be relaxed by requiring to serve a job per slot on average as proposed in [11], which is further approached by Lagrangian methods [23] and can be decomposed into a singlejob pricebased parameterized optimization problem. For a price or Lagrangiar parameter , we will therefore study the user subproblem: The Lagrangiar parameter can be interpreted as the perslot cost of serving.
Let us define serving set , which prescribes to serve a user if , while not to serve this user if . We will refer to states as active and as passive. This way, it is possible to rewrite problem (7) as where and are, respectively, the expected total reward and expected total work in state.
Conjecture 1. Problem (8) is indexable.
Due to the complexity of our model, proving the indexability of Problem (8) is not an easy task. In this way, we assume Conjecture 1, and therefore, we suppose that for each state of user the Whittle index exists. From now on, we omit user label . In spite of considering the undiscounted case, the mathematical analysis is carried out for the discounted case so as to avoid infinite values caused by terms in denominators for , and then, we obtain the undiscounted index values in the limit .
Following Whittle index definition in [10], we formally write the Whittle index for Problem (8), , as
Whittle index, , represents the rate between marginal reward and marginal work, where the marginal reward (work) is the difference of the expected reward earned (work required) by serving and not serving at the initial state and employing policy afterwards.
Lemma 2. For any state and under any policy we obtain
Proof. From the definition of reward and work, respectively, we have
By substituting the expressions (11) and (12) in formula (9), we achieve a more complete (10) formulation.
Nevertheless, in order to obtain an analytically tractable Whittle index expression for (10) in the state , it is necessary to determine the policy of the future states that have influence on the index computation of this state. As can be observed from (10), (11), and (12) expressions, in these future states the value of attained service is maintained or increased, whereas delay grows over time. Since increasing attained service decreases instant completion probability for a Pareto distribution and increasing delay degrades QoE, it would seem intuitively correct to suppose that future states with the same or worse channel conditions are passive. In such a way, we could say that in a state with the same or worse channel condition being attained service and delay higher that the Whittle index value is lower. However, for future states with better channel conditions their passive or active property is not trivial, because the improvement in channel condition does not guarantee being a better state due to the aforementioned behavior of the increasing attained service and delay. Therefore, finding the active set of all the influential states for Whittle index computation becomes a challenging task.
In such situation, aimed at achieving the characteristics of the optimal policy and the subsequent active set building evolution we will use an algorithm called AdaptiveGreedy, shortly algorithm (see [10] for a survey). This algorithm computes Whittle indices numerically taking into account optimal active sets. This way, we have implemented algorithm for our case study in order to verify the fundamental properties of the Whittle index for our problem, which will be useful to determine the active set, and consequently, a closedform index expression. We provide the operation of the applied algorithm in Algorithm 1.

We have performed several numerical experiments for case, in which we have reduced the state space of both delay and channel states so as to do the execution of the algorithm computationally and in time viably. We have only considered two channel conditions and eleven levels of delay in the scenarios. In such a way, from experimental results we conjecture the following fundamental properties.(1)The Whittle index value depends on attained service, delay, and channel state. Moreover, as can be seen in graphs from Figures 3 and 4, we can clearly distinguish the ranges of MOS function; above the upper delay threshold, the index value is null.(2)For the same channel condition and the same delay level, Whittle index values are decreasing with attained service: (illustrated in Figures 3 and 4).(3)For the same channel condition and the same attained service level, Whittle index values are decreasing with delay: . (See Figures 3 and 4).(4)In a better channel condition, for the same attained service level and the same delay level, the Whittle index value is higher: . (Observe Figure 3).(5)Opposite to the previous channelaware work [6], where Whittle indices in the best channel condition are infinite having absolute priority in this channel state, in our case being in the best channel condition does not guarantee that the Whittle index value is higher than in the rest of the channel conditions: it is possible to fulfil . If we compare the two graphs in Figure 3, it can be easily seen how for several attained services and/or delayed values the Whittle index values in the worst channel condition are higher than in the best channel condition (taking also into account nonnull values).(6)Dependency on channel condition probabilities: hile for channel condition the Whittle index value remains invariant with channel state probability, for the rest of channel conditions its value decreases as long as channel state probabilities increase in better channel conditions (shown in Figure 4).
(a)
(b)
(a)
(b)
Once we analyzed the properties of the Whittle index, we set out to define the structure of the active set in order to obtain an analytically and computationally tractable Whittle index expression. In such a way, to derive a closedform characterization of the Whittle index of state, first of all it is necessary to specify the activity of and states. As concluded from Whittle index properties, it is known that states are passive. However, it is not trivial to specify above which attained service level and/or delay level if states with better channel condition are active or passive.
Nonetheless, in the best channel the active set is totally specified since all the states that compute in the index are passive, because of nonexisting channel improvement and being in the future attained service equal or higher and delay higher. As shown in Proposition 3, for the best channel condition the Whittle index has a closedform expression, where the index value is the multiplication between completion probability and normalized MOS function. Note that the obtained Whittle index expression in the best channel state is equivalent to the GM index achieved for constant channel [16].
Proposition 3. The Whittle index for Problem (8) in the best channel condition is given by
Proof. It is known that for index computation that . In this way, using (11) and (12): . Hence, analogolously, , , and . Substituting these null terms in expression (10), and being , we achieve (13).
Besides, as stated before, obtaining a computationally tractable Whittle index expression for any channel condition is not possible. In such situation, in order to achieve a tractable Whittlebased QoEaware index rule for timevarying channels, we propose a Whittlebased approximate heuristic in the next subsection.
4.1. Approximate Whittle Index Proposal
In the following we propose an approximate solution for the QoEaware Whittle index in a timevarying channel context (10).
First of all, we will determine the structure of the active set. For the Whittle index computation of state, as previously concluded, it is known that states are passive, but we need to define the unknown activity of states. Thus, so as to determine if these states are active or passive we apply the following two considerations.(1) approximation: (2) approximation:
On the one hand, we assume that for any value that (14). It is known that states are active, and being , we could consider that states are active.
On the other hand, we assume (15), since QoE degradation in a higher delay level is negligible in the millisecond scale that delay evolves. Taking into account the approximation about MOS function, being states active, we can suppose that states are also active. This way, combining approximations (14) and (15), we consider that states are active.
Then, suppose that the aforementioned and MOS approximations cause and . So, using these reward and work approximations, the first simplified expression for the Whittle index formulation (10) is
In such a way, applying the previous simplifications, we obtain the Whittle index approximation presented in Proposition 4, what we call attained service PI MOS (ASPIM). The methodology to compute the closedform expression of this QoEaware Whittlebased index is presented in the Appendix.
Proposition 4. The formulation of the ASPIM index is:
If we analyze expression (17), we observe that this QoEaware and channelaware Whittlebased ASPIM index is implementable. This proposal depends on QoE, size, and channel properties. Moreover, the obtained index for the best channel condition is equal to the original Whittle index (13), since the summation in the denominator becomes null. However, as shown in (18), when channel state is not the best the index value goes to zero:
Hence, a user in its best channel state with nonnull instantaneous normalized MOS function has priority over a user which is not in its best channel condition. Thus, we summarize the proposed QoEaware Whittlebased ASPIM index rule in Definition 5.
Definition 5. The ASPIM index rule consists of: at every decision slot,(i)serving the user in its best channel condition with nonnull instantaneous normalized MOS function with the highest value of ;(ii)if there is no user in its best channel condition with nonnull instantaneous normalized MOS function, serving the user with the highest value of (using (18) normalized by ).
So, if channel condition is not the best, the index value is the ratio of the multiplication between the actual service completion and the instantaneous normalized MOS function respect to the expected potential improvement of the completion probability.
Therefore, in the context of timevarying channels aimed at maximizing average MOS, we have proposed ASPIM scheduling algorithm, which is a simple and tractable QoEaware, channelaware, and sizebased index rule. Nevertheless, in order to verify the correct behavior of this heuristic, we analyze its performance in the subsequent section.
5. Performance Evaluation
In this section, we evaluate the performance of the proposed ASPIM index rule (see proposal in Definition 5). To that end, we compare the QoE performance of our heuristic with the one achieved by wellknown scheduling strategies in several simulation scenarios. Below we provide a brief description of the scheduling priority rules used in our experiments.(i)Max rate (MR): this opportunistic policy consists in serving the user with the best channel capacity.(ii)Proportional fair (PF): this channelaware discipline serves the user with the highest ratio between instantaneous transmission rate and current attained throughput.(iii)Cost and attained service dependant (AS): although it might seem counterintuitive, minimizing mean delay does not result in maximizing average MOS [24]. This way, in order to compare the performance with delayaware solutions, we propose to adapt the classical rule aimed at minimizing mean delay [25] to attained service. This adaptation consists of serving the user with the highest value of , being the waiting cost paid for every slot while flow transmission is uncompleted.(iv)GM [16]: this QoEaware pioneering index rule, originally proposed for constant channel capacity, chooses the user with the highest value of .
In case of ties, we use random tiebreaking rule.
Concerning traffic flows in the system, we use the Pareto size distribution given in (1) in simulations, whose and parameters determine mean flow size. We consider typical mean sizes used in wireless networks: 0.5 Mb (small traffic: a web page, email), 5 Mb (medium traffic: a PDF document, a picture, etc.), and 50 Mb (MP3 audio, a group of pictures from a video sequence, etc.). User traffic flows arrive to the scheduler according to a Poisson process with rate.
Referring to network characteristics, we use transmission rates that depend on CQIs used in real wireless networks, which are associated with 4G modulation and coding schemes. We show the mapping between CQI values and transmission rates in Table 1. The presented values are adapted and extrapolated from [26] for a cell of 20 MHz.

Moreover, we consider different network loads in order to analyze the behavior of scheduling strategies under different network conditions. Flow arrival rate determines network load, , where . We consider seven network states in this study: low load (), medium load (), and high load ().
Regarding QoE aspects, as typically utilized in perceived quality studies [27], we use a logshaped MOS function for QoE degradation, presented in (19). This way, employing the MOS curve modelling provided in [21] as starting point and taking into account the flow size, the network type and the user expertise, we have suitably defined and delay thresholds used in simulations, and consequently and MOS coefficients. Consider
We have implemented the whole network environment for the scheduling of network resources in MATLAB. In relation to the performed experiments, for each scenario, for each combination of scheduling discipline and network load, we have carried out a set of 10 rounds of 10000 s length simulations. These rounds differ in the randomly pregenerated traces (input vector) of sizes and arrivals. Thus, not only the overall average MOS values will be provided but we will also include their 95% confidence intervals.
Besides, in order to guarantee that the obtained performance results are generally valid, we will analyze different simulation scenarios. We have chosen five relevant settings, which differ in QoE, size, channel, or/and cost characteristics. The parameters of these scenarios are summarized in Table 2. Note that in the first four scenarios, we restrict ourselves to the case of two channel conditions so as to get fundamental performance results.

We have carefully selected the scenarios used in simulations. First of all, we define a setting (Scenario 1) in which the values of all the parameters are typical or standard; the values of these parameters introduce low error in the approximations considered for our ASPIM scheduler proposal. Then, in order to evaluate worst performance cases for ASPIM, we have chosen three scenarios (Scenario 2, Scenario 3, and Scenario 4) in which the error of the approximation used for ASPIM is high due to a different factor (gradient of MOS, rate and size, resp.).
Apart from that, it is worth mentioning that for obtaining the ASPIM index rule that we have considered an isolated MDP model per user. Hence, the achieved solution is based on a unique user type or single class. However, it would be interesting to examine the proposal when there are different types of users, which is analyzed in Scenario 5. Furthermore, this last setting resembles a real 4G wireless network. Note that, even though this technology allows the simultaneous transmission of multiple flows per TTI, we assume that a single flow transmits in each TTI.
Next we describe the results we have achieved in the aforementioned scenarios.
5.1. Scenario 1: Basic Case
In this first family of simulations we consider a basic or typical scenario, which takes into account the equiprobable channel case and mediumsized selfsimilar flows. Figure 5 shows average MOS results for this setting employing different scheduling policies and network loads. The main conclusion achieved from these results is that ASPIM outperforms the rest of disciplines, which is better than the QoEaware GM rule proposed for constant channel capacity. Moreover, contrary to the other policies, it holds average MOS values in the same level for all . Referring to the traditional channelaware disciplines, both PF and MR give bad behavior comparing to the rest of scheduling strategies.
In order to justify the previous performance results we now analyze delay and MOS statistics for the highest network load considered. In such a way, as regarded in Figure 6(a) for the Probability Density Function (PDF) of delay, we highlight the probability mass or peak that appears for some scheduling disciplines before the low delay threshold (0.5 s). The height and weight of a peak belonging to a policy have direct impact on its average MOS. ASPIM discipline shows the highest and the widest peak among the presented ones, followed by GM but with a notable difference. Therefore, our novel QoEaware policy seems to try to finish the transmission of most of the flows before , which is reflected in the aforementioned peaks. Note that the cAS rule aimed at minimizing delay behaves in different manner, which shows more fairness among delays in the range shown compared to the QoEaware ones.
(a)
(b)
(c)
What is previously concluded about delay statistics has a direct reflection on QoE statistics. This way, as can be seen in the PDF of MOS presented in Figure 6(b), if the objective is to maximize average QoE it is important, on the one hand, to be the height and the weight of the peaks around the maximum MOS high, and on the other hand, to be the value around the minimum MOS low. ASPIM fulfils these positive properties, and consequently, it is the best in average QoE terms. Besides, if we analyze the Complementary Cumulative Distribution Function (CCDF) of MOS provided in Figure 6(c), we observe that ASPIM gives the highest probability of being MOS value higher. Considering, for example, a QoE satisfaction threshold, , of , in this point ASPIM improves GM a 10%, whereas the deterioration of MR and PF is higher than a 50%.
5.2. Scenario 2: High Demanding Users Case
This scenario reflects the case of users with low tolerance to delay, in which the MOS function in the QoE degradation range shows a great slope. Thus, in this setting the error of the MOS approximation (15) is high since the gradient of MOS among consecutive delay levels becomes significant. In Figure 7 we compare the results of this second scenario and the basic scenario. From the illustrated graph we could say that our QoEaware proposal adapts to any MOS function. This is, whereas for the rest of disciplines (excluding GM) the negative effects due to the abrupt degradation of the QoE utility function are appreciable, ASPIM manages to avoid this, adapting to any MOS function.
5.3. Scenario 3: Higher Error Case due to Rate
During the simplification process carried out for obtaining the ASPIM index rule in Section 4.1, the approximation used (see (14)) introduces a rate dependant error. Note that increasing the rate in the good channel state makes bigger the error of the approximation used. Therefore, we will analyze the obtained performance for our new QoEaware proposal under this higher error due to rate. In the previous scenarios we have assumed that the transmission rate in the good channel condition is twice the rate in the bad channel state. In this setting we consider that the transmission rate in the good channel is four times the transmission rate in the bad channel. This way, according to the QoE results presented in Figure 8, in spite of a higher error caused by rate, ASPIM is superior to the rest of policies.
5.4. Scenario 4: High Error Case due to Size
Besides, it is known that for a distribution with a decreasing hazard rate such as Pareto that the gradient of completion probability increases as long as the mean size decreases, which makes the error of the approximation (14) for the ASPIM rule achievement bigger. In this way, in this scenario we verify that ASPIM gives satisfactory behavior under a more notable error due to size, where we consider a mean size ten times lower. We conclude from Figure 9 that even under a high error caused by size that ASPIM shows good performance, being the best alternative.
5.5. Scenario 5: A Real Network Case
This last scenario reflects a real 4G wireless network context, but with the simplification that a single user transmits in each TTI. To that end, we use CQI traces obtained from a systemlevel radio access simulator [28], considering users moving at 5 km/h in the cell. For this setting with dynamic users, we define channel propagation by the Extended Pedestrian A model. As in current 4G networks, we consider 16 channel conditions. The state probabilities of these channel states are collected in Table 3, which are decreasing with the improvement of channel quality. In this scenario we consider two classes of users, which differ in size, QoE and cost. The mean size in class 2 is ten times bigger than in class 1, and we give priority to the smallsized class by having bigger cost and QoE requirement.

Under this mixture of classes and realistic channel models, in Figure 10(a) we observe that scheduling strategies show the same tendency as in the previous scenarios. It is remarkable that the difference between ASPIM and GM is bigger, especially for high values. Moreover, under the highest network load ASPIM achieves MOS values above QoE satisfaction threshold (). Furthermore, average MOS values for MR and PF dramatically decrease comparing with the previous settings.
(a)
(b)
(c)
Concerning class behaviour, we show perclass results in Figures 10(b) and 10(c). As can be concluded from the presented results, QoEaware policies guarantee fairness among classes, and ASPIM is the best choice also inside class. Aside from that, for the cost dependant cAS discipline the deterioration of the class with no priority is significant.
6. Conclusions
This paper goes into detail about obtaining a scheduling algorithm aimed at maximizing users’ perception of quality in channels with timevarying capacity. Furthermore, this work considers traffic flows with nonexponential size distributions, contrary to previous approaches that assume memoryless distributions.
As first contribution, we provide a MDPbased model for the scheduling problem presented. This model combines QoEawareness, sizeawareness and channelawareness. Nevertheless, since this model can not be solved either analytically or computationally, from the proposed MDP model we derive a simple, tractable and implementable wellperforming scheduling index rule by applying a methodology based on Whittle approach.
The proposed scheduling solution still combines QoE, size and channelawareness. Our scheduling proposal gives priority to users in their best channel condition that are out of the service unavoidable range. As verified from simulation results, the proposed scheduling strategy shows suitable subjective quality performance in a wide range of scenarios, including the case of a simplified 4G network with heterogeneous users.
Therefore, the results of this work present a relevant mathematical basis for developing scheduling algorithms aimed at maximizing QoE in current and future networks with timevarying channel capacity. This way, the proposed channelaware ASPIM discipline will be useful for network operators in order to guarantee to their customers service satisfaction in timevarying wireless networks. Apart from that, the approximations and simplifications used in Whittle method are applicable to any area when the gradient of system components is considerably small among consecutive decision slots.
As future research, we will extend our work to a multiuser approach, considering the simultaneous transmission of multiple flows per TTI.
Appendix
Proof of Proposition 4
Lemma A.1. Suppose that states are passive and that states are active, and Along with this, if we assume that and for any value, for the undiscounted case:
Proof. Referring to work elements, using (12):
And by (A.3) expressions,
This way, isolating ,
Analogously, for reward elements, using (11):
And thus, from (A.6),
Isolating ,
Substituting (A.5) and (A.8) in (A.1), and by simplifying, we obtain
And therefore, for the undiscounted case, the limit for expression (A.9) is (A.2).
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
Acknowledgment
The research leading to these results has received funding from the European Union Seventh Framework Programme (FP7/2007–2013) under Grant agreement 284863 (FP7 SEC GERYON).
References
 ITUT, P.10/G.100, Amendment 1 (01/07): New Appendix Idefinition of Quality of Experience (QoE), 2006.
 S. Aalto and P. Lassila, “Flowlevel stability and performance of channelaware prioritybased schedulers,” in Proceeedings of the 6th Euro NF Conference on Next Generation Internet (NGI '10), IEEE, June 2010. View at: Publisher Site  Google Scholar
 L. Kleinrock, Queueing Systems, vol. 2 of Computer Applications, John Wiley & Sons, 1976.
 J. Gittins, “Bandit processes and dynamic allocation indices,” Journal of the Royal Statistical Society B: Methodological, vol. 41, no. 2, pp. 148–177, 1979. View at: Google Scholar  MathSciNet
 K. Avrachenkovt, U. Ayesta, P. Brown, and E. Nyberg, “Differentiation between short and long TCP ows: predictability of the response time,” in Proceedings of the 23rd Annual Joint Conference of the IEEE Computer and Communications Societies (INFOCOM '04), vol. 2, p. 762, IEEE, 2004. View at: Google Scholar
 U. Ayesta, M. Erausquin, and P. Jacko, “A modeling framework for optimizing the flowlevel scheduling with timevarying channels,” Performance Evaluation, vol. 67, no. 11, pp. 1014–1029, 2010. View at: Publisher Site  Google Scholar
 S. Aalto, U. Ayesta, and R. Righter, “Properties of the Gittins index with application to optimal scheduling,” Probability in the Engineering and Informational Sciences, vol. 25, no. 3, pp. 269–288, 2011. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 M. Puterman, “Markov decision processes,” in Handbooks in Operations Research and Management Science, vol. 2, pp. 331–434, 1990. View at: Google Scholar
 C. H. Papadimitriou and J. N. Tsitsiklis, “The complexity of Markov decision processes,” Mathematics of Operations Research, vol. 12, no. 3, pp. 441–450, 1987. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 J. NiñoMora, “Dynamic priority allocation via restless bandit marginal productivity indices,” TOP, vol. 15, no. 2, pp. 161–198, 2007. View at: Publisher Site  Google Scholar  MathSciNet
 P. Whittle, “Restless bandits: activity allocation in a changing world,” Journal of Applied Probability, vol. 25, pp. 287–298, 1988. View at: Google Scholar
 U. Ayesta, M. Erausquin, M. Jonckheere, and I. M. Verloop, “Scheduling in a random environment: stability and asymptotic optimality,” IEEE/ACM Transactions on Networking, vol. 21, no. 1, pp. 258–271, 2013. View at: Publisher Site  Google Scholar
 S. Khan, S. Thakolsri, E. Steinbach, and W. Kellerer, “QoEbased crosslayer optimization for wireless multiuser systems,” in Proceedings of the 18th ITC Specialist Seminar on Quality of Experience, 2008. View at: Google Scholar
 S. Thakolsri, S. Khan, E. Steinbach, and W. Kellerer, “QoEdriven crosslayer optimization for high speed downlink packet access,” Journal of Communications, vol. 4, no. 9, pp. 669–680, 2009. View at: Publisher Site  Google Scholar
 P. Ameigeiras, J. J. RamosMunoz, J. NavarroOrtiz, P. Mogensen, and J. M. LopezSoler, “QoE oriented crosslayer design of a resource allocation algorithm in beyond 3G systems,” Computer Communications, vol. 33, no. 5, pp. 571–582, 2010. View at: Publisher Site  Google Scholar
 I. Taboada, F. Liberal, J. O. Fajardo, and U. Ayesta, “QoEaware optimization of multimedia flow scheduling,” Computer Communications, vol. 36, no. 1516, pp. 1629–1638, 2013. View at: Publisher Site  Google Scholar
 U. Ayesta, M. Erausquin, and P. Jacko, “Resourcesharing in a single server with timevarying capacity,” in Proceeding of the 49th Annual Allerton Conference on Communication, Control, and Computing (Allerton ’11), pp. 377–383, Monticello, Ill, USA, September 2011. View at: Publisher Site  Google Scholar
 K. Thompson, G. J. Miller, and R. Wilder, “Widearea internet traffic patterns and characteristics,” IEEE Network, vol. 11, no. 6, pp. 10–23, 1997. View at: Publisher Site  Google Scholar
 M. Fiedler, T. Hossfeld, and P. TranGia, “A generic quantitative relationship between quality of experience and quality of service,” IEEE Network, vol. 24, no. 2, pp. 36–41, 2010. View at: Publisher Site  Google Scholar
 A. Richards, M. Antoniades, V. Witana, and G. Rogers, “Mapping user level QoS from a single parameter,” in Proceedings of the International Conference on Multimedia Networks and Services, Citeseer, 1998. View at: Google Scholar
 E. Ibarrola, F. Liberal, I. Taboada, and R. Ortega, “Web QoE evaluation in multiagent networks: validation of ITUT G.1030,” in Proceedings of the 5th International Conference on Autonomic and Autonomous Systems (ICAS '09), pp. 289–294, Valencia, Spain, April 2009. View at: Publisher Site  Google Scholar
 S. Aalto and U. Ayesta, “Optimal scheduling of jobs with a DHR tail in the M/G/1 queue,” in Proceedings of the 3rd International Conference on Performance Evaluation Methodologies and Tools (ValueTools '08), article 50, Institute for Computer Sciences, SocialInformatics and Telecommunications Engineering, Athens, Greece, 2008. View at: Publisher Site  Google Scholar
 V. Visweswaran, “Decomposition techniques for MILP: lagrangian relaxation,” in Encyclopedia of Optimization, pp. 632–638, Springer, New York, NY, USA, 2009. View at: Publisher Site  Google Scholar
 I. Taboada, J. O. Fajardo, and F. Liberal, “Performance analysis of scheduling algorithms for web QoE optimization in wireless networks,” Network Protocols and Algorithms, vol. 4, pp. 27–43, 2012. View at: Google Scholar
 J. A. van Mieghem, “Dynamic scheduling with convex delay costs: the generalized $c\mu $ rule,” The Annals of Applied Probability, vol. 5, no. 3, pp. 809–833, 1995. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 S. Sesia, I. Toufik, and M. Baker, LTEThe UMTS Long Term Evolution: From Theory to Practice, Wiley, 2011.
 P. Reichl, B. Tuffin, and R. Schatz, “Logarithmic laws in service quality perception: where microeconomics meets psychophysics and quality of experience,” Telecommunication Systems, vol. 52, no. 2, pp. 587–600, 2013. View at: Publisher Site  Google Scholar
 J. C. Ikuno, M. Wrulich, and M. Rupp, “System level simulation of LTE networks,” in Proceedings of the IEEE 71st Vehicular Technology Conference, Taipei, Taiwan, May 2010. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2014 Ianire Taboada and Fidel Liberal. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.