Source author record

Mattia Zorzi

Mattia Zorzi appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

33works
9topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

33 published item(s)

preprint2023arXiv

A Well-Posed Multidimensional Rational Covariance and Generalized Cepstral Extension Problem

In the present paper we consider the problem of estimating the multidimensional power spectral density which describes a second-order stationary random field from a finite number of covariance and generalized cepstral coefficients. The latter can be framed as an optimization problem subject to multidimensional moment constraints, i.e., to search a spectral density maximizing an entropic index and matching the moments. In connection with systems and control, such a problem can also be posed as finding a multidimensional shaping filter (i.e., a linear time-invariant system) which can output a random field that has identical moments with the given data when fed with a white noise, a fundamental problem in system identification. In particular, we consider the case where the dimension of the random field is greater than two for which a satisfying theory is still missing. We propose a multidimensional moment problem which takes into account a generalized definition of the cepstral moments, together with a consistent definition of the entropy. We show that it is always possible to find a rational power spectral density matching exactly the covariances and approximately the generalized cepstral coefficients, from which a shaping filter can be constructed via spectral factorization. In plain words, our theory allows to construct a well-posed spectral estimator for any finite dimension.

preprint2023arXiv

On the Statistical Consistency of a Generalized Cepstral Estimator

We consider the problem to estimate the generalized cepstral coefficients of a stationary stochastic process or stationary multidimensional random field. It turns out that a naive version of the periodogram-based estimator for the generalized cepstral coefficients is not consistent. We propose a consistent estimator for those coefficients. Moreover, we show that the latter can be used in order to build a consistent estimator for a particular class of cascade linear stochastic systems.

preprint2022arXiv

A second-order generalization of TC and DC kernels

Kernel-based methods have been successfully introduced in system identification to estimate the impulse response of a linear system. Adopting the Bayesian viewpoint, the impulse response is modeled as a zero mean Gaussian process whose covariance function (kernel) is estimated from the data. The most popular kernels used in system identification are the tuned-correlated (TC), the diagonal-correlated (DC) and the stable spline (SS) kernel. TC and DC kernels admit a closed form factorization of the inverse. The SS kernel induces more smoothness than TC and DC on the estimated impulse response, however, the aforementioned property does not hold in this case. In this paper we propose a second-order extension of the TC and DC kernel which induces more smoothness than TC and DC, respectively, on the impulse response and a generalized-correlated kernel which incorporates the TC and DC kernels and their second order extensions. Moreover, these generalizations admit a closed form factorization of the inverse and thus they allow to design efficient algorithms for the search of the optimal kernel hyperparameters. We also show how to use this idea to develop higher oder extensions. Interestingly, these new kernels belong to the family of the so called exponentially convex local stationary kernels: such a property allows to immediately analyze the frequency properties induced on the estimated impulse response by these kernels.

preprint2022arXiv

Distributed Kalman filtering with event-triggered communication: a robust approach

We consider the problem of distributed Kalman filtering for sensor networks in the case there is a limit in data transmission and there is model uncertainty. More precisely, we propose a distributed filtering strategy with event-triggered communication in which the state estimators are computed according to the least favorable model. The latter belongs to a ball (in Kullback-Leibler topology) about the nominal model. We also present a preliminary numerical example in order to test the performance of the proposed strategy.

preprint2022arXiv

Mean-square consistency of the $f$-truncated $\text{M}^2$-periodogram

The paper deals with the problem of estimating the M$^2$ (i.e. multivariate and multidimensional) spectral density function of a stationary random process or random field. We propose the $f$-truncated periodogram, i.e. a truncated periodogram where the truncation point is a suitable function $f$ of the sample size. We discuss the asymptotic consistency of the estimator and we provide three concrete problems that can be solved using the proposed approach. Simulation results show the effectiveness of the procedure.

preprint2022arXiv

Robust Distributed Kalman filtering with Event-Triggered Communication

We consider the problem of distributed Kalman filtering for sensor networks in the case there are constraints in data transmission and there is model uncertainty. More precisely, we propose two distributed filtering strategies with event-triggered communication where the state estimators are computed according to the least favorable model. The latter belongs to a ball about the nominal model. We also show that both the methods are stable in the sense that the mean-square of the state estimation error is bounded in all the nodes.

