Preprint
Article

Randomized Nonuniform Sampling for Random Signals Bandlimited in the Special Affine Fourier Transform Domain

Altmetrics

Downloads

79

Views

20

Comments

0

A peer-reviewed article of this preprint also exists.

This version is not peer-reviewed

Submitted:

12 January 2024

Posted:

12 January 2024

You are already at the latest version

Alerts
Abstract
The special affine Fourier transform (SAFT) is a useful and powerful analyzing tool in signal processing, optics and communications. In this paper, we mainly discuss the randomized nonuniform sampling and reconstruction for random signals bandlimited in the SAFT domain. First, we show that the nonuniform sampling is identical to the uniform sampling after a pre-filter in the sense of second order statistic characters. Then, we propose an approximate reconstruction based on sinc interpolation for the nonuniform sampling of random signals bandlimited in the SAFT domain. Finally, we give the mean square error estimate for the proposed approximate recovery approach.
Keywords: 
Subject: Computer Science and Mathematics  -   Analysis

1. Introduction

The special affine Fourier transform (SAFT) was firstly proposed in [1] to model optical systems. It offers a unified viewpoint of known signal processing transforms, such as Fourier transform (FT), fractional Fourier transform (FrFT), linear canonical transform (LCT), Laplace transform (LT) and so on. It can also include some optical operations on light waves, such as rotation, magnification, hyperbolic transformation, free space propagation, Lens transformation and so on. The SAFT is a six-parameters linear integral transform which is defined by offsetting two extra parameters on the basis of the LCT, so SAFT is also known as the offset linear canonical transform (OLCT). It has been proved that the SAFT is a useful tool for signal processing, communications, quantum mechanics and optics [12,16,23,26]. Many classical results such as Zak transform, Poisson summation formula and convolution theorems are established in the SAFT domain [6,24,33].
Let
A = a b u 0 c d ω 0
be a matrix with six real parameters satisfying a d b c = 1 . The continuous-time SAFT associated with the parameter matrix A of a signal f ( t ) is defined as in [1],
F A ( u ) = S A F T [ f ] ( u ) = + f ( t ) K A ( t , u ) d t , b 0 , d e j c d u u 0 2 2 + j w 0 u f d u u 0 , b = 0 ,
where the kernel function K A ( t , u ) is given by
K A ( t , u ) = 1 2 π j b e j d u 0 2 2 b e j 2 b a t 2 + 2 t u 0 u 2 u d u 0 b w 0 + d u 2 .
It is noted that when b = 0 , the SAFT of a signal is essentially a chirp multiplication. Therefore, we shall confine our attention to the case of b 0 . The inverse SAFT is expressed as
f ( t ) = C + F A ( u ) K A 1 ( u , t ) d u ,
where C = e j 2 c d u 0 2 2 a d u 0 w 0 + a b w 0 2 and
A 1 : = d b b ω 0 d u 0 c a c u 0 a ω 0 .
Sampling is one of the most fundamental process in digital signal processing which provides a bridge between the continuous physical signals and the discrete digital signals. Beginning with the Shannon’s sampling theorem of bandlimited signals [15], various sampling such as nonuniform sampling, average sampling, dynamic sampling, random sampling, mobile sampling, timing sampling and multi-channel sampling have been generally studied for signals bandlimited in the FT domain [2,3,5,9]. With the appearance and developments of the more general transforms, the corresponding sampling theories are extended to the signals bandlimited in the FrFT, LCT and SAFT domains [6,12,14,18,19,21,22,23,25,26,27,30,32].
Signals in the real world often presents random characteristics and sampling for random signals bandlimited in the FT domain has been generally studied [5,7,8,17]. In recent years, there have existed many researches for sampling of random signals bandlimited in the FrFT and LCT domains [10,11,20,28,31]. The uniform sampling theorems in [10] was extended to the SAFT domain as in [29]. Nonuniform sampling is a more realistic sampling scheme due to the limitations of data acquisition and processing ability. In fact, the nonuniform sampling theories including the periodic nonuniform sampling model, N-order recurrent nonuniform sampling model, nonuniform sampling due to migration of a finite number of uniform samples and the general nonuniform sampling have been given for signals bandlimited in the LCT domain [31] and signals bandlimited in the SAFT domain [4,30], respectively. In particular, the nonuniform sampling problem was also considered in [11] for random signals bandlimited in the LCT domain, where a randomized nonuniform sampling method and a class of approximate recovery approaches by using sinc interpolation functions were studied. In this paper, we will further study the randomized nonuniform sampling for random signals bandlimited in the SAFT domain and also give an approximate recovery method based on the sinc interpolation.
The paper is organized as follows. In section 2, we give the definition of the power spectral density in the SAFT domain. In section 3, we study the nonuniform sampling scheme and propose an approximate recovery approach. In section 4, the mean square error estimate for the proposed approximate recovery method is demonstrated.

2. Power Spectral Density in the SAFT Domain

