MPE Mathematical Problems in Engineering 1563-5147 1024-123X Hindawi 10.1155/2017/1354129 1354129 Research Article Poisson Summation Formulae Associated with the Special Affine Fourier Transform and Offset Hilbert Transform http://orcid.org/0000-0002-8277-8160 Zhuo Zhi-Hai 1 Lo Schiavo Alessandro Beijing Information Science & Technology University Beijing China bistu.edu.cn 2017 1582017 2017 07 05 2017 16 07 2017 1582017 2017 Copyright © 2017 Zhi-Hai Zhuo. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

This paper investigates the generalized pattern of Poisson summation formulae from the special affine Fourier transform (SAFT) and offset Hilbert transform (OHT) points of view. Several novel summation formulae are derived accordingly. Firstly, the relationship between SAFT (or OHT) and Fourier transform (FT) is obtained. Then, the generalized Poisson sum formulae are obtained based on above relationships. The novel results can be regarded as the generalizations of the classical results in several transform domains such as FT, fractional Fourier transform, and the linear canonical transform.

National Natural Science Foundation of China 61402044 Beijing City Board of Education Science and Technology Plan KM201711232009
1. Introduction

The classical Poisson summation formula, which demonstrates that the sum of infinite samples in the time domain of a signal f(t) is equivalent to the sum of infinite samples of F(·) in the Fourier domain, is of importance in theories and applications of signal processing . The traditional Poisson sum formula can be represented as follows : (1)k=-+ft+kτ=1τn=-+Fnτejnt/τor (2)k=-+fkτ=1τn=-+Fnτ,t=0,where F(·) denotes the Fourier transform (FT) of a signal f(t) and τ stands for the sampling interval. Not only does Poisson summation formula play a key role in various branches of the mathematics, but also it finds numerous applications in lots of fields, for example, mechanics, signal processing community, and many other scientific fields. The Poisson summation formula is related to the Fourier transform, and, with the development of modern signal processing technologies, there are many other kinds of transforms that have been proposed, it is therefore worthwhile and interesting to investigate the Poisson sum formula in deep associated with these kinds of new integral transforms.

The special affine Fourier transform (SAFT) [3, 4], also known as the offset linear canonical transform [5, 6] or the inhomogeneous canonical transform , is a six-parameter (a,b,c,d,u0,w0) class of linear integral transform. Many well-known transforms in signal processing and optic systems are its special cases such as Fourier transform (FT), fractional Fourier transform (FRFT), the linear canonical transform (LCT), time shifting and scaling, frequency modulation, pulse chirping, and others [7, 8]. SAFT can be interpreted as a time shifting and frequency modulated version of LCT , that is much more flexible because of its extra parameters (u0,w0). Recently, it has been widely noticed in many practical applications along with the rapid development of LCT . Thus, developing relevant theorems for SAFT is of importance and necessary in optical systems and many signal processing applications as well.

In addition, the generalized Hilbert transform closely related to SAFT, called offset Hilbert transform (OHT), is another powerful tool in the fields of optics and signal processing community . It has been presented recently and widely used for image processing, especially for edge detection and enhancement, because it can emphasize the derivatives of the image [16, 17]. In recent decades, many essential theories and useful applications of SAFT and OHT have been derived from in-depth researching on it [8, 15, 18, 19]. To the best of our knowledge, Poisson sum formula has been generalized with many transforms such as FRFT, LCT, fractional Laplace transform and fractional Hilbert transform [1, 2, 20, 21]. However, none of the research papers throw light on the study of the traditional Poisson sum formula associated with the SAFT and OHT yet. Based on the existing results, the motivation of this paper is to generalize the above-mentioned Poisson sum formula into SAFT and OHT domains.

The rest of this paper is organized as follows. Section 2 gives some fundamental knowledge of SAFT and OHT. In Section 3, we give the relationships between SAFT/OHT and FT in detail. Some novel Poisson summation formulae associated with SAFT are presented in Section 4. Section 5 concludes the paper.

2. Preliminaries 2.1. The Special Affine Fourier Transform