preprint2021arXiv

Learning the tuned liquid damper dynamics by means of a robust EKF

The tuned liquid dampers (TLD) technology is a feasible and cost-effective seismic design. In order to improve its efficiency it is fundamental to find accurate models describing their dynamic. A TLD system can be modeled through the Housner model and its parameters can be estimated by solving a nonlinear state estimation problem. We propose a robust extended Kalman filter which alleviates the model discretization and the fact that the noise process is not known. We test the effectiveness of the proposed approach by using some experimental data corresponding to two classical seismic waves, namely the El Centro wave and the Hachinohe wave.

preprint2021arXiv

Optimal Transport between Gaussian Stationary Processes

We consider the optimal transport problem between multivariate Gaussian stationary stochastic processes. The transportation effort is the variance of the filtered discrepancy process. The main contribution of this technical note is to show that the corresponding solution leads to a weighted Hellinger distance between multivariate power spectral densities. Then, we propose a spectral estimation approach in the case of indirect measurements which is based on this distance.

preprint2020arXiv

A new kernel-based approach for spectral estimation

The paper addresses the problem to estimate the power spectral density of an ARMA zero mean Gaussian process. We propose a kernel based maximum entropy spectral estimator. The latter searches the optimal spectrum over a class of high order autoregressive models while the penalty term induced by the kernel matrix promotes regularity and exponential decay to zero of the impulse response of the corresponding one-step ahead predictor. Moreover, the proposed method also provides a minimum phase spectral factor of the process. Numerical experiments showed the effectiveness of the proposed method.

preprint2020arXiv

Autoregressive Identification of Kronecker Graphical Models

We address the problem to estimate a Kronecker graphical model corresponding to an autoregressive Gaussian stochastic process. The latter is completely described by the power spectral density function whose inverse has support which admits a Kronecker product decomposition. We propose a Bayesian approach to estimate such a model. We test the effectiveness of the proposed method by some numerical experiments. We also apply the procedure to urban pollution monitoring data.

preprint2020arXiv

Data-driven Link Prediction over Graphical Models

The positive link prediction (PLP) problem is formulated in a system identification framework: we consider dynamic graphical models for auto-regressive moving-average (ARMA) Gaussian random processes. For the identification of the parameters, we model our network on two different time scales: a quicker one, over which we assume that the process representing the dynamics of the agents can be considered to be stationary, and a slower one in which the model parameters may vary. The latter accounts for the possible appearance of new edges. The identification problem is cast into an optimization framework which can be seen as a generalization of the existing methods for the identification of ARMA graphical models. We prove the existence and uniqueness of the solution of such an optimization problem and we propose a procedure to compute numerically this solution. Simulations testing the performances of our method are provided.

preprint2020arXiv

Distributed Kalman Filtering under Model Uncertainty

We study the problem of distributed Kalman filtering for sensor networks in the presence of model uncertainty. More precisely, we assume that the actual state-space model belongs to a ball, in the Kullback-Leibler topology, about the nominal state-space model and whose radius reflects the mismatch modeling budget allowed for each time step. We propose a distributed Kalman filter with diffusion step which is robust with respect to the aforementioned model uncertainty. Moreover, we derive the corresponding least favorable performance. Finally, we check the effectiveness of the proposed algorithm in the presence of uncertainty through a numerical example.

preprint2020arXiv

Image compression by means of the multidimensional circulant covariance extension problem -- Revisited

We revisit the image compression problem using the framework introduced by Ringh, Karlsson and Lindquist. More precisely, we explore the possibility to consider a family of objective functions and a different way to design the prior in the corresponding multidimensional circulant covariance extension problem. The latter leads to refined compression paradigms.

preprint2020arXiv

Link Prediction: A Graphical Model Approach

We consider the problem of link prediction in networks whose edge structure may vary (sufficiently slowly) over time. This problem, with applications in many important areas including social networks, has two main variants: the first, known as positive link prediction or PLP consists in estimating the appearance of a link in the network. The second, known as negative link prediction or NLP consists in estimating the disappearance of a link in the network. We propose a data-driven approach to estimate the appearance/disappearance of edges. Our solution is based on a regularized optimization problem for which we prove existence and uniqueness of the optimal solution.

