Next Article in Journal
On the Roots of a Family of Polynomials
Next Article in Special Issue
Fast Linear Canonical Transform for Nonequispaced Data
Previous Article in Journal
On Coupled System of Langevin Fractional Problems with Different Orders of μ-Caputo Fractional Derivatives
Previous Article in Special Issue
2D Linear Canonical Transforms on Lp and Applications
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Brief Report

Unlimited Sampling Theorem Based on Fractional Fourier Transform

1
School of Mathematics and Statistics, Beijing Institute of Technology, Beijing 100081, China
2
Beijing Key Laboratory on MCAACI, Beijing Institute of Technology, Beijing 100081, China
*
Author to whom correspondence should be addressed.
Fractal Fract. 2023, 7(4), 338; https://doi.org/10.3390/fractalfract7040338
Submission received: 18 February 2023 / Revised: 13 April 2023 / Accepted: 15 April 2023 / Published: 18 April 2023
(This article belongs to the Special Issue Recent Advances in Fractional Fourier Transforms and Applications)

Abstract

:
The recovery of bandlimited signals with high dynamic range is a hot issue in sampling research. The unlimited sampling theory expands the recordable range of traditional analog-to-digital converters (ADCs) arbitrarily, and the signal is folded back into a low dynamic range measurement, avoiding the saturation problem. Since the non-bandlimited signal in the Fourier domain cannot be directly applied to its existing theory, the non-bandlimited signal in the Fourier domain may be bandlimited in the fractional Fourier domain. Therefore, this brief report studies the unlimited sampling problem of high dynamic non-bandlimited signals in the Fourier domain based on the fractional Fourier transform. Firstly, a mathematical signal model for unlimited sampling is proposed. Secondly, based on this mathematical model, the annihilation filtering method is used to estimate the arbitrary folding time. Finally, a novel fractional Fourier domain unlimited sampling theorem is obtained. The theory proves that, based on the folding characteristics of the self-reset ADC, the number of samples is not affected by the modulo threshold, and any folding time can be handled.

1. Introduction

In signal processing, sampling [1,2,3] is the primary task faced in the process of digitizing the signal. Since Shannon’s sampling theorem was proposed [4], sampling theory has been developed for more than 70 years, and its theoretical results [5,6,7,8] are so rich that it has become one of the research hotspots in the field of signal processing. From a practical standpoint, point-wise samples of the function are obtained using the analog-to-digital converter (ADC), but the ADC has a limited dynamic range [ λ , λ ] [9]. Whenever the signal exceeds a certain preset threshold λ , the ADC will saturate, and the aliases signal will be clipped due to clipping, we refer to these ADCs as clipped ADCs or traditional ADCs (C-ADC) [10,11]. Since most signals in practical applications are not limited by broadband, the dynamic range is very wide. Therefore, recovering the signal from the tailored version of the C-ADC is an inaccurate inverse problem. Generally, recovery methods almost alleviate the effects of clipping at the expense of oversampling. To solve these problems, so self-reset ADC (S-ADC) was proposed [12,13,14]. Each time the input signal reaches the upper (lower) saturation limit, these S-ADCs will be reset to the other corresponding thresholds, which allows the S-ADC to reset rather than saturate, resulting in analog sampling. When the signal reaches the upper (lower) threshold point, it will fold backward (front) by an integer multiple of 2 λ . This phenomenon is equivalent to modulo arithmetic on the input signal, which is very helpful for processing high dynamic range signals.
Because of the S-ADC’s ability to process high dynamic range signals, Bandari et al. recently made the first pioneering contribution [15]. He proposed the unlimited sampling theorem and developed the first provable reflector of the guaranteed algorithm. Similar to Shannon’s sampling theorem, the unlimited sampling theorem proves that the bandlimited signal can be recovered from analog sampling as long as it meets a certain sampling density criterion and is not affected by the ADC threshold. In this way, the results allow the perfect recovery of the bandlimited function, whose amplitude exceeds the ADC threshold by orders of magnitude. The results [15] have led to a lot of follow-up work, and the theoretical research of unlimited sampling has gradually enriched [16,17,18,19,20]. In this modular sampling framework, many scholars have studied the sampling and reconstruction of bandlimited functions and smooth functions under different backgrounds. The paper is shown that the bandlimited function is uniquely characterized by modular samples under certain conditions [21]. Ordentlich et al. studied the recovery of quantization modulus samples by using edge information [22]. Musa et al. gave the modulus sampling theory of the S-ADC sparse signal [23]. The unlimited sampling method based on wavelet is suitable for general smooth signals, not limited to bandlimited signals [24]. The paper is mainly applicable to bandlimited signals in the Fourier domain on unlimited sampling method [25]. Most of the unlimited sampling frameworks are based on the bandlimited signals, but there are few articles on non-bandlimited signals. For various applications of the non-bandlimited signal models, the original results are not directly applicable. Therefore, it is very necessary to study the unlimited sampling theory of non-bandlimited signals in the Fourier domain.
As a general form of the Fourier transform (FT), the fractional Fourier transform (FRFT) extends dimensionality of traditional FT-based spectral analysis [6,26,27]. The FRFT has an additional degree of freedom compared to the FT, which makes it more flexible and suitable for non-stationary signals. By transforming α from 0 to π 2 , the FRFT will be able to fully characterize the signal at the transition from time to frequency, overcoming the limitation of the FT only to perform frequency analysis. Furthermore, unlike conventional FT, which uses a complex exponential signal as the basis function, the FRFT utilizes a chirp signal as the basis function, which implies that a signal that is not bandlimited in the Fourier domain can be bandlimited in the fractional Fourier domain. Therefore, it is very meaningful to study the sampling theorem under the unlimited sampling framework based on the FRFT.
The FRFT can expand the signal range applicable to traditional sampling theory. Since the non-bandlimited signal in the Fourier domain cannot be directly applied to its existing theory, the non-bandlimited signal in the Fourier domain may be bandlimited in the fractional Fourier domain. Therefore, we propose a method based on the fractional Fourier domain bandlimited signal reconstruction theory under the framework of unlimited sampling to make up for the shortcomings of the Fourier domain reconstruction theory. To enrich the content of the unlimited sampling framework, we extend the theory of the unlimited sampling framework in the Fourier domain to the fractional Fourier domain. Our main work is to perform a modular operation on folds introduced by modular nonlinearity in the fractional Fourier domain and to estimate the fractional spectrum of instantaneous folding time to obtain a new sampling theorem in the fractional Fourier domain.
The main contributions of this brief report are summarized as follows. Firstly, the folding formula of non-linear and modulus mapping in the fractional Fourier domain is introduced, and a mathematical signal model of unlimited sampling for the FRFT is proposed. Secondly, based on this mathematical model, the fractional spectrum of the arbitrary folding time is estimated using the annihilation filtering method. Finally, a novel fractional Fourier domain unlimited sampling theorem is obtained. It is shown that the sampling theorem is independent of the modulo threshold and that it can be applied to arbitrary folding times.
This report is organized as follows. In Section 2, we briefly describe the FRFT and unlimited sampling theorem in the Fourier domain. In Section 3, a mathematical signal model for unlimited sampling is proposed in the fractional Fourier domain. In Section 4, a novel unlimited sampling theorem in the fractional Fourier domain is proposed. In Section 5, relevant applications of the proposed method are provided. We conclude this report in Section 6.

