Fast, accurate, and automatic extraction of parameters of nuclear magnetic resonance free induction decay (FID) signal for chemical spectroscopy is a challenging problem. Recently, the Steiglitz-McBride algorithm has been shown to exhibit superior performance in terms of speed, accuracy, and automation when applied to the extraction of T relaxation parameters for myelin water imaging of brain. Applying it to FID data reveals that it falls short of the second objective, the accuracy. Especially, it struggles with the issue of missed spectral peaks when applied to chemical samples with relatively dense frequency spectra. To overcome this issue, a preprocessing stage of subband decomposition is proposed before the application of Steiglitz-McBride algorithm to the FID signal. It is demonstrated that by doing so, a considerable improvement in accuracy is achieved. But this is not gained at the cost of the first objective, the speed. An adaptive subband decomposition is employed in conjunction with the Bayesian information criteria to carry out an efficient decomposition according to spectral content of the signal under investigation. Furthermore, adaptive subband decomposition and the Bayesian information criteria also serve to make the resulting algorithm independent of user input, which also fulfills the third objective, the automation. This makes the proposed algorithm favorable for fast, accurate, and automatic extraction of FID signal parameters.

We report a novel approach for inversion of large random matrices in massive Multiple-Input Multiple Output (MIMO) systems. It is based on the concept of inverse vectors in which an inverse vector is defined for each column of the principal matrix. Such an inverse vector has to satisfy two constraints. Firstly, it has to be in the null-space of all the remaining columns. We call it the null-space problem. Secondly, it has to form a projection of value equal to one in the direction of selected column. We term it as the normalization problem. The process essentially decomposes the inversion problem and distributes it over columns. Each column can be thought of as a node in the network or a particle in a swarm seeking its own solution, the inverse vector, which lightens the computational load on it. Another benefit of this approach is its applicability to all three cases pertaining to a linear system: the fully-determined, the over-determined, and the under-determined case. It eliminates the need of forming the generalized inverse for the last two cases by providing a new way to solve the least squares problem and the Moore and Penrose's pseudoinverse problem. The approach makes no assumption regarding the size, structure or sparsity of the matrix. This makes it fully applicable to much in vogue large random matrices arising in massive MIMO systems. Also, the null-space problem opens the door for a plethora of methods available in literature for null-space computation to enter the realm of matrix inversion. There is even a flexibility of finding an exact or approximate inverse depending on the null-space method employed. We employ the Householder's null-space method for exact solution and present a complete exposition of the new approach. A detailed comparison with well-established matrix inversion methods in literature is also given.

A novel approach for solving linear estimation problem in multi-user massive MIMO systems is proposed. In this approach, the difficulty of matrix inversion is attributed to the incomplete definition of the dot product. The general definition of dot product implies that the columns of channel matrix are always orthogonal whereas, in practice, they may be not. If the latter information can be incorporated into dot product, then the unknowns can be directly computed from projections without inverting the channel matrix. By doing so, the proposed method is able to achieve an exact solution with a 25% reduction in computational complexity as compared to the QR method. Proposed method is stable, offers an extra flexibility of computing any single unknown, and can be implemented in just twelve lines of code.

Rapid, accurate, and automatic quantitation of two‐dimensional nuclear magnetic resonance(2D‐NMR) data is a challenging problem. Recently, a Bayesian information criterion based subband Steiglitz–McBride algorithm has been shown to exhibit superior performance on all three fronts when applied to the quantitation of one‐dimensional NMR free induction decay data. In this paper, we demonstrate that the 2D Steiglitz–McBride algorithm, in conjunction with 2D subband decomposition and the 2D Bayesian information criterion, also achieves excellent results for 2D‐NMR data in terms of speed, accuracy, and automation—especially when compared in these respects to the previously published analysis techniques for 2D‐NMR data.

Impaired tissue perfusion underlies many chronic disease states and aging. Diffusion-weighted imaging (DWI) is a noninvasive MRI technique that has been widely used to characterize tissue perfusion. Parametric models based on DWI measurements can characterize microvascular perfusion modulated by functional and microstructural alterations in the skeletal muscle. The intravoxel incoherent motion (IVIM) model uses a biexponential form to quantify the incoherent motion of water molecules in the microvasculature at low b-values of DWI measurements. The fractional Fickian diffusion (FFD) model is a parsimonious representation of anomalous superdiffusion that uses the stretched exponential form and can be used to quantify the microvascular volume of skeletal muscle. Both models are established measures of perfusion based on DWI, and the prognostic value of model parameters for identifying pathophysiological processes has been studied. Although the mathematical properties of individual models have been previously reported, quantitative connections between IVIM and FFD models have not been examined. This work provides a mathematical framework for obtaining a direct, one-way transformation of the parameters of the stretched exponential model to those of the biexponential model. Numerical simulations are implemented, and the results corroborate analytical results. Additionally, analysis of in vivo DWI measurements in skeletal muscle using both biexponential and stretched exponential models is shown and compared with analytical and numerical models. These results demonstrate the difficulty of model selection based on goodness of fit to experimental data. This analysis provides a framework for better interpreting and harmonizing perfusion parameters from experimental results using these two different models.