Subspace averaging and order determination for source enumeration

This paper was presented in part at the 2016 Workshop on Statistical Signal Processing, Palma de Mallorca, Spain, June 2016, and in part at the 2018 Workshop on Statistical Signal Processing, Freiburg, Germany, June 2018.In this paper, we address the problem of subspace averaging, with special emphasis placed on the question of estimating the dimension of the average. The results suggest that the enumeration of sources in a multi-sensor array, which is a problem of estimating the dimension of the array manifold, and as a consequence the number of radiating sources, may be cast as a problem of averaging subspaces. This point of view stands in contrast to conventional approaches, which cast the problem as one of identifiying covariance models in a factor model. We present a robust formulation of the proposed order fitting rule based on majorization-minimization algorithms. A key element of the proposed method is to construct a bootstrap procedure, based on a newly proposed discrete distribution on the manifold of projection matrices, for stochastically generating subspaces from a function of experimentally determined eigenvalues. In this way, the proposed subspace averaging (SA) technique determines the order based on the eigenvalues of an average projection matrix, rather than on the likelihood of a covariance model, penalized by functions of the model order. By means of simulation examples, we show that the proposed SA criterion is especially effective in high-dimensional scenarios with low sample support.The work of V. Garg and I. Santamaría was supported in part by the Ministerio de Economía y Competitividad (MINECO) of Spain, and in part by the AEI/FEDER funds of the E.U., under Grants TEC2016-75067-C4-4-R (CARMEN), TEC2015-69648-REDC, and BES-2017-080542. The work of D. Ramírez was supported in part by the Ministerio de Ciencia, Innovación y Universidades under Grant TEC2017-92552-EXP (aMBITION), in part by the Ministerio de Ciencia, Innovación y Universidades, jointly with the European Commission (ERDF), under Grants TEC2015-69868- C2-1-R (ADVENTURE) and TEC2017-86921-C2-2-R (CAIMAN), and in part by The Comunidad de Madrid under Grant Y2018/TCS-4705 (PRACTICOCM). The work of L.L Scharf was supported in part by the U.S. NSF under Contract CISE-1712788.Publicad