preprint2020arXiv

Low-rank Kalman filtering under model uncertainty

We consider a robust filtering problem where the nominal state space model is not reachable and different from the actual one. We propose a robust Kalman filter which solves a dynamic game: one player selects the least-favorable model in a given ambiguity set, while the other player designs the optimum filter for the least-favorable model. It turns out that the robust filter is governed by a low-rank risk sensitive-like Riccati equation. Finally, simulation results show the effectiveness of the proposed filter.

preprint2020arXiv

M$^2$-Spectral Estimation: A Relative Entropy Approach

This paper deals with M$^2$-signals, namely multivariate (or vector-valued) signals defined over a multidimensional domain. In particular, we propose an optimization technique to solve the covariance extension problem for stationary random vector fields. The multidimensional Itakura-Saito distance is employed as an optimization criterion to select the solution among the spectra satisfying a finite number of moment constraints. In order to avoid technicalities that may happen on the boundary of the feasible set, we deal with the discrete version of the problem where the multidimensional integrals are approximated by Riemann sums. The spectrum solution is also discrete, which occurs naturally when the underlying random field is periodic. We show that a solution to the discrete problem exists, is unique and depends smoothly on the problem data. Therefore, we have a well-posed problem whose solution can be tuned in a smooth manner. Finally, we have applied our theory to the target parameter estimation problem in an integrated system of automotive modules. Simulation results show that our spectral estimator has promising performance.

preprint2016arXiv

Online semi-parametric learning for inverse dynamics modeling

This paper presents a semi-parametric algorithm for online learning of a robot inverse dynamics model. It combines the strength of the parametric and non-parametric modeling. The former exploits the rigid body dynamics equa- tion, while the latter exploits a suitable kernel function. We provide an extensive comparison with other methods from the literature using real data from the iCub humanoid robot. In doing so we also compare two different techniques, namely cross validation and marginal likelihood optimization, for estimating the hyperparameters of the kernel function.

preprint2016arXiv

Robust Kalman Filtering under Model Perturbations

We consider a family of divergence-based minimax approaches to perform robust filtering. The mismodeling budget, or tolerance, is specified at each time increment of the model. More precisely, all possible model increments belong to a ball which is formed by placing a bound on the Tau-divergence family between the actual and the nominal model increment. Then, the robust filter is obtained by minimizing the mean square error according to the least favorable model in that ball. It turns out that the solution is a family of Kalman like filters. Their gain matrix is updated according to a risk sensitive like iteration where the risk sensitivity parameter is now time varying. As a consequence, we also extend the risk sensitive filter to a family of risk sensitive like filters according to the Tau-divergence family.

preprint2016arXiv

Sparse plus Low rank Network Identification: A Nonparametric Approach

Modeling and identification of high-dimensional stochastic processes is ubiquitous in many fields. In particular, there is a growing interest in modeling stochastic processes with simple and interpretable structures. In many applications, such as econometrics and biomedical sciences, it seems natural to describe each component of that stochastic process in terms of few factor variables, which are not accessible for observation, and possibly of few other components of the stochastic process. These relations can be encoded in graphical way via a structured dynamic network, referred to as "sparse plus low-rank (S+L) network" hereafter. The problem of finding the S+L network as well as the dynamic model can be posed as a system identification problem. In this paper, we introduce two new nonparametric methods to identify dynamic models for stochastic processes described by a S+L network. These methods take inspiration from regularized estimators based on recently introduced kernels (e.g. "stable spline", "tuned-correlated" etc.). Numerical examples show the benefit to introduce the S+L structure in the identification procedure.

preprint2015arXiv

An Interpretation of the Dual Problem of the THREE-like Approaches

Spectral estimation can be preformed using the so called THREE-like approach. Such method leads to a convex optimization problem whose solution is characterized through its dual problem. In this paper, we show that the dual problem can be seen as a new parametric spectral estimation problem. This interpretation implies that the THREE-like solution is optimal in terms of closeness to the correlogram over a certain parametric class of spectral densities, enriching in this way its meaningfulness.

