It is with great pleasure that we introduce this special issue of Signal Processing, published to mark the 45th anniversary of the journal. Since its inception in 1979, Signal Processing has played a foundational role in shaping the evolution of the field. As the flagship journal of the European Association for Signal Processing (EURASIP), it has served not only as a premier venue for scientific dissemination but also as a pillar of the signal processing research community. [...]

In the paper, we consider the line spectral estimation problem in an unlimited sensing framework (USF), where a modulo analog-to-digital converter (ADC) is employed to fold the input signal back into a bounded interval before quantization. Such an operation is mathematically equivalent to taking the modulo of the input signal with respect to the interval. To overcome the noise sensitivity of higher-order difference-based methods, we explore the properties of the first-order difference of modulo samples, and develop two line spectral estimation algorithms based on the first-order difference, which are robust against noise. Specifically, we show that, with a high probability, the first-order difference of the original samples is equivalent to that of the modulo samples. By utilizing this property, line spectral estimation is solved via a robust sparse signal recovery approach. The second algorithms is built on our finding that, with a sufficiently high sampling rate, the first-order difference of the original samples can be decomposed as a sum of the first-order difference of the modulo samples and a sequence whose elements are confined to three possible values. This decomposition enables us to formulate the line spectral estimation problem as a mixed integer linear program that can be efficiently solved. Simulation results show that both proposed methods are robust against noise and achieve a significant performance improvement over the higher-order difference-based method. methods.

This paper investigates the positioning of the pilot symbols, as well as the power distribution between the pilot and the communication symbols for the orthogonal time frequency space (OTFS) modulation scheme. We analyze the pilot placements that minimize the mean squared error (MSE) in estimating the channel taps. This allows us to identify two new pilot allocations for OTFS that save approximately 50% of the pilot overhead compared to existing allocations. In addition, we optimize the average channel capacity by adjusting the power distribution. We show that this leads to a significant increase in average capacity. The results provide valuable guidance for designing the OTFS parameters to achieve maximum capacity. Numerical simulations are performed to validate the findings.

Developing methods to process irregularly structured data is crucial in applications like gene-regulatory, brain, power, and socioeconomic networks. Graphs have been the go-to algebraic tool for modeling the structure via nodes and edges capturing their interactions, leading to the establishment of the fields of graph signal processing (GSP) and graph machine learning (GML). Key graph-aware methods include Fourier transform, filtering, sampling, as well as topology identification and spatiotemporal processing. Although versatile, graphs can model only pairwise dependencies in the data. To this end, topological structures such as simplicial and cell complexes have emerged as algebraic representations for more intricate structure modeling in data-driven systems, fueling the rapid development of novel topological-based processing and learning methods. This paper first presents the core principles of topological signal processing through the Hodge theory, a framework instrumental in propelling the field forward thanks to principled connections with GSP-GML. It then outlines advances in topological signal representation, filtering, and sampling, as well as inferring topological structures from data, processing spatiotemporal topological signals, and connections with topological machine learning. The impact of topological signal processing and learning is finally highlighted in applications dealing with flow data over networks, geometric processing, statistical ranking, biology, and semantic communication.

This paper addresses graph topology identification for applications where the underlying structure of systems like brain and social networks is not directly observable. Traditional approaches based on signal matching and spectral templates have limitations, particularly in handling scale issues and sparsity assumptions. We introduce a novel covariance matching methodology that efficiently reconstructs the graph topology using observable data. For the structural equation model (SEM) using an undirected graph, we demonstrate that our method can converge to the correct result under relatively soft conditions. Furthermore, we extend our methodology to polynomial models and any known distribution of latent variables, broadening its applicability and utility in diverse graph-based systems.

CANDECOMP/PARAFAC (CP) decomposition is the mostly used model to formulate the received tensor signal in a massive MIMO system, as the receiver generally sums the components from different paths or users. To achieve accurate and low-latency channel estimation, good and fast CP decomposition (CPD) algorithms are desired. The CP alternating least squares (CPALS) is the workhorse algorithm for calculating the CPD. However, its performance depends on the initializations, and good starting values can lead to more efficient solutions. Existing initialization strategies are decoupled from the CPALS and are not necessarily favorable for solving the CPD. This paper proposes a deep-learning-aided CPALS (DL-CPALS) method that uses a deep neural network (DNN) to generate favorable initializations. The proposed DL-CPALS integrates the DNN and CPALS to a model-based deep learning paradigm, where it trains the DNN to generate an initialization that facilitates fast and accurate CPD. Moreover, benefiting from the CP low-rankness, the proposed method is trained using noisy data and does not require paired clean data. The proposed DL-CPALS is applied to millimeter wave MIMO-OFDM channel estimation. Experimental results demonstrate the significant improvements of the proposed method in terms of both speed and accuracy for CPD and channel estimation.