2. Preliminaries

2.1. Fractional Fourier Transform

Definition 1. 
The FRFT of a signal x ( t ) L 2 ( R ) with an angle α is defined as [27,28,29,30]
X α ( u ) = F α x ( t ) ( u ) + x ( t ) K α ( u , t ) d t ,
 where F α is the FRFT operator, u stands for fractional frequency, K α ( u , t ) denotes the kernel function of the FRFT
K α ( u , t ) = A α e i ( cot   α 2 t 2 csc   α u t + cot   α 2 u 2 ) , α k π δ ( t u ) , α = 2 k π δ ( t + u ) , α = ( 2 k 1 ) π
where A α 1 i   cot   α 2 π , the rotation angle of the FRFT is expressed as α = p π 2 , and p is the order of the FRFT. The domain 0 < α < π 2 are called fractional Fourier domains in [28], and this definition is also adopted in this report.
The FRFT can be understood as the rotation of the time-frequency plane. The essence of the FRFT of a signal is to decompose the signal with the chirp signal K α ( u , t ) as the basis function. According to the FRFT of the signal x ( t ) , it can be determined whether it is bandlimited in the fractional Fourier domain.
The FRFT has linear transform additivity, namely
F α + β x ( t ) ( u ) = F α x ( t ) ( u ) · F β x ( t ) ( u ) = X α ( u ) · X β ( u ) .
It can be seen that the inverse transform of the FRFT relative to the α angle is the FRFT with the parameter α angle, we have
x ( t ) = F α X α ( u ) = + X α ( u ) K α ( u , t ) d u ,
when α = π 2 , the FRFT degenerates to the traditional inverse FT; when α = π 2 , the FRFT degenerates to traditional FT, X π 2 ( u ) = + x ( t ) e i 2 π u t d t ; when α = 0 , the FRFT degenerates to an identity transformation, X 0 ( u ) = x ( t ) ; when α = π , the FRFT degenerates to the inversion of the signal with respect to the time axis, X π ( u ) = x ( t ) .
Definition 2. 
A signal x ( t ) is called Ω α bandlimited signal in the fractional Fourier domain, which means
X α ( u ) = 0 , u > Ω α ,
where Ω α is called the bandwidth of signal x ( t ) in the fractional Fourier domain. It has been shown that if a nonzero signal is bandlimited in the αth fractional Fourier domain, it can’t be bandlimited in the fractional Fourier domain with another angle β, where β ± α + n π for any integer n [26].

2.2. Unlimited Sampling Theorem in the Fourier Domain

