BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Jun Fang

Jun Fang appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
Information Theorymath.ITastro-ph.HEeess.SPMachine Learningphysics.chem-phphysics.soc-phInformation Retrievalmath.NANumerical Analysisphysics.comp-phArtificial IntelligenceComputer Visioncond-mat.stat-mechcond-mat.str-elnlin.CG

Catalog footprint

What is connected

59works
16topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

59 published item(s)

preprint2026arXiv

D-Artemis: A Deliberative Cognitive Framework for Mobile GUI Multi-Agents

Graphical User Interface (GUI) agents aim to automate a wide spectrum of human tasks by emulating user interaction. Despite rapid advancements, current approaches are hindered by several critical challenges: data bottleneck in end-to-end training, high cost of delayed error detection, and risk of contradictory guidance. Inspired by the human cognitive loop of Thinking, Alignment, and Reflection, we present D-Artemis -- a novel deliberative framework in this paper. D-Artemis leverages a fine-grained, app-specific tip retrieval mechanism to inform its decision-making process. It also employs a proactive Pre-execution Alignment stage, where Thought-Action Consistency (TAC) Check module and Action Correction Agent (ACA) work in concert to mitigate the risk of execution failures. A post-execution Status Reflection Agent (SRA) completes the cognitive loop, enabling strategic learning from experience. Crucially, D-Artemis enhances the capabilities of general-purpose Multimodal large language models (MLLMs) for GUI tasks without the need for training on complex trajectory datasets, demonstrating strong generalization. D-Artemis establishes new state-of-the-art (SOTA) results across both major benchmarks, achieving a 75.8% success rate on AndroidWorld and 96.8% on ScreenSpot-V2. Extensive ablation studies further demonstrate the significant contribution of each component to the framework.

0 citations0 reviewsNo code linkedOpen access
preprint2025arXiv

Ultrahigh-Energy Gamma-ray Emission Associated with Black Hole-Jet Systems

