Abstract

We consider opportunistic access to spectrum resources in cognitive wireless networks. The users equipment, or the network nodes in general are able to sense the spectrum and adopt a subset of available resources (the spectrum and the power) individually and independently in a distributed manner, that is, based on their local channel quality information and not knowing the Channel State Information (CSI) of the other nodes' links in the considered network area. In such a network scenery, the competition of nodes for available resources is observed, which can be modeled as a game. To obtain spectrally efficient and fair spectrum allocation in this competitive environment with the nodes having no information on the other players, taxation of resources is applied to coerce desired behavior of the competitors. In the paper, we present mathematical formulation of the problem of finding the optimal taxation rate (common for all nodes) and propose a reduced-complexity algorithm for this optimization. Simulation results for these derived optimal values in various scenarios are also provided.

1. Introduction

Opportunistic spectrum access and flexible and efficient spectrum allocation procedures as well are considered as measures to increase the utilization of the scarce radio resources in future wireless communication networks. Apart from the spectral efficiency, the Quality of Experience (QoE), and the associated fairness in resources distribution are in the focus of research towards the cognitive, opportunistic, and dynamic spectrum access. The spectrum allocation procedures are usually centralized, require the Channel State Information (CSI) of all links in the network, and involve the overhead traffic, which in turn occupies the scarce radio resources. For the future communication concepts, such as cognitive or opportunistic radio, the nodes are expected to take intelligent decisions on the amount of resources to be utilized in a distributed way, thus minimizing or eliminating the overhead traffic.

In this paper, we consider opportunistic acquisition of orthogonal frequency channels by the network nodes. An example of the multiple-access technique using such orthogonal channels is the well-known Orthogonal Frequency Division Multiple Access (OFDMA). In the opportunistic OFDMA, the network nodes are able to adopt a subset of accessible subcarriers (SCs) individually, as well as the transmission rate and power allocated to these SCs [1]. Below, we consider a more general scenario of the opportunistic access to frequency channels of any bandwidth, limited centralized management, and very limited control traffic, that is, there is no central frequency-channel scheduler, and no CSI exchange between the network nodes. Our approach to opportunistic spectrum allocation is related to noncooperative game theory and to the concept of pricing.

The game-theoretic scheduling for OFDMA has been considered in the literature as centralized and distributed SCs allocation. The centralized schemes allow for more efficient and fair spectrum utilization; however they require centralized management and a considerable amount of control traffic related to the CSI of all possible links in all considered frequency channels and to the information on the allocated channels. This information has to be exchanged or to be available at a central unit (e.g., at a base station of a cellular network) every time the channels qualities change for the nodes in the network area. Newest results for such centralized solutions based on cooperative complete information game models have been presented in [2–4]. Distributed decision making, on the contrary, deploys noncooperative games and seeks for Nash Equilibrium (NE) as a game solution. However, for the spectrum allocation, only the complete-information games have been considered in the literature so far. We believe that such models cannot be considered for practical applications in dynamically changing wireless networks, since the complete knowledge of the CSI related to all links to be available at every other node would require a lot of control traffic between the nodes in dedicated control channels. This information would have to be sent every time the channels change, so in mobile environment, the control traffic would be comparable to the information-data traffic. Thus, noncooperative complete-information game models are only suitable for multicell environment, where the players are the base stations, which have the CSI of all links in their cell areas [5, 6], or in static wireless scenarios.

The concept of resource pricing (or coercive taxation) has been considered in the literature extensively for power allocation, for example, for OFDM and OFDMA in [7–9]. There, the resource that is taxed is the power used by the network nodes, and the goal is to maximize the sum throughput given the total allowable transmission power in the network. In these papers however, it has been assumed that the SCs distribution among the users has already been done somehow centrally.

Pricing has been also applied for distributed power and interference management in the network, for example, in [10] for the code-division multiple access. For this purpose the complete-information noncooperative game models have been formulated. Such game-theoretic problem formulation has practical application for interference management, due to the fact that the complete information on the interference level can be available for each player, since the nodes can measure it locally. However, for the spectrum management the problem is different.

Some papers, for example, [11–13] consider distributed allocation of resources based on pricing in a multicell scenario, where the base stations act as players. The pricing concepts developed for the multicell scenario, cannot be considered for the distributed resource allocation in decentralized networks, because contrary to a base station that may have the CSI of all links, a network node may only have the CSI of its own link. Another approach to price based spectrum management is based on iterative water filling, which allows all players to use the same frequency channels and adjust their power levels in these channels based on pricing function [14]. Definition of the pricing function for each player requires the CSI knowledge and exchange between the neighboring players, as well as the number of iterations. Similarly, in [15] the information exchange between the secondary and primary users is assumed for the spectrum leasing. In [16] the spectrum sharing is modeled as the oligopoly competition between the primary users and the Bertrand game model, which again requires the knowledge of the secondary users’ CSI by the players (primary users). Although the above-mentioned works have contributed to significant advance in the game-theoretic price-based models for spectrum sharing, they all make an assumption on the complete information available for all players that relate to their links CSI or they narrow to the power and interference management.

In our earlier work [17], we have presented distributed SCs allocation interference-free method in a network of the OFDMA-based opportunistic radios. It aimed at rational and efficient spectrum utilization in both the uplink or downlink transmission. Rationality in our case means that apart from maximizing the spectral efficiency, the network and each individual node aim at lowering the cost of this efficiency (resulting from taxation of resources) and at increasing the QoE (resulting from the number of served nodes). These rationality measures have been reflected in the definition of the noncooperative game model with complete information, and in the utility function defined for each player (the network node). The definition of this game involved aggregation of the players, in such a way that each player (the network node) can view all other players as one, named the network-nodes community (NNC). The complete information required in this game does not include the individual CSI of the other network nodes, but only the local (single-link) CSI and the taxation rate. This way noncooperative game with full information is reasonable and practically applicable in the dynamically changing network scenarios.

Here, in this paper, we present a generalized framework for the taxation-based allocation of orthogonal frequency channels in the opportunistic radio. First, we show inappropriateness of the complete-information game models in the considered framework. Then, we consider selfish and social behavior of the players by appropriate definitions of the utility functions reflecting such behaviors. These utility functions include the linear-taxation summand dependent on the amount of the acquired spectrum resources. We aim at finding the optimal taxation rate to come up with a high overall network efficiency defined in two ways. We provide the mathematical description of the problem of finding the optimal tax rate, show that the problem complexity is NP hard, and present and examine the reduced-complexity algorithm for solving it.

In Section 2, we present the main idea of the proposed game-theoretic approach to distributed spectrum allocation and provide formal definitions of the considered games. In Section 3, we mathematically derive the amount of bandwidth each player is inclined to acquire. In Section 4, we present the reduced-complexity algorithm to obtain the optimal taxation, where the optimality is defined in a number of ways. The simulation results are presented in Section 5, and the work is concluded in Section 6.

2. Taxation-Based Models of Distributed Spectrum Allocation

We consider the scenery of multiple cognitive-radio nodes (or users) appearing in the opportunistic network area, that make use of the orthogonal frequency channels, for example, OFDMA subcarriers. It implies that the nodes do not have to apply any guard frequency bands to limit the out-of-band interference. The frequency correction at the receiver is also assumed to be perfect. This scenery of the opportunistic and cognitive radio network is presented in Figure 1. The nodes are able to sense the radio environment, detect the parameter called tax rate available in a given area, detect available spectrum resources, and acquire a subset of these resources usable for their intended transmission, for example, peer-to-peer communication, an access to a cellular network or to any wireless network in general. The goal of each node is to make the best use of these resources, that is, obtain high data rate at the lowest cost. As a proof of our concept, we consider the freedom in the spectrum allocation, that is, theoretically even the smallest part of the spectrum can be used by a player. This theoretical assumption can be refined for practical applications, if we assume that the nodes demand the spectrum only if they can make use of it, that is, if there is a minimal contiguous part of the spectrum available for their intended communication (such as one OFDM subcarrier band) that may include the protection band to mitigate the interference to other transmissions.

