On the Duality of Discrete and Periodic Functions

: Although versions of Poisson’s Summation Formula (PSF) have already been studied extensively, there seems to be no theorem that relates discretization to periodization and periodization to discretization in a simple manner. In this study, we show that two complementary formulas , both closely related to the classical Poisson Summation Formula, are needed to form a reciprocal Discretization-Periodization Theorem on generalized functions. We deﬁne discretization and periodization on generalized functions and show that the Fourier transform of periodic functions are discrete functions and, vice versa, the Fourier transform of discrete functions are periodic functions.


Introduction
Versions of Poisson's summation formula have already been the subject of many publications and it is clear that it plays a very central role in mathematics, physics and electrical engineering.Many publications describe the interaction of discretization (or sampling) and periodization.In some publications we find a clear statement that the Poisson's Summation Formula (PSF) is linked to both discretization and periodization [1][2][3], and in rare cases we find that the PSF actually is a discretization-periodization theorem [4].Despite the circumstance that PSF versions are commonly understood as discretization -periodization relationships, it appears to be difficult to find a simple statement on the special interrelation between discretization and periodization.In this paper, we show that two complementary formulas, very closely related to Poisson's classical summation formula, are needed to form a complete, i.e., reciprocal, discretization-periodization theorem in the tempered distribution sense.
Most publications about discretization and periodization relationships are so far application driven.A comprehensive derivation of interactions between sampling (discretization) and periodization has already been performed for Gabor systems.Orr and Janssen initially covered the transition from continuous settings to discrete settings [5,6] and later Kaiblinger [7] covered the transition from discrete settings to continuous settings.Finally, Søndergaard summarized and complemented these results in [8] and, amongst other things, wrote a Poisson Summation Formula without sums but with sampling and periodization operators instead ( [4], Theorem A.1), very similar to the derivation in this paper but in S 0 (R), Feichtinger's algebra [9].
Another comprehensive treatment of interactions between discretization and periodization is given in Ozaktas [10], whereby discreteness and periodicity are characterized as Fourier duals such as translation and modulation.Amongst others, the notion of discrete functions is clarified as being Fourier duals of periodic functions.Other publications establish a connection between sampling and periodization, on the one hand, and Shannon's Sampling Theorem on the other [11] without compelling necessity to mention Poisson's summation formula in this context.There are even excellent signal processing textbooks that can manage without Poisson's summation formula but not without describing its impacts, e.g., [12].Poisson's summation formula also plays an important role in wavelet theory [13][14][15], although only in proofs, because it provides transitions from continuous to discrete settings.
Lastly, Benedetto and Zimmermann investigated the validity of Poisson's summation formula in several classical spaces in analysis including the space of tempered distributions [3,16].They define sampling and periodization operators and, amongst others, characterize uniform sampling as the Fourier transform of periodization.
The Poisson Summation Formula is also honored as being the link between the classical Fourier transform and Fourier series in Girgensohn [17], Hunter and Nachtergaele [18] and Strichartz [19].Indeed, we will see that this link is discretization, formally correct only in the sense of generalized functions, called distributions.Hence, in order to formally treat discrete functions correctly, we will need the theory of distributions, basically developed by Laurent Schwartz in [20,21], with contributions by Temple [22] and Lighthill [2], Gelfand and Schilow [23,24] and Zemanian [25].For our purposes, it is especially convenient to show everything within the space of tempered distributions S ′ , because the Fourier transform of a tempered distribution is again a tempered distribution [21,26].Note that the Fourier transform definition can even be more generalized to spaces larger than S ′ , see [27] for example, but extending these results accordingly is beyond the scope of this paper.
Conversely, results presented in this paper may also hold in less general settings.For engineering purposes, it often suffices to apply the definitions of both sampling and periodization and use their symbols, but neglect their tempered distribution nature, in order to simplify the already established formulas found in textbooks.A first impression of this technique is given in Section 5 that uses almost no generalized functions theory.Section 2 will first familiarize the reader with the notations used throughout this article.Section 3 explains the proof method; readers familiar with distribution theory may skip this section.In Section 4 we define discretization and periodization as operations on suitable subspaces of tempered distributions and in Section 5 we hope to motivate the reader to apply them to already known results in textbooks.Section 6 provides calculation rules for both discretization and periodization, such that the proof of our theorem in Section 7 reduces to two lines.Finally, we summarize our results in Section 8 and provide an outlook in Section 9.

