Stability estimates for phase retrieval from discrete Gabor measurements

TL;DR

This paper investigates the stability of phase retrieval from discrete Gabor measurements, demonstrating that semi-global reconstruction can achieve polynomially bounded stability constants in finite-dimensional, bandlimited signal settings.

Contribution

It provides the first evidence that semi-global phase retrieval from discrete Gabor measurements in finite dimensions can have stability constants that grow only polynomially with dimension.

Findings

Stability constant scales polynomially with dimension

Semi-global regime enables stable phase retrieval for bandlimited signals

Utilizes recent reconstruction formulae for analysis

Abstract

Phase retrieval refers to the problem of recovering some signal (which is often modelled as an element of a Hilbert space) from phaseless measurements. It has been shown that in the deterministic setting phase retrieval from frame coefficients is always unstable in infinite-dimensional Hilbert spaces [7] and possibly severely ill-conditioned in finite-dimensional Hilbert spaces [7]. Recently, it has also been shown that phase retrieval from measurements induced by the Gabor transform with Gaussian window function is stable under a more relaxed semi-global phase recovery regime based on atoll functions [1]. In finite dimensions, we present first evidence that this semi-global reconstruction regime allows one to do phase retrieval from measurements of bandlimited signals induced by the discrete Gabor transform in such a way that the corresponding stability constant only scales like a low order polynomial in the space dimension. To this end, we utilise reconstruction formulae which have become common tools in recent years [6,12,18,20].

Equations

$$f_{\lambda }^{+}(t):=g(t-\lambda )+g(t+\lambda )\text{mbox}and{f}_{\lambda }^{-}(t):=g(t-\lambda )-g(t+\lambda )$$

$${\mathcal{V}}_{\varphi }[x](m,n):=\frac{1}{\sqrt{L}}\cdot \sum _{\ell =0}^{L-1}x(\ell )\overline{\varphi (\ell -m)}{\mathrm{e}}^{-2\pi \mathrm{i}\frac{\ell n}{L}},m,n=0,\dots ,L-1.$$

$${\mathcal{V}}_{\varphi }[x](m,n)=\mathcal{F}\left[x_m\right](n),m,n=0,\dots ,L-1,$$

$$\mathcal{F}[x](k):=\frac{1}{\sqrt{L}}\cdot \sum _{\ell =0}^{L-1}x(\ell ){\mathrm{e}}^{-2\pi \mathrm{i}\frac{\ell k}{L}},k=0,\dots ,L-1,$$

$${\mathcal{F}}^{-1}[x](\ell )=\frac{1}{\sqrt{L}}\cdot \sum _{k=0}^{L-1}x(k){\mathrm{e}}^{2\pi \mathrm{i}\frac{k\ell }{L}},\ell =0,\dots ,L-1.$$

$$M_{\varphi }[x](m,n):={\left|{\mathcal{V}}_{\varphi }[x](m,n)\right|}^2,m,n=0,\dots ,L-1.$$

$$\underset{\alpha \in \mathbb{R}}{\inf }\Vert x-{\mathrm{e}}^{\mathrm{i}\alpha }y{\Vert }_2\lesssim \Vert \left|{\mathcal{V}}_{\varphi }[x]\right|-\left|{\mathcal{V}}_{\varphi }[y]\right|{\Vert }_{\mathrm{F}},$$