Let us consider the resource acquisition procedure as a game, which each network node plays against the other nodes (the players). Let us assume that there are 𝐾 players, and the available bandwidth is 𝐡. (For the simplicity of our considerations, in the remainder of this paper, we assume that 𝐡 and 𝐾 are fixed, although in dynamically changing network environment, the number of players, their demands, and the available bandwidth vary.) A single player decides what portion of the available bandwidth she is going to use. (Note that such a personification and female pronouns are established in the game-theoretic convention.) Selfish player aiming at her throughput maximization would occupy the whole available spectrum; however, such a behavior decreases the spectral efficiency and the capacity of the network, as well as the QoE of other players who cannot access the network. The problem is known as the Tragedy of Commons described in [18].

As a countermeasure for the problem of commons sharing and utilization, taxation of the resources is introduced. In our network scenery presented in Figure 1, the taxation rate is the same for all network nodes or users (the players). This tax rate is known in the considered area. It can be stored in an area database (among other parameters, required for the efficient operation of the cognitive or opportunistic users, e.g., the spectrum masks for available bands in a given location and time, primary-users detection thresholds, etc.) or transmitted by the elected master node in case, when the considered network operates independently from the area database in an ad hoc manner. It is being updated periodically and broadcasted in this area as one parameter among many other ones in a typical Broadcast Control CHannel (BCCH), or specifically defined Cognitive Pilot Channel (CPC) [19]. Let us stress that this broadcast transmission of a single parameter occupies really minor resources, contrary to the situation of transmitting full CSI of all links in the considered frequency band using dedicated channels.

2.1. Inappropriateness of the All Links Complete-Information Game Model

Let us first show that the complete-information game model which makes use of the CSI of all involved links is not suitable for our scenario. The utility function for player π‘˜ in such a game, in which the concept of resource taxation is applied, reflects her throughput (revenue) and its related cost, and is defined as (Note that from this point, the mathematical analysis in this paper is performed in a continuous space for the sake of generality, however it can be easily translated into a discrete orthogonal channel scenario.) πœξ€·π‘€π‘˜ξ€Έ=ξ€œ(𝑓)𝐡2𝐡1log2ξ€Ί1+π›Ώπ‘˜π‘ƒπ‘˜(𝑓)π›Ύπ‘˜(𝑓)π‘€π‘˜ξ€»(𝑓)π‘‘π‘“βˆ’π‘Ÿ0π‘π‘˜,(1) where π‘€π‘˜(𝑓) is the function indicating the occupancy of the frequencies by player π‘˜ (π‘€π‘˜(𝑓)=1 if frequency 𝑓 is assigned to user π‘˜ and π‘€π‘˜(𝑓)=0 otherwise), 𝐡1 and 𝐡2 are the lower and the upper bounds of the available spectrum (𝐡2βˆ’π΅1=𝐡), π‘π‘˜ is the amount of bandwidth the player acquires, π›Ύπ‘˜(𝑓)=|π»π‘˜(𝑓)|2/𝒩0 is the Carrier-to-Noise Ratio (CNR) measured at the frequency 𝑓, π»π‘˜(𝑓) and π‘ƒπ‘˜(𝑓) are the π‘˜th user’s channel characteristic and the power spectral density allocated to this frequency, respectively, and 𝒩0 is the noise power spectral density. Let us note that for the case of orthogonal channels, interference that is usually added to noise equals zero. Moreover, in (1), π›Ώπ‘˜ is the factor (often called the Signal-to-Noise Ratio (SNR)-gap) depending on the assumed player’s Bit Error Probability (BEP) π‘ƒπ‘’π‘˜. (In case of M-QAM, π›Ώπ‘˜=βˆ’1,5/ln(0,5π‘ƒπ‘’π‘˜), while for 𝑀β‰₯4 and the SNR in the range of 0–30 dB it can be set more precisely as π›Ώπ‘˜=βˆ’1,5/ln(5π‘ƒπ‘’π‘˜) [20].) Finally, in the above equation, π‘Ÿ0 is a linear tax rate.

Since every channel can be used by a sole player: βˆ€π‘“βˆƒ!π‘˜βˆΆπ‘€π‘˜π‘(𝑓)=1,π‘˜=ξ€œπ΅2𝐡1π‘€π‘˜(𝑓)𝑑𝑓.(2)

To find the NE in such a game, we shall be solving this problem numerically. This numerical representation of finding the NE is the binary linear programming problem, that is, for a given set of π›Ύπ‘˜(𝑓) values for all 𝐾 nodes and for a considered number of available channels we shall find binary values π‘€π‘˜(𝑓) for a given π‘Ÿ0. Then, looking for optimum π‘Ÿ0 to maximize some goal function, for example, the network sum throughput, would be even more complex. Note, that solution of such a defined problem would require the knowledge of all π›Ύπ‘˜(𝑓) by every node, and thus, as mentioned before such a model for resource allocation is not practical.

2.2. Game Models Using Only the Local CSI

To narrow the space of this analysis and to eliminate the necessity for complete information concerning other players’ strategies and payoffs, we propose to treat the rest of the players as a whole (the NNC). Note that NNC is not formally organized in any way. It is only viewed as such by a single player. Moreover, we let the players take decisions independently and subsequently one after another, as they appear in the network. It is a usual case in all wireless networks that some collision avoidance mechanism is implemented, or an access to decision-making entity (e.g., base-station) makes use of the random access channel to avoid taking decisions at the same time. We also limit the players in the maximum amount of bandwidth they can take at a time as a countermeasure for their greedy behavior, and let this maximum allowed bandwidth be 𝐡𝐼. The rationale behind such a limitation is proved in [21].

Let us first look at the utility defined for the π‘˜th player so, as to reflect the player’s throughput: πœπ‘˜ξ€·π‘π‘˜ξ€Έ=ξ€œπ‘“βˆˆπ’π‘˜log2ξ€Ί1+π›Ώπ‘˜π‘ƒπ‘˜(𝑓)π›Ύπ‘˜(𝑓)π‘‘π‘“βˆ’π‘Ÿ0π‘π‘˜,(3) where π’π‘˜ is the (noncountable) set of frequencies player π‘˜ occupies. Let us note that before a decision is taken by player π‘˜, she senses the available spectrum resources and knows which frequency bands are already occupied (the amount of available bandwidth for the π‘˜th player is π΅π‘˜, and there is no possibility that different players acquire the same frequencies: βˆ€π‘˜β‰ jβˆΆπ‘†π‘˜βˆ©π‘†π‘—= Ø). Thus, the game model is dynamic. Moreover, the game for each player is two-dimensional only, what noticeably simplifies the problem.

The model defined above is not really a commonly understood game, since the decision of the NNC is not represented in the utility function (unless we call it a game against the radio environment, or a game against the nature as in [22]). This model reflects selfish behavior of the players, who care only for their own throughput. In the reminder of the paper we will call this model Selfish-Behavior Model.

In [17] we have proposed a noncooperative complete-information two-dimensional game, where the payoffs of both players (player π‘˜ and the NCC at the π‘˜th game stage) reflect their decisions. In such a case, being the alternative to the framework discussed above, the strategy space of the considered π‘˜th player is the same (consists of possible number of the orthogonal channels the player acquires). The strategies of the NNC are the numbers of channels this community may occupy all together apart from the considered player. The utility function of the considered π‘˜th player is defined as follows: πœ‰π‘˜ξ€·π‘π‘˜,π‘π‘˜ξ€Έ=1π΅βˆ«π‘“βˆˆπ’π‘˜log2ξ€Ί1+π›Ώπ‘˜π‘ƒπ‘˜(𝑓)π›Ύπ‘˜ξ€»ξ‚‡β‹…ξ€½π΅(𝑓)π‘‘π‘“π‘˜βˆ’π‘π‘˜βˆ’π‘π‘˜ξ€Ύβˆ’π‘Ÿ0π‘π‘˜,(4) where π‘π‘˜ is the amount of bandwidth that will be occupied by NNC, and π΅π‘˜ is the amount of bandwidth available at the π‘˜th game stage. One may interpret formula (4) as the total normalized throughput (throughput per frequency unit of the total available bandwidth 𝐡) which could be obtained by the new incoming players in case they occupied the remaining bandwidth and had the same average spectral efficiency as the considered player. This way, in the decision making on how much of the spectrum to occupy, the players factor the social aspect of the network (to serve multiple nodes) and not just their own benefit. In the reminder of the paper the game model with the π‘˜th player’s utility function defined by (4) will be called the Social-Behavior Model.

