The frame sets of the B-splines

Dengfeng LI

Front. Math. China ›› 2023, Vol. 18 ›› Issue (4) : 223 -233.

PDF (499KB)
Front. Math. China ›› 2023, Vol. 18 ›› Issue (4) : 223 -233. DOI: 10.3868/s140-DDD-023-0021-x
SURVEY ARTICLE
SURVEY ARTICLE

The frame sets of the B-splines

Author information +
History +
PDF (499KB)

Abstract

It is not completely clear which elements constitute the frame sets of the B-splines currently, but some considerable results have been obtained. In this paper, firstly, the background of frame set is introduced. Secondly, the main progress of the frame sets of the B-splines in the past more than twenty years are reviewed, and particularly the progress for the frame set of the 2 order B-spline and the frame set of the 3 order B-spline are explained, respectively.

Keywords

B-spline / Gabor frame / frame set

Cite this article

Download citation ▾
Dengfeng LI. The frame sets of the B-splines. Front. Math. China, 2023, 18(4): 223-233 DOI:10.3868/s140-DDD-023-0021-x

登录浏览全文

4963

注册一个新账户 忘记密码

1 Introduction

Suppose that gL2(R) and (a,b)R+2, where R+2={(x,y)|x>0,y>0}. Let

(g,a,b)={ei2πmbxg(xna)}m,nZ={EmbTnag(x)}m,nZ.

Then (g,a,b) is called a Gabor system for L2(R), where Z is the integer set, Ta denotes shift operator:

Ta:L2(R)L2(R),(Taf)(x)=f(xa),aR,

and Eb is the modulation operator:

Eb:L2(R)L2(R),(Ebf)(x)=ei2πbxf(x),bR.

Defintion 1.1 If there exists constants A>0 and B< such that

fL2(R),Af||2m,nZ|f,EmbTnag|2Bf2,

then (g,a,b) is called a Gabor frame for L2(R) with bounds A and B, and (g,a,b) is also called a Weyl-Heisenberg frame, correspondingly, where g is called the window function or the generating function. If only the right inequality holds in (1.1), then (g,a,b) is called a Bessel sequence for L2(R); if A=B in (1.1), then the frame is called a tight Gabor frame for L2(R).

The existence, construction, and properties of the generating function g and time-frequency parameters (a,b) in Gabor frame (g,a,b) constitute the part of main contents of Gabor analysis [9, 20-22, 26, 37, 41]. The research on Gabor analysis predates wavelet analysis historically, but the research on Gabor frame has gained great attention with the rapid development of wavelet analysis [15, 16, 20-22], this is mainly due to the widespread applications of Gabor frame in fields such as image processing, digital communication and others [20-22, 30, 45, 46]. As a tool, Gabor frame can be used to study quasi-differential operator, complex analysis, non-commutative geometry, etc. [23, 38, 39, 43]. Specifically, Gabor frame is a fundamental research topic in time-frequency analysis [18, 20-22].

Obviously, whether a Gabor system forms a frame depends on the window function and time-frequency parameters. Ideally, the characteristic description of the window function g and time-frequency parameters (a,b) should be established so that the Gabor system (g,a,b) is a frame for L2(R), but this is a very challenging problem, which has not been completely solved yet. However, for given gL2(R), the search for time-frequency parameters (a,b) in the system (g,a,b) forming a frame for L2(R) has obtained some results, and the study of this situation is currently one of the basic research contents of Gabor frame. Therefore, the concept of frame set was born.

Definition 1.2 Assume that gL2(R), the following set

F(g)={(a,b)R+2|(g,a,b)isaGaborframeforL2(R)}

is called the frame set associated with g.

It is easy to know that if (a,b)F(g), then ab1, which comes from the density theorem of Gabor frame as follows (see [9, 15-17, 29, 37]).

Theorem 1.1  Let gL2(R) and (a,b)R+2.

(1) If ab>1, then (g,a,b) is not a frame for L2(R).

(2) The system (g,a,b) is an orthonormal basis for L2(R) if and only if the (g,a,b) is a tight frame for L2(R), g=1, and ab=1.

In general, for fixed gL2(R), fully giving its frame set is an unresolved problem. But an important result of the frame set was established in [19].