Sparse array design is used to help reduce computational, hardware, and power requirements compared to uniform arrays while maintaining acceptable performance. Although minimizing the Cramér-Rao bound has been adopted previously for sparse sensing, it did not consider multiple targets and unknown target directions. To handle the unknown target directions when optimizing the Cramér-Rao bound, we propose to use the worst-case Cramér-Rao bound of two uncorrelated equal power sources with arbitrary angles. This new worst-case two-target Cramér-Rao bound metric has some resemblance to the peak sidelobe level metric which is commonly used in unknown multi-target scenarios. We cast the sensor selection problem for 3-D arrays using the worst-case two-target Cramér-Rao bound as a convex semi-definite program and obtain the binary selection by randomized rounding. We illustrate the proposed method through numerical examples, comparing it to solutions obtained by minimizing the single-target Cramér-Rao bound, minimizing the Cramér-Rao bound for known target angles, the concentric rectangular array and the boundary array. We show that our method selects a combination of edge and center elements, which contrasts with solutions obtained by minimizing the single-target Cramér-Rao bound. The proposed selections also exhibit lower peak sidelobe levels without the need for sidelobe level constraints.

Ultrasonography could allow operator-independent examination and continuous monitoring of the carotid artery (CA) but normally requires complex and expensive transducers, especially for 3-D. By employing computational ultrasound imaging (cUSi), using an aberration mask and model-based reconstruction, a monitoring device could be constructed with a more affordable simple transducer design comprising only a few elements. We aim to apply the cUSi concept to create a CA monitoring system. The system’s possible configurations for the 2-D imaging case were explored using a linear array setup emulating a cUSi device in silico, followed by in vitro testing and in vivo CA imaging. Our study shows enhanced reconstruction performance with the use of an aberrating mask, improved lateral resolution through proper choice of the mask delay variation, and more accurate reconstructions using least-squares with QR (LSQR) decomposition compared to matched filtering (MF). Together, these advancements enable B-mode reconstruction and power Doppler imaging (PDI) of the CA with sufficient quality for monitoring using a configuration of 12 transceivers coupled with a random aberration mask with a maximum delay variation of four wave periods (WPs).

Neural networks on simplicial complexes (SCs) can learn representations from data residing on simplices such as nodes, edges, triangles, etc. However, existing works often overlook the Hodge theorem that decomposes simplicial data into three orthogonal characteristic subspaces, such as the identifiable gradient, curl and harmonic components of edge flows. This provides a universal tool to understand the machine learning models on SCs, thus, allowing for better principled and effective learning. In this paper, we study the effect of this data inductive bias on learning on SCs via the principle of convolutions. Particularly, we present a general convolutional architecture that respects the three key principles of uncoupling the lower and upper simplicial adjacencies, accounting for the inter-simplicial couplings, and performing higher-order convolutions. To understand these principles, we first use Dirichlet energy minimizations on SCs to interpret their effects on mitigating simplicial oversmoothing. Then, we show the three principles promote the Hodge-aware learning of this architecture, through the lens of spectral simplicial theory, in the sense that the three Hodge subspaces are invariant under its learnable functions and the learning in two nontrivial subspaces is independent and expressive. Third, we investigate the learning ability of this architecture in optic of perturbation theory on simplicial topologies and prove that the convolutional architecture is stable to small perturbations. Finally, we corroborate the three principles by comparing with methods that either violate or do not respect them. Overall, this paper bridges learning on SCs with the Hodge theorem, highlighting its importance for rational and effective learning from simplicial data, and provides theoretical insights to convolutional learning on SCs.