Definition 3. 
The central modulo operation is defined by the mapping [15]
M λ : g 2 λ g 2 λ + 1 2 1 2 , g = d e f g g ,
where g defines the fractional part of input signal g and λ > 0 is the ADC threshold. Note that Equation (6) is a nonlinear modulus mapping, which converts a smooth function into a discontinuous function. It is equivalent to a centered modulo operation since M λ ( g ) g mod 2 λ . By implementing the mapping Equation (6), it is clear that out-of-range amplitudes are folded back into the dynamic range [ λ , λ ] .
Let’s review some important conclusions in [16,25].
Lemma 1 
((Modular decomposition property) [16]). Let g B Ω , where B Ω denotes the space of σ-bandlimited functions, and M λ ( · ) be defined in Equation (6) with λ is a fixed, positive constant. Then, the bandlimited function g ( t ) admits a decomposition
g ( t ) = z ( t ) + ε g ( t ) ,
ε g ( t ) = 2 λ m Z e m 1 D m ( t ) , e m Z ,
where z ( t ) = M λ g ( t ) , ε g ( t ) is a simple function, and m Z D m = R is a partition of the real line into intervals D m .
The process of solving discontinuities is very critical. Lemma 1 just proves this problem. Each bandlimited function, whether continuous or discrete, can be decomposed into the sum of the modular function and the stepwise residual of the simple function. Observe that the output function z ( t ) is the difference between g ( t ) and a piecewise constant function ε g ( t ) .
Theorem 1 
((Unlimited sampling theorem in the Fourier domain) [25]). Let g B Ω be a τ-periodic function. Suppose that we are given Q modulo samples of where y [ k ] = M λ ( g ( k T ) ) folded at most M times. Then a sufficient condition for recovery of g ( t ) from y [ k ] (up to a constant) is that, T τ Q and Q 2 Ω σ 2 π + M + 1 .
Unlike the traditional FT, which uses a complex exponential signal as the basis function, the FRFT uses a chirp signal as the basis function. This connotation determines that a non-bandlimited signal in the Fourier domain may be bandlimited in the fractional Fourier domain. Therefore, the unlimited sampling study of non-bandlimited signals in the Fourier domain can be transformed into the theoretical study of bandlimited signals in the fractional Fourier domain. The next step is to study the unlimited sampling theory of bandlimited signals in the fractional Fourier domain.

3. Mathematical Model for Unlimited Sampling with FRFT

In this Section, we will study the unlimited sampling theory of bandlimited signals in the fractional Fourier domain. Here we make the following symbolic regulations: the sets of real, integer, and complex-valued numbers are denoted by R , Z , and C , respectively.

3.1. Mathematical Signal Model

Based on the periodic signal model in the Fourier domain, this report proposes a periodic signal model in the fractional Fourier domain, so mathematically x ( t ) can be represented as follows.
Let x ( t ) be a Ω α bandlimited function satisfies x ( t ) = x ( t + σ ) , t R . Then x ( t ) has a fractional Fourier series (FRFS) expansion
x ( t ) = | w | R X ^ α ( w ) Φ α ( w , t ) ,
where X ^ α ( w ) is FRFS coefficient, and
Φ α ( w , t ) = sin α i cos α σ e i ( cot   α 2 t 2 csc ( α ) w u 0 t + cot   α 2 w 2 u 0 2 ) ,
where R = Ω α u 0 , u 0 = 2 π sin α σ .
The discrete-time representation x ( n T s ) , n Z of the signal x ( t ) can be obtained by uniformly sampling at intervals of T s . The discrete-time FRFT of the α angle of the discrete-time signal x ( t ) is defined as follows
X α , s = F α x ( n T s ) ( u ) n = + x ( n T s ) K α ( u , n T s ) ,
where K α is given by Equation (2), and T s is the sampling period.
Let sampling function x ( t ) obtains Q modular samples at the sampling rate T in the interval t σ 2 , σ 2 , then FRFS coefficients of x ( t ) in the fractional Fourier domain have a form
X ^ α ( w ) = σ 2 σ 2 x ( t ) Φ α ( w , t ) d t , w E R , Q , 0 , w I Q \ E R , Q , | w | > Ω α ,
where I Q = 0 , 1 , , Q 1 denote the set of Q contiguous integers, and E R , Q is given by
E R , Q = 0 , R Q R , Q 1 , | E R , Q = 2 R + 1 .
Remark 1. 
The well-known Fourier series (FS) is just a special case of the FRFS for α = π 2 , please see [25]. In order to solve for X ^ α ( w ) in (12), we must require Q 2 R + 1 . Because of Q T = σ , so T σ Q σ 2 R + 1 .
The hypothesis of periodic functions in our report only provides a practical method for recovering signals from the folding measurements below. However, when the signal is aperiodic, the theoretical reconstruction guarantees that the aperiodic signal can also be expanded by the discrete-time FRFT, but additional requirements are required for sampling samples, and this report will not expand in detail.

3.2. Nonlinear Modulus Mapping