Notation
We mainly follow established notations as in Schwartz [21], including his Fourier transform definition Equation (1), Benedetto [3], Strichartz [26] and Walter [27], and we substantially follow Woodward [1] and Bracewell [28].One of the primary goals of this paper is to extend Bracewell's idea of establishing symbolic calculation rules in signal processing for educational purposes.His textbook Fourier Transform and its Applications is a highly recognized standard work in electrical engineering literature.By replacing Bracewell's Shah symbol III by two Shah-like symbols ⊥ ⊥ ⊥ and △ △ △ we achieve three things: (1) we are able to express the interrelation between discretization and periodization, (2) we are able to express the duality of both symbols and (3) we overcome difficulties that arise when conventional dissimilar symbols are used.
The application of a distribution f ∈ S ′ to some test function φ ∈ S is denoted as ⟨f, φ⟩ and, in particular, if f is a locally integrable function such that the following integral exists, we define it by For singular distributions, e.g., δ ∈ S ′ , it is required to write down their definition explicitly, i.e., ⟨δ, φ⟩ := φ(0).The Fourier Transform F for integrable functions φ ∈ L 1 (R n ) is defined by where i = √ −1 and t • σ is the usual scalar product in R n .For generalized functions f ∈ S ′ we define its Fourier transform by ⟨Ff, φ⟩ := ⟨f, Fφ⟩ i.e., F is "rolled off" to test functions φ, a standard technique for definitions in distribution theory.Defining the Fourier transform in that way implies that Fδ = 1 and F1 = δ, see for example [21].Also note that F −1 Ff = f and FF −1 f = f for any f ∈ S ′ , see [21,26], and FFφ = φ where φ(t) := φ(−t) for test functions φ ∈ S. Furthermore, we define Fourier Series F S for integrable periodic functions The translation operator τ for some function f : R n → l C and some constant a ∈ R n is defined by (τ a f )(t) := f (t − a).For generalized functions it is defined by ⟨τ a f, φ⟩ := ⟨f, τ −a φ⟩.
Additionally, we abbreviate τ a δ := δ a for translated Dirac impulses δ.Often in the literature we find a Dirac impulse train, also Dirac comb as Instead of δ k , we often use δ kT where By kT we mean component-wise multiplication with k ∈ Z n .Note that III ∈ S ′ , see [2,26,29] for example.
As mentioned already, the Fourier transform of a tempered distribution is again a tempered distribution.So, for convenience, we stay in S ′ .Additionally, we sometimes need to restrict ourselves to several subspaces of S ′ and would like to remind the reader that the following continuous embeddings of distribution spaces have been investigated in [21].We will need the space of test functions E := C ∞ , its dual space E ′ , the space of infinitely differentiable functions with compact support D := C ∞ c and the space of bounded test functions B and its dual space D ′ L 1 , we need the Schwartz space S of rapidly decreasing, infinitely differentiable functions and its dual space, the space of tempered distributions S ′ and exploit that F is an automorphism of S, i.e., F(S) = S and that this implies that F is an automorphism of S ′ , i.e., F(S ′ ) = S ′ .
Additionally, we require the space of slowly increasing test functions O M and the space of rapidly decreasing distributions to be the image of O M in S ′ with respect to F, i.e., F(O M ) = O C ′ and Recall that O M and O C ′ are the spaces of multiplication operators in S ′ and convolution operators in S ′ , respectively [21].Finally, all those spaces X(R n ) we abbreviate to X, for functions φ : R n → l C depending on t ∈ R n we write φ instead of φ(t).For test functions in O M we use α, β, for test functions in S we use φ, ψ and for generalized functions in O C ′ or S ′ we use f, g.