This paper investigates jointly optimal array geometry and waveform designs for active sensing. Specifically, we focus on minimizing the Cramér-Rao lower bound (CRB) of the angle of a single target in white Gaussian noise. We first find that several array-waveform pairs can yield the same CRB by virtue of sequences with equal sums of squares, i.e., solutions to certain Diophantine equations. Furthermore, we show that under physical aperture and sensor number constraints, the CRB-minimizing receive array geometry is unique, whereas the transmit array can be chosen flexibly. We leverage this freedom to design a novel sparse array geometry that not only minimizes the single-target CRB given an optimal waveform, but also has a nonredundant and contiguous sum co-array—a desirable property when launching independent waveforms, with relevance also to the multi-target case.

In this letter, we propose an analytical solution for recovering a low-rank positive semi-definite (PSD) matrix from its rank-one measurements. We show that by utilizing a set of structured measurement vectors, we can analytically determine the null space of this low-rank PSD matrix. Based on the result, the PSD matrix can be efficiently recovered. Our analysis shows that the proposed method only requires (N - K)(2K + 1)+ K ² measurements to guarantee exact recovery of the PSD matrix, where N and K respectively denote the dimension and the rank of the PSD matrix. Numerical results show that the proposed method achieves a considerable improvement over existing state-of-the-art methods in terms of both sample complexity and computational efficiency. Specifically, the proposed method helps improve the computational efficiency by an order of magnitude as compared with existing methods.

Four-dimensional ultrasound imaging of complex biological systems such as the brain is technically challenging because of the spatiotemporal sampling requirements. We present computational ultrasound imaging (cUSi), an imaging method that uses complex ultrasound fields that can be generated with simple hardware and a physical wave prediction model to alleviate the sampling constraints. cUSi allows for high-resolution four-dimensional imaging of brain hemodynamics in awake and anesthetized mice.

In this paper, we present a novel convolution theorem which encompasses the well known convolution theorem in (graph) signal processing as well as the one related to time-varying filters. Specifically, we show how a node-wise convolution for signals supported on a graph can be expressed as another node-wise convolution in a frequency domain graph, different from the original graph. This is achieved through a parameterization of the filter coefficients following a basis expansion model. After showing how the presented theorem is consistent with the already existing body of literature, we discuss its implications in terms of non-stationarity. Finally, we propose a data-driven algorithm based on subspace fitting to learn the frequency domain graph, which is then corroborated by experimental results on synthetic and real data.

Identifying overlapping communities from data is crucial for grasping the complex structure and dynamics of networks, amongst others in fields such as computational neuroscience. Research using fMRI has demonstrated that brain regions can change their functional network membership over time using temporal independent component analysis (tICA). However, reproducibility of such overlapping communities remains a challenge. Recently, several alternative approaches have been proposed to identify such overlapping communities. While results are promising, less is known about the model and assumptions that underlie these approaches. This paper shows that the bilinear model, combined with the assumption of quasi-stationary and uncorrelated sources, underlies novel methods for identifying overlapping brain networks. Furthermore, we propose a new algorithm, and through simulations, we investigate the robustness of our algorithm and several existing methods to solve the problem in noisy conditions with few available data samples. We conclude that quasi-stationary blind source separation-based techniques can have a promising advantage over tICA in terms of identifiability of overlapping communities and thus have the potential to contribute towards greater reproducibility of results.

In this work, we deal with the problem of reconstructing a complete bandlimited graph signal from partially sampled noisy measurements. For a known graph structure, an efficient greedy algorithm is presented to partition the graph nodes into disjoint subsets such that sampling the graph signal from any subset leads to a sufficiently accurate reconstruction. Furthermore, we consider a scenario where the graph is massive and data processing centrally is no longer practical. To overcome this issue, a distributed framework is proposed that allows us to implement partitioning algorithms in a parallelized fashion. Finally, we provide numerical simulation results on synthetic and real-world data to show that our proposals outperform the state-of-the-art.

Hearing impairment is a prevalent problem with daily challenges like impaired speech intelligibility and sound localisation. One of the shortcomings of spatial filtering in hearing aids is that speech intelligibility is often not optimised directly, meaning that different auditory processes contributing to intelligibility are often not considered. One example is the perceptual phenomenon known as spatial release from masking (SRM). This paper develops a signal model that explicitly considers SRM in the beamforming design, achieved by transforming the binaural intelligibility prediction model (BSIM) into a signal processing framework. The resulting extended signal model is used to analyse the performance of reference beamformers and design a novel beamformer that more closely considers how the auditory system perceives binaural sound. It can be shown that the binaural minimum variance distortionless response (BMVDR) beamformer is also an optimal solution for the extended, perceived model, suggesting that SRM does not play a significant role in intelligibility enhancement after optimal beamforming. However, the optimal beamformer is no longer unique in the extended signal model. The additional secondary degrees of freedom can be used to preserve binaural cues of interfering sources while still achieving the same perceived performance of the BMVDR beamformer, though with a possible high sensitivity to intelligibility model mismatch errors.