Given a probability space ( Ω , F , p ) , a stochastic process x ( t ) is said to be wide sense stationary if it has zero mean and its auto-correlation function
R x x ( t + τ , t ) = E x ( t + τ ) x * ( t )
is independent of t R , i.e., R x x ( t + τ , t ) = R x x ( τ ) , where E [ · ] denotes mathematical expectation and the superscript * stands for the complex conjugate. Two stochastic processes x ( t ) and y ( t ) are said to be jointly stationary, if x ( t ) and y ( t ) are both stationary and their cross-correlation function
R x y ( t + τ , t ) = E x ( t + τ ) y * ( t )
is independent of t R , i.e., R x y ( t + τ , t ) = R x y ( τ ) .
We next introduce the SAFT auto-correlation function, the SAFT cross-correlation function, the SAFT auto-power spectral density and the SAFT cross-power spectral density as in [29]. For two random signals x ( t ) and y ( t ) , the SAFT auto-correlation function of x ( t ) is defined as
R x x A ( τ ) = lim T + 1 2 T T T R x x ( t + τ , t ) e j a b t τ d t .
Similarly, the SAFT cross-correlation function of x ( t ) and y ( t ) is defined as
R x y A ( τ ) = lim T + 1 2 T T T R x y ( t + τ , t ) e j a b t τ d t .
Remark 2.1.If the random signal x ˜ ( t ) = x ( t ) e j a 2 b t 2 is stationary, then x 1 ( t ) = x ˜ ( t ) e j u 0 b t is also stationary. In fact,
R x 1 x 1 ( t + τ , t ) = e j u 0 τ b R x ˜ x ˜ ( t + τ , t ) .
Moreover, one has
R x ˜ x ˜ ( t + τ , t ) = E x ˜ ( t + τ ) x ˜ * ( t ) = E x ( t + τ ) e j a 2 b ( t + τ ) 2 x * ( t ) e j a 2 b t 2 = E x ( t + τ ) x * ( t ) e j a 2 b τ 2 e j a t τ b = R x x ( t + τ , t ) e j a t τ b e j a 2 b τ 2 .
Therefore, R x x ( t + τ , t ) e j a t τ b must be independent of t. In such case, we have
R x x A ( τ ) = R x 1 x 1 ( τ ) e j a 2 b τ 2 e j u 0 b τ .
Define the SAFT auto-power spectral density of the random signal x ( t ) by
P x x A ( u ) = 1 j 2 π b e j d 2 b u 2 e j d u 0 2 2 b e j u b d u 0 b w 0 F A R x x A ( τ ) ( u )
and the SAFT cross-power spectral density of the random signals x ( t ) and y ( t ) as
P x y A ( u ) = 1 j 2 π b e j d 2 b u 2 e j d u 0 2 2 b e j u b d u 0 b w 0 F A R x y A ( τ ) ( u ) .
It follows from (1.1) and (2.6) that
R x x A ( τ ) = C · + P x x A u 1 / 1 j 2 π b   e j d 2 b u 2 e j d 2 b u 0 2 e j u b d u 0 b w 0 1 j 2 π b   × e j a 2 b ( b w 0 d u 0 ) 2 e j d 2 b u 2 e j u b b w 0 d u 0 τ e j u 0 b τ e   j a 2 b τ 2 d u   = C · + P x x A u e j d 2 b u 0 2 e j a 2 b ( b w 0 d u 0 ) 2 e j u b τ e j u 0 b τ e   j a 2 b τ 2 d u   = e j 2 c d u 0 2 2 a d u 0 w 0 + a b w 0 2 + P x x A u e j d 2 b u 0 2 e j a 2 b ( b w 0 d u 0 ) 2 e j u b τ e j u 0 b τ e j a 2 b τ 2 d u   = + P x x A u e j a 2 b τ 2 e j 1 b ( u u 0 ) τ d u .
Multiplicative filtering in the SAFT domain is showed in Figure 1, which has been introduced in [29]. More specifically, we first obtain the SAFT of the input signal f 1 ( t ) and apply the multiplicative filter H ( u ) in the SAFT domain. Then the output signal f 2 ( t ) in the time domain is obtained by the inverse SAFT. Mathematically, the output f 2 ( t ) is given by
f 2 ( t ) = F A 1 F 2 ( u ) ( t ) = F A 1 F 1 ( u ) H ( u ) ( t ) ,
where F 1 ( u ) = F A f 1 ( t ) ( u ) and F 2 ( u ) = F 1 ( u ) H ( u ) .
Define normalized convolution
( f Θ g ) ( t ) = 1 2 π R f ( x ) g ( t x ) e j a 2 b t 2 x 2 d x
for f , g L 2 ( R ) [23]. Then we have the following conclusion.
Proposition 2.2. 
Let
H ( u ) = 1 2 π + h ( t ) e j u u 0 b t d t .
Then the multiplicative filter in Figure 1 is equivalent to
f 2 ( t ) = f 1 Θ h ( t ) .
Proof we only need to prove
F A f 1 Θ h ( t ) ( u ) = F 1 ( u ) H ( u ) .
It follows from the definition of the SAFT that
F A f 1 Θ h ( t ) ( u ) = R f 1 Θ h ( t ) 1 2 π j b e j d u 0 2 2 b e j 2 b a t 2 + 2 t u 0 u 2 u d u 0 b w 0 + d u 2 d t = 1 2 π R R f 1 ( x ) h ( t x ) e j a 2 b t 2 x 2 1 2 π j b e j d u 0 2 2 b e j 2 b a t 2 + 2 t u 0 u 2 u d u 0 b w 0 + d u 2 d x d t = R f 1 ( x ) 1 2 π j b e j d u 0 2 2 b e j 2 b a x 2 + 2 x u 0 u 2 u d u 0 b w 0 + d u 2 d x 1 2 π R h ( t x ) e j b t u 0 u e j 2 b 2 x u 0 u d t = R f 1 ( x ) 1 2 π j b e j d u 0 2 2 b e j 2 b a x 2 + 2 x u 0 u 2 u d u 0 b w 0 + d u 2 d x 1 2 π R h ( t ) e j u u 0 b t d t = F 1 ( u ) H ( u ) .
Lemma 2.3. 
[29] Suppose that the random signals x ( t ) and y ( t ) are the input and the output in Figure 1, then
P x y A ( u ) = H ( u ) P x x A ( u )
and
P y y A ( u ) = | H ( u ) | 2 P x x A ( u ) .