Finally, the payoff for the NNC at the π‘˜th game stage can be defined as πœ€π‘˜=π‘π‘˜, where π‘π‘˜βˆˆ[0,max(0,π΅π‘˜βˆ’π΅πΌ)], that is, the number of resources that can be potentially occupied by NNC.

3. Optimal Choice of Spectrum Resources: The Players’ Perspective

Let us now consider how the players choose their strategic options, and how to coerce their desired behavior to obtain specifically defined network benefit.

3.1. Selfish-Behavior Model

As the practical approach to the Selfish-Behavior Model (SelBM) described by the utility function (3), we propose to eliminate dominated (disadvantageous) strategies of the players. Taking the π‘˜th player’s strongest frequencies (the ones that have the highest π›Ύπ‘˜(𝑓) values) into account is in fact the elimination of dominated strategies of the considered player. Note, that from both the individual node and the whole network perspective, making use of the strongest channels results in higher spectral efficiency. Thus, the strategies of a single node are all possible channels of the strongest frequencies (from 0 to 𝐡𝐼). In order to find the optimum taxation, let us consider rewriting formula (3): Μ‡πœπ‘˜ξ€·π‘π‘˜ξ€Έ=ξ€œπ‘π‘˜0log2ξ€Ί1+π›Ώπ‘˜Μ‡π‘ƒπ‘˜ξ€·π‘“π‘˜ξ€ΈΜ‡π›Ύπ‘˜ξ€·π‘“π‘˜ξ€Έξ€»π‘‘π‘“π‘˜βˆ’π‘Ÿ0π‘π‘˜,(5) where Μ‡π›Ύπ‘˜(π‘“π‘˜) is the function resulting from ordering (in the descending manner) of the continuous values of π›Ύπ‘˜(𝑓): (Ordering operation described by (6) is a hypothetical bijective mapping, which cannot be generically defined for any continuous space of values and depends on π›Ύπ‘˜(𝑓). This operation involves mapping of both the domain arguments and the codomain images of π›Ύπ‘˜(𝑓) to new arguments and images of Μ‡π›Ύπ‘˜(π‘“π‘˜) belonging to the same domain and codomain, respectively. For a discrete set of values ordering can be done by a standard sorting algorithm.) Μ‡π›Ύπ‘˜ξ€·π‘“π‘˜ξ€ΈβŸ·π›Ύπ‘˜(𝑓)βˆΆβˆ€π‘›,π‘šπ‘“π‘˜π‘›β‰€π‘“π‘˜π‘šβŸΉΜ‡π›Ύπ‘˜ξ€·π‘“π‘˜π‘›ξ€Έβ‰₯Μ‡π›Ύπ‘˜ξ€·π‘“π‘˜π‘šξ€Έ,(6) and Μ‡π‘ƒπ‘˜(π‘“π‘˜) is the power spectral density resulting from the optimal power allocation for a total-power constraint, that is, from the water filling. Note that this ordering takes place at every stage, so in (5), the lower integration endpoint is always equal to zero. Moreover, as mentioned before π‘π‘˜ξ€½π΅β‰€min𝐼,π΅π‘˜ξ€Ύ,(7) where 𝐡I is the maximum allowable bandwidth one player can take at a time, and π΅π‘˜ is the bandwidth available at the π‘˜th game stage (for the π‘˜th player). We also limit our choices of π‘π‘˜ to only useful frequencies, that is, π‘π‘˜β‰€π΅WFπ‘˜,(8) where 𝐡WFπ‘˜ is the useful bandwidth after the water-filling, that is, the bandwidth, in which: Μ‡π‘ƒπ‘˜(π‘“π‘˜)>0 for all π‘“π‘˜βˆˆ[0,𝐡WFπ‘˜]. In such a case: βˆ€π‘“π‘˜β‰€π‘π‘˜βˆΆΜ‡π‘ƒπ‘˜ξ€·π‘“π‘˜ξ€Έ=1𝛾0π‘˜ξ€·π‘π‘˜ξ€Έβˆ’1π›Ώπ‘˜Μ‡π›Ύπ‘˜ξ€·π‘“π‘˜ξ€Έ,(9) where 1/𝛾0π‘˜(π‘π‘˜) is the water level obtained in the water-filling algorithm over the acquired bandwidth π‘π‘˜ of player π‘˜. The utility function Μ‡πœπ‘˜(π‘π‘˜) can be thus expressed as: Μ‡πœπ‘˜ξ€·π‘π‘˜ξ€Έ=ξ€œπ‘π‘˜0log2ξ€Ίπ›Ώπ‘˜Μ‡π›Ύπ‘˜ξ€·π‘“π‘˜ξ€Έ/𝛾0π‘˜ξ€·π‘π‘˜ξ€Έξ€»π‘‘π‘“π‘˜βˆ’π‘Ÿ0π‘π‘˜.(10) It can be easily shown that the above-defined function is concave, (For all π‘π‘˜ defined by (7) and (8), the first summand is concave and monotonically increasing because π›Ώπ‘˜Μ‡π›Ύπ‘˜(π‘“π‘˜)/𝛾0π‘˜(π‘π‘˜)>1 for all π‘“π‘˜, and the second summand is linearly decreasing with π‘π‘˜) so we may find its maximum (as each rational player would do) by solving the following equation: πœ•πœ•π‘π‘˜Μ‡πœπ‘˜ξ€·π‘π‘˜ξ€Έ=0.(11) As derived in Appendix A (formula (A.10)): πœ•πœ•π‘π‘˜Μ‡πœπ‘˜ξ€·π‘π‘˜ξ€Έ=β„±π‘˜ξ€·π‘π‘˜ξ€Έβˆ’π‘Ÿ0,(12) where β„±π‘˜(π‘π‘˜) can be defined in a number of ways depending on the player’s CNR characteristic π›Ύπ‘˜(𝑓), and its resulting sorted values at the π‘˜th game stage Μ‡π›Ύπ‘˜(π‘“π‘˜). As derived in Appendix A in (A.11) for the two-path propagation model it can be approximated as β„±ξ€·π‘π‘˜ξ€Έβ‰ƒlog2ξ€½π΄π‘˜π›Ώπ‘˜ξ€Ίξ€·1+cos2πœ‹π‘π‘˜πœπ‘˜βˆ’πœ‘π‘˜ξ€Έξ€»ξ€Ύ+log2ξ‚»π‘π‘˜βˆ’1ξ‚ΈΞ“π‘˜+12πœ‹πœπ‘˜π΄π‘˜π›Ώπ‘˜ξ‚ƒξ‚€tgπœ‹π‘π‘˜πœπ‘˜βˆ’πœ‘π‘˜2ξ‚ξ‚€πœ‘βˆ’tgπ‘˜2+π‘ξ‚ξ‚„ξ‚„ξ‚‡π‘˜Γ—ξ‚»Ξ“ln2π‘˜+12πœ‹πœπ‘˜π΄π‘˜π›Ώπ‘˜ξ‚ƒξ‚€tgπœ‹π‘π‘˜πœπ‘˜βˆ’πœ‘π‘˜2ξ‚ξ‚€πœ‘βˆ’tgπ‘˜2ξ‚Όξ‚ξ‚„βˆ’1⋅2π΄π‘˜π›Ώπ‘˜cos2ξ‚€πœ‹π‘π‘˜πœπ‘˜βˆ’πœ‘π‘˜2ξ‚ξ‚‡βˆ’1,(13) where the π΄π‘˜, πœ‘π‘˜, and πœπ‘˜ are the parameters of the considered multipath propagation model depending on the signal attenuation, average phase difference between the arriving multipath signal components, and the multipath delay spread (see Appendix A). Thus, by solving (11) we obtain the amount of bandwidth π‘βˆ—π‘˜ the π‘˜th player is inclined to acquire π‘π‘˜=π‘βˆ—π‘˜ξ€·π‘βˆΆβ„±βˆ—π‘˜ξ€Έ=π‘Ÿ0.(14) Note, that for a given π‘Ÿ0, user π‘˜ can find π‘βˆ—π‘˜, not knowing other players CSI, and this finding is independent from other players’ choices.

