
In this paper, symbolbysymbol maximum likelihood (ML) detection is proposed
for a cooperative diffusionbased molecular communication (MC) system. In this
system, a fusion center (FC) chooses the transmitter's symbol that is more
likely, given the likelihood of the observations from multiple receivers (RXs).
We propose three different ML detection variants according to different
constraints on the information available to the FC, which enables us to
demonstrate tradeoffs in their performance versus the information available.
The system error probability for one variant is derived in closed form.
Numerical and simulation results show that the ML detection variants provide
lower bounds on the error performance of the simpler cooperative variants and
demonstrate that majority rule detection has performance comparable to ML
detection when the reporting is noisy.

In this paper, the error performance achieved by cooperative detection among
K distributed receivers in a diffusionbased molecular communication (MC)
system is analyzed and optimized. In this system, the receivers first make
local hard decisions on the transmitted symbol and then report these decisions
to a fusion center (FC). The FC combines the local hard decisions to make a
global decision using an NoutofK fusion rule. Two reporting scenarios,
namely, perfect reporting and noisy reporting, are considered. Closedform
expressions are derived for the expected global error probability of the system
for both reporting scenarios. New approximated expressions are also derived for
the expected error probability. Convex constraints are then found to make the
approximated expressions jointly convex with respect to the decision thresholds
at the receivers and the FC. Based on such constraints, suboptimal convex
optimization problems are formulated and solved to determine the optimal
decision thresholds which minimize the expected error probability of the
system. Numerical and simulation results reveal that the system error
performance is greatly improved by combining the detection information of
distributed receivers. They also reveal that the solutions to the formulated
suboptimal convex optimization problems achieve nearoptimal global error
performance.

In this paper, symbolbysymbol maximum likelihood (ML) detection is proposed
for a cooperative diffusionbased molecular communication (MC) system. In this
system, a fusion center (FC) chooses the transmitter (TX) symbol that is more
likely, given the likelihood of its observations from multiple receivers (RXs),
where the TX sends a common information symbol to all RXs. The transmission of
a sequence of binary symbols and the resultant intersymbol interference are
considered in the cooperative MC system. We propose five ML detection variants
according to different constraints on the knowledge at the FC. These five
variants demonstrate tradeoffs between their performance and the information
available. The system error probabilities for three ML detector variants are
derived, two of which are in closed form. Numerical and simulation results show
that the ML detection variants provide lower bounds on the error performance of
the simpler cooperative variants and demonstrate that these simpler cooperative
variants have error performance comparable to ML detectors when the reporting
from RXs to the FC is noisy.

This paper considers the communication for omniscience (CO) problem: A set of
users observe a discrete memoryless multiple source and want to recover the
entire multiple source via noisefree broadcast communications. We study the
problem of how to determine an optimal rate vector that attains omniscience
with the minimum sumrate, the total number of communications. The results
cover both asymptotic and nonasymptotic models where the transmission rates
are real and integral, respectively. We propose a modified decomposition
algorithm (MDA) and a sumrate increment algorithm (SIA) for the asymptotic and
nonasymptotic models, respectively, both of which determine the value of the
minimum sumrate and a corresponding optimal rate vector in polynomial time.
For the coordinate saturation capacity (CoordSatCap) algorithm, a nesting
algorithm in MDA and SIA, we propose to implement it by a fusion method and
show by experimental results that this fusion method contributes to a reduction
in computation complexity. Finally, we show that the separable convex
minimization problem over the optimal rate vector set in the asymptotic model
can be decomposed by the fundamental partition, the optimal partition of the
user set that determines the minimum sumrate, so that the problem can be
solved more efficiently.

For the accurate representation and reconstruction of bandlimited signals on
the sphere, an optimaldimensionality sampling scheme has been recently
proposed which requires the optimal number of samples equal to the number of
degrees of freedom of the signal in the spectral (harmonic) domain. The
computation of the spherical harmonic transform (SHT) associated with the
optimaldimensionality sampling requires the inversion of a series of linear
systems in an iterative manner. The stability of the inversion depends on the
placement of isolatitude rings of samples along colatitude. In this work, we
have developed a method to place these isolatitude rings of samples with the
objective of improving the wellconditioning of the linear systems involved in
the computation of the SHT. We also propose a multipass SHT algorithm to
iteratively improve the accuracy of the SHT of bandlimited signals.
Furthermore, we review the changes in the computational complexity and
improvement in accuracy of the SHT with the embedding of the proposed methods.
Through numerical experiments, we illustrate that the proposed variations and
improvements in the SHT algorithm corresponding to the optimaldimensionality
sampling scheme significantly enhance the accuracy of the SHT.