3. Nonuniform Sampling and Approximate Recovery

In this section, we will study the nonuniform sampling and reconstruction of random signals which are bandlimited in the SAFT domain.
Definition 3.1. 
[29] We say that a random signal x ( t ) is bandlimited in the SAFT domain if its SAFT power spectral density P x x A ( u ) satisfies
P x x A ( u ) = 0 , | u | > u r ,
where u r is called the bandwidth of the random signal x ( t ) in the SAFT domain.
Lemma 3.2. 
Assume that a random signal x ( t ) is bandlimited in the SAFT domain with bandwidth u r and x ˜ ( t ) = x ( t ) e j a 2 b t 2 is stationary. Then x 1 ( t ) is bandlimited in the FT domain with bandwidth u r b and the power spectral density satisfies supp { P x 1 x 1 ( u ) } [ u r b , u r b ] .
Proof Since x 1 ( t ) is stationary, it follows from (2.5) and (2.6) that
P x x A u = 1 j 2 π b e j d 2 b u 2 e   j d 2 b u 0 2 e j u b d u 0 b w 0 F A R x x A ( τ ) u = 1 j 2 π b e j d 2 b u 2 e   j d 2 b u 0 2 e j u b d u 0 b w 0 F A R x 1 x 1 τ e j a 2 b τ 2 e j u 0 b τ u = 1 j 2 π b e j d 2 b ( u 2 + u 0 2 ) e j u b d u 0 b w 0 + R x 1 x 1 τ e j a 2 b τ 2 e j u 0 b τ × 1 2 π j b e j d u 0 2 2 b e j 2 b a τ 2 + 2 τ u 0 u 2 u d u 0 b w 0 + d u 2 d τ = 1 2 π b + R x 1 x 1 τ e j u b τ d τ
= 1 2 π b P x 1 x 1 u b .
Note that P x x A ( u ) = 0 , | u | > u r . Then the desired result is proved. First, we will show that the nonuniform sampling is identical to uniform sampling after a pre-filter in the sense of second order statistic characters.
Theorem 3.3. 
Suppose that the random signal x ( t ) is bandlimited in the SAFT domain with bandwidth u r and x ˜ ( t ) = x ( t ) e j a 2 b t 2 is stationary. Then in the sense of second order statistic characters, the nonuniform sampling of x ( t ) at the sampling points t n = n T + ξ n (Figure 2 ) is identical to the uniform sampling after a SAFT filter h 1 ( t ) as in Figure 3, where T T N = π b u r , { ξ n } is a sequence of independent identically distributed random variables with zero mean in the interval ( T / 2 , T / 2 ) . Moreover,
H 1 ( u ) = ϕ ξ u b = 1 2 π + h 1 ( t ) e j u u 0 b t d t
and ϕ ξ ( u ) denotes the characteristic function of ξ n .
Proof Note that y ( t ) = x Θ h 1 ( t ) . Then it follows from Lemma 2.2 that
P y y A ( u ) = H 1 ( u ) 2 P x x A ( u ) .
Moreover, one has
y 1 ( t ) = y ( t ) e j a 2 b t 2 e j u 0 b t   = 1 2 π e j u 0 b t R x ( s ) h 1 ( t s ) e j a 2 b s 2 d s = 1 2 π e j u 0 b t R x ˜ ( s ) h 1 ( t s ) d s .
Hence, we have
R y 1 y 1 ( t + τ , t ) = 1 2 π e j u 0 b τ R R h 1 * ( s ) R x ˜ x ˜ τ s + s h 1 s d s d s ,
which is independent of t and y 1 ( t ) is stationary. It follows from (2.8) and (3.3) that
R y y A ( k T ) = u r u r P y y A ( u ) e j a 2 b ( k T ) 2 + j 1 b u u 0 k T d u = u r u r H 1 ( u ) 2 P x x A ( u ) e j a 2 b ( k T ) 2 + j 1 b u u 0 k T d u .
This together with (2.5) obtains
R y 1 y 1 ( n T , ( n k ) T ) = R y 1 y 1 ( k T )   = R y y A ( k T ) e j a 2 b ( k T ) 2 e j u 0 b k T   = e j a 2 b ( k T ) 2 e j u 0 b ( k T ) u r u r H 1 ( u ) 2 P x x A ( u ) e j a 2 b ( k T ) 2 + j 1 b u u 0 k T d u = u r u r H 1 ( u ) 2 P x x A ( u ) e j u b k T d u .
Combining (2.5) and (2.8), we have
E R x 1 x 1 k T + ξ n ξ n k = E R x x A k T + ξ n ξ n k e j a 2 b k T + ξ n ξ n k 2 e j u 0 b k T + ξ n ξ n k   = E u r u r P x x A ( u ) e j a 2 b k T + ξ n ξ n k 2 + j 1 b u u 0 k T + ξ n ξ n k d u · e j a 2 b k T + ξ n ξ n k 2 e j u 0 b k T + ξ n ξ n k   = E u r u r P x x A ( u ) e j u b k T + ξ n ξ n k d u = u r u r P x x A ( u ) e j u b k T E e j u b ξ n ξ n k d u .
Let Z = ξ n ξ n k and f Z ( η ) be the probability density function of Z. Note that ξ n and ξ n k are independent and have identical distributions. Let f ξ ( η ) be their common probability density function. Then we have
f Z ( η ) = f ξ ( · ) * f ξ ( · ) ( η ) ,
where * denotes the convolution operator. Moreover, one has
E e j u b ξ n ξ n k = + f Z ( η ) e j u b η d η   = + f ξ ( · ) * f ξ ( · ) ( η ) e j u b η d η   = + f ξ ( η ) e j u b η d η · + f ξ ( η ) e j u b η d η = ϕ ξ u b 2 ,
where
ϕ ξ ( u ) = + f ξ ( η ) e j u η d η .
Substituting (3.9) into (3.7) obtains
E R x 1 x 1 k T + ξ n ξ n k = u r u r ϕ ξ u b 2 P x x A ( u ) e j u b k T d u .
This together with H 1 ( u ) = ϕ ξ u b and (3.6) proves the desired result.
In the following, we will give an approximate recovery method for bandlimited signals in the SAFT domain based on randomized nonuniform samples.
Lemma 3.4. 
[13] Suppose that the random signal x ( t ) is bandlimited in the Fourier transform domain with bandwidth u r b , { ξ n } and { ζ n } are two sequences of independent identically distributed random variables with zero mean. Then an approximate recovery formula of nonuniform sampling for the random signal x ( t ) can be represented by
x ( t ) = T T N n = + x t n h 2 t t ˜ n ,
where h 2 ( t ) = s i n c u r t b , s i n c ( x ) sin x x , t n = n T + ξ n and t ˜ n = n T + ζ n .
Theorem 3.5. 
Suppose that the random signal x ( t ) is bandlimited in the SAFT domain with bandwidth u r and x ˜ ( t ) = x ( t ) e j a 2 b t 2 is stationary. Then x ( t ) can be approximated from its nonuniform samples by utilizing the sinc interpolation function as
x ^ ( t ) = T T N e j u 0 b t e j a 2 b t 2 n = + x t n e j a 2 b t n 2 e j u 0 b t n h 2 t t ˜ n ,
where t n and t ˜ n are as in Lemma 3.4.
Proof It follows from Lemma 3.2 that x 1 ( t ) is bandlimited in the FT domain with bandwidth u r b . By (3.11), we know that
x 1 ( t ) = T T N n = + x 1 ( t n ) h 2 t t ˜ n = T T N n = + x t n e j a 2 b t n 2 e j u 0 b t n h 2 t t ˜ n
is an approximation of x 1 ( t ) . Note that x ( t ) = e j u 0 b t e j a 2 b t 2 x 1 ( t ) . Then x ^ ( t ) in (3.12) is an approximate recovery approach of x ( t ) and the proof is completed.
From Theorem 3.5, one can see that the approximate recovery approach using the sinc interpolation for a random signal that is bandlimited in the SAFT domain can be expressed in Figure 4.

