Distributed Game Theoretic Optimization Of Frequency Selective Interference Channels: A Cross Layer Approach

Similar documents
The success of unlicensed broadband communication

COMPUTING the capacity region of the interference

Game theory and the frequency selective interference channel - A tutorial

Resource Allocation Challenges in Future Wireless Networks

Stable matching for channel access control in cognitive radio systems

Dynamic Fair Channel Allocation for Wideband Systems

Low-Complexity OFDMA Channel Allocation With Nash Bargaining Solution Fairness

Cooperative Spectrum Sharing in Cognitive Radio Networks: A Game-Theoretic Approach

Power Minimization for Multi-Cell OFDM Networks Using Distributed Non-cooperative Game Approach

Resource Allocation in Energy-constrained Cooperative Wireless Networks

Cognitive Radios Games: Overview and Perspectives

Multiple Antennas. Mats Bengtsson, Björn Ottersten. Basic Transmission Schemes 1 September 8, Presentation Outline

Pareto Optimization for Uplink NOMA Power Control

ENERGY EFFICIENT WATER-FILLING ALGORITHM FOR MIMO- OFDMA CELLULAR SYSTEM

3432 IEEE TRANSACTIONS ON INFORMATION THEORY, VOL. 53, NO. 10, OCTOBER 2007

Joint User Selection and Beamforming Schemes for Inter-Operator Spectrum Sharing

Channel Capacity Estimation in MIMO Systems Based on Water-Filling Algorithm

Adaptive Wireless. Communications. gl CAMBRIDGE UNIVERSITY PRESS. MIMO Channels and Networks SIDDHARTAN GOVJNDASAMY DANIEL W.

Degrees of Freedom of the MIMO X Channel

Optimal Power Allocation over Fading Channels with Stringent Delay Constraints

Modeling the Dynamics of Coalition Formation Games for Cooperative Spectrum Sharing in an Interference Channel

THE emergence of multiuser transmission techniques for

Keywords: Wireless Relay Networks, Transmission Rate, Relay Selection, Power Control.

Energy Efficient Power Control for the Two-tier Networks with Small Cells and Massive MIMO

Frequency and Power Allocation for Low Complexity Energy Efficient OFDMA Systems with Proportional Rate Constraints

Sequential Multi-Channel Access Game in Distributed Cognitive Radio Networks

Performance Analysis of Multiuser MIMO Systems with Scheduling and Antenna Selection

Analysis and Improvements of Linear Multi-user user MIMO Precoding Techniques

How (Information Theoretically) Optimal Are Distributed Decisions?

Jamming Games for Power Controlled Medium Access with Dynamic Traffic

On the Value of Coherent and Coordinated Multi-point Transmission

UNEQUAL POWER ALLOCATION FOR JPEG TRANSMISSION OVER MIMO SYSTEMS. Muhammad F. Sabir, Robert W. Heath Jr. and Alan C. Bovik

QoS Optimization For MIMO-OFDM Mobile Multimedia Communication Systems

Simple, Optimal, Fast, and Robust Wireless Random Medium Access Control

Joint Relaying and Network Coding in Wireless Networks

OPTIMAL POWER ALLOCATION FOR MULTIPLE ACCESS CHANNEL

University of Bristol - Explore Bristol Research. Peer reviewed version. Link to published version (if available): /VETECF.2011.

IEEE TRANSACTIONS ON SIGNAL PROCESSING, VOL. 59, NO. 1, JANUARY B. Related Works

The Performance Loss of Unilateral Interference Cancellation

International Journal of Advance Engineering and Research Development. Channel Estimation for MIMO based-polar Codes

SPECTRUM SHARING IN CRN USING ARP PROTOCOL- ANALYSIS OF HIGH DATA RATE

Auction-Based Optimal Power Allocation in Multiuser Cooperative Networks

Joint Transmitter-Receiver Adaptive Forward-Link DS-CDMA System

A Two-Layer Coalitional Game among Rational Cognitive Radio Users

Performance Analysis of Cognitive Radio based on Cooperative Spectrum Sensing

Fig.1channel model of multiuser ss OSTBC system

Delay Constrained Point to Multi-Point Scheduling in Wireless Fading Channels

A Game-Theoretic Analysis of Uplink Power Control for a Non-Orthogonal Multiple Access System with Two Interfering Cells