Idea
In signal processing it is well-known that sampling a function means periodizing its spectrum.And, vice versa, periodizing a function means sampling its spectrum.Our goal is now, first, to formulate this in the simplest way and, secondly, to prove it the widest possible way.We are especially interested under which conditions these statements hold in the widest sense.
Denoting a generalized function as f and its spectrum as Ff , we additionally require a symbol for discretization, say ⊥ ⊥ ⊥, and one for periodization, say △ △ △.We will then write a sampled function as f s = ⊥ ⊥ ⊥f and a periodized function as f p = △ △ △f .This notation is deliberately chosen such that it reminds us of linear algebra notations like y = Ax.Having a symbol for discretization, we are now able to express the relationship between Fourier series F s and the Fourier transform as F s (f ) = F( ⊥ ⊥ ⊥f ), i.e., first applying discretization and then F to f corresponds to what we know as the Fourier series of f , shortly F s = F • ⊥ ⊥ ⊥ where • is concatenation.Now, coming back to the first two sentences in this section, our goal is to prove that But what do these properties look like, why do we need distribution theory and does this really have applications?
Properties.We will see that f must be an element in O M ⊂ S ′ to fulfill the first equation, and an element in F(O M ) ⊂ S ′ to fulfill the second.

Distribution Theory
Indeed, the only way to prove these equations is using the theory of generalized functions (distribution theory) because sampling some function φ(t), often explained via definition is mathematically wrong.Either this integral would be zero according to the definition of δ or it contradicts the Lebesgue integral.Distribution theory has been developed to overcome this problem.Instead of the integral expression, we now write to apply δ to a function φ.The novelty is that arbitrarily defined "functions" f can now be applied to other functions φ without contradicting the Lebesgue integral.Hence, ⟨•, •⟩ can be considered as a "generalized integral" but, equivalently, we must take care that ⟨f, φ⟩ < ∞ holds for all "functions" f operating on functions φ.This is secured by using spaces X of φ having certain properties and spaces X ′ of f with complementary properties.So, for example, if f is slowly growing then φ must be rapidly decreasing.Vice versa, if f is rapidly decreasing then φ is allowed to be slowly growing.Coming back to our equations, we must show that if f ∈ X ′ then the above equations, e.g., for every test function φ ∈ X.Also, note that spaces of test functions are continuously embedded in spaces of generalized functions X ⊂ X ′ , such that test functions themselves are generalized functions.Moreover, in distribution theory every test function φ is required being infinitely differentiable (φ ∈ C ∞ ) such that derivatives of generalized functions can always be found no matter whether their derivatives do exist in the conventional sense or not.Distribution theory is therefore also regarded as the completion of Fourier analysis [2] or of differential calculus [26].

Applications
As usually done in distribution theory, we will claim in the following, e.g., Equations ( 5), ( 7), ( 10), ( 17), ( 21) and (23), that all functions involved must be test functions, i.e., they are in C ∞ .This is required in order to be able to roll-off the multiplication of functions with generalized functions to the multiplication of functions with test functions in definition ⟨ αf , φ⟩ := ⟨ f , αφ⟩ such that αφ ∈ S is secured for all φ ∈ S. The property of being infinitely differentiable is a pure technical requirement in distribution theory that does not restrict the theorem we prove for two reasons: (1) In distribution theory, every locally integrable function (considered as generalized function) is already infinitely differentiable and (2) distribution theory is connected to its applications via several approximation theorems (see e.g., [27,30]).For example, every continuous function but also every generalized function can be approximated by test functions, i.e., functions in the classical sense being infinitely differentiable.We therefore always operate (as usual) on either spaces of generalized functions or spaces of test functions.

