Sampling formulas involving differences in
shift-invariant subspaces: a unified approach

arXiv:1709.03893v1 [math.FA] 12 Sep 2017

Antonio G. Garcı́a∗ and Marı́a J. Muñoz-Bouzo†.

* Departamento de Matemáticas, Universidad Carlos III de Madrid, Avda. de la
Universidad 30, 28911 Leganés-Madrid, Spain.
† Departamento de Matemáticas Fundamentales, U.N.E.D., Senda del Rey 9, 28040
Madrid, Spain.
Abstract
Successive differences on a sequence of data help to discover some smoothness
features of this data. This was one of the main reasons for rewriting the classical
interpolation formula in terms of such data differences. The aim of this paper is
to mimic them to a sequence of regular samples of a function in a shift-invariant
subspace allowing its stable recovery. A suitable expression for the functions in
the shift-invariant subspace by means of an isomorphism with the L2 (0, 1) space
is the key to identify the simple pattern followed by the dual Riesz bases involved
in the derived formulas. The paper contains examples illustrating different nonexhaustive situations including also the two-dimensional case.

Keywords: Forward and backward differences; Averages; Shift-invariant subspaces;
Dual Riesz bases; Sampling formulas.
AMS: 42C15; 94A20; 97N50.

1

Statement of the problem

P
2
2
Consider the shift-invariant space Vϕ2 =
n∈Z an ϕ(t − n) : {an } ∈ ` (Z) of L (R)
2
generated by ϕ ∈ L (R). The goal of this work is to investigate the existence of stable
sampling formulas for any f ∈ Vϕ2 from differences (or averages) obtained from a regular
sequence of samples {f (a + n)}n∈Z . That is, the involved samples are p subsequences
taken from
 k
 k

∆+ f (a + pn) n∈Z ; and/or
∆− f (a + pn) n∈Z ,
(1)
f (a + pn) n∈Z ;
k−1
where ∆k+ f (a + pn) := ∆k−1
+ f (a + pn + 1) − ∆+ f (a + pn) denotes the k-th forward
k−1
difference operator, ∆k− f (a + pn) := ∆k−1
− f (a + pn) − ∆− f (a + pn − 1) denotes the
k-th backward difference operator, 1 ≤ k ≤ p − 1, and a is a fixed number in [0, 1);
∗
†

E-mail:agarcia@math.uc3m.es
E-mail:mjmunoz@mat.uned.es

1


here ∆0± = I. Successive differences on samples f (a + n) n∈Z can help us to discover
some smoothness features on the original function f .
It is a well-known that, under some appropriate hypotheses (see Section 2.2 below),
a Shannon-type sampling formula as
X
f (t) =
f (a + n) Sa (t − n) , t ∈ R ,
(2)
n∈Z


holds in Vϕ2 for some sampling function Sa ∈ Vϕ2 such that the sequence Sa (t − n) n∈Z
is a Riesz basis for Vϕ2 . As a consequence, the above sampling formula is an interpolation
formula: Indeed, for each n ∈ Z, the function Sa satisfies the interpolation property
Sa (a + n) = δn .
The sampling formula (2) was first proved by Walter in [22]; since then, sampling
theory in shift-invariant subspaces has been largely studied; see, for instance, Refs.
[1, 3, 7, 10, 15, 16, 20, 23]. For the general theory of shift-invariant subspaces and their
applications, see, for instance, Refs. [2, 4, 5, 6, 18, 21].
Roughly speaking, the aim of the paper is to rewrite formula (2) by replacing some
samples by any other information given by differences, averages, etc. Let us consider
an easy example. Initially, we write formula (2) as
X
f (t) =
f (a + 2n) Sa (t − 2n) + f (a + 2n + 1) Sa (t − 2n − 1) , t ∈ R .
n∈Z

Then, adding and subtracting f (a + 2n) Sa (t − 2n − 1) in each summand yields
X


f (t) =
f (a+2n) Sa (t−2n)+Sa (t−2n−1) +∆+ f (a+2n) Sa (t−2n−1) , t ∈ R .
n∈Z



The stable recovery of any f ∈ Vϕ2 from the samples f (a+2n) n∈Z ∪ ∆+ f (a+2n)
follows from the fact that the sequence


Sa (t − 2n) + Sa (t − 2n − 1) n∈Z ∪ Sa (t − 2n − 1) n∈Z

n∈Z

turns out to be a Riesz basis for Vϕ2 , or equivalently, there exist constants 0 < A ≤ B
such that
X
Akf k2 ≤
|f (a + 2n)|2 + |∆+ f (a + 2n)|2 ≤ Bkf k2 for each f ∈ Vϕ2 .
n∈Z

In this paper we derive, in a systematic way, stable sampling formulas involving
forward and/or backward differences of arbitrary order or averages of the sequence
of samples {f (a + n)}n∈Z . The used technique will allow to deal also with central
differences and averages; this is done in Section 3. Thus we obtain a sort of counterpart
to the family of finite interpolation formulas in the infinite setting of sampling in shiftinvariant subspaces of L2 (R). The existence of an isomorphism Tϕ between L2 (0, 1)
and Vϕ2 and an associated representation of the functions in Vϕ2 (see (4) below) make
the mathematical technique very simple: we identify the available samples as inner
products with respect to a suitable Riesz basis in L2 (0, 1) for which we obtain its
corresponding dual Riesz basis; then we derive the corresponding sampling formula in
Vϕ2 by using Tϕ . Finally, the two-dimensional case is also treated in Section 4.
2

2

Sampling theorem in a shift-invariant subspaces of L2 (R)

For the sake of completeness we include the preliminaries on sampling in shift-invariant
subspaces needed in the sequel.

2.1

Some preliminaries on shift-invariant subspaces

Assume that for ϕ ∈ L2 (R) the sequence {ϕ(t − n)}n∈Z is a Riesz sequence for L2 (R),
i.e., a Riesz basis for its closed span. A Riesz basis for a (separable) Hilbert space is
the image of an orthonormal basis by means of an invertible bounded operator. Thus,
the (principal) shift-invariant subspace Vϕ2 := span{ϕ(t − n)}n∈Z can be described as
Vϕ2 =

X

an ϕ(t − n) : {an } ∈ `2 (Z) .

n∈Z

We also assume that the generator ϕ is continuous on R and the 1-periodic function
P
2
2
n∈Z |ϕ(t − n)| is bounded on [0, 1]; this is equivalently to assume that Vϕ is a shiftinvariant subspace of continuous functions on R [23]. As a consequence, the functions
2
in Vϕ2 are
P continuous on R and any f ∈ Vϕ is described on R as the pointwise sum
f (t) = n∈Z an ϕ(t − n), t ∈ R.
On the other hand, the shift-invariant space Vϕ2 is the image of L2 (0, 1) by means
of the isomorphism:
Tϕ :

L2 (0, 1)
−→ Vϕ2
{e−2πinx }n∈Z 7−→ {ϕ(t − n)}n∈Z ,

which maps the orthonormal basis {e−2πinw }n∈Z for L2 (0, 1) onto the Riesz basis {ϕ(t−
n)}n∈Z for Vϕ2 . This isomorphism Tϕ satisfies, for each F ∈ L2 (0, 1) and n ∈ Z, the
following shifting property