Theorem 1.2 [19]  If gL2(R) and gM1(R), then F(g) is a non-empty open set in R+2, where the modulation space M1(R) (see [7]) is defined by

M1(R)={f|RR|WFff(x,ξ)|dxdξ<},

and WFg is the windowed Fourier transform, i.e.,

WFgf(x,ξ)=Rf(t)g(tx)¯ei2πtξdt.

For any gL2(R), it is not known which elements F(g) are composed of, but for the special function g, F(g) has been fully characterized.

Theorem 1.3  If g is any of the following functions, then

F(g)={(a,b)R+2|ab<1}.

(1) The Guassian function g(x)=eπx2 (see [39, 42, 44]).

(2) The hyperbolic secant function g(x)=1coshx (see [34]).

(3) The one-sided exponential decay function g(x)=exχ[0,)(x) (see [31]).

(4) The two-sided exponential decay function g(x)=e|x| (see [32]).

(5) The total positive function of finite type g(x), i.e., the Fourier transform of g(x) has the form g^(ξ)=j=1K(1+i2πβjξ)1, where KZ+, Z+ is the positive integer set, and jk,βjβkR (see [27]).

(6) The total positive function of Gaussian type g(x), i.e., the Fourier transform of g(x) has the form g^(ξ)=j=1K(1+i2πβjξ)1eaξ2, where KZ+ and jk,βjβkR,a>0 (see [25]).

Remark 1.1 It is easy to see that a two-sided exponential decay function is a total positive function of finite type and the Gaussian function is a total positive function of Gaussian type.

In addition, the structure of F(g) of g(x)=χ[0,c](x)(c>0) is quite complex, and the F(g) has been fully described in [13]. In history, earlier relevant results of the F(g) were established in [28, 33].

The structure of F(g) of Haar function g(x)=χ[12,0)(x)+χ[0,12)(x) was completely given in [14]. We have recently seen that the profound results on the frame sets of Hergoltz functions and the rational functions were established in [6].

It is well known that the spline functions are a type of piecewise smooth functions with certain smoothness at the connections of each segment. It is a powerful tool for function approximation and a commonly used method for interpolation and fitting of scattered data. It is not only widely used in data processing, differential equations, numerical calculus and other mathematical branches, but also closely related to functional analysis, computational geometry, variational problems, optimal control, statistics and other directions. At the same time, the spline functions have profound applications in many engineering sciences such as aerospace, computers, and materials. Even the simplest B-spline in such a class of functions, the specific composition of its frame set is not completely clear. In the past twenty years, there have been many results in the study of the B-spline frame sets, but there are still many problems that need to be addressed. The main advances in this area of research will be reviewed in this paper.

The paper is organized as follows: Section 2 reviews the existing results of the B-spline frame set so far. In Sections 3 and 4, the results of the 2 order B-spline and the 3 order B-spline frame sets are described in detail, respectively.

2 Results on the frame sets of the B-splines

The definition of the B-splines is as follows.

g1(x)=χ[12,12](x),gN(x)=(g1gN1)(x),N2.

For the discontinuous function g1, the structure of F(g1) is very complex, as shown in [13], one of the reason for this situation is that g1 is the only B-spline that does not belong to M1(R). In fact, F(g1) is not an open set in R+2, and (a,b)(0,1]×(0,1]F(g1). It was shown in [24] that the window function other than M1(R) has poor properties, so its frame set is difficult to determine. Notice that the non-smooth spline does not belong to M1(R), hence determining the frame set of the non-smooth spline is more difficult.

For N2, gNM1(R), thus by Theorem 1.2, F(gN) is an open set in R+2. In addition, the paper [8] lists the characterization of F(gN) as one of the six open questions in frame theory.

In the following, we assume that N2. In [10-12, 35], the necessary and sufficient conditions for the duality of Gabor frames and the properties of Zak transforms are used to discuss in depth the frame sets of different classes of window functions with compact support and certain smoothness, and gN is a special case of these different types of window functions, so the characterizations of the following F(gN) comes from the literatures [10-12, 35]. However, let's first review the number pairs that are not in F(gN). Note that supp gN=[N2,N2], hence for a>N, the support union of TnagN(nZ) cannot cover the real number set, thus b>0, (g,a,b) cannot be a frame for L2(R).