The special affine Fourier transform (SAFT) with real parameters A=(a,b,c,d,u0,w0) of a signal f(t) is defined by the following [5, 22]: (3)FAu=OLAftu=-+fthAt,udtb0dejcd/2u-u02+jw0ufdu-u0b=0,where (4)hAt,u=KAej/2bat2+2tu0-u-2udu0-bw0+du2,KA=1j2πbejd/2bu02and ad-bc=1. Note that, for b=0, the SAFT of a signal is essentially a chirp multiplication and it is of no particular interest for our objective in this work. Hence, without loss of generality, we set b0 in the following section unless stated otherwise. The inverse of an SAFT with parameter A=(a,b,c,d,u0,w0) is given by an SAFT with parameter A-1=(d,-b,-c,a,bw0-du0,cu0-aw0), which is (5)ft=OLA-1FAut=C-+FAuhA-1u,tdu,where C=ej1/2(cdu02-2adu0w0+abw02). This can be verified by the definition of SAFT. Most of important transforms can be its special cases when parameter A is replaced with specific parameters. For example, when A=(0,1,-1,0,0,0), SAFT coincides with FT; when A=(cosα,sinα,-cosα,sinα,0,0), SAFT is FRFT; when A=(a,b,c,d,0,0), SAFT equals LCT. Furthermore, many important theories on SAFT have been investigated [8, 15, 23, 24].

2.2. Offset Hilbert Transform

The offset Hilbert transform (OHT) of a signal f(t) is defined as follows : (6)fOAt=HOAft=p.v.e-jat2+2u0t/2bπ-+fxt-xejax2+2u0x/2bdx.It should be noted that the above definition uses the Cauchy principal value of the integral (denoted here by p.v.). To obtain the relationship between the stand and HT and OHT, we can rewrite (6) as (7)fOAt=HOAft=e-jat2+2u0t/2bftejat2+2u0t/2bht.Notice that computing the OHT of a signal f is equivalent to multiplying it by a chirp, ejat2+2u0t/2b, then passing the product through a standard Hilbert filter and finally multiplying the output by the chirp, e-jat2+2u0t/2b. This relationship between OHT and the classical HT can be shown in Figure 1.

Offset Hilbert filter.

3. The Relationships between SAFT/OHT and Fourier Transform

In order to derive novel Poisson summation formulae based on SAFT and OHT, some relationships between SAFT/OHT and FT are obtained in this section firstly.

Lemma 1.

Suppose the SAFT of a signal f(t) with parameters A=(a,b,c,d,u0,w0) is FA(·), and set g(t)=f(t)eja/2bt2, and then the following relations hold: (8)FAu=1j2πbej1/2bdu02+u2-2udu0-bw0Gu-u0b,Gv=j2πbe-j1/2dbv2+2bvw0+2u0w0FAbv+u0,where G(·) is the FT of signal g(t).

Proof.

It is easy to verify Lemma 1 by the definitions of SAFT and FT.

Lemma 2.

Suppose the SAFT of a signal f(t) with parameters A=(a,b,c,d,u0,w0) is FA(·), and set q(x)=f(t-x)/xeja/2bx2, and then the following relations hold: (9)fOAt=1πQat+u0b,Qv=πfOAbv-u0a,where Q(·) is the FT of signal q(x).

Proof.

According to the definition of OHT, (10)fOAt=e-jat2+2u0t/2bπ-+fxejax2+2u0x/2bt-xdx=e-jat2+2u0t/2bπ-+ft-xejat-x2+2u0t-x/2bxdx=1π-+ft-xxeja/2bx2e-jat+u0/bxdx=1π-+qxe-jat+u0/bxdx=1πQat+u0b.By replacing v=at+u0/b, (10) can be rewritten as (11)fOAbv-u0a=1πQv.This completes the proof of Lemma 2.

4. Main Results

Based on the relationships in Lemmas 1 and 2, the generalized Poisson summation formulae associated with SAFT and OHT are derived in following subsections, respectively.

4.1. The Poisson Sum Formula Based on SAFT Theorem 3.