Tϕ F (x) e−2πinx (t) = Tϕ [F ](t − n) , t ∈ R .
(3)
After some straightforward calculations (see [9, 12] for the details), any f ∈ Vϕ2 can be
expressed as
f (t) = F, Kt L2 (0,1) , t ∈ R ,
(4)
P
where f = Tϕ F and Kt (x) = n∈Z ϕ(t − n) e−2πinx ∈ L2 (0, 1). Notice that Kt (x) =
P
−2πinx is just the Zak transform of the
Zϕ(t, x), where Zϕ(t, x) :=
n∈Z ϕ(t + n)e
function ϕ (see [8] for properties and uses of the Zak transform). For t ∈ R and m ∈ Z,
the following property holds: Kt+m (x) = e−2πimx Kt (x).
From expression (4) we easily deduce that evaluation functionals are bounded on
Vϕ2 , i.e., it is a reproducing kernel Hilbert space (RKHS). Thus, convergence in the
L2 (R)-norm implies
pointwise convergence which is uniform on R since the function
P
t 7→ kKt k2 = n∈Z |ϕ(t − n)|2 is assumed to be bounded on R. As a consequence,
from now on we are only concerned with L2 (R)-norm convergence. See more details in
[12, 13].

3

2.2

Shannon-type sampling formula in Vϕ2

Having in mind expression (4), for a ∈ [0, 1) fixed and n ∈ Z we have
f (a + n) = hF, Ka+n iL2 (0,1) = hF, e−2πinx Ka iL2 (0,1) ,

F = Tϕ−1 f .

Thus, the sequence of samples {f
 (a + n)}n∈Z forms a stable sampling set for any f ∈
Vϕ2 if and only if the sequence e−2πinx Ka (x) n∈Z is a Riesz basis for L2 (0, 1). This
occurs if and only if the inequalities 0 < kKa k0 ≤ kKa k∞ < ∞ hold, where kKa k0 :=
ess inf x∈(0,1) |Ka (x)| and kKa k∞ := ess supx∈(0,1) |Ka (x)|. Moreover, its dual Riesz basis

is e−2πinx /Ka (x) n∈Z (see [9, 12] for the details). Therefore, for each F ∈ L2 (0, 1) we
have the expansion
F =

X

hF, e−2πinx Ka i

n∈Z

e−2πinx
Ka (x)

=

X
n∈Z

f (a + n)

e−2πinx
Ka (x)

in L2 (0, 1) .

Applying the isomorphism Tϕ and its shifting property we obtain in Vϕ2 the Shannontype sampling formula
X
f (t) =
f (a + n) Sa (t − n) , t ∈ R ,
n∈Z


where Sa := Tϕ 1/Ka ∈ Vϕ2 . The convergence of the series is absolute due to the
unconditional convergence of a Riesz basis expansion and it is also uniform on R since
the shift-invariant space Vϕ2 is a RKHS.
Some well-known examples
A famous example of shift-invariant subspace in L2 (R) is the Paley-Wiener space of
band-limited functions to an interval, say [−π, π], i.e.,

P Wπ := f ∈ L2 (R) : supp fb ⊆ [−π, π] .
2 whose generator ϕ is the sinc funcIt coincides with the shift-invariant subspace Vsinc
tion, i.e., sinc t := sin πt/πt. In P Wπ , for any a ∈ [0, 1), the sampling formula

f (t) =

X
n∈Z

f (n + a)

sin π(t − n − a)
,
π(t − n − a)

t ∈ R,

holds. In particular, taking a = 0 we recover the famous Shannon’s sampling formula.
Although Shannon’s sampling theory has had an enormous impact, it has a number
of problems as pointed out in [21]. From a practical point of view, the most important
examples of shift-invariant spaces Vϕ2 are those generated by B-splines. Consider ϕ :=
Nm where Nm is the B-spline of order m − 1, i.e., Nm := N1 ∗ N1 ∗ · · · ∗ N1 (m times)
where N1 := χ[0,1] denotes the characteristic function of the interval [0, 1]. It is known
that the sequence Nm (t − n) n∈Z is a Riesz basis for VN2m (see, for instance, [8]). For
the quadratic and cubic splines the following sampling formulas hold:

4

Figure 1: Quadratic S1/2 and cubic S0 sampling funtions
(1) Quadratic Spline N3 . Here, taking a = 1/2, for any f ∈ VN23 we obtain the sampling
formula
X
1
f (t) =
f (n + ) S1/2 (t − n) , t ∈ R ,
2
n∈Z
√
√ P
where S1/2 (t) = 2 n∈Z (2 2 − 3)|n+1| N3 (t − n) (see Fig.1). For a = 0 the sampling
result fails since kK0 k0 = 0; see [12] for the details.
(2) Cubic Spline N4 . Here, taking a = 0, for any f ∈ VN24 we obtain the sampling
formula
X
f (t) =
f (n) S0 (t − n) , t ∈ R ,
n∈Z

√ P
√
where S0 (t) = 3 n∈Z (−1)n (2 − 3)|n| N4 (t − n + 2) (see Fig.1). See [12] for the
details. In Ref. [19] one can find more information on cardinal spline interpolation.

2.3

The involved Riesz bases in L2 (0, 1) and their duals

Having in mind expression (4), for the forward difference operator ∆+ of f ∈ Vϕ2 , i.e.,
∆+ f (t) := f (t + 1) − f (t), t ∈ R, we have
∆+ f (t) = F, Kt+1 − Kt

L2 (0,1)

= F, (e−2πix − 1) Kt

L2 (0,1)

,

t ∈ R.

In general, for the k-th forward difference operator we have the expression
∆k+ f (t) = F, (e−2πix − 1)k Kt

L2 (0,1)

,

t ∈ R.

(5)

For the backward difference operator ∆− of f ∈ Vϕ2 , i.e., ∆− f (t) := f (t) − f (t − 1),
t ∈ R, we have
∆− f (t) = F, Kt − Kt−1

L2 (0,1)

= F, (1 − e2πix ) Kt

L2 (0,1)

,

t ∈ R.

In general, for the k-th backward difference operator we have the expression
∆k− f (t) = F, (1 − e2πix )k Kt

5

L2 (0,1)

,

t ∈ R.

(6)

Similar expressions can be also obtained from iterates of the central difference operator
∆0 , i.e., ∆0 f (t) := f (t + 1) − f (t − 1), or from the average operators µ+ , µ− or µ0
defined as:
1
1
1
µ+ f (t) := [f (t+1)+f (t)] , µ− f (t) := [f (t)+f (t−1)] , µ0 f (t) := [f (t+1)+f (t−1)] .
2
2
2
All these cases
P can be englobed by considering generalized finite differences defined as
∆am,n f (t) := N
k=M ak f (t + k), where a = (ak ) ⊂ C and M, N ∈ Z.
Going back to our previous example in Section 1, we have for the samples f (a + 2n)
and ∆+ f (a + 2n), n ∈ Z:
f (a + 2n) = F, e−2πi(2n)x Ka
∆+ f (a + 2n) =

,

L2 (0,1)
F, e−2πi(2n+1)x Ka − e−2πi(2n)x Ka L2 (0,1) .