We present the generalized iterative residual fitting (IRF) for the
computation of the spherical harmonic transform (SHT) of bandlimited signals
on the sphere. The proposed method is based on the partitioning of the subspace
of bandlimited signals into orthogonal subspaces. There exist sampling schemes
on the sphere which support accurate computation of SHT. However, there are
applications where samples~(or measurements) are not taken over the predefined
grid due to nature of the signal and/or acquisition setup. To support such
applications, the proposed IRF method enables accurate computation of SHTs of
signals with randomly distributed sufficient number of samples. In order to
improve the accuracy of the computation of the SHT, we also present the
socalled multipass IRF which adds multiple iterative passes to the IRF. We
analyse the multipass IRF for different sampling schemes and for different
size partitions. Furthermore, we conduct numerical experiments to illustrate
that the multipass IRF allows sufficiently accurate computation of SHTs.

In this paper, we develop a new method for the fast and memoryefficient
computation of Slepian functions on the sphere. Slepian functions, which arise
as the solution of the Slepian concentration problem on the sphere, have
desirable properties for applications where measurements are only available
within a spatially limited region on the sphere and/or a function is required
to be analyzed over the spatially limited region. Slepian functions are
currently not easily computed for large bandlimits for an arbitrary spatial
region due to high computational and large memory storage requirements. For the
special case of a polar cap, the symmetry of the region enables the
decomposition of the Slepian concentration problem into smaller subproblems and
consequently the efficient computation of Slepian functions for large
bandlimits. By exploiting the efficient computation of Slepian functions for
the polar cap region on the sphere, we develop a formulation, supported by a
fast algorithm, for the approximate computation of Slepian functions for an
arbitrary spatial region to enable the analysis of modern datasets that support
large bandlimits. For the proposed algorithm, we carry out accuracy analysis
of the approximation, computational complexity analysis, and review of memory
storage requirements. We illustrate, through numerical experiments, that the
proposed method enables faster computation, and has smaller storage
requirements, while allowing for sufficiently accurate computation of the
Slepian functions.

In this paper, we develop an analytical formulation for the Slepian
spatialspectral concentration problem on the sphere for a limited
colatitudelongitude spatial region on the sphere, defined as the Cartesian
product of a range of positive colatitudes and longitudes. The solution of the
Slepian problem is a set of functions that are optimally concentrated and
orthogonal within a spatial or spectral region. These properties make them
useful for applications where measurements are taken within a spatially limited
region of the sphere and/or a signal is only to be analyzed within a region of
the sphere. To support localized spectral/spatial analysis, and estimation and
sparse representation of localized data in these applications, we exploit the
expansion of spherical harmonics in the complex exponential basis to develop an
analytical formulation for the Slepian concentration problem for a limited
colatitudelongitude spatial region. We also extend the analytical formulation
for spatial regions that are comprised of a union of rotated limited
colatitudelongitude subregions. By exploiting various symmetries of the
proposed formulation, we design a computationally efficient algorithm for the
implementation of the proposed analytical formulation. Such a reduction in
computation time is demonstrated through numerical experiments. We present
illustrations of our results with the help of numerical examples and show that
the representation of a spatially concentrated signal is indeed sparse in the
Slepian basis.

This paper proposes a multishell sampling scheme and corresponding
transforms for the accurate reconstruction of the diffusion signal in diffusion
MRI by expansion in the spherical polar Fourier (SPF) basis. The sampling
scheme uses an optimal number of samples, equal to the degrees of freedom of
the bandlimited diffusion signal in the SPF domain, and allows for
computationally efficient reconstruction. We use synthetic data sets to
demonstrate that the proposed scheme allows for greater reconstruction accuracy
of the diffusion signal than the multishell sampling schemes obtained using
the generalised electrostatic energy minimisation (gEEM) method used in the
Human Connectome Project. We also demonstrate that the proposed sampling scheme
allows for increased angular discrimination and improved rotational invariance
of reconstruction accuracy than the gEEM schemes.

We propose a multishell sampling grid and develop corresponding transforms
for the accurate reconstruction of the diffusion signal in diffusion MRI by
expansion in the spherical polar Fourier (SPF) basis. The transform is exact in
the radial direction and accurate, on the order of machine precision, in the
angular direction. The sampling scheme uses an optimal number of samples equal
to the degrees of freedom of the diffusion signal in the SPF domain.