Generality
The theorem we prove corresponds to the Fourier transform in S ′ , the widest possible comprehension of the Fourier transform as automorphism (according to Laurent Schwartz' distribution theory [20,21]) and is, hence, most general in this sense.It comprises, moreover, several classical definitions of the Fourier transform.We note, for example, that there is no need to have a special Fourier transform F P designed for periodic functions.Also constants, in particular also being periodic functions, can be Fourier transformed without difficulties.Finally, all results presented in this paper can be used in pure symbolic calculations once they have been proven (in S ′ ).Our theorem shows that discretization and periodization are intertwined by the Fourier transform, symbolically where • is an operator argument placeholder.

Definitions
We often find the Shah symbol III in the literature (Bracewell [28], Feichtinger [31] and Lecture Notes, e.g., [32]) as an abbreviation for the Dirac impulse train Equation (4).It is convenient to only have one symbol because the Dirac impulse train can be used to both discretize and periodize functions: Multiplication with it yields discrete functions, convolution with it yields periodic functions [28].Furthermore, having only one symbol avoids introducing a theory of operations on functions which is especially of advantage for reasons of brevity in engineering books.However, any interaction between discretization and periodization cannot be expressed via one symbol.

Discretization
A major difficulty in S ′ is that the product of two tempered distributions is generally not defined.We will therefore use Laurent Schwartz' [21] space of multiplication operators O M in the following definition.
Definition 1 (Discretization).For slowly increasing functions α ∈ O M ⊂ S ′ and increments T ∈ R n + we claim that the tempered distribution defined by converges in S ′ (Figure 1).The operation ⊥ ⊥ ⊥ T is called uniform discretization (or sampling) with increments T .We call ⊥ ⊥ ⊥ T α the discrete (sampled) function of α discretized (sampled) with increments T .Discretization (sampling) is a linear continuous operator Discretization will be well-defined in this sense and it is clear that these series indeed converge in S ′ , see [26] for example.This definition corresponds in fact to the multiplication of function α with a Dirac comb III T ∈ S ′ .It is today already a common engineering practice to model the sampling process by the multiplication with the Dirac comb [33].However, an operator approach is much more advantageous as we will see later.It can already be found as comb T (α) in [1], later re-published in [34].

Periodization
Analogous to discretization where the space of multiplication operators O M in S ′ is required, we will now need the space of convolution operators O C ′ in S ′ , again, introduced by Laurent Schwartz in [21].
Definition 2 (Periodization).For rapidly decreasing tempered distributions f ∈ O C ′ ⊂ S ′ and increments T ∈ R n + we claim that the tempered distribution defined by converges in S ′ (Figure 2).The operation △ △ △ T is called periodization (replicating) or periodic continuation of f with periods T .Periodization (replicating) is a linear continuous operator △ △ △ T : Again, it is clear that periodization is well-defined in this sense and that these series do converge in S ′ .
In some textbooks, this operator can be found as rep T f [1,34,35], as f per [15], as • f T [16], as P f [26], as f [36] as T f [37] or as f T [38].However, △ △ △ is similar to ⊥ ⊥ ⊥ and can, hence, be used for symbolic calculations as we will see later.We will in particular see that, by using ⊥ ⊥ ⊥ and △ △ △, especially Poisson Summation Formulas reduce to something that is much easier to understand.

Motivation
The interaction of discretization and periodization is ubiquitous in many calculations that we find in today's textbooks.In order to reveal these relations, in this section we use Equations ( 5) and ( 6), for the sake of simplicity, neglect considerations in the tempered distribution sense.We give three examples: first, we show that Fourier Series result from Fourier transforming discretized functions.Secondly, we develop periodized functions into Fourier series seeing that it is related to discretization.And, thirdly, we apply the definitions of Section 4 to known versions of Poisson's summation formula found in the literature.

Fourier Series and Fourier Transform
Using our definition of discretization we verify that Fourier series actually correspond to a concatenation of discretization followed by Fourier transformation.Hence, having F and ⊥ ⊥ ⊥, the definition of F S as given in Equation ( 3) is obsolete; we may now write F ⊥ ⊥ ⊥ instead of introducing Equation (3).

Lemma 1 (Fourier Series vs Fourier Transform). Let α ∈ O M , then
Proof.We may apply ⊥ ⊥ ⊥ to α ∈ O M and, accordingly, ⊥ ⊥ ⊥α converges in S ′ such that F( ⊥ ⊥ ⊥α) ∈ S ′ .Due to the convergence of ⊥ ⊥ ⊥α in S ′ we may calculate term by term.Symbolically, we obtain with equalities in the tempered distributions sense.
A consequence is the following equality in the operations sense where • is the concatenation of operations in S ′ .We will later come back to this result in Equation (25).

Fourier Series of Periodized Functions
For further motivation we show that a Discretization -Periodization Theorem for Schwartz functions φ ∈ S can already be derived from the above definitions: We construct the Fourier series of a periodically continued function φ ∈ S and see that it is indeed related to discretization.Note that every periodic function f can be written as f = △ △ △ϕ for some function ϕ.Because φ ∈ S ⊂ O C ′ it follows that △ △ △φ ∈ S ′ is a tempered distribution.As every periodic function can be developed into a sum of e k (t) = e 2πi t•k , k ∈ Z n , we will now do this for △ △ △φ.Starting from we finally arrive at (see [27,39,40] for example) Inserting now the definition of discretization yields this important relationship, a formula being closely related to Poisson's summation formula, between periodization and discretization so far for Schwartz functions φ ∈ S. We note that Fφ ∈ S ⊂ O M such that ⊥ ⊥ ⊥(F φ) is meaningful as well because ⊥ ⊥ ⊥ requires functions in O M according to Definition 1. Indeed, Equation ( 10) can be more generalized and an inversion formula can be found.It will, moreover, be surprising to learn at this point that both operators in Equations ( 5) and ( 6) have already been defined in Woodward's 1953 textbook [1] leading to Equations ( 28) and ( 29), meant on functions in S. In Section 7, we will extend these results to generalized functions f ∈ O C ′ ⊃ S and S ⊂ O M ∋ Ff (see also Appendix B).

Versions of Poisson Summation Formulas
In general, each Poisson Summation Formula is only one of two possible statements: It either turns discretization into periodization resulting in Equations ( 23) and (28) or it turns periodization into discretization resulting in Equations ( 24) and ( 29) via the Fourier transform.And special cases lead to Equations ( 16) and (27).This can be seen by applying Equations ( 5) and ( 6) to PSF versions found in the literature.In the following we distinguish a general version, a classical version, a symmetrical version and, finally, a special version.
For simplicity, let φ ∈ S and T ∈ R + (one-dimensional case) for now.Recall that PSF are valid on φ ∈ S, see e.g., [3,9].In Section 7, we will extend these results to generalized functions f ∈ O C ′ ⊃ S and S ⊂ O M ∋ Ff (see also Appendix B).We note that Equation (11) corresponds to Equations ( 24) and ( 29) in Section 7.
A general version of Poisson's summation formula can be found for example as in Benedetto and Zimmermann [3] but also in [16,37,38,[41][42][43].Inserting Equations ( 5) and ( 6) and using Equation ( 9), this formula yields (29).A variation of (11) emerges if T = 1 , e.g., in Gröchenig [39].Now, let t ≡ 0, then Equation ( 11) reduces to the classical version of Poisson's summation formula found in [1][2][3][43][44][45] for example.Now, let T ≡ 1, then Equation ( 12) reduces to the symmetrical version of Poisson's summation formula where φ := Fφ, e.g., found in Feichtinger and Strohmer [9], also in [37,[46][47][48].We note that the latter two Equations ( 12) and ( 13), especially the last one, look most elegant and one might guess that there would be some loss of information from top to bottom.But surprisingly, the deduction could also be done from bottom to top, first via argument substitution x → x T , i.e., using Fourier transform rule F(φ(xT )) = (1/T ) (Fφ)(x/T ) and then using argument substitution x → t + x, i.e., using Fourier transform rule F(φ(t + x)) = (Fφ)(x) e 2πit•x .Hence, Equations (11), ( 12) and ( 13) can be transferred into one another for φ ∈ S. Now, let us assume that instead of φ ∈ S ⊂ S ′ , we indeed have some suitable f ∈ S ′ , let f ≡ δ, then Equation ( 11) is valid in S ′ (known result) and reduces to the special version of Poisson's summation formula also found very often in the literature, for example in [16,26,39,49] stating that the Fourier transform of a Dirac comb is a Dirac comb.Thus, Equation ( 11) with f ≡ δ yields (27).Additionally having T = 1 yields (16), found for example in [3,16,50,51].Note that in Equation ( 14) we already use identity . Such rules will be treated in the next section.Finally, by looking at the from bottom-to-top derivation above, we note that there is a certain degree of freedom on what side to move the Fourier transform.This allows us, in contrast to what is already done, to periodize Ff instead of f , resulting in the other two Equations ( 23) and ( 28) we will prove in Section 7.These formulas are, in contrast to Equations ( 24) and (29), very difficult to find in the literature.However, [52] is one of the rare sources where beside the classical Poisson Summation Equation (37.1) in [52] another dual formula Equation (37.2) in [52] can be found instead of only one; Equation (37.2) in [52] corresponds to Equations ( 23) and (28) and Equation (37.1) in [52] corresponds to Equations ( 24) and (29).

Calculation Rules
In order to present a concise proof of Theorem 1, we need to compile few important operator properties.We note that the Fourier transform of a Dirac comb within S ′ is commonly known (e.g., [16,39]) whereas the (asymmetric) duality of multiplication and convolution within S ′ is a more generalized case of commonly known (symmetric) versions (e.g., in S), it can be found in [21].
Dirac Comb Identity.The first, maybe surprising observation by looking at the two definitions above is that both operations generate a Dirac impulse train by either sampling the function that is constantly 1 or by periodically repeating δ, symbolically Moreover, 1 and δ are unitary elements with respect to multiplication and convolution in S ′ , respectively.They are also a Fourier pair, more precisely F1 = δ and Fδ = 1.Also note that the Shah symbol III in Equation ( 4) is called the sampling or replicating symbol in [28] because it can be used for either sampling (discretizing) or replicating (periodizing) a function.However, instead of III we now have either ⊥ ⊥ ⊥1 or, alternatively, △ △ △δ which allows us to switch from discretization to periodization or from periodization to discretization (needed in proof of Theorem 1).Dirac Comb Fourier Transform.We also need that F( ⊥ ⊥ ⊥1) = ⊥ ⊥ ⊥1 and, accordingly, Omitting Brackets.A very convenient property of both operators is that brackets can be omitted.For reasons of brevity, we move this proof to Appendix A.
In particular, with ⊥ ⊥ ⊥1 = III = △ △ △δ, additionally let us identify ⊥ ⊥ ⊥1 ≡ ⊥ ⊥ ⊥ and △ △ △δ ≡ △ △ △, with Equations ( 17) and ( 18) we obtain relationships published in Bracewell [28] It means that discretization is realized via multiplication and periodization is realized via convolution with III.Note that replacing III by operations ⊥ ⊥ ⊥, △ △ △ is the central proposal in this paper.Multiplication and Convolution.Recall that multiplication or convolution products do not necessarily exist in S ′ .We must rather restrict our operator arguments to suitable subspaces of S ′ (see operator definitions), i.e., to O M and O C ′ respectively.

Lemma 3 (Multiplication and Convolution
A detailed treatment of multiplication and convolution in S ′ including a proof of this lemma is given in [21], Théorème XV, p.124.Notice that O M can rarely be found in the literature ( [30,[53][54][55][56]) while O C ′ is better known.A description of interactions of O M with O C ′ , aside from [21], is even more seldom in the literature, e.g., [57,58] describing the Fourier transform as an algebra isomorphism

The Discretization-Periodization Theorem
For the sake of simplicity, we first consider only increments of T ν = 1 in all dimensions 1 ⩽ ν ⩽ n and then generalize these results in a corollary below.

Unitary Increments
Everything is prepared now to present a discretization -periodization theorem in S ′ .It consists of two formulas very closely related to Poisson's Summation Formula starting from discretization and from periodization, respectively.Compared to classical versions of Poisson's summation formula we simply hide one sum in ⊥ ⊥ ⊥ and the other in △ △ △, such that these formulas become better understandable.
Theorem 1 (Discretization vs. Periodization, unitary case).Let α ∈ O M and f ∈ O C ′ .Then Furthermore, recall that F( ⊥ ⊥ ⊥α) and △ △ △(F α) are the Fourier series F S of α according to Equation (7).It is, in other words, the Fourier transform of discrete functions F S .And the second equation corresponds to the Fourier transform of periodic functions F P (Figure 3).See one-page summaries: Figure 3.31 in [12] or Prolog II in [16], for example, embracing all four different Fourier transforms and Fourier series including their (either integral or discrete) inversion formulas.Proof.Formally, according to the calculation rules shown above the following equalities hold ′ is the unitary element with respect to convolution in Then we apply Equations ( 18), ( 22), ( 16), ( 15) and (17), in this order.Finally, let α := Ff ∈ O M , we use α = 1 • α because 1 ∈ O M is the unitary element with respect to multiplication in O M to show the second formula in an analogous manner.
To demonstrate the usefulness of the theorem, we use Equations ( 23), ( 24) and ( 15) to prove Equation ( 16) by and Finally, we use the theorem to conclude that Equation ( 8) is complemented having in the operations sense • : O M → S ′ .Hence, Fourier series F S either correspond to discretization followed by the Fourier transform or, equivalently, they correspond to the Fourier transform followed by periodization.Accordingly, for the second Poisson Summation Formula ( 24) we obtain in the operations sense

Arbitrary Increments
In our Equations ( 1), ( 3) and ( 4) as well as in Equations ( 16), ( 23) and ( 24) we used periodic functions with periods of T ν = 1 in all dimensions 1 ⩽ ν ⩽ n, which is convenient but not the general case.The following lemma helps to generalize our results to arbitrary periods and increments.
A proof of this lemma is straight forward from what is known in the literature (see for example [16,32,40,49]).We will omit it here for the sake of conciseness.Let us now see the theorem in the arbitrary increments case. where This is instantly verified replacing Equation ( 16) by Equation (27) in the above proof.

Conclusions
We have seen that, in the generalized functions sense, there is a (Fourier transform) duality between discrete and periodic functions.The theorem states that the generalized function being discretized must not grow faster than polynomials and the generalized function being periodized must rapidly decay to zero in order to obtain convergence on both sides of the equations in the generalized functions sense.Furthermore, it is shown that both formulas in this theorem are closely related to the classical Poisson Summation Formula.This means that these conditions will have impact on the validity of Poisson Summation Formula versions.We have, furthermore, seen that one formula in this theorem corresponds to Fourier series and the other formula corresponds to the Fourier transform of periodic functions.Finally, by choosing two similar symbols for discretization and periodization, we showed that an easy-to-use symbolic calculation scheme is obtained.

Outlook
The theorem above is indisputably connected to the classical sampling theorem, see, e.g., Woodward [1] or Benedetto and Zimmermann [3].An all-embracing theory including the Sampling Theorem will, however, require to have two more operations, one reversing discretization and one reversing periodization, i.e., interpolation of discrete functions and restriction of a periodic function to one period.Surprisingly, S ′ has been suitable to model discretization and periodization but will not be suitable to model their inverses.Further studies are required.for all test functions φ ∈ S, i.e., β ⊥ ⊥ ⊥ T α = ⊥ ⊥ ⊥ T (βα) in S ′ .Finally, the roles of α and β can be exchanged due to commutativity βα = αβ in O M .
Proof of Lemma 2. Periodization, △ △ △.We know that △ △ △ T g ∈ S ′ and with f ∈ O C ′ it follows that f * △ △ △ T g ∈ S ′ such that it can be applied to test functions φ ∈ S. Due to the convergence of △ △ △ T g in S ′ and the continuity of convolution * : O C ′ × S ′ → S ′ , (f, g) → f * g we may calculate term by term, we use Equations (A3), (A4), (A5), (A3), (A4), in this order, obtaining for all test functions φ ∈ S, i.e., f Finally, the roles of f and g can be exchanged due to commutativity The following definitions and calculation rules are needed for the above proof.
Definition A1.Let f ∈ S ′ and φ ∈ S. Then Proof.We may write The last expression is a function in S because with Ff ∈ O M and Fφ ∈ S it follows that Ff • Fφ ∈ S. Finally,  6) and (A1), exchanges in the sequence of translations and again Equation (A1), in this order.

B. Summary
Let us now have a look at the beauty of "sortedness" of functions in S ′ (Figure 4).It is this "sortedness" that led us from constant function 1 via the Gaussian function g(t) = e −πt 2 up to the discovery of δ.Also note that discrete functions are being injected on the smooth side of S ′ (left) whereas periodic functions are being injected on the rough side of S ′ (right).But both are actually located on the same side of S ′ (opposite of S) because left hand side and right hand side are being identical, recall that ⊥ ⊥ ⊥1 ≡ △ △ △δ.So, S ′ reminds us of a torus.

)
Note that, in contrast to α ∈ O M and f ∈ O C ′ , both formulas hold for Schwartz functions φ ∈ S understood as rapidly decreasing tempered distributions because φ ∈ S ⊂ O M ⊂ S ′ on the one hand, and φ ∈ S ⊂ F(O M ) = O C ′ ⊂ S ′ on the other.Hence, it proves Equation (10) in particular (see Appendix B).

Figure 4 .
Figure 4.The space of tempered distributions and important subspaces.