Inducing Cooperation for Optimal Coexistence in Cognitive Radio Networks: A Game Theoretic Approach

Imperfect Monitoring in Multi-agent Opportunistic Channel Access

Multi-Band Spectrum Allocation Algorithm Based on First-Price Sealed Auction

EasyChair Preprint. A User-Centric Cluster Resource Allocation Scheme for Ultra-Dense Network

Capacity and Optimal Resource Allocation for Fading Broadcast Channels Part I: Ergodic Capacity

SPECTRUM resources are scarce and fixed spectrum allocation

Maximising Average Energy Efficiency for Two-user AWGN Broadcast Channel

DEGRADED broadcast channels were first studied by

Dirty Paper Coding vs. TDMA for MIMO Broadcast Channels

IEEE JOURNAL OF SELECTED TOPICS IN SIGNAL PROCESSING, VOL. 6, NO. 2, APRIL

Decentralized and Fair Rate Control in a Multi-Sector CDMA System

On Multiple Users Scheduling Using Superposition Coding over Rayleigh Fading Channels

IN RECENT years, wireless multiple-input multiple-output

Chapter 12. Cross-Layer Optimization for Multi- Hop Cognitive Radio Networks

A Distributed Opportunistic Access Scheme for OFDMA Systems

ORTHOGONAL Frequency Division Multiplexing Access. Non-Cooperative Resource Competition Game by Virtual Referee in Multi-Cell OFDMA Networks

Communication over MIMO X Channel: Signalling and Performance Analysis

IN a large wireless mesh network of many multiple-input

Game Theory and Economics of Contracts Lecture 4 Basics in Game Theory (2)

Coalitional Games in Cooperative Radio Networks

DYNAMIC POWER ALLOCATION SCHEME USING LOAD MATRIX TO CONTROL INTERFERENCE IN 4G MOBILE COMMUNICATION SYSTEMS

ADAPTIVE RESOURCE ALLOCATION FOR WIRELESS MULTICAST MIMO-OFDM SYSTEMS

Research Collection. Multi-layer coded direct sequence CDMA. Conference Paper. ETH Library

Section Notes 6. Game Theory. Applied Math 121. Week of March 22, understand the difference between pure and mixed strategies.

1366 IEEE TRANSACTIONS ON COMMUNICATIONS, VOL. 53, NO. 8, AUGUST 2005

Beamforming with Imperfect CSI

Stability Analysis for Network Coded Multicast Cell with Opportunistic Relay

A Distributed Merge and Split Algorithm for Fair Cooperation in Wireless Networks

DOWNLINK BEAMFORMING AND ADMISSION CONTROL FOR SPECTRUM SHARING COGNITIVE RADIO MIMO SYSTEM

Nonstationary Resource Sharing with Imperfect Binary Feedback: An Optimal Design Framework for Cost Minimization

Degrees of Freedom of Multi-hop MIMO Broadcast Networks with Delayed CSIT

ISSN Vol.03,Issue.17 August-2014, Pages:

Resource Allocation on the MISO Interference Channel

Random Beamforming with Multi-beam Selection for MIMO Broadcast Channels

Distributed Approaches for Exploiting Multiuser Diversity in Wireless Networks

MULTICARRIER communication systems are promising

Opportunistic Beamforming Using Dumb Antennas

Multi-user Space Time Scheduling for Wireless Systems with Multiple Antenna

Competitive Distributed Spectrum Access in QoS-Constrained Cognitive Radio Networks

Diversity Techniques

Resource Allocation for OFDM and Multi-user. Li Wei, Chathuranga Weeraddana Centre for Wireless Communications

MULTIPATH fading could severely degrade the performance

Optimal Spectrum Management in Multiuser Interference Channels

ELEC E7210: Communication Theory. Lecture 11: MIMO Systems and Space-time Communications

Multiuser Scheduling and Power Sharing for CDMA Packet Data Systems

4740 IEEE TRANSACTIONS ON INFORMATION THEORY, VOL. 57, NO. 7, JULY 2011

Low-Complexity Beam Allocation for Switched-Beam Based Multiuser Massive MIMO Systems

LTE in Unlicensed Spectrum

Gradient-based scheduling and resource allocation in OFDMA systems

Measured propagation characteristics for very-large MIMO at 2.6 GHz

Optimal Utility-Based Resource Allocation for OFDM Networks with Multiple Types of Traffic