The system that we study in this paper contains a set of users that observe a
discrete memoryless multiple source and communicate via noisefree channels
with the aim of attaining omniscience, the state that all users recover the
entire multiple source. We adopt the concept of successive omniscience (SO),
i.e., letting the local omniscience in some user subset be attained before the
global omniscience in the entire system, and consider the problem of how to
efficiently attain omniscience in a successive manner. Based on the existing
results on SO, we propose a CompSetSO algorithm for determining a complimentary
set, a user subset in which the local omniscience can be attained first without
increasing the sumrate, the total number of communications, for the global
omniscience. We also derive a sufficient condition for a user subset to be
complimentary so that running the CompSetSO algorithm only requires a lower
bound, instead of the exact value, of the minimum sumrate for attaining global
omniscience. The CompSetSO algorithm returns a complimentary user subset in
polynomial time. We show by example how to recursively apply the CompSetSO
algorithm so that the global omniscience can be attained by multistages of SO.

We propose a GaussLegendre quadrature based sampling on the rotation group
for the representation of a bandlimited signal such that the Fourier transform
(FT) of a signal can be exactly computed from its samples. Our figure of merit
is the sampling efficiency, which is defined as a ratio of the degrees of
freedom required to represent a bandlimited signal in harmonic domain to the
number of samples required to accurately compute the FT. The proposed sampling
scheme is asymptotically as efficient as the most efficient scheme developed
very recently. For the computation of FT and inverse FT, we also develop fast
algorithms of complexity similar to the complexity attained by the fast
algorithms for the existing sampling schemes. The developed algorithms are
stable, accurate and do not have any precomputation requirements. We also
analyse the computation time and numerical accuracy of the proposed algorithms
and show, through numerical experiments, that the proposed Fourier transforms
are accurate with errors on the order of numerical precision.

We develop a method for the accurate reconstruction of nonbandlimited finite
rate of innovation signals on the sphere. For signals consisting of a finite
number of Dirac functions on the sphere, we develop an annihilating filter
based method for the accurate recovery of parameters of the Dirac functions
using a finite number of observations of the bandlimited signal. In comparison
to existing techniques, the proposed method enables more accurate
reconstruction primarily due to better conditioning of systems involved in the
recovery of parameters. For the recovery of $K$ Diracs on the sphere, the
proposed method requires samples of the signal bandlimited in the spherical
harmonic~(SH) domain at SH degree equal or greater than $ K + \sqrt{K +
\frac{1}{4}}  \frac{1}{2}$. In comparison to the existing stateofthe art
technique, the required bandlimit, and consequently the number of samples, of
the proposed method is the same or less. We also conduct numerical experiments
to demonstrate that the proposed technique is more accurate than the existing
methods by a factor of $10^{7}$ or more for $2 \le K\le 20$.

This paper studies a wireless network consisting of multiple
transmitterreceiver pairs sharing the same spectrum where interference is
regarded as noise. Previously, the throughput region of such a network was
characterized for either one time slot or an infinite time horizon. This work
aims to close the gap by investigating the throughput region for transmissions
over a finite time horizon. We derive an efficient algorithm to examine the
achievability of any given rate in the finitehorizon throughput region and
provide the rateachieving policy. The computational efficiency of our
algorithm comes from the use of A* search with a carefully chosen heuristic
function and a tree pruning strategy. We also show that the celebrated
maxweight algorithm which finds all achievable rates in the infinitehorizon
throughput region fails to work for the finitehorizon throughput region.

We propose a novel algorithm to design user loadachieving pilot sequences
that mitigate pilot contamination in multicell massive multipleinput
multipleoutput (MIMO) networks. To this end, we first derive expressions for
the user load and the load region of the network considering both smallscale
and largescale propagation effects. We then develop the pilot sequence
algorithm for multicell massive MIMO networks as per the rules of generalized
Welch bound equality design. Notably, we find that our algorithm and the
corresponding downlink power allocation ensure that the user load is achieved
when the signaltointerferenceplusnoise ratio (SINR) requirements for the
users lie within the load region. Furthermore, we demonstrate the performance
advantage of our proposed design relative to the existing designs, in terms of
a larger load region and a higher maximum permitted SINR. Finally, we show that
our proposed design can satisfy the predefined SINR requirements for users
with a finite number of antennas at the base station (BS), while the existing
designs cannot satisfy the same requirements even with an infinite number of
antennas at the BS.

