System Identification in Dynamical Sampling
Abstract
We consider the problem of spatiotemporal sampling in a discrete infinite dimensional spatially invariant evolutionary process x(n)=Anx to recover an unknown convolution operator A given by a filter a ∈ 1(Z) and an unknown initial state x modeled as avector in 2(Z). Traditionally, under appropriate hypotheses, any x can be recovered from its samples on Z and A can be recovered by the classical techniques of deconvolution. In this paper, we will exploit the spatiotemporal correlation and propose a new spatiotemporal sampling scheme to recover A and x that allows to sample the evolving states x,Ax, ·s, AN-1x on a sub-lattice of Z, and thus achieve the spatiotemporal trade off. The spatiotemporal trade off is motivated by several industrial applications Lv09. Specifically, we show that \x(mZ), Ax(mZ), ·s, AN-1x(mZ): N ≥ 2m\ contains enough information to recover a typical "low pass filter" a and x almost surely, in which we generalize the idea of the finite dimensional case in AK14. In particular, we provide an algorithm based on a generalized Prony method for the case when both a and x are of finite impulse response and an upper bound of their support is known. We also perform the perturbation analysis based on the spectral properties of the operator A and initial state x, and verify them by several numerical experiments. Finally, we provide several other numerical methods to stabilize the method and numerical example shows the improvement.