# Spectral Analysis of Sampled Signals in the Linear Canonical Transform Domain

Authors: Bing-Zhao Li; Tian-Zhou Xu
Journal: Mathematical Problems in Engineering
 (2012)

---

## Abstract

The spectral analysis of uniform or nonuniform sampling signal is one of the hot topics in digital signal processing community. Theories and applications of uniformly and nonuniformly sampled one-dimensional or two-dimensional signals in the traditional Fourier domain have been well studied. But so far, none of the research papers focusing on the spectral analysis of sampled signals in the linear canonical transform domain have been published. In this paper, we investigate the spectrum of sampled signals in the linear canonical transform domain. Firstly, based on the properties of the spectrum of uniformly sampled signals, the uniform sampling theorem of two dimensional signals has been derived. Secondly, the general spectral representation of periodic nonuniformly sampled one and two dimensional signals has been obtained. Thirdly, detailed analysis of periodic nonuniformly sampled chirp signals in the linear canonical transform domain has been performed.


---

## Body

## 1. Introduction

The sampling process is one of the fundamental concepts of digital signal processing, which serves as a bridge between the continuous physical signals and discrete signals. The sampling process can be marked as uniform sampling and nonuniform sampling according to the sampling offsets. Theories and applications of uniformly and nonuniformly sampled one and two-dimensional signals in traditional Fourier domain have been well studied [1–3]. The spectral analysis of uniform and nonuniform samples has also been investigated in the Fourier domain [4–6]. However, most of the theories and methods derived in the literature are only suitable for bandlimited signals of Fourier transform (FT) domain, they will derive the wrong (or at least suboptimal) conclusions about signals and systems that not bandlimited in the Fourier domain [7]. For example, signals such as gravity waves, broad-band chirp signals, and radar and sonar signals may behave as nonbandlimited in the Fourier domain. It is therefore worthwhile and interesting to explore the new nonstationary signal processing tools and derive novel properties of nonbandlimited signals of the Fourier domain in the new transform domain.With the efforts of the scientists of signal processing community, many useful tools such as the wavelet transform, the fractional Fourier transform [8], the Hilbert-Huang transform [9], and the linear canonical transform (LCT) have been proposed to process nonstationary signals. The LCT, which was introduced during the 1970s with four parameters [10], has been proven to be one of the most powerful tools for nonstationary signal processing. The signal processing transforms, such as the Fourier transform (FT), the fractional Fourier transform (FrFT), the Fresnel transform, and the scaling operations are all special cases of the LCT [8]. The well-known concepts associated with the traditional Fourier transform, such as the uncertainty principles [11], the convolution and product theorem [12, 13], the Hilbert transform [14–16], and the Poisson summation formula [17] are well studied and extended in the LCT domain. The expansion of the classical uniform or nonuniform sampling theorem for the band-limited or time-limited signal in the LCT domain has recently been studied [18–20]. Recently, the authors of [21, 22] address the multichannel sampling problems and their novel results can be looked as the generalization of the well-known sampling theorem for the LCT and FT domain. These sampling theorems establish the fact that a band-limited or time-limited continuous signal in the LCT domain can be completely reconstructed by a set of equidistantly spaced signal samples, but so far none of the research papers covering the spectral analysis of sampled signals in the LCT domain have been published. Therefore, exploring the spectral properties for sampled signals in the LCT domain is worthwhile and interesting.Focusing on the spectral analysis of uniformly or nonuniformly sampled signals in the LCT domain, this paper investigates the spectrum of signals based on the uniform and periodic nonuniform samples in the LCT domain. The paper is organized as follows. The Preliminary is proposed in Section2. In Section 3, the uniform sampling theorem of two-dimensional LCT is derived based on the definition of LCT and the traditional sampling theorem. The spectral analysis of one and two dimensional periodic nonuniform sampling signals is derived in Section 4. Section 5 is the conclusion of the paper.


## 2. The Preliminary

### 2.1. The Linear Canonical Transform (LCT)