3.2. Social-Behavior Model

Let us now consider formula (4) reflecting the Social Behavior Model (SocBM) in the form with ordered values of π›Ύπ‘˜(𝑓) (similarly as in the previous section): Μ‡πœ‰π‘˜ξ€·π‘π‘˜,π‘π‘˜ξ€Έ=1π΅ξ‚»ξ€œπ‘π‘˜0log2ξ€Ί1+π›Ώπ‘˜Μ‡π‘ƒπ‘˜ξ€·π‘“π‘˜ξ€ΈΜ‡π›Ύπ‘˜ξ€·π‘“π‘˜ξ€Έξ€»π‘‘π‘“π‘˜ξ‚Όβ‹…ξ€½π΅π‘˜βˆ’π‘π‘˜βˆ’π‘π‘˜ξ€Ύβˆ’π‘Ÿ0π‘π‘˜,(15) Based on (6)–(9), the above formula can be easily (again similarly as in the previous section) converted to Μ‡πœ‰π‘˜ξ€·π‘π‘˜,π‘π‘˜ξ€Έ=1π΅ξ‚»ξ€œπ‘π‘˜0log2ξ€Ίπ›Ώπ‘˜Μ‡π›Ύπ‘˜ξ€·π‘“π‘˜ξ€Έ/𝛾0π‘˜ξ€·π‘π‘˜ξ€Έξ€»π‘‘π‘“π‘˜ξ‚Όβ‹…ξ€½π΅π‘˜βˆ’π‘π‘˜βˆ’π‘π‘˜ξ€Ύβˆ’π‘Ÿ0π‘π‘˜.(16) It can be shown that function (16) is concave, (For all π‘π‘˜ defined by (7) and (8), the first factor in the first summand is concave and monotonically increasing, the second factor in the first summand is linearly decreasing with π‘π‘˜, and so is the second summand.) so we may find its maximum (as each rational player would do) by solving the following equation: πœ•πœ•π‘π‘˜Μ‡πœ‰π‘˜ξ€·π‘π‘˜,π‘π‘˜ξ€Έ=0.(17) The derivative (πœ•/πœ•π‘π‘˜)Μ‡πœ‰π‘˜(π‘π‘˜,π‘π‘˜) is defined by formula (B.1) obtained in Appendix B, whose simplified form is the following: πœ•πœ•π‘π‘˜Μ‡πœ‰π‘˜ξ€·π‘π‘˜,π‘π‘˜ξ€Έ=π’’π‘˜ξ€·π‘π‘˜,π‘π‘˜ξ€Έβˆ’π‘Ÿ0,(18) where π’’π‘˜(π‘π‘˜,π‘π‘˜) can be defined as in (B.3) for rural areas. (We do not repeat its long formula here. See Appendix B for its definition.) Thus, by solving (16) we obtain the amount of bandwidth π‘βˆ—π‘˜ which the π‘˜th player is inclined to acquire taking into account the considered amount of bandwidth to be occupied by the NNC at the π‘˜th game stage: π‘π‘˜ξ€·π‘π‘˜ξ€Έ=π‘βˆ—π‘˜ξ€·π‘π‘˜ξ€Έξ€·π‘βˆΆπ’’βˆ—π‘˜ξ€Έ=π‘Ÿ0.(19) In other words, π‘βˆ—π‘˜(π‘π‘˜) is the best-response function in the considered two-dimensional game.

Let us recall that the payoff of the NNC is defined as πœ€π‘˜=π‘π‘˜, where π‘π‘˜βˆˆ[0,max(0,π΅π‘˜βˆ’B𝐼)] and is not dependent on the strategy of the π‘˜th player π‘π‘˜. Thus, the NNC would always choose the strategy π‘βˆ—π‘˜=max(0,π΅π‘˜βˆ’π΅πΌ) resulting in its highest possible payoff. For this NNC strategy (π‘βˆ—π‘˜) and for π‘βˆ—π‘˜(π‘βˆ—π‘˜) we obtain the NE. Thus, for the calculated equilibrium strategies the players acquire a portion of bandwidth for their transmission.

4. Optimal Tax Rates

To obtain the desired behavior of the players and high overall network efficiency, the tax rate for the considered games (presenting either the SelBM or SocMB) should be properly chosen to obtain the maximum benefit for the whole network in the considered framework, for example, the maximum sum throughput reflecting the efficiency of the spectrum distribution. We can define our objective function as πœ‚ST=1π΅πΎβˆ’1ξ“π‘˜=0ξ€œπ‘βˆ—π‘˜0log2ξ€Ίπ›Ώπ‘˜Μ‡π›Ύπ‘˜ξ€·π‘“π‘˜ξ€Έ/𝛾0π‘˜ξ€·π‘βˆ—π‘˜ξ€Έξ€»π‘‘π‘“π‘˜,(20) which is the sum throughput (ST) of all players averaged over the total available bandwidth 𝐡. Alternatively, we may look at maximizing the actual spectral efficiency (SE) of the transmission in the network (the sum throughput averaged over the actually used frequency bandwidth): πœ‚SE=1βˆ‘πΎβˆ’1π‘˜=0π‘βˆ—π‘˜πΎβˆ’1ξ“π‘˜=0ξ€œπ‘βˆ—π‘˜0log2ξ€Ίπ›Ώπ‘˜Μ‡π›Ύπ‘˜ξ€·π‘“π‘˜ξ€Έ/𝛾0π‘˜ξ€·π‘βˆ—π‘˜ξ€Έξ€»π‘‘π‘“π‘˜.(21) The next step would be to find the optimum value of π‘Ÿ0 to maximize either function (20) or (21). Note that many other definitions of the objective function are possible, that could reflect the fairness or proportional fairness in the distribution of resources, as well as other factors, for example, the percentage of used bandwidth or the percentage of served players. Below, we will examine the two objective functions defined above by (20) and (21) and show that some fairness in resource distribution is also achieved with a tax-rate optimal for (20).

The values of πœ‚ST or πœ‚SE depend on π‘Ÿ0, and on the values of Μ‡π›Ύπ‘˜(π‘“π‘˜). (This is because π‘Ÿ0 and Μ‡π›Ύπ‘˜(π‘“π‘˜) have the implication on π‘βˆ—π‘˜ and on the throughput obtained by player π‘˜.) Moreover, the order of appearance of the players in the game matters, since Μ‡π›Ύπ‘˜(π‘“π‘˜) depends on ̇𝛾𝑗(𝑓𝑗) for all 𝑗<π‘˜ (frequencies allocated to players taking decisions before player π‘˜ must be excluded from this player considerations). Unfortunately, both πœ‚ST and πœ‚SE are neither strictly concave or convex functions of π‘Ÿ0. In general, for low π‘Ÿ0 it pays off for all players to acquire the highest possible amount of bandwidth, irrespective to their channel qualities. As π‘Ÿ0 increases, it becomes affordable to acquire some bandwidth only for the players with good channel conditions (high π›Ύπ‘˜(𝑓) and Μ‡π›Ύπ‘˜(π‘“π‘˜)). Thus, in such a case, the spectral efficiency πœ‚SE is increased, and the average sum-throughput πœ‚ST may be decreased. However, when π‘Ÿ0 is too high, that is, close to the barrage tax rate, only very few players can afford some small portion of the frequency band, so a lot of available bandwidth is not used, and thus both the average sum-throughput πœ‚ST and the spectral efficiency of the network πœ‚SE are low. Thus, there exist some optimum values for π‘Ÿ0 to maximize either πœ‚ST or πœ‚SE: π‘Ÿ0=π‘Ÿβˆ—ST0βˆΆπœ‚STξ‚€π‘Ÿβˆ—ST0=maxπœ‚STξ€·π‘Ÿ0ξ€Έ,π‘Ÿ0=π‘Ÿβˆ—SE0βˆΆπœ‚SEξ‚€π‘Ÿβˆ—SE0=maxπœ‚SEξ€·π‘Ÿ0ξ€Έ.(22) However, it is not straightforward to determine these optimum values. (As mentioned before, πœ‚ST and πœ‚ST depend on Μ‡π›Ύπ‘˜(π‘“π‘˜) functions, which have different arguments π‘“π‘˜ for different players, and on the order of players’ appearance.) To find this optimum, even numerically, is a complex NP-hard problem, and the optimization procedure has to be performed every time the users’ channels change. Some simplifications can be obtained in finding this optimum in a proper time span (not necessarily shorter than the coherence time of the tracked processes: π›Ύπ‘˜(𝑓)), because the value of the optimum in the next time instant should be found close to the optimum value in the previous time instant. For this purpose we may apply a method that actually traces the variations of the objective function (πœ‚ST or πœ‚SE) rather than the variations of the players’ channel conditions.