This paper studies a wireless network consisting of multiple
transmitterreceiver pairs where interference is treated as noise. Previously,
the throughput region of such networks was characterized for either one time
slot or an infinite time horizon. We aim to fill the gap by investigating the
throughput region for transmissions over a finite time horizon. Unlike the
infinitehorizon throughput region, which is simply the convex hull of the
throughput region of one time slot, the finitehorizon throughput region is
generally nonconvex. Instead of directly characterizing all achievable
ratetuples in the finitehorizon throughput region, we propose a metric termed
the rate margin, which not only determines whether any given ratetuple is
within the throughput region (i.e., achievable or unachievable), but also tells
the amount of scaling that can be done to the given achievable (unachievable)
ratetuple such that the resulting ratetuple is still within (brought back
into) the throughput region. Furthermore, we derive an efficient algorithm to
find the rateachieving policy for any given ratetuple in the finitehorizon
throughput region.

We propose a novel pilot sequence design to mitigate pilot contamination in
multicell multiuser massive multipleinput multipleoutput networks. Our
proposed design generates pilot sequences in the multicell network and devises
power allocation at base stations (BSs) for downlink transmission. The pilot
sequences together with the power allocation ensure that the user capacity of
the network is achieved and the predefined signaltointerferenceplusnoise
ratio (SINR) requirements of all users are met. To realize our design, we first
derive new closedform expressions for the user capacity and the user capacity
region. Built upon these expressions, we then develop a new algorithm to obtain
the required pilot sequences and power allocation. We further determine the
minimum number of antennas required at BSs to achieve certain SINR requirements
of all users. Numerical results are presented to corroborate our analysis and
to examine the impact of key parameters, such as the pilot sequence length and
the total number of users, on the network performance. A pivotal conclusion is
reached that our design achieves a larger user capacity region than the
existing designs and needs less antennas at the BS to fulfill the predefined
SINR requirements of all users in the network than the existing designs.

We propose a modified decomposition algorithm (MDA) to solve the asymptotic
communication for omniscience (CO) problem where the communication rates could
be real or fractional. By starting with a lower estimation of the minimum
sumrate, the MDA algorithm iteratively updates the estimation by the optimizer
of a Dilworth truncation problem until the minimum is reached with a
corresponding optimal rate vector. We also propose a fusion method
implementation of the coordinatewise saturation capacity algorithm
(CoordSatCapFus) for solving the Dilworth truncation problem, where the
minimization is done over a fused user set with a cardinality smaller than the
original one. We show that the MDA algorithm is less complex than the existing
ones. In addition, we show that the nonasymptotic CO problem, where the
communication rates are integral, can be solved by one more call of the
CoordSatCapfus algorithm. By choosing a proper linear ordering of the user
indices in the MDA algorithm, the optimal rate vector is also the one with the
minimum weighted sumrate.

We study the adaptive modulation (AM) problem in a networkcoded twoway
relay channel (NCTWRC), where each of the two users controls its own bit rate
in the $m$ary quadrature amplitude modulation ($m$QAM) to minimize the
transmission error rate and enhance the spectral efficiency. We show that there
exists a strategic complementarity, one user tends to transmit while the other
decides to do so in order to enhance the overall spectral efficiency, which is
beyond the scope of the conventional singleagent AM scheduling method. We
propose a twoplayer game model parameterized by the signaltonoise ratios
(SNRs) of two usertouser channels and prove that it is a supermodular game
where there always exist the extremal pure strategy Nash equilibria (PSNEs),
the largest and smallest PSNEs. We show by simulation results that the extremal
PSNEs incur a similar bit error rate (BER) as the conventional singleagent AM
scheme, but significantly improve the spectral efficiency in the NCTWRC
system. The study also reveals the Pareto order of the extremal PSNEs: The
largest and smallest PSNEs are Pareto worst and best PSNEs, respectively.
Finally, we derive the sufficient conditions for the extremal PSNEs to be
symmetric and monotonic in channel SNRs. We also discuss how to utilize the
symmetry and monotonicity to relieve the complexity in the PSNE learning
process.

In this paper, a cooperative diffusionbased molecular communication system
is considered where distributed receivers collaboratively determine a
transmitter's signal. In this system, the receivers first make local hard
decisions about the current transmitted bit and then report these decisions to
a fusion center (FC). The FC combines the local hard decisions to make a global
decision using an NoutofK fusion rule. Asymmetric and symmetric topologies
are considered and for each topology, two reporting scenarios, namely, perfect
reporting and noisy reporting, are addressed. Closedform analytical
expressions for the expected global error probability are derived for all
considered topologies and scenarios. Numerical and simulation results show that
system reliability can be greatly improved by combining the detection
information of distributed receivers.