The Poisson summation formulae of a signal f(t) in the SAFT domain with parameter A are (12)k=-+ft+kτeja/2b2ktτ+k2τ2=j2πbτe-ju0w0+a/2bt2n=-+e-jdbn2/2τ2+bnw0/τFAbnτ+u0ejnt/τ,(13)k=-+fkτeja/2bk2τ2=j2πbτe-ju0w0n=-+e-jdbn2/2τ2+bnw0/τFAbnτ+u0,t=0.

Proof.

If we set g(t)=f(t)eja/2bt2, from the traditional Poisson sum formula for g(t) in the Fourier domain, that is, (1), we obtain (14)k=-+gt+kτ=1τn=-+Gnτejnt/τ.By directly using Lemma 1, we derive that (15)k=-+ft+kτeja/2bt+kτ2=1τn=-+j2πbe-j1/2bdb2n/τ2+2b2n/τw0+2bu0w0FAbnτ+u0ejnt/τ.Theorem 3 is proved by simple calculation on (15).

Equations (13) and (15) can be regarded as the generalization of classical Poisson sum formula based on SAFT. It should be noticed that when the parameters of the SAFT are chosen to be the special cases of the SAFT, the derived results reduce to the classical results of Fourier transform domain, fractional Fourier transform domain, and linear canonical transform domains. It clearly demonstrates that the infinite sum of periodic phase-shifted replica of a signal f(t) in the time domain is equivalent to the infinite sum of periodic phase-shifted replica FA(·) in the SAFT domain.

In addition, it is of importance to investigate the Poisson sum formula of signals with compact support in SAFT domain. A signal f(t) is said to have compact support in SAFT domain if its SAFT FA(u)>ΩA, where ΩA>0 is some real number. Without loss of generality, let a>0,b>0,u0>0 in the following analysis.

Corollary 4.

Suppose a signal f(t) is band-limited in SAFT domain with a compact support ΩA; then the Poisson sum formula derived in Theorem 3 can be rewritten as the following forms according to the replica period τ:

When b/τ>ΩA+u0 and ΩA>u0, (16)k=-+ft+kτeja/2b2ktτ+k2τ2=j2πbτe-ja/2bt2e-ju0w0FAu0.

When ΩA+u0/2<b/τ<ΩA-u0 and ΩA>3u0,(17)k=-+ft+kτeja/2b2ktτ+k2τ2=j2πbτe-ja/2bt2e-ju0w0×FAu0+e-j1/2db/τ2+2bw0/τFAbτ+u0ejt/τ+e-j1/2db/τ2-2bw0/τFA-bτ+u0e-jt/τ.

When ΩA+u0/m+1<b/τ<ΩA-u0/m and ΩA>(2m+1)u0, (18)k=-+ft+kτeja/2b2ktτ+k2τ2=j2πbτe-ja/2bt2e-ju0w0n=-mme-jdbn2/2τ2+bnw0/τFAbnτ+u0ejnt/τ.

Proof.

(a) Since f(t) is a band-limited signal in SAFT domain with a compact support ΩA, it is easy to derive the right hand FAbn/τ+u0 of (12) that is equal to zeros when n0. That is, from b/τ>ΩA+u0, we derive that b/τ+u0>b/τ-u0>ΩA and -2b/τ+u0<-b/τ+u0<-ΩA. Thus, it is easy to derive that (19)j2πbτe-ju0w0+a/2bt2n=-+e-jdbn2/2τ2+bnw0/τFAbnτ+u0ejnt/τ=j2πbτe-ja/2bt2e-ju0w0FAu0.Substituting (19) into (12) yields the final results.

(b) It is easy to prove that only when n=-1,0,+1, FAbn/τ+u0 is nonzero. The right hand of (12) is (20)j2πbτe-ju0w0+a/2bt2n=-+e-jdbn2/2τ2+bnw0/τFAbnτ+u0ejnt/τ=j2πbτe-ja/2bt2e-ju0w0FAu0+e-j1/2db/τ2+2bw0/τFAbτ+u0ejt/τ+e-j1/2db/τ2-2bw0/τFA-bτ+u0e-jt/τ.Substituting (20) into (12) yields the final results.