Black holes (BH), one of the most intriguing objects in the universe, can manifest themselves through electromagnetic radiation initiated by the accretion flow. Some stellar-mass BHs drive relativistic jets when accreting matter from their companion stars, forming microquasars. Non-thermal emission from the radio to tera-electronvolt (TeV) gamma-ray band has been observed from microquasars, indicating the acceleration of relativistic particles. Here we report detection of four microquasars (SS 433, V4641 Sgr, GRS 1915+105, MAXI J1820+070) of spectrum extending to the ultrahigh-energy (UHE; photon energy $E>100$ TeV) band and one microquasar (Cygnus X-1) of spectrum approaching 100 TeV, using the Large High Altitude Air Shower Observatory (LHAASO). Notably, the total emission associated with SS 433 cannot be interpreted with a single leptonic component. In the UHE band, its emission is in spatial coincidence with a giant atomic cloud, which is consistent with a hadronic origin. An elongated source is discovered from V4641 Sgr with the spectrum continuing up to 800 TeV. The detection of UHE gamma rays demonstrates that accreting BHs and their environments can operate as extremely efficient accelerators of particles out of 1 peta-electronvolt (PeV), suggesting microquasars to be important contributors to Galactic cosmic rays especially around the `knee' region.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

A Variational Bayesian Inference-Inspired Unrolled Deep Network for MIMO Detection

The great success of deep learning (DL) has inspired researchers to develop more accurate and efficient symbol detectors for multi-input multi-output (MIMO) systems. Existing DL-based MIMO detectors, however, suffer several drawbacks. To address these issues, in this paper, we develop a model-driven DL detector based on variational Bayesian inference. Specifically, the proposed unrolled DL architecture is inspired by an inverse-free variational Bayesian learning framework which circumvents matrix inversion via maximizing a relaxed evidence lower bound. Two networks are respectively developed for independent and identically distributed (i.i.d.) Gaussian channels and arbitrarily correlated channels. The proposed networks, referred to as VBINet, have only a few learnable parameters and thus can be efficiently trained with a moderate amount of training samples. The proposed VBINet-based detectors can work in both offline and online training modes. An important advantage of our proposed networks over state-of-the-art MIMO detection networks such as OAMPNet and MMNet is that the VBINet can automatically learn the noise variance from data, thus yielding a significant performance improvement over the OAMPNet and MMNet in the presence of noise variance uncertainty. Simulation results show that the proposed VBINet-based detectors achieve competitive performance for both i.i.d. Gaussian and realistic 3GPP MIMO channels.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Beam Training and Alignment for RIS-Assisted Millimeter Wave Systems:State of the Art and Beyond

Reconfigurable intelligent surface (RIS) has recently emerged as a promising paradigm for future cellular networks. Specifically, due to its capability in reshaping the propagation environment, RIS was introduced to address the blockage issue in millimeter Wave (mmWave) or even Terahertz (THz) communications. The deployment of RIS, however, complicates the system architecture and poses a significant challenge for beam training (BT)/ beam alignment (BA), a process that is required to establish a reliable link between the transmitter and the receiver. In this article, we first review several state-of-the-art beam training solutions for RIS-assisted mmWave systems and discuss their respective advantages and limitations. We also present a new multi-directional BT method, which can achieve a decent BA performance with only a small amount of training overhead. Finally, we outline several important open issues in BT for RIS-assisted mmWave systems.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Compressed Channel Estimation for IRS-Assisted Millimeter Wave OFDM Systems: A Low-Rank Tensor Decomposition-Based Approach

We consider the problem of downlink channel estimation for intelligent reflecting surface (IRS)-assisted millimeter Wave (mmWave) orthogonal frequency division multiplexing (OFDM) systems. By exploring the inherent sparse scattering characteristics of mmWave channels, we show that the received signals can be expressed as a low-rank third-order tensor that admits a tensor rank decomposition, also known as canonical polyadic decomposition (CPD). A structured CPD-based method is then developed to estimate the channel parameters. Our analysis reveals that the training overhead required by our proposed method is as low as O(U^2), where U denotes the sparsity of the cascade channel. Simulation results are provided to illustrate the efficiency of the proposed method.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Investigating the radiative properties of LHAASO J1908+0621

LHAASO J1908+0621 has been recently detected as a source emitting $γ$-rays with energies above 100 TeV, and the multiband observations show that a break around 1 TeV appears in the $γ$-ray spectrum. We have reanalyzed the GeV $γ$-ray properties for the 100 TeV source using 14 years of data recorded by the Fermi Large Area Telescope (Fermi-LAT). The spectrum in the energy range range 30-500 GeV has an index of 1.50 $\pm$ 0.26, which is much smaller than that detected in the TeV $γ$-rays. Additionally, the radiation properties of this source are investigated based on a one-zone time-dependent model. In the model, LHAASO J1908+0621 is associated with a pulsar wind nebula (PWN) powered by the pulsar PSR J1907$+$0602. High-energy particles composed of electrons and positrons are injected into the nebula. Multiband nonthermal emission is produced via synchrotron radiation and inverse Compton scattering (ICS). Taking the effect of radiative energy losses and adiabatic cooling into account, the spectral energy distribution from the model with a broken power-law for the distribution of the injected particles can explain the detected fluxes in the $γ$-ray bands. The results support that LHAASO J1908+0621 originates from the PWN powered by PSR J1907$+$0602, and the $γ$-rays with energy above 100 TeV are produced by the electrons/positrons in the nebula via ICS.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Kinetic Simulation on Electron, Proton and Helium Acceleration in a Nonrelativistic Quasiparallel shock

In addition to electrons and protons, nonrelativistic quasiparallel shocks are expected to possess the ability to accelerate heavy ions. The shocks in supernova remnants are generally supposed to be accelerators of the Galactic cosmic rays, which consist of many species of particles. We investigate diffusive shock acceleration (DSA) of electrons, protons and helium ions in a nonrelativistic quasiparallel shock through 1D particle-in-cell (PIC) simulation with a helium-to-proton number density ratio of $0.1$, which is relevant for the Galactic cosmic rays. The simulation indicates that waves can be excited by the flow of the energetic protons and helium ions upstream of the nonrelativistic quasiparallel shock with a sonic Mach number of 14 and an alfvén Mach number of 19.5 in the shock rest frame, and the charged particles are scattered by the self-generated waves and accelerated gradually. Moreover, the spectra of the charged particles downstream of the shock are thermal plus a nonthermal tail, and the acceleration is efficient with about $7\%$ and $5.4\%$ of the bulk kinetic energy transferred into the nonthermal protons and helium ions in the near downstream region at the end of the simulation, respectively.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Over-the-Air Federated Multi-Task Learning via Model Sparsification and Turbo Compressed Sensing

To achieve communication-efficient federated multitask learning (FMTL), we propose an over-the-air FMTL (OAFMTL) framework, where multiple learning tasks deployed on edge devices share a non-orthogonal fading channel under the coordination of an edge server (ES). In OA-FMTL, the local updates of edge devices are sparsified, compressed, and then sent over the uplink channel in a superimposed fashion. The ES employs over-the-air computation in the presence of intertask interference. More specifically, the model aggregations of all the tasks are reconstructed from the channel observations concurrently, based on a modified version of the turbo compressed sensing (Turbo-CS) algorithm (named as M-Turbo-CS). We analyze the performance of the proposed OA-FMTL framework together with the M-Turbo-CS algorithm. Furthermore, based on the analysis, we formulate a communication-learning optimization problem to improve the system performance by adjusting the power allocation among the tasks at the edge devices. Numerical simulations show that our proposed OAFMTL effectively suppresses the inter-task interference, and achieves a learning performance comparable to its counterpart with orthogonal multi-task transmission. It is also shown that the proposed inter-task power allocation optimization algorithm substantially reduces the overall communication overhead by appropriately adjusting the power allocation among the tasks.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Quasi-periodic behaviour in the $γ$-ray light curve of the blazar PKS 0405-385

We analyze the quasi-periodic oscillation (QPO) of the historical light curve of FSRQs PKS 0405-385 detected by the Fermi LAT from August 2008 to November 2021. To identify and determine the QPO signal of PKS 0405-385 in the $γ$-ray light curve, we use four time series analysis techniques based on frequency and time domains, i.e., the Lomb-Scargle periodogram (LSP), the weighted wavelet z-transform (WWZ), the REDFIT and the epoch folding. The results show that PKS 0405-385 has a quasi-periodic behavior of $\sim$2.8 yr with the significance of $\sim$4.3$σ$ in Fermi long-term monitoring. Remarkably, we also performed QPO analysis in the G-band light curve observed from October 2014 to October 2021 using LSP and WWZ technology, and the results ($\sim$4$σ$ of significance) are consistent with the periodic detection in $γ$-ray. This may imply that the optical emission is radiated by an electron population same as the \gr\ emission. In discussing the possible mechanism of quasi-periodic behavior, either the helical motion within a jet or the supermassive black hole binary system provides a viable explanation for the QPO of 2.8 yr, and the relevant parameters have been estimated.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Spatial Channel Covariance Estimation and Two-Timescale Beamforming for IRS-Assisted Millimeter Wave Systems

We consider the problem of spatial channel covariance matrix (CCM) estimation for intelligent reflecting surface (IRS)-assisted millimeter wave (mmWave) communication systems. Spatial CCM is essential for two-timescale beamforming in IRS-assisted systems; however, estimating the spatial CCM is challenging due to the passive nature of reflecting elements and the large size of the CCM resulting from massive reflecting elements of the IRS. In this paper, we propose a CCM estimation method by exploiting the low-rankness as well as the positive semi-definite (PSD) 3-level Toeplitz structure of the CCM. Estimation of the CCM is formulated as a semidefinite programming (SDP) problem and an alternating direction method of multipliers (ADMM) algorithm is developed. Our analysis shows that the proposed method is theoretically guaranteed to attain a reliable CCM estimate with a sample complexity much smaller than the dimension of the CCM. Thus the proposed method can help achieve a significant training overhead reduction. Simulation results are presented to illustrate the effectiveness of our proposed method and the performance of two-timescale beamforming scheme based on the estimated CCM.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Investigating the energy distribution of the high-energy particles in the Crab nebula

The Crab nebula is a prominent pulsar wind nebula (PWN) detected in multiband observations ranging from radio to very high-energy (VHE) $γ$-rays. Recently, $γ$-rays with energies above $1 \mathrm{PeV}$ had been detected by the Large High Altitude Air Shower Observatory (LHAASO), and the energy of the most energetic particles in the nebula can be constrained. In this paper, we investigate the broadest spectral energy distribution of the Crab nebula and the energy distribution of the electrons emitting the multiwavelength nonthermal emission based on a one-zone time-dependent model. The nebula is powered by the pulsar, and high-energy electrons/positrons with a broken power-law spectrum are continually injected in the nebula as the pulsar spins down. Multiwavelength nonthermal emission is generated by the leptons through synchrotron radiation and inverse Compton scattering. Using appropriate parameters, the detected fluxes for the nebula can be well reproduced, especially for the $γ$-rays from $10^2\,\mathrm{MeV}$ to $1\,\mathrm{PeV}$. The results show that the detected $γ$-rays can be produced by the leptons via the inverse Compton scattering, and the lower limit of the Lorentz factor of the most energetic leptons is $\sim 8.5\times10^{9}$. It can be concluded that there are electrons/positrons with energies higher than $4.3$\,PeV in the Crab nebula.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Joint Active and Passive Beamforming for IRS-Assisted Radar

Intelligent reflecting surface (IRS) is a promising technology being considered for future wireless communications due to its ability to control signal propagation. This paper considers the joint active and passive beamforming problem for an IRS-assisted radar, where multiple IRSs are deployed to assist the surveillance of multiple targets in cluttered environments. Specifically, we aim to maximize the minimum target illumination power at multiple target locations by jointly optimizing the active beamformer at the radar transmitter and the passive phase-shift matrices at the IRSs, subject to an upperbound on the clutter power at each clutter scatterer. The resulting optimization problem is nonconvex and solved with a sequential optimization procedure along with semedefinite relaxation (SDR). Simulation results show that IRSs can help create effective line-of-sight (LOS) paths and thus substantially improve the radar robustness against target blockage.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Session-based Recommendation with Self-Attention Networks

Session-based recommendation aims to predict user's next behavior from current session and previous anonymous sessions. Capturing long-range dependencies between items is a vital challenge in session-based recommendation. A novel approach is proposed for session-based recommendation with self-attention networks (SR-SAN) as a remedy. The self-attention networks (SAN) allow SR-SAN capture the global dependencies among all items of a session regardless of their distance. In SR-SAN, a single item latent vector is used to capture both current interest and global interest instead of session embedding which is composed of current interest embedding and global interest embedding. Some experiments have been performed on some open benchmark datasets. Experimental results show that the proposed method outperforms some state-of-the-arts by comparisons.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Terahertz Multi-User Massive MIMO with Intelligent Reflecting Surface: Beam Training and Hybrid Beamforming

Terahertz (THz) communications open a new frontier for the wireless network thanks to their dramatically wider available bandwidth compared to the current micro-wave and forthcoming millimeter-wave communications. However, due to the short length of THz waves, they also suffer from severe path attenuation and poor diffraction. To compensate the THz-induced propagation loss, this paper proposes to combine two promising techniques, viz., massive multiple input multiple output (MIMO) and intelligent reflecting surface (IRS), in THz multi-user communications, considering their significant beamforming and aperture gains. Nonetheless, channel estimation and low-cost beamforming turn out to be two main obstacles to realizing this combination, due to the passivity of IRS for sending/receiving pilot signals and the large-scale use of expensive RF chains in massive MIMO. In view of these limitations, this paper first develops a cooperative beam training scheme to facilitate the channel estimation with IRS. In particular, we design two different hierarchical codebooks for the proposed training procedure, which are able to balance between the robustness against noise and searching complexity. Based on the training results, we further propose two cost-efficient hybrid beamforming (HB) designs for both single-user and multi-user scenarios, respectively. Simulation results demonstrate that the proposed joint beam training and HB scheme is able to achieve close performance to the optimal fully digital beamforming (FDB) which is implemented even under perfect channel state information (CSI).

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Beamforming Optimization for Intelligent Reflecting Surface Assisted MIMO: A Sum-Path-Gain Maximization Approach

Recently, intelligent reflecting surface (IRS) has emerged as an appealing technique that enables wireless communications with low hardware cost and low power consumption. In this letter, we consider an IRS-assisted point-to-point multi-input multi-output (MIMO) system, where a source communicates with its destination with the help of an IRS. Our goal is to maximize the spectral efficiency of this system by jointly optimizing the (active) precoding at the source and the (passive) phase shifters (PSs) at the IRS. However, this turns out to be an intractable mixed integer non-convex optimization problem. To circumvent the intractability, we propose a new sum-path-gain maximization (SPGM) criterion to obtain a high-quality and efficient suboptimal solution to this problem. Specifically, the PSs are first designed based on a simplified optimization problem, which aims to maximize the sum-gains of the spatial paths between the source and the destination. Then, a low-complexity alternating direction method of multipliers (ADMM) algorithm is utilized to solve this simplified problem. Finally, with the above obtained PSs, the source precoding is derived by performing the singular value decomposition (SVD) on the effective channel between the source and the destination. Numerical results demonstrate that the proposed scheme can achieve near-optimal performance.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Compressed Channel Estimation and Joint Beamforming for Intelligent Reflecting Surface-Assisted Millimeter Wave Systems

In this paper, we consider channel estimation for intelligent reflecting surface (IRS)-assisted millimeter wave (mmWave) systems, where an IRS is deployed to assist the data transmission from the base station (BS) to a user. It is shown that for the purpose of joint active and passive beamforming, the knowledge of a large-size cascade channel matrix needs to be acquired. To reduce the training overhead, the inherent sparsity in mmWave channels is exploited. By utilizing properties of Katri-Rao and Kronecker products, we find a sparse representation of the cascade channel and convert cascade channel estimation into a sparse signal recovery problem. Simulation results show that our proposed method can provide an accurate channel estimate and achieve a substantial training overhead reduction.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Post-Training Piecewise Linear Quantization for Deep Neural Networks

Quantization plays an important role in the energy-efficient deployment of deep neural networks on resource-limited devices. Post-training quantization is highly desirable since it does not require retraining or access to the full training dataset. The well-established uniform scheme for post-training quantization achieves satisfactory results by converting neural networks from full-precision to 8-bit fixed-point integers. However, it suffers from significant performance degradation when quantizing to lower bit-widths. In this paper, we propose a piecewise linear quantization (PWLQ) scheme to enable accurate approximation for tensor values that have bell-shaped distributions with long tails. Our approach breaks the entire quantization range into non-overlapping regions for each tensor, with each region being assigned an equal number of quantization levels. Optimal breakpoints that divide the entire range are found by minimizing the quantization error. Compared to state-of-the-art post-training quantization methods, experimental results show that our proposed method achieves superior performance on image classification, semantic segmentation, and object detection with minor overhead.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Reconfigurable Intelligent Surface Aided Constant-Envelope Wireless Power Transfer

By reconfiguring the propagation environment of electromagnetic waves artificially, reconfigurable intelligent surfaces (RISs) have been regarded as a promising and revolutionary hardware technology to improve the energy and spectrum efficiency of wireless networks. In this paper, we study a RIS aided multiuser multiple-input single-output (MISO) wireless power transfer (WPT) system, where the transmitter is equipped with a constant-envelope analog beamformer. We formulate a novel problem to maximize the total received power of all the users by jointly optimizing the beamformer at transmitter and the phase shifts at the RISs, subject to the individual minimum received power constraints of users. We further solve the problem iteratively with a closed-form expression for each step. Numerical results show the performance gain of deploying RIS and the effectiveness of the proposed algorithm.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Early acceleration of electrons and protons at the nonrelativistic quasiparallel shocks with different obliquity angles

The early acceleration of protons and electrons in the nonrelativistic collisionless shocks with three obliquities are investigated through 1D particle-in-cell simulations. In the simulations, the charged particles possessing a velocity of $0.2\, c$ flow towards a reflecting boundary, and the shocks with a sonic Mach number of $13.4$ and a Alfven Mach number of $16.5$ in the downstream shock frame are generated. In these quasi-parallel shocks with the obliquity angles $θ= 15^\circ$, $30^\circ$, and $45^\circ$, some of the protons and the electrons can be injected into the acceleration processes, and their downstream spectra in the momentum space show a power law tail at a time of $1.89\times10^5 ω_{\rm pe}^{-1}$, where $ω_{\rm pe}$ is the electron plasma frequency. Moreover, the charged particles reflected at the shock excite magnetic waves upstream of the shock. The shock drift acceleration is more prominent with a larger obliquity angle for the shocks, but the accelerated particles diffuse parallel to the shock propagation direction more easily to participate in the diffusive shock acceleration. At the time still in the early acceleration stage, more energetic protons and electrons appear in the downstream of the shock for $θ= 15^\circ$ compared with the other two obliquities; moreover, in the upstream region, the spectrum of the accelerated electrons is the hardest for $θ_{\rm nB} = 45^\circ$ among the three obliquities, whereas the proton spectra for $θ_{\rm nB} = 15^\circ$ and $45^\circ$ are similar as a result of the competition of the effectiveness of the shock drift acceleration and the diffusive shock acceleration.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Numerically investigating the morphology of the supernova remnant SN 1006 in the ambient medium with a density discontinuity

Multiband observations on the type Ia supernova remnant SN 1006 indicate peculiar properties in its morphologies of emission in the radio, optical and X-ray bands. In the hard X-rays, the remnant is bilateral with two opposite bright limbs with prominent protrusions. Moreover, a filament has been detected at the radio, optical and soft X-ray wavelengthes. The reason for these peculiar features in the morphologies of the remnant is investigated using 3D HD simulations. With the assumption that the supernova ejecta is evolved in the ambient medium with a density discontinuity, the radius of the remnant's boundary is smaller in the tenuous medium, and the shell consists of two hemispheres with different radiuses. Along particular line of sights, protrusions appear on the periphery of the remnants since the emission from the edge of the hemisphere with a larger radius is located outside that from the shell of the small hemisphere. Furthermore, the northwest filament of SN 1006 arises as a result of the intersection of the line of sight and the shocked material near the edges of the two hemispheres. It can be concluded that the features that the protrusions on the northeast and southwest limbs and the northwest filament in the morphologies of SN 1006 can be reproduced as the remnants interacting with the medium with a density discontinuity.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Channel Estimation for Millimeter Wave Multiuser MIMO Systems via PARAFAC Decomposition

We consider the problem of uplink channel estimation for millimeter wave (mmWave) systems, where the base station (BS) and mobile stations (MSs) are equipped with large antenna arrays to provide sufficient beamforming gain for outdoor wireless communications. Hybrid analog and digital beamforming structures are employed by both the BS and the MS due to hardware constraints. We propose a layered pilot transmission scheme and a CANDECOMP/PARAFAC (CP) decomposition-based method for joint estimation of the channels from multiple users (i.e. MSs) to the BS. The proposed method exploits the sparse scattering nature of the mmWave channel and the intrinsic multi-dimensional structure of the multiway data collected from multiple modes. The uniqueness of the CP decomposition is studied and sufficient conditions for essential uniqueness are obtained. The conditions shed light on the design of the beamforming matrix, the combining matrix and the pilot sequences, and meanwhile provide general guidelines for choosing system parameters. Our analysis reveals that our proposed method can achieve a substantial training overhead reduction by employing the layered pilot transmission scheme. Simulation results show that the proposed method presents a clear advantage over a compressed sensing-based method in terms of both estimation accuracy and computational complexity.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Dominant itinerant ferromagnetism in Eu0.5Sr0.5CoO3: evidences from both critical analysis and Rhodes-Wohlfarth's criterion

The critical behaviors of ferromagnet Eu0.5Sr0.5CoO3 arround TC=140.5K have been comprehensively investigated by analyzing a series of isothermal magnetization M(H) curves. Both Modified Arrott plot and Kouvel-Fisher methods give nearly the same critical exponents, which scale nicely the M(H) curves into two different branches below and above TC. The exponents γ=1.044 and δ=3.06 demonstrate the relevance of mean-field characters for this material. The conclusion of mean-field behavior proves a dominant itinerant ferromagnetism (FM) due to a long range exchange interaction in the system. Meanwhile, by using Rhodes-Wohlfarth's criterion [P. Rhodes and E. P. Wohlfarth, Proc. R. Soc. Lond. A 273, 247 (1963)], it is further confirmed that the itinerant FM dominates in the system.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

GSVD-Based Precoding in MIMO Systems With Integrated Services

This letter considers a two-receiver multiple-input multiple-output (MIMO) Gaussian broadcast channel model with integrated services. Specifically, we combine two sorts of service messages, and serve them simultaneously: one multicast message intended for both receivers and one confidential message intended for only one receiver. The confidential message is kept perfectly secure from the unauthorized receiver. Due to the coupling of service messages, it is intractable to seek capacity-achieving transmit covariance matrices. Accordingly, we propose a suboptimal precoding scheme based on the generalized singular value decomposition (GSVD). The GSVD produces several virtual orthogonal subchannels between the transmitter and the receivers. Subchannel allocation and power allocation between multicast message and confidential message are jointly optimized to maximize the secrecy rate in this letter, subject to the quality of multicast service (QoMS) constraints. Since this problem is inherently complex, a difference-of-concave (DC) algorithm, together with an exhaustive search, is exploited to handle the power allocation and subchannel allocation, respectively. Numerical results are presented to illustrate the efficacy of our proposed strategies.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Improving Wireless Physical Layer Security via Exploiting Co-Channel Interference

This paper considers a scenario in which a source-destination pair needs to establish a confidential connection against an external eavesdropper, aided by the interference generated by another source-destination pair that exchanges public messages. The goal is to compute the maximum achievable secrecy degrees of freedom (S.D.o.F) region of a MIMO two-user wiretap network. First, a cooperative secrecy transmission scheme is proposed, whose feasible set is shown to achieve all S.D.o.F. pairs on the S.D.o.F. region boundary. In this way, the determination of the S.D.o.F. region is reduced to a problem of maximizing the S.D.o.F. pair over the proposed transmission scheme. The maximum achievable S.D.o.F. region boundary points are obtained in closed form, and the construction of the precoding matrices achieving the maximum S.D.o.F. region boundary is provided. The obtained analytical expressions clearly show the relation between the maximum achievable S.D.o.F. region and the number of antennas at each terminal.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Learning Dominant Wave Directions For Plane Wave Methods For High-Frequency Helmholtz Equations

We present a ray-based finite element method (ray-FEM) by learning basis adaptive to the underlying high-frequency Helmholtz equation in smooth media. Based on the geometric optics ansatz of the wave field, we learn local dominant ray directions by probing the medium using low-frequency waves with the same source. Once local ray directions are extracted, they are incorporated into the finite element basis to solve the high-frequency Helmholtz equation. This process can be continued to further improve approximations for both local ray directions and the high frequency wave field iteratively. The method requires a fixed number of grid points per wavelength to represent the wave field and achieves an asymptotic convergence as the frequency $ω\rightarrow \infty$ without the pollution effect. A fast solver is developed for the resulting linear system with an empirical complexity $\mathcal{O}(ω^d)$ up to a poly-logarithmic factor. Numerical examples in 2D are presented to corroborate the claims.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Linear Precoder Design for a MIMO Gaussian Wiretap Channel with Full-Duplex Source and Destination Nodes

We consider linear precoder design for a multiple-input multiple-output (MIMO) Gaussian wiretap channel, which comprises two legitimate nodes, i.e., Alice and Bob, operating in Full-Duplex (FD) mode and exchanging confidential messages in the presence of a passive eavesdropper. Using the sum secrecy degrees of freedoms (sum S.D.o.F.) as reliability measure, we formulate an optimization problem with respect to the precoding matrices. In order to solve this problem, we first propose a cooperative secrecy transmission scheme, and prove that its feasible set is sufficient to achieve the maximum sum S.D.o.F.. Based on that feasible set, we then determine the maximum achievable sum S.D.o.F. in closed form, and provide a method for constructing the precoding matrix pair which achieves the maximum sum S.D.o.F.. Results show that, the FD based network provides an attractive secrecy transmission rate performance.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Low-Rank Covariance-Assisted Downlink Training and Channel Estimation for FDD Massive MIMO Systems

We consider the problem of downlink training and channel estimation in frequency division duplex (FDD) massive MIMO systems, where the base station (BS) equipped with a large number of antennas serves a number of single-antenna users simultaneously. To obtain the channel state information (CSI) at the BS in FDD systems, the downlink channel has to be estimated by users via downlink training and then fed back to the BS. For FDD large-scale MIMO systems, the overhead for downlink training and CSI uplink feedback could be prohibitively high, which presents a significant challenge. In this paper, we study the behavior of the minimum mean-squared error (MMSE) estimator when the channel covariance matrix has a low-rank or an approximate low-rank structure. Our theoretical analysis reveals that the amount of training overhead can be substantially reduced by exploiting the low-rank property of the channel covariance matrix. In particular, we show that the MMSE estimator is able to achieve exact channel recovery in the asymptotic low-noise regime, provided that the number of pilot symbols in time is no less than the rank of the channel covariance matrix. We also present an optimal pilot design for the single-user case, and an asymptotic optimal pilot design for the multi-user scenario. Lastly, we develop a simple model-based scheme to estimate the channel covariance matrix, based on which the MMSE estimator can be employed to estimate the channel. The proposed scheme does not need any additional training overhead. Simulation results are provided to verify our theoretical results and illustrate the effectiveness of the proposed estimated covariance-assisted MMSE estimator.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Low-Rank Tensor Decomposition-Aided Channel Estimation for Millimeter Wave MIMO-OFDM Systems

We consider the problem of downlink channel estimation for millimeter wave (mmWave) MIMO-OFDM systems, where both the base station (BS) and the mobile station (MS) employ large antenna arrays for directional precoding/beamforming. Hybrid analog and digital beamforming structures are employed in order to offer a compromise between hardware complexity and system performance. Different from most existing studies that are concerned with narrowband channels, we consider estimation of wideband mmWave channels with frequency selectivity, which is more appropriate for mmWave MIMO-OFDM systems. By exploiting the sparse scattering nature of mmWave channels, we propose a CANDECOMP/PARAFAC (CP) decomposition-based method for channel parameter estimation (including angles of arrival/departure, time delays, and fading coefficients). In our proposed method, the received signal at the BS is expressed as a third-order tensor. We show that the tensor has the form of a low-rank CP decomposition, and the channel parameters can be estimated from the associated factor matrices. Our analysis reveals that the uniqueness of the CP decomposition can be guaranteed even when the size of the tensor is small. Hence the proposed method has the potential to achieve substantial training overhead reduction. We also develop Cramer-Rao bound (CRB) results for channel parameters, and compare our proposed method with a compressed sensing-based method. Simulation results show that the proposed method attains mean square errors that are very close to their associated CRBs, and presents a clear advantage over the compressed sensing-based method in terms of both estimation accuracy and computational complexity.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Multiple Staggered Mesh Ewald: Boosting the Accuracy of the Smooth Particle Mesh Ewald Method

The smooth particle mesh Ewald (SPME) method is the standard method for computing the electrostatic interactions in the molecular simulations. In this work, the multiple staggered mesh Ewald (MSME) method is proposed to boost the accuracy of the SPME method. Unlike the SPME that achieves higher accuracy by refining the mesh, the MSME improves the accuracy by averaging the standard SPME forces computed on, e.g. $M$, staggered meshes. We prove, from theoretical perspective, that the MSME is as accurate as the SPME, but uses $M^2$ times less mesh points in a certain parameter range. In the complementary parameter range, the MSME is as accurate as the SPME with twice of the interpolation order. The theoretical conclusions are numerically validated both by a uniform and uncorrelated charge system, and by a three-point-charge water system that is widely used as solvent for the bio-macromolecules.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

On the existence of the optimal order for wavefunction extrapolation in Born-Oppenheimer molecular dynamics

Wavefunction extrapolation greatly reduces the number of self-consistent field (SCF) iterations and thus the overall computational cost of Born-Oppenheimer molecular dynamics (BOMD) that is based on the Kohn-Sham density functional theory. Going against the intuition that the higher order of extrapolation possesses a better accuracy, we demonstrate, from both theoretical and numerical perspectives, that the extrapolation accuracy firstly increases and then decreases with respect to the order, and an optimal extrapolation order in terms of minimal number of SCF iterations always exists. We also prove that the optimal order tends to be larger when using larger MD time steps or more strict SCF convergence criteria. By example BOMD simulations of a solid copper system, we show that the optimal extrapolation order covers a broad range when varying the MD time step or the SCF convergence criterion. Therefore, we suggest the necessity for BOMD simulation packages to open the user interface and to provide more choices on the extrapolation order. Another factor that may influence the extrapolation accuracy is the alignment scheme that eliminates the discontinuity in the wavefunctions with respect to the atomic or cell variables. We prove the equivalence between the two existing schemes, thus the implementation of either of them does not lead to essential difference in the extrapolation accuracy.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Perturbation-Assisted PAPR Reduction for Large-Scale MIMO-OFDM Systems via ADMM

We consider the problem of peak-to-average power ratio (PAPR) reduction for orthogonal frequency-division multiplexing (OFDM) based large-scale multiple-input multipleoutput (MIMO) systems. A novel perturbation-assisted scheme is developed to reduce the PAPRs of the transmitted signals by exploiting the redundant degrees-of-freedom (DoFs) inherent in the large-scale antenna array. Specifically, we introduce artificial perturbation signals to the frequency-domain precoded signals, with the aim of reducing the PAPRs of their time-domain counterpart signals. Meanwhile, the additive perturbation signal associated with each tone is constrained to lie in the null-space of its associated channel matrix, such that it does not cause any multi-user inference or out-of-band radiations. Such a problem is formulated as a convex optimization problem, and an efficient algorithm is developed by resorting to the variable splitting and alterative direction method of multipliers (ADMM) techniques. Simulation results show that the proposed method has a fast convergence rate and achieves substantial PAPR reduction within only tens of iterations. In addition, unlike other precoding-based PAPR reduction methods, our proposed method which introduces perturbation signals to the precoded signals is independent of the precoding stage and thus could be more suitable for practical systems.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Robust Bayesian Compressed sensing

We consider the problem of robust compressed sensing whose objective is to recover a high-dimensional sparse signal from compressed measurements corrupted by outliers. A new sparse Bayesian learning method is developed for robust compressed sensing. The basic idea of the proposed method is to identify and remove the outliers from sparse signal recovery. To automatically identify the outliers, we employ a set of binary indicator hyperparameters to indicate which observations are outliers. These indicator hyperparameters are treated as random variables and assigned a beta process prior such that their values are confined to be binary. In addition, a Gaussian-inverse Gamma prior is imposed on the sparse signal to promote sparsity. Based on this hierarchical prior model, we develop a variational Bayesian method to estimate the indicator hyperparameters as well as the sparse signal. Simulation results show that the proposed method achieves a substantial performance improvement over existing robust compressed sensing techniques.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Sampling the isothermal-isobaric ensemble by Langevin dynamics

We present a new method of conducting molecular dynamics simulation in isothermal-isobaric ensemble based on Langevin equations of motion. The stochastic coupling to all particle and cell degrees of freedoms is introduced in a correct way, in the sense that the stationary configurational distribution is proved to be in consistent with that of the isothermal-isobaric ensemble. In order to apply the proposed method in computer simulations, a second order symmetric numerical integration scheme is developed by Trotter's splitting of the single-step propagator. Moreover, a practical guide of choosing working parameters is suggested for user specified thermo- and baro-coupling time-scales. The method and software implementation are carefully validated by a numerical example.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

An Efficient Bayesian PAPR Reduction Method for OFDM-Based Massive MIMO Systems

We consider the problem of peak-to-average power ratio (PAPR) reduction in orthogonal frequency-division multiplexing (OFDM) based massive multiple-input multiple-output (MIMO) downlink systems. Specifically, given a set of symbol vectors to be transmitted to K users, the problem is to find an OFDM-modulated signal that has a low PAPR and meanwhile enables multiuser interference (MUI) cancellation. Unlike previous works that tackled the problem using convex optimization, we take a Bayesian approach and develop an efficient PAPR reduction method by exploiting the redundant degrees-of-freedom of the transmit array. The sought-after signal is treated as a random vector with a hierarchical truncated Gaussian mixture prior, which has the potential to encourage a low PAPR signal with most of its samples concentrated on the boundaries. A variational expectation-maximization (EM) strategy is developed to obtain estimates of the hyperparameters associated with the prior model, along with the signal. In addition, the generalized approximate message passing (GAMP) is embedded into the variational EM framework, which results in a significant reduction in computational complexity of the proposed algorithm. Simulation results show our proposed algorithm achieves a substantial performance improvement over existing methods in terms of both the PAPR reduction and computational complexity.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

An Iterative Reweighted Method for Tucker Decomposition of Incomplete Multiway Tensors

We consider the problem of low-rank decomposition of incomplete multiway tensors. Since many real-world data lie on an intrinsically low dimensional subspace, tensor low-rank decomposition with missing entries has applications in many data analysis problems such as recommender systems and image inpainting. In this paper, we focus on Tucker decomposition which represents an Nth-order tensor in terms of N factor matrices and a core tensor via multilinear operations. To exploit the underlying multilinear low-rank structure in high-dimensional datasets, we propose a group-based log-sum penalty functional to place structural sparsity over the core tensor, which leads to a compact representation with smallest core tensor. The method for Tucker decomposition is developed by iteratively minimizing a surrogate function that majorizes the original objective function, which results in an iterative reweighted process. In addition, to reduce the computational complexity, an over-relaxed monotone fast iterative shrinkage-thresholding technique is adapted and embedded in the iterative reweighted process. The proposed method is able to determine the model complexity (i.e. multilinear rank) in an automatic way. Simulation results show that the proposed algorithm offers competitive performance compared with other existing algorithms.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Computationally Efficient Sparse Bayesian Learning via Generalized Approximate Message Passing

The sparse Beyesian learning (also referred to as Bayesian compressed sensing) algorithm is one of the most popular approaches for sparse signal recovery, and has demonstrated superior performance in a series of experiments. Nevertheless, the sparse Bayesian learning algorithm has computational complexity that grows exponentially with the dimension of the signal, which hinders its application to many practical problems even with moderately large data sets. To address this issue, in this paper, we propose a computationally efficient sparse Bayesian learning method via the generalized approximate message passing (GAMP) technique. Specifically, the algorithm is developed within an expectation-maximization (EM) framework, using GAMP to efficiently compute an approximation of the posterior distribution of hidden variables. The hyperparameters associated with the hierarchical Gaussian prior are learned by iteratively maximizing the Q-function which is calculated based on the posterior approximation obtained from the GAMP. Numerical results are provided to illustrate the computational efficacy and the effectiveness of the proposed algorithm.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Global and Local Structure Preserving Sparse Subspace Learning: An Iterative Approach to Unsupervised Feature Selection

As we aim at alleviating the curse of high-dimensionality, subspace learning is becoming more popular. Existing approaches use either information about global or local structure of the data, and few studies simultaneously focus on global and local structures as the both of them contain important information. In this paper, we propose a global and local structure preserving sparse subspace learning (GLoSS) model for unsupervised feature selection. The model can simultaneously realize feature selection and subspace learning. In addition, we develop a greedy algorithm to establish a generic combinatorial model, and an iterative strategy based on an accelerated block coordinate descent is used to solve the GLoSS problem. We also provide whole iterate sequence convergence analysis of the proposed iterative algorithm. Extensive experiments are conducted on real-world datasets to show the superiority of the proposed approach over several state-of-the-art unsupervised feature selection approaches.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

On the Secrecy Capacity of a MIMO Gaussian Wiretap Channel with a Cooperative Jammer

We study the secrecy capacity of a helper-assisted Gaussian wiretap channel with a source, a legitimate receiver, an eavesdropper and an external helper, where each terminal is equipped with multiple antennas. Determining the secrecy capacity in this scenario generally requires solving a nonconvex secrecy rate maximization (SRM) problem. To deal with this issue, we first reformulate the original SRM problem into a sequence of convex subproblems. For the special case of single-antenna legitimate receiver, we obtain the secrecy capacity via a combination of convex optimization and one-dimensional search, while for the general case of multi-antenna legitimate receiver, we propose an iterative solution. To gain more insight into how the secrecy capacity of a helper-assisted Gaussian wiretap channel behaves, we examine the achievable secure degrees of freedom (s.d.o.f.) and obtain the maximal achievable s.d.o.f. in closed-form. We also derive a closed-form solution to the original SRM problem which achieves the maximal s.d.o.f.. Numerical results are presented to illustrate the efficacy of the proposed schemes.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Real-Valued Khatri-Rao Subspace Approaches on the ULA and a New Nested Array

In underdetermined direction-of-arrival (DOA) estimation using the covariance-based signal models, the computational complexity turns into a noticeable issue because of the high dimension of the virtual array manifold. In this paper, real-valued Khatri-Rao (KR) approaches are developed on the uniform linear array (ULA) and the nested array. The complexities of subspace decomposition and spectral search are reduced compared with the complex-valued KR approach. By designing a special transformation matrix, the influence of the noise is removed in the mean time while the data is transformed from the complex domain to the real domain. Deploying the sensors with nonuniform spacings can raise the degree of freedom (DOF) and hence help detect more sources in the underdetermined situation. To increase the DOF further, a new nested array geometry is designed. The real-valued denoising KR approach developed on the new nested array can resolve more sources with reduced complexities. The performance improvement is demonstrated by numerical studies.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Sparse Bayesian Dictionary Learning with a Gaussian Hierarchical Model

We consider a dictionary learning problem whose objective is to design a dictionary such that the signals admits a sparse or an approximate sparse representation over the learned dictionary. Such a problem finds a variety of applications such as image denoising, feature extraction, etc. In this paper, we propose a new hierarchical Bayesian model for dictionary learning, in which a Gaussian-inverse Gamma hierarchical prior is used to promote the sparsity of the representation. Suitable priors are also placed on the dictionary and the noise variance such that they can be reasonably inferred from the data. Based on the hierarchical model, a variational Bayesian method and a Gibbs sampling method are developed for Bayesian inference. The proposed methods have the advantage that they do not require the knowledge of the noise variance \emph{a priori}. Numerical results show that the proposed methods are able to learn the dictionary with an accuracy better than existing methods, particularly for the case where there is a limited number of training signals.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Spectral Compressed Sensing via CANDECOMP/PARAFAC Decomposition of Incomplete Tensors

We consider the line spectral estimation problem which aims to recover a mixture of complex sinusoids from a small number of randomly observed time domain samples. Compressed sensing methods formulates line spectral estimation as a sparse signal recovery problem by discretizing the continuous frequency parameter space into a finite set of grid points. Discretization, however, inevitably incurs errors and leads to deteriorated estimation performance. In this paper, we propose a new method which leverages recent advances in tensor decomposition. Specifically, we organize the observed data into a structured tensor and cast line spectral estimation as a CANDECOMP/PARAFAC (CP) decomposition problem with missing entries. The uniqueness of the CP decomposition allows the frequency components to be super-resolved with infinite precision. Simulation results show that the proposed method provides a competitive estimate accuracy compared with existing state-of-the-art algorithms.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Two-Dimensional Pattern-Coupled Sparse Bayesian Learning via Generalized Approximate Message Passing

We consider the problem of recovering two-dimensional (2-D) block-sparse signals with \emph{unknown} cluster patterns. Two-dimensional block-sparse patterns arise naturally in many practical applications such as foreground detection and inverse synthetic aperture radar imaging. To exploit the block-sparse structure, we introduce a 2-D pattern-coupled hierarchical Gaussian prior model to characterize the statistical pattern dependencies among neighboring coefficients. Unlike the conventional hierarchical Gaussian prior model where each coefficient is associated independently with a unique hyperparameter, the pattern-coupled prior for each coefficient not only involves its own hyperparameter, but also its immediate neighboring hyperparameters. Thus the sparsity patterns of neighboring coefficients are related to each other and the hierarchical model has the potential to encourage 2-D structured-sparse solutions. An expectation-maximization (EM) strategy is employed to obtain the maximum a posterior (MAP) estimate of the hyperparameters, along with the posterior distribution of the sparse signal. In addition, the generalized approximate message passing (GAMP) algorithm is embedded into the EM framework to efficiently compute an approximation of the posterior distribution of hidden variables, which results in a significant reduction in computational complexity. Numerical results are provided to illustrate the effectiveness of the proposed algorithm.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Prior Support Knowledge-Aided Sparse Bayesian Learning with Partly Erroneous Support Information

It has been shown both experimentally and theoretically that sparse signal recovery can be significantly improved given that part of the signal's support is known \emph{a priori}. In practice, however, such prior knowledge is usually inaccurate and contains errors. Using such knowledge may result in severe performance degradation or even recovery failure. In this paper, we study the problem of sparse signal recovery when partial but partly erroneous prior knowledge of the signal's support is available. Based on the conventional sparse Bayesian learning framework, we propose a modified two-layer Gaussian-inverse Gamma hierarchical prior model and, moreover, an improved three-layer hierarchical prior model. The modified two-layer model employs an individual parameter $b_i$ for each sparsity-controlling hyperparameter $α_i$, and has the ability to place non-sparsity-encouraging priors to those coefficients that are believed in the support set. The three-layer hierarchical model is built on the modified two-layer prior model, with a prior placed on the parameters $\{b_i\}$ in the third layer. Such a model enables to automatically learn the true support from partly erroneous information through learning the values of the parameters $\{b_i\}$. Variational Bayesian algorithms are developed based on the proposed hierarchical prior models. Numerical results are provided to illustrate the performance of the proposed algorithms.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Super-Resolution Compressed Sensing: A Generalized Iterative Reweighted L2 Approach

Conventional compressed sensing theory assumes signals have sparse representations in a known, finite dictionary. Nevertheless, in many practical applications such as direction-of-arrival (DOA) estimation and line spectral estimation, the sparsifying dictionary is usually characterized by a set of unknown parameters in a continuous domain. To apply the conventional compressed sensing technique to such applications, the continuous parameter space has to be discretized to a finite set of grid points, based on which a "presumed dictionary" is constructed for sparse signal recovery. Discretization, however, inevitably incurs errors since the true parameters do not necessarily lie on the discretized grid. This error, also referred to as grid mismatch, may lead to deteriorated recovery performance or even recovery failure. To address this issue, in this paper, we propose a generalized iterative reweighted L2 method which jointly estimates the sparse signals and the unknown parameters associated with the true dictionary. The proposed algorithm is developed by iteratively decreasing a surrogate function majorizing a given objective function, resulting in a gradual and interweaved iterative process to refine the unknown parameters and the sparse signal. A simple yet effective scheme is developed for adaptively updating the regularization parameter that controls the tradeoff between the sparsity of the solution and the data fitting error. Extension of the proposed algorithm to the multiple measurement vector scenario is also considered. Numerical results show that the proposed algorithm achieves a super-resolution accuracy and presents superiority over other existing methods.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Super-Resolution Compressed Sensing: An Iterative Reweighted Algorithm for Joint Parameter Learning and Sparse Signal Recovery

In many practical applications such as direction-of-arrival (DOA) estimation and line spectral estimation, the sparsifying dictionary is usually characterized by a set of unknown parameters in a continuous domain. To apply the conventional compressed sensing to such applications, the continuous parameter space has to be discretized to a finite set of grid points. Discretization, however, incurs errors and leads to deteriorated recovery performance. To address this issue, we propose an iterative reweighted method which jointly estimates the unknown parameters and the sparse signals. Specifically, the proposed algorithm is developed by iteratively decreasing a surrogate function majorizing a given objective function, which results in a gradual and interweaved iterative process to refine the unknown parameters and the sparse signal. Numerical results show that the algorithm provides superior performance in resolving closely-spaced frequency components.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

The structure of TeV-bright shell-type supernova remnants

Aims. Two-dimensional MHD simulations are used to model the emission properties of TeV-bright shell-type supernova remnants (SNRs) and to explore their nature. Methods. In the leptonic scenario for the TeV emission, the $γ$-ray emission is produced via Inverse Compton scattering of background soft photons by high-energy electrons accelerated by the shocks of the SNRs. The TeV emissivity is proportional to the magnetic field energy density and MHD simulations can be used to model the TeV structure of such remnants directly. 2D MHD simulations for SNRs are then performed under the assumption that the ambient interstellar medium is turbulent with the magnetic field and density fluctuations following a Kolmogorov-like power-law spectrum. Results. (1) As expected, these simulations confirm early 1D and 2D modelings of these sources, namely the hydrodynamical evolution of the shock waves and amplification of magnetic field by Rayleigh-Taylor convective flows and by shocks propagating in a turbulent medium; (2) We reproduce rather complex morphological structure for $γ$-rays, suggesting intrinsic variations of the source morphology not related to the structure of the progenitor and environment; (3)The observed radial profile of several remnants are well reproduced with an ambient medium density of $0.1-1$ cm$^{-3}$. An even lower ambient density leads to a sharper drop of the TeV brightness with radius than what is observed near the outer edge of these remnants. Conclusions. In a turbulent background medium, we can reproduce the observed characteristics of several shell-type TeV SNRs with reasonable parameters except for a higher ambient density than that inferred from X-ray observations.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

On Feasibility of Generalized Interference Alignment with Partial Interference Cancelation

We study a new IA strategy which is referred to as "Partial Interference Cancelation-based Interference Alignment" (PIC-IA). Unlike the conventional IA strategy, PIC-IA does not strive to eliminate interference from all users. Instead, it aims to remove the most significant interference signals. This PIC-IA strategy generalizes the conventional IA concept by addressing partial, instead of complete, interference cancelation. The feasibility of this new strategy is studied in this paper. Our results show that for a symmetric, single-stream system with $N_t$ transmit antennas and $N_r$ receive antennas, the PIC-IA is feasible when the number of significant interference signals to be removed at each receiver is no more than $N_t+N_r-2$, no matter how many users are in the network. This is in sharp contrast to the conventional IA whose feasibility is severely limited by the number of users $K$.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

One-Bit Quantization Design and Adaptive Methods for Compressed Sensing

There have been a number of studies on sparse signal recovery from one-bit quantized measurements. Nevertheless, little attention has been paid to the choice of the quantization thresholds and its impact on the signal recovery performance. This paper examines the problem of one-bit quantizer design for sparse signal recovery. Our analysis shows that the magnitude ambiguity that ever plagues conventional one-bit compressed sensing methods can be resolved, and an arbitrarily small reconstruction error can be achieved by setting the quantization thresholds close enough to the original data samples without being quantized. Note that unquantized data samples are unaccessible in practice. To overcome this difficulty, we propose an adaptive quantization method that adaptively adjusts the quantization thresholds in a way such that the thresholds converges to the optimal thresholds. Numerical results are illustrated to collaborate our theoretical results and the effectiveness of the proposed algorithm.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Pattern-Coupled Sparse Bayesian Learning for Recovery of Block-Sparse Signals

We consider the problem of recovering block-sparse signals whose structures are unknown \emph{a priori}. Block-sparse signals with nonzero coefficients occurring in clusters arise naturally in many practical scenarios. However, the knowledge of the block structure is usually unavailable in practice. In this paper, we develop a new sparse Bayesian learning method for recovery of block-sparse signals with unknown cluster patterns. Specifically, a pattern-coupled hierarchical Gaussian prior model is introduced to characterize the statistical dependencies among coefficients, in which a set of hyperparameters are employed to control the sparsity of signal coefficients. Unlike the conventional sparse Bayesian learning framework in which each individual hyperparameter is associated independently with each coefficient, in this paper, the prior for each coefficient not only involves its own hyperparameter, but also the hyperparameters of its immediate neighbors. In doing this way, the sparsity patterns of neighboring coefficients are related to each other and the hierarchical model has the potential to encourage structured-sparse solutions. The hyperparameters, along with the sparse signal, are learned by maximizing their posterior probability via an expectation-maximization (EM) algorithm. Numerical results show that the proposed algorithm presents uniform superiority over other existing methods in a series of experiments.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

A Fast Iterative Algorithm for Recovery of Sparse Signals from One-Bit Quantized Measurements

This paper considers the problem of reconstructing sparse or compressible signals from one-bit quantized measurements. We study a new method that uses a log-sum penalty function, also referred to as the Gaussian entropy, for sparse signal recovery. Also, in the proposed method, sigmoid functions are introduced to quantify the consistency between the acquired one-bit quantized data and the reconstructed measurements. A fast iterative algorithm is developed by iteratively minimizing a convex surrogate function that bounds the original objective function, which leads to an iterative reweighted process that alternates between estimating the sparse signal and refining the weights of the surrogate function. Connections between the proposed algorithm and other existing methods are discussed. Numerical results are provided to illustrate the effectiveness of the proposed algorithm.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

A Symmetry-based Decomposition Approach to Eigenvalue Problems: Formulation, Discretization, and Implementation

In this paper, we propose a decomposition approach for eigenvalue problems with spatial symmetries, including the formulation, discretization as well as implementation. This approach can handle eigenvalue problems with either Abelian or non-Abelian symmetries, and is friendly for grid-based discretizations such as finite difference, finite element or finite volume methods. With the formulation, we divide the original eigenvalue problem into a set of subproblems and require only a smaller number of eigenpairs for each subproblem. We implement the decomposition approach with finite elements and parallelize our code in two levels. We show that the decomposition approach can improve the efficiency and scalability of iterative diagonalization. In particular, we apply the approach to solving Kohn--Sham equations of symmetric molecules consisting of hundreds of atoms.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

A Two-Dimensional CA Traffic Model with Dynamic Route Choices Between Residence and Workplace

The Biham, Middleton and Levine (BML) model is extended to describe dynamic route choices between the residence and workplace in cities. The traffic dynamic in the city with a single workplace is studied from the velocity diagram, arrival time probability distribution, destination arrival rate and convergence time. The city with double workplaces is also investigated to compared with a single workplace within the framework of four modes of urban growth. The transitional region is found in the velocity diagrams where the system undergoes a continuous transition from a moving phase to a completely jamming phase. We perform a finite-size scaling analysis of the critical density from a statistical point of view and the order parameter of this jamming transition is estimated. It is also found that statistical properties of urban traffic are greatly influenced by the urban area, workplace area and urban layout.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Jamming Transition of Point-to-Point Traffic Through Cooperative Mechanisms

We study the jamming transition of two-dimensional point-to-point traffic through cooperative mechanisms using computer simulation. We propose two decentralized cooperative mechanisms which are incorporated into the point-to-point traffic models: stepping aside (CM-SA) and choosing alternative routes (CM-CAR). Incorporating CM-SA is to prevent a type of ping-pong jumps from happening when two objects standing face-to-face want to move in opposite directions. Incorporating CM-CAR is to handle the conflict when more than one object competes for the same point in parallel update. We investigate and compare four models mainly from fundamental diagrams, jam patterns and the distribution of cooperation probability. It is found that although it decreases the average velocity a little, the CM-SA increases the critical density and the average flow. Despite increasing the average velocity, the CM-CAR decreases the average flow by creating substantially vacant areas inside jam clusters. We investigate the jam patterns of four models carefully and explain this result qualitatively. In addition, we discuss the advantage and applicability of decentralized cooperation modeling.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Optimal Precoding Design and Power Allocation for Decentralized Detection of Deterministic Signals

We consider a decentralized detection problem in a power-constrained wireless sensor networks (WSNs), in which a number of sensor nodes collaborate to detect the presence of a deterministic vector signal. The signal to be detected is assumed known \emph{a priori}. Given a constraint on the total amount of transmit power, we investigate the optimal linear precoding design for each sensor node. More specifically, in order to achieve the best detection performance, shall sensor nodes transmit their raw data to the fusion center (FC), or transmit compressed versions of their original data? The optimal power allocation among sensors is studied as well. Also, assuming a fixed total transmit power, we examine how the detection performance behaves with the number of sensors in the network. A new concept "detection outage" is proposed to quantify the reliability of the overall detection system. Finally, decentralized detection with unknown signals is studied. Numerical results are conducted to corroborate our theoretical analysis and to illustrate the performance of the proposed algorithm.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

The fundamental Diagram of Pedestrian Model with Slow Reaction

The slow-to-start models are a classical cellular automata model in simulating vehicle traffic. However, to our knowledge, the slow-to-start effect has not considered in modeling pedestrian dynamic. We verify the similar behavior between pedestrian and vehicle, and propose an new lattice gas (LG) model called the slow reaction (SR) model to describe the pedestrian's delayed reaction in single-file movement. We simulate and reproduce the Seyfried's field experiments at the research centre Julich, and use its empirical data to validate our SR model. We compare the SR model with the standard LG model. We test different probability of slow reaction ps in SR model and found the simulation data of ps=0.3 fit the empirical data best. The RMS error of mean velocity of SR model is smaller than that of standard LG model. In the range of ps=0.1~0.3, our fundamental diagram between velocity and density by simulation coincides with field experiments. The distribution of individual velocity in fundamental diagram in SR model agrees with the empirical data better than that of standard LG model. In addition, we observe the stop-and-go waves and phase separation in pedestrian flow by simulation. We reproduced the phenomena of uneven distribution of interspaces by SR model while the standard LG model did not implement. The SR model can reproduce the evolution of spatio-temporal structures of pedestrian flow with higher fidelity to Seyfried's experiments than the standard LG model.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Recovery of Block-Sparse Representations from Noisy Observations via Orthogonal Matching Pursuit

We study the problem of recovering the sparsity pattern of block-sparse signals from noise-corrupted measurements. A simple, efficient recovery method, namely, a block-version of the orthogonal matching pursuit (OMP) method, is considered in this paper and its behavior for recovering the block-sparsity pattern is analyzed. We provide sufficient conditions under which the block-version of the OMP can successfully recover the block-sparse representations in the presence of noise. Our analysis reveals that exploiting block-sparsity can improve the recovery ability and lead to a guaranteed recovery for a higher sparsity level. Numerical results are presented to corroborate our theoretical claim.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Multiband Nonthermal Radiative Properties of HESS J1813-178

The source HESS J1813-178 was detected in the survey of the inner Galaxy in TeV gamma-rays, and a SNR G12.8-0.0 was identified in the radio band to be associated with it. The PWN embedded in the SNR is powered by an energetic pulsar PSR J1813-1749, which was recently discovered. Whether the TeV gamma-rays originate from the SNR shell or the PWN is uncertain now. We investigate theoretically the multiwavelength nonthermal radiation from the composite SNR G12.8-0.0. The emission from the particles accelerated in the SNR shell is calculated based on a semianalytical method to the nonlinear diffusive shock acceleration mechanism. In the model, the magnetic field is self-generated via resonant streaming instability, and the dynamical reaction of the field on the shock is taken into account. Based on a model which couples the dynamical and radiative evolution of a PWN in a non-radiative SNR, the dynamics and the multi-band emission of the PWN are investigated. The particles are injected with a spectrum of a relativistic Maxwellian plus a power law high-energy tail with an index of -2.5. Our results indicate that the radio emission from the shell can be well reproduced as synchrotron radiation of the electrons accelerated by the SNR shock; with an ISM number density of 1.4 cm^{-3} for the remnant, the gamma-ray emission from the SNR shell is insignificant, and the observed X-rays and VHE gamma-rays from the source are consistent with the emission produced by electrons/positrons injected in the PWN via synchrotron radiation and IC scattering, respectively; the resulting gamma-ray flux for the shell is comparable to the detected one only with a relatively larger density of about 2.8 cm^{-3}. The VHE gamma-rays of HESS J1813-178 can be naturally explained to mainly originate from the nebula although the contribution of the SNR shell becomes significant with a denser ambient medium.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

On the Nonthermal Emission from the Supernova Remnant W51C

The middle-aged supernova remnant (SNR) W51C is an interesting source for the interaction of the shell with a molecular cloud. The shell emits intense radio synchrotron photons, and high-energy gamma-rays from the remnant have been detected using the {\it Fermi} Large Area Telescope (LAT), the H.E.S.S. telescope, and the Milagro gamma-ray observatory. Based on a semi-analytical approach to the nonlinear shock acceleration process, we investigate the multiband nonthermal emission from W51C. The result shows that the radio emission from the remnant can be explained as synchrotron radiation of the electrons accelerated by a part of the shock flowing into the ambient medium. On the other hand, the high-energy gamma-rays detected by the {\it Fermi} LAT are mainly produced via proton-proton collisions of the high-energy protons with the ambient matter in the molecular cloud overtaken by the other part of the shock. We propose a possible explanation of the multiband nonthermal emission from W51C, and it can be concluded that a molecular cloud overtaken by a shock wave can be an important emitter in GeV $γ$-rays.

0 citations0 reviewsNo code linkedOpen access