The considerations presented in this section for continuous orthogonal channels can be easily translated to discrete orthogonal channels, that is, to the case of having a set of 𝑁 available channels (e.g., OFDM subcarriers) to be acquired by the players. In such a case, the integrations in (5), (10), (15), (16), (20), and (21) should be replaced by summations, and the value of π‘π‘˜ and π‘βˆ—π‘˜ should be approximated by the discrete number of channels π‘–π‘˜ of a particular bandwidth Δ𝑓 (with a particular resolution). Moreover, as shown in [17], there exists the NE for the discrete orthogonal channels (like in OFDMA). Below, for such a case of the available bandwidth discretization, we present the optimal tax-rate-searching algorithm with reduced complexity tracing the instantaneous variations of the network-objective function πœ‚ST around its maximal value.

Step 1. Initialize algorithm:
Determine 𝑁 available channels of Δ𝑓 bandwidth,
Determine the range of π‘Ÿ0∈[π‘Ÿmin,π‘Ÿmax],
Determine the increment of π‘Ÿ0: Ξ”π‘Ÿ0,
Determine acceptable value of πœ‚ST: πœ‚STmin,
Determine 𝐾, 𝐡𝐼, measure the players’ π›Ύπ‘˜(𝑓),
Determine the order of players’ appearance.

Step 2. Find the π‘βˆ—π‘˜ values for all considered tax-rates π‘Ÿ0.

Step 3. Calculate πœ‚ST for all r0.

Step 4. Find optimal tax-rate π‘Ÿβˆ—ST0 that maximizes πœ‚ST.

Loop 1:
 Step L1.1. Monitor the network-objective function πœ‚ST. Step L1.2. If πœ‚ST<πœ‚STmin, go to Loop 2, else go to Loop 1.

Step 5. Update (increase) π‘Ÿβˆ—ST0: π‘Ÿβˆ—ST+∢=π‘Ÿβˆ—ST0+Ξ”π‘Ÿ0,

Step 6. Calculate the resulting πœ‚ST+=πœ‚ST(π‘Ÿβˆ—ST+),

Step 7. If πœ‚ST+>πœ‚ST, go to Step 9
else go to Loop 3.

Loop 2:
Repeat Step L2.1. π‘Ÿβˆ—ST0∢=π‘Ÿβˆ—ST+, Step L2.2. πœ‚ST∢=πœ‚ST+, Step L2.3. π‘Ÿβˆ—ST+∢=π‘Ÿβˆ—ST0+Ξ”π‘Ÿ0, Step L2.4. Calculate πœ‚ST+=πœ‚ST(π‘Ÿβˆ—ST+), Until πœ‚ST+>πœ‚ST.

Step 8. Go to Step 12.

Step 9. Update (decrease) π‘Ÿβˆ—ST0: π‘Ÿβˆ—STβˆ’βˆΆ=π‘Ÿβˆ—ST0βˆ’Ξ”π‘Ÿ0,

Step 10. Calculate the resulting πœ‚STβˆ’=πœ‚ST(π‘Ÿβˆ—STβˆ’),

Step 11. If πœ‚STβˆ’>πœ‚ST, go to Loop 3
else go to Step 12.

Loop 3:
Repeat Step L3.1. π‘Ÿβˆ—ST0∢=π‘Ÿβˆ—STβˆ’, Step L3.2. πœ‚ST∢=πœ‚STβˆ’,  Step L3.3. π‘Ÿβˆ—STβˆ’βˆΆ=π‘Ÿβˆ—ST0βˆ’Ξ”π‘Ÿ0, Step L3.4. Calculate πœ‚STβˆ’=πœ‚ST(π‘Ÿβˆ—STβˆ’), Until πœ‚STβˆ’>πœ‚ST.

Step 12. If πœ‚STβ‰₯πœ‚STmin
communicate new tax-rate and go to Loop 1,
else Warn: β€œNo tax-rate meeting the objectives”

As it will be shown in the next section, πœ‚ST is more appropriate as the network-objective function for the fairness of resource distribution among the players in the case of both SelBM and SocBM, and always has a maximum when 𝐡I is properly chosen (not too small) for a given 𝐾 and 𝐡. Analogous algorithm to the one presented below can be performed for searching the maximum of πœ‚SE. The presented algorithm has reduced complexity due to the application of the following methods: optimum tax-rate searching around the previous optimum and optimization procedure running only when πœ‚ST drops below the required value: πœ‚STmin.

Alternatively, to reduce the rate of necessary calculations to solve the optimization problem, we may maximize the expected values of (20) or (21) over the set of random variables Μ‡π›Ύπ‘˜(π‘“π‘˜): πœ‚ST=π”ΌΜ‡π›Ύπ‘˜{πœ‚ST}, πœ‚SE=π”ΌΜ‡π›Ύπ‘˜{πœ‚SE}. Such definitions of the objective functions could be useful if we were able to approximate the expectation values with the average values and use them in a static (or slowly changing) environment. The resulting tax rate would approximate the optimum one (either π‘Ÿβˆ—ST0 or π‘Ÿβˆ—SE0) with unknown accuracy, while the optimization procedure can be performed off-line. This option is to be investigated in the future.

5. Numerical Results

Our simulation setup is the following. We have considered an available bandwidth 𝐡 with the resolution Δ𝑓=𝐡/256, where Δ𝑓 can be considered as the smallest spectrum unit, that can be occupied by orthogonal signals, for example, OFDM subcarriers. In our considered scenario, the total transmission power has been fixed. The power constraint for each link results from the distance between the transmitter and the receiver and from the power-control mechanism. (Usually this mechanism is applied to combat the near-far effect and the interference between the users; however, here, we assume orthogonal frequency channels, so this mechanism is only used to assure the appropriate quality of the link, i.e., the required average SNR, which in our case has been set to 30 dB). For our simulation purposes, it has been assumed that the order of appearance of the players in a game is random. Furthermore, we assume that the power control mechanism has a tolerance of 3 dB, so that random deviation from the average SNR is possible for any node (average SNR =(30Β±3) dB). This average SNR deviation (which also reflects the accuracy of power-control in modern radio systems) has been chosen to differentiate possible link qualities. Moreover, two example channel models have been compared. The first one is the two-path Rayleigh-fading channel with the delay spread ranging from 0 to 1/64 of 1/𝐡, and the average power of the second path being βˆ’3 dB relative to the first path (such a model can be considered as suitable for rural environment). The second considered model is the six-path channel, with paths having the same power, and delays uniformly spread between 0 to 1/𝐡. (This is a test-channel model often used for the test of equalizers that reflects particularly hostile environment with very small coherence bandwidth and very deep fading.) We have observed 1000 channel realizations and assumed the target BEP π‘ƒπ‘’π‘˜=10βˆ’3 for all links (for all π‘˜).

For the comparison purposes, we present results of our proposed framework and the reduced-complexity algorithm of finding the optimal tax rate together with the results of the greedy algorithm (that assigns the frequencies to the players with the highest CNR values at these frequencies) and Round-Robin algorithm of resource distribution. Although both of these algorithms can be only implemented in a centralized manner, they give the two opposite extremes: either maximum spectral efficiency or maximum fairness for the case of the whole used bandwidth.

Let us first analyze the network performance in the case of the users’ SelBM and the influence of the tax rate π‘Ÿ0, the restricted amount of bandwidth 𝐡I, and the number of players 𝐾 on the network behavior. In Figures 2 and 3 we observe the averaged sum throughput πœ‚ST defined by (20) and the network spectral efficiency πœ‚SE defined by (21), respectively, for the two-path channel model. As we can see, there is some optimum tax rate π‘Ÿβˆ—ST0 that maximizes πœ‚ST when 𝐡I is not too small for a given 𝐾. Otherwise, πœ‚ST is constant for low tax rates, and then, for higher tax rates decreases to zero. The optimal tax rate π‘Ÿβˆ—SE0 that maximizes πœ‚SE is close to the barrage tax.

