Source author record

Jameson Cahill

Jameson Cahill 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

18works
6topics
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

18 published item(s)

preprint2022arXiv

Group-invariant max filtering

Given a real inner product space $V$ and a group $G$ of linear isometries, we construct a family of $G$-invariant real-valued functions on $V$ that we call max filters. In the case where $V=\mathbb{R}^d$ and $G$ is finite, a suitable max filter bank separates orbits, and is even bilipschitz in the quotient metric. In the case where $V=L^2(\mathbb{R}^d)$ and $G$ is the group of translation operators, a max filter exhibits stability to diffeomorphic distortion like that of the scattering transform introduced by Mallat. We establish that max filters are well suited for various classification tasks, both in theory and in practice.

preprint2016arXiv

Connectivity and Irreducibility of Algebraic Varieties of Finite Unit Norm Tight Frames

In this paper, we settle a long-standing problem on the connectivity of spaces of finite unit norm tight frames (FUNTFs), essentially affirming a conjecture first appearing in [Dykema and Strawn, 2003]. Our central technique involves continuous liftings of paths from the polytope of eigensteps to spaces of FUNTFs. After demonstrating this connectivity result, we refine our analysis to show that the set of nonsingular points on these spaces is also connected, and we use this result to show that spaces of FUNTFs are irreducible in the algebro-geometric sense, and also that generic FUNTFs are full spark.

preprint2016arXiv

The gap between the null space property and the restricted isometry property

The null space property (NSP) and the restricted isometry property (RIP) are two properties which have received considerable attention in the compressed sensing literature. As the name suggests, NSP is a property that depends solely on the null space of the measurement procedure and as such, any two matrices which have the same null space will have NSP if either one of them does. On the other hand, RIP is a property of the measurement procedure itself, and given an RIP matrix it is straightforward to construct another matrix with the same null space that is not RIP. %Furthermore, RIP is known to imply NSP and therefore RIP is a strictly stronger assumption than NSP. We say a matrix is RIP-NSP if it has the same null space as an RIP matrix. We show that such matrices can provide robust recovery of compressible signals under Basis pursuit which in many applicable settings is comparable to the guarantee that RIP provides. More importantly, we constructively show that the RIP-NSP is stronger than NSP with the aid of this robust recovery result, which shows that RIP is fundamentally stronger than NSP.

preprint2015arXiv

Phase retrieval

We answer a number of open problems concerning phase retrieval and phase retrieval by projections. In particular, one main theorem classifies phase retrieval by projections via collections of sequences of vectors allowing norm retrieval. Another key result computes the minimal number of vectors needed to add to a frame in order for it to possess the complement property and hence allow phase retrieval. In furthering this idea, in a third main theorem we show that when a collection of subspaces is one subspace short from allowing phase retrieval, then any partition of orthonormal bases from these subspaces into two sets which fail to span, then each spans a hyperplane. We offer many more results in this area as well as provide a large number of examples showing the limitations of the theory.

preprint2014arXiv

Phase retrieval and norm retrieval

Phase retrieval has become a very active area of research. We will classify when phase retrieval by Parseval frames passes to the Naimark complement and when phase retrieval by projections passes to the orthogonal complements. We introduce a new concept we call norm retrieval and show that this is what is necessary for passing phase retrieval to complements. This leads to a detailed study of norm retrieval and its relationship to phase retrieval. One fundamental result: a frame $\{φ_i\}_{i=1}^M$ yields phase retrieval if and only if $\{Tφ_i\}_{i=1}^M$ yields norm retrieval for every invertible operator $T$.

preprint2014arXiv

Robust width: A characterization of uniformly stable and robust compressed sensing

Compressed sensing seeks to invert an underdetermined linear system by exploiting additional knowledge of the true solution. Over the last decade, several instances of compressed sensing have been studied for various applications, and for each instance, reconstruction guarantees are available provided the sensing operator satisfies certain sufficient conditions. In this paper, we completely characterize the sensing operators which allow uniformly stable and robust reconstruction by convex optimization for many of these instances. The characterized sensing operators satisfy a new property we call the robust width property, which simultaneously captures notions of widths from approximation theory and of restricted eigenvalues from statistical regression. We provide a geometric interpretation of this property, we discuss its relationship with the restricted isometry property, and we apply techniques from geometric functional analysis to find random matrices which satisfy the property with high probability.

preprint2013arXiv

A note on scalable frames

We study the problem of determining whether a given frame is scalable, and when it is, understanding the set of all possible scalings. We show that for most frames this is a relatively simple task in that the frame is either not scalable or is scalable in a unique way, and to find this scaling we just have to solve a linear system. We also provide some insight into the set of all scalings when there is not a unique scaling. In particular, we show that this set is a convex polytope whose vertices correspond to minimal scalings.

preprint2013arXiv

Phase Retrieval By Projections

The problem of recovering a vector from the absolute values of its inner products against a family of measurement vectors has been well studied in mathematics and engineering. A generalization of this phase retrieval problem also exists in engineering: recovering a vector from measurements consisting of norms of its orthogonal projections onto a family of subspaces. There exist semidefinite programming algorithms to solve this problem, but much remains unknown for this more general case. Can families of subspaces for which such measurements are injective be completely classified? What is the minimal number of subspaces required to have injectivity? How closely does this problem compare to the usual phase retrieval problem with families of measurement vectors? In this paper, we answer or make incremental steps toward these questions. We provide several characterizations of subspaces which yield injective measurements, and through a concrete construction, we prove the surprising result that phase retrieval can be achieved with $2M-1$ projections of arbitrary rank in $\HH_M$. Finally we present several open problems as we discuss issues unique to the phase retrieval problem with subspaces.

