On Reconstructing Gaussian Mixtures from the Distance Between Two Samples: An Algebraic Perspective

Kindyl King, Purdue University

Abstract

This thesis is concerned with the problem of characterizing the orbits of certain probability density functions under the action of the Euclidean group. Our motivating application is the recognition of a point configuration where the coordinates of the points are measured under noisy conditions. Consider a random variable X in R¢ with probability density function p(.e). Let x and rz be independent random samples following p(x). Define A as the squared Euelidean distance between x; and zz. It has previously been shown that two distributions p(x) and p(x) consisting of Dirac delta distributions in generic positions that have the same respective distributions of A are necessarily related by a rigid motion. That is, there exists some rigid motion g in the Euclidean group p(x) such that p(x) = P(g - x) for all x € Rd. To account for noise in the measurements. we assume X is a random variable in R¢ whose density is a k-component mixture of Gaussian distributions with means in generic position. We further assume that the covariance matrices of the Gaussian components are equal and of the form © = 971, with o? € Roo. In Theorem 3.1.1 and Theorem 3.2.1, we prove that, when o? is known, generic k- component Gaussian mixtures are uniquely reconstructible up to a rigid motion from the density of ∆. A more general formulation is proven in Theorem 3.2.3. Similarly, when is unknown, we prove in Theorem 4.1.1 and Theorem 4.1.2 that generic equally-weighted k-component Gaussian mixtures with & = 1 and k = 2 are uniquely reconstructible up to a rigid motion from the distribution of ∆. There are at most three non-equivalent equally- weighted 3-component Gaussian mixtures up to a rigid motion having the same distribution of A, as proven in Theorem 4.1.3. In Theorem 1.1.4, we present. a test to check if, for a given & and d, the number of non-equivalent equally-weighted k-component Gaussian mixtures in R¢ having the same distribution of ∆ is at most (3) +1. Numerical computations showed that distributions with k = 4,5.6,7 such that d < k—2 and (k,d) = (8,1) pass the test, and thus have a finite number of reconstructions up to a rigid motion. When a2 is unknown and the mixture weights are also unknown, we prove in Theorem 1.2.) that there are at most four non-equivalent 2-component Gaussian mixtures up to a rigid motion having the same distribution of ∆.

Degree

Ph.D.

Advisors

Boutin, Purdue University.

Subject Area

Mathematics

COinS