To better understand the mechanism of increasing tax rates in the network, already discussed in the previous section, let us analyze Figures 4 and 5. There, the percentage of served nodes and the percentage of used bandwidth are shown versus the tax rate π‘Ÿ0. (We assume that the node is served if it is able to acquire a portion of bandwidth satisfying her target BEP.) As we can see, low taxes allow to utilize most of the bandwidth and serve most of the users, again when 𝐡I is properly chosen. In general, restricting the players in the amount of bandwidth they can take at their turn has negative influence on πœ‚ST (the maximum is always achieved for 𝐡𝐼=𝐡) and on the percentage of used bandwidth, and positive influence on πœ‚SE and the percentage of served nodes, but only for relatively low tax rates. For higher tax rates, for which πœ‚ST and the percentage of used bandwidth dramatically drop, both πœ‚SE and the percentage of served nodes are not dependent on 𝐡𝐼. Thus, our first conclusion is that for the fairness of the resource distribution, it is better to apply 𝐡𝐼=𝐡, and calibrate just the tax rate to optimize πœ‚ST rather than πœ‚SE. Similar (analogous) conclusions can be derived for the SocBM and for the other channel model.

In Figures 6 and 7 we can observe the tax rates π‘Ÿβˆ—ST0 optimizing πœ‚ST for both the SelBM and the SocBM and in the case of two-path and six-path channel models. These tax rates have been found using the algorithm defined in the previous section. As we can see, for the SelBM the optimum tax rates π‘Ÿβˆ—ST0 for different 𝐡𝐼, converge to the same value as 𝐾 increases. It is not the case for the SocBM. Moreover, for the two-path channel model, the values of the optimum tax rate π‘Ÿβˆ—ST0 are higher than for the six-path channel model.

The tax-rate that optimizes πœ‚SE does not depend on 𝐾. For the SelBM, it also does not depend on 𝐡𝐼, but only varies for different channel models. For the SocBM, π‘Ÿβˆ—SE0 depends on both, 𝐡𝐼 and the channel model. This can be observed in Table 1.

In Figures 8 and 9 one can observe the average sum throughput πœ‚ST resulting from the optimal taxation versus the number of competing players 𝐾, for both the SelBM and the SocBM and in the case of two-path (Figure 8) and six-path (Figure 9) channel models. In Figures 10 and 11 the network spectral efficiency πœ‚SE is presented versus 𝐾 for the same cases. Note, that for a given channel model the achievable average sum throughput πœ‚ST is exactly the same for both behavior models: either SelBM or SocBM (although the respective optimal tax rates are different). The same holds for the achievable spectral efficiency πœ‚SE. The difference between the plots occurs for different channel models. The achievable πœ‚ST as well as πœ‚SE are higher for the six-path channel model than for the two-path channel model due to the fact that this six-path model presents higher diversity in the subbands qualities for each player, so the players can make better choices. Finally, we have observed that when the optimum tax rate π‘Ÿβˆ—ST0 is applied in either scenario, 99–100% of nodes are served.

Note, that our framework cannot result in the maximal achievable sum throughput, which can be only obtained when the problem described by (1) is solved, which assumes complete-information of all links CSI and 𝐾-dimensional game that can be solved in either cooperative or noncooperative manner.

6. Conclusions

Here above, we have presented a game-theory-related framework for distributed allocation of spectrum resources in the opportunistic radio access networks. Contrary to the methods presented in the literature so far, in our game models, we do not assume the complete knowledge of the players CSI. Each player has the information on her own CSI only. Additionally, the taxation-rate parameter available in a data base and mandatory in the considered area is made known to the players through the broadcast channel (BCCH or the CPC). This significantly reduces the amount of control traffic in the network when compared with the frequent exchange of the all links CSI in the dedicated channels. Above, we have proposed a reduced-complexity algorithm of finding and tracing the optimum tax-rate value maximizing the network objective function. Our presented framework and the algorithm of finding the optimal taxation-rate result in high network benefit reflected in the sum throughput, but also in fairness of resource distribution (understood as the number of served nodes). The simulation results show that it is more beneficial for the network and for the individual players to use taxation with the tax rate maximizing the network sum throughput rather than to additionally limit the users in the maximum bandwidth they can acquire at the time. It is also more beneficial than maximization of the network spectral efficiency due to better utilization of the spectrum resources and higher percentage of served nodes. Simulation results also show that in the considered scenarios, when the optimal tax rate is applied the achievable sum throughput per frequency unit is as high as 5.5–6 [bits/s/Hz] (depending on the considered propagation model) for sufficiently high number of players. Moreover, in such a case, 99-100% of nodes are served in the network, that is, are able to acquire some resources satisfying their target BEP.

Appendices

A.

Below, we calculate the derivative of Μ‡πœπ‘˜(π‘π‘˜): πœ•πœ•π‘π‘˜Μ‡πœπ‘˜ξ€·π‘π‘˜ξ€Έ=πœ•πœ•π‘π‘˜ξ‚»ξ€œπ‘π‘˜0log2ξ€Ίπ›Ώπ‘˜Μ‡π›Ύπ‘˜ξ€·π‘“π‘˜ξ€Έ/𝛾0π‘˜ξ€·π‘π‘˜ξ€Έξ€»π‘‘π‘“π‘˜ξ‚Όβˆ’πœ•πœ•π‘π‘˜ξ€½π‘Ÿ0π‘π‘˜ξ€Ύ=πœ•πœ•π‘π‘˜ξ‚»ξ€œπ‘π‘˜0log2ξ€Ίπ›Ώπ‘˜Μ‡π›Ύπ‘˜ξ€·π‘“π‘˜ξ€Έξ€»π‘‘π‘“π‘˜ξ‚Όβˆ’πœ•πœ•π‘π‘˜ξ‚»ξ€œπ‘π‘˜0log2𝛾0π‘˜ξ€·π‘π‘˜ξ€Έξ€»π‘‘π‘“π‘˜ξ‚Όβˆ’π‘Ÿ0.(A.1)