preprint2013arXiv

Saving phase: Injectivity and stability for phase retrieval

Recent advances in convex optimization have led to new strides in the phase retrieval problem over finite-dimensional vector spaces. However, certain fundamental questions remain: What sorts of measurement vectors uniquely determine every signal up to a global phase factor, and how many are needed to do so? Furthermore, which measurement ensembles lend stability? This paper presents several results that address each of these questions. We begin by characterizing injectivity, and we identify that the complement property is indeed a necessary condition in the complex case. We then pose a conjecture that 4M-4 generic measurement vectors are both necessary and sufficient for injectivity in M dimensions, and we prove this conjecture in the special cases where M=2,3. Next, we shift our attention to stability, both in the worst and average cases. Here, we characterize worst-case stability in the real case by introducing a numerical version of the complement property. This new property bears some resemblance to the restricted isometry property of compressed sensing and can be used to derive a sharp lower Lipschitz bound on the intensity measurement mapping. Localized frames are shown to lack this property (suggesting instability), whereas Gaussian random measurements are shown to satisfy this property with high probability. We conclude by presenting results that use a stochastic noise model in both the real and complex cases, and we leverage Cramer-Rao lower bounds to identify stability with stronger versions of the injectivity characterizations.

preprint2013arXiv

Tight and random nonorthogonal fusion frames

First we show that tight nonorthogonal fusion frames a relatively easy to com by. In order to do this we need to establish a classification of how to to wire a self adjoint operator as a product of (nonorthogonal) projection operators. We also discuss the link between nonorthogonal fusion frames and positive operator valued measures, we define and study a nonorthogonal fusion frame potential, and we introduce the idea of random nonorthogonal fusion frames.

preprint2012arXiv

Full Spark Frames

Finite frame theory has a number of real-world applications. In applications like sparse signal processing, data transmission with robustness to erasures, and reconstruction without phase, there is a pressing need for deterministic constructions of frames with the following property: every size-M subcollection of the M-dimensional frame elements is a spanning set. Such frames are called full spark frames, and this paper provides new constructions using the discrete Fourier transform. Later, we prove that full spark Parseval frames are dense in the entire set of Parseval frames, meaning full spark frames are abundant, even if one imposes an additional tightness constraint. Finally, we prove that testing whether a given matrix is full spark is hard for NP under randomized polynomial-time reductions, indicating that deterministic full spark constructions are particularly significant because they guarantee a property which is otherwise difficult to check.

preprint2011arXiv

Constructing finite frames of a given spectrum and set of lengths

When constructing finite frames for a given application, the most important consideration is the spectrum of the frame operator. Indeed, the minimum and maximum eigenvalues of the frame operator are the optimal frame bounds, and the frame is tight precisely when this spectrum is constant. Often, the second-most important design consideration is the lengths of frame vectors: Gabor, wavelet, equiangular and Grassmannian frames are all special cases of equal norm frames, and unit norm tight frame-based encoding is known to be optimally robust against additive noise and erasures. We consider the problem of constructing frames whose frame operator has a given spectrum and whose vectors have prescribed lengths. For a given spectrum and set of lengths, the existence of such frames is characterized by the Schur-Horn Theorem---they exist if and only if the spectrum majorizes the squared lengths---the classical proof of which is nonconstructive. Certain construction methods, such as harmonic frames and spectral tetris, are known in the special case of unit norm tight frames, but even these provide but a few examples from the manifold of all such frames, the dimension of which is known and nontrivial. In this paper, we provide a new method for explicitly constructing any and all frames whose frame operator has a prescribed spectrum and whose vectors have prescribed lengths. The method itself has two parts. In the first part, one chooses eigensteps---a sequence of interlacing spectra---that transform the trivial spectrum into the desired one. The second part is to explicitly compute the frame vectors in terms of these eigensteps; though nontrivial, this process is nevertheless straightforward enough to be implemented by hand, involving only arithmetic, square roots and matrix multiplication.

preprint2011arXiv

The Paulsen Problem in Operator Theory

The Paulsen Problem in Hilbert space frame theory has proved to be one of the most intractable problems in the field. We will help explain why by showing that this problem is equivalent to a fundamental, deep problem in operator theory. Along the way we will give a new exact computation for chordal distances, we will give a generalization of these problems and we will spell out exactly the complementary versions of the problem.

preprint2010arXiv

Non-orthogonal fusion frames and the sparsity of fusion frame operators

Fusion frames have become a major tool in the implementation of distributed systems. The effectiveness of fusion frame applications in distributed systems is reflected in the efficiency of the end fusion process. This in turn is reflected in the efficiency of the inversion of the fusion frame operator $S_{\cW}$, which in turn is heavily dependent on the sparsity of $S_{\cW}$. We will show that sparsity of the fusion frame operator naturally exists by introducing a notion of {\it non-orthogonal fusion frames}. We show that for a fusion frame $\{W_i,v_i\}_{i\in I}$, if $\text{dim}(W_i)=k_i$, then the matrix of the non-orthogonal fusion frame operator $\cSw$ has in its corresponding location at most a $k_i\times k_i$ block matrix. We provide necessary and sufficient conditions for which the new fusion frame operator $\cSw$ is diagonal and/or a multiple of an identity. A set of other critical questions are also addressed. A scheme of {\it multiple fusion frames} whose corresponding fusion frame operator becomes an diagonal operator is also examined.