4. Error estimate for Nonuniform Sampling

Since the reconstruction with randomized sinc interpolation is an approximate method, we will estimate the approximation error in this section.
Lemma 4.1. 
Let random signals x 1 ( t ) and y 1 ( t ) be the input and output of the FT multiplicative filter as in Figure 5. Then
P y 1 y 1 ( u ) = h ^ 3 ( u ) 2 P x 1 x 1 ( u ) ,
where h ^ 3 ( u ) is the FT of h 3 ( t ) , that is,
h ^ 3 ( u ) = R h 3 ( t ) e j u t d t .
Proof Note that y 1 ( t ) = + x 1 ( t u ) h 3 ( u ) d u . Then
R y 1 x 1 t + τ , t = E y 1 ( t + τ ) x 1 * ( t ) = + R x 1 x 1 τ u h 3 ( u ) d u ,
which is independent of t. Moreover, one has
R y 1 y 1 t + τ , t = E y 1 ( t + τ ) y 1 * ( t ) = + R y 1 x 1 τ + u h 3 * ( u ) d u .
Taking FT on both sides of (4.1) and (4.2) obtains
P y 1 x 1 u = h ^ 3 ( u ) P x 1 x 1 u
and
P y 1 y 1 u = h ^ 3 * ( u ) P y 1 x 1 u .
Combining (4.3) and (4.4) gives
P y 1 y 1 ( u ) = h ^ 3 ( u ) 2 P x 1 x 1 ( u ) .
Theorem 4.2. 
Suppose that the random signal x ( t ) is bandlimited in the SAFT domain with bandwidth u r and x ˜ ( t ) = x ( t ) e j a 2 b t 2 is stationary. Let v ( t ) be an additive noise with zero mean, which is stationary, uncorrelated with x ( t ) and has the power spectral density
P v v ( u ) = T u r u r P x x A ( u 1 ) 1 ϕ ξ ζ u 1 b , u 2 d u 1 , | u | u r b ,
where ϕ ξ ζ ( s , t ) is the joint characteristic function of the random variables ξ n and ζ n . If ϕ ξ ζ ( u , u ) is the frequency response of the filter h 3 ( t ) , then the model described in Figure 5 is identical to the procedure represented in Figure 4 in the sense of second order statistic characters. Moreover, we have
E x ^ ( t ) x ( t ) 2 = u r u r P x x A ( u ) 1 ϕ ξ ζ u b , u b 2 d u   + T 2 π b u r u r P x x A u u r u r 1 ϕ ξ ζ u b , u 1 b 2 d u 1 d u .
Proof It follows from Theorem 3.5 that
x ¯ ( t ) = x ^ ( t ) e j u 0 b t e j a 2 b t 2 = T T N n = + x ( t n ) e j a 2 b t n 2 e j u 0 t n b h 2 ( t t ˜ n ) .
Then one has
R x ¯ x ¯ ( t , t τ ) = T T N 2 E n = + x ( n T + ξ n ) e j a 2 b ( n T + ξ n ) 2 e j u 0 ( n T + ξ n ) b h 2 ( t n T ζ n ) ·   k = + x * ( k T + ξ k ) e j a 2 b ( k T + ξ k ) 2 e j u 0 ( k T + ξ k ) b h 2 * ( t τ k T ζ k )   = T T N 2 n = + k = + E R x 1 x 1 ( n T k T + ξ n ξ k ) h 2 ( t n T ζ n ) h 2 * ( t τ k T ζ k ) .
Moreover, it can be represented by two terms as
R x ¯ x ¯ ( t , t τ ) = T T N 2 R x 1 x 1 ( 0 ) n = + E h 2 ( t n T ζ n ) h 2 * ( t τ n T ζ n )     + T T N 2 n k   E R x 1 x 1 ( n T k T + ξ n ξ k ) h 2 ( t n T ζ n ) h 2 * ( t τ k T ζ k )   = Δ I + I I .
Note that n e j ( u 2 u 1 ) n T = 2 π k δ ( u 2 u 1 ) T 2 π k and
h 2 ( t ) = 1 2 π R H 2 ( u ) e j u b t d u = b 2 π R H 2 ( u b ) e j u t d u .
These together with the fact that H 2 ( u ) = π 2 π u r χ u r , u r ( u ) show that
I = 1 2 π T b T N 2 R x 1 x 1 ( 0 ) R R H 2 ( b u 1 ) H 2 * ( b u 2 ) e j ( u 1 u 2 ) t e j u 2 τ n = + e j ( u 2 u 1 ) n T E e j ( u 2 u 1 ) ζ n d u 1 d u 2 = T b T N 2 R x 1 x 1 ( 0 ) u r b u r b 1 T H 2 ( b u ) 2 e j u τ d u   = T b T N 2 1 2 π u r b u r b P x 1 x 1 ( u 1 ) d u 1 u r b u r b H 2 ( b u ) 2 e j u τ d u   = T 4 π 2 u r b u r b e j u τ u r b u r b P x 1 x 1 ( u 1 ) d u 1   d u .
Moreover, we have
I I = b T 2 π T N 2 n k   E u r b u r b P x 1 x 1 ( u ) e j u ( n T k T + ξ n ξ k ) d u R H 2 ( b u 1 ) e j u 1 ( t n T ζ n ) d u 1 · R H 2 * ( b u 2 ) e j u 2 ( t τ k T ζ k ) d u 2 ] = b T 2 π T N 2 n k   u r b u r b R R P x 1 x 1 ( u ) H 2 ( b u 1 ) H 2 * ( b u 2 ) e j u 2 τ e j ( u 1 u 2 ) t e j ( u u 1 ) n T e j ( u u 2 ) k T · E e j u ξ n e j u ξ k e j u 1 ζ n e j u 2 ζ k d u 1 d u 2 d u = b T 2 π T N 2 u r b u r b R R P x 1 x 1 ( u ) H 2 ( b u 1 ) H 2 * ( b u 2 ) ϕ ξ ζ u , u 1 ϕ ξ ζ * u , u 2 e j u 2 τ e j ( u 1 u 2 ) t · n e j ( u u 1 ) n T k e j ( u u 2 ) k T d u 1 d u 2 d u b T 2 π T N 2 u r b u r b R R P x 1 x 1 ( u ) H 2 ( b u 1 ) H 2 * ( b u 2 ) · ϕ ξ ζ u , u 1 ϕ ξ ζ * u , u 2 e j u 2 τ e j ( u 1 u 2 ) t n e j ( u 2 u 1 ) n T d u 1 d u 2 d u = b T N 2 u r b u r b P x 1 x 1 ( u ) ϕ ξ ζ u , u 2 H 2 ( b u ) 2 e j u τ d u T 2 π b T N 2 u r b u r b u r b u r b P x 1 x 1 ( u 1 ) ϕ ξ ζ u 1 , u 2 H 2 ( b u ) 2 e j u τ d u 1 d u   = b T N 2 u r b u r b H 2 ( b u ) 2 e j u τ P x 1 x 1 ( u ) ϕ ξ ζ u , u 2 T 2 π u r b u r b P x 1 x 1 ( u 1 ) ϕ ξ ζ u 1 , u 2 d u 1 du = 1 2 π u r b u r b e j u τ P x 1 x 1 ( u ) ϕ ξ ζ u , u 2 T 2 π u r b u r b P x 1 x 1 ( u 1 ) ϕ ξ ζ u 1 , u 2 d u 1 d u .
Substituting (4.9) and (4.10) into (4.7) obtains
R x ¯ x ¯ ( t , t τ ) = T 4 π 2 u r b u r b u r b u r b P x 1 x 1 ( u 1 ) 1 ϕ ξ ζ u 1 , u 2 d u 1 e j u τ d u   +   1 2 π u r b u r b e j u τ P x 1 x 1 ( u ) ϕ ξ ζ u , u 2 d u .  
Similarly, we can obtain
R x ¯ x 1 ( t , t τ ) = 1 2 π u r b u r b P x 1 x 1 ( u ) e j u τ ϕ ξ ζ ( u , u ) d u .
Therefore, we have
P x ¯ x ¯ ( u ) = P x 1 x 1 ( u ) ϕ ξ ζ u , u 2 + T 2 π u r b u r b P x 1 x 1 ( u 1 ) 1 ϕ ξ ζ u 1 , u 2 d u 1  
and
P x ¯ x 1 ( u ) = P x 1 x 1 ( u ) ϕ ξ ζ ( u , u ) .
It follows from Lemma 4.1 that the first term P x 1 x 1 ( u ) ϕ ξ ζ u , u 2 in (4.13) is the FT power spectral density of y 1 ( t ) in Figure 5. Furthermore, since x ¯ ( t ) = y 1 ( t ) + v ( t ) and v ( t ) is uncorrelated with x ( t ) , then
R x ¯ x ¯ ( t + τ , t ) = E y 1 ( t + τ ) + v ( t + τ ) y 1 ( t ) + v ( t ) *   = R y 1 y 1 ( t + τ , t ) + R y 1 v ( t + τ , t ) + R v y 1 ( t + τ , t ) + R v v ( t + τ , t )   = R y 1 y 1 ( t + τ , t ) + R v v ( t + τ , t ) .
Moreover, one has
P x ¯ x ¯ ( u ) = P y 1 y 1 ( u ) + P v v ( u ) ,
which shows that the second term in (4.13) is just the power spectral density of v ( t ) , that is,
P v v ( u ) = T 2 π u r b u r b P x 1 x 1 ( u 1 ) 1 ϕ ξ ζ u 1 , u 2 d u 1   = T u r u r P x x A ( u 1 ) 1 ϕ ξ ζ u 1 b , u 2 d u 1 , | u | u r b .
Therefore, the model described in Figure 5 is identical to the procedure represented in Figure 4 in the sense of second order statistic characters.
Next, we will estimate the error E x ^ ( t ) x ( t ) 2 . Let ε ( t ) = x ^ ( t ) x ( t ) . Combining (3.2) and (4.13), we get
P x ^ x ^ A u = 1 2 π b P x ¯ x ¯   ( u b )   = 1 2 π b P x 1 x 1 ( u b ) ϕ ξ ζ u b , u b 2 + T 4 π 2 b u r b u r b P x 1 x 1 ( u 1 ) 1 ϕ ξ ζ u 1 , u b 2   d u 1   = P x x A u ϕ ξ ζ u b , u b 2 + T 2 π b u r u r P x x A u 1 1 ϕ ξ ζ u 1 b , u b 2   d u 1 .  
Similarly, we can obtain
P x ^ x A u = P x x A u ϕ ξ ζ u b , u b .
In fact, it is easy to see that
R x ¯ x 1 t + τ , t = E x ^ ( t + τ ) e j u 0 b ( t + τ ) e j a 2 b ( t + τ ) 2 x * ( t ) e j u 0 b t e j a 2 b t 2   = R x ^ x t + τ , t e j u 0 b τ e j a b t τ e j a 2 b τ 2 .
Therefore, R x ^ x t + τ , t e j a b t τ is independent of t due to (4.12). Then
R x ^ x A τ = R x ^ x t + τ , t e j a b t τ   = R x ¯ x 1 t + τ , t e j u 0 b τ e j a 2 b τ 2 .
Moreover, it follows from (2.7) that
P x ^ x A u = 1 j 2 π b e j d 2 b u 2 e   j d 2 b u 0 2 e j u b d u 0 b w 0 F A R x ¯ x 1 τ e j u 0 b τ e j a 2 b τ 2 u   = 1 j 2 π b e j d 2 b ( u 2 + u 0 2 ) e j u b d u 0 b w 0 + R x ¯ x 1 τ e j u 0 b τ e j a 2 b τ 2 1 2 π j b e j d u 0 2 2 b   · e j 2 b a τ 2 + 2 τ u 0 u 2 u d u 0 b w 0 + d u 2 d τ   = 1 2 π b + R x ¯ x 1 τ e j u b τ d τ   = 1 2 π b P x ¯ x 1 u b   = 1 2 π b P x 1 x 1 u b ϕ ξ ζ u b , u b   = P x x A u ϕ ξ ζ u b , u b .  
Hence, the SAFT auto-power spectral density of the reconstruction error ε ( t ) is
P ε ε A u = P x ^ x ^ A u P x ^ x A u P x x ^ A u + P x x A u   = P x x A u ϕ ξ ζ u b , u b 2 + T 2 π b u r u r P x x A u 1 1 ϕ ξ ζ u 1 b , u b 2   d u 1     P x x A u ϕ ξ ζ u b , u b P x x A u ϕ ξ ζ u b , u b * + P x x A u = P x x A u 1 ϕ ξ ζ u b , u b 2 + T 2 π b u r u r P x x A u 1 1 ϕ ξ ζ u 1 b , u b 2   d u 1 ,
where we have used the fact that P x x A u is real due to (3.2). Note that
ε 1 ( t ) = ε ( t ) e j u 0 b t e j a 2 b t 2 = ( x ^ ( t ) x ( t ) ) e j u 0 b t e j a 2 b t 2 = x ¯ ( t ) x 1 ( t ) .
Then ε 1 ( t ) is stationary. Moreover, it follows from (2.5) and (2.8) that
E ε ( t ) 2 = R ε 1 ε 1 ( 0 ) = R ε ε A ( 0 ) = u r u r P ε ε A u d u = u r u r P x x A u 1 ϕ ξ ζ u b , u b 2 d u + T 2 π b u r u r P x x A u u r u r 1 ϕ ξ ζ u b , u 1 b 2 d u 1 d u .
This completes the proof.
()f ξ n and ζ n equal to zero, then the nonuniform sampling studied in this paper reduces to the uniform sampling. In such case, ϕ ξ ζ ( s , t ) 1 . Then it follows from Theorem 4.2 that
E | x ^ ( t ) x ( t ) | 2 = u r u r P x x A ( u ) 1 ϕ ξ ζ u b , u b 2 d u + T 2 π b u r u r P x x A ( u ) u r u r 1 ϕ ξ ζ u b , u 1 b 2 d u 1 d u = 0 .
That is to say, x ( t ) is equal to its approximation x ^ ( t ) in the mean square sense. From Theorem 3.5, one can see that for T = T N = π b u r , the approximation of x ( t ) obtained in (3.12) becomes
x ^ ( t ) = e j a 2 b t 2 n = + x ( n T ) e j a 2 b ( n T ) 2 e j u 0 b ( n T t ) s i n c u r ( t n T ) b ,
which coincides with Theorem 3 in [29]. Therefore, the result of uniform sampling proposed in [29] is a special case of Theorems 3.5 and 4.2 in this paper.