This report uses the definition and properties of generalized modular non-linear mapping in Equation (6); this phenomenon is equivalent to modulo arithmetic on the input function. According to Lemma 1, it gives the following form
v x ( t ) = x ( t ) M λ x ( t ) = : m Z c m 1 [ t m , t m + 1 ] t ,
where c m R , 1 [ t a , t b ] is the indicator function on [ t a , t b ] , and t m σ 2 , σ 2 denotes the folding instants with t a < t b . Obviously, the output function M λ x ( t ) is the difference between x ( t ) and a residual function v x ( t ) . The sampling process to obtain modulo samples of a function is outlined in Figure 1. The [16] requires that the correlation coefficient of the residual function v x ( t ) is an integer multiple of 2 λ , while [25] does not need to make assumptions about its correlation coefficient.
Without loss of the generality, we make the following symbolic regulations:
  • Let f k = def x k T , h k = def M λ x k T , v k = def v x k T , then
    f k = h k + v k .
  • Let Δ N f = Δ N 1 Δ f denote the Nth difference operator with Δ f = f [ k + 1 ] f [ k ] , f ¯ [ k ] = def Δ f [ k ] , h ¯ [ k ] = def Δ h [ k ] , and v ¯ [ k ] = def Δ v x ( k T ) , then
    f ¯ [ k ] = h ¯ [ k ] + v ¯ [ k ] = h ¯ [ k ] + m M c m δ ( k T t m ) , k I Q ,
    where δ denotes the Dirac distribution, c [ m ] are unknown weights, t m are unknown fold instant, and the size of the set M depends on the dynamic range of the signal relative to the threshold λ .
Similar to the phase unwrapping theory, we can obtain the following fact from Itoh’s condition [31]. When the max-norm of the first-order finite difference of the samples is bounded by 2 λ or | f [ k + 1 ] f [ k ] | 2 λ , the first-order finite difference operator on the modular sequence can be reversed operation to restore.
Equation (16) is written as the fractional Fourier domain
H ¯ α [ n ] = F ¯ α [ n ] V ¯ α [ n ] , n E R , Q 1 V ¯ α [ n ] , n I Q 1 \ E R , Q 1
where F ¯ α , H ¯ α , and V ¯ α are the FRFT of f ¯ , h ¯ , and v ¯ , respectively. At the same time, the discrete FRFT form of h ¯ [ k ] is given
H ¯ α [ n ] = k I Q 1 h ¯ [ k ] K α ( n , k ) = k I Q 1 A α h ¯ [ k ] e i cot   α 2 k 2 csc ( α ) u ¯ 0 k n + cot   α 2 u ¯ 0 2 n 2 ,
where A α 1 i   cot   α 2 π and u ¯ 0 = 2 π sin α Q 1 . When α = π 2 , this transform is the discrete FT, see [25] for details.
If we want to recover f [ k ] , we must solve v [ k ] , then Equation (15) is transformed into solving Equation (16), and the key to solving Equation (17) is to find the value of the unknown folding instant c [ m ] , t m m Z . Using Equations (10) and (16), we can obtain
V ¯ α [ n ] = k I Q 1 m M c [ m ] δ k T t m Φ α ( n , k T ) = m M c [ m ] A α e i cot   α 2 T 2 t m 2 csc ( α ) u ¯ 0 n T t m + cot   α 2 u ¯ 0 2 n 2 ,
where A α 1 i   cot   α 2 π , and M denotes the size of the set depends on the dynamic range of the signal relative to the threshold λ . When α = π 2 , this transform is the discrete FT [25]. The estimation of the unknown parameters in Equation (19) is called the spectral estimation problem [32,33].

4. Unlimited Sampling Theorem in the Fractional Fourier Domain

4.1. Computing the Folding Instants