(c) The proof of this situation is similar to the proof of (a) and (b), and we omit it here.

4.2. The Poisson Sum Formula Based on OHT Theorem 5.

The Poisson sum formula of a signal f(t) in the OHT domain with parameter A is as follows: (21)k=-+ft-y-kτy+kτeja/2b2kyτ+k2τ2=πτe-ja/2by2n=-+fOAbn-u0τaτejny/τ.

Proof.

If we set q(x)=f(t-x)/xeja/2bx2, it is easy to verify Theorem 5 via (1) and Lemma 2: (22)k=-+qy+kτ=1τn=-+Qnτejny/τ.That is (23)k=-+ft-y-kτy+kτeja/2by+kτ2=πτn=-+fOAbn-u0τaτejny/τBy simple calculation, Theorem 5 is completed.

Equation (21) can be seen as the Poisson sum formula associated with offset Hilbert transform. Furthermore, it is worthwhile and interesting to study the signals with compact support in offset Hilbert transform domain. Let fOA(y) be the OHT of a signal f(t). Then f(t) is said to have compact support in OHT domain, if fOA(y)=0 for y>ΩA, where ΩA>0 is some real number.

Corollary 6.

Suppose signal f(t) is band-limited in OHT domain with a compact support ΩA; then the Poisson sum formula derived in Theorem 5 can be rewritten as the following forms according to the replica period τ:

When b/aτ>ΩA+u0/a and ΩA>u0/a, (24)k=-+ft-y-kτy+kτeja/2by+kτ2=πτfOA-u0a.

When ΩA+u0/a/2<b/τ<ΩA-u0/a and ΩA>3u0/a, (25)k=-+ft-y-kτy+kτeja/2by+kτ2=πτfOA-b-u0τaτe-jy/τ+fOA-u0a+fOAb-u0τaτejy/τ.

When ΩA+u0/a/m+1<b/τ<ΩA-u0/a/m and ΩA>(2m+1)u0/a, (26)k=-+ft-y-kτy+kτeja/2by+kτ2=πτn=-mmfOAbn-u0τaτejny/τ.

Proof.

It is easy to verify this corollary using Theorem 5 and the similar method in Corollary 4.

5. Conclusion

In this paper, the traditional Poisson summation formula has been generalized into SAFT and OHT domain. Theorems 3 and 5 are the generalizations of Poisson summation formulae based on SAFT and OHT, respectively. In addition, signals with compact support are mostly used in signal processing and considered in this paper as well. Some novel results associated with Poisson summation formula have been derived in the form of Corollaries 4 and 6.

Conflicts of Interest

The author declares that there are no conflicts of interest regarding the publication of this paper.

Acknowledgments

This work was supported in part by the National Natural Science Foundation of China (no. 61402044) and Beijing City Board of Education Science and Technology Plan (no. KM201711232009).

