Missing data are a problem in geophysical surveys, and interpolation and reconstruction of missing data is part of the data processing and interpretation. Based on the sparseness of the geophysical data or the transfo...Missing data are a problem in geophysical surveys, and interpolation and reconstruction of missing data is part of the data processing and interpretation. Based on the sparseness of the geophysical data or the transform domain, we can improve the accuracy and stability of the reconstruction by transforming it to a sparse optimization problem. In this paper, we propose a mathematical model for the sparse reconstruction of data based on the LO-norm minimization. Furthermore, we discuss two types of the approximation algorithm for the LO- norm minimization according to the size and characteristics of the geophysical data： namely, the iteratively reweighted least-squares algorithm and the fast iterative hard thresholding algorithm. Theoretical and numerical analysis showed that applying the iteratively reweighted least-squares algorithm to the reconstruction of potential field data exploits its fast convergence rate, short calculation time, and high precision, whereas the fast iterative hard thresholding algorithm is more suitable for processing seismic data, moreover, its computational efficiency is better than that of the traditional iterative hard thresholding algorithm.展开更多
Using the inversion of the auto correlation function Toeplitz matrix of pseudo random binary sequence (PRBS) derived in this paper and the theorem of partitioned matrix inversion, a fast multistage least squares (FM...Using the inversion of the auto correlation function Toeplitz matrix of pseudo random binary sequence (PRBS) derived in this paper and the theorem of partitioned matrix inversion, a fast multistage least squares (FMLS) method is developed. Its performances are theoretically analyzed and digital simulation is made to compare FMLS with multistage least squares (MSLS), correlation least squares(COR LS) and LS for their computer speed and identification accuracy. Finally, FMLS is applied to identifying the heat excharger dynamics. It is shown that FMLS is a good and effective identification technique.展开更多
基金supported by the National Natural Science Foundation of China (Grant No.41074133)
文摘Missing data are a problem in geophysical surveys, and interpolation and reconstruction of missing data is part of the data processing and interpretation. Based on the sparseness of the geophysical data or the transform domain, we can improve the accuracy and stability of the reconstruction by transforming it to a sparse optimization problem. In this paper, we propose a mathematical model for the sparse reconstruction of data based on the LO-norm minimization. Furthermore, we discuss two types of the approximation algorithm for the LO- norm minimization according to the size and characteristics of the geophysical data： namely, the iteratively reweighted least-squares algorithm and the fast iterative hard thresholding algorithm. Theoretical and numerical analysis showed that applying the iteratively reweighted least-squares algorithm to the reconstruction of potential field data exploits its fast convergence rate, short calculation time, and high precision, whereas the fast iterative hard thresholding algorithm is more suitable for processing seismic data, moreover, its computational efficiency is better than that of the traditional iterative hard thresholding algorithm.
文摘Using the inversion of the auto correlation function Toeplitz matrix of pseudo random binary sequence (PRBS) derived in this paper and the theorem of partitioned matrix inversion, a fast multistage least squares (FMLS) method is developed. Its performances are theoretically analyzed and digital simulation is made to compare FMLS with multistage least squares (MSLS), correlation least squares(COR LS) and LS for their computer speed and identification accuracy. Finally, FMLS is applied to identifying the heat excharger dynamics. It is shown that FMLS is a good and effective identification technique.
