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 China61402044Beijing City Board of Education Science and Technology PlanKM2017112320091. 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 [1]. The traditional Poisson sum formula can be represented as follows [2]: (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 [5], 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 [9–11], 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 [12–14]. 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 [15]. 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. Preliminaries2.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,udtb≠0dejcd/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 b≠0 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 [15]: (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/2b∗ht.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 SAFTTheorem 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/2bt2∑n=-∞+∞e-jdbn2/2τ2+bnw0/τFAbnτ+u0ejnt/τ,(13)∑k=-∞+∞fkτeja/2bk2τ2=j2πbτe-ju0w0∑n=-∞+∞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-ju0w0∑n=-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 n≠0. 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/2bt2∑n=-∞+∞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/2bt2∑n=-∞+∞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 OHTTheorem 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/2by2∑n=-∞+∞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).
LiB.-Z.TaoR.XuT.-Z.WangY.The poisson sum formulae associated with the fractional fourier transform200989585185610.1016/j.sigpro.2008.10.0302-s2.0-58849116271ZhangJ.-F.HouS.-P.The generalization of the poisson sum formula associated with the linear canonical transform201220121910203910.1155/2012/1020392-s2.0-84872769004AbeS.SheridanJ. T.Generalization of the fractional Fourier transformation to an arbitrary linear lossless transformation: an operator approach1994271241794187MR128261710.1088/0305-4470/27/12/0232-s2.0-21344493264AbeS.SheridanJ. T.Optical operations on wave functions as the abelian subgroups of the special affine fourier transformation199419221801180310.1364/OL.19.0018012-s2.0-0028546432PeiS.-C.DingJ.-J.Eigenfunctions of the offset fourier, fractional fourier, and linear canonical transforms200320352253210.1364/JOSAA.20.000522MR1966689PeiS.-C.DingJ.-J.Eigenfunctions of fourier and fractional fourier transforms with complex offsets and parameters20075471599161110.1109/TCSI.2007.900182MR2369822CaiL. Z.Special affine Fourier transformation in frequency-domain20001854-62712762-s2.0-003431675110.1016/S0030-4018(00)01005-1XiangQ.QinK.Convolution, correlation, and sampling theorems for the offset linear canonical transform2014834334422-s2.0-8489427241110.1007/s11760-012-0342-0HealyJ. J.OzaktasH. M.Sampling and discrete linear canonical transforms2016198New York, NY, USASpringer24125610.1007/978-1-4939-3028-9_8MR3467308XuT. Z.LiB. Z.2013Beijing, BJ, ChinaScience PressMoshinskyM.QuesneC.Linear canonical transformations and their unitary representations1971121772178010.1063/1.1665805MR02863852-s2.0-36849112629TaoR.UrynbassarovaD.Li BZ.The wigner-ville distribution in the linear canonical transform domain2016464559563LiB.-Z.TaoR.WangY.New sampling formulae related to linear canonical transform20078759839902-s2.0-3384657149810.1016/j.sigpro.2006.09.008Zbl1186.94201FengQ.LiB.Convolution theorem for fractional cosine-sine transform and its application201740103651366510.1002/mma.4251XiangQ.QinK.-Y.HuangQ.-Z.Multichannel sampling of signals band-limited in offset linear canonical transform domains20133252385240610.1007/s00034-013-9575-6MR31007212-s2.0-84883322977LohmannA. W.MendlovicD.ZalevskyZ.Fractional Hilbert transform19962142812832-s2.0-000038006110.1364/OL.21.000281DavisJ. A.McNamaraD. E.CottrellD. M.Analysis of the fractional Hilbert transform19983729691169132-s2.0-000046110210.1364/AO.37.006911GoelN.SinghK.Convolution and correlation theorems for the offset fractional Fourier transform and its application20167021381502-s2.0-8495363590810.1016/j.aeue.2015.10.009BhandariA.ZayedA. I.Shift-invariant and sampling spaces associated with the fractional Fourier transform domain20126041627163710.1109/TSP.2011.2177260MR2934185SheikhA.GudadheA.Poisson summation formulae associated with the generalized fractional Hilbert transform2016GudadheS. A.DeshmukhR. P.Poisson summation formula associated with the fractional Laplace transform2012SternA.Sampling of compact signals in offset linear canonical transform domains2007143593672-s2.0-3814900689010.1007/s11760-007-0029-0Zbl1122.94023ZhiX.WeiD.ZhangW.A generalized convolution theorem for the special affine Fourier transform and its application to filtering201612752613261610.1016/j.ijleo.2015.11.2112-s2.0-84959016712XuS.ChaiY.HuY.JiangC.LiY.Reconstruction of digital spectrum from periodic nonuniformly sampled signals in offset linear canonical transform domain201534859652-s2.0-8492530305310.1016/j.optcom.2015.03.027