Distributed Power Control in Cellular and Wireless Networks - A Comparative Study

Transcription:

2010 IEEE 26-th Convention of Electrical and Electronics Engineers in Israel Distributed Game Theoretic Optimization Of Frequency Selective Interference Channels: A Cross Layer Approach Amir Leshem and Ephraim Zehavi School of Engineering Bar Ilan University Abstract In this paper we propose a smart multichannel extension of CSMA (SXMA with collision avoidance as a mean to perform distributed spectrum allocation using the Nash Bargaining solution. The smart carrier sensing replaces explicit CSI signalling required by standard approaches to cooperative game theoretic spectrum allocation in cognitive radio networks. It allows fast convergence of the network to a cooperative Pareto optimal point. Keywords: Cognitive radio, distributed spectrum management, game theory, Nash bargaining solution, interference channel, multiple access channel, carrier sensing. I. INTRODUCTION In recent years cooperative approaches derived from game theory have been used for efficient radio resource allocation. The most popular approach is the Nash Bargaining Solution (NBS [1],[2]. The NBS is based on four axioms that lead to a unique solution for the bargaining problem [3]. The NBS was applied for a flat fading channel (SISO [4], MISO [5] and MIMO [6], as well as for a frequency selective channel under two types of constraints: an average power constraint [1], [7] and a power mask constraint [2]. Most published results are based on a central controller that provides the allocation for NBS. The central controller collects gets from each player the rate vector and the feasible rate with competition. Base on this information the central controller allocates the subchannels to the player. Recently, Gao et al. [8] proposed a semi distributed algorithm with a coordinator using the Lagrange dual decomposition technique. In this approach each player solves the allocation problem independently based on the thresholds that are provided by the coordinator, and transmits his allocation to the coordinator. The coordinator updates the threshold and broadcasts them. It is highly desirable to achieve the game theoretic solution using a purely distributed approach. In recent years a new family of opportunistic medium access control protocols has emerged in the context of wireless sensor networks, where energy and bandwidth constraints are the dominant limitations [9], [10], [11]. These protocols use the channel response of each sensor in order to determine the back-off function. Recently, more advanced techniques have been proposed, where the back-off function depends on a network utility function [12]. However these techniques were single channel techniques. In this paper we establish a novel approach for distributed computation of the NBS with no central coordinator and with no explicit transmission of the achievable rates. The proposed approach is based on a novel multichannel smart carrier sensing and can be applied to other cooperative solutions like stable matching [13] and to generalized NBS [7]. The paper is organized as follows. In section II we define the Nash axioms and the Nash Bargaining solution. In section III we describe the cooperative bargaining over frequency selective channel under mask constraint. In section IV we presents a novel multichannel smart carrier sensing approach that solves the NBS without a central controller. Simulation results are shown in section V followed by conclusions. II. THE NASH BARGAINING SOLUTION In this section we describe the Nash bargaining solution (NBS and describe the axioms used in order to derive this cooperative solution. While, Nash equilibria [14] are inevitable whenever a non-cooperative zero sum game is played, they can lead to a substantial loss to all players, compared to a cooperative strategy in the non-zero sum case, where players can cooperate. An example of this situation is the well known prisoner s dilemma. The main issue in this case is how to achieve the cooperation in a stable manner and what rates can be achieved through cooperation. One of the earliest solutions to cooperative N player games is the Nash bargaining solution [15], [3]. The underlying structure for a Nash bargaining problem in an N player game is a set of outcomes of the bargaining process S R N which is compact and convex and a designated disagreement outcome d (which represents the agreement to disagree and solve the problem competitively. S can be considered as a set of outcomes of the possible joint strategies or states, Alternatively, some authors consider S as a set of states, d a disagreement state and a multiuser utility function U : S {d} R N. such that U (S {d} is compact and convex. The two approaches are identical and the first is obtained from the second by defining the game by the set of utilities of the possible outcomes. We will use the first formulation since it simplifies notation. However, in some cases we will define the outcomes of the game in terms of strategies. The set S in the first definition is then obtained by identifying it with U (S {d} of the second definition. 978-1-4244-8682-3/10/$26.00 c 2010 IEEE 000826

The Nash bargaining solution is a function F which assigns to each bargaining problem S {d} as above an element of S {d}, satisfying the following four axioms: Linearity. Assume that we consider the bargaining problem S {d } obtained from the problem S {d} by transformations: s i = α is i + β i, i = ( 1,..., N. d i = α i d i + β i. Then the solution satisfies F i S {d } = α i F i (S {d}+β i, for all i =1,..., N. Independence of irrelevant alternatives. This axiom states that if the bargaining solution of a large game T {d} is obtained in a small set S. Then the bargaining solution assigns the same solution to the smaller game, i.e., The irrelevant alternatives in T \S do not affect the outcome of the bargaining. Symmetry. If two players i<jare identical in the sense that S is symmetric with respect to changing the i th and the j th coordinates, then F i (S {d} =F j (S {d}. Equivalently, players which have identical bargaining preferences, get the same outcome at the end of the bargaining process. Pareto optimality. Ifs is the outcome of the bargaining then no other state t exists such that s < t (coordinate wise. A good discussion of these axioms can be found in [16]. Nash proved that there exists a unique solution to the bargaining problem satisfying these four axioms. The solution is obtained by solving the following problem: N s Nash = arg max (s n d n. (1 s S {d} Typically, one assumes that there exist at least one feasible s S such that d < s coordinatewise, but otherwise we can assume that the bargaining solution is d. We also define the Nash function F (s :S {d} R N F (s = (s n d n, (2 where s = [s 1,..., s N ] T The Nash bargaining solution is obtained by maximizing the Nash function over all possible states. Since the set of possible outcomes S {d} is compact and convex F (s has a unique maximum on the boundary of S {d}. Whenever the disagreement situation can be decided by a competitive game, it is reasonable to assume that the disagreement state is given by a Nash equilibrium of the relevant competitive game. In some cases there are other possibilities for the disagreement point. When the utility for user n is given by the rate R n, and d is the competitive Nash equilibrium, it is obtained by iterative waterfilling for general ISI channels. For the case of mask constraints the competitive solution is simply given by all users using the maximal PSD at all tones. III. BARGAINING OVER FREQUENCY SELECTIVE CHANNELS UNDER MASK CONSTRAINT In this section we define the cooperative game corresponding to the joint FDM/TDM achievable rate region for the frequency selective N user interference channel. We limit ourselves to the PSD mask constrained case since this case is actually the more practical one. In real applications, the regulator limits the PSD mask and not only the total power constraint. Let the K channel matrices at frequencies k =1,..., K be given by H k : k =1,..., K. Each player is allowed to transmit at maximum power p (k in the k th frequency bin. In non-cooperative scenario, under mask constraint, all players transmit at the maximal power they can use. Thus, all players choose the PSD, p i = p i (k :1 k K. The payoff for user i in the non-cooperative game is therefore given by: K Ri C (p i = log 2 (1 + SINR i (k. (3 k=1 Here, Ri C is the capacity available to player i given a PSD mask constraint distributions p, and SINR i (k is the received signal to interference-plus-noise ratio. Note that without loss of generality, and in order to simplify notation, we assume that the width of each bin is normalized to 1. We now define the cooperative game G TF (N,K,p. Definition 3.1: The FDM/TDM game G TF (N,K,p is a game between N players transmitting over K frequency bins under common PSD mask constraint. Each user has full knowledge of the channel matrices H k. The following conditions hold: 1 Player i transmits using a PSD limited by p i (k : k =1,..., K. 2 Strategies for player i are vectors α = [α i (1,..., α i (K] T where α i (k is the proportion of time player i uses the k th frequency channel. This is the TDM part of the strategy. 3 The utility of the i th player is given by R i = K k=1 R i(k = ( K k=1 α i(k log 2 1+ hii(k 2 p i(k (4 σi 2(k. Note that interference is avoided by time sharing at each frequency band, i.e only one player transmits at a given frequency bin at any time. Furthermore, since at each time instance each frequency is used by a single user, each user can transmit using maximal power. The Nash bargaining can be posed as an optimization problem max N ( i=1 Ri (α i Ri C subject to: k N i=1 α i(k =1, (5 i, k α i (k 0, i Ri C R i (α i, where, R i (α i = ( K k=1 α i(k log 2 1+ hii(k 2 P max(k σi 2(k = K k=1 α i(kr i (k. (6 This problem is convex and therefore can be solved efficiently using convex optimization techniques. Furthermore, the following theorem holds [2]: Theorem 3.1: Assume that for all i j the values {L ij (k :k =1,..., K} are all distinct. Then in the optimal 000827

solution at most ( N 2 frequencies are shared between different users. This theorem suggests, that when ( N 2 << K the optimal FDM NBS is very close to the joint FDM/TDM solution. It is obtained by allocating the common frequencies to one of the users. Therefore, we will turn now to a pure FDM game, where each frequency is allocated to a single user. IV. SMART CARRIER SENSING WITH COLLISION AVOIDANCE In this section we present a multichannel carrier sensing scheme with a time varying back-off function. The back-off function at each channel is derived based on the utility function of each user, the achievable rate at that channel and the total rate currently available for this user. Similarly to standard CSMA protocols, each user performs sensing on each of its intended channels up to a time defined by a back-off function. if the the back-off period expired and he did not sense any activity on the channel he is allowed to transmit its own packet. However, in contrast to standard CSMA, the smart carrier sensing has the following properties: 1 The purpose of the smart carrier sensing is not to provide a random access MAC, but to allow a distributed convergence of the different users to the Nash Bargaining solution for the resource allocation problem. As such convergence time should be short compared to the time scale of channel statistics variations. 2 The sensing is performed simultaneously on all the channels. This is reasonable in view of recent advances in cognitive radio technologies. 3 The back-off function depends on the achievable rate of the user and the total rate it has currently allocated on all channels. This allows to optimize complicated network utility functions. 4 The back-off function is both frequency and time dependent, since after a channel has been allocated, the back-off function, might allocate it to a different user at the next packet. 5 In contrast to our work on stable matching [13], users are allowed to transmit on multiple frequencies simultaneously. The specific choice of the back-off function at each time and each frequency is the important part in the smart carrier sensing. It should be derived from the total utility optimization. The smart carrier sensing algorithm is performed as described in Table I. Proper choice of the back-off function b n (k for each user and each carrier can lead to optimization of various network utility functions. In the next section we will discuss the solution of the Nash Bargaining problem using smart carrier sensing. TABLE I SMART CARRIER SENSING MULTIPLE ACCESS (SXMA Initialization: For each user n: Estimates the channels h n(k,k =1,..., K. l =0. Set α n(k =1. Set R tot (l (n = K k=1 αn(krn(k. Until Convergence For each user n (parallel: For k =1,..., K ( For each user n set b n(k =f R n(k,r tot (l (n Try to access to channel k using carrier sensing and back-off function b n(k. If user n won the channel or lost it, update R tot (l+1 (n and α n(k. l := l +1 end V. DISTRIBUTED COMPUTATION OF THE NASH BARGAINING SOLUTION We would like to optimize a total utility function of the form ( N K U (α 1,..., α N = log α n (kr n (k R nc. (7 k=1 Note that log(r is a monotonic function of R, and by Pareto optimality N α n (k =1 for all k. (8 and 0 α n (k 1 Computing R nc is easy to do without cooperation, since each user can estimate its own channel, and the interference it experiences during the channel acquisition and sensing phase. However, in random access environment we can use R nc =0 since a collision results in a packet loss. To simplify notation we assume R nc =0. We propose to use the carrier sensing function given by: b n (k = R(l+1 tot (n R n (k 2 (9 Theorem 5.1: Smart carrier sensing with the function b n (k converges to the the Nash Bargaining Solution. The proof of the theorem is postponed to the full version of this paper. However, we hint that the function b n (k was chosen, such that at each step the user with largest gradient of the Nash function at the given channel, receives the access to the channel. The ensures, that we optimize a first order approximation of the Nash Bargaining at each step. VI. SIMULATIONS We now present several examples for the performance of the smart carrier sensing for optimizing the Nash Bargaining function. In the first example we assume that two users share a channel composed of 32 i.i.d Rayleigh sub-channels. The 000828

transmit power to receiver noise ratio at each sub channel was 10 db and 13 db for users 1 and 2 respectively. Figure 1 describes the utility as a function of number, where each includes serial update of all frequency channels. Figure 2 demonstrates the rate of each user. We can clearly see the convergence within few s. We repeated the experiment with 16 users and 512 frequency channels. The SNR of the users was set from 0dB to 45dB in steps of 3dB. The results are depicted in Figures 3 and 4. We can clearly see the rapid convergence as well as the final rate of each user. Figure 5 demonstrates the rates achieved after convergence. Utility 80.9 80.85 80.8 80.75 80.7 80.65 80.6 80.55 80.5 80.45 Total utility (NBS vs. 80.4 0 20 40 60 80 100 Fig. 3. Total utility as a function of. 16 users, 512 frequency channels. SNR n =3(n 1dB. 3500 Total rate per user (NBS vs. 3000 2500 Rate 2000 1500 Fig. 1. Total utility as a function of. Two users, 32 frequency channels. SNR 1 =10dB, SNR 2 =13dB. 1000 500 0 0 20 40 60 80 100 Fig. 4. Rate for each user as a function of. 16 users, 512 frequency channels. SNR n =3(n 1dB. 700 Total rate per user (NBS vs. 600 500 Rate 400 300 Fig. 2. Rate for each user as a function of. Two users, 32 frequency channels. SNR 1 =10dB, SNR 2 =13dB. 200 100 VII. CONCLUSIONS In this paper we propose a new smart multichannel CSMA (SXMA protocol that can optimize the spectrum allocation under a total PSD mask constraint. The proposed protocol 0 0 200 400 600 800 1000 Fig. 5. Rate for each user as a function of. 16 users, 512 frequency channels. SNR n =3(n 1dB. 000829

achieves rapid convergence to the Nash Bargaining Solution (up to second order terms. The technique is easy to implement and fully distributed. REFERENCES [1] Z. Han, Z. Ji, and K. Liu, Fair multiuser channel allocation for OFDMA networks using the Nash bargaining solutions and coalitions, IEEE Trans. on Communications, vol. 53, pp. 1366 1376, Aug. 2005. [2] A. Leshem and E. Zehavi, Cooperative game theory and the Gaussian interference channel, IEEE Journal on Selected Areas in Communications, vol. 26, pp. 1078 1088, Sept. 2008. [3] J. Nash, Two-person cooperative games, Econometrica, vol. 21, pp. 128 140, Jan. 1953. [4] A. Leshem and E. Zehavi, Bargaining over the interference channel, in Proc. ISIT, pp. 2225 2229, July 2006. [5] E. Jorswieck and E. Larsson, The MISO interference channel from a game-theoretic perspective: A combination of selfishness and altruism achieves pareto optimality, Proc. ICASSP, pp. 5364 5367, 31 2008- April 4 2008. [6] M. Nokleby, A. Swindlehurst, Y. Rong, and Y. Hua, Cooperative power scheduling for wireless MIMO networks, Proc. GLOBECOM, pp. 2982 2986, Nov. 2007. [7] E. Zehavi and A. Leshem, Bargaining over the interference channel with total power constraints, in Proc. of GameNets, May 2009. [8] J. Gao, S. A. Vorobyov, and H. Jiang, Game theoretic solutions for precoding strategies over the interference channel, Proc. GLOBECOM, 2008. [9] Y. Chen and Q. Zhao, An integrated approach to energy aware medium access for wireless sensor networks, IEEE Trans. Signal Process., vol. 55, pp. 3429 3444, July 2007. [10] Y. Chen, Q. Zhao, V. Krishnamurthy, and D. Djonin, Transmission scheduling for optimizing sensor network lifetime: a stochastic shortest path approach, IEEE Trans. Signal Process., vol. 55, pp. 2294 2309, May 2007. [11] K. Cohen and A. Leshem, A time-varying opportunistic approach to lifetime maximization of wireless sensor networks, IEEE Transactions on Signal Processing, vol. 58, pp. 5307 5319, oct. 2010. [12] L. Jiang and J. Walrand, A distributed csma algorithm for throughput and utility maximization in wireless networks, in 46th Annual Allerton Conference on Communication, Control, and Computing, pp. 1511 1519, sep. 2008. [13] A. L. Y. Yaffe and E. Zehavi, Stable matching for channel access control for cognitive radio, in Proc. of Cognitive Information Processing CIP 2010, pp. 470 475, jun 2010. [14] J. Nash, The bargaining problem, Econometrica, vol. 18, pp. 155 162, Apr. 1950. [15] J. Nash, Non-cooperative games, The Annals of Mathematics, vol. 54, pp. 286 295, Sept. 1951. [16] G. Owen, Game theory. Academic Press, third ed., 1995. 000830