English

System Identification in Dynamical Sampling

Information Theory 2017-06-19 v2 math.IT

Abstract

We consider the problem of spatiotemporal sampling in a discrete infinite dimensional spatially invariant evolutionary process x(n)=Anxx^{(n)}=A^nx to recover an unknown convolution operator AA given by a filter a1(Z)a \in \ell^1(\mathbb{Z}) and an unknown initial state xx modeled as avector in 2(Z)\ell^2(\mathbb{Z}). Traditionally, under appropriate hypotheses, any xx can be recovered from its samples on Z\mathbb{Z} and AA 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 AA and xx that allows to sample the evolving states x,Ax,,AN1xx,Ax, \cdots, A^{N-1}x on a sub-lattice of Z\mathbb{Z}, and thus achieve the spatiotemporal trade off. The spatiotemporal trade off is motivated by several industrial applications \cite{Lv09}. Specifically, we show that {x(mZ),Ax(mZ),,AN1x(mZ):N2m}\{x(m\mathbb{Z}), Ax(m\mathbb{Z}), \cdots, A^{N-1}x(m\mathbb{Z}): N \geq 2m\} contains enough information to recover a typical "low pass filter" aa and xx almost surely, in which we generalize the idea of the finite dimensional case in \cite{AK14}. In particular, we provide an algorithm based on a generalized Prony method for the case when both aa and xx 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 AA and initial state xx, and verify them by several numerical experiments. Finally, we provide several other numerical methods to stabilize the method and numerical example shows the improvement.

Keywords

Cite

@article{arxiv.1502.02741,
  title  = {System Identification in Dynamical Sampling},
  author = {Sui Tang},
  journal= {arXiv preprint arXiv:1502.02741},
  year   = {2017}
}

Comments

25 pages, 7 figures

R2 v1 2026-06-22T08:26:07.185Z