Acknowledgments

The project is partially supported by the National Natural Science Foundation of China (No.12261025), Center for Applied Mathematics of Guangxi (No. AD23023002), Guangxi Colleges and Universities Key Laboratory of Data Analysis and Computation.

Conflicts of Interest

The authors declare that they have no competing interests.

References

  1. S. Abe, J. T. Sheridan. Optical operations on wave functions as the Abelian subgroups of the special affine Fourier transformation. Optics Letters, 1994, 19: 1801–1803.
  2. K. Adam, A. Scholefield, M. Vetterli. Sampling and reconstruction of bandlimited signals with multi-channel time encoding. IEEE Transactions on Signal Processing, 2020, 69: 1105-1119.
  3. A. Aldroubi, J. Davis, I. Krishtal. Dynamic sampling: Time-space trade-off. Appl. Comput. Harmon. Anal., 2013, 34: 495–503.
  4. Y. Bao, Y. Zhang, Y. Song, B. Li, P. Dang. Nonuniform sampling theorems for random signals in the offset linear canonical transform domain. Proceedings of APSIPA Annual Summit and Conference, 2017, Kuala Lumpur, Malaysia, 94–99.
  5. R. F. Bass, K. Gröchenig. Random sampling of bandlimited functions. Isral J. Math., 2010, 177(1): 1-28.
  6. A. Bhandari, A. I. Zayed. Shift-invariant and sampling spaces associated with the special affine Fourier transform. Applied and Computational Harmonic Analysis, 2019, 47: 30–52.
  7. H. Boche, U. J. Mönich. Approximation of wide-sense stationary stochastic processes by Shannon sampling series. IEEE Transactions on Information Theory, 2010, 56(12): 6459–6469.
  8. F. Gilles, K. Sinuk. Average sampling of band-limited stochastic processes. Applied and Computational Harmonic Analysis, 2013, 35(3): 527–534.
  9. K. Gröchenig, J. L. Romero, J. Unnikrishnan, M. Vetterli. On minimal trajectories for mobile sampling of bandlimited fields. Appl. Comput. Harmon. Anal., 2015, 39: 487–510.
  10. H. Huo, W. Sun. Sampling theorems and error estimates for random signals in the linear canonical transform domain. Signal Processing, 2015, 111: 31–38.
  11. H. Huo, W. Sun. Nonuniform sampling for random signals bandlimited in the linear canonical transform domain. Multidimensional Systems and Signal Processing, 2020, 31: 927–950.
  12. X. Liu, J. Shi, W. Xiang, Q. Zhang, N. Zhang. Sampling expansion for irregularly sampled signals in fractional Fourier transform domain. Digital Signal Processing, 2014, 34: 74–81.
  13. S. Maymon, A. V. Oppenheim. Sinc interpolation of nonuniform samples. IEEE Trans. Signal Process., 2011, 59 (10): 4745–4758.
  14. F. Oktem, H. Ozaktas. Exact relation between continuous and discrete linear canonical transforms. IEEE Signal Processing Letters, 2009, 16(8): 727–730.
  15. C. E. Shannon. Communication in the presence of noise. Proc. IRE, 1949, 37: 10–12.
  16. K. Sharma, S. Joshi. Signal separation using linear canonical and fractional Fourier transforms. Optics Communications, 2006, 265(2): 454–460.
  17. Z. Song, W. Sun, X. Zhou, Z. Hou. An average sampling theorem for bandlimited stochastic processes. IEEE Transactions on Information Theory, 2007, 53(12): 4798–4800.
  18. A. Stern. Sampling of compact signals in the offset linear canonical transform domains. Signal, Image and Video Processing, 2007, 1(4): 359–367.
  19. R. Tao, B. Li, Y. Wang, G. Aggrey. On sampling of band-limited signals associated with the linear canonical transform. IEEE Transactions on Signal Processing, 2008, 56(11): 5454–5464.
  20. R. Tao, F. Zhang, Y. Wang. Sampling random signals in a fractional Fourier domain. Signal Processing, 2011, 91(6): 1394–1400.
  21. D. Wei, Q. Ran, Y. Li. Multichannel sampling and reconstruction of bandlimited signals in the linear canonical transform domain. IET Signal Processing, 2011, 5(8): 717–727.
  22. H. Zhao, R. Wang, D. Song. Recovery of bandlimited signals in linear canonical transform domain from noisy samples. Circuits Systems and Signal Processing, 2014, 33(6): 1997–2008.
  23. H. Zhao, R. Wang, D. Song, T. Zhang, D. Wu. Extrapolation of discrete bandlimited signals in linear canonical transform domain. Signal Processing, 2014, 94: 212–218.
  24. Q. Xiang, K. Qin. Convolution, correlation, and sampling theorems for the offset linear canonical transform. Signal, Image and Video Processing, 2014, 8(3): 433–442.
  25. Q. Xiang, K. Qin, Q. Huang. Multichannel sampling of signals band-limited in offset linear canonical transform domains. Circuits Systems and Signal Processing, 2013, 32(5): 2385–2406.
  26. S. Xu, Y. Chai, Y. Hu. Spectral analysis of sampled band-limited signals in the offset linear canonical transform domain. Circuits Systems and Signal Processing, 2015, 34(12): 3979–3997.
  27. S. Xu, Y. Chai, Y. Hu, C. Jiang, Y. Li. Reconstruction of digital spectrum from periodic nonuniformly sampled signals in offset linear canonical transform domain. Optics Communications, 2015, 348: 59–65.
  28. S. Xu, L. Feng, Y. Chai, Y. He. Analysis of A-stationary random signals in the linear canonical transform domain. Signal Processing, 2018, 146: 126–132.
  29. S. Xu, L. Feng, Y. Chai, Y. Hu, L. Huang. Sampling theorems for bandlimited random signals in the offset linear canonical transform domain. AEU-International Journal of Electronics and Communications, 2017, 81: 114–119.
  30. S. Xu, L. Huang, Y. Chai, Y. He. Nonuniform sampling theorems for bandlimited signals in the offset linear canonical transform. Circuits Systems and Signal Processing, 2018, 37: 3227–3244.
  31. S. Xu, C. Jiang, Y. Chai, Y. Hu, L. Huang. Nonuniform sampling theorems for random signals in the linear canonical transform domain. International Journal of Electronics, 2018, 105(6): 1051–1062.
  32. L. Xu, F. Zhang, R. Tao. Randomized nonuniform sampling and reconstruction in fractional Fourier domain. Signal Processing, 2016, 120: 311–322.
  33. X. Zhi, D. Wei, W. Zhang. A generalized convolution theorem for the special affine Fourier transform and its application to filtering. Optik, 2016, 127: 2613–2616.