We consider the problem of how to fairly distribute the minimum sumrate
among the users in communication for omniscience (CO). We formulate a problem
of minimizing a weighted quadratic function over a submodular base polyhedron
which contains all achievable rate vectors, or transmission strategies, for CO
that have the same sumrate. By solving it, we can determine the rate vector
that optimizes the Jain's fairness measure, a more commonly used fairness index
than the Shapley value in communications engineering. We show that the
optimizer is a lexicographically optimal (lexoptimal) base and can be
determined by a decomposition algorithm (DA) that is based on submodular
function minimization (SFM) algorithm and completes in strongly polynomial
time. We prove that the lexoptimal minimum sumrate strategy for CO can be
determined by finding the lexoptimal base in each user subset in the
fundamental partition and the complexity can be reduced accordingly.

We propose a coalition game model for the problem of communication for
omniscience (CO). In this game model, the core contains all achievable rate
vectors for CO with sumrate being equal to a given value. Any rate vector in
the core distributes the sumrate among users in a way that makes all users
willing to cooperate in CO. We give the necessary and sufficient condition for
the core to be nonempty. Based on this condition, we derive the expression of
the minimum sumrate for CO and show that this expression is consistent with
the results in multivariate mutual information (MMI) and coded cooperative data
exchange (CCDE). We prove that the coalition game model is convex if the
sumrate is no less than the minimal value. In this case, the core is nonempty
and a rate vector in the core that allocates the sumrate among the users in a
fair manner can be found by calculating the Shapley value.

This paper considers the discrete convexity of a crosslayer onoff
transmission control problem in wireless communications. In this system, a
scheduler decides whether or not to transmit in order to optimize the longterm
quality of service (QoS) incurred by the queueing effects in the data link
layer and the transmission power consumption in the physical (PHY) layer
simultaneously. Using a Markov decision process (MDP) formulation, we show that
the optimal policy can be determined by solving a minimization problem over a
set of queue thresholds if the dynamic programming (DP) is submodular. We prove
that this minimization problem is discrete convex. In order to search the
minimizer, we consider two discrete stochastic approximation (DSA) algorithms:
discrete simultaneous perturbation stochastic approximation (DSPSA) and
Lnaturalconvex stochastic approximation (Lnaturalconvex SA). Through
numerical studies, we show that the two DSA algorithms converge significantly
faster than the existing continuous simultaneous perturbation stochastic
approximation (CSPSA) algorithm in multiuser systems. Finally, we compare the
convergence results and complexity of two DSA and CSPSA algorithms where we
show that DSPSA achieves the best tradeoff between complexity and accuracy in
multiuser systems.

This paper considers a crosslayer adaptive modulation system that is modeled
as a Markov decision process (MDP). We study how to utilize the monotonicity of
the optimal transmission policy to relieve the computational complexity of
dynamic programming (DP). In this system, a scheduler controls the bit rate of
the mquadrature amplitude modulation (mQAM) in order to minimize the
longterm losses incurred by the queue overflow in the data link layer and the
transmission power consumption in the physical layer. The work is done in two
steps. Firstly, we observe the Lnaturalconvexity and submodularity of DP to
prove that the optimal policy is always nondecreasing in queue occupancy/state
and derive the sufficient condition for it to be nondecreasing in both queue
and channel states. We also show that, due to the Lnaturalconvexity of DP,
the variation of the optimal policy in queue state is restricted by a bounded
marginal effect: The increment of the optimal policy between adjacent queue
states is no greater than one. Secondly, we use the monotonicity results to
present two low complexity algorithms: monotonic policy iteration (MPI) based
on Lnaturalconvexity and discrete simultaneous perturbation stochastic
approximation (DSPSA). We run experiments to show that the time complexity of
MPI based on Lnaturalconvexity is much lower than that of DP and the
conventional MPI that is based on submodularity and DSPSA is able to adaptively
track the optimal policy when the system parameters change.

This paper considers how to accurately estimate the minimum sumrate so as to
reduce the complexity of solving cooperative data exchange (CDE) problems. The
CDE system contains a number of geographically close clients who send packets
to help the others recover an entire packet set. The minimum sumrate is the
minimum value of total number of transmissions that achieves universal recovery
(the situation when all the clients recover the whole packet set). Based on a
necessary and sufficient condition for a supermodular base polyhedron to be
nonempty, we show that the minimum sumrate for a CDE system can be determined
by a maximization over all possible partitions of the client set. Due to the
high complexity of solving this maximization problem, we propose a
deterministic algorithm to approximate a lower bound on the minimum sumrate.
We show by experiments that this lower bound is much tighter than those lower
bounds derived in the existing literature. We also show that the deterministic
algorithm prevents from repetitively running the existing algorithms for
solving CDE problems so that the overall complexity can be reduced accordingly.