Because the one integrand in the first summand does not depend on π‘π‘˜, and the one integrand in the second summand does not depend on π‘“π‘˜, further derivation of the above formula is the following: πœ•πœ•π‘π‘˜Μ‡πœπ‘˜ξ€·π‘π‘˜ξ€Έ=log2ξ€Ίπ›Ώπ‘˜Μ‡π›Ύπ‘˜ξ€·π‘π‘˜βˆ’πœ•ξ€Έξ€»πœ•π‘π‘˜ξ€½π‘π‘˜log2𝛾0π‘˜ξ€·π‘π‘˜ξ€Έξ€»ξ€Ύβˆ’π‘Ÿ0=log2ξ€Ίπ›Ώπ‘˜Μ‡π›Ύπ‘˜ξ€·π‘π‘˜ξ€Έξ€»βˆ’π‘π‘˜πœ•πœ•π‘π‘˜ξ€½log2𝛾0π‘˜ξ€·π‘π‘˜ξ€Έξ€»ξ€Ύβˆ’log2𝛾0π‘˜ξ€·π‘π‘˜ξ€Έξ€»βˆ’π‘Ÿ0=log2ξ€Ίπ›Ώπ‘˜Μ‡π›Ύπ‘˜ξ€·π‘π‘˜βˆ’1𝑏ln2π‘˜π›Ύ0π‘˜ξ€·π‘π‘˜ξ€Έπœ•πœ•π‘π‘˜ξ€½π›Ύ0π‘˜ξ€·π‘π‘˜ξ€Έξ€Ύβˆ’log2𝛾0π‘˜ξ€·π‘π‘˜ξ€Έξ€»βˆ’π‘Ÿ0.(A.2) Now, this formula does not have a closed form, that is, (πœ•/πœ•π‘π‘˜)Μ‡πœπ‘˜(π‘π‘˜) depends on the particular shapes of Μ‡π›Ύπ‘˜(π‘π‘˜) and of 𝛾0π‘˜(π‘π‘˜) (where the later function in turn depends on the shape of the former). As an example, let us consider the rural area with the two-path propagation. In such a case, the channel power characteristic can be described as π›Ύπ‘˜||𝛼(𝑓)=π‘˜||2||1βˆ’π›½π‘˜ξ€½ξ€·expβˆ’π‘—2πœ‹π‘“πœπ‘˜βˆ’πœ‘π‘˜||ξ€Έξ€Ύ2,(A.3) where π›Όπ‘˜ is the complex amplitude attenuation of the first path, π›½π‘˜ is the attenuation of the second path relative to the first path, πœπ‘˜ is the multipath delay spread, and πœ‘π‘˜ is average phase difference between the arriving multipath signal components. Let us denote πœ™π‘˜=2πœ‹π‘“πœπ‘˜βˆ’πœ‘π‘˜. For the sufficient distance between the transmitting and receiving antennas π›½π‘˜β‰ˆ1. Consequently, π›Ύπ‘˜(𝑓)≃𝛼2π‘˜ξ€·1βˆ’cosπœ™π‘˜+𝑗sinπœ™π‘˜ξ€Έξ€·1βˆ’cosπœ™π‘˜βˆ’π‘—sinπœ™π‘˜ξ€Έ=2𝛼2π‘˜ξ€·1βˆ’cosπœ™π‘˜ξ€Έ=2𝛼2π‘˜ξ€Ίξ€·1βˆ’cos2πœ‹π‘“πœπ‘˜βˆ’πœ‘π‘˜.ξ€Έξ€»(A.4) The above function is periodic and monotonically increasing for πœ™βˆˆ[0,πœ‹]. Therefore, its sorted (in a descending order) version Μ‡π›Ύπ‘˜(π‘“π‘˜) can be approximated as Μ‡π›Ύπ‘˜ξ€·π‘“π‘˜ξ€Έβ‰ƒπ΄π‘˜ξ€Ίξ€·1+cos2πœ‹π‘“π‘˜πœπ‘˜βˆ’πœ‘π‘˜ξ€Έξ€»,(A.5) where π΄π‘˜ is the proportionality constant.

Now, for our propagation model, let us find the formula for the water-level 𝛾0π‘˜(π‘π‘˜) dependent on the shape of function Μ‡π›Ύπ‘˜(π‘“π‘˜) and on the considered bandwidth π‘π‘˜. To this end we will integrate both sides of (9) resulting in π‘π‘˜π›Ύ0π‘˜ξ€·π‘π‘˜ξ€Έ=Ξ“π‘˜+ξ€œπ‘π‘˜01π›Ώπ‘˜Μ‡π›Ύπ‘˜ξ€·π‘“π‘˜ξ€Έπ‘‘π‘“π‘˜,(A.6) where Ξ“π‘˜ is the power limit for player π‘˜. Consequently 𝛾0π‘˜ξ€·π‘π‘˜ξ€Έβ‰ƒπ‘π‘˜ξƒ―Ξ“π‘˜+ξ€œπ‘π‘˜01π΄π‘˜π›Ώπ‘˜ξ€Ίξ€·1+cos2πœ‹π‘“π‘˜πœπ‘˜βˆ’πœ‘π‘˜ξ€Έξ€»π‘‘π‘“π‘˜ξƒ°βˆ’1,(A.7) and because ∫(1/cosπœ™π‘˜)π‘‘πœ™π‘˜=tgπœ™π‘˜/2, we obtain 𝛾0π‘˜ξ€·π‘π‘˜ξ€Έβ‰ƒπ‘π‘˜ξ‚»Ξ“π‘˜+12πœ‹πœπ‘˜π΄π‘˜π›Ώπ‘˜ξ‚ƒξ‚€tgπœ‹π‘π‘˜πœπ‘˜βˆ’πœ‘π‘˜2ξ‚ξ‚€πœ‘βˆ’tgπ‘˜2ξ‚Όξ‚ξ‚„βˆ’1.(A.8) Using the above expression we will derive the second term in (A.2): π‘π‘˜π›Ύ0π‘˜ξ€·π‘π‘˜ξ€Έπœ•πœ•π‘π‘˜ξ€½π›Ύ0π‘˜ξ€·π‘π‘˜β‰ƒξ‚»Ξ“ξ€Έξ€Ύπ‘˜+12πœ‹πœπ‘˜π΄π‘˜π›Ώπ‘˜ξ‚ƒξ‚€tgπœ‹π‘π‘˜πœπ‘˜βˆ’πœ‘π‘˜2ξ‚ξ‚€πœ‘βˆ’tgπ‘˜2ξ‚Όβ‹…ξ€·ξ‚ξ‚„βˆ’π‘π‘˜ξ€Έξ‚»Ξ“π‘˜+12πœ‹πœπ‘˜π΄π‘˜π›Ώπ‘˜Γ—ξ‚ƒξ‚€tgπœ‹π‘π‘˜πœπ‘˜βˆ’πœ‘π‘˜2ξ‚ξ‚€πœ‘βˆ’tgπ‘˜2ξ‚Όξ‚ξ‚„βˆ’2⋅2π΄π‘˜π›Ώπ‘˜cos2ξ‚€πœ‹π‘π‘˜πœπ‘˜βˆ’πœ‘π‘˜2ξ‚ξ‚‡βˆ’1=ξ€·βˆ’π‘π‘˜ξ€ΈΓ—ξ‚»Ξ“π‘˜+12πœ‹πœπ‘˜π΄π‘˜π›Ώπ‘˜ξ‚ƒξ‚€tgπœ‹π‘π‘˜πœπ‘˜βˆ’πœ‘π‘˜2ξ‚ξ‚€πœ‘βˆ’tgπ‘˜2ξ‚Όξ‚ξ‚„βˆ’1⋅2π΄π‘˜π›Ώπ‘˜cos2ξ‚€πœ‹π‘π‘˜πœπ‘˜βˆ’πœ‘π‘˜2ξ‚ξ‚‡βˆ’1.(A.9) We can now substitute (A.5)–(A.9) to (A.2), which results in πœ•πœ•π‘π‘˜Μ‡πœπ‘˜ξ€·π‘π‘˜ξ€Έ=β„±π‘˜ξ€·π‘π‘˜ξ€Έβˆ’π‘Ÿ0,(A.10) where the first term on the right-hand side of the above equation β„±π‘˜(π‘π‘˜) does not depend on π‘Ÿ0, and for the considered channel model can be approximated as β„±π‘˜ξ€·π‘π‘˜ξ€Έβ‰ƒlog2ξ€½π΄π‘˜π›Ώπ‘˜ξ€Ίξ€·1+cos2πœ‹π‘π‘˜πœπ‘˜βˆ’πœ‘π‘˜ξ€Έξ€»ξ€Ύ+log2ξ‚»π‘π‘˜βˆ’1ξ‚ΈΞ“π‘˜+12πœ‹πœπ‘˜π΄π‘˜π›Ώπ‘˜ξ‚ƒξ‚€tgπœ‹π‘π‘˜πœπ‘˜βˆ’πœ‘π‘˜2ξ‚ξ‚€πœ‘βˆ’tgπ‘˜2+π‘ξ‚Ήξ‚Ήξ‚Όπ‘˜ξ‚»Ξ“ln2π‘˜+12πœ‹πœπ‘˜π΄π‘˜π›Ώπ‘˜ξ‚ƒξ‚€tgπœ‹π‘π‘˜πœπ‘˜βˆ’πœ‘π‘˜2ξ‚ξ‚€πœ‘βˆ’tgπ‘˜2ξ‚ξ‚Ήξ‚Όβˆ’1×2π΄π‘˜π›Ώπ‘˜cos2ξ‚€πœ‹π‘π‘˜πœπ‘˜βˆ’πœ‘π‘˜2ξ‚ξ‚‡βˆ’1.(A.11) The above formula can be further simplified, when we assume that the phase difference between the arriving two-path waveform components is negligible due to similar distance that both waves travel, that is, when πœ‘π‘˜β‰ˆ0. Note, that formula (A.10) is very general, and function β„±(π‘π‘˜) can be defined in a number of ways depending on the assumed propagation environment models.