Figure 1. Multiplicative filtering in the SAFT
Figure 1. Multiplicative filtering in the SAFT
Preprints 96227 g001
Figure 2. The nonuniform sampling process
Figure 2. The nonuniform sampling process
Preprints 96227 g002
Figure 3. The equivalent system of the nonuniform sampling, where the filtering through filter h 1 ( t ) means that y ˜ ( t ) = 1 2 π R x ˜ ( s ) h 1 ( t s ) d s
Figure 3. The equivalent system of the nonuniform sampling, where the filtering through filter h 1 ( t ) means that y ˜ ( t ) = 1 2 π R x ˜ ( s ) h 1 ( t s ) d s
Preprints 96227 g003
Figure 4. The approximate reconstruction with sinc interpolation function, where x ¯ ( t ) = T T N n = + x 1 ( t n ) h 2 ( t t ˜ n )
Figure 4. The approximate reconstruction with sinc interpolation function, where x ¯ ( t ) = T T N n = + x 1 ( t n ) h 2 ( t t ˜ n )
Preprints 96227 g004
Figure 5. An equivalent nonuniform sampling and reconstruction system to Figure 4, where v ( t ) is an additive noise with zero mean which is uncorrelated with x ( t ) and has power spectral density as (4.5).
Figure 5. An equivalent nonuniform sampling and reconstruction system to Figure 4, where v ( t ) is an additive noise with zero mean which is uncorrelated with x ( t ) and has power spectral density as (4.5).
Preprints 96227 g005
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.
Copyright: This open access article is published under a Creative Commons CC BY 4.0 license, which permit the free download, distribution, and reuse, provided that the author and preprint are cited in any reuse.
Prerpints.org logo

Preprints.org is a free preprint server supported by MDPI in Basel, Switzerland.

Subscribe

© 2024 MDPI (Basel, Switzerland) unless otherwise stated