Li B.-Z. Tao R. Xu T.-Z. Wang Y. The poisson sum formulae associated with the fractional fourier transform Signal Processing 2009 89 5 851 856 10.1016/j.sigpro.2008.10.030 2-s2.0-58849116271 Zhang J.-F. Hou S.-P. The generalization of the poisson sum formula associated with the linear canonical transform Journal of Applied Mathematics 2012 2012 1 9 102039 10.1155/2012/102039 2-s2.0-84872769004 Abe S. Sheridan J. T. Generalization of the fractional Fourier transformation to an arbitrary linear lossless transformation: an operator approach Journal of Physics. A. Mathematical and General 1994 27 12 4179 4187 MR1282617 10.1088/0305-4470/27/12/023 2-s2.0-21344493264 Abe S. Sheridan J. T. Optical operations on wave functions as the abelian subgroups of the special affine fourier transformation Optics Letters 1994 19 22 1801 1803 10.1364/OL.19.001801 2-s2.0-0028546432 Pei S.-C. Ding J.-J. Eigenfunctions of the offset fourier, fractional fourier, and linear canonical transforms Journal of the Optical Society of America A. Optics, Image Science, and Vision 2003 20 3 522 532 10.1364/JOSAA.20.000522 MR1966689 Pei S.-C. Ding J.-J. Eigenfunctions of fourier and fractional fourier transforms with complex offsets and parameters IEEE Transactions on Circuits and Systems. I. Regular Papers 2007 54 7 1599 1611 10.1109/TCSI.2007.900182 MR2369822 Cai L. Z. Special affine Fourier transformation in frequency-domain Optics Communications 2000 185 4-6 271 276 2-s2.0-0034316751 10.1016/S0030-4018(00)01005-1 Xiang Q. Qin K. Convolution, correlation, and sampling theorems for the offset linear canonical transform Signal, Image and Video Processing 2014 8 3 433 442 2-s2.0-84894272411 10.1007/s11760-012-0342-0 Healy J. J. Ozaktas H. M. Sampling and discrete linear canonical transforms Linear canonical transforms 2016 198 New York, NY, USA Springer 241 256 10.1007/978-1-4939-3028-9_8 MR3467308 Xu T. Z. Li B. Z. Linear Canonical Transform and Its Applications 2013 Beijing, BJ, China Science Press Moshinsky M. Quesne C. Linear canonical transformations and their unitary representations Journal of Mathematical Physics 1971 12 1772 1780 10.1063/1.1665805 MR0286385 2-s2.0-36849112629 Tao R. Urynbassarova D. Li B Z. The wigner-ville distribution in the linear canonical transform domain IAENG International Journal of Applied Mathematics 2016 46 4 559 563 Li B.-Z. Tao R. Wang Y. New sampling formulae related to linear canonical transform Signal Processing 2007 87 5 983 990 2-s2.0-33846571498 10.1016/j.sigpro.2006.09.008 Zbl1186.94201 Feng Q. Li B. Convolution theorem for fractional cosine-sine transform and its application Mathematical Methods in the Applied Sciences 2017 40 10 3651 3665 10.1002/mma.4251 Xiang Q. Qin K.-Y. Huang Q.-Z. Multichannel sampling of signals band-limited in offset linear canonical transform domains Circuits, Systems, and Signal Processing 2013 32 5 2385 2406 10.1007/s00034-013-9575-6 MR3100721 2-s2.0-84883322977 Lohmann A. W. Mendlovic D. Zalevsky Z. Fractional Hilbert transform Optics Letters 1996 21 4 281 283 2-s2.0-0000380061 10.1364/OL.21.000281 Davis J. A. McNamara D. E. Cottrell D. M. Analysis of the fractional Hilbert transform Applied Optics 1998 37 29 6911 6913 2-s2.0-0000461102 10.1364/AO.37.006911 Goel N. Singh K. Convolution and correlation theorems for the offset fractional Fourier transform and its application AEU - International Journal of Electronics and Communications 2016 70 2 138 150 2-s2.0-84953635908 10.1016/j.aeue.2015.10.009 Bhandari A. Zayed A. I. Shift-invariant and sampling spaces associated with the fractional Fourier transform domain IEEE Transactions on Signal Processing 2012 60 4 1627 1637 10.1109/TSP.2011.2177260 MR2934185 Sheikh A. Gudadhe A. Poisson summation formulae associated with the generalized fractional Hilbert transform Asian Journal of Mathematics and Computer Research 2016 Gudadhe S. A. Deshmukh R. P. Poisson summation formula associated with the fractional Laplace transform Journal of Science and Arts 2012 Stern A. Sampling of compact signals in offset linear canonical transform domains Signal, Image and Video Processing 2007 1 4 359 367 2-s2.0-38149006890 10.1007/s11760-007-0029-0 Zbl1122.94023 Zhi X. Wei D. Zhang W. A generalized convolution theorem for the special affine Fourier transform and its application to filtering Optik-International Journal for Light and Electron Optics 2016 127 5 2613 2616 10.1016/j.ijleo.2015.11.211 2-s2.0-84959016712 Xu S. Chai Y. Hu Y. Jiang C. Li Y. Reconstruction of digital spectrum from periodic nonuniformly sampled signals in offset linear canonical transform domain Optics Communications 2015 348 59 65 2-s2.0-84925303053 10.1016/j.optcom.2015.03.027