In other words, denoting xn := e−2πinx Ka (x), n ∈ Z, the questions are whether the
sequence {x2n }n∈Z ∪ {x2n+1 − x2n }n∈Z is a Riesz basis for L2 (0, 1) and how to obtain, in
the affirmative case, its dual basis. As it was mentioned in the introduction section, this
is true: Indeed, denoting yn := e−2πinx /Ka (x), n ∈ Z, the dual Riesz basis of {xn }n∈Z ,
the dual Riesz basis of {x2n }n∈Z ∪ {x2n+1 − x2n }n∈Z is {y2n + y2n+1 }n∈Z ∪ {y2n+1 }n∈Z .
A general answer which involves the aforementioned generalized finite differences is
given in the following result:
Lemma 1. Let {xn }n∈Z be a Riesz basis for a separable Hilbert space H with dual
Riesz basis {yn }n∈Z . For a fixed p ∈ N, consider
a partition {x1n }n∈Z ∪ {x2n }n∈Z ∪

· · · ∪ {xpn }n∈Z of {xn }n∈Z and let M = ajk ∈ Cp×p be a matrix of scalars. Then, the
new sequence defined by using the rows of M as
p
X

a1k xkn

k=1

∪
n∈Z

p
X

a2k xkn

k=1

∪ ··· ∪
n∈Z

p
X

apk xkn

n∈Z

(7)

k=1

is a Riesz basis for H if and only if det M 6= 0. In this case, its dual Riesz basis is
given by
p
p
p
X
X
X
bk1 ykn n∈Z ∪
bk2 ykn n∈Z ∪ · · · ∪
bkp ykn n∈Z
k=1

k=1

k=1

>

where b1k b2k . . . bpk
denotes the k-th column of the inverse matrix M −1 , k =
1, 2, . . . , p. Moreover, if a sequence as (7) with det M 6= 0 is a Riesz basis for H
then the sequence {xn }n∈Z is also a Riesz basis for H.
P
P
P
Proof. If det M = 0 then { pk=1 a1k xkm , pk=1 a2k xkn , . . . , pk=1 apk xkn } is a linearly
dependent system for each m ∈ Z, so that the sequence defined in (7) is not a Riesz
basis.
Conversely, assume that det M 6= 0. Denoting the new sequence defined in (7) as
{zn }n∈Z and its partition as {z1n }n∈Z ∪ {z2n }n∈Z ∪ · · · ∪ {zpn }n∈Z , we rewrite (7) using
the symbolic notation
>
>
M x1n , x2n , . . . , xpn = z1n , z2n , . . . , zpn , n ∈ Z .
6

It is easily deduced what we also symbolically write as
x1n , x2n , . . . , xpn

>

= M −1 z1n , z2n , . . . , zpn

>

,

n ∈ Z.

Thus span{zn }n∈Z = span{xn }n∈Z and consequently the sequence {zn }n∈Z is complete
in H.
P
In addition, for every finite scalar sequence {cm } and for every {zm } = { pk=1 ajm k xknm }
with nm ∈ Z and jm ∈ {1, 2, . . . , p}, one has
X

cm zm

2

=

X

m

cm

m

= B

X

p
X

ajm k xknm

2

≤B

p
XX

X

|ajm k |2 ≤ BD

|cm |

|cm |2

m

k=1

where B is an upper Riesz bound of {xn }n∈Z and D = maxpj=1
X

cm zm

2

=

X

m

cm

m

= A

X

p
X

ajm k xknm

k=1
p
X
2

|cm |

m

|cm ajm k |2

m k=1

k=1
p
X
2

m





2

≥A

Pp

2
k=1 |ajk | .

p
XX

Similarly,

|cm ajm k |2

m k=1

|ajm k |2 ≥ AC

X

|cm |2

m

k=1

P
where A is a lower Riesz bound of the Riesz basis {xn }n∈Z and C = minpj=1 pk=1 |ajk |2 .
Since det M 6= 0, then C > 0 and consequently, the sequence {zn }n∈Z satisfies the Riesz
sequence condition (see [8, Theorem 3.6.6])
X
X
X
2
A0
|cm |2 ≤
cm zm ≤ B 0
|cm |2 .
m

m

m

with constants A0 = AC > 0 and B 0 =
> 0.
PBD
p
0
0
Finally, since hxq , yq i = δq,q and k=1 ajk bkj 0 = δj,j 0
p
X
k=1

ajk xkn ,

p
X

bk0 j 0 yk0 m

k0 =1

 Pp
proving that
k=1 ajk xkn
pair of Riesz bases for H.

p
X