The LCT with parameters(abcd) of a signal f(t) is defined as [8]:
(2.1)LfA(u)=LfA[f(t)](u)={1j2πb∫-∞+∞f(t)ej(1/2)[(a/b)t2-(2/b)tu+(d/b)u2]dt,b≠0,dej(1/2)cdu2f(du),b=0,
where
(2.2)A≡(abcd),det⁡(A)=ad-bc=1.
The LCT of two-dimensional signal f(x,y) is defined as [16]:
(2.3)LfA,B{f(x,y)}=LfA,B(u,v)=∬-∞+∞f(x,y)KA,B(x,y,u,v)dxdy={12πj1b1b2⋅ej[(d1u2/2b1)+(d2v2/2b2)]×∬-∞+∞e-j[(ux/b1)+(vy/b2)]⋅ej[(a1x2/2b1)+(a2y2/2b2)]f(x,y)dxdy,b1b2≠0,|A|=|B|=1d1d2⋅ej[(c1d1u2+c2d2v2/2)]f(d1u,d2v),b12+b22=0,
where
(2.4)KA,B(x,y,u,v)=KA(x,u)KB(y,v),KA(x,u)=12πjb1⋅ej(d1u2/2b1)e-j(ux/b1)ej(a1x2/2b1);KB(y,v)=12πjb2⋅ej(d2v2/2b2)e-j(vy/b2)ej(a2y2/2b2),A=(a1b1c1d1),B=(a2b2c2d2).
In this paper we restrict ourselves to the class of the LCT with real parameters, that is, the parameters in matrix A,B are real numbers.


### 2.2. The Sampling Theorem of One-Dimensional Signal in the LCT Domain

It is shown in [18, 19] that if a signal is nonbandlimited in the Fourier domain, then it can be bandlimited in the LCT domain. Therefore, the common sampling theorem in Fourier domain may not suitable for these nonbandlimited signals. In order to solve these problems, the uniform sampling theorems for one-dimensional signal in the LCT domain have been derived as shown below [18, 19].Lemma 2.1.
Assume a signalf(t) bandlimited to ΩA in the LCT domain with parameter A=(abcd) and b>0; then the sampling theorem for signal f(t) can be expressed as:
(2.5)f(t)=e-j(a/2b)t2∑n=-∞+∞f(nT)ej(a/2b)(nT)2sin⁡[ΩA(t-nT)/b]ΩA(t-nT)/b,
where T is the sampling period and satisfies T=πb/ΩA;  and the Nyquist rate of sampling theorem associated with the LCT is fs=ΩA/πb.Proof.
see [18, 19].


### 2.3. The Periodic Nonuniform Sampling Model

The periodic nonuniform sampling [4–6, 21–23], which also named as recurrent nonuniform sampling, arises in a broad range of applications. For example [4], we might consider converting a continuous time signal to a discrete time signal using a series of A/D converters, and each is operating at a rate lower than the Nyquist rate such that the average sampling rate is equal to the Nyquist rate. This may be beneficial in applications where high-rate A/D converters are required. Typically, the cost and complexity of a converter will increase (more than linearly) with the rate. In such cases, we can benefit from converting a continuous time signal to a discrete time signal using N  A/D converters [23], each is operating at one Nth of the Nyquist rate. Since the converters are typically not synchronized, the resulting discrete time signal is a combination of N sequences of uniform samples. Thus, the resulting discrete time signal corresponds to recurrent nonuniform samples of the continuous time signal.In this form of sampling, the sampling points are divided into several groups of points. The groups have a recurrent period, which is denoted byNT and equal to N times the Nyquist period. Each period consists of N nonuniform sampling points. Denoting the points in one period by tn,n=0,1,…,N-1, the complete set of sampling points are tn+kNT, n=0,1,2,…,N-1,k∈Z. Periodic nonuniform samples can be regarded as a combination of N sequences of uniform samples taken at one Nth of the Nyquist rate. An example of this periodic nonuniform sampling distribution is depicted in Figure 1 [5].Figure 1
The recurrent nonuniform sampling model [5].Besides the one-dimensional recurrent nonuniform sampling model, in real applications we often meet the multidimensional nonuniform sampling. Jenq investigates the two-dimensional periodic nonuniform sampling model and obtains a perfect spectral reconstruction method in the Fourier domain [5, 6], Feuer and Coodwin introduce the multi-dimensional recurrent nonuniform sampling model in [4] and derive the reconstruction methods based on the filterbanks. It is shown that a specific example of this situation arises when one utilizes multiple identical digital cameras on the same scene [4].The signal recovery from one or multidimensional recurrent nonuniform sampling points are well studied in the Fourier domain; however, there are no paper published about the spectral analysis from recurrent nonuniform sampling points in the LCT domain. It is worthwhile and interesting to investigate the spectral analysis and reconstruction in the LCT domain.


## 2.1. The Linear Canonical Transform (LCT)

The LCT with parameters(abcd) of a signal f(t) is defined as [8]:
(2.1)LfA(u)=LfA[f(t)](u)={1j2πb∫-∞+∞f(t)ej(1/2)[(a/b)t2-(2/b)tu+(d/b)u2]dt,b≠0,dej(1/2)cdu2f(du),b=0,
where
(2.2)A≡(abcd),det⁡(A)=ad-bc=1.
The LCT of two-dimensional signal f(x,y) is defined as [16]:
(2.3)LfA,B{f(x,y)}=LfA,B(u,v)=∬-∞+∞f(x,y)KA,B(x,y,u,v)dxdy={12πj1b1b2⋅ej[(d1u2/2b1)+(d2v2/2b2)]×∬-∞+∞e-j[(ux/b1)+(vy/b2)]⋅ej[(a1x2/2b1)+(a2y2/2b2)]f(x,y)dxdy,b1b2≠0,|A|=|B|=1d1d2⋅ej[(c1d1u2+c2d2v2/2)]f(d1u,d2v),b12+b22=0,
where
(2.4)KA,B(x,y,u,v)=KA(x,u)KB(y,v),KA(x,u)=12πjb1⋅ej(d1u2/2b1)e-j(ux/b1)ej(a1x2/2b1);KB(y,v)=12πjb2⋅ej(d2v2/2b2)e-j(vy/b2)ej(a2y2/2b2),A=(a1b1c1d1),B=(a2b2c2d2).
In this paper we restrict ourselves to the class of the LCT with real parameters, that is, the parameters in matrix A,B are real numbers.


## 2.2. The Sampling Theorem of One-Dimensional Signal in the LCT Domain

It is shown in [18, 19] that if a signal is nonbandlimited in the Fourier domain, then it can be bandlimited in the LCT domain. Therefore, the common sampling theorem in Fourier domain may not suitable for these nonbandlimited signals. In order to solve these problems, the uniform sampling theorems for one-dimensional signal in the LCT domain have been derived as shown below [18, 19].Lemma 2.1.
Assume a signalf(t) bandlimited to ΩA in the LCT domain with parameter A=(abcd) and b>0; then the sampling theorem for signal f(t) can be expressed as:
(2.5)f(t)=e-j(a/2b)t2∑n=-∞+∞f(nT)ej(a/2b)(nT)2sin⁡[ΩA(t-nT)/b]ΩA(t-nT)/b,
where T is the sampling period and satisfies T=πb/ΩA;  and the Nyquist rate of sampling theorem associated with the LCT is fs=ΩA/πb.Proof.
see [18, 19].


## 2.3. The Periodic Nonuniform Sampling Model

The periodic nonuniform sampling [4–6, 21–23], which also named as recurrent nonuniform sampling, arises in a broad range of applications. For example [4], we might consider converting a continuous time signal to a discrete time signal using a series of A/D converters, and each is operating at a rate lower than the Nyquist rate such that the average sampling rate is equal to the Nyquist rate. This may be beneficial in applications where high-rate A/D converters are required. Typically, the cost and complexity of a converter will increase (more than linearly) with the rate. In such cases, we can benefit from converting a continuous time signal to a discrete time signal using N  A/D converters [23], each is operating at one Nth of the Nyquist rate. Since the converters are typically not synchronized, the resulting discrete time signal is a combination of N sequences of uniform samples. Thus, the resulting discrete time signal corresponds to recurrent nonuniform samples of the continuous time signal.In this form of sampling, the sampling points are divided into several groups of points. The groups have a recurrent period, which is denoted byNT and equal to N times the Nyquist period. Each period consists of N nonuniform sampling points. Denoting the points in one period by tn,n=0,1,…,N-1, the complete set of sampling points are tn+kNT, n=0,1,2,…,N-1,k∈Z. Periodic nonuniform samples can be regarded as a combination of N sequences of uniform samples taken at one Nth of the Nyquist rate. An example of this periodic nonuniform sampling distribution is depicted in Figure 1 [5].Figure 1
The recurrent nonuniform sampling model [5].Besides the one-dimensional recurrent nonuniform sampling model, in real applications we often meet the multidimensional nonuniform sampling. Jenq investigates the two-dimensional periodic nonuniform sampling model and obtains a perfect spectral reconstruction method in the Fourier domain [5, 6], Feuer and Coodwin introduce the multi-dimensional recurrent nonuniform sampling model in [4] and derive the reconstruction methods based on the filterbanks. It is shown that a specific example of this situation arises when one utilizes multiple identical digital cameras on the same scene [4].The signal recovery from one or multidimensional recurrent nonuniform sampling points are well studied in the Fourier domain; however, there are no paper published about the spectral analysis from recurrent nonuniform sampling points in the LCT domain. It is worthwhile and interesting to investigate the spectral analysis and reconstruction in the LCT domain.


## 3. The Uniform Sampling Theorems

### 3.1. The Uniform Sampling Theorem of Two-Dimensional Signals

This section focuses on the uniform sampling theorem for two-dimensional signals in the LCT domain. Letf(x,y) be an analog two-dimensional signal with its continuous LCT LfA,B(u,v) bandlimited to (Ωx,Ωy) in the LCT domain. In other words, LfA,B(u,v)=0, when |u|≥Ωx and |v|≥Ωy.Supposefs(x,y) is obtained by sampling the signal f(x,y) at uniformly spaced grids Tx and Ty, Tx and Ty are uniform sampling period on X and Y-axis, respectively, we have
(3.1)fs(x,y)=∑n=-∞∞∑k=-∞∞f(x,y)δ(x-nTx)δ(y-kTy).
Then, the spectrum of these samples function fs(x,y) in the LCT domain can be shown in the following theorem.Theorem 3.1.
Supposefs(x,y) is obtained by sampling the continuous signal f(x,y) at uniformly spaced grids Tx and Ty, then the spectrum of these uniformly samples can be represented as:
(3.2)DLfA,B(u,v)=1TxTyej[(d1u2/2b1)+(d2v2/2b2)]×∑k=-∞+∞∑l=-∞+∞LfA,B(u-2πkTxb1,v-2πlTyb2)×e-j(d1/2b1)(u-(2πk/Tx)b1)2e-j(d2/2b2)(v-(2πl/Ty)b2)2.Proof.
Applying the LCT to both sides of (3.1), we obtain
(3.3)DLfA,B(u,v)=∬-∞+∞fs(x,y)KA,B(x,y,u,v)dxdy=12πj1b1b2⋅ej[(d1u2/2b1)+(d2v2/2b2)]⋅∬-∞+∞e-j[(ux/b1)+(vy/b2)]ej[(a1x2/2b1)+(a2y2/2b2)]×∑n=-∞+∞∑k=-∞+∞f(x,y)δ(x-nTx)δ(y-kTy)dxdy.
Rearranging (3.3) and using the identity ∑k=-∞+∞δ(x-kTx)∑k=-∞+∞δ(x-kTx)=(1/Tx)∑k=-∞+∞ej(2kπx/Tx)(1/Ty)∑k=-∞+∞ej(2kπy/Ty), (3.3), can be rewritten as
(3.4)DLfA,B=12πj1b1b2⋅ej[(d1u2/2b1)+(d2v2/2b2)]⋅∬-∞+∞f(x,y)e-j[(ux/b1)+(vy/b2)]ej[(a1x2/2b1)+(a2y2/2b2)]×{1Tx∑k=-∞+∞e(j2πkx/Tx)⋅1Ty∑l=-∞+∞e(j2πly/Ty)}dxdy=12πj1TxTy1b1b2⋅ej[(d1u2/2b1)+(d2v2/2b2)]⋅∑k=-∞+∞∑l=-∞+∞∬-∞+∞f(x,y)e-j[((u-(2πk/Tx)b1)x/b1)+(v-(2πk/Ty)b2y/b2)]ej[(a1x2/2b1)+(a2y2/2b2)].
The final result can be obtained by the definition of the LCT.From (3.2), the spectra of two-dimensional uniform sampling signal in the LCT domain can be seen as the repetition of LfA,B(u,v)e-j(d1/2b1)u2e-j(d2/2b2)v2, and these repetitions do not overlap if we chose the sampling interval Tx≤πb1/Ωx,Ty≤πb2/Ωy. In order to single out just one copy, we apply a low-pass filter to DLfA,B(u,v). 								(3.5)H(u,v)={TxTy,|u|≤Ωx,|v|≤Ωy,0,others.
So the original signal f(x,y) can be derived by the inverse LCT transform of DLfA,B(u,v)H(u,v). The reconstruction formula is presented in the following theorem.Theorem 3.2.
supposef(x,y) be an analog two-dimensional signal with its continuous LCT LfA,B(u,v) bandlimited to (Ωx,Ωy) in the LCT domain. Then the original signal f(x,y) can be reconstructed by the following uniform sampling formula
(3.6)f(x,y)=e-j(a1/2b1)x2-j(a2/2b2)y2×∑n=-∞+∞∑k=-∞+∞f(nTx,kTy)ej(a1/2b1)(nTx)2+j(a2/2b2)(kTy)2⋅sin⁡[Ωx(x-nTx)/b1]Ωx(x-nTx)/b1sin⁡[Ωy(y-kTy)/b2]Ωy(y-kTy)/b2,
where Tx,Ty is the sampling period in the X-axis and Y-axis, respectively, and satisfies Tx≤πbx/Ωx,Ty≤πby/Ωy.Proof.
From the definition ofDLfA,B(u,v)H(u,v), the original signal can be derived by the inverse LCT transform of DLfA,B(u,v)H(u,v). In other words, the original signal can be rewritten as:
(3.7)f(x,y)=LA-1,B-1{DLfA,B(u,v)H(u,v)}(x,y)=1TxTy∬-∞+∞KA-1(u,x)KB-1(v,y)×∬-∞+∞f(xn,k,yn,k)KA(xn,k,u)KB(yn,k,v)H(u,v)dxn,kdyn,kdudv=∫-∞+∞f(xn,k,yn,k)×∫-∞+∞KA-1(u,x)KA(xn,k,u)KB-1(v,y)KB(yn,k,v)H(u,v)dudvdxn,kdyn,k=(12π)21b1b2e-(a1/2b1)x2-(a2/2b2)y2×∬-∞+∞f(xn,k,yn,k)e(a1/2b1)xn,k2e(a2/2b2)yn,k2×∬-∞+∞e(1/b1)u(x-xn,k)e(1/b2)v(y-yn,k)H(u,v)dudvdxn,kdyn,k.
Because
(3.8)∬-∞+∞e(1/b1)u(x-xn,k)e(1/b2)v(y-yn,k)H(u,v)dudv=TxTy∫-Ωy+Ωy∫-Ωx+Ωxej(1/b1)u(x-xl)e(1/b2)v(y-yn,k)dudv=TxTy2sin⁡[Ωx(x-xn,k)/b1](x-xn,k)/b12sin⁡[Ωy(y-yn,k)/b2](y-yn,k)/b2.
Substitute (3.8) into (3.7) we obtain
(3.9)f(x,y)=e-(a1/2b1)x2-(a2/2b2)y2×∬-∞+∞f(xn,k,yn,k)e(a1/2b1)xn,k2e(a2/2b2)yn,k2×Txsin⁡[Ωx(x-xn,k)/b1]π(x-xn,k)Tysin⁡[Ωy(y-yn,k)/b2]π(y-yn,k)dxn,kdyn,k.
The final result can be obtained by substituting (3.1) into above equation and letting the sample interval to be Tx=πbx/Ωx,Ty=πby/Ωy.


### 3.2. Special Cases of the Uniform Sampling Theorem of Two-Dimensional Signals

(1) Let A=(01-10) and B=(01-10), then (3.2) and (3.6) reduce to
(3.10)DLfA,B(u,v)=1TxTy∑k=-∞+∞∑l=-∞+∞LfA,B(u-k2πTx,v-l2πTy),(3.11)f(x,y)=∑n=-∞+∞∑k=-∞+∞f(nTx,kTy)sin⁡[Ωx(x-nTx)]Ωx(x-nTx)sin⁡[Ωy(y-kTy)]Ωy(y-kTy),
which are the well known digital spectral representation and the reconstruction formula of uniformly sampled two-dimensional signals in traditional Fourier transform domain.(2) Let A=(cosθ1sinθ1-sinθ1cosθ1) and B=(cosθ2sinθ2-sinθ2cosθ2), then (3.2) and (3.6) reduce to
(3.12)DLfA,B(u,v)=1TxTyej[(cotθ1u2/2)+(cotθ2v2/2)]×∑k=-∞+∞∑l=-∞+∞LfA,B(u-2πkTxsin⁡θ1,v-2πlTysin⁡θ2)×e-j(cotθ1/2)(u-(2πk/Tx)sin⁡θ1)2e-j(cotθ2/2)(v-(2πl/Ty)sin⁡θ2)2,(3.13)f(x,y)=e-j(cotθ1/2)x2-j(cotθ2/2)y2×∑n=-∞+∞∑k=-∞+∞f(nTx,kTy)ej(cotθ1/2)(nTx)2+j(cotθ2/2)(kTy)2⋅sin⁡[cscθ1Ωx(x-nTx)]cscθ1Ωx(x-nTx)sin⁡[cscθ2Ωy(y-kTy)]cscθ2Ωy(y-kTy).
They can be looked as the digital spectrum of a uniformly sampled signal and the reconstruction formula in fractional Fourier transform domain, respectively.


## 3.1. The Uniform Sampling Theorem of Two-Dimensional Signals

This section focuses on the uniform sampling theorem for two-dimensional signals in the LCT domain. Letf(x,y) be an analog two-dimensional signal with its continuous LCT LfA,B(u,v) bandlimited to (Ωx,Ωy) in the LCT domain. In other words, LfA,B(u,v)=0, when |u|≥Ωx and |v|≥Ωy.Supposefs(x,y) is obtained by sampling the signal f(x,y) at uniformly spaced grids Tx and Ty, Tx and Ty are uniform sampling period on X and Y-axis, respectively, we have
(3.1)fs(x,y)=∑n=-∞∞∑k=-∞∞f(x,y)δ(x-nTx)δ(y-kTy).
Then, the spectrum of these samples function fs(x,y) in the LCT domain can be shown in the following theorem.Theorem 3.1.
Supposefs(x,y) is obtained by sampling the continuous signal f(x,y) at uniformly spaced grids Tx and Ty, then the spectrum of these uniformly samples can be represented as:
(3.2)DLfA,B(u,v)=1TxTyej[(d1u2/2b1)+(d2v2/2b2)]×∑k=-∞+∞∑l=-∞+∞LfA,B(u-2πkTxb1,v-2πlTyb2)×e-j(d1/2b1)(u-(2πk/Tx)b1)2e-j(d2/2b2)(v-(2πl/Ty)b2)2.Proof.
Applying the LCT to both sides of (3.1), we obtain
(3.3)DLfA,B(u,v)=∬-∞+∞fs(x,y)KA,B(x,y,u,v)dxdy=12πj1b1b2⋅ej[(d1u2/2b1)+(d2v2/2b2)]⋅∬-∞+∞e-j[(ux/b1)+(vy/b2)]ej[(a1x2/2b1)+(a2y2/2b2)]×∑n=-∞+∞∑k=-∞+∞f(x,y)δ(x-nTx)δ(y-kTy)dxdy.
Rearranging (3.3) and using the identity ∑k=-∞+∞δ(x-kTx)∑k=-∞+∞δ(x-kTx)=(1/Tx)∑k=-∞+∞ej(2kπx/Tx)(1/Ty)∑k=-∞+∞ej(2kπy/Ty), (3.3), can be rewritten as
(3.4)DLfA,B=12πj1b1b2⋅ej[(d1u2/2b1)+(d2v2/2b2)]⋅∬-∞+∞f(x,y)e-j[(ux/b1)+(vy/b2)]ej[(a1x2/2b1)+(a2y2/2b2)]×{1Tx∑k=-∞+∞e(j2πkx/Tx)⋅1Ty∑l=-∞+∞e(j2πly/Ty)}dxdy=12πj1TxTy1b1b2⋅ej[(d1u2/2b1)+(d2v2/2b2)]⋅∑k=-∞+∞∑l=-∞+∞∬-∞+∞f(x,y)e-j[((u-(2πk/Tx)b1)x/b1)+(v-(2πk/Ty)b2y/b2)]ej[(a1x2/2b1)+(a2y2/2b2)].
The final result can be obtained by the definition of the LCT.From (3.2), the spectra of two-dimensional uniform sampling signal in the LCT domain can be seen as the repetition of LfA,B(u,v)e-j(d1/2b1)u2e-j(d2/2b2)v2, and these repetitions do not overlap if we chose the sampling interval Tx≤πb1/Ωx,Ty≤πb2/Ωy. In order to single out just one copy, we apply a low-pass filter to DLfA,B(u,v). 								(3.5)H(u,v)={TxTy,|u|≤Ωx,|v|≤Ωy,0,others.
So the original signal f(x,y) can be derived by the inverse LCT transform of DLfA,B(u,v)H(u,v). The reconstruction formula is presented in the following theorem.Theorem 3.2.
supposef(x,y) be an analog two-dimensional signal with its continuous LCT LfA,B(u,v) bandlimited to (Ωx,Ωy) in the LCT domain. Then the original signal f(x,y) can be reconstructed by the following uniform sampling formula
(3.6)f(x,y)=e-j(a1/2b1)x2-j(a2/2b2)y2×∑n=-∞+∞∑k=-∞+∞f(nTx,kTy)ej(a1/2b1)(nTx)2+j(a2/2b2)(kTy)2⋅sin⁡[Ωx(x-nTx)/b1]Ωx(x-nTx)/b1sin⁡[Ωy(y-kTy)/b2]Ωy(y-kTy)/b2,
where Tx,Ty is the sampling period in the X-axis and Y-axis, respectively, and satisfies Tx≤πbx/Ωx,Ty≤πby/Ωy.Proof.
From the definition ofDLfA,B(u,v)H(u,v), the original signal can be derived by the inverse LCT transform of DLfA,B(u,v)H(u,v). In other words, the original signal can be rewritten as:
(3.7)f(x,y)=LA-1,B-1{DLfA,B(u,v)H(u,v)}(x,y)=1TxTy∬-∞+∞KA-1(u,x)KB-1(v,y)×∬-∞+∞f(xn,k,yn,k)KA(xn,k,u)KB(yn,k,v)H(u,v)dxn,kdyn,kdudv=∫-∞+∞f(xn,k,yn,k)×∫-∞+∞KA-1(u,x)KA(xn,k,u)KB-1(v,y)KB(yn,k,v)H(u,v)dudvdxn,kdyn,k=(12π)21b1b2e-(a1/2b1)x2-(a2/2b2)y2×∬-∞+∞f(xn,k,yn,k)e(a1/2b1)xn,k2e(a2/2b2)yn,k2×∬-∞+∞e(1/b1)u(x-xn,k)e(1/b2)v(y-yn,k)H(u,v)dudvdxn,kdyn,k.
Because
(3.8)∬-∞+∞e(1/b1)u(x-xn,k)e(1/b2)v(y-yn,k)H(u,v)dudv=TxTy∫-Ωy+Ωy∫-Ωx+Ωxej(1/b1)u(x-xl)e(1/b2)v(y-yn,k)dudv=TxTy2sin⁡[Ωx(x-xn,k)/b1](x-xn,k)/b12sin⁡[Ωy(y-yn,k)/b2](y-yn,k)/b2.
Substitute (3.8) into (3.7) we obtain
(3.9)f(x,y)=e-(a1/2b1)x2-(a2/2b2)y2×∬-∞+∞f(xn,k,yn,k)e(a1/2b1)xn,k2e(a2/2b2)yn,k2×Txsin⁡[Ωx(x-xn,k)/b1]π(x-xn,k)Tysin⁡[Ωy(y-yn,k)/b2]π(y-yn,k)dxn,kdyn,k.
The final result can be obtained by substituting (3.1) into above equation and letting the sample interval to be Tx=πbx/Ωx,Ty=πby/Ωy.


## 3.2. Special Cases of the Uniform Sampling Theorem of Two-Dimensional Signals

(1) Let A=(01-10) and B=(01-10), then (3.2) and (3.6) reduce to
(3.10)DLfA,B(u,v)=1TxTy∑k=-∞+∞∑l=-∞+∞LfA,B(u-k2πTx,v-l2πTy),(3.11)f(x,y)=∑n=-∞+∞∑k=-∞+∞f(nTx,kTy)sin⁡[Ωx(x-nTx)]Ωx(x-nTx)sin⁡[Ωy(y-kTy)]Ωy(y-kTy),
which are the well known digital spectral representation and the reconstruction formula of uniformly sampled two-dimensional signals in traditional Fourier transform domain.(2) Let A=(cosθ1sinθ1-sinθ1cosθ1) and B=(cosθ2sinθ2-sinθ2cosθ2), then (3.2) and (3.6) reduce to
(3.12)DLfA,B(u,v)=1TxTyej[(cotθ1u2/2)+(cotθ2v2/2)]×∑k=-∞+∞∑l=-∞+∞LfA,B(u-2πkTxsin⁡θ1,v-2πlTysin⁡θ2)×e-j(cotθ1/2)(u-(2πk/Tx)sin⁡θ1)2e-j(cotθ2/2)(v-(2πl/Ty)sin⁡θ2)2,(3.13)f(x,y)=e-j(cotθ1/2)x2-j(cotθ2/2)y2×∑n=-∞+∞∑k=-∞+∞f(nTx,kTy)ej(cotθ1/2)(nTx)2+j(cotθ2/2)(kTy)2⋅sin⁡[cscθ1Ωx(x-nTx)]cscθ1Ωx(x-nTx)sin⁡[cscθ2Ωy(y-kTy)]cscθ2Ωy(y-kTy).
They can be looked as the digital spectrum of a uniformly sampled signal and the reconstruction formula in fractional Fourier transform domain, respectively.


## 4. The Spectral Analysis of Nonuniformly Sampled Signals

It is shown in the preliminary Section that the recurrent nonuniform sampling occurs frequently in real applications, and there are many published works on the spectral analysis methods in the Fourier domain [4–6]. However, these results are only suitable for the bandlimited signals in the Fourier domain and may obtain incorrect results for nonbanlimited signals of Fourier domain [19]. Because the nonbandlimited signals can be bandlimited in the LCT domain [8], so it is worthwhile and interesting to investigate the spectral analysis problems in the LCT domain.


### 4.1. Spectral Analysis of One-Dimensional Signal from Recurrent Nonuniform Samples

In above-mentioned, one-dimensional periodic nonuniform sampling model, themth sampling sequence can be seen as uniformly sampling the signal with sampling period MT:
(4.1)sm=[f(tm),f(tM+m),f(t2M+m),…],
the total sampling sequence can be represented as:
(4.2)s={f(tkm)∣tkm=kMT+tm,m=0,1,2,…,M-1,k∈Z}.
Let s̅m=[f(tm),0,…,(M-1)zeros,f(tM+m),0,0,…],m=0,1,…,M-1 and shift s̅mmT positions to the right, for m=0,1,…,M-1 to obtain
(4.3)s̅mz-m=[(mzeros)f(tm),0,…,(M-1)zeros,f(tM+m),0,0,…],
where z-1 is the unit delay operator. Finally, summing up all the subsequences to obtain the original sequence s=∑m=0M-1s̅mz-m, therefore the discrete LCT of the periodic nonuniformly sampled signal sequence s can be derived by the summation of the LCT of M subsequence s̅mz-m(m=0toM-1). The result can be represented as the following theorem.Theorem 4.1.
Letf(t) be an analog signal with its CLCT LfA(u) bandlimited to (-Ωα,Ωα) in the LCT domain; f(t) is sampled nonuniformly with recurrent period MT. Then, the digital spectrum of these periodic nonuniformly sampled points can be represented as:
(4.4)DLfA[f(t)](u)=1MT×∑n=-∞∞∑m=0M-1LfA[u-2πnNTb+armT]⋅ej(1/2)((a-a2d)/b)T2rm2ej(2πn/MT)d(u-(πn/MT)b)⋅e-j[((1-ad)/b)u-(1-ad)(2πn/MT)]rmTe-j(2πnm/M),
where T is the nominal uniform sampling period, and rm=(mT-tm)/T,m=0,1,…,M-1 are sampling time offsets.Proof.
The result can be derived by using the similar methods as in [5] and the properties of LCT.


### 4.2. The Spectral Analysis of Two-Dimensional Signal

Firstly, consider the functionfs(x,y) which are obtained by sampling a two-dimensional signal f(x,y) at nonuniformly spaced grids. The sampling offsets on both X-axis and Y-axis are periodic and have the period M and N, respectively. That means
(4.5)xu,v=uTx+Δxu,v,yu,v=vTy+Δyu,v,
where Tx and Ty are nominal uniform sampling periods on the X-axis and Y-axis, respectively. Δxu,v and Δyu,v are periodic sequence with period M on u and N on v. Because of the periodic structure of Δxu,v and Δyu,v, we can let u=kM+m,v=lN+n, where -∞<k,l<∞, 0≤m≤M-1 and 0≤n≤N-1. Then
(4.6)xu,v=uTx+Δxu,v=(kM+m)Tx+ΔxkM+m,lN+n=kMTx+mTx+Δxm,n=kMTx+mTx+rm,nTx,yu,v=vTy+Δyu,v=(lN+n)Ty+ΔykM+m,lN+n=lNTy+nTy+Δym,n=lNTy+nTy+sm,nTy,
where rm,n=Δxm,n/Tx and sm,n=Δym,n/Ty are the ratios of sampling offsets to the nominal uniform sampling periods. The sampling signal fs(x,y) can be represented as:
(4.7)fs(x,y)=f(x,y)×δ[x-(Mk+m+rm,n)Tx]×δ[y-(Nl+n+sm,n)Ty].
The DLCT of this two-dimensional signal f(xu,v,yu,v) can be obtained as:
(4.8)DLfA,B(u′,v′)=∑u=-∞+∞∑v=-∞+∞f(xu,v,yu,v)×12πj1b1b2ej[(d1u′2/2b1)+(d2v′2/2b2)]e-j[u′uTx/b1+v′vTy/b2]×ej[a1(uTx)2/2b1+a2(vTy)2/2b2],
where b1b2≠0,|A|=|B|=1.Based on the properties of the DLCT and (4.8), the digital spectrum of the nonuniformly sampled two-dimensional signals can be represented in the following Theorem 4.2.Theorem 4.2.
Letf(x,y) be an analog two-dimensional signal with its continuous LCT bandlimited in the LCT domain, f(x,y) is sampled periodically with period M and N in the x and y directions. Then the digital spectrum of these periodic nonuniformly sampled points can be presented as:
(4.9)DLfA,B(u′,v′)=1TxTyej[(d1u′2/2b1)+(d2v′2/2b2)]×∑k=-∞+∞∑l=-∞+∞H(k,l,u′,v′)×e-j[(d1/2b1)(u′+a1rm,nTx-(2πk/MTx)b1)2+(d2/2b2)(v′+a2sm,nTy-(2πl/NTy)b2)2]⋅LfA,B(u′+a1rm,nTx-2πkMTxb1,v′+a2sm,nTy-2πlNTyb2),
where
(4.10)H(k,l,u′,v′)=(1MN)×∑m=0M-1∑n=0N-1ej[(u′/b1+a1rm,nTx/2b1)rm,nTx-(2πk(m+rm,n)/M)]×ej[(v′/b2+a2sm,nTy/2b2)sm,nTy-(2πl(n+sm,n)/N)].Proof.
Following the definition of the DLCT, the digital spectrum of the samples can be written as:(4.11)DLfA,B(u′,v′)=∑u=-∞+∞∑v=-∞+∞fs(x,y)12πj1b1b2×ej[(d1u′2/2b1)+(d2v′2/2b2)]×e-j[u′uTx/b1+v′vTy/b2]×ej[a1(uTx)2/2b1+a2(vTy)2/2b2]=∑k=-∞+∞∑l=-∞+∞∑m=0M-1∑n=0N-112πj1b1b2⋅ej[(d1u′2/2b1)+(d2v′2/2b2)]⋅e-j[u′(kM+m)Tx/b1+v′(lN+n)Ty/b2]⋅ej[a1(kM+m)2Tx2/2b1+a2(lN+n)2Ty2/2b2]×∬-∞+∞f(x,y)×δ[x-(Mk+m+rm,n)Tx]×δ[y-(Nl+n+sm,n)Ty]dxdy.
By the Fourier series of δ[x-(Mk+m+rm,n)Tx] and δ[y-(Nl+n+sm,n)Ty], (4.11) can be rearranged as:
(4.12)DLfA,B(u′,v′)=12πj1b1b2⋅ej[(d1u′2/2b1)+(d2v′2/2b2)]⋅∑m=0M-1∑n=0N-1∬-∞+∞f(x,y)×ej[u′rm,nTx/b1+v′sm,nTy/b2]⋅ej[a1(rm,nTx)2/2b1+a2(sm,nTy)2/2b2]×e-j[u′x/b1+v′y/b2]ej[a1(x2-2rm,nTxx)/2b1+a2(y2-2sm,nTyy)/2b2]×{1MTx∑k=-∞+∞ej2πk[x-(m+rm,n)Tx]/MTx×1NTy∑l=-∞+∞ej2πl[y-(n+sm,n)Ty]/NTy}dxdy.
The final result can be obtained by rearranging (4.12) and following the definition of the LCT.


### 4.3. Special Cases of the Derived Results

Before introducing the new results, let us investigate the special cases of the derived results.(1) When the parameters of LCT reduce toA=(01-10) and B=(01-10), (4.9) can be reduced to
(4.13)DLfA,B(u′,v′)=1TxTy∑k=-∞+∞∑l=-∞+∞H(k,l,u′,v′)LfA,B(u′-k2πMTx,v′-l2πNTy),
where
(4.14)H(k,l,u′,v′)=(1MN)∑m=0M-1∑n=0N-1ej[(u′+(k2π/MTx))rm,nTx-mk(2π/M)]×ej[(v′+(l2π/NTy))sm,nTy-nl(2π/N)],
which is the well known digital spectrum of periodic nonuniformly sampled two-dimensional signal in Fourier transform domain [6].(2) Furthermore, whenrm,n and sm,n are zeros, In this case, the periodic nonuniform sampling model reduces to uniform sampling model. It is easy to verify that H(k,l,u′,v′)=1 when k is an integer multiple of M and l is an integer multiple of N; otherwise H(k,l,u′,v′)=0. Therefore,
(4.15)DLfA,B(u,v)=1TxTyej[(d1u2/2b1)+(d2v2/2b2)]×∑k=-∞+∞∑l=-∞+∞LfA,B(u-2πkTxb1,v-2πlTyb2)e-j(d1/2b1)(u-(2πk/Tx)b1)2×e-j(d2/2b2)(v-(2πk/Tx)b2)2,
which is the digital spectrum of a uniformly sampled two-dimensional signal in linear canonical transform domain derived in Section 3 (Theorem 3.1).(3) When the parameter of LCT reduces toA=(cosθ1sinθ1-sinθ1cosθ1) and B=(cosθ2sinθ2-sinθ2cosθ2), then (4.6) reduces to
(4.16)DLfA,B(u′,v′)=1TxTyej[(d1u′2/2b1)+(d2v′2/2b2)]×∑k=-∞+∞∑l=-∞+∞H(k,l,u′,v′)×e-j[(d1/2b1)(u′+a1rm,nTx-(2πk/MTx)b1)2+(d2/2b2)(v′+a2sm,nTy-(2πl/NTy)b2)2]×LfA,B(u′+a1rm,nTx-2πkMTxb1,v′+a2sm,nTy-2πlNTyb2),
where
(4.17)H(k,l,u′,v′)=(1MN)×∑m=0M-1∑n=0N-1ej[(u′cscθ1+cotθ1rm,nTx/2+(k2π/MTx))rm,nTx-mk(2π/M)]×ej[(v′cscθ2+cotθ2sm,nTy/2+(l2π/NTy))sm,nTy-nl(2π/N)],
this result can be looked as the digital spectral of two-dimensional signals in the fractional Fourier domain.(4) For one-dimensional case, when the parameter of LCT reduces toA=(cosθsinθ-sinθcosθ), the results of Theorem 4.1 reduce to the results of fractional Fourier transform domain of [24], and reduce to the results of Fourier domain of [5] when θ=π/2.From the above-mentioned special cases of the Theorem4.2, we conclude that the formula we derived is more general and more suitable for nonstationary signal processing. Because of chirp signal is one of the mostly used nonstationary signal in the Radar signal processing society. Therefore, the properties and potential applications of this result for chirp signal processing are investigated in the following Sections.


### 4.4. The Spectrum of Chirp Signals in the LCT Domain

The chirp signal can be looked as one of the typical nonstationary signals in modern signal processing community. It is shown that the LCT is one of the most important tools for chirp signal processing. It is therefore worthwhile and interesting to investigate the spectrum of chirp signal in the LCT domain. It is easy to derive that when(a1/b1)+m1=0,(a2/b2)+m2=0, the continuous LCT of a two-dimensional chirp signal f(x,y)=Cej(ω1x+(1/2)m1x2)ej(ω2y+(1/2)m2y2) can be obtained as:
(4.18)LfA,B(u,v)=C2πj1b1b2ej[(d1u2/2b1)+d2v2/2b2]⋅δ(ω1-ub1)⋅δ(ω2-vb2).
Substitute (4.18) in (4.9), the digital spectrum of two-dimensional chirp signal can be represented as following theorem.Theorem 4.3.
Letf(x,y)=Cej(ω1x+(1/2)m1x2)ej(ω2y+(1/2)m2y2), then the digital spectrum from the periodic nonuniformly sampled signal points with period M and N in the x and y directions, respectively, can be represented by the following:
(4.19)DLfA,B(u′,v′)=C2πj1b1b21TxTyej[(d1u′2/2b1)+(d2v′2/2b2)]∑k=-∞+∞∑l=-∞+∞H(k,l)⋅δ[ω1-(u′b1+a1b1rm,nTx-2πkMTx)]δ[ω2-(v′b2+a2b2sm,nTy-2πlNTy)],
where
(4.20)H(k,l)=1MN×∑m=0M-1∑n=0N-1ej[(ω1-(a1/2b1)rm,nTx+(4πk/MTx))rm,nTx-mk(2π/M)]×ej[(ω2-(a2/2b2)sm,nTy+(4πl/MTy))sm,nTy-nl(2π/N)].Proof.
This formula can be derived by substituting (4.18) in (4.9), and applying the properties of the discrete Fourier transform.This theorem indicates that the LCT of a two-dimensional chirp signal from periodic nonuniform sampling points contains line spectra components in the LCT domain, with the magnitude at(k,l) being H(k,l). The difference between the uniform and periodic nonuniform sampling points is that for periodic nonuniform sampling signals, it will produce some additional pseudo spectra in the same LCT domain besides the main spectral component.


### 4.5. The Simulation Results

In order to show the importance of spectrum analysis of the sampling points associated with the LCT, we choose the recurrent nonuniform sampling model as an example to perform the simulation in the LCT domain.Firstly, we choose the chirp signal asf(t)=exp{j(5t-6t2)}, the sampling period T=0.2s in a period [-10,10]. The spectral of uniformly sampled signal is plotted in Figure 2, it is shown that the uniform sampling in the time domain will produce the periodic spectrum in the LCT domain as reflected in Lemma 2.1. Figures 3, 4, and 5 plot the spectrum of nonuniformly sampling points for M = 3, 4, and 5, respectively. From these simulation results, we can see that it will produce additional M-1 pseudospectra in the same LCT domain besides the main spectral component. This was coinciding with the derived results shown in Theorem 4.1.Figure 2
The spectrum of uniform sampling points.Figure 3
The spectrum of nonuniform sampling points withM=3.Figure 4
The spectrum of nonuniform sampling points withM=4.Figure 5
The spectrum of nonuniform sampling points withM=5.Secondly, we choose the two-dimensional chirp signal asf(x,y)=exp{j(5x-6x2)}exp{j(y-y2)}, the sampling period of x and y is equal and choose to be T=0.2s in a period [-10,10]. The spectral of uniformly sampled signal is plotted in Figure 6, it is shown in the simulation that the spectra of two-dimensional uniform sampling signal in the LCT domain are periodic as reflected in Theorem 3.1 of the paper. The spectral of periodic nonuniformly sampled signal with N=4, M=3 is plotted in Figure 7.Figure 6
The spectrum of uniform sampling points.Figure 7
The spectral of nonuniformly sampled signal withN=4, M=3.From Figures6 and 7, we can conclude that the LCT spectrum of periodic nonuniform sampling points for two-dimensional chirp signals will produce the pseudospectra in the same LCT domain besides the main spectral component, as reflected in Theorem 4.3 of the paper.


## 4.1. Spectral Analysis of One-Dimensional Signal from Recurrent Nonuniform Samples

In above-mentioned, one-dimensional periodic nonuniform sampling model, themth sampling sequence can be seen as uniformly sampling the signal with sampling period MT:
(4.1)sm=[f(tm),f(tM+m),f(t2M+m),…],
the total sampling sequence can be represented as:
(4.2)s={f(tkm)∣tkm=kMT+tm,m=0,1,2,…,M-1,k∈Z}.
Let s̅m=[f(tm),0,…,(M-1)zeros,f(tM+m),0,0,…],m=0,1,…,M-1 and shift s̅mmT positions to the right, for m=0,1,…,M-1 to obtain
(4.3)s̅mz-m=[(mzeros)f(tm),0,…,(M-1)zeros,f(tM+m),0,0,…],
where z-1 is the unit delay operator. Finally, summing up all the subsequences to obtain the original sequence s=∑m=0M-1s̅mz-m, therefore the discrete LCT of the periodic nonuniformly sampled signal sequence s can be derived by the summation of the LCT of M subsequence s̅mz-m(m=0toM-1). The result can be represented as the following theorem.Theorem 4.1.
Letf(t) be an analog signal with its CLCT LfA(u) bandlimited to (-Ωα,Ωα) in the LCT domain; f(t) is sampled nonuniformly with recurrent period MT. Then, the digital spectrum of these periodic nonuniformly sampled points can be represented as:
(4.4)DLfA[f(t)](u)=1MT×∑n=-∞∞∑m=0M-1LfA[u-2πnNTb+armT]⋅ej(1/2)((a-a2d)/b)T2rm2ej(2πn/MT)d(u-(πn/MT)b)⋅e-j[((1-ad)/b)u-(1-ad)(2πn/MT)]rmTe-j(2πnm/M),
where T is the nominal uniform sampling period, and rm=(mT-tm)/T,m=0,1,…,M-1 are sampling time offsets.Proof.
The result can be derived by using the similar methods as in [5] and the properties of LCT.


## 4.2. The Spectral Analysis of Two-Dimensional Signal

Firstly, consider the functionfs(x,y) which are obtained by sampling a two-dimensional signal f(x,y) at nonuniformly spaced grids. The sampling offsets on both X-axis and Y-axis are periodic and have the period M and N, respectively. That means
(4.5)xu,v=uTx+Δxu,v,yu,v=vTy+Δyu,v,
where Tx and Ty are nominal uniform sampling periods on the X-axis and Y-axis, respectively. Δxu,v and Δyu,v are periodic sequence with period M on u and N on v. Because of the periodic structure of Δxu,v and Δyu,v, we can let u=kM+m,v=lN+n, where -∞<k,l<∞, 0≤m≤M-1 and 0≤n≤N-1. Then
(4.6)xu,v=uTx+Δxu,v=(kM+m)Tx+ΔxkM+m,lN+n=kMTx+mTx+Δxm,n=kMTx+mTx+rm,nTx,yu,v=vTy+Δyu,v=(lN+n)Ty+ΔykM+m,lN+n=lNTy+nTy+Δym,n=lNTy+nTy+sm,nTy,
where rm,n=Δxm,n/Tx and sm,n=Δym,n/Ty are the ratios of sampling offsets to the nominal uniform sampling periods. The sampling signal fs(x,y) can be represented as:
(4.7)fs(x,y)=f(x,y)×δ[x-(Mk+m+rm,n)Tx]×δ[y-(Nl+n+sm,n)Ty].
The DLCT of this two-dimensional signal f(xu,v,yu,v) can be obtained as:
(4.8)DLfA,B(u′,v′)=∑u=-∞+∞∑v=-∞+∞f(xu,v,yu,v)×12πj1b1b2ej[(d1u′2/2b1)+(d2v′2/2b2)]e-j[u′uTx/b1+v′vTy/b2]×ej[a1(uTx)2/2b1+a2(vTy)2/2b2],
where b1b2≠0,|A|=|B|=1.Based on the properties of the DLCT and (4.8), the digital spectrum of the nonuniformly sampled two-dimensional signals can be represented in the following Theorem 4.2.Theorem 4.2.
Letf(x,y) be an analog two-dimensional signal with its continuous LCT bandlimited in the LCT domain, f(x,y) is sampled periodically with period M and N in the x and y directions. Then the digital spectrum of these periodic nonuniformly sampled points can be presented as:
(4.9)DLfA,B(u′,v′)=1TxTyej[(d1u′2/2b1)+(d2v′2/2b2)]×∑k=-∞+∞∑l=-∞+∞H(k,l,u′,v′)×e-j[(d1/2b1)(u′+a1rm,nTx-(2πk/MTx)b1)2+(d2/2b2)(v′+a2sm,nTy-(2πl/NTy)b2)2]⋅LfA,B(u′+a1rm,nTx-2πkMTxb1,v′+a2sm,nTy-2πlNTyb2),
where
(4.10)H(k,l,u′,v′)=(1MN)×∑m=0M-1∑n=0N-1ej[(u′/b1+a1rm,nTx/2b1)rm,nTx-(2πk(m+rm,n)/M)]×ej[(v′/b2+a2sm,nTy/2b2)sm,nTy-(2πl(n+sm,n)/N)].Proof.
Following the definition of the DLCT, the digital spectrum of the samples can be written as:(4.11)DLfA,B(u′,v′)=∑u=-∞+∞∑v=-∞+∞fs(x,y)12πj1b1b2×ej[(d1u′2/2b1)+(d2v′2/2b2)]×e-j[u′uTx/b1+v′vTy/b2]×ej[a1(uTx)2/2b1+a2(vTy)2/2b2]=∑k=-∞+∞∑l=-∞+∞∑m=0M-1∑n=0N-112πj1b1b2⋅ej[(d1u′2/2b1)+(d2v′2/2b2)]⋅e-j[u′(kM+m)Tx/b1+v′(lN+n)Ty/b2]⋅ej[a1(kM+m)2Tx2/2b1+a2(lN+n)2Ty2/2b2]×∬-∞+∞f(x,y)×δ[x-(Mk+m+rm,n)Tx]×δ[y-(Nl+n+sm,n)Ty]dxdy.
By the Fourier series of δ[x-(Mk+m+rm,n)Tx] and δ[y-(Nl+n+sm,n)Ty], (4.11) can be rearranged as:
(4.12)DLfA,B(u′,v′)=12πj1b1b2⋅ej[(d1u′2/2b1)+(d2v′2/2b2)]⋅∑m=0M-1∑n=0N-1∬-∞+∞f(x,y)×ej[u′rm,nTx/b1+v′sm,nTy/b2]⋅ej[a1(rm,nTx)2/2b1+a2(sm,nTy)2/2b2]×e-j[u′x/b1+v′y/b2]ej[a1(x2-2rm,nTxx)/2b1+a2(y2-2sm,nTyy)/2b2]×{1MTx∑k=-∞+∞ej2πk[x-(m+rm,n)Tx]/MTx×1NTy∑l=-∞+∞ej2πl[y-(n+sm,n)Ty]/NTy}dxdy.
The final result can be obtained by rearranging (4.12) and following the definition of the LCT.


## 4.3. Special Cases of the Derived Results

Before introducing the new results, let us investigate the special cases of the derived results.(1) When the parameters of LCT reduce toA=(01-10) and B=(01-10), (4.9) can be reduced to
(4.13)DLfA,B(u′,v′)=1TxTy∑k=-∞+∞∑l=-∞+∞H(k,l,u′,v′)LfA,B(u′-k2πMTx,v′-l2πNTy),
where
(4.14)H(k,l,u′,v′)=(1MN)∑m=0M-1∑n=0N-1ej[(u′+(k2π/MTx))rm,nTx-mk(2π/M)]×ej[(v′+(l2π/NTy))sm,nTy-nl(2π/N)],
which is the well known digital spectrum of periodic nonuniformly sampled two-dimensional signal in Fourier transform domain [6].(2) Furthermore, whenrm,n and sm,n are zeros, In this case, the periodic nonuniform sampling model reduces to uniform sampling model. It is easy to verify that H(k,l,u′,v′)=1 when k is an integer multiple of M and l is an integer multiple of N; otherwise H(k,l,u′,v′)=0. Therefore,
(4.15)DLfA,B(u,v)=1TxTyej[(d1u2/2b1)+(d2v2/2b2)]×∑k=-∞+∞∑l=-∞+∞LfA,B(u-2πkTxb1,v-2πlTyb2)e-j(d1/2b1)(u-(2πk/Tx)b1)2×e-j(d2/2b2)(v-(2πk/Tx)b2)2,
which is the digital spectrum of a uniformly sampled two-dimensional signal in linear canonical transform domain derived in Section 3 (Theorem 3.1).(3) When the parameter of LCT reduces toA=(cosθ1sinθ1-sinθ1cosθ1) and B=(cosθ2sinθ2-sinθ2cosθ2), then (4.6) reduces to
(4.16)DLfA,B(u′,v′)=1TxTyej[(d1u′2/2b1)+(d2v′2/2b2)]×∑k=-∞+∞∑l=-∞+∞H(k,l,u′,v′)×e-j[(d1/2b1)(u′+a1rm,nTx-(2πk/MTx)b1)2+(d2/2b2)(v′+a2sm,nTy-(2πl/NTy)b2)2]×LfA,B(u′+a1rm,nTx-2πkMTxb1,v′+a2sm,nTy-2πlNTyb2),
where
(4.17)H(k,l,u′,v′)=(1MN)×∑m=0M-1∑n=0N-1ej[(u′cscθ1+cotθ1rm,nTx/2+(k2π/MTx))rm,nTx-mk(2π/M)]×ej[(v′cscθ2+cotθ2sm,nTy/2+(l2π/NTy))sm,nTy-nl(2π/N)],
this result can be looked as the digital spectral of two-dimensional signals in the fractional Fourier domain.(4) For one-dimensional case, when the parameter of LCT reduces toA=(cosθsinθ-sinθcosθ), the results of Theorem 4.1 reduce to the results of fractional Fourier transform domain of [24], and reduce to the results of Fourier domain of [5] when θ=π/2.From the above-mentioned special cases of the Theorem4.2, we conclude that the formula we derived is more general and more suitable for nonstationary signal processing. Because of chirp signal is one of the mostly used nonstationary signal in the Radar signal processing society. Therefore, the properties and potential applications of this result for chirp signal processing are investigated in the following Sections.


## 4.4. The Spectrum of Chirp Signals in the LCT Domain

The chirp signal can be looked as one of the typical nonstationary signals in modern signal processing community. It is shown that the LCT is one of the most important tools for chirp signal processing. It is therefore worthwhile and interesting to investigate the spectrum of chirp signal in the LCT domain. It is easy to derive that when(a1/b1)+m1=0,(a2/b2)+m2=0, the continuous LCT of a two-dimensional chirp signal f(x,y)=Cej(ω1x+(1/2)m1x2)ej(ω2y+(1/2)m2y2) can be obtained as:
(4.18)LfA,B(u,v)=C2πj1b1b2ej[(d1u2/2b1)+d2v2/2b2]⋅δ(ω1-ub1)⋅δ(ω2-vb2).
Substitute (4.18) in (4.9), the digital spectrum of two-dimensional chirp signal can be represented as following theorem.Theorem 4.3.
Letf(x,y)=Cej(ω1x+(1/2)m1x2)ej(ω2y+(1/2)m2y2), then the digital spectrum from the periodic nonuniformly sampled signal points with period M and N in the x and y directions, respectively, can be represented by the following:
(4.19)DLfA,B(u′,v′)=C2πj1b1b21TxTyej[(d1u′2/2b1)+(d2v′2/2b2)]∑k=-∞+∞∑l=-∞+∞H(k,l)⋅δ[ω1-(u′b1+a1b1rm,nTx-2πkMTx)]δ[ω2-(v′b2+a2b2sm,nTy-2πlNTy)],
where
(4.20)H(k,l)=1MN×∑m=0M-1∑n=0N-1ej[(ω1-(a1/2b1)rm,nTx+(4πk/MTx))rm,nTx-mk(2π/M)]×ej[(ω2-(a2/2b2)sm,nTy+(4πl/MTy))sm,nTy-nl(2π/N)].Proof.
This formula can be derived by substituting (4.18) in (4.9), and applying the properties of the discrete Fourier transform.This theorem indicates that the LCT of a two-dimensional chirp signal from periodic nonuniform sampling points contains line spectra components in the LCT domain, with the magnitude at(k,l) being H(k,l). The difference between the uniform and periodic nonuniform sampling points is that for periodic nonuniform sampling signals, it will produce some additional pseudo spectra in the same LCT domain besides the main spectral component.


## 4.5. The Simulation Results

In order to show the importance of spectrum analysis of the sampling points associated with the LCT, we choose the recurrent nonuniform sampling model as an example to perform the simulation in the LCT domain.Firstly, we choose the chirp signal asf(t)=exp{j(5t-6t2)}, the sampling period T=0.2s in a period [-10,10]. The spectral of uniformly sampled signal is plotted in Figure 2, it is shown that the uniform sampling in the time domain will produce the periodic spectrum in the LCT domain as reflected in Lemma 2.1. Figures 3, 4, and 5 plot the spectrum of nonuniformly sampling points for M = 3, 4, and 5, respectively. From these simulation results, we can see that it will produce additional M-1 pseudospectra in the same LCT domain besides the main spectral component. This was coinciding with the derived results shown in Theorem 4.1.Figure 2
The spectrum of uniform sampling points.Figure 3
The spectrum of nonuniform sampling points withM=3.Figure 4
The spectrum of nonuniform sampling points withM=4.Figure 5
The spectrum of nonuniform sampling points withM=5.Secondly, we choose the two-dimensional chirp signal asf(x,y)=exp{j(5x-6x2)}exp{j(y-y2)}, the sampling period of x and y is equal and choose to be T=0.2s in a period [-10,10]. The spectral of uniformly sampled signal is plotted in Figure 6, it is shown in the simulation that the spectra of two-dimensional uniform sampling signal in the LCT domain are periodic as reflected in Theorem 3.1 of the paper. The spectral of periodic nonuniformly sampled signal with N=4, M=3 is plotted in Figure 7.Figure 6
The spectrum of uniform sampling points.Figure 7
The spectral of nonuniformly sampled signal withN=4, M=3.From Figures6 and 7, we can conclude that the LCT spectrum of periodic nonuniform sampling points for two-dimensional chirp signals will produce the pseudospectra in the same LCT domain besides the main spectral component, as reflected in Theorem 4.3 of the paper.


## 5. Conclusion

Based on the spectral properties of uniform and nonuniform samples in the LCT domain, this paper investigates the spectrum of one and two dimensional periodic nonuniformly sampled signals in the LCT domain. The digital spectral representation of one- and two-dimensional periodic nonuniformly sampled signals in the LCT domain has been obtained, and a detailed analysis for chirp signals has been carried out. The derived results can be looked as the generalization of the traditional results in the Fourier domain and fractional Fourier transform domain and more suitable for the nonstationary signal analysis and processing. The future work direction will be the spectrum analysis ofNth order periodic nonuniform sampling [23] points and the spectrum reconstruction from the periodic or Nth order periodic nonuniform sampling signals in the LCT domain.



---

*Source: 536464-2012-04-04.xml*