B.

Below, we calculate the derivative of Μ‡πœ‰(π‘π‘˜,π‘π‘˜): πœ•πœ•π‘π‘˜Μ‡πœ‰π‘˜ξ€·π‘π‘˜,π‘π‘˜ξ€Έ=πœ•πœ•π‘π‘˜ξ€œξ‚»ξ‚Έπ‘π‘˜0log2ξ€·π›Ώπ‘˜Μ‡π›Ύπ‘˜ξ€·π‘“π‘˜ξ€Έ/𝛾0π‘˜ξ€·π‘π‘˜ξ€Έξ€Έπ‘‘π‘“π‘˜ξ‚ΉΓ—ξ€·π΅π‘˜βˆ’π‘π‘˜βˆ’π‘π‘˜ξ€Έξ‚Όβˆ’πœ•πœ•π‘π‘˜ξ€½π‘Ÿ0π‘π‘˜ξ€Ύ=πœ•πœ•π‘π‘˜ξ‚»ξ€œπ‘π‘˜0log2ξ€·π›Ώπ‘˜Μ‡π›Ύπ‘˜ξ€·π‘“π‘˜ξ€Έ/𝛾0π‘˜ξ€·π‘π‘˜ξ€Έξ€Έπ‘‘π‘“π‘˜ξ‚ΌΓ—ξ€·π΅π‘˜βˆ’π‘π‘˜βˆ’π‘π‘˜ξ€Έβˆ’ξ€œπ‘π‘˜0log2ξ€·π›Ώπ‘˜Μ‡π›Ύπ‘˜ξ€·π‘“π‘˜ξ€Έ/𝛾0π‘˜ξ€·π‘π‘˜ξ€Έξ€Έπ‘‘π‘“π‘˜βˆ’π‘Ÿ0=πœ•πœ•π‘π‘˜ξ‚»ξ€œπ‘π‘˜0log2ξ€Ίπ›Ώπ‘˜Μ‡π›Ύπ‘˜ξ€·π‘“π‘˜ξ€Έξ€»π‘‘π‘“π‘˜ξ‚Όξ€·π΅π‘˜βˆ’π‘π‘˜βˆ’π‘π‘˜ξ€Έβˆ’πœ•πœ•π‘π‘˜ξ‚»ξ€œπ‘π‘˜0log2𝛾0π‘˜ξ€·π‘π‘˜ξ€Έξ€»π‘‘π‘“π‘˜ξ‚Όβ‹…ξ€·π΅π‘˜βˆ’π‘π‘˜βˆ’π‘π‘˜ξ€Έβˆ’ξ€œπ‘π‘˜0log2ξ€·π›Ώπ‘˜Μ‡π›Ύπ‘˜ξ€·π‘“π‘˜ξ€Έ/𝛾0π‘˜ξ€·π‘π‘˜ξ€Έξ€Έπ‘‘π‘“π‘˜βˆ’π‘Ÿ0=ξ‚»log2ξ€Ίπ›Ώπ‘˜Μ‡π›Ύπ‘˜ξ€·π‘π‘˜βˆ’πœ•ξ€Έξ€»πœ•π‘π‘˜ξ€½π‘π‘˜log2𝛾0π‘˜ξ€·π‘π‘˜ξ‚ΌΓ—ξ€·π΅ξ€Έξ€»ξ€Ύπ‘˜βˆ’π‘π‘˜βˆ’π‘π‘˜ξ€Έβˆ’ξ€œπ‘π‘˜0log2ξ€Ίπ›Ώπ‘˜Μ‡π›Ύπ‘˜ξ€·π‘“π‘˜ξ€Έξ€»π‘‘π‘“π‘˜+ξ€œπ‘π‘˜0log2𝛾0π‘˜ξ€·π‘π‘˜ξ€Έξ€»π‘‘π‘“π‘˜βˆ’π‘Ÿ0=ξ€½log2ξ€Ίπ›Ώπ‘˜Μ‡π›Ύπ‘˜ξ€·π‘π‘˜ξ€Έξ€»βˆ’log2𝛾0π‘˜ξ€·π‘π‘˜Γ—ξ€·π΅ξ€Έξ€»ξ€Ύπ‘˜βˆ’π‘π‘˜βˆ’π‘π‘˜ξ€Έβˆ’ξƒ―1𝑏ln2π‘˜π›Ύ0π‘˜ξ€·π‘π‘˜ξ€Έπœ•πœ•π‘π‘˜ξ€½π›Ύ0π‘˜ξ€·π‘π‘˜ξƒ°Γ—ξ€·π΅ξ€Έξ€Ύπ‘˜βˆ’π‘π‘˜βˆ’π‘π‘˜ξ€Έβˆ’ξ€œπ‘π‘˜0log2ξ€Ίπ›Ώπ‘˜Μ‡π›Ύπ‘˜ξ€·π‘“π‘˜ξ€Έξ€»π‘‘π‘“π‘˜+π‘π‘˜log2𝛾0π‘˜ξ€·π‘π‘˜ξ€Έξ€»βˆ’π‘Ÿ0.(B.1) We can write the above expression in a simpler form: πœ•πœ•π‘π‘˜Μ‡πœ‰π‘˜ξ€·π‘π‘˜,π‘π‘˜ξ€Έξ€·π‘=π’’π‘˜,π‘π‘˜ξ€Έβˆ’π‘Ÿ0,(B.2) and we can substitute expressions (A.5)–(A.9) to formula (B.1) to obtain the derivative of function Μ‡πœ‰π‘˜(π‘π‘˜,π‘π‘˜) for the rural channel model and the expression for π’’π‘˜(π‘π‘˜): π’’π‘˜ξ€·π‘π‘˜ξ€Έβ‰ƒξ€½log2ξ€½π΄π‘˜π›Ώπ‘˜ξ€Ίξ€·1+cos2πœ‹π‘π‘˜πœπ‘˜βˆ’πœ‘π‘˜β‹…ξ€·π΅ξ€Έξ€»ξ€Ύξ€Ύπ‘˜βˆ’π‘π‘˜βˆ’π‘π‘˜ξ€Έ+log2ξ‚»π‘π‘˜βˆ’1ξ‚ΈΞ“π‘˜+12πœ‹πœπ‘˜π΄π‘˜π›Ώπ‘˜ξ‚ƒξ‚€tgπœ‹π‘π‘˜πœπ‘˜βˆ’πœ‘π‘˜2ξ‚ξ‚€πœ‘βˆ’tgπ‘˜2ξ‚Γ—ξ€·π΅ξ‚Ήξ‚Ήξ‚Όπ‘˜βˆ’π‘π‘˜βˆ’π‘π‘˜ξ€Έ+π‘π‘˜Γ—ξ‚»Ξ“ln2π‘˜+12πœ‹πœπ‘˜π΄π‘˜π›Ώπ‘˜ξ‚ƒξ‚€tgπœ‹π‘π‘˜πœπ‘˜βˆ’πœ‘π‘˜2ξ‚ξ‚€πœ‘βˆ’tgπ‘˜2ξ‚Όξ‚ξ‚„βˆ’1×2π΄π‘˜π›Ώπ‘˜cos2ξ‚€πœ‹π‘π‘˜πœπ‘˜βˆ’πœ‘π‘˜2ξ‚ξ‚‡βˆ’1ξ€·π΅π‘˜βˆ’π‘π‘˜βˆ’π‘π‘˜ξ€Έβˆ’ξ€œπ‘π‘˜0log2ξ€½π΄π‘˜π›Ώπ‘˜ξ€Ίξ€·1+cos2πœ‹π‘“π‘˜πœπ‘˜βˆ’πœ‘π‘˜ξ€Έξ€»ξ€Ύπ‘‘π‘“π‘˜βˆ’π‘π‘˜log2ξ‚»1π‘π‘˜ξ‚ΈΞ“π‘˜+12πœ‹πœπ‘˜π΄π‘˜π›Ώπ‘˜ξ‚ƒξ‚€tgπœ‹π‘π‘˜πœπ‘˜βˆ’πœ‘π‘˜2ξ‚ξ‚€πœ‘βˆ’tgπ‘˜2.ξ‚Ήξ‚Ήξ‚Ό(B.3)