(
1
=
ajk bk0 j 0 hxkn , yk0 m i =
0
k,k0 =1
j=1,2,...,p, n∈Z

and

 Pp

if n 6= m and j = j 0
otherwise

k=1 bkj ykn j=1,2,...p, n∈Z

are a dual

Observe that in the case of one forward difference, the partition
 of the
 Riesz basis
1 0
{xn }n∈Z is {x2n }n∈Z ∪ {x2n+1 }n∈Z , the associated matrix is M =
; its inverse
−1 1


1 0
gives the dual Riesz basis {y2n + y2n+1 }n∈Z ∪ {y2n+1 }n∈Z .
1 1
Although we will not derive overcomplete expansions in the sequel, it is worth to
mention the frame counterpart to the above lemma. Under the circumstances of Lemma
1 we have:
7


Lemma 2. Let M = ajk ∈ Cq×p be a matrix
Ppof scalars such that q > p and rank M =
p. By using the rows of M we define zjn = k=1 ajk xkn , 1 ≤ j ≤ q and n ∈ Z. Then
the sequence {z1n }n∈Z ∪ {z2n }n∈Z ∪ · · · ∪ {zqn }n∈Z is a frame for H. Moreover, a dual
frame {e
z1n }n∈Z ∪{e
zqn }n∈ZP
is obtained from the columns of a left inverse
 z2n }n∈Z ∪· · ·∪{e
matrix N = bkj ∈ Cp×q of M as zejn = pk=1 bkj ykn , 1 ≤ j ≤ q and n ∈ Z.
Proof. Clearly, both sequences
are
P
P Bessel sequences. Following [8, Lemma 5.6.2], it
suffices to prove that x = n∈Z qj=1 hx, zejn i zjn for each x ∈ H. Indeed,
q
XX

x,

n∈Z j=1

=

X

p
X
k=1

p
X

n∈Z k,k0 =1

since

Pq

j=1 bkj ajk0

bkj ykn

p
X

ajk xkn =

k=1

δk,k0 hx, ykn ixk0 n =

p
X X

n∈Z k,k0 =1
p
XX

q
X


bkj ajk0 hx, ykn ixk0 n

j=1

hx, ykn ixkn =

n∈Z k=1

X

hx, yn ixn = x ,

n∈Z

= δk,k0 , and {xn }n∈Z and {yn }n∈Z form a pair of dual Riesz bases.

Observe
that the
as

 set of all left inverse matrices N of M can be expressed
−1N =∗
†
†
∗
Iq −M M where U denotes an arbitrary p×q matrix and M = M M
M
is the Moore-Penrose pseudo-inverse of the matrix M .
Observe that the proof of the above lemma is also valid in case the sequences
{xn }n∈Z and {yn }n∈Z form a pair of dual frames for H.

M † +U

3

The resulting sampling formulas

Along this section we fix a ∈ [0, 1) such that 0 < kKa k0 ≤ kKa k∞ < ∞, and we will
denote
xn := e−2πinx Ka (x) and yn := e−2πinx /Ka (x) , n ∈ Z ,
the associated pair of dual Riesz bases for L2 (0, 1) which appeared in Section 2.2. The
general procedure is easy: For any f ∈ Vϕ2 , let F ∈ L2 (0, 1) be such that Tϕ F = f ; once
we identify the Riesz basis for L2 (0, 1) associated with successive differences considered
on the data samples, we should compute its dual Riesz basis by using Lemma 1. Then
we expand F with respect to this dual basis in L2 (0, 1). Finally, by applying the
isomorphism Tϕ we obtain a sampling formula in Vϕ2 with L2 (R)-convergence; as Vϕ2
is a RKHS we obtain the pointwise convergence properties of the sampling expansion:
the convergence will be uniform on R. It will be also absolute due the unconditional
convergence of an expansion with respect to a Riesz basis.
Next we will exhibit a variety of examples involving different types of samples obtained from forward, backward or central difference operators, or from average operators.

8

3.1

The case involving forward differences

Using first and second forward differences
For the samples f (a + 3n), ∆+ f (a + 3n), and ∆2+ f (a + 3n), n ∈ Z, we have (see (5)):
f (a + 3n) = F, e−2πi(3n)x Ka
∆+ f (a + 3n) =
∆2+ f (a + 3n) =

,

L2 (0,1)
−2πi(3n+1)x
F, e
Ka − e−2πi(3n)x Ka L2 (0,1) ,
F, e−2πi(3n+2)x Ka − 2e−2πi(3n+1)x Ka + e−2πi(3n)x Ka L2 (0,1) .

If we consider the partition {x3n }n∈Z ∪ {x3n+1 }n∈Z ∪ {x3n+2 }n∈Z of {xn }n∈Z , according
to Lemma 1 the sequence
{x3n }n∈Z ∪ {x3n+1 − x3n }n∈Z ∪ {x3n+2 − 2x3n+1 + x3n }n∈Z
is a Riesz basis for L2 (0, 1) and its dual Riesz basis is
{y3n + y3n+1 + y3n+2 }n∈Z ∪ {y3n+1 + 2y3n+2 }n∈Z ∪ {y3n+2 }n∈Z ,
since


−1 
1 0 0
1
0 0
−1 1 0 = 1 1 0
1 2 1
1 −2 1


Expanding F with respect to the above dual basis we obtain
Xn




F =
f (a + 3n) y3n + y3n+1 + y3n+2 + ∆+ f (a + 3n) y3n+1 + 2y3n+2 +
n∈Z

+ ∆2+ f (a + 3n) y3n+2

o

in L2 (0, 1) .

Applying the isomorphism Tϕ and its shifting property (3), we obtain in Vϕ2 the sampling
formula

f (t) =

Xn


f (a + 3n) Sa (t − 3n) + Sa (t − 3n − 1) + Sa (t − 3n − 2)
n∈Z



+ ∆+ f (a + 3n) Sa (t − 3n − 1) + 2Sa (t − 3n − 2)
o
+ ∆2+ f (a + 3n) Sa (t − 3n − 2) , t ∈ R .



Observe that, due to the shifting property (3), we have Tϕ ypn+k = Sa t − (pn + k) .

9

Using p − 1 forward differences
Proceeding as in the above example, for the samples f (a + pn), ∆+ f (a + pn), . . . , and
∆p−1
+ f (a + pn), n ∈ Z, we have (see (5)):
f (a + pn) = F, e−2πi(pn)x Ka

L2 (0,1)

,

∆+ f (a + pn) = F, e−2πi(pn+1)x Ka − e−2πi(pn)x Ka

L2 (0,1)

,

∆2+ f (a + pn) = F, e−2πi(pn+2)x Ka − 2e−2πi(pn+1)x Ka + e−2πi(pn)x Ka

L2 (0,1)

······
∆p−1
+ f (a

+ pn) = F, e

−2πi(pn+p−1)x

+ · · · + (−1)p−2



−2πi(pn+p−2)x



p − 1 −2πi(pn+p−3)x
Ka +
e
Ka
2

Ka − (p − 1)e

p − 1 −2πi(pn+1)x
e
Ka + (−1)p−1 e−2πi(pn)x Ka
p−2

L2 (0,1)

If we consider the partition {xpn }n∈Z ∪ {xpn+1 }n∈Z ∪ · · · ∪ {xpn+p−1 }n∈Z of {xn }n∈Z ,
according to Lemma 1 the sequence defined by using the rows of the lower triangular
matrix:


1

 −1
1



 1
2
− 1
1




3
3
−
1

 −1
1
2


.
.
.
.
.


..
..
..
..
..







p−2
p−2
p−2
p−2
p−3
p−4
p−5
(−1)
(−1)
(−1)
(−1)
.
.
.
1


1 
2 
3 

p−1
p−4 p−1
p−3 p−1
1
.
.
.
−
(−1)
(−1)
(−1)p−1 (−1)p−2 p−1
p−2
3
2
1
is a Riesz basis for L2 (0, 1) and its dual Riesz basis and, consequently, the sampling
functions are obtained from the columns of its inverse which is:


1
1

1


1

2
1


1



3
3
1
1


1
2
.

..
..
..
..
 ..

.
.
.
.







p−2
p−2
p−2
.
.
.
1
1

1 
2 
3 

p−1
p−1
p−1
1 p−1
.
.
.
1
1
2
3
p−2
The resulting sampling formula in Vϕ2 reads:

f (t) =

Xn
f (a + pn) Ta,1 (t − pn) + ∆+ f (a + pn) Ta,2 (t − pn)
n∈Z

o
+ ∆2+ f (a + pn) Ta,3 (t − pn) + · · · + ∆p−1
+ f (a + pn) Ta,p (t − pn) ,
where, for 1 ≤ k ≤ p, the sampling function Ta,k is given by
10

t ∈ R,






k−1
k
Ta,k (t) =
Sa (t − (k − 1)) +
Sa (t − k)+
k−1
k−1




k+1
p−1
+
Sa (t − (k + 1)) + · · · +
Sa (t − (p − 1)) ,
k−1
k−1

t ∈ R.

Observe that the sampling function Ta,k has p − k + 1 summands, 1 ≤ k ≤ p.

3.2

The general case involving p − 1 backward differences

Here, for the samples f (a + pn), ∆− f (a + pn), . . . , and ∆p−1
− f (a + pn), n ∈ Z, we have
(see (6)):
f (a + pn) = F, e−2πi(pn)x Ka
∆− f (a + pn) = F, e

−2πi(pn)x

L2 (0,1)

,

−2πi(pn−1)x

Ka − e

Ka

L2 (0,1)

,

∆2− f (a + pn) = F, e−2πi(pn)x Ka − 2e−2πi(pn−1)x Ka + e−2πi(pn−2)x Ka

L2 (0,1)

······


p − 1 −2πi(pn−2)x
−2πi(pn)x
−2πi(pn−1)x
∆p−1
f
(a
+
pn)
=
F,
e
K
−
(p
−
1)e
K
+
e
Ka
a
a
−
2


p − 1 −2πi(pn−p+2)x
e
Ka + (−1)p−1 e−2πi(pn−p+1)x Ka
+ · · · + (−1)p−2
p−2

L2 (0,1)

If we consider the partition {xpn−p+1 }n∈Z ∪{xpn−p+2 }n∈Z ∪· · ·∪{xpn−1 }n∈Z ∪{xpn }n∈Z
of {xn }n∈Z , according to Lemma 1 the sequence defined by using the rows of the upper
triangular matrix:



 p−1
 
p−1
p−1
p−3 p−1
(−1)p−1 (−1)p−2 p−1
(−1)
.
.
.
−
−
1
2 
3 
2 
1  1


p−2
p−2
p−2
p−2
p−2
p−3
(−1)
(−1)
... − 3 

1
2  − 1  1


p−3
(−1)p−3
. . . − p−3
− p−3
1

3
2
1

..
..
..
.. 
..


.
.
.
. 
.



3
− 31 1
−1

2



1
1
− 21



−1
1
1
is a Riesz basis for L2 (0, 1) and its dual Riesz basis and, consequently, the sampling
functions are obtained from the columns of its inverse which here coincides with itself.
Therefore, writing from the last column to the first one, the resulting sampling formula
in Vϕ2 reads:

f (t) =

Xn
f (a + pn) Tea,1 (t − pn) + ∆− f (a + pn) Tea,2 (t − pn)
n∈Z

o
ea,p (t − pn) ,
+ ∆2− f (a + pn) Tea,3 (t − pn) + · · · + ∆p−1
f
(a
+
pn)
T
−

where, for 1 ≤ k ≤ p, the sampling function Tea,k is given by
11

t ∈ R,



h p − 1 
p−2
k−1
e
Sa (t + p − 1) +
Sa (t + p − 2)+
Ta,k (t) = (−1)
k−1
k−1




i
p−3
k−1
+
Sa (t + p − 3) + · · · +
Sa (t + k − 1) , t ∈ R .
k−1
k−1
Observe that the sampling function Tea,k has p − k + 1 summands, 1 ≤ k ≤ p.

3.3

Two representative examples involving forward and backward differences

Using first forward and backward differences
For the samples f (a + 3n), ∆− f (a + 3n), and ∆+ f (a + 3n), n ∈ Z, we have:
∆− f (a + 3n) = F, e−2πi(3n)x Ka − e−2πi(3n−1)x Ka
f (a + 3n) = F, e−2πi(3n)x Ka
∆+ f (a + 3n) =

L2 (0,1)

,

,

L2 (0,1)
−2πi(3n+1)x
F, e
Ka − e−2πi(3n)x Ka L2 (0,1) .

If we consider the partition {x3n−1 }n∈Z ∪ {x3n }n∈Z ∪ {x3n+1 }n∈Z of {xn }n∈Z , according
to Lemma 1 the sequence
{−x3n−1 + x3n }n∈Z ∪ {x3n }n∈Z ∪ {−x3n + x3n+1 }n∈Z
is a Riesz basis for L2 (0, 1) and its dual Riesz basis is
{−y3n−1 }n∈Z ∪ {y3n−1 + y3n + y3n+1 }n∈Z ∪ {y3n+1 }n∈Z
since

−1 

−1 1 0
−1 1 0
0
1 0 =  0 1 0 
0 −1 1
0 1 1


The resulting sampling formula in Vϕ2 reads:

f (t) =

Xn

− ∆− f (a + 3n) Sa (t − 3n + 1)

n∈Z



+ f (a + 3n) Sa (t − 3n + 1) + Sa (t − 3n) + Sa (t − 3n − 1)
o
+ ∆+ f (a + 3n) Sa (t − 3n − 1) , t ∈ R .

12

Using first backward and first and second forward differences
For the samples f (a + 4n), ∆+ f (a + 4n), ∆2+ f (a + 4n), and ∆− f (a + 4n), n ∈ Z, we
have (see (5) and (6)):
∆− f (a + 4n) = F, e−2πi(4n)x Ka − e−2πi(4n−1)x Ka
f (a + 4n) = F, e−2πi(4n)x Ka

L2 (0,1)

L2 (0,1)

,

L2 (0,1)

,

,

∆+ f (a + 4n) = F, e−2πi(4n+1)x Ka − e−2πi(4n)x Ka

∆2+ f (a + 4n) = F, e−2πi(4n+2)x Ka − 2e−2πi(4n+1)x Ka + e−2πi(4n)x Ka

L2 (0,1)

.

If we consider the partition {x4n−1 }n∈Z ∪ {x4n }n∈Z ∪ {x4n+1 }n∈Z ∪ {x4n+2 }n∈Z of
{xn }n∈Z , according to Lemma 1 the sequence
{−x4n−1 + x4n }n∈Z ∪ {x4n }n∈Z ∪ {−x4n + x4n+1 }n∈Z ∪ {x4n − 2x4n+1 + x4n+2 }n∈Z
is a Riesz basis for L2 (0, 1) and its dual Riesz basis is
{−y4n−1 }n∈Z ∪ {y4n−1 + y4n + y4n+1 + y4n+2 }n∈Z ∪ {y4n+1 + 2y4n+2 }n∈Z ∪ {y4n+2 }n∈Z
since

−1 1
0
0
1
0

 0 −1 1
0
1 −2

−1 
0
−1


0
0
=
0
0
1
0

1
1
1
1

0
0
1
2


0
0

0
1

The resulting sampling formula in Vϕ2 reads:

f (t) =

Xn


− ∆− f (a + 4n)Sa (t − 4n + 1) + f (a + 4n) Sa (t − 4n + 1) + Sa (t − 4n)

n∈Z


+ Sa (t − 4n − 1) + Sa (t − 4n − 2)


+ ∆+ f (a + 4n) Sa (t − 4n − 1) + 2Sa (t − 4n − 2)
o
+ ∆2+ f (a + 4n) Sa (t − 4n − 2) , t ∈ R .

From the above example and Sections 3.1 and 3.2 the general case involving p
forward and q − 1 backward differences with sampling period p + q easily comes out.

3.4

Some examples involving averages and differences

Using forward averages and differences
For the samples µ+ f (a + 2n) and ∆+ f (a + 2n), n ∈ Z, we have:
1
µ+ f (a + 2n) = F, (e−2πi(2n+1)x Ka + e−2πi(2n)x Ka L2 (0,1) ,
2
∆+ f (a + 2n) = F, e−2πi(2n+1)x Ka − e−2πi(2n)x Ka L2 (0,1) .
13

If we consider the partition {x2n }n∈Z ∪ {x2n+1 }n∈Z of {xn }n∈Z , according to Lemma 1
the sequence
1
1
{ x2n + x2n+1 }n∈Z ∪ {−x2n + x2n+1 }n∈Z
2
2
is a Riesz basis for L2 (0, 1) and its dual Riesz basis is
1
1
{y2n + y2n+1 }n∈Z ∪ {− y2n + y2n+1 }n∈Z
2
2
since


−1 

1/2 1/2
1 −1/2
=
−1
1
1 1/2

The resulting sampling formula in Vϕ2 reads:

f (t) =

Xn


µ+ f (a + 2n) Sa (t − 2n) + Sa (t − 2n − 1)
n∈Z

+ ∆+ f (a + 2n)

o
1
− Sa (t − 2n) + Sa (t − 2n − 1) ,
2

t ∈ R.

Using central averages and differences
For the samples f (a + 3n), µ0 f (a + 3n), and ∆0 f (a + 3n), n ∈ Z, we have:
1
µ0 f (a + 3n) = F, (e−2πi(3n+1)x Ka + e−2πi(3n−1)x Ka
2
f (a + 3n) = F, e−2πi(3n)x Ka L2 (0,1) ,
∆0 f (a + 3n) = F, e−2πi(3n+1)x Ka − e−2πi(3n−1)x Ka

L2 (0,1)

L2 (0,1)

,

.

If we consider the partition {x3n−1 }n∈Z ∪ {x3n }n∈Z ∪ {x3n+1 }n∈Z of {xn }n∈Z , according
to Lemma 1 the sequence
1
1
{x3n }n∈Z ∪ { x3n−1 + x3n+1 }n∈Z ∪ {−x3n−1 + x3n+1 }n∈Z
2
2
is a Riesz basis for L2 (0, 1) and its dual Riesz basis is
1
1
{y3n }n∈Z ∪ {y3n−1 + y3n+1 }n∈Z ∪ {− y3n−1 + y3n+1 }n∈Z
2
2
since


−1 
0 1 0
0 1 −1/2
1/2 0 1/2 = 1 0
0 
−1 0 1
0 1 1/2


The resulting sampling formula in Vϕ2 reads:
14

f (t) =

Xh


f (a + 3n) Sa (t − 3n) + µ0 f (a + 3n) Sa (t − 3n + 1) + Sa (t − 3n − 1)
n∈Z

+ ∆0 f (a + 3n)

4

i
1
− Sa (t − 3n + 1) + Sa (t − 3n − 1) ,
2

t ∈ R.

Sampling results for the two-dimensional case

Firstly, we recall the usual sampling formulas in the two-dimensional setting:

4.1

Sampling formulas for two-dimensional shift-invariant subspaces

The general two-dimensional case

Consider a generator Φ ∈ L2 (R2 ) such that the sequence Φ(t − n, s − m) n,m∈Z forms
a Riesz sequence for L2 (R2 ). Thus,
n X
o
anm Φ(t − n, s − m) : {anm } ∈ `2 (Z2 ) .
VΦ2 := span{Φ(t − n, s − m)}n,m∈Z =
n,m∈Z

P
Assuming that Φ is continuous in R2 and that n,m∈Z |Φ(t − n, s − m)|2 is bounded in
[0, 1]2 one can easily deduce the same results that in Section 2 for the one dimensional
case. For instance, the isomorphism TΦ : L2 (0, 1)2 −→ VΦ2 which maps the orthonormal
basis {e−2πinx e−2πimy }n,m∈Z onto the Riesz basis {Φ(t − n, s − m)}n,m∈Z satisfies for
any F ∈ L2 (0, 1)2 the shifting property


TΦ e−2πinx e−2πimy F (x, y) (t, s) = TΦ [F ](t − n, s − m) , t, s ∈ R , n, m ∈ Z .
Moreover, any f = TΦ F ∈ Vϕ2 can be expressed as
D
E
f (t, s) = F (x, y), Kt,s (x, y) 2

L (0,1)2

,

t, s ∈ R ,

(8)

P
where Kt,s (x, y) = n,m∈Z Φ(t − n, s − m) e−2πinx e−2πimy ∈ L2 (0, 1)2 . The details can
be found in [11]. Thus, for fixed a, b ∈ [0, 1) and n, m ∈ Z we get the expression for the
samples
f (a + n, b + m) = F, Ka+n,b+m L2 (0,1)2 = F, e−2πinx e−2πimy Ka,b (x, y) L2 (0,1)2 .

The sequence e−2πinx e−2πimy Ka,b (x, y) n,m∈Z is a Riesz basis for L2 (0, 1)2 if and only
if the inequalities
0 < kKa,b k0 ≤ kKa,b k∞ < ∞ hold. Besides, its dual Riesz basis is
 −2πinx
the sequence e
e−2πimy /Ka,b (x, y) n,m∈Z .
Proceeding as in Section 2.2 we derive for any f ∈ VΦ2 the sampling formula
X
f (t, s) =
f (a + n, b + m) Sa,b (t − n, s − m) , t, s ∈ R ,
(9)
n,m∈Z

15


where Sa,b := TΦ 1/Ka,b ∈ VΦ2 . The convergence of the series is absolute due to the
unconditional convergence of a Riesz basis expansion and it is also uniform on R2 since
the shift-invariant space VΦ2 is a RKHS contained in L2 (R2 ).
The case of two separable variables
Consider two generators ϕ(x) and ψ(y) in L2 (R) satisfying the same hypotheses as in
Section 2, and define Φ(x, y) := ϕ(x)ψ(y) in L2 (R2 ) and its associated shift-invariant
2 := V 2 . It can be seen as the tensor product V 2 ⊗ V 2 . For a brief on
subspace Vϕψ
ϕ
Φ
ψ
tensor products see, for instance, [14] and references therein.
Straightforward, Tϕψ := TΦ = Tϕ ⊗ Tψ holds and, consequently, for H, G ∈ L2 (0, 1)
and n, m ∈ Z, the shifting property reads:


Tϕψ e−2πinx e−2πimy H(x)G(y) (t, s) = Tϕ [H](t − n)Tψ [G](s − m) , t, s ∈ R .
2 , where F ∈ L2 (0, 1)2 , can be expressed as
Besides, any f = Tϕψ F ∈ Vϕψ
D
E
e s (y)
f (t, s) = F (x, y), Kt (x)K
, t, s ∈ R ,
2
2

(10)

L (0,1)

P
−2πimy belong to
e s (y) = P
where Kt (x) = n∈Z ϕ(t − n) e−2πinx and K
m∈Z ψ(s − m) e
2
L (0, 1). Thus, for fixed a, b ∈ [0, 1) and n, m ∈ Z we get the following expression for
the samples
e b+m
f (a + n, b + m) = F, Ka+n K

L2 (0,1)2

e b (y)
= F, e−2πinx Ka (x) e−2πimy K

L2 (0,1)2

.

e b k0 ≤ kK
e b k∞ < ∞, the
Assumingthat 0 < kKa k0 ≤ kKa k∞ < ∞ and 0 < kK
−2πinx
−2πimy
2
e
sequence e
e
Ka (x)Kb (y) n,m∈Z is a Riesz basis for L (0, 1)2 with dual
 −2πinx −2πimy
e b (y))
basis e
e
/(Ka (x)K
. By using the notation in Section 3 the
n,m∈Z


sequences xn ⊗ x
em n,m∈Z and yn ⊗ yem n,m∈Z form a pair of dual Riesz bases in
2
2
L (0, 1) .
2 the sampling formula
Proceeding as in Section 2.2 we derive for any f ∈ Vϕψ
X
f (t, s) =
f (a + n, b + m) Sa (t − n) Seb (s − m) , t, s ∈ R ,
(11)
n,m∈Z



e b ∈ V 2 . The convergence of the series
where Sa := Tϕ 1/Ka ∈ Vϕ2 and Seb := Tψ 1/K
ψ
is absolute due to the unconditional convergence of a Riesz basis expansion and it is
2 is a RKHS contained in L2 (R2 ).
also uniform on R2 since the shift-invariant space Vϕψ

4.2

The resulting sampling formulas

Throughout this section we will use the forward differences defined as:
0
0 −1 
1,0
k−1,k0 
0,1
∆k,k
f = ∆+
∆k,k
f ,
+ f = ∆+ ∆+
+
where
1,0
∆+
f (t, s) = f (t + 1, s) − f (t, s)

and
16

∆0,1
+ f (t, s) = f (t, s + 1) − f (t, s) ;

we adopt the convention ∆0,0
+ = I.
2
For any f ∈ Vϕψ (such that Tϕψ F = f ) or f ∈ VΦ2 (such that TΦ F = f ), we
respectively get
0

0

−2πix
e s (y)
∆k,k
− 1)k (e−2πiy − 1)k Kt (x) K
+ f (t, s) = F, (e

or

0

L2 (0,1)2

0

−2πix
∆k,k
− 1)k (e−2πiy − 1)k Kt,s (x, y)
+ f (t, s) = F, (e

L2 (0,1)2

,

.

(12)

(13)

Next we exhibit a couple of examples involving forward differences up to second
order. The obtained results could be generalized to high order forward differences and
even to higher dimensions; besides, backward differences or averages could be used as
in the previous section.
0

0
2
Using the samples ∆k,k
+ f (a + 2n, b + 3m), (k = 0, 1 ; k = 0, 1, 2) in Vϕψ
2 we have the expression for the samples (see (12)):
For any f = Tϕψ F in Vϕψ

D
E
0
0
−2πix
e b (y) ,
∆k,k
−1)k (e−2πiy −1)k e−2πi(2n)x Ka (x)e−2πi(3m)y K
+ f (a+2n, b+3m) = F, (e
where n, m ∈ Z and k = 0, 1 ;k 0 = 0, 1, 2.
Consider the Riesz basis xn ⊗ x
em n,m∈Z for L2 (0, 1)2 , where xn := e−2πinx Ka (x)
e b (y), n, m ∈ Z. According to Lemma 1, for the partition
and x
em := e−2πimy K
{x2n ⊗ x
e3m }n,m∈Z ∪ {x2n ⊗ x
e3m+1 }n,m∈Z ∪ {x2n ⊗ x
e3m+2 }n,m∈Z
∪ {x2n+1 ⊗ x
e3m }n,m∈Z ∪ {x2n+1 ⊗ x
e3m+1 }n,m∈Z ∪ {x2n+1 ⊗ x
e3m+2 }n,m∈Z ,
the matrix of the Riesz basis associated with the above samples will be the Kronecker
product (tensor product) of the corresponding matrices derived, for the one-dimensional
case, in Section 3. That is:


1
0
0
0
0 0

 −1 1
0
0
0 0




1
0 0


1 0
1
−2
1
0
0
0
.
⊗ −1 1 0 = 

−1 1
−1 0
0
1
0 0


1 −2 1
 1 −1 0 −1 1 0
−1

2

−1

Its inverse is the Kronecker product of their respective
[17]); namely:

1 0 0

 1 1 0



1 0 0
1 2 1
1 0


⊗ 1 1 0 =
1 0 0
1 1

1 2 1
1 1 0
1 2 1

1

−2 1

inverses (see, for instance, Ref.
0
0
0
1
1
1

0
0
0
0
1
2


0
0

0
.
0

0
1

2 reads:
Thus, the corresponding sampling formula for each f ∈ Vϕψ

17

(14)

f (t, s) =

1 X
2
X nX
n,m∈Z

o
0
ek,k0
∆k,k
+ f (a + 2n, b + 3m) Ta,b (t − 2n, s − 3m) ,

t, s ∈ R ,

k=0 k0 =0

k,k0
where the sampling functions Tea,b
(t, s) are obtained from the columns of the inverse

matrix (14) as
0,0
Tea,b
(t, s) = Sa (t) Seb (s) + Sa (t) Seb (s − 1) + Sa (t) Seb (s − 2) + Sa (t − 1) Seb (s)

+ Sa (t − 1) Seb (s − 1) + Sa (t − 1) Seb (s − 2) ;
Te0,1 (t, s) = Sa (t) Seb (s − 1) + 2Sa (t) Seb (s − 2) + Sa (t − 1) Seb (s − 1) + 2Sa (t − 1) Seb (s − 2) ;
a,b

0,2
Tea,b
(t, s) = Sa (t) Seb (s − 2) + Sa (t − 1) Seb (s − 2) ;
1,0
Tea,b
(t, s) = Sa (t − 1) Seb (s) + Sa (t − 1) Seb (s − 1) + Sa (t − 1) Seb (s − 2) ;
1,1
Tea,b
(t, s) = Sa (t − 1) Seb (s − 1) + 2Sa (t − 1) Seb (s − 2) ;
1,2
Tea,b
(t, s) = Sa (t − 1) Seb (s − 2) ,

t, s ∈ R .

Observe that Tϕψ y2n+k ⊗ ye3m+k0 = Sa (t − (2n + k))Seb (s − (3m + k 0 )), where
e b (y), n, m ∈ Z.
yn := e−2πinx /Ka (x) and yem := e−2πimy /K
0

0
2
Using the samples ∆k,k
+ f (a + 3n, b + 3m), (k, k = 0, 1, 2) in VΦ

For any f = TΦ F in VΦ2 we have the expression for the samples (see (13)):
D
E
0
−2πix
k −2πiy
k0 −2πi(3n)x −2πi(3m)y
∆k,k
f
(a
+
3n,
b
+
3m)
=
F,
(e
−
1)
(e
−
1)
e
e
K
(x,
y)
,
a,b
+

where n, m ∈ Z and k, k 0 = 0, 1, 2. Consider the Riesz basis xn,m n,m∈Z for L2 (0, 1)2 ,
where xn,m := e−2πinx e−2πimy Ka,b (x, y), n, m ∈ Z. According to Lemma 1, for the
partition
{x3n,3m }n,m∈Z ∪ {x3n,3m+1 }n,m∈Z ∪ {x3n,3m+2 }n,m∈Z
∪ {x3n+1,3m }n,m∈Z ∪ {x3n+1,3m+1 }n,m∈Z ∪ {x3n+1,3m+2 }n,m∈Z
∪ {x3n+2,3m }n,m∈Z ∪ {x3n+2,3m+1 }n,m∈Z ∪ {x3n+2,3m+2 }n,m∈Z ,
the Riesz basis associated to the above samples has as matrix the Kronecker product:


1
0
0
0
0
0
0
0 0
−1 1
0
0
0
0
0
0 0


 1 −2 1

0
0
0
0
0
0


 
 

1
0 0
1
0 0
0
1
0
0
0
0 0
−1 0

−1 1 0 ⊗ −1 1 0 =  1 −1 0 −1 1
0
0
0 0

,

0
0 0
1 −2 1
1 −2 1
−1 2 −1 1 −2 1

1
0
0 −2 0
0
1
0 0


−1 1
0
2 −2 0 −1 1 0
1 −2 1 −2 4 −2 1 −2 1
18

whose inverse is:

1
1



 
 1
1
1 0 0
1 0 0

1 1 0 ⊗ 1 1 0 = 1

1
1 2 1
1 2 1

1

1
1

0
1
2
0
1
2
0
1
2

0
0
1
0
0
1
0
0
1

0
0
0
1
1
1
2
2
2

0
0
0
0
1
2
0
2
4

0
0
0
0
0
1
0
0
2

0
0
0
0
0
0
1
1
1

0
0
0
0
0
0
0
1
2


0
0

0

0

0
.
0

0

0

(15)

1

Thus, the corresponding sampling formula for each f ∈ VΦ2 reads:

f (t, s) =

2
X n X
n,m∈Z

o
0
k,k0
∆k,k
+ f (a + 3n, b + 3m) Ta,b (t − 3n, s − 3m) ,

t, s ∈ R ,

k,k0 =0

0

k,k
(t, s) are obtained from the columns of the inverse
where the sampling functions Ta,b

matrix (15) as:
0,0
(t, s) = Sa,b (t, s) + Sa,b (t, s − 1) + Sa,b (t, s − 2) + Sa,b (t − 1, s) + Sa,b (t − 1, s − 1)
Ta,b

+ Sa,b (t − 1, s − 2) + Sa,b (t − 2, s) + Sa,b (t − 2, s − 1) + Sa,b (t − 2, s − 2) ;
0,1
(t, s)
Ta,b

= Sa,b (t, s − 1) + 2Sa,b (t, s − 2) + Sa,b (t − 1, s − 1) + 2Sa,b (t − 1, s − 2)
+ Sa,b (t − 2, s − 1) + 2Sa,b (t − 2, s − 2) ;

0,2
Ta,b
(t, s) = Sa,b (t, s − 2) + Sa,b (t − 1, s − 2) + Sa,b (t − 2, s − 2) ;
1,0
Ta,b
(t, s) = Sa,b (t − 1, s) + Sa,b (t − 1, s − 1) + Sa,b (t − 1, s − 2) + 2Sa,b (t − 2, s − 2)

+ 2Sa,b (t − 2, s − 1) + 2Sa,b (t − 2, s − 2) ;
1,1
Ta,b
(t, s) = Sa,b (t − 1, s − 1) + 2Sa,b (t − 1, s − 2) + 2Sa,b (t − 2, s − 1) + 4Sa,b (t − 2, s − 2) ;
1,2
Ta,b
(t, s) = Sa,b (t − 1, s − 2) + 2Sa,b (t − 2, s − 2) ;
2,0
Ta,b
(t, s) = Sa,b (t − 2, s) + Sa,b (t − 2, s − 1) + Sa,b (t − 2, s − 2) ;
2,1
Ta,b
(t, s) = Sa,b (t − 2, s − 1) + 2Sa,b (t − 2, s − 2) ;
2,2
Ta,b
(t, s) = Sa,b (t − 2, s − 2) ,

t, s ∈ R .

Observe that TΦ y3n+k,3m+k0 = Sa,b (t − (3n + k), s − (3m + k 0 )), where yn,m :=
e−2πinx e−2πimy /Ka,b (x, y), n, m ∈ Z.

5

Conclusions

Starting from a Shannon-type sampling formula valid in a shift-invariant subspace,
we have introduced a systematic way to derive stable sampling formulas involving
19

differences and/or averages from the data available in the former sampling formula. The
mathematical technique is very simple and relies on a representation of the elements
of the shift-invariant subspace as an inner product in the auxiliary space L2 (0, 1); this
allows us to identify an involved dual pair of Riesz bases by means of an easy matrices
calculation. Besides, a suitable isomorphism between L2 (0, 1) and the shift-invariant
subspace gives the resulting sampling formulas. These formulas can be viewed as a
multichannel filter processing.
Thus, a non-exhaustive variety of examples in the one-dimensional case are exhibited
throughout the work. All these sampling formulas can be understood as a change of
basis in the original Shannon-type sampling formula. The two-dimensional case is also
treated: it reduces to the one-dimensional case via the Kronecker product of matrices.
The same procedure could be generalized to higher dimensions.
Although the obtained formulas are Riesz bases expansions, as it was pointed out
the used technique still remains valid for overcomplete expansions i.e., the frame setting. Furthermore, this technique also works in other types of sampling formulas:
for instance, those involving samples of filtered versions of the functions in the shiftinvariant subspace.
Acknowledgments: This work has been supported by the grant MTM2014-54692-P
from the Spanish Ministerio de Economı́a y Competitividad (MINECO)

References
[1] A. Aldroubi and K. Gröchenig. Non-uniform sampling and reconstruction in shiftinvariant spaces. SIAM Rev., 43:585–620, 2001.
[2] A. Aldroubi and M. Unser. Sampling procedures in function spaces and asymptotic
equivalence with Shannon sampling theorem. Numer. Funct. Anal. Optimiz., 15:1–
21, 1994.
[3] A. Aldroubi, Q. Sun, and W-S. Tang. Convolution, average sampling, and a
Calderon resolution of the identity for shift-invariant spaces. J. Fourier Anal.
Appl., 11(2):215–244, 2005.
[4] A. Aldroubi, M. Unser and M. Eden. Cardinal spline filters: Stability and convergence to the ideal sinc interpolator. Signal Process., 28:127–138, 1992.
[5] C. Boor, R. DeVore, and A. Ron. Approximation from shift-invariant subspaces
in L2 (Rd ). Trans. Amer. Math. Soc., 341:787–806, 1994.
[6] C. Boor, R. DeVore, and A. Ron. The structure of finitely generated shift-invariant
spaces in L2 (Rd ). J. Funct. Anal., 119:37–78, 1994.
[7] W. Chen, S. Itoh, and J. Shiki. On sampling in shift invariant spaces. IEEE Trans.
Signal Processing, 48(10):2802–2810, 2002.
[8] O. Christensen. An Introduction to Frames and Riesz Bases. Birkhäuser, Boston,
2003.
20

[9] A. G. Garcı́a, G. Pérez-Villalón, and A. Portal. Riesz bases in L2 (0, 1) related to
sampling in shift-invariant spaces. J. Math. Anal. Appl., 308(2):703–713, 2005.
[10] A. G. Garcı́a and G. Pérez-Villalón. Dual frames in L2 (0, 1) connected with
generalized sampling in shift-invariant spaces. Appl. Comput. Harmon. Anal.,
20(3):422–433, 2006.
[11] A. G. Garcı́a and G. Pérez-Villalón. Multivariate generalized sampling in shiftinvariant spaces and its approximation properties. J. Math. Anal. Appl., 355:397–
413, 2009.
[12] A. G. Garcı́a, M. A. Hernández-Medina and M. J. Muñoz-Bouzo. The Kramer
sampling theorem revisited. Acta Appl. Math., 133:87–111, 2014.
[13] A. G. Garcı́a. Sampling theory and reproducing kernel Hilbert spaces. In: Operator
Theory, D. Alpay Editor, Springer, Basel, 2015.
[14] A. G. Garcı́a, A. Ibort and M. J. Muñoz-Bouzo. Modeling sampling in tensor
products of unitary invariant subspaces. J. Funct. Spaces, Vol. 2016, Article ID
4573940, 14 pp., 2016.
[15] S. Kang and K. H. Kwon. Generalized average sampling in shift invariant spaces.
J. Math. Anal. Appl., 377:70–78, 2011.
[16] J. M. Kim and K. H. Kwon. Sampling expansion in shift invariant spaces. Int. J.
Wavelets Multiresolut. Inf. Process., 6:223–248, 2008.
[17] A. J. Laub. Matrix Analysis. Siam, Philadelfia, 2005.
[18] J. J. Lei, R. Q. Jia, and E. W. Cheney. Approximation from shift-invariant spaces
by integral operators. SIAM J. Math. Anal., 28(2):481–498, 1997.
[19] I. J. Schoenberg. Cardinal Spline Interpolation. Siam, Philadelfia, 1993.
[20] W. Sun and X. Zhou. Average sampling in spline subspaces. Appl. Math. Lett.,
15:223–237, 2002.
[21] M. Unser. Sampling 50 Years After Shannon. Proc. IEEE, 88(4):569–587, 2000.
[22] G. G. Walter. A sampling theorem for wavelet subspaces. IEEE Trans. Inform.
Theory, 38:881–884, 1992.
[23] X. Zhou and W. Sun. On the Sampling Theorem for Wavelet Subspaces. J. Fourier
Anal. Appl., 5(4):347–354, 1999.

21

