arXiv:1706.08740v1 [math.CA] 27 Jun 2017 n a lowrite also can One ucin ftecniuu ore rnfr operator o transform basis Fourier continuous orthonormal complete the a of form functions functions Hermite that known enda follows as defined [ Grunbaum Keywords h edrcnfidsm xmlso uhcntutosi h pap the in constructions eigenbas such eige the of if this examples beneficial for some be ask find also can might would reader we it The so properties”, have numerically; would desirable computable and “other easily functions of Hermite to terms analogous In be would which DFT, [email protected] eerhltrtr:Hwcnw osrc notoomlbasis orthonormal an construct (DFT) we transform can How Fourier results discrete literature: the research main about question the following The and Introduction 1 Classification Subject Mathematics 2010 principle rclw oad mi:[email protected] Email: Wroc Poland. law, ∗ † e ssaetemi rpriso emt ucin htwl eus be will that functions Hermite of properties main the state us Let et fMteaisadSaitc,Yr nvriy 70KeeS Keele 4700 University, York Statistics, and Mathematics of Dept. aut fPr n ple ahmtc,Wo a nvriyo Tec of University Wroc law Mathematics, Applied and Pure of Faculty htte lostsyatretr eurnerlto and relation recurrence direct three-term are a DFT as satisfy the also of approp eigenvectors they an of optimal that sense these that the that in show show optimal we is we that paper DFT this the of In eigenbasis (DFT). transform Fourier discrete iia emt-yeegnai ftedsrt Fourier discrete the of eigenbasis Hermite-type Minimal N hr xs aywy obida rhnra ai of basis orthonormal an build to ways many exist There ievcos iceeFuirtasom rhgnlbasis, orthogonal transform, Fourier discrete eigenvectors, : nrae oinfinity. to increases 3 ,Mha[ Mehta ], ψ n ( 7 x n e n hn [ Chang and Pei and ] ( = ) ψ lxyKuznetsov Alexey n ( √ x =( := ) π 2 n n ( F !) − − g 1) 1 )( / rmr 23,Scnay65T50 Secondary 42A38, Primary : n 2 x ( e √ = ) ue2,2017 28, June − transform x π 2 8 2 / Abstract √ ]. 2 n H n ∗ 1 2 aes Kwa´snicki Mateusz , !) n π 1 ( − x Z 1 ,where ), / R 2 e e F − x endas defined , 2 i xy / 2 g d htte ovret emt functions Hermite to converge they that R d ( x n y N re,Trno N 3 P,Cnd.Email: Canada. 1P3, M3J ON, Toronto, treet, of n H it netit rnil.Moreover, principle. uncertainty riate )d ossigo h ievcoso the of eigenvectors the of consisting , e nlg,Wbz˙eWsi´sig 7 50-370 Wyspia´nskiego 27, Wybrze˙ze hnology, n nlge fteHriefunctions, Hermite the of analogues hr sol n uhorthonormal such one only is there R − y, x r emt oyoil.I swell- is It polynomials. Hermite are N 2 f a eevdalto teto nthe in attention of lot a received has L ossigo h ievcoso the of eigenvectors the of consisting r yDcisnadSegiz[ Steiglitz and Dickinson by ers bsst eepii r tleast, at or, explicit be to nbasis n 2 ( R emt ucin,uncertainty functions, Hermite ≥ dltr emt ucin are functions Hermite later. ed eohrdsrbeproperties? desirable other me n htte r h eigen- the are they that and ) swsi oesneunique. sense some in was is 0 , † (1) (2) 1 ], so that ( ψ )(x)=( i)nψ (x). F n − n Many constructions of the eigenbasis of the DFT start either with the fact that the Hermite functions are the eigenfunctions of the operator := d2/dx2 x2, or that they satisfy the recursion relations L −
ψn′ (x)+ xψn(x)= √2nψn 1(x), ψn′ (x) xψn(x)= 2(n + 1)ψn+1(x). − − − The goal, typically, is to find a discrete counterpart of the operator por or the operators d/dx + x and d/dx x and to use them to construct discrete analogues of HermiteL functions. A major problem with − this approach is that there are many ways to approximate a differential operator by a matrix, and this leads to many different constructions and results in lack of uniqueness of the “canonical eigenbasis” of the DFT. However, we do not need to be restricted to using differential operators if we want to characterize Hermite functions. It turns out that a much more fruitful approach, in terms of finding discrete analogues of Hermite functions, is to use Hardy’s uncertainty principle. This result was derived back in 1933 and m x2/2 it states that if both functions f and f are O( x e− ) for large x and some m, then f and f are finite linear combinations of HermiteF functions| (see| [4, Theorem 1]). The following theorem followsF easily from Hardy’s result. Theorem 1.1 (Uncertainty principle characterization of Hermite functions). Assume that the functions f : R R satisfy the following three conditions: n 7→ (i) f are the eigenfunctions of the Fourier transform: ( f )(x)=( i)nf (x) for n 0 and x R; n F n − n ≥ ∈ R (ii) fn n 0 form an orthonormal set in L2( ), that is R fn(x)fm(x)dx = δn,m; { } ≥ n 1 x2/2 (iii) for every n 0 we have x− − e f (x) 0 as x R . ≥ n → | |→∞ Then for all n 0 we have f ψ or f ψ . ≥ n ≡ n n ≡ − n The main goal of our paper is to find a set of vectors Tn 0 n N Theorem 1.3. For every N 2 there exist unique (up to change of sign) vectors Tn 0 n Definition 1.4. The basis Tn 0 n The vectors Tn 0 n 3 dim(E0) dim(E1) dim(E2) dim(E3) N =4L L +1 L L L 1 − N =4L +1 L +1 L L L N =4L +2 L +1 L L +1 L N =4L +3 L +1 L +1 L +1 L Table 1: Dimensions of the eigenspaces of the DFT. where Km := (N +2+ m)/4 . The result (4) is usually presented in the form of a table, by considering different residue⌊ classes of N modulo⌋ 4, see Table 1. The dimensions of the eigenspaces of the DFT were first found by Schur in 1921, however they can also be easily obtained from a much earlier result of Gauss on the law of quadratic reciprocity. By using this result and the fact that the vectors Tn 0 n N Proposition 1.5. Let Tn 0 n We have argued above that the minimal Hermite-type eigenvectors Tn 0 n L a(k)= a(k +1)+ a(k 1) + 2cos(ωk)a(k). (5) − Note that here we interpret vectors a and L a as N-periodic functions on Z, see the discussion on page 2. It is clear that L is a self-adjoint operator and it is also known that it commutes with the DFT (see [1]). N Theorem 1.6. Let Tn 0 n