If we want to recover v [ k ] , we must find the value of the unknown folding instant c [ m ] , t m m Z . Equation (19) is the spectral estimation problem. The commonly used spectrum estimation methods are annihilation filter (AF) [32,34], ESPRIT [35], MUSIC [36], etc. Among them, the AF is the most commonly used method in many theoretical analyses and practical applications. In principle, the signal reconstruction process is to use the obtained set of moments or fractional Fourier coefficients of the input signal to solve a spectrum problem to achieve an accurate estimation of the unknown parameter c [ m ] , t m m Z . For convenience, Equation (19) is written as follows
V ¯ α [ n ] = A α e i cot   α 2 u ¯ 0 2 n 2 κ ( n ) m M c [ m ] e i cot   α 2 T 2 t m 2 χ m · e i   csc ( α ) u ¯ 0 T n t m ς m n ( n ) ,
where ( n ) = m M χ m ς m n . Since the formula is very complicated, we rewrite the formula V ¯ α [ n ] = A α κ ( n ) ( n ) . Because the part of κ ( n ) does not contain unknown parameters, we separately perform an annihilation filter to get c [ m ] , t m , and finally bring in Equation (20), and get v ¯ through inverse FRFT.
First, let’s analyze ( n ) in detail below,
( n ) = m M χ m ς m n .
Equation (21) is a classic spectrum estimation problem, which can be handled by an annihilation filter. It is known from the [32] that we can accurately estimate the unknown parameters χ m and ς m from 2 K continuous non-zero measured values ( n ) . The following is divided into two parts to solve separately.
(1) Construct the filter Γ [ ϑ ] ϑ = 0 , 1 , , M so that its zero point is the parameter
ς m = e i csc ( α ) u ¯ 0 T t ϑ ϑ = 0 ϑ = M 1 ,
then the z transform of the filter can be expressed as
Γ [ z ] = m = 0 M 1 ( 1 ς m z 1 ) = ϑ = 0 M Γ [ ϑ ] z ϑ .
It can be seen that the root of the polynomial is the parameter ς m . Therefore, this report convolutes it directly, so it has
Γ * [ n ] = ϑ = 0 M Γ [ ϑ ] [ n ϑ ] = ϑ = 0 M m = 0 M 1 c [ m ] e i cot   α 2 T 2 t m 2 · Γ [ ϑ ] · e i csc ( α ) · ( n ϑ ) u ¯ 0 T t m = m = 0 M 1 c [ m ] e i cot   α 2 T 2 t m 2 ϑ = 0 M Γ [ ϑ ] e i csc ( α ) · ϑ u ¯ 0 T t m Γ [ ς m ] e i csc ( α ) · n u ¯ 0 T t m = 0 .
We write Equation (23) in the form of matrix-vector to obtain
[ M 1 ] [ M 2 ] [ 0 ] [ M ] [ M 1 ] [ 1 ] [ N 1 ] [ N 2 ] [ N M ] Γ [ 1 ] Γ [ 2 ] Γ [ M ] = [ M ] [ M + 1 ] [ N ] ,
where = [ 0 ] , [ 1 ] , , [ M ] T , [ M ] = 1 . The unique solution can be obtained Γ [ ϑ ] , ϑ = 1 , 2 , , M , and the instantaneous folding time t m m Z can be obtained.
(2) To estimate the amplitude parameter χ m , extract M continuous values from the known coefficient [ n ] , that is to say, m = 0 , 1 , , M , and write the ( n ) = m = 0 M 1 χ m ς m n in the form of a matrix-vector
U χ = ,
1 1 1 ς 0 ς 1 ς M 1 ς 0 M 1 ς 1 M 1 ς M 1 M 1 e i cot   α 2 T 2 t 0 2 e i cot   α 2 T 2 t 1 2 e i cot   α 2 T 2 t M 1 2 · c [ 0 ] c [ 1 ] c [ M 1 ] = [ 0 ] [ 1 ] [ M 1 ] ,
where U is a vandermonde matrix, it is a matrix whose columns are geometric series. For any integer a , b = 0 , 1 , , M 1 ( a b ) satisfy U a U b , and U is non-singular, at this time, Equation (25) has a unique solution. It needs to be emphasized here that we generally use the least squares method to obtain an estimate of the amplitude information.
Through the above two steps, the instantaneous folding time c [ m ] , t m can be obtained. we bring c [ m ] , t m into Equation (20), and get v ¯ through the inverse FRFT.

4.2. Unlimited Sampling Theorem in the Fractional Fourier Domain

Through the above research content, we have found the value of the unknown folding instant c [ m ] , t m m Z . If v [ k ] is known, we can infer v [ k ] from h [ k ] and recover f [ k ] from h [ k ] . The unlimited sampling theorem in the fractional Fourier domain is given below.
Theorem 2 
(Unlimited sampling theorem in the fractional Fourier domain). Let x ( t ) be a Ω α bandlimited function satisfies x ( t ) = x ( t + σ ) , t R , and h k = M λ x k T folded at most M times. Then a sufficient condition for recovery of x ( t ) from h [ k ] , is that T σ Q and Q 2 Ω α σ 2 π sin α + M + 1 , where M is known.
Proof of Theorem 2. 
From the unlimited sampling part of the fractional Fourier domain in Equations (17) and (20), we can get the
[ M 1 ] [ M 2 ] [ 0 ] [ M ] [ M 1 ] [ 1 ] [ N 1 ] [ N 2 ] [ N M ] Γ [ 1 ] Γ [ 2 ] Γ [ M ] = [ M ] [ M + 1 ] [ N ] ,
that is, Γ [ ϑ ] .
From
Γ * [ n ] = 0 ,
we can get the root of ς m .
Bring Equation (19) into Equation (17), we obtain
F ¯ α [ n ] = H ¯ α [ n ] + V ¯ α [ n ] = k I Q 1 A α h ¯ [ k ] e i cot   α 2 k 2 csc ( α ) · u ¯ 0 k n + cot   α 2 u ¯ 0 2 n 2 + m M c [ m ] A α e i cot   α 2 T 2 t m 2 csc ( α ) · u ¯ 0 n T t m + cot   α 2 u ¯ 0 2 n 2 .
Using the least square method to estimate c [ m ] , and get v k in Equation (15). we develop a method that allows for inferring v [ k ] from h [ k ] .
Based on Theorem 1, we get the sampling density criterion, and the following conclusions can be drawn
| I Q 1 \ E R , Q 1 | = Q 2 R 2 2 M ,
where M is known.
Due Q T = σ , R = Ω α u 0 , u 0 = 2 π sin α σ , we have
T = T F R F T σ 2 ( R + M + 1 ) = σ 2 ( Ω α σ / 2 π sin α + M + 1 ) .
Proof completed. ☐
Remark 2. 
After performing M folds, Equation (29) can guarantee the restoration and reconstruction of folding moment c [ m ] , t m m = 0 M 1 . Theorem 2 turns out that the unlimited sampling theorem has nothing to do with the modulus threshold and can handle arbitrary folding time. When α = π 2 , see [25].
Remark 3. 
Indeed, the choice of the optimal angle α of the FRFT must take into account the actual application: if the modulation frequency is known, the optimal estimate can be obtained under the angle α determined by the modulation frequency; if the modulation frequency is unknown, it can be estimated in advance using the instants frequency estimation method, and then the estimated modulation frequency is used to determine the angle α.
Remark 4. 
The proposed method only needs to compute the first-order difference and then isolate the nonlinearity-induced folding in the FRFT domain. The [16] needs to compute the N-order difference and recover the high-order difference from the model sample. In comparison to [16], the number of differentiations is much lower in the proposed method.
The unlimited sampling theorem proves that the non-bandlimited signal in the Fourier domain based on the FRFT can be recovered from analog sampling as long as it meets Equation (29), whose amplitude exceeds the ADC threshold by orders of magnitude. It is particularly important that the signal is not affected by the ADC threshold.