Theorem 2.1  For a>N and b>0, (a,b)F(gN).

In addition, using the necessary conditions of Gabor frame (see [37], Theorem 8.2.2), it is easy to prove

Theorem 2.2  If a>0,b=2,3,, then (a,b)F(gN).

For the critical values a and b that satisfy ab=1, the following was obtained by the Zak transform properties of gN and the Zibulski Zeevi matrix in [40].

Theorem 2.3  If ab=1, then (a,b)F(gN).

Combining Theorem 2.1, Theorem 2.2 and Theorem 2.3, we have

F(gN){(a,b)R+2|ab<1,a<N,b2,3,}.

Based on (2.1), Gröchenig [24] proposed the following conjecture.

Conjecture 2.1 (Gröchenig)  N2,

F(gN)={(a,b)R+2|ab<1,a<N,b2,3,}.

Subsequently, it was proved that the following result holds by fully using the Zak transform properties of g2 and the Zibulski Zeevi matrix in [36, 40].

Theorem 2.4  Let ab=56,b[73,83], then (a,b)F(g2).

Theorem 2.4 negates Gröchenig conjecture, one specific example is that if a=13 and b=52, then (a,b)F(g2). In fact, the other values of a and b were provided in [36] such that (a,b)F(gN), these values of a and b were included in Theorem 2.4 for N=2. At the same, paper [36] proposed the following conjecture based on Theorem 2.4.

Conjecture 2.2 If a=12n+1,b=2l+12,n,lZ+,l>n and ab<1, then (a,b)F(g2). Furthermore, if a~b~=2l+12(2n+1),b~[baln2,b+aln2], then (a~,b~)F(g2).

It is obvious to see that when n=1,l=2 in conjecture 2.2, a=13,b=52, a~b~=56,b~[73,83], which is just Theorem 2.4.

In Conjecture 2.2, if a is the reciprocal of a positive even number, then it is believed that this conjecture is also correct by the analysis of number values in [40], that is,

Conjecture 2.3 If a=12n,b=2l+12,n,lZ+,l>n,ab<1, then (a,b)F(g2).

To our knowledge, the above two conjectures has been unsolved yet.

Let's discuss the number pairs in the set F(gN) below. Based on the results of the literatures [10-12, 35], it can be concluded that

(0,N)×(0,4N+3a]F(gN),

[N2,N)×(0,1a)F(gN).

Indeed, formulas (2.2) and (2.3) come from literatures [12] and [11], respectively. In addition, the following results were provided in [11].

Theorem 2.5  LetN2,a>0,b>0 and ab1.

(1) If a<N and bN1, then (a,b)F(gN).

(2) If there exists a positive integer l such that N1bN12 and 21N<al<b1, then (a,b)F(gN).

(3) If b{1,21,,(N1)1}, then (a,b)F(gN).

(4) If a=p1l, where l=1,2,,N1 and p is a positive integer, and b<l1, then (a,b)F(gN).

Remark 2.1 It is easy to see that formula (2.2) implies (1) in Theorem 2.5, and in the case l=1, (2) in Theorem 2.5 is just (2.3). Additionally, when l=1, the results (3) and (4) in Theorem 2.5 all stem from [35].

We also noticed that paper [5] used the partitioning technique to obtain

[N3,N2]×[4N+3a,2N]F(gN).

In fact, the frame sets of a class of functions (referred to as generalized B-splines) containing B-spline functions discussed in [12] were considered in [5], and the obtained result is to give a new pair of numbers (a,b)F(gN) on the basis of the original results. Combining the formulas (2.2), (2.3), and (2.4) immediately yields

ENF(gN),

where

EN={(a,b)R+2|ab<1,0<a<N,0<bmax{4N+3a,2N}}.

As far as we know, the set EN in (2.5) is the largest subset of the frame set F(gN) of the B-spline, except for N=2 and N=3. When N=2 and N=3, F(gN) contain more values of a and b than in the set EN, so the following two sections summarize the results of the 2 order B-spline and the 3 order B-spline frame sets, respectively.

3 The case of the 2 order B-spline