$${\mathcal{F}}^{-1}\left[M_{\varphi }[x](m,\cdot )\right](n)=\frac{1}{\sqrt{L}}\cdot \left(x_m\ast x_m^{\#}\right)(n),m,n=0,\dots ,L-1.$$

$$\left(x_m\ast x_m^{\#}\right)(n)=\sum _{\ell =0}^{L-1}x(\ell )\overline{x(\ell -n)}\varphi (\ell -n-m)\overline{\varphi (\ell -m)},m,n=0,\dots ,L-1.$$

$$\mathcal{F}\left[M_{\varphi }[x]\right](m,n)=\mathcal{A}[x](-n,m)\overline{\mathcal{A}[\varphi ](-n,m)},\text{mbox}form,n=0,\dots ,L-1.$$

$$\mathcal{A}[\varphi ](m,n)\ne 0,\text{mbox}form,n=0,\dots ,L-1.$$

$$\mathcal{A}[\varphi ](m,n)\ne 0,\text{mbox}form=0,1,n=0,\dots ,L-1.$$

$$\mathcal{A}[\varphi ](m,n)\ne 0,\text{mbox}form\in \{0,\dots ,\Delta +1\},n\in \{0,\dots ,L-1\}.$$

$$(\ell ,k)\in E\iff |\ell -k|\in (0,\Delta +1]\cup [L-\Delta -1,L),$$

$$x(\ell )={\mathrm{e}}^{\mathrm{i}{\alpha }_k}y(\ell ),\ell \in V_k.$$

$$\mathrm{supp}\varphi =\{n_0,\dots ,n_0+{\ell }_{\varphi }\}$$

$$|\ell -k|\in \{{\ell }_{\varphi },L-{\ell }_{\varphi }\}$$

$$V_{\delta }=\{\ell \in \{0,\dots ,L-1\}||x(\ell )|\ge \delta \}$$

$$|\ell -k|\in (0,\Delta +1]\cup [L-\Delta -1,L),$$

$$\mathrm{supp}\mathcal{F}[x]\subset \{-B,\dots ,B\}\mathrm{mod}L.$$

$$\mathrm{supp}\mathcal{A}[x](m,\cdot )\subset \{-2B,\dots ,2B\}\mathrm{mod}L.$$

$$\Vert x{\Vert }_{{\ell }^2(S)}:={\left(\sum _{\ell \in S}|x(\ell ){|}^2\right)}^{\frac{1}{2}}.$$

$$\underset{n\in \{-2B,\dots ,2B\}}{\min }|\mathcal{A}[\varphi ](0,n)|\ge \frac{1}{c},$$

$${\Vert |x|-|y|\Vert }_{{\ell }^2(V_{\delta })}\le \frac{c}{\delta }{\Vert M_{\varphi }[x]-M_{\varphi }[y]\Vert }_{\mathrm{F}}.$$

$$\left||x(\ell )|-|y(\ell )|\right|=\frac{\left||x(\ell ){|}^2-|y(\ell ){|}^2\right|}{\left||x(\ell )|+|y(\ell )|\right|}\le \frac{1}{\delta }\left||x(\ell ){|}^2-|y(\ell ){|}^2\right|.$$

$${\Vert |x|-|y|\Vert }_{{\ell }^2(V_{\delta })}\le \frac{1}{\delta }{\Vert |x{|}^2-|y{|}^2\Vert }_{{\ell }^2(V_{\delta })}\le \frac{1}{\delta }{\Vert |x{|}^2-|y{|}^2\Vert }_2.$$

$${\Vert |x{|}^2-|y{|}^2\Vert }_2^2$$

$$=\sum _{n=-2B}^{2B}{\left|\mathcal{A}[x](0,n)-\mathcal{A}[y](0,n)\right|}^2=\sum _{n=-2B}^{2B}{\left|\frac{\mathcal{F}\left[M_{\varphi }[x]-M_{\varphi }[y]\right](n,0)}{\mathcal{A}[\varphi ](0,n)}\right|}^2$$

$$\le c^2\sum _{n=-2B}^{2B}{\left|\mathcal{F}\left[M_{\varphi }[x]-M_{\varphi }[y]\right](n,0)\right|}^2\le c^2{\Vert \mathcal{F}\left[M_{\varphi }[x]-M_{\varphi }[y]\right]\Vert }_2^2$$

$$=c^2{\Vert M_{\varphi }[x]-M_{\varphi }[y]\Vert }_2^2.$$

$${\Vert |x{|}^2-|y{|}^2\Vert }_2\le c{\Vert M_{\varphi }[x]-M_{\varphi }[y]\Vert }_2.$$

Full text

Recently, it has also been shown that phase retrieval from measurements induced by the Gabor transform with Gaussian window function is stable under a more relaxed semi-global phase recovery regime based on atoll functions [1].

In finite dimensions, we present first evidence that this semi-global reconstruction regime allows one to do phase retrieval from measurements of bandlimited signals induced by the discrete Gabor transform in such a way that the corresponding stability constant only scales like a low order polynomial in the space dimension. To this end, we utilise reconstruction formulae which have become common tools in recent years [6, 13, 19, 21].

Keywords Phase retrieval, Gabor frames, Time-frequency analysis

Mathematics Subject Classification 42C15, 42A38, 94A12, 65T50

1 Introduction

Phase retrieval generally alludes to the non-linear inverse problem of recovering some signal (which in this paper will be modelled by $x\in\mathbb{C}^{L}$) from phaseless measurements. Some of its more well-known applications include ptychography for coherent diffraction imaging [16, 20, 24, 29] and audio processing [14, 18, 12]. It has been shown that the phase retrieval problem for frames in finite-dimensional Hilbert spaces [7] and a forteriori in finite-dimensional reflexive Banach spaces [2] is always stable, which elicits the question: Why are we concerned with stability estimates for phase retrieval from discrete Gabor measurements at all? The reason is that phase retrieval for frames in infinite-dimensional spaces is always unstable [7, 2] and in addition one can construct sequences of finite-dimensional subspaces of infinite-dimensional Hilbert spaces along with frames for which the stability constant of phase retrieval increases exponentially in the dimension of the constructed subspaces [7]. Recent research [1] into the infinite-dimensional phase retrieval problem has however led us to believe that the instability of phase retrieval is not an insurmountable obstacle to reconstruction. It was shown that stability can be restored for examples that exhibit a disconnectedness in the measurements by only reconstructing the phase semi-globally or in an atoll sense. Furthermore, it was shown in [15] that such disconnectedness in the measurements is the only source of instabilities for phase retrieval.

A simple example of instability can be obtained by considering the Gaussian functions $g(t):=\mathrm{e}^{-\pi t^{2}}$ in conjunction with the signals

[TABLE]

depicted in figure 1. When $\lambda$ increases, the Gaussian bumps in the signals $f_{\lambda}^{\pm}$ start to move further apart effectively generating what we call a time gap whose length depends linearly on $\lambda$. It can be shown, see [3], that the measurements generated by the continuous Gabor transform with Gaussian window of the signals $f_{\lambda}^{\pm}$ have distance on the order of $\mathrm{e}^{-\lambda^{2}}$ in the standard Sobolev space $W^{1,2}(\mathbb{R}^{2})$ and that one can therefore not stably retrieve $f_{\lambda}^{\pm}$ from continuous Gabor transform measurements. Similar phenomena can be observed for the discrete setting considered in this paper and we do therefore propose a similar paradigm as in [1] and try to recover signals in a semi-global fashion that is not common in the phase retrieval literature up to this point. Note that in audio processing, it is natural to consider (audio) signals up to semi-global phase as human listeners are not able to distinguish between two signals which differ by semi-global phase [1].

One should note that in recent years a variety of stability result for phase retrieval have been proven. Some highlights of this research include:

1. i.

The PhaseLift method [8, 11] which guarantees stable recovery from $\mathcal{O}(L)$ randomly chosen Gaussian measurements with high probability [9]. 2. ii.

The research on polarisation for phase retrieval [4, 5, 22, 26] in which the authors supplement an existing measurement ensemble in order to obtain a phase retrieval problem that is efficiently and stably solvable. 3. iii.

Wirtinger flow and related methods [10, 27, 28] which offer stability guarantees for sufficiently many randomly chosen Gaussian measurements. 4. iv.

The eigenvector-based angular synchronisation approach [17] which relies on a certain weak form of invertibility of the phase retrieval problem to prove a stability result for deterministic measurement systems. 5. v.

The very recent work [25] in which the stability of phase retrieval from (random) frames whose frame vectors are uniformly distributed on the unit sphere (but not necessarily independent) is considered.

In some way or another, all of these results are based on different setups than ours: As opposed to the papers referenced in item i., iii. and v. we will not work with a probabilistic measurement system but with a deterministic one. We will also not supplement our measurement ensemble as is done in the results referenced in item ii. and we will not work with the weak form of invertibility that is present in the paper referenced in item iv. In fact, we will consider the two well-known formulae (1) and (2) presented in section 2 which are heavily used to develop methods for exact phase retrieval from Gabor measurements in the literature [13, 19, 21]. We show that through further analysis of the formulae (1) and (2), one can derive stability results for some of those methods and therefore also for phase retrieval in general. Our stability results are designed for bandlimited signals and come with constants that scale in the square root of the space dimension at the cost of relaxing the notion of stability to resemble the one proposed in [1].

Outline

In section 2, we present the reader with the uniqueness results and the formulae on which our stability results hinge. In section 3, we utilise the ambiguity function relation (2) in order to show that phase retrieval can be stably done for bandlimited signals based on the considerations in [6, 23]. In section 4, we use the autocorrelation relation (1) in order to show that phase retrieval can be done stably for bandlimited signals utilising results from [13, 19]. As the proofs of our main results are a bit technical, they appear separately in section 5.

2 Prerequisites

Throughout this paper, we fix the dimension $L\in\mathbb{N}$ and let $x,\phi\in\mathbb{C}^{L}$. We define the discrete Gabor transform (DGT) of $x$ with window function $\phi$ to be

[TABLE]

Here and throughout this paper, the indexing is understood to be periodic. In particular, we use the convention $\phi(\ell)=\phi(\ell\mod L)$, for $\ell\in\mathbb{Z}$. A helpful way of looking at the DGT is to view it as a collection of windowed Fourier transforms. For this purpose, we denote $x_{m}(\ell):=x(\ell)\overline{\phi(\ell-m)}$, for $\ell,m\in\{0,\dots,L-1\}$, and obtain

[TABLE]

where $\mathcal{F}:\mathbb{C}^{L}\to\mathbb{C}^{L}$ denotes the discrete Fourier transform (DFT)

[TABLE]

with inverse

[TABLE]

We will frequently use the two-dimensional discrete Fourier transform which is the composition of two DFTs as defined above. Additionally, we define the ambiguity function of a signal $x$ via $\mathcal{A}[x]:=\mathcal{V}_{x}[x]$. We are interested in the recovery of signals $x\in\mathbb{C}^{L}$ from the measurements

[TABLE]

It is immediately obvious that $x\in\mathbb{C}^{L}$ and any signal $\mathrm{e}^{\mathrm{i}\alpha}x$, with $\alpha\in\mathbb{R}$, yield the same measurements $M_{\phi}[\mathrm{e}^{\mathrm{i}\alpha}x]=M_{\phi}[x]$. Therefore, to have any chance of recovery, we will actually view $M_{\phi}$ as an operator defined on the quotient space $\mathbb{C}^{L}/\mathcal{S}^{1}$, where $\mathcal{S}^{1}$ denotes the unit circle. Under various assumptions, which we will lay out in the following, one can show that $M_{\phi}:\mathbb{C}^{L}/\mathcal{S}^{1}\to\mathbb{R}_{+}^{L\times L}$ is an injective operator and that phase retrieval is therefore possible up to a global phase factor. In addition, it was shown in [7] that

[TABLE]

for all $x,y\in\mathbb{C}^{L}$, where $\lVert\cdot\rVert_{\mathrm{F}}$ denotes the Frobenius norm and the estimate depends on a constant which might increase exponentially in the space dimension $L$. Our phase retrieval problem is therefore possibly ill-conditioned.

As mentioned before, the number of known uniqueness results has seen a stark rise in the past few years. In the following, we want to mention those that inspired our stability estimates. Let us start by remarking that almost all uniqueness results can be traced back to two consequential formulae which are well-known in the literature. The first of these relates the Gabor measurements to the autocorrelation of $x_{m}$. In what follows, time-reversal of a signal will be denoted by $x^{\#}(\ell)=\overline{x(-\ell)}$.

Lemma 2.1**.**

For any $x\in\mathbb{C}^{L}$,

[TABLE]

Proof.

See appendix B. ∎

The right-hand side in the above result is the aforementioned autocorrelation of $x_{m}$:

[TABLE]

The second of these formulae relates the Gabor measurements to the ambiguity function of $x$ and the ambiguity function of $\phi$.

Lemma 2.2**.**

For any $x\in\mathbb{C}^{L}$, the following holds:

[TABLE]

Proof.

See appendix B. ∎

Next, we will consider the uniqueness results from [6, 23] which are based on equation (2).

Corollary 2.3** (Theorem 2.2 in [6], p. 547).**

Suppose that $\phi$ satisfies

[TABLE]

Then, $x$ is uniquely determined by the measurements $M_{\phi}[x]$ up to global phase.

While this result is exceptionally nice in the sense that it does not impose any requirements on the signal, it is quite restrictive in its requirements on the window function $\phi$. For instance, windows $\phi$ with support length $\lvert\operatorname{supp}\phi\rvert$ smaller than $L/2$ will always have zero entries in their ambiguity function.

Corollary 2.4** (Theorem 2.4 in [6], p. 549).**

Let $x\in\mathbb{C}^{L}$ be nowhere-vanishing, i.e. $\operatorname{supp}x=\{0,\dots,L-1\}$, and

[TABLE]

Then, $x$ is uniquely determined by the measurements $M_{\phi}[x]$ up to global phase.

This result is in some sense orthogonal to corollary 2.3: Its requirements on the window function are moderate while its requirements on the signal are rather restrictive. Of course, we might also infer a variety of results that are based on different trade-offs between restrictions on the window and restrictions on the signal. For this purpose, we introduce the parameter $\Delta\in\mathbb{N}_{0}$. It corresponds to the maximum number of adjacent zeroes across which we may propagate phase in the reconstruction procedure used in the proof of the following corollary. Stated a bit more precisely: If $x$ is a signal of which we only know its measurements $M_{\phi}[x]$, then it follows from $\mathcal{A}[\phi](0,\cdot)$ being nowhere-vanishing (and the use of the ambiguity function relation) that we can reconstruct the magnitudes of $x$. Therefore, it suffices to propagate phase between the entries of $x$ to reconstruct $x$ up to global phase. When we assume that $\mathcal{A}[\phi](m,\cdot)$ is nowhere-vanishing, for some $m$, then we allow (according to the ambiguity function relation) the phase propagation from the entry with index $\ell$ to the entry with index $\ell+m$ (and to the entry with index $\ell-m$). Whether this allows us to reconstruct $x$ up to global phase depends on the set of $m$ for which $\mathcal{A}[\phi](m,\cdot)$ is nowhere-vanishing and on the support set of $x$. This is the central idea on which the following corollary is built:

Corollary 2.5**.**

Let $\Delta\in\mathbb{N}_{0}$ and let $x,y,\phi\in\mathbb{C}^{L}$ be such that $M_{\phi}[x]=M_{\phi}[y]$ and

[TABLE]

Furthermore, let $G=(V,E)$ denote the graph with vertex set $V=\operatorname{supp}x$ and edge set $E\subset V\times V$ such that

[TABLE]

i.e. two vertices are connected if and only if they are at most $\Delta+1$ apart. If $\{V_{k}\}_{k=1}^{K}$ constitute the vertex sets of the connected components of $G$, then for each $k\in\{1,\dots,K\}$ there exists an $\alpha_{k}\in\mathbb{R}$ such that

[TABLE]

Proof.

See section 5. ∎

Remark 2.6**.**

The corollary above is more general than corollaries 2.3 and 2.4. Indeed, if $\Delta\geq\tfrac{L}{2}-1$, then $G$ is connected. In fact, one can see from the definition of the edge set that $G$ is the complete graph on the vertex set $\operatorname{supp}x$. In particular, corollary 2.3 follows. If $\Delta=0$ and $x$ is nowhere-vanishing, then $G$ is the circle graph on $L$ vertices and is thus connected. In this way, we recover corollary 2.4.

Finally, we will work with a uniqueness result first proven in [13] and later generalised in [19] based mostly on equation (1). Consider the following statement.

Corollary 2.7** (Theorem 1 in [13], p. 639).**

Let $n_{0},\ell_{\phi}\in\{0,\dots,L-1\}$ be such that $\ell_{\phi}<L/2$ and suppose that $\ell_{\phi}-1$ and $L$ are coprime. If

[TABLE]

and $\mathcal{F}[\lvert\phi\rvert^{2}]$ and $x$ are nowhere-vanishing, then $x$ is uniquely determined by the measurements $M_{\phi}[x]$ up to global phase.

The work in [19] shows that one can also derive this result as part of a graph-theoretical formulation for phase retrieval.

Corollary 2.8** (Theorem 3.1 in [19], p. 373).**

Let $n_{0},\ell_{\phi}\in\{0,\dots,L-1\}$ be such that $\ell_{\phi}<L/2$ and suppose that

[TABLE]

and $\mathcal{F}[\lvert\phi\rvert^{2}]$ is nowhere-vanishing. Let the graph $G=(V,E)$ defined by having the vertex set $V=\operatorname{supp}x$ and an edge between $\ell,k\in V$ if

[TABLE]

be connected. Then, $x$ is uniquely determined by the measurements $M_{\phi}[x]$ up to global phase.

3 Stability estimates based on the ambiguity function relation

3.1 Stability for a single island

First, we derive stability estimates by employing equation (2) and corollaries 2.3–2.5. In doing this, we want to start with the very simple setup of corollary 2.4.

One can immediately see that there are some intricacies to the phase retrieval problem for signals $x\in\mathbb{C}^{L}$. One of those is dealing with entries $x(\ell)$ of $x$ which have small (or even vanishing) magnitude. For these entries, extracting the phase of $x(\ell)$ is unstable (or even impossible). See figure 2 for a depiction of this situation. Because of this, we will mostly work with a graph capturing only the larger entries of the signals.

Definition 3.1**.**

Let $\Delta\in\mathbb{N}_{0}$ and $\delta>0$. We call the graph $G_{\delta}=(V_{\delta},E)$ defined by having the vertex set

[TABLE]

and an edge between $\ell,k\in V$ if

[TABLE]

the essential support graph of $x$ with time separation parameter $\Delta$. We will also simply call the essential support graph of $x$ with time separation parameter zero the essential support graph of $x$.

The stability estimates we derive hold for bandlimited signals defined as follows:

Definition 3.2**.**

Let $B\in\mathbb{N}_{0}$. We call $x\in\mathbb{C}^{L}$ $B$-bandlimited if

[TABLE]

One important property of bandlimited signals is that their ambguity function does not have full support (when $4B<L$).

Proposition 3.3**.**

Let $x\in\mathbb{C}^{L}$ be $B$-bandlimited, for some $B\in\mathbb{N}_{0}$. Then, it holds for all $m\in\{0,\dots,L-1\}$ that

[TABLE]

We included a proof of this basic proposition in appendix B. In the following, we will work with the $\ell^{2}$-norm on subsets $S$ of $\{0,\dots,L-1\}$. For such sets, we define

[TABLE]

We may now prove the following result on the stability of magnitude retrieval:

Lemma 3.4** (Stability of magnitude retrieval).**

Let $\delta>0$ and let $x,y\in\mathbb{C}^{L}$ be $B$-bandlimited, for some $B\in\mathbb{N}_{0}$. Define $G_{\delta}=(V_{\delta},E)$ to be the essential support graph of $x$ and let $\phi\in\mathbb{C}^{L}$ be such that

[TABLE]

for some $c>0$. Then,

[TABLE]

Proof.

Let us start with the simple fact that for $\ell\in V_{\delta}$ (cf. definition 3.1),

[TABLE]

Thus, we have

[TABLE]

Next, suppose that $4B<L$. Employing Plancherel’s theorem, the ambiguity function relation and proposition 3.3, we find

[TABLE]

The proof for $4B\geq L$ is even simpler: In this case $\mathcal{A}[\phi](0,\cdot)$ is lower bounded everywhere. ∎

Remark 3.5**.**

Note that we need to restrict our stability result to the essential support $V_{\delta}$ of $x$ because the square root $t\mapsto\sqrt{t}$ is not Lipschitz continuous. For this reason, we obtain the dependence of our stability result on $\delta$. We note that in the above proof, we derive

[TABLE]

Hence, magnitude retrieval is, in fact, stable even when we consider small entries as long as we compare the squared magnitudes of the signals with the squared absolute values of the Gabor transform.

Next, we turn to the retrieval of the phases. First, in accordance with corollary 2.4, we will only use the entries of $\mathcal{A}[x](1,\cdot)$ (and of $\mathcal{A}[x](0,\cdot)$) for our recovery which allows us to do phase propagation on adjacent entries. To be precise, we can propagate the phase from $x(\ell)$ to $x(\ell+1)$ (or back), for any $\ell\in\{0,\dots,L-1\}$. Mathematically this fact can be captured with the help of the essential support graph $G_{\delta}$ of $x$ with time-separation parameter zero, i.e. $\Delta=0$. In the following, we will call the connected components of $G_{\delta}$ temporal islands.

Theorem 3.6** (Stability of phase retrieval on a single temporal island).**

Let $\delta>0$ and let $x,y\in\mathbb{C}^{L}$ be $B$-bandlimited, for $B\in\mathbb{N}_{0}$. For $G_{\delta}=(V_{\delta},E)$ the essential support graph of $x$, assume that the subgraph induced by the vertex set $S_{\delta}=V_{\delta}\cap\operatorname{supp}y$ is connected. If $\phi\in\mathbb{C}^{L}$ is such that

[TABLE]

for some $c>0$, then

[TABLE]

Proof.

See section 5. ∎

Remark 3.7**.**

The stability constant derived in the above result is

[TABLE]

and consists of a contribution from the magnitude retrieval estimate in lemma 3.4 and the phase retrieval estimate presented in section 5. The contribution from the phase retrieval estimate contains a factor $\sqrt{2\lvert S_{\delta}\rvert}$ which can be interpreted as a mild ill-conditioning of phase retrieval as it might scale like $L^{\frac{1}{2}}$. Additionally, the phase retrieval estimate contains a factor $1/\delta$ which captures the dependence on the size of the magnitudes of $x$.

For a visualisation, we plot the magnitudes ambiguity functions of four commonly used window functions $\phi\in\mathbb{C}^{L}$ in figure 3. For reference, we use $L=1024$ and the windows

[TABLE]

Example 3.8**.**

We want to present an example to clarify the statement of Theorem 3.6. For this purpose, we let $L\in\mathbb{N}$, with $L\geq 6$, be arbitrary but fixed and consider the rectangular window of length two

[TABLE]

Note that the choice of the rectangular window of length two is rather arbitrary. One could, in fact, perform similar calculations with most other common window functions as long as one picks the window in such a way that its time support is small enough. We observe that, for large $L$, the rectangular window of length two will have a small time support (by which we mean that $\lvert\operatorname{supp}\phi\rvert=2$ is small compared to $L$) and a large frequency support (by which we mean that $\lvert\operatorname{supp}\mathcal{F}[\phi]\rvert$, which is readily seen to be $L$ or $L-1$ depending on whether $L$ is even or odd, is comparable to $L$). This property of the rectangular window of length two will carry over to its ambiguity function in the sense that even for large $B$, we find that

[TABLE]

is rather large and thereby the constant $c$ in Theorem 3.6 is rather small. So let $\widetilde{c}>1$ and $B\in\mathbb{N}_{0}$ be such that $\tfrac{L}{\widetilde{c}}\leq 6B\leq L$. Then, we find that

[TABLE]

Therefore, it follows that Theorem 3.6 holds with $c=\sqrt{L}$, i.e. for $\delta>0$ and $x,y\in\mathbb{C}^{L}$ $B$-bandlimited such that the subgraph of the essential support graph of $x$ induced by the vertex set $S_{\delta}=V_{\delta}\cap\operatorname{supp}y$ is connected, we have

[TABLE]

It follows that for the rectangular window of size two and for $6B\leq L$, the stability constant scales linearly in $L$. The dimension of the space of $B$-bandlimited signals is $d:=2B+1$ and therefore it follows from $\tfrac{L}{\widetilde{c}}\leq 6B$ that

[TABLE]

Example 3.9**.**

We also want to give an example which elaborates on the dependence of the stability constant on $\delta^{-1}$. On first glance, it might look like this dependence is only due to our analysis (the propagation of phases between adjacent entries.) However, this is not the case. Consider the rectangular window $\phi$ of length two along with the signals

[TABLE]

where $\delta\in(0,1)$. Then, we have

[TABLE]

In addition, we have

[TABLE]

Therefore, it follows that

[TABLE]

is a lower bound for the stability constant. Note that this example is independent of the reconstruction technique we have chosen.

3.2 Multiple islands and frequency gaps

The phase propagation procedure presented as part of the proof of theorem 3.6 carries over quite naturally to the case where the graph $G_{\delta}=(V_{\delta},E)$ is disconnected rather than connected. We say the graph $G_{\delta}$ has multiple temporal islands. It is of course interesting to consider this case, as there is a wide range of signals for which $G$ will be disconnected. For instance, recordings of human speech will typically consist of multiple temporal islands as speakers tend to leave short gaps (i.e. modes of silence) in between words. In addition, a discretisation of the signal $f_{\lambda}^{+}$ from the introduction (see figure 4) will yield two temporal islands.

Theorem 3.10** (Stability of phase retrieval on multiple temporal islands).**

Let $\delta>0$ and let $x,y\in\mathbb{C}^{L}$ be $B$-bandlimited, for $B\in\mathbb{N}_{0}$. For the essential support graph $G_{\delta}=(V_{\delta},E)$ of $x$, assume that the subgraph induced by the vertex set $S_{\delta}=V_{\delta}\cap\operatorname{supp}y$ has $K$ connected components whose vertex sets are denoted by $\{S_{k}\}_{k=1}^{K}$. If $\phi\in\mathbb{C}^{L}$ is such that

[TABLE]

for some $c>0$, then

[TABLE]

Proof.

See theorem 3.11. ∎

Furthermore, one should note that until now we have only worked with minimal restrictions on the ambiguity function $\mathcal{A}[\phi]$ of the window $\phi$, i.e. we have only utilised the ambiguity function for $m=0,1$. In the following, we want to generalise our result to be able to use $\mathcal{A}[\phi](m,n)$ for $m=0,\dots,\Delta+1$. In particular, we may be able to harness corollary 2.5 in order to propagate phase stably across a section of the signal in which the entries consistently fall below a threshold $\delta$. To precisely describe this phase propagation procedure, we make use of the essential support graph of signals $x$ with time-separation parameter $\Delta$.

Theorem 3.11** (Main theorem).**

Let $\Delta\in\mathbb{N}_{0}$, let $\delta>0$ and suppose that $x,y\in\mathbb{C}^{L}$ are $B$-bandlimited, for $B\in\mathbb{N}_{0}$. Let $G_{\delta}=(V_{\delta},E)$ be the essential support graph of $x$ with time-separation parameter $\Delta$ and assume that the subgraph induced by the vertex set $S_{\delta}=V_{\delta}\cap\operatorname{supp}y$ has $K$ connected components whose vertex sets are denoted by $\{S_{k}\}_{k=1}^{K}$. If $\phi\in\mathbb{C}^{L}$ is such that

[TABLE]

for some $c>0$, then

[TABLE]

Proof.

See section 5. ∎

Remark 3.12**.**

Alternatively, one can show

[TABLE]

under the assumptions laid out in theorem 3.11. We prefer the result above as it is more compact. Note that neither of these results is stronger or weaker than the other.

We note that one can obtain dual results to the above by considering $\mathcal{F}[x]$ instead of $x$. A straight-forward calculation (see proof of proposition 3.3) yields

[TABLE]

In this light, it is not surprising that we can derive stability results for recovering $\mathcal{F}[x]$ from the measurements $M_{\phi}[x]$ resembling the theorems derived above. Note that in this way, one can also show the following dual of the ambguity function relation:

[TABLE]

Theorem 3.13** (Stability for frequency gaps).**

Let $B,\Delta\in\mathbb{N}_{0}$, let $\delta>0$ and suppose that $x,y\in\mathbb{C}^{L}$ have their support contained in $\{-B,\dots,B\}\mod L$. Let $G_{\delta}=(V_{\delta},E)$ be the essential support graph of $\mathcal{F}[x]$ with time-separation parameter $\Delta$ and assume that the subgraph induced by the vertex set $S_{\delta}=V_{\delta}\cap\operatorname{supp}\mathcal{F}[y]$ has $K$ connected components whose vertex sets are denoted by $\{S_{k}\}_{k=1}^{K}$. If $\phi\in\mathbb{C}^{L}$ is such that

[TABLE]

for some $c>0$, then

[TABLE]

Remark 3.14**.**

In the preceding pages, we have presented approaches for phase retrieval for signals with multiple temporal or frequency islands. Unfortunately, it is not so clear how to extend this work to the more general case of time-frequency atolls considered in [1, 15]. It is likely that one has to come up with a different approach that allows one to do phase propagation in frequency and time direction simultaneously to actually handle time-frequency atolls.

We want to end this section by remarking that from our proof strategy for the frequency result a straight-forward dual version of corollary 2.5 follows.

Corollary 3.15**.**

Let $\Delta\in\mathbb{N}_{0}$ and let $x,y,\phi\in\mathbb{C}^{L}$ be such that $M_{\phi}[x]=M_{\phi}[y]$. Assume that the window $\phi$ satisfies

[TABLE]

Define $G=(V,E)$ as the graph with vertex set $V=\operatorname{supp}\mathcal{F}[x]$ and edge set $E\subset\{0,\dots,L-1\}\times\{0,\dots,L-1\}$ given by

[TABLE]

If $\{V_{k}\}_{k=1}^{K}$ are the vertex sets of the connected components of $G$, then for all $k\in\{1,\dots,K\}$ there exists an $\alpha_{k}\in\mathbb{R}$ such that

[TABLE]

4 Stability estimates based on the autocorrelation relation

The goal of this section is to apply the techniques we have developed thus far to the setup proposed in [19]. Our approach will be designed to work for bandlimited signals $x\in\mathbb{C}^{L}$ which potentially have very small entries. In doing so, we do not need to require that the Fourier transform of the absolute value squared of the window is nowhere-vanishing. We emphasise these particularities as the stability result developed in [19] (theorem 4.2 on p. 375) relies on window functions for which $\mathcal{F}[\lvert\phi\rvert^{2}]$ is nowhere-vanishing, and while one may apply it to signals $x\in\mathbb{C}^{L}$ with very small entries, the resulting stability constant will be ill-behaved as it depends inversely on $\min_{\ell\in\operatorname{supp}x}\lvert x(\ell)\rvert^{2}$.

Recall from section 2 that the authors of [19] consider the graph $G=(V,E)$ with vertex set $V:=\operatorname{supp}(x)$ and an edge between $\ell,k\in V$ if

[TABLE]

where $\ell_{\phi}+1$ denotes the support length of the window function. They do then propose to reconstruct the magnitude of a signal $x\in\mathbb{C}^{L}$ using the ambiguity function relation (2) (as in section 3) and propagate phase from one entry of $x$ to another if their indices have distance $\ell_{\phi}$ using equation (1). As before, we will work with local lower bounds on the ambiguity function of the window and the signal in order to ensure that all the aforedescribed steps can be carried out stably.

In order to introduce the local lower bounds on the signal, we will have to modify the graphs presented in [19] slightly: Let us consider a signal $x\in\mathbb{C}^{L}$, a tolerance parameter $\delta>0$ and a window function $\phi\in\mathbb{C}^{L}$ such that

[TABLE]

for $n_{0},\ell_{\phi}\in\{0,\dots,L-1\}$ such that $2\ell_{\phi}<L$. Introduce the graph $G_{\delta}=(V_{\delta},E)$ with vertex set

[TABLE]

and an edge between $\ell,k\in V$ if

[TABLE]

Under these assumptions, we can state the following stability estimate whose proof is inspired by the proof of corollary 2.8:

Theorem 4.1**.**

Let $\ell_{\phi},n_{0}\in\mathbb{N}_{0}$ such that $2\ell_{\phi}<L$. Furthermore, let $\delta>0$ and let $x,y\in\mathbb{C}^{L}$ be $B$-bandlimited, for $B\in\mathbb{N}_{0}$. Suppose that the subgraph induced by the vertex set $S_{\delta}=V_{\delta}\cap\operatorname{supp}y$ is connected and that the window $\phi$ satisfies

[TABLE]

as well as

[TABLE]

for some $c>0$. Then,

[TABLE]

Proof.

See section 5. ∎

Remark 4.2**.**

The stability constant in this result is

[TABLE]

The part of the constant due to magnitude retrieval is exactly the same as in the results in section 3. The part of the constant stemming from phase retrieval is a slight modification from the constants in section 3. It is mostly the term $\lvert\phi(n_{0})\phi(n_{0}+\ell_{\phi})\rvert$ in the denominator that deserves some attention. It is clear that phase propagation based on the relation

[TABLE]

will be unstable whenever the ends of the window $\phi(n_{0})$ and $\phi(n_{0}+\ell_{\phi})$ are close to zero. In particular, the reconstruction method proposed by the authors of [13, 19] benefits from windows for which $\lvert\phi(n_{0})\phi(n_{0}+\ell_{\phi})\rvert$ is large such as the Hamming or rectangular window.

As remarked before, in contrast to the stability result in [19], our result is applicable even when the Fourier transform of the magnitude squared of the window function $\mathcal{F}[\lvert\phi\rvert^{2}]$ has vanishing entries at the cost of only being applicable to bandlimited signals. We should also note that in [19] the stability constant for the phase retrieval estimate scales like $\sqrt{\lvert\operatorname{supp}x\rvert}\cdot L^{3}$ (which becomes $L^{7/2}$ for nowhere-vanishing signals) whereas our stability constant merely scales like $\sqrt{\lvert S_{\delta}\rvert L}$ (which becomes $L$ for signals whose entries have absolute values in excess of the threshold $\delta$).

Remark 4.3** (On disconnected graphs and duality results).**

Finally, we would like to note that one can prove a result resembling theorem 4.1 in the case where $S_{\delta}$ has $K\in\mathbb{N}$ connected components whose vertex sets are denoted by $S_{1},\dots,S_{K}\subset S_{\delta}$.

Similarly, we may utilise lemma 2.1 in order to deduce that

[TABLE]

holds, for $x,\phi\in\mathbb{C}^{L}$ and $m,n=0,\dots,L-1$. This in turn can be used to deduce a stability result which is essentially the Fourier-dual of theorem 4.1.

5 Proofs of the main results

Proof of corollary 2.5.

By the ambiguity function relation and our assumptions, we find that

[TABLE]

Therefore,

[TABLE]

and in particular $x$ and $y$ have the same magnitudes. Let us now consider $k\in\{1,\dots,K\}$ as well as some $\ell_{0}$ in $V_{k}$. We have $\lvert x(\ell_{0})\rvert=\lvert y(\ell_{0})\rvert$ and hence $x(\ell_{0})=\mathrm{e}^{\mathrm{i}\alpha_{k}}y(\ell_{0})$, for some $\alpha_{k}\in\mathbb{R}$. As $V_{k}$ is the vertex set of a connected component of $G$, it follows that for all $\ell\in V_{k}\setminus\{\ell_{0}\}$, there exists a (simple) path from $\ell_{0}$ to $\ell$. Therefore, we can consider $\ell\in V_{k}\setminus\{\ell_{0}\}$ and let $(u_{0},\dots,u_{n})$ be the vertex sequence of the path from $u_{0}=\ell_{0}$ to $u_{n}=\ell$. For $j\in\{0,\dots,n-1\}$ one has, by definition of the edge set, that

[TABLE]

Thus, there exists an $m_{j}\in\{1,\dots,\Delta+1\}$ such that $u_{j+1}-u_{j}=m_{j}\mod L$ or $u_{j}-u_{j+1}=m_{j}\mod L$. In either case, it follows from equation (5) that $x(u_{j+1})\overline{x(u_{j})}=y(u_{j+1})\overline{y(u_{j})}$. By induction on $j$, we find that $x(\ell)=\mathrm{e}^{\mathrm{i}\alpha_{k}}y(\ell)$. ∎

Proof of theorem 3.6.

The case $4B\leq L$ is similar to the case $4B>L$ but simpler: Indeed consider $4B\leq L$. In this case, we have

[TABLE]

by assumption. Therefore, we can replace all sums over $\{-2B,\dots,2B\}$ by sums over $\{0,\dots,L-1\}$ in this proof. So let us consider $4B>L$: Let $\alpha\in\mathbb{R}$ be arbitrary. Employing proposition A.1, we have

[TABLE]

The magnitude difference is estimated as in lemma 3.4. For the estimate of the phase difference, we develop inequalities in the following. Let $\ell,k\in S_{\delta}$. According to proposition A.2, we have

[TABLE]

Using the above inequality recursively, one obtains that for all $M\in\mathbb{N}$ and $u_{0},u_{1},\dots,u_{M}\in S_{\delta}$:

[TABLE]

Suppose now that $\ell_{0}$ is chosen such that any other vertex $\ell\in S_{\delta}$ has graph distance (in the induced subgraph) at most $\lvert S_{\delta}\rvert/2$ from $\ell_{0}$. Then, for any $\ell\in S_{\delta}\setminus\{\ell_{0}\}$, there exists $M(\ell)\in\mathbb{N}$, with $M(\ell)\leq\lvert S_{\delta}\rvert/2$, and a sequence $u_{0}^{\ell}=\ell_{0},u_{1}^{\ell},\dots,u_{M(\ell)}^{\ell}=\ell$ in $S_{\delta}$ such that (cf. definition 3.1 with $\Delta=0$)

[TABLE]

Therefore, there exists a sequence $\sigma_{1}^{\ell},\dots,\sigma_{M(\ell)}^{\ell}$ in $\{-1,1\}$ such that

[TABLE]

Now, let $\alpha\in\mathbb{R}$ be such that

[TABLE]

Then, we have for any $\ell\in S_{\delta}$, according to the above considerations (and inequality (7)), that

[TABLE]

For the second term of the right-hand side of inequality (6) this yields

[TABLE]

Applying Jensen’s inequality on the square of the inner sum and noting that $M(\ell)\leq\lvert S_{\delta}\rvert/2$, we obtain

[TABLE]

Since $u^{\ell}_{j}\in S_{\delta}$, for $j\in\{1,\dots,M(\ell)\}$, we can further estimate

[TABLE]

and with $\sigma^{\ell}_{j}\in\{\pm 1\}$, for $j\in\{1,\dots,M(\ell)\}$, we also get

[TABLE]

There are no repetitions in the sequences $u^{\ell}_{1},u^{\ell}_{2},\dots,u^{\ell}_{M(\ell)}$ and hence

[TABLE]

Therefore, we have

[TABLE]

Suppose now that $4B<L$. By Plancherel’s theorem, it holds that

[TABLE]

It follows from the ambiguity function relation and the lower bound on the ambiguity function of the window on $\{0,1\}\times\{-2B,\dots,2B\}$ that

[TABLE]

where we have used Plancherel’s theorem in the last equality. ∎

Proof of theorem 3.11.

Let $k\in\{1,\dots,K\}$ and $\alpha_{k}\in\mathbb{R}$. As in the proof of theorem 3.6, we start by splitting the estimate into a phase and a magnitude estimate using proposition A.1:

[TABLE]

As the connected components $\{S_{k}\}_{k=1}^{K}$ are disjoint subsets of $V_{\delta}$, we can use Jensen’s inequality to see that

[TABLE]

Employing lemma 3.4, we obtain for the magnitude retrieval estimate

[TABLE]

The phase difference is estimated just like in theorem 3.6: First, observe that there must exist a vertex $\ell_{0}\in S_{k}$ such that any other vertex $\ell\in S_{k}$ has graph distance $M(\ell)$ at most $\tfrac{L+\Delta}{2+\Delta}$ from $\ell_{0}$. Indeed, consider the following argument: The worst case which could happen is that we need to connect the vertex [math] to the vertex $\lfloor L/2\rfloor$. By definition of the graph, it will take us exactly one step to go from [math] to any $\ell\in\{1,\dots,\Delta+1\}\cap S_{k}$; it will take us exactly two steps to go from [math] to $\Delta+2$, if the latter is in $S_{k}$; and it will take us at most three steps to go from [math] to $\ell\in\{\Delta+3,2\Delta+3\}\cap S_{k}$. Following this logic, it is not too hard to see that it will take us at most $2n$ steps to go from [math] to $n(\Delta+2)$, if the latter is in $S_{k}$, and it will take us at most $2n+1$ steps to go from [math] to $\ell\in\{n(\Delta+2)+1,\dots,(n+1)(\Delta+2)-1\}\cap S_{k}$. So, if there exists an element $n\in\mathbb{N}$ such that $n(\Delta+2)=\lfloor L/2\rfloor$, then it will take us at most

[TABLE]

steps to connect [math] and $\lfloor L/2\rfloor$. Similarly, if there is an element $n\in\mathbb{N}_{0}$ such that $\lfloor L/2\rfloor\in\{n(\Delta+2)+1,\dots,(n+1)(\Delta+2)-1\}$, then there is a $\beta\in\{1,\dots,\Delta+1\}$ such that

[TABLE]

So for any $\ell\in S_{k}\setminus\{\ell_{0}\}$, there exists a path $u_{0}^{\ell}=\ell_{0},u_{1}^{\ell},\dots,u_{M(\ell)}^{\ell}=\ell$ from $\ell_{0}$ to $\ell$. By definition, this path satisfies

[TABLE]

Therefore, there exist sequences $\sigma_{1}^{\ell},\dots,\sigma_{M(\ell)}^{\ell}\in\{\pm 1\}$ and $\Delta_{1}^{\ell},\dots,\Delta_{M(\ell)}^{\ell}\in\{1,\dots,\Delta+1\}$ such that

[TABLE]

We let $\alpha_{k}\in\mathbb{R}$ be chosen in such a way that

[TABLE]

Proceeding as in the proof of theorem 3.6 (now with $M(\ell)\leq(L+\Delta)/(2+\Delta)$), we derive

[TABLE]

We first treat the case $2\Delta<L-2$ and use that $\sigma_{j}^{\ell}\Delta_{j}^{\ell}\in\{-\Delta-1,\dots,\Delta+1\}$ to estimate

[TABLE]

We may use that there are no repetitions in the sequences $u_{1}^{\ell},\dots,u_{M(\ell)}^{\ell}$ to obtain

[TABLE]

Suppose furthermore that $4B>L$. According to Plancherel’s theorem, we find that

[TABLE]

Next, we use the ambiguity function relation, inequality (4) as well as the symmetry of the ambiguity function of the window around the origin to derive

[TABLE]

Note that we may once again use Jensen’s inequality to show that

[TABLE]

Thus, combining the phase and the magnitude estimates yields

[TABLE]

The cases in which $2\Delta\geq L-2$ or $4B\geq L$ are dealt with similarly. ∎

Proof of theorem 4.1.

Let $\alpha\in\mathbb{R}$. As in the proof of theorem 3.6, we start by splitting the estimate into a phase and a magnitude estimate

[TABLE]

First, noting that

[TABLE]

we can apply lemma 3.4 to obtain the estimate

[TABLE]

For the estimate on the phase retrieval part, we need to consider a new strategy based on equation (1). First, we find that there must exist a vertex $\ell_{0}\in S_{\delta}$ such that any other vertex $\ell\in S_{\delta}$ has graph distance $M(\ell)$ at most $\lvert S_{\delta}\rvert/2$ from $\ell_{0}$. So for any $\ell\in S_{\delta}\setminus\{\ell_{0}\}$, there exists a path $u_{0}^{\ell}=\ell_{0},u_{1}^{\ell},\dots,u_{M(\ell)}^{\ell}=\ell$ from $\ell_{0}$ to $\ell$. By definition, this path satisfies

[TABLE]

Therefore, there exists a sequence $\sigma_{1}^{\ell},\dots,\sigma_{M(\ell)}^{\ell}\in\{\pm 1\}$ such that

[TABLE]

With this at hand, we proceed similarly to the proof of theorem 3.6. We let $\alpha\in\mathbb{R}$ be chosen in such a way that

[TABLE]

Then, we have that for any $\ell\in S_{\delta}$

[TABLE]

and

[TABLE]

by Jensen’s inequality, $M(\ell)\leq\lvert S_{\delta}\rvert/2$ and the fact that $u_{j}^{\ell},u_{j}^{\ell}-\sigma_{j}^{\ell}\ell_{\phi}\in S_{\delta}$, for all $j\in\{1,\dots,M(\ell)\}$. As before, we can further employ a crude estimate to derive

[TABLE]

since $\sigma_{j}^{\ell}\in\{\pm 1\}$, for all $j\in\{1,\dots,M(\ell)\}$ and all $\ell\in S_{\delta}$, and because for fixed $\ell\in S_{\delta}$, the $u_{j}^{\ell}$ are all distinct. Next, we note that due to $\operatorname{supp}\phi=\{n_{0},\dots,n_{0}+\ell_{\phi}\}\mod L$, $2\ell_{\phi}<L$, and the autocorrelation relation, we have for $u\in S_{\delta}$

[TABLE]

as well as

[TABLE]

This implies

[TABLE]

Yet another crude estimate results in

[TABLE]

∎

Appendix A Pointwise estimates for the stability results

Let us start by the typical splitting of signal differences into phase and magnitude part.

Proposition A.1**.**

Let $x,y\in\mathbb{C}^{L}$. Then, we have for all $\alpha\in\mathbb{R}$ that

[TABLE]

holds for all $\ell\in\{0,\dots,L-1\}$ for which $x(\ell)\neq 0$ and $y(\ell)\neq 0$. Moreover,

[TABLE]

holds for all $\ell\in\{0,\dots,L-1\}$ for which $x(\ell)=0$ or $y(\ell)=0$.

Proof.

Consider the case $x(\ell)\neq 0$ and $y(\ell)\neq 0$:

[TABLE]

If $x(\ell)=0$, then

[TABLE]

∎

In addition, a result about phase propagation will be handy.

Proposition A.2**.**

Let $x,y\in\mathbb{C}^{L}$ and let $\alpha\in\mathbb{R}$. Then, we have that

[TABLE]

holds, for all $\ell,k\in\{0,\dots,L-1\}$ for which $x(k),x(\ell),y(k)$ and $y(\ell)$ are all different from zero.

Proof.

We compute

[TABLE]

The claim follows by noting that using the triangle inequality and the reverse triangle inequality, one has

[TABLE]

for $z_{0},z_{1}\in\mathbb{C}$. ∎

Appendix B Proofs of the most fundamental formulae

For convenience of the reader, we present the proofs of the two formulae presented in section 2. We note that these formulae are well-known in the literature and have repeatedly been used to prove uniqueness results in recent years [6, 13, 19, 21]. We start with lemma 2.1.

Proof of lemma 2.1.

For $m,n\in\{0,\dots,L-1\}$, one can calculate

[TABLE]

using the convolution theorem for the DFT. Applying the inverse Fourier transform yields the statement. ∎

Proof of lemma 2.2.

According to lemma 2.1, we have

[TABLE]

for $m,n\in\{0,\dots,L-1\}$, and therefore

[TABLE]

Taking the DFT in $m$ yields

[TABLE]

∎

Proof of proposition 3.3.

Let us use the notation $x^{\prime}_{m}(\ell)=\overline{x(\ell-m)}$, for $\ell,m\in\{0,\dots,L-1\}$, and denote the entry-wise product of two vectors via $\odot$. For $m,n\in\{0,\dots,L-1\}$ we apply the convolution theorem for the DFT to see that

[TABLE]

Next, we compute

[TABLE]

for $k\in\{0,\dots,L-1\}$. Therefore,

[TABLE]

Let us assume that $4B<L$ and consider $k\in\{-B,\dots,B\}$. First, suppose that $n\in(-\tfrac{L}{2},-2B)\cap\mathbb{Z}$. Then, it follows that

[TABLE]

and since $L-B\geq\tfrac{L}{2}+B$, it follows that $\mathcal{F}[x](k-n)=0$. Therefore,

[TABLE]

Secondly, suppose that $n\in(2B,\tfrac{L}{2}]\cap\mathbb{Z}$. Then, one has

[TABLE]

and $\tfrac{L}{2}-B>B$ implies $\mathcal{F}[x](k-n)=0$. Thus,

[TABLE]

What remains is the case $4B\geq L$: Note that by $\{-2B,\dots,2B\}=\{0,\dots,L-1\}\mod L$ the proposition is trivial in this case. ∎