5. Potential Application

5.1. Self-Reset ADC

The standard process of digitizing a signal involves bandwidth limiting by an annihilation filter, followed by sampling with an ADC. In practice, the C-ADC has a limited dynamic range [ λ , λ ] ; if the input signal exceeds this range, the signal is clipped. The clipping problem is a serious problem that manifests itself in the form of non-linear artifacts in audio-visual data and applications involving the sampling of biomedical data. The S-ADCs are rapidly developing with advances in ADC design. The theoretical concept of the S-ADCs has been in the literature since as early as the late 1970s, and physical implementations have only been in development since the early 2000s. One of the earliest references was by Rhee and Joo [37], where the authors proposed the S-ADC in the context of CMOS image sensors. When the upper or lower saturation threshold is reached, it resets to the corresponding threshold so that subsequent changes are captured even if the saturation limit is exceeded. The conceptual difference between the C-ADC and S-ADC is shown in Figure 2. The aim of developing the S-ADC is to enable the dynamic range of natural images in real-world applications to exceed the range that the C-ADC can handle. This capability is critical not only for consumer photography, but also for life sciences and bioimaging.
The work presented in this study revolves around the theoretical aspects of unlimited sampling. Based on the good properties of the FRFT, we propose a novel fractional Fourier domain sampling theorem that isolates the folds introduced by the modulus nonlinearity in the fractional Fourier domain, which can deal with the non-idealities and uncertainties introduced by the hardware while ensuring a low sampling rate. It is shown in theory that, based on the folding properties of the S-ADC, we can recover some signals up to multiples of the ADC threshold λ and can deal with arbitrarily long folding times.
  • In the practical application of the S-ADC, the proposed method can satisfy the condition that λ is unknowable, so it can resist any non-ideality.
  • The proposed method only needs to calculate the first-order differential, which is especially useful for the S-ADC in case of errors.
  • The proposed mathematical model in this study has a certain possibility in S-ADC. The proposed signals in this report are periodic bandlimited signals in the fractional Fourier domain, and this limitation reflects a practical limitation. Typically, instead of sampling on an ideal real line, the signal is sampled at finite intervals.
Based on the folding characteristics of the S-ADC hardware, the core of the proposed method is a mathematical model of the FRFT unlimited sampling signal. The folding introduced by the modulus nonlinearity can be isolated in the fractional Fourier domain, leading to frequency estimation problems. An annihilating filter estimation method is used to deal with arbitrarily close folding instants. In future work, we are committed to combining theory and practice, and have achieved breakthrough results in both hardware and algorithms.

5.2. Future Directions

  • Based on the mathematical conclusions obtained in this report, we will study a series of simulation experiments with the proposed theory. Some practical application ADC examples using the proposed unlimited sampling theorem will also be investigated in the future.
  • In the future, we will study the optimal fractional order angle α selection method in sampling theory and tell readers how to determine the optimal parameter α of the FRFT in the reconstruction of the analog signal from its sampled signal.
  • The broader signal transformation domain is also a topic worthy of our attention in the future. Extending our results to a wider domain of transforms, such as linear canonical transform, linear canonical wavelet transform, and canonical S-transform, etc., is a very interesting follow-up question.

6. Conclusions

In this report, we study the sampling theorem of bandlimited signals in the fractional Fourier domain based on the unlimited sampling framework of modulo measurement. Our main work is to perform modular operations in the fractional Fourier domain with the folding introduced by modular nonlinearity, and then to deal with the problem of fractional spectrum estimation. It turns out that the unlimited sampling theorem has nothing to do with the modulus threshold and can handle arbitrary folding time. In future work, we are committed to combining theory and practice and have achieved breakthrough results in both hardware and algorithms.

Author Contributions