preprint2015arXiv

Sparse plus low-rank autoregressive identification in neuroimaging time series

This paper considers the problem of identifying multivariate autoregressive (AR) sparse plus low-rank graphical models. Based on the corresponding problem formulation recently presented, we use the alternating direction method of multipliers (ADMM) to efficiently solve it and scale it to sizes encountered in neuroimaging applications. We apply this decomposition on synthetic and real neuroimaging datasets with a specific focus on the information encoded in the low-rank structure of our model. In particular, we illustrate that this information captures the spatio-temporal structure of the original data, generalizing classical component analysis approaches.

preprint2014arXiv

AR Identification of Latent-variable Graphical Models

The paper proposes an identification procedure for autoregressive gaussian stationary stochastic processes wherein the manifest (or observed) variables are mostly related through a limited number of latent (or hidden) variables. The method exploits the sparse plus low-rank decomposition of the inverse of the manifest spectral density and the efficient convex relaxations recently proposed for such decomposition.

preprint2014arXiv

Multivariate Spectral Estimation based on the concept of Optimal Prediction

In this technical note, we deal with a spectrum approximation problem arising in THREE-like multivariate spectral estimation approaches. The solution to the problem minimizes a suitable divergence index with respect to an a priori spectral density. We derive a new divergence family between multivariate spectral densities which takes root in the prediction theory. Under mild assumptions on the a priori spectral density, the approximation problem, based on this new divergence family, admits a family of solutions. Moreover, an upper bound on the complexity degree of these solutions is provided.

preprint2013arXiv

A Contraction Analysis of the Convergence of Risk-Sensitive Filters

A contraction analysis of risk-sensitive Riccati equations is proposed. When the state-space model is reachable and observable, a block-update implementation of the risk-sensitive filter is used to show that the N-fold composition of the Riccati map is strictly contractive with respect to the Riemannian metric of positive definite matrices, when N is larger than the number of states. The range of values of the risk-sensitivity parameter for which the map remains contractive can be estimated a priori. It is also found that a second condition must be imposed on the risk-sensitivity parameter and on the initial error variance to ensure that the solution of the risk-sensitive Riccati equation remains positive definite at all times. The two conditions obtained can be viewed as extending to the multivariable case an earlier analysis of Whittle for the scalar case.

preprint2013arXiv

A new family of high-resolution multivariate spectral estimators

In this paper, we extend the Beta divergence family to multivariate power spectral densities. Similarly to the scalar case, we show that it smoothly connects the multivariate Kullback-Leibler divergence with the multivariate Itakura-Saito distance. We successively study a spectrum approximation problem, based on the Beta divergence family, which is related to a multivariate extension of the THREE spectral estimation technique. It is then possible to characterize a family of solutions to the problem. An upper bound on the complexity of these solutions will also be provided. Simulations suggest that the most suitable solution of this family depends on the specific features required from the estimation problem.

preprint2013arXiv

Minimum Relative Entropy for Quantum Estimation: Feasibility and General Solution

We propose a general framework for solving quantum state estimation problems using the minimum relative entropy criterion. A convex optimization approach allows us to decide the feasibility of the problem given the data and, whenever necessary, to relax the constraints in order to allow for a physically admissible solution. Building on these results, the variational analysis can be completed ensuring existence and uniqueness of the optimum. The latter can then be computed by standard, efficient standard algorithms for convex optimization, without resorting to approximate methods or restrictive assumptions on its rank.

preprint2013arXiv

Rational approximations of spectral densities based on the Alpha divergence

We approximate a given rational spectral density by one that is consistent with prescribed second-order statistics. Such an approximation is obtained by minimizing a suitable distance from the given spectrum and under the constraints corresponding to imposing the given second-order statistics. Here, we consider the Alpha divergence family as a distance measure. We show that the corresponding approximation problem leads to a family of rational solutions. Secondly, such a family contains the solution which generalizes the Kullback-Leibler solution proposed by Georgiou and Lindquist in 2003. Finally, numerical simulations suggest that this family contains solutions close to the non-rational solution given by the principle of minimum discrimination information.