On the Generalization of the Random Projection Method for Problems of the Recovery of Object Signal Described by Models of Convolution Type

Abstract

Introduction. In technical systems, there is a common situation when transformation input-output is described by the integral equation of convolution type. This situation accurses if the object signal is recovered by the results of remote measurements. For example, in spectrometric tasks, for an image deblurring, etc. Matrices of the discrete representation for the output signal and the kernel of convolution are known. We need to find a matrix of the discrete representation of a signal of the object. The well known approach for solving this problem includes the next steps. First, the kernel matrix has to be represented as the Kroneker product. Second, the input-output transformation has to be presented with the usage of Kroneker product matrices. Third, the matrix of the discrete representation of the object has to be found. The object signal matrix estimation obtained with the help of pseudo inverting of Kroneker decomposition matrices is unstable. The instability of the object signal estimation in the case of usage of Kroneker decomposition matrices is caused by their discrete ill posed matrix properties (condition number is big and the series of the singular numbers smoothly decrease to zero). To find solutions of discrete ill-posed problems we developed methods based on the random projection and the random projection with an averaging by the random matrices. These methods provide a stable solutions with a small computational complexity. We consider the problem of object signals recovering in the systems where an input-output transformation is described by the integral equation of a convolution. To find a solution for these problems we need to build a generalization for two-dimensional signals case of the random projection method. Purpose. To develop a stable method of the recovery of object signal for the case in which an input-output transformation is described by the integral equation of a convolution. Results and conclusions. We developed the method of a stable recovery of object signal for the case in which an input-output transformation is described by the integral equation of a convolution. The stable estimation of the object signal is provided by Kroneker decomposition of the kernel matrix of convolution, computation of random projections for Kroneker factorization matrices, and a selection of the optimal dimension of a projector matrix. The method is illustrated by its application in technical problems.