Conceptualization, H.Z. and B.-Z.L.; formal analysis, H.Z. and B.-Z.L.; investigation, H.Z. and B.-Z.L.; writing—original draft preparation, H.Z. and B.-Z.L.; funding acquisition, B.-Z.L. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Natural Science Foundation of China [No. 62171041].

Data Availability Statement

Not applicables.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Shannon, C.E. Communication in the presence of noise. Proc. IRE 1949, 37, 447–457. [Google Scholar] [CrossRef]
  2. Nyquist, H. Introduction to Certain topics in telegraph transmission theory. Proc. IEEE 2002, 90, 276–279. [Google Scholar] [CrossRef]
  3. Butzer, P.L.; Stens, R.L. Sampling theory for not necessarily band-limited functions: A historical overview. SIAM Rev. 1992, 34, 40–53. [Google Scholar] [CrossRef]
  4. Unser, M. Sampling-50 years after Shannon. Proc. IEEE 2000, 88, 567–587. [Google Scholar] [CrossRef]
  5. Blu, T.; Unser, M. Quantitative Fourier analysis of approximation techniques. I. Interpolators and projectors. IEEE Trans. Signal Process. 1999, 47, 2796–2806. [Google Scholar] [CrossRef]
  6. Jing, R.M.; Feng, Q.; Li, B.Z. Higher order derivative sampling associated with fractional Fourier transform. Circuits Syst. Signal Process. 2019, 38, 1751–1774. [Google Scholar] [CrossRef]
  7. Li, B.Z.; Tao, R.; Wang, Y. New sampling formula related to linear canonical transform. Signal Process. 2007, 87, 983–990. [Google Scholar] [CrossRef]
  8. Butzer, P.L.; Dodson, M.M.; Ferreira, P.J.S.G.; Higgins, J.R.; Schmeisser, G.; Stens, R.L. Seven pivotal theorems of Fourier analysis, signal analysis, numerical analysis and number theory: Their interconnections. Bull. Math. Sci. 2014, 4, 481–525. [Google Scholar] [CrossRef]
  9. Kwon, U.K.; Kim, D.; Im, G.H. Amplitude clipping and iterative reconstruction of MIMO-OFourier domainM signals with optimum equalization. IEEE Trans. Wirel. Commun. 2009, 8, 268–277. [Google Scholar] [CrossRef]
  10. Kollár, Z.; Gazda, J.; Horvxaxth, P.; Varga, L.; Kocur, D. Iterative signal reconstruction of deliberately clipped SMT signals. Sci. China Inf. Sci. 2014, 56, 1–13. [Google Scholar] [CrossRef]
  11. Abel, J.S.; Smith, J.O. Restoring a clipped signal. In Proceedings of the IEEE International Conference on Acoustics, Speech, and Signal Processing, Toronto, ON, Canada, 14–17 April 1991; pp. 1745–1748. [Google Scholar]
  12. Kester, W. MT-025 tutorial ADC architectures VI: Folding ADCs. Analog Devices Tech. Rep. 2014, 27, 623–656. [Google Scholar]
  13. Park, D.; Rhee, J.; Joo, Y. A wide dynamic-range CMOS image sensor using self-reset technique. IEEE Electron Device Lett. 2007, 28, 890–892. [Google Scholar] [CrossRef]
  14. Park, D.; Joo, Y.; Koppa, S. A simple and robust self-reset CMOS image sensor. In Proceedings of the Midwest Symposium on Circuits and Systems, Seattle, WA, USA, 1–4 August 2010; pp. 680–683. [Google Scholar]
  15. Bhandari, A.; Krahmer, F.; Raskar, R. On unlimited sampling. In Proceedings of the 2017 International Conference on Sampling Theory and Applications, Tallinn, Estonia, 3–7 July 2017; pp. 1–26. [Google Scholar]
  16. Bhandari, A.; Krahmer, F.; Raskar, R. On unlimited sampling and reconstruction. IEEE Trans. Signal Process. 2020, 69, 3827–3839. [Google Scholar] [CrossRef]
  17. Shah, V.; Hegde, C. Sparse signal recovery from modulo observations. EURASIP J. Adv. Signal Process. 2021, 15, 1–17. [Google Scholar] [CrossRef]
  18. Bhandari, A.; Krahmer, F.; Raskar, R. Unlimited sampling of sparse sinusoidal mixtures. In Proceedings of the Symposium on Information Theory, Vail, CO, USA, 17–22 June 2018; pp. 336–340. [Google Scholar]
  19. Rencker, L.; Bach, F.; Wang, W.; Plumbley, M.D. Sparse recovery and dictionary learning from nonlinear compressive measurements. IEEE Trans. Signal Process. 2019, 67, 5659–5670. [Google Scholar] [CrossRef]
  20. Chen, W.; Fu, Z.; Grafakos, L.; Wu, Y. Fractional Fourier transforms on Lp and applications. Comput. Appl. Math. 2021, 55, 71–96. [Google Scholar]
  21. Romanov, E.; Ordentlich, O. Above the Nyquist rate, modulo folding does not hurt. IEEE Signal Process. Lett. 2019, 26, 1167–1171. [Google Scholar] [CrossRef]
  22. Ordentlich, O.; Tabak, G.; Hanumolu, P.K.; Singer, A.C.; Wornell, G.W. A modulo-based architecture for analog-to-digital conversion. IEEE J. Sel. Top. Signal Process. 2018, 12, 825–840. [Google Scholar] [CrossRef]
  23. Musa, O.; Jung, P.; Goertz, N. Generalized approximate message passing for unlimited sampling of sparse signals. In Proceedings of the 2018 IEEE Global Conference on Signal and Information Processing, Anaheim, CA, USA, 26–29 November 2018; pp. 336–340. [Google Scholar]
  24. Rudresh, S.; Adiga, A.; Shenoy, B.A.; Seelamantula, C.S. Wavelet-based reconstruction for unlimited sampling. In Proceedings of the 2018 IEEE International Conference on Acoustics, Speech and Signal Processing, Calgary, AB, Canada, 15–20 April 2018; pp. 4584–4588. [Google Scholar]
  25. Bhandari, A.; Krahmer, F.; Poskitt, T. Unlimited sampling from theory to practice: Fourier-Prony recovery and prototype ADC. IEEE Trans. Signal Process. 2021, 70, 1131–1141. [Google Scholar] [CrossRef]
  26. Xia, X.G. On bandlimited signals with fractional Fourier transform. IEEE Signal Process. Lett. 1996, 3, 72–74. [Google Scholar] [CrossRef]
  27. Almeida, L.B. The fractional Fourier transform and time-frequency representation. IEEE Trans. Signal Process. 1994, 42, 3084–3091. [Google Scholar] [CrossRef]
  28. Ozaktas, H.M.; Aytür, O. Fractional Fourier domains. Signal Process. 1995, 46, 119–124. [Google Scholar] [CrossRef]
  29. Pei, S.C.; Yeh, M.H.; Luo, T.L. Fractional Fourier series expansion for finite signals and dual extension to discrete-time fractional Fourier transform. IEEE Trans. Signal Process. 1999, 47, 2883–2888. [Google Scholar]
  30. Bhandari, A.; Marziliano, P. Sampling and reconstruction of sparse signals in fractional Fourier domain. IEEE Signal Process. Lett. 2010, 17, 221–224. [Google Scholar] [CrossRef]
  31. Itoh, K. Analysis of the phase unwrapping algorithm. Appl. Opt. 1982, 21, 2470. [Google Scholar] [CrossRef] [PubMed]
  32. Stoica, P.; Moses, R. Introduction to Spectral Analysis; Prentice-Hall: Englewood Cliffs, NJ, USA, 2000. [Google Scholar]
  33. Stoica, P.; Moses, R. Spectral Analysis of Signals; Prentice Hall: London, UK, 2005. [Google Scholar]
  34. Vetterli, M.; Marziliano, P.; Blu, T. Sampling signal with finite rate of innovation. IEEE Trans. Signal Process. 2002, 50, 1417–1428. [Google Scholar] [CrossRef]
  35. Roy, R.; Kailath, T. Esprit-estimation of signal parameters via rotational invariance technique. IEEE Trans. Acoust. Speech. Signal Process. 1989, 37, 984–995. [Google Scholar] [CrossRef]
  36. Wang, Y.Y.; Chen, J.T.; Fang, W.H. TST-MUSIC for joint DOA-delay estimation. IEEE Trans. Signal Process. 2001, 49, 721–729. [Google Scholar] [CrossRef]
  37. Rhee, J.; Joo, Y. Wide dynamic range CMOS image sensor with pixel level ADC. Electron. Lett. 2003, 39, 360–361. [Google Scholar] [CrossRef]