We consider the problem of recovering complex-valued block sparse signals with unknown borders. Such signals arise naturally in numerous applications. Several algorithms have been developed to solve the problem of unknown block partitions. In pattern-coupled sparse Bayesian learning (PCSBL), each coefficient involves its own hyperparameter and those of its immediate neighbors to exploit the block sparsity. Extended block sparse Bayesian learning (EBSBL) assumes the block sparse signal consists of correlated and overlapping blocks to enforce block correlations. We propose a simpler alternative to EBSBL and reveal the underlying relationship between the proposed method and a particular case of EBSBL. The proposed algorithm uses the fact that immediate neighboring sparse coefficients are correlated. The proposed model is similar to classical sparse Bayesian learning (SBL). However, unlike the diagonal correlation matrix in conventional SBL, the unknown correlation matrix has a tridiagonal structure to capture the correlation with neighbors. Due to the entanglement of the elements in the inverse tridiagonal matrix, instead of a direct closed-form solution, an approximate solution is proposed. The alternative algorithm avoids the high dictionary coherence in EBSBL, reduces the unknowns of EBSBL, and is computationally more efficient. The sparse reconstruction performance of the algorithm is evaluated with both correlated and uncorrelated block sparse coefficients. Simulation results demonstrate that the proposed algorithm outperforms PCSBL and correlation-based methods such as EBSBL in terms of reconstruction quality. The numerical results also show that the proposed correlated SBL algorithm can deal with isolated zeros and nonzeros as well as block sparse patterns.

The main focus of this paper is an active sensing application that involves selecting transmit and receive sensors to optimize the Cramér-Rao bound (CRB) on target parameters. Although the CRB is non-convex in the transmit and receive selection, we demonstrate that it is convex in the virtual array weight vector, which describes the multiplicity of the virtual array elements. Based on this finding, we propose a novel algorithm that optimizes the virtual array weight vector first and then finds a matching transceiver array. This greatly enhances the efficiency of the transmit and receive sensor selection problem.

Computational ultrasound imaging (cUSi) offers high-resolution 3D imaging with simpler hardware by relying on computational power. Central to cUSi is a large model matrix that stores all pulse-echo signals. For 3D imaging this matrix easily surpasses 1 terabyte, hindering in-memory storage and real-time processing. This paper presents a solution for cUSi through an aberrating layer by introducing a virtual array concept, which uses transfer functions to map data from the real to a virtual array, enabling the use of conventional reconstruction techniques like delay-and-sum (DAS). We demonstrate the mathematical similarity of this approach to using a full model matrix and validate it with promising imaging results.

In this paper, we propose a new method for joint ranging and Phase Offset (PO) estimation of multiple transponder-equipped aviation vehicles (TEAVs), including Manned Aerial Vehicles (MAVs) and Unmanned Aerial Vehicles (UAVs). The proposed method employs the overlapping uncoordinated Automatic Dependent Surveillance-Broadcast (ADS-B) packets broadcasted by the TEAVs for joint range and PO estimation prior to ADS-B packet decoding; thus, it can improve air safety when packet decoding is infeasible due to packet collision. Moreover, it enables coherent detection of ADS-B packets, which can result in more reliable multiple target tracking in aviation systems using cooperative sensors for sense and avoid. By minimizing the Kullback-Leibler Divergence (KLD), we show that the received complex baseband signal, coming from K uncoordinated TEAVs, which is corrupted by Additive White Gaussian Noise (AWGN) at a single antenna receiver can be approximated by an independent and identically distributed (i.i.d.) Gaussian Mixture (GM) with 2^K mixture components in the two-dimensional plane. The proposed estimator employs the Expectation-Maximization (EM) algorithm to estimate the modes of the 2D Gaussian mixture followed by a reordering estimation technique to jointly estimate range and PO. Simulation results show that the proposed joint estimator outperforms excising methods, such as the time segmentation method and the blind adaptive beamforming.