Computation over WideBand MAC: Improved Achievable Rate through SubFunction Allocation
Abstract
Future networks are expected to connect an enormous number of nodes wirelessly using wideband transmission. This brings great challenges. To avoid collecting a large amount of data from the massive number of nodes, computation over multiaccess channel (CoMAC) is proposed to compute a desired function over the air utilizing the signalsuperposition property of MAC. Due to frequency selective fading, wideband CoMAC is more challenging and has never been studied before. In this work, we propose the use of orthogonal frequency division multiplexing (OFDM) in wideband CoMAC to transmit functions in a similar way to bit sequences through division, allocation and reconstruction of function. An achievable rate without any adaptive resource allocation is derived. To prevent a vanishing computation rate from the increase of the number of nodes, a novel subfunction allocation of subcarriers is derived. Furthermore, we formulate an optimization problem considering power allocation. A spongesqueezing algorithm adapted from the classical waterfilling algorithm is proposed to solve the optimal power allocation problem. The improved computation rate of the proposed framework and the corresponding allocation has been verified through both theoretical analysis and simulation.
I Introduction
Current 5G and internet of things envision the interconnections of up to 1 trillion products, machines and devices through wideband transmission[fettweis20145g, al2015internet, gupta2015survey]. With such an enormous number of nodes, it is impractical to transmit data using conventional orthogonal multiaccess schemes for the future networks because this would result in excessive network latency and low spectrum utilization efficiency.
To solve this problem, computation over multiaccess channel (CoMAC) was proposed to exploit the signalsuperposition property of MAC to compute the desired function through concurrent node transmissions that combines computation with communication[goldenbaum2015nomographic, goldenbaum2013robust, abari2016over, goldenbaum2014channel, kortke2014analog, zhu2018over, wu2019experimental, nazer2007computation, appuswamy2014computing, erez2005lattices, wu2018experimental, nazer2011compute, jeon2014computation, Jeon2016Opportunistic, wang2015interactive, goldenbaum2014computation]. A straightforward application of CoMAC to attain transmission and computation in parallel is wireless sensor networks whose goal of communication is typically for a fusion center to obtain a function value of the sensor readings (e.g., arithmetic mean, polynomial or number of active nodes), rather than store all readings from all sensors. Apart from this, CoMAC can also be applied in networks that focus on computing a class of so called nomographic functions of distributed data via concurrent node transmissions[goldenbaum2015nomographic].
A simple computation framework, called analog CoMAC, has been studied in [goldenbaum2013robust, abari2016over, goldenbaum2014channel, kortke2014analog, zhu2018over, wu2019experimental]. In analog CoMAC, all nodes participate in the transmission. The authors used preprocessing at each node and postprocessing at the fusion center to deal with fading and compute or estimate functions[goldenbaum2013robust]. The designs of preprocessing and postprocessing were used to compute linear and nonlinear functions discussed in[abari2016over], and the effect of channel estimation error was characterized in [goldenbaum2014channel]. In order to verify whether the analog CoMAC is feasible in practice, an implementation using software defined radio was built in [kortke2014analog]. A multifunction computation method utilizing a multiantenna fusion center to collect data transmitted by a cluster of multiantenna multimodal sensors has been presented in [zhu2018over]. Besides, using the multiplexing gain to compute multifunction has been discussed in [wu2019experimental]. In summary, the simple analog CoMAC has led to an active area focusing on the design and implementation techniques for receiving a desired function of concurrent signals.
The main limitation of analog CoMAC is the lack of robust design against noise. Hence, digital CoMAC was proposed to use joint sourcechannel coding in [nazer2007computation, appuswamy2014computing, erez2005lattices, nazer2011compute, jeon2014computation, Jeon2016Opportunistic, wang2015interactive, goldenbaum2014computation, wu2018experimental] to reduce the noise effect. The potential of linear source coding was discussed in [nazer2007computation], and its application was presented in [appuswamy2014computing] for the function computation over MAC. Compared with linear source coding, nested lattice coding can approach the performance of standard random coding [erez2005lattices]. This latticebased CoMAC was extended to a general framework in [nazer2011compute] for relay networks with linear channels and additive white Gaussian noise. In order to combat nonuniform fading, a uniformforcing transceiver design was given in [wu2018experimental]. Available computation rates were given in [nazer2011compute, jeon2014computation, Jeon2016Opportunistic] for digital CoMAC through theoretical analysis. In [Jeon2016Opportunistic], Jeon and Bang found that the computation rates achieved by the above computation techniques decrease as the number of nodes increases, and eventually go to zero due to fading MAC. To handle this problem, they proposed an opportunistic CoMAC where a subset of nodes opportunistically participate in the transmission at each time, which provided a nonvanishing computation rate even when the number of nodes in the network goes to infinity.
To the best of our knowledge, the above mentioned works only consider a narrowband CoMAC. A wideband CoMAC has not been discussed before. Wideband signal transmission has to face frequency selective fading because the coherence bandwidth of the channel is smaller than the signal bandwidth. In the time domain, the multipath interference is not easy to remove[tse2005fundamentals]. Orthogonal frequency division multiplexing (OFDM) is an effective solution to intersymbol interference caused by a dispersive channel[hasan2009energy, mohanram2007joint, chen2008fast], and has been applied in many communication fields.
Motivated by the above observations, we employ OFDM in the implementation of CoMAC. However, the conventional OFDM systems which load bit sequences of different lengths into different subcarriers cannot be used directly because CoMAC transmits a desired function over the air, not a bit sequence. Thus, the desired function as a whole cannot be allocated into several subcarriers. To deal with the issue, we propose a method of subfunction allocation for the division, allocation and reconstruction of the desired function. This method divides the desired function to some subfunctions, allocates these subfunctions into different subcarriers, and reconstructs the desired function at the fusion center. The theoretical expression of achievable computation rate is derived based on the classical results of nested lattice coding[jeon2014computation, Jeon2016Opportunistic, wang2015interactive, goldenbaum2014computation]. It is shown to provide a nonvanishing computation rate as the number of nodes increases. An optimization considering power allocation is further discussed, and a spongesqueezing algorithm adapted from the classical waterfilling algorithm is proposed to solve the optimal power allocation problem. Owing to subfunction allocation and power allocation, the computation rate is improved significantly. Our contributions can be summarized as follows:

Novel CoMACOFDM. In order to tackle frequency selective fading, we exploit the mechanism of OFDM to implement wideband CoMAC. Unlike the conventional CoMAC schemes, CoMACOFDM divides a desired function to subfunctions, allocates subfunctions to each subcarrier and reconstructs the desired function at the fusion center.

Improved computation rate. A novel subfunction allocation is designed to assign subfunctions into subcarriers, which provides a nonvanishing computation rate as the number of nodes increases. After that, the theoretical expression of improved computation rate is derived.

Optimal power allocation. An optimal power allocation is considered to improve computation rate. Since the solution of our optimization is not suitable to utilize the classical waterfilling algorithm, we propose a spongesqueezing algorithm adapted from the classical waterfilling algorithm to carry out the power allocation, which provides a simple and concise interpretation of the necessary optimality conditions.
The paper is organized as follows. Section II introduces the definitions of CoMAC and the system model for CoMACOFDM. In Section III, we summarize the main results of this paper and compare them with the existing results. Section IV presents the proposed CoMACOFDM with subfunction allocation in detail and analyzes the computation rate. Section V focuses on the improvement of the computation rate using power allocation. Simulation results are presented in Section VI and conclusions are given in Section VII.
Ii Preliminaries & Network Model
In this section, we give some definitions before introducing CoMAC. Throughout this paper, we denote , and . For a set , denotes the cardinality of . is the ceiling function as . Let the entropy of a random variable be and denote the diagonal matrix of which the diagonal elements are from to . A set is written as for short. represents the transpose of a vector or matrix.
Iia CoMAC
A classical framework of CoMAC is depicted in Fig. 1. In this framework, the fusion center is designed to compute a desired function from nodes. Each node obtains data from a corresponding random source for times, then giving a data vector of length .
Definition 1 (Data Matrix)
Assuming that is a random source from which the th node gets data, let be a random vector associated with a joint probability mass function . A data matrix represents the data from nodes during time slots. It is expressed as
(1) 
where , , is the th data of the th node from the random source , is the th data of all nodes and is the data vector of node . Note that belongs to , which means it is mapped to a number between 0 and through quantization. And is independently drawn from .
Using Definition 1, the desired function determined by the random source vector can be expressed as , and its definition is in the following.
Definition 2 (Desired Function)
For all , every function
(2) 
that has to be computed at the fusion center is called a desired function where is independently drawn from (See Definition 1). Every function can be seen as a realization of . Thus, it has functions when each node gets data from each random source for times.
Remark 1
As studied in [giridhar2005computing, wang2015interactive, goldenbaum2014computation, jeon2014computation], CoMAC can be designed to compute different types of desired function. There are two typical functions that we focus on. A function whose values are in the set is called the arithmetic sum function where is the weighting factor for node and belongs to . The arithmetic sum function is a weighted sum function, which includes the mean function for all nodes and the function for the active node only as special cases. Apart from this, a function with values in the set of is regarded as the type function where denotes the indicator function. As pointed out in [giridhar2005computing], any symmetric function such as mean, variance, maximum, minimum and median can be attained from the type function.
Definition 2 shows that the desired function is computed by all nodes participating in the computation. In some cases, we only need to choose nodes to participate in the computation each time. Then the desired function has to be divided into parts. The function which only uses part of the nodes is called subfunction, and its definition is given as follows.
Definition 3 (SubFunction)
Let
(3) 
denote a set where each element is the index from the chosen nodes. Suppose that and for all , A function is said to be a subfunction if and only if there exists a function satisfying .
Thus, the fusion center reconstructs the desired function based on subfunctions, when nodes are divided into parts.
In order to achieve reliable computations against noise, we apply a block code named sequences of nested lattice codes[nazer2011compute] throughout this paper. Based on the block coding, the definitions of encoding and decoding are given as follows.
Definition 4 (Encoding & Decoding)
Let denote the data vector for the th node whose length is (see Definition 1). Denote as the length transmitted vector for node . The received vector whose length is is given by at fusion center. Assuming there is a block code with length , the encoding and decoding functions can be expressed as follows.