Figure 1. Unlimited sampling architecture for obtaining modulo samples.
Figure 1. Unlimited sampling architecture for obtaining modulo samples.
Fractalfract 07 00338 g001
Figure 2. Conceptual difference between the C−ADC and S−ADC. (a) Input signal and output signal in C−ADC; (b) input signal and output signal in S−ADC.
Figure 2. Conceptual difference between the C−ADC and S−ADC. (a) Input signal and output signal in C−ADC; (b) input signal and output signal in S−ADC.
Fractalfract 07 00338 g002
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Zhao, H.; Li, B.-Z. Unlimited Sampling Theorem Based on Fractional Fourier Transform. Fractal Fract. 2023, 7, 338. https://doi.org/10.3390/fractalfract7040338

AMA Style

Zhao H, Li B-Z. Unlimited Sampling Theorem Based on Fractional Fourier Transform. Fractal and Fractional. 2023; 7(4):338. https://doi.org/10.3390/fractalfract7040338

Chicago/Turabian Style

Zhao, Hui, and Bing-Zhao Li. 2023. "Unlimited Sampling Theorem Based on Fractional Fourier Transform" Fractal and Fractional 7, no. 4: 338. https://doi.org/10.3390/fractalfract7040338

APA Style

Zhao, H., & Li, B. -Z. (2023). Unlimited Sampling Theorem Based on Fractional Fourier Transform. Fractal and Fractional, 7(4), 338. https://doi.org/10.3390/fractalfract7040338

Article Metrics

Back to TopTop