The definition of the 2 order B-spline is as follows.

g2(x)={1+x,x[1,0],1x,x[0,1].

Of course, F(g2) is an open set in R+2, and which elements exactly constitute F(g2) is still an open question. The proof ideas of existing results all rely on the necessary and sufficient conditions for the dual Gabor frame established by Janssen (see [9, 20, 37]).

Theorem 3.1  Let g,hL2(R) and (a,b)R+2. Then two Bessel sequences (g,a,b) and (h,a,b) are dual Gabor frames if and only if

1blZg(xlakb)h(xla)¯=δk,0,a.e.,x[0,a].

For the Bessel sequences (g,a,b) and (h,a,b), if

fL2(R),f(x)=m,nZf,EmbTnagEmbTnah(x),

then (g,a,b) and (h,a,b) are called the dual Gabor frames, alternatively, h is called a dual of g. In fact, the two are mutually dual.

In order to reflect the historical process of obtaining the elements in F(g2) based on Theorem 3.1, we summarize as follows.

Theorem 3.2  For g2, the following statements hold.

(1) If 0<a<2 and 0<b22+a, then (a,b)F(g2), and there exists a unique dual hL2(R)L(R) such that supp h[a2,a2] (see [12]).

(2) If 0<a<2 and 22+a<b42+3a, then (a,b)F(g2), and there is a unique dual hL2(R)L(R) such that supp h[3a2,3a2] (see [36]).

(3) If 0<a<12 and 42+3a<b21+a, then (a,b)F(g2), and there is a unique dual hL2(R)L(R) such that supp h[5a2,5a2] (see [3]).

(4) If 12a45, and 42+3a<b62+5a and b>1, then (a,b)F(g2), and there exists a unique dual hL2(R)L(R) such that supp h[5a2,5a2] (see [3]).

(5) If 23a1 and 42+3a<b<1, then (a,b)F(g2), and there is a unique compactly supported dual hL2(R)L(R) (see [4]).

Certainly, the formula (2.5) also holds for N=2, that is,

E2F(g2).

Recently, some new elements of the set F(g2) were found in [2] based on the linear algebra method used in [4].

Theorem 3.3  Γ3ΛF(g2), where

Γ3={(a,b)R+2|a(0,29](27,12),b(42+3a,21+a]},

and

Λ=m=3Λm.

Here

Λ3={(a,b)R+2|a[12,45],b(42+3a,62+5a),b>1}.

And for m4,

Λm={(a,b)R+2|a[m3m2,2(m1)2m1],b(2(m1)2+(2m3)a,γm,a],b>1},

γm,a=min{2m2+(2m1)a,21+a}.

The proof idea of Theorem 3.3 is still based on the necessary and sufficient condition of the dual Gabor frames, but this condition is a concretization of the special cases of g and h in Theorem 3.1, that is,

Theorem 3.4  Let N2,(a,b)R+2,a<N,m1, and h be a bounded real function with supp h[2m12a,2m12a]. Then Gabor systems (gN,a,b) and (h,a,b) are dual Gabor frames if and only if

1bl=1mm1gN(xlakb)h(xla)=δk,0,|k|m1,a.e.,x[a2,a2].

Remark 3.1 For the proof of Theorem 3.4, the cases m=1 and m=2 were given in [36] and [12], respectively, and the cases m3 were established in [4].

After carefully analyzing the proof methods and processes of Theorems 3.2 and 3.3 and formula 3.1, Okoudjou [41] proposed the following unsolved conjecture.

Conjecture 3.1 [41]

{(a,b)R+2|12a<1,62+5ab<21+a,b>1}F(g2).

4 The case of the 3 order B-spline

The 3 order B-spline is defined by

g3(x)={12x2+32x+98,x[32,12],x2+34,x[12,12],12x232x+98,x[12,32].

Naturely, F(g3) is an open set in R+2, and by (2.5), the following set E3 is a subset of F(g3),

E3={(a,b)R+2|ab<1,0<a<3,0<bmaxa{43+3a,23}}.

Some new elements of F(g3) were founded by making full use of Theorem 3.4 and the invertibility of block matrix in [1], see the following theorem.

Theorem 4.1  Let

Δ3={(a,b)R+2|a[34,65],b(43+3a,63+5a),b>1}.

For m4,

Δm={(a,b)R+2|a[3(m3)(2(m2),3(m1)2m1],b(2(m1)3+(2m3)a,θm,a],b>23},

θm,a=min{2m3+(2m1)a,43+2a}.

If (a,b)Δm, then (a,b)F(g3), and g3 exists a unique dual hL2(R) such that

supph[(2m1)a2,(2m1)a2].

Theorem 4.1 shows that

m=3ΔmF(g3).

References

[1]

Atindehou A G D. On the frame set for the 3-spline. Appl Math 2022; 13: 377–400

[2]

Atindehou A G D, Frederick C, Kouagou Y B, Okoudjou K A. On the frame set of the second-order B-spline. Appl Comput Harmon Anal 2023; 62: 237–250

[3]

AtindehouA G DKouagouY BOkoudjouK A. Frame sets for generalized B-splines. 2018, arXiv: 1804-02450

[4]

Atindehou A G D, Kouagou Y B, Okoudjou K A. Frame sets for a class of compactly supported continuous functions. Asian-Eur J Math 2020; 13(5): 2050093

[5]

AtindehouA G DYebeniBKouagouY BOkoudjouK. A. Frame sets for generalized B-splines. 2018, arXiv: 1804.02450v2

[6]

Belov Y, Kulikov A, Lyubarskii Y. Gabor frames for rational functions. Invent Math 2023; 23(2): 431–466

[7]

BényiÁOkoudjouK A. Modulation Spaces: With Applications to Pseudodifferential Operators and Nonlinear Schrödinger Equations. Applied and Numerical Harmonic Analysis. New York: Birkhaüser/Springer, 2020

[8]

ChristensenO. Six (seven) problems in frame theory. In: New Perspectives on Approximation and Sampling Theory. Applied and Numerical Harmonic Analysis. Cham: Birkhaüser/Springer, 2014, 337–358

[9]

ChristensenO. An Introduction to Frames and Riesz Bases, 2nd ed. Applied and Numerical Harmonic Analysis. Cham: Birkhaüser/Springer, 2016

[10]

Christensen O, Kim H O, Kim R Y. Gabor windows supported on [−1,1] and compactly supported dual windows. Appl Comput Harmon Anal 2010; 28(1): 89–103

[11]

Christensen O, Kim H O, Kim R Y. On Gabor frames generated by sign-changing windows and B-splines. Appl Comput Harmon Anal 2015; 39(3): 534–544

[12]

Christensen O, Kim H O, Kim R Y. On Gabor frame set for compactly supported continuous functions. J Inequal Appl 2016; 94

[13]

Dai X R, Sun Q Y. The abc-problem for Gabor systems. Mem Amer Math Soc 2016; 244(1152): ix+99 pp

[14]

DaiX RZhuM. Frame sets for Gabor systems with Haar window. 2022, arXiv: 2205.06479v1

[15]

Daubechies I. The wavelet transform, time-frequency localization and signal analysis. IEEE Trans Inform Theory 1990; 36(5): 961–1005

[16]

DaubechiesI. Ten Lectures on Wavelets. CBMS-NSF Regional Conference Series in Applied Mathematics, Vol 61. Philadelphia, PA: SIAM, 1992

[17]

Daubechies I, Landau H, Landau Z. Gabor time-frequency lattices and the Wexler-Raz identity. J Fourier Anal Appl 1995; 1(4): 437–478

[18]

Feichtinger H G, Gröchenig K. Gabor frames and time-frequency analysis of distributions. J Funct Anal 1997; 146(2): 464–495

[19]

Feichtinger H G, Kaiblinger N. Varying the time-frequency lattices of Gabor frames. Trans Amer Math Soc 2004; 356(5): 2001–2023

[20]

FeichtibgerH GStrohmerT (eds). Gabor Analysis and Algorithms: Theory and Applications. Applied and Numerical Harmonic Analysis. Boston, MA: Birkhaüser Boston, Inc, 1998

[21]

FeichtingerH GStrohmerT (eds). Advances in Gabor Analysis. Applied and Nu- merical Harmonic Analysis. Boston, MA: Birkhaüser Boston, Inc, 2003

[22]

GröchenigK. Foundations of Time-Frequency Analysis. Applied and Numerical Har- monic Analysis. Boston, MA: Birkhaüser Boston, Inc, 2001

[23]

Gröchenig K. Time-frequency analysis of Sjöstrand’s class. Rev Mat Iberoam 2006; 22(2): 703–724

[24]

Gröchenig K. The mystery of Gabor frames. J Fourier Anal Appl 2014; 20(4): 865–895

[25]

Gröchenig K, Romero J L, Stöckler J. Sampling theorems for shift-invariant spaces, Gabor frames, and totally positive functions. Invent Math 2018; 211(3): 1119–1148

[26]

Gröchenig K, Romero J L, Stöckler J. Sharp results on sampling with derivatives in shift-invariant spaces and multi-window Gabor frames. Constr Approx 2020; 51(1): 1–25

[27]

Gröchenig K, Stöckler J. Gabor frames and totally positive functions. Duke Math J 2013; 162(6): 1003–1031

[28]

Gu Q, Han D G. When a characteristic function generates a Gabor frame?. Appl Comput Harmon Anal 2008; 24(3): 290–309

[29]

Heil C. History and evolution of the density theorem for Gabor frames. J Fourier Anal Appl 2007; 13(2): 113–166

[30]

HlawatschFMatzG. Wireless communications over rapidly time-varying channels. Amsterdam: Academic Press, 2011

[31]

Janssen A J E M. Some Weyl-Heisenberg frame bound calculations. Indag Math (N S) 1996; 7(2): 165–183

[32]

Janssen A J E M. On generating tight Gabor frames at critical density. J Fourier Anal Appl 2003; 9(2): 175–214

[33]

JanssenA J E M. Zak transforms with few zeros and the tie. In: Advances in Gabor Analysis. App Numer Harmon Anal. Boston, MA: Birkhaüser Boston, Inc, 2003, 31–70

[34]

Janssen A J E M, Strohmer T. Hyperbolic secants yield Gabor frames. Appl Comput Harmon Anal 2002; 12(2): 259–267

[35]

Kloos T, Stöckler J. Zak transforms and Gabor frames of totally positive functions and exponential B-splines. J Approx Theory 2014; 184: 209–237

[36]

Lemvig J, Nielsen K H. Counterexamples to the B-spline conjecture for Gabor frames. J Fourier Anal Appl 2016; 22(6): 1440–1451

[37]

LiD F. Mathematical Theory of Wavelet Analysis. Beijing: Science Press, 2017 (in Chinese)

[38]

Luef F. Projective modules over noncommutative tori are multi-window Gabor frames for modulation spaces. J Funct Anal 2009; 257(6): 1921–1946

[39]

LyubarskiiY I. Frames in the Bargmann space of entire function functions. In: Entire and Subharmonic Functions. Adv Sovit Math, Vol 11. Providence, RI: AMS, 1992, 167–180

[40]

NielsenK H. The frame set of Gabor systems with B-spline generators. Master Thesis. København: Technical University of Denmark, 2015

[41]

Okoudjou K A. An invitation to Gabor analysis. Notices Amer Math Soc 2019; 66(6): 808–819

[42]

Seip K. Density theorems for sampling and interpolation in the Bargmann-Fock space, I. J Reine Angew Math 1992; 429: 91–106

[43]

SeipK. Interpolation and sampling in spaces of analytic functions. University Lecture Series, Vol 33. Providence, RI: AMS, 2004

[44]

Seip K, Wallstén R. Density theorems for sampling and interpolation in the Bargmann-Fock space, II. J Reine Angew Math 1992; 429: 107–113

[45]

Strohmer T. Approximation of dual Gabor frames, window decay, and wireless communications. Appl Comput Harmon Anal 2001; 11(2): 243–262

[46]

Strohmer T. Pseudodifferential operators and Banach algebras in mobile communications. Appl Comput Harmon Anal 2006; 20(2): 237–249

RIGHTS & PERMISSIONS

Higher Education Press 2023

AI Summary AI Mindmap
PDF (499KB)

743

Accesses

0

Citation

Detail

Sections
Recommended

AI思维导图

/