Encoding Functions: the univariate function which generates is an encoding function of node . It means with length is mapped to a transmitted vector with length for node .

Decoding Functions: the decoding function is used to estimate the th desired function , which satisfies . It implies the fusion center obtains desired functions depending on the whole received vector with length .
Considering the block code with length , the definition of computation rate[Jeon2016Opportunistic, jeon2014computation, goldenbaum2015nomographic, goldenbaum2014computation] can be given as follows.
Definition 5 (Computation rate)
The computation rate specifies how many function values can be computed per channel use within a predefined accuracy. It can be written as where is the number of function values (see Definition 2), is the length of the block code and is the entropy of . Apart from this, is achievable only if there is a length block code in order that the probability as increases.
IiB Fading MAC with Wideband signal
Wideband signal transmission has to face frequency selective fading because the coherence bandwidth of the channel is smaller than the signal bandwidth. OFDM is wellsuited to provide robustness against fading. Thus, we describe a CoMACOFDM system with subcarriers during OFDM symbols while the length of the block code is . Then, the th received OFDM symbol at the fusion center in the frequency domain can be given by
(4) 
where , is the power allocation matrix of node where each diagonal element is the power allocated to the th subcarrier, is a diagonal matrix of which the diagonal element is the frequency response of each subcarrier for node , is the transmitted diagonal matrix of node at the th OFDM symbol, the diagonal element in represents an element whose index is in the transmitted vector (See Definition 4) for node and the diagonal element of can be regarded as identically and independently distributed (i.i.d.) complex Gaussian random noise following .
According to the above expression, the received signal in the th subcarrier at the th OFDM symbol can be expressed as
(5) 
where , , is the th diagonal element from , from is the frequency response of the th subcarrier of node at the th OFDM symbol, is i.i.d. complex Gaussian random noise following .
The above equations assume perfect synchronization and perfect removal of intercarrier interference.
Iii Problem Statement & Main Results
In this section, we first introduce the conventional CoMAC schemes, their features and limitations. Thus, we describe our main results and the improvements by comparing them with conventional CoMAC schemes. For easy presentation in following sections, we define the ordered indexes to describe the indexes of nodes which are sorted by the cosponsoring channel gain. For the th subcarrier, let be the th element in the set of ordered indexes of such that , e.g., . Further, can be omitted if the number of subcarriers is .
Iiia Previous Works
In [jeon2014computation], the authors provided the computation rate for the arithmetic sum and type functions (See Remark 1) in fading MAC, i.e., the channel coefficients are i.i.d. and vary independently over time.
Theorem 1 (Conventional CoMAC Rate)
Based on [jeon2014computation, Theorem 3], the average computation rate with average power constraint at each time slot is given as
(6) 
where is the gain from the average power constraint, each channel response of the th node and are i.i.d. random variables and the condition follows on the ordered indexes.
Unfortunately, for wireless networks with a massive number of nodes, the computation rate approaches 0 as increases, which is a vanishing computation rate. To resolve the problem, an opportunistic CoMAC system has been discussed in [Jeon2016Opportunistic]. In each time slot, only nodes which have the largest channel gains participate in the transmission such that the indexes of the chosen nodes is in the set .
Theorem 2 (Opportunistic CoMAC Rate)
As shown in [Jeon2016Opportunistic, Theorem 1], for any satisfying , the computation rate of the opportunistic CoMAC is given by
(7) 
where is the th element of the set of ordered indexes .
Theorem 2 improved the computation rate in fading environments and provided a nonvanishing computation rate as increases.
However, the wideband CoMAC has not been discussed and the improvement of computation rate is still limited in [Jeon2016Opportunistic, giridhar2005computing, goldenbaum2014computation, goldenbaum2015nomographic, jeon2014computation, wang2015interactive, appuswamy2014computing, zhan2013linear, ma2012interactive]. Therefore, Our paper focuses on designing a framework of CoMAC for wideband transmission, which can be robust against both frequency selective fading and the vanishing computation rate. Apart from this, under our subfunction allocation and power allocation, we achieve a higher computation rate and provide a improved nonvanishing computation rate compared with the conventional CoMAC schemes.
IiiB Main Results
We employ OFDM design into wideband CoMAC and divide a desired function into subfunctions. By allocating these subfunctions into different subcarriers, the computation rate based on CoMACOFDM is presented in Theorem 3.
Theorem 3 (General CoMACOFDM Rate)
For any satisfying , the computation rate of wideband CoMACOFDM over fading MAC is given by
(8)  
(9) 
where the condition follows while , is the number of OFDM symbols as shown in Eqs. (4) and (5), is the number of nodes for a desired function (See Definition 2), is the number of the chosen nodes for a subfunction (See Definition 3), is the number of the subcarriers for each OFDM symbol, is the frequency response of the th subcarrier of node , is the power allocated to the th node in the th subcarrier at the th OFDM symbol.
Proof:
We refer to Section IV for proof. \qed
Remark 2
A simple computation rate for wideband CoMACOFDM system without subfunction allocation and adaptive power allocation can be easily obtained from Theorem 3.
Corollary 1 (Direct CoMACOFDM Rate)
For a straightforward wideband CoMACOFDM system, each desired function is assigned to each subcarrier directly and the node has an identical average power constraint in the th subcarrier, i.e., for all and . The computation rate can be given as
(10) 
Proof:
We refer to Section VA for proof. \qed
By setting in Corollary 1, the computation rate is the same as Theorem 1. Comparing (6) and (10), the computation rate for wideband CoMACOFDM can be improved by increasing the number of subcarriers . However, even though can keep increasing, it still has upper limit due to electronic devices. Thus, cannot be infinity in practice. The computation rates in Eq. (6) and Eq. (10) with a fixed will become vanishing as the number of nodes increases, and converge to 0 in fading MAC.
In order to deal with this, we propose subfunction allocation to avoid the vanishing computation rate even if the number of nodes in the network increases. Based on Theorem 3, the computation rate of CoMACOFDM with subfunction allocation can be expressed as follows.
Corollary 2 (CoMACOFDM Rate with SubFunction Allocation)
Suppose that a desired function is divided into subfunctions, and each subfunction is allocated to a subcarrier where the channel gains of the nodes are the largest. The computation rate of wideband CoMACOFDM with subfunction allocation and average power constraint can be given as
(11) 
Proof:
We refer to Section VA for proof. \qed
By setting in Corollary 2, the computation rate is the same as that in Theorem 2. Compared with Eq. (7), it shows that the proposed scheme can not only provide a nonvanishing computation rate but also improve the computation rate when the number of subcarriers increases.
In addition to average power constraint, OFDM can achieve a higher rate with adaptive power allocation for each node. The instantaneous computation rate for a OFDM symbol is obtained as follows.
Corollary 3 (CoMACOFDM Rate with SubFunction Allocation and Optimal Power Allocation)
The optimal computation rate with subfunction allocation at the th OFDM symbol is given as
(12) 
where is the optimal level for the th subcarrier at the th OFDM symbol, which will be further explained in Section VB .
Proof:
We refer to Section VB1 for the proof. \qed
Depending on the proposed subfunction allocation, we consider an optimization problem to achieve a higher computation rate in OFDM with total power constraint for each node. Based on the solution to the optimization problem, we propose a spongesqueezing algorithm perform the optimal power allocation method.
Iv SubFunction Allocation for CoMACOFDM
In this section, we explain the framework of CoMACOFDM, which is designed to transmit wideband signals and improve the computation rate. In Section IVA, we present subfunction allocation and expound some technical issues that need to be addressed. Then, we resolve these issues and provide the computation rate in Theorem 3 from Section IVB to Section IVC.
Iva SubFunction Allocation
As shown in Fig. 2, subfunction allocation includes three main parts, function division, function allocation and function reconstruction. We provide a simplified description on subfunction allocation based on CoMACOFDM system as follows.

Function Division. In Fig. 2, Node gets data from the corresponding random source and obtains the data vector . Then, desired functions is divided into subfunctions, which means only chosen nodes participate in the computation for a subfunction.

Function Allocation. Each subfunction is only allocated to the subcarrier where the channel gains of the chosen nodes are the largest. For a subfunction in Fig. 2, the chosen nodes are Node 1 and Node 4. When the channel gains of Nodes 1 and 4 are the largest ones in the th subcarrier at the th OFDM symbol, e.g., , the subfunction should be assigned to this subcarrier. After function allocation, all the nodes transmit the corresponding OFDM symbols over the air.

Function Reconstruction. At the fusion center in Fig. 2, two OFDM symbols are received, and each subfunction can be obtained in each subcarrier. Using the received subfunctions, desired functions are reconstructed by using the relationship between the subfunctions and the desired functions.
However, the above description of subfunction allocation is ideal, and the following practical issues should be addressed.

The order of the channel gains at each subcarrier is random and depends on the channel realizations in practice. This means that it is important to establish the rules for subfunction allocation.

The above description ignores the encoding and decoding of nested lattice code. It needs employ sequences of nested lattice codes into subfunction allocation.
Further discussion on solving these issues is in Section IVC.
IvB Computation Rate for SubFunctions
In the subfunction allocation shown in Section IVA, the fusion center computes one subfunction at a subcarrier. Each subfunction is only allocated into the subcarrier where the channel gains of the chosen nodes are the largest ones. The indexes of those chosen nodes are in a set for the th subcarrier at the th OFDM symbol. For the sake of fairness, we assume the bandwidth of CoMACOFDM with N subcarriers is the same as the bandwidth of the conventional CoMAC system which uses single carrier for transmission. It means the bandwidth of the subcarrier is of the bandwidth of the single carrier. Then, the computation rate of a subfunction at the th OFDM symbol can be given as follows.
Lemma 1
When the length of block , the instantaneous computation rate of single subfunction of the th subcarrier at the th OFDM symbol with additive white Gaussian noise and a variance of is given as
(13) 
where is the frequency response of the th subcarrier for node (See Eq. (5)), is the allocated power of the th node in the th subcarrier at the th OFDM symbol.
Proof:
[nazer2011compute, Theorem 4, Eq. (85)] shows that the computation rate with a noise variance of is given as
(14) 
where is the equation coefficient vector, is channel coefficient vector and is the power of each transmitted symbol.
By setting and using [nazer2011compute, Theorems 4 and 11], the computation rate can be written as
(15) 
Then by combining Eq. (15) with [jeon2014computation, Theorem 3], the instantaneous computation rate in fading MAC with additive white Gaussian noise whose variance is can be expressed as
(16) 
at the th time slot. Furthermore, the symbol of a subcarrier in OFDM lasts time slots while the symbol in the single carrier lasts 1 time slot. According Definition 5, the number of the function values in a subcarrier is of the number of the function values in single carrier, which means the computation rate of a subcarrier is of the computation rate in Eq. (16). In conclusion, Lemma 1 has been proved. \qed
IvC Reconstruction of SubFunctions
Only nodes are chosen from all nodes for a subfunction, which means there need be subfunctions to reconstruct a desired function. denotes a set where the elements are the indexes of the chosen nodes for a subfunction (See Definition 3). Then, the set
(17) 
includes all the possible subfunctions^{1}^{1}1For easy presentation, we use the element stands for the subfunction which is computed by these nodes in ., and the cardinality of is . According to Definition 3, the subfunctions can be reconstructed into the desired function if and only if and for all and . We define
(18) 
which contains all possible combinations that can be used to reconstruct the desired function, and .
Under our design, the length block code is transmitted through OFDM symbols. In order to simplify the derivation, the number of OFDM symbols instead of in Eq. (4), which means the number of all the subcarriers during OFDM symbols is . Based on subfunction allocation, all the subcarriers should be assigned to every subfunction. The set for all includes those subcarriers are assigned to subfunction . And let a set include the subcarriers that are from and allocated in the combination that contains the corresponding subfunction .
As mentioned in Issue i, Section IVA, the order of the indexes for nodes in each subcarrier is random depending on channel realizations in practice. It means the number of the subcarriers in is random. Thus, we use the following lemma to characterize the the minimum deterministic values of and .
Lemma 2
When is large, for all and , each set contains
(19) 
subcarriers.
Proof:
The probability, , for all is since channel gains are i.i.d. in the th subcarrier at the th OFDM symbol. Thus, holds depending on [csiszar2011information, Lemma 2.12]. It shows for all as n increases. Besides, in the set , there are combinations which include as an element for all . Therefore, each combination including subfunction is assigned subcarriers. \qed
According to the definition of , we find that there are subcarriers which are used to transmit subfunction in the combination over OFDM symbols. It also implies that the length of the transmitted signal vector of each node for subfunction is . To simplify the derivation, the transmitted signal vector is sent centrally in OFDM symbols in wideband CoMACOFDM. From Eq. (5), the received vector for subfunction based on the combination at the th OFDM symbol in the frequency domain can be given as
(20) 
Then from Eq. (5), the received signal in th subcarrier can be given as
(21) 
where is . From the conclusion of Lemma 1 and Eqs. (20)(21) , the average rate for computing a subfunction during OFDM symbols is given as
(22) 
where the condition follows because .
Note that we employ nested lattice code in the proposed system and the length of the transmitted signal vector of each node in is . In order to have this length for the transmitted signal vector, the corresponding length of data vector should be depending on Definition 5, which addresses Issue ii, Section IVA. Thus, a desired function reconstructed by combination can be described as (See Definition 3), and the number of desired function values reconstructed by combination is . Then, the number of desired function values for all is
(23) 
during OFDM symbols. As a result, the computation rate based on Definition 5 for computing the desired function via CoMACOFDM
is achievable as increases where the condition follows because of Lemma 2 and Eq. (23), the condition follows from Eq. (22), the cardinality of is from Eq. (18), the condition follows while and the condition follows due to the increase of . Hence, Theorem 3 is proved.
V Power Allocation for SubFunction Allocation
In this section, we propose two power allocation methods to implement CoMACOFDM with subfunction allocation. For average power allocation, the achievable rate is give in Section VA. We further consider an optimization problem and propose a spongesqueezing algorithm which provides a simple and concise interpretation of the necessary optimality conditions based on the solution in Section VB.
Va Average Power Allocation
Based on subfunction allocation discussed in Section IVA, we proposed a power allocation method with average power constraint.
Because only nodes participate in computing a subfunction in each subcarrier depending on Theorem 3, let represent the transmitted power in th subcarrier at the th OFDM symbol for the node , which can be given as
(24) 
In order to compare it with the conventional CoMAC which uses the single carrier to transmit signal in a fair way, we assume the bandwidth of OFDM symbols with subcarriers is the same as the bandwidth of the conventional single carrier symbols. Then, a complete OFDM symbol should use time slots to be transmitted while the number of transmitted symbols should be in the single carrier system. Hence, the average power in the th subcarrier at the th OFDM symbol is . From
(25) 
with , we can calculate and set
(26) 
We further put from Eq. (26) into Eq. (8) in Theorem 3, the computation rate for CoMACOFDM with subcarrier loading scheme can be written as
(27) 
where the condition follows because the channel gains in different subcarriers are i.i.d. random variables. In conclusion, Eq. (VA) is achievable as increases (See Theorem 3), which completes the proof of Corollary 2. Particularly, Corollary 1 can be proved by setting , which means each subcarrier is allocated a desired function.
VB Optimal Power allocation
Based on the above analysis, we demonstrate the computation rate based on CoMACOFDM in Section IVC. We further find that the computation rate using on the average power constraint for each node in each subcarrier is not optimal. Therefore, we consider the optimal power allocation for CoMACOFDM, discuss an optimization problem and propose an algorithm called spongesqueezing based on the solution.
VB1 Problem Formulation
Corollary 3 shows that the computation rate in fading MAC is the mean of the instantaneous computation rates, which means Eq. (8) can be rewritten as
(28) 
where the instantaneous computation rate at the th OFDM symbol is
(29) 
It implies that our power allocation should focus on the current OFDM symbol depending on the channel gains of each subcarrier. Thus, the optimization problem based on subfunction loading is given as
(30) 
From Eq. (29), we find that the problem is difficult to solve, because the function is nested in the function in the objective function and the ordered indexes also create difficulty. In order to make it tractable, we introduce a subfunction allocation matrix , and each element is given as
(31) 
Furthermore, in Eq. (29) implies that the computation rate is only determined by the node which has the minimum product for the chosen nodes in the th subcarrier. It means other nodes need not put much power into this subcarrier except for the node which has the minimum product, and the unused power can be allocated to the other subcarriers to improve the rate. Thus, we use the level to replace the function by making for all in the th subcarrier. Then, the optimization objective can be written as
(32) 
Therefore, the problem can be expressed as:
(33) 
where and are necessary, , , is the element of the subfunction allocation matrix at the th OFDM symbol.
Note that the objective function in Eq. (VB1) is concave since it is a sum of functions, which are concave themselves. Besides, the constraint set is convex as it is composed of linear constraints. Hence, the above optimization problem is a convex optimization problem, and has a unique maximizer. The Lagrangian function {dmath} L=MKN∑_g=1^N C^+( NM+Nη_g[m] )∑_i=1^K μ_i(∑_g=1^NG_i,g[m]η_g[m]ω_i,g[m]P ) with the complimentary slackness condition