spectral graph theory examines relationships between the structure of a graph and the eigenvalues ( or spectrum ) of a matrix associated with that graph .
different matrices are able to give different information , but all the common matrices have limitations .
this is because there are graphs which have the same spectrum for a certain matrix but different structure  such graphs are called cospectral with respect to that matrix .
the following are some of the matrices studied in spectral graph theory :    * the adjacency matrix @xmath1 .
this is defined by @xmath2 when @xmath3 and @xmath4 are adjacent and @xmath5 otherwise .
the spectrum of the adjacency matrix can determine the number of edges and if a graph is bipartite , but it can not determine if a graph is connected ( see @xcite for more information ) .
see figure  [ fig:1a ] , the `` saltire pair '' @xcite , for an example of two graphs which are cospectral with respect to the adjacency matrix but not any of the other matrices we will discuss . *
the laplacian @xmath6 .
this matrix is also known as the combinatorial laplacian and is found by taking the difference of the diagonal degree matrix @xmath7 and the adjacency matrix .
the spectrum of the laplacian can determine the number of edges and the number of connected components , but it can not determine if a graph is bipartite ( see @xcite for more information ) .
see figure  [ fig:1b ] @xcite , for an example of two graphs which are cospectral with respect to the laplacian but not any of the other matrices we will discuss .
* the signless laplacian @xmath8 .
this matrix is found by taking the sum of the diagonal degree matrix @xmath7 and the adjacency matrix .
the spectrum of the signless laplacian can determine the number of edges and the number of connected components which are bipartite , but it can not determine if a graph is bipartite or connected ( see @xcite for more information ) .
see figure  [ fig:1c ] @xcite , for an example of two graphs which are cospectral with respect to the signless laplacian but not any of the other matrices we will discuss . * the normalized laplacian @xmath9 .
this matrix is defined by @xmath10 , where by convention if we have an isolated vertex then it will contribute @xmath5 to the spectrum .
the spectrum of the normalized laplacian is closely related to the spectrum of the probability transition matrix of a random walk .
this spectrum can determine if a graph is bipartite and the number of connected components , but it can not determine the number of edges ( see @xcite for more information ) .
see figure  [ fig:1d ] for an example of two graphs which are cospectral with respect to the normalized laplacian but not any of the other matrices we will discuss .
one way to understand what structure the spectrum of a matrix can not identify is to study cospectral graphs .
cospectral graphs for the adjacency matrix @xcite and the laplacian matrix @xcite have been studied , particularly for graphs with few vertices .
cospectral graphs for the signless laplacian have been little studied beyond their enumeration and rules which apply to all of these matrices @xcite .
little is also known about cospectral graphs for the normalized laplacian compared to other matrices .
previously , the only cospectral graphs were bipartite ( complete bipartite graphs @xcite and bipartite graphs found by `` unfolding '' a small bipartite graph in two ways @xcite ) or were regular and cospectral for @xmath1 ( since a regular graph is cospectral for all of these four matrices if it is cospectral for any one of them ) .    in table  [ table ] ,
we have listed the number of graphs with a cospectral mate for the normalized laplacian matrix for graphs on nine or fewer vertices , counted using sage ( see appendix  [ sagecode ] for example code ) .
the number of cospectral graphs for the normalized laplacian have also been counted by wilson and zhu @xcite , though they give percentages and not the count of how many have cospectral mates .
we have also included similar counts for the other three matrices which come from @xcite .
@xmath11    given the relatively small number of graphs with a cospectral mate with respect to the normalized laplacian , it is surprising that so little is known about forming cospectral graphs for that matrix .
the problem is that some of the main tools that are used to form cospectral graphs for other matrices do not generalize to the normalized laplacian .
one such example is a technique known as switching , which is accomplished by replacing edges by non - edges and non - edges by edges between two subsets ( see @xcite ; a simple example is shown in figure  [ fig : switch ] ) .
given some basic assumptions , this is an easy way to construct cospectral graphs for the adjacency matrix .
however , switching does not in general work for the normalized laplacian , in particular , it will only be guaranteed to work when the degrees are unchanged ( see @xcite ) .
so the graphs shown in figure  [ fig : switch ] are _ not _ cospectral with respect to the normalized laplacian .
in section  [ sec : swap ] we will introduce a method to construct cospectral graphs for the normalized laplacian .
this construction will work similarly to switching in that we will make a small local change to the graph by swapping in one bipartite graph with a cospectral mate and show the two graphs still share the same eigenvalues . in section  [ sec : other ]
we will show that if we add additional constraints to the bipartite graphs which are swapped then the resulting graphs are also cospectral with respect to @xmath1 , @xmath12 and @xmath13 . in section  [ sec :
applications ] we will show how to use this construction to produce large families of graphs which are mutually cospectral . finally , in section  [
sec : comments ] we will give some concluding remarks .
the method of finding cospectral graphs for the adjacency matrix by switching reduces to making a local change of the graph from one bipartite subgraph to its complement and showing that the spectrum is unchanged .
we will consider something similar , namely a local change , but instead of replacing a bipartite subgraph with its complement we will swap out a bipartite subgraph for a cospectral mate .    for a subset @xmath14 of the vertices @xmath15 of @xmath16 , we will let @xmath17 $ ]
be the induced subgraph of @xmath16 on the vertex set @xmath14 .
[ thm : biswitch ] let @xmath18 and @xmath19 be bipartite , cospectral graphs with respect to the normalized laplacian on the vertex set @xmath20 such that all edges go between @xmath21 and @xmath22 and where all vertices in @xmath21 have degree @xmath23 for both graphs .
let @xmath24 be a graph on the vertices @xmath25 where @xmath26 $ ] is an arbitrary graph ; @xmath27 $ ] and @xmath28 $ ] have no edges ; there are no edges going between @xmath1 and @xmath21 , between @xmath1 and @xmath22 , or between @xmath29 and @xmath22 ; @xmath30 $ ] is a complete bipartite graph ; and @xmath31=p_1 $ ] .
the graph @xmath32 is defined similarly except that @xmath33=p_2 $ ]
. then @xmath24 and @xmath32 are cospectral with respect to the normalized laplacian .
if the dimension of the eigenspace associated with @xmath34 intersected with the subspace of vectors that are nonzero only on @xmath21 is the same for both @xmath18 and @xmath19 , then the graphs @xmath35 and @xmath36 are also cospectral with respect to the normalized laplacian , where @xmath35 only differs from @xmath24 in that @xmath37 $ ] is the complete graph and similarly @xmath36 only differs from @xmath32 in that @xmath38 $ ] is the complete graph .
an example of the construction described in theorem  [ thm : biswitch ] is shown in figure  [ fig : exswitch ] .    in this case @xmath18 and @xmath19 are the graphs @xmath39 and @xmath40 where we have put the degree @xmath41 vertices into @xmath21 .
it is well known that the spectrum of a complete bipartite graph @xmath42 is @xmath43},1^{[p+q-2]},2^{[1]}$ ] , where the exponent indicates multiplicity , hence these graphs are easily seen to be cospectral .
moreover it is easy to check that the dimension of the eigenspace for the eigenvalue @xmath41 intersected with the subspace of vectors which are nonzero only on @xmath21 is @xmath44 for both graphs , i.e. , we could also have put a complete graph on the vertices of @xmath21 in the graphs shown in figure  [ fig : exswitch ] and still had a cospectral pair .
additional examples of graphs @xmath18 and @xmath19 satisfying ( both ) the conditions on theorem  [ thm : biswitch ] are @xmath45 and @xmath46 where the vertices of degree @xmath47 are all placed in @xmath21 and @xmath48 .
two further examples are shown in figures  [ fig : other1 ] and  [ fig : other2 ] .    to prove theorem  [ thm : biswitch ]
we will find it convenient to work with the harmonic eigenvectors of the normalized laplacian .
namely , if @xmath49 is an eigenvector associated with the eigenvalue @xmath50 , i.e. , @xmath51 , then the corresponding harmonic eigenvector is @xmath52 .
this translates the relationship @xmath51 into @xmath53 , which at a vertex @xmath4 becomes @xmath54 where @xmath55 is the degree of the vertex @xmath4 .
we will say that two harmonic eigenvectors @xmath56 and @xmath57 are orthogonal if @xmath58 , i.e. , if the corresponding eigenvectors are orthogonal .    before we begin the proof of theorem  [ thm : biswitch ] it will be useful to make some observations about bipartite graphs that will come up in the proof .    [
lem : help ] let @xmath59 be a bipartite graph on the vertices @xmath20 where all edges go between @xmath21 and @xmath22 and the vertices in @xmath21 all have degree @xmath23 .
further , if @xmath60 is an eigenvector for the normalized laplacian then we can write it as @xmath61 where @xmath62 is the vector @xmath60 restricted to the vertices of @xmath21 , and similarly for @xmath63 . then the following hold :    * if @xmath64 , i.e. , @xmath60 is an eigenvector for the eigenvalue @xmath41 , then @xmath65 and @xmath66 . in other words
, we can divide the eigenspace for the eigenvalue @xmath41 between those which are nonzero only on the vertices of @xmath21 and nonzero only on the vertices of @xmath22 . *
if @xmath61 is an eigenvector associated with @xmath50 , then @xmath67 is an eigenvector associated with @xmath68 . *
if @xmath69 are a set of orthogonal harmonic eigenvectors of @xmath59 , then for @xmath70 we have @xmath71 is orthogonal to @xmath72 and @xmath73 is orthogonal to @xmath74 . in particular , for @xmath70 the sum of the entries of @xmath71 is @xmath5 .
parts ( a ) and ( b ) easily follow from . for ( c ) we note that @xmath75 will be orthogonal to the vectors @xmath76 and @xmath77 , also since the vertices in @xmath21 are regular the result on the sum of the entries of @xmath71 follows by the definition of orthogonality .
let @xmath69 be a full set of orthogonal harmonic eigenvectors of @xmath18 associated with @xmath78 respectively .
let @xmath79 denote the degree of a vertex in @xmath21 for the graph @xmath24 and @xmath23 the degree in @xmath18 .
then for @xmath70 we have that @xmath80 is a harmonic eigenvector for @xmath24 associated with eigenvalue @xmath81 . to see this we need to examine what happens for each vertex @xmath4 in @xmath24 using .
* for @xmath82 : then all of the entries of the vertex and its neighbors are all @xmath5 and so trivially holds . * for @xmath83 : all of the neighbors in @xmath1 and @xmath29
are @xmath5 while the sum of the entries in @xmath21 will be @xmath5 by lemma  [ lem : help](c ) , so both sides again are @xmath5 and holds . * for @xmath84
: we note that all of the nonzero elements of @xmath85 adjacent to @xmath4 are in @xmath22 , and further we have that for the harmonic eigenvector @xmath75 for @xmath18 that @xmath86 so we have for @xmath24 that @xmath87 so again holds .
* for @xmath88 : we proceed similarly as we did for vertices in @xmath21 and we get @xmath89    next we note by lemma  [ lem : help](c ) that @xmath72 and @xmath74 are orthogonal to each of the vectors @xmath85 for @xmath70 . on @xmath18 ,
the dimension of the space orthogonal to the @xmath85 for @xmath70 restricted to @xmath18 is two . a basis for the orthogonal complement of @xmath90 is @xmath91 ,
so any harmonic vector orthogonal to all of these _ must _ be a linear combination of @xmath72 and @xmath74 .
in particular , if we let @xmath92 be a harmonic eigenvector of @xmath24 that is orthogonal to all of the @xmath85 for @xmath70 , then @xmath92 restricted to @xmath21 and @xmath22 is orthogonal to the @xmath85 , so we have @xmath93 for some constants @xmath94 and @xmath95 .
in other words , the harmonic eigenvector @xmath92 is constant on the vertices of @xmath21 and constant on the vertices of @xmath22 .
everything that we have done for @xmath24 carries over to @xmath32 .
in particular , since @xmath18 and @xmath19 are cospectral , then the @xmath96 found by generalizing the harmonic eigenvectors of @xmath18 and @xmath19 remain the same .
furthermore , every other harmonic eigenvector orthogonal to the @xmath85 must be constant on the vertices of @xmath21 and the vertices of @xmath22 .
to finish off the first claim of the theorem , we now only need to observe that _ any _ harmonic eigenvector which is orthogonal to all of the @xmath85 in @xmath24 is also a harmonic eigenvector for @xmath32 for the same eigenvalue . again to see this we only need to consider what happens for each vertex .    * for @xmath97 : then all of the entries of the vertex and its neighbors are the same for both graphs , so trivially holds . * for @xmath84 :
all of the neighbors in @xmath29 are the same , and we have the exact same number of neighbors in @xmath22 as before with the same value and so again holds . * for @xmath88 : for the graph @xmath24 , let @xmath98 be the fixed value of the vertices in @xmath21 .
we have @xmath99 of course the @xmath55 terms cancel and we are left with @xmath100 .
so even though the degree of the vertex in @xmath22 might change , it will have no effect on this relationship and so again holds .    in summary
, we were able to find @xmath101 harmonic eigenvectors for @xmath24 and @xmath32 that gave the same set of eigenvalues .
for any other harmonic eigenvector orthogonal to these , the same harmonic eigenvector worked for both graphs and so the remaining set of eigenvalues also agreed .
so we can conclude that @xmath24 and @xmath32 are cospectral with respect to the normalized laplacian .
we now turn to the second statement of the theorem .
let @xmath79 be the degree of a vertex in @xmath21 in the graph @xmath35 .
so let @xmath69 be a full set of orthogonal harmonic eigenvectors of @xmath18 associated with @xmath78 respectively .
first , we consider the harmonic eigenvectors of the graph @xmath18 in the eigenspace corresponding to @xmath34 . by lemma  [ lem : help](a )
we can assume that the vectors @xmath75 are either of the form @xmath71 or @xmath73 .
let the corresponding harmonic eigenvector of @xmath35 be @xmath102 ,
i.e. , we simply expand @xmath75 to be zero outside of @xmath18 .
then we have two cases .    *
if @xmath103 : in this case , we know there are no problems for the vertices in @xmath1 ( every term in is @xmath5 ) or @xmath22 ( since reduces to what was done in @xmath18 ) . for the vertices in @xmath29 we can use lemma  [ lem : help](c ) to see that both sides of are @xmath5 .
finally , suppose that @xmath4 is a vertex in @xmath21 .
then we have @xmath104 in particular , we have that @xmath105 is an eigenvalue , i.e. , that @xmath106 is an eigenvalue , for this harmonic eigenvector . * if @xmath107 : in this case we similarly know there are no problems for vertices in @xmath1 , @xmath29 and @xmath21 . for vertices in @xmath22 , reduces to what we had in @xmath18 and so we can conclude that this is a harmonic eigenvector for the eigenvalue @xmath108 .    now suppose that @xmath109 is associated with @xmath110 for @xmath18 .
then we now will create two harmonic eigenvectors for @xmath35 , namely @xmath111 these will be associated with the eigenvalues of @xmath112 respectively , for @xmath35 .
some simple computations show that the following relationships hold : @xmath113    ( before we proceed to the next step of showing that these are indeed harmonic eigenvectors and eigenvalues we first make the observation that we are not creating more harmonic eigenvectors than we had before .
this is because in lemma  [ lem : help](b ) we can pair up harmonic eigenvectors and eigenvalues . in particular , we would have generated the same new harmonic eigenvectors and eigenvalues if we had used @xmath114 and @xmath115 .
so really we are taking pairs of eigenvectors and eigenvalues to new pairs of eigenvectors and eigenvalues . )    to verify that these are harmonic eigenvectors with the specified eigenvalues we need to examine what happens for each vertex @xmath4 in @xmath35 using .
we will step through @xmath116 and @xmath117 , the arguments for @xmath118 and @xmath119 are the same .    * for @xmath82 : then all of the entries of the vertex and its neighbors are all @xmath5 and so trivially holds . * for @xmath83 : all of the neighbors in @xmath1 and @xmath29 are @xmath5 while the sum of the entries in @xmath21 will be @xmath5 by lemma  [ lem : help](c ) , so both sides again are @xmath5 and holds . * for @xmath84
: we note that all of the nonzero elements of @xmath85 adjacent to @xmath4 are in @xmath21 and @xmath22 , and further we have that for the harmonic eigenvector @xmath75 for @xmath18 that @xmath120 so we have for @xmath35 that @xmath121 so again holds . * for @xmath88 : we proceed similarly as we did for vertices in @xmath21 and we get @xmath122    the remainder of this case now proceeds as before .
namely , everything that we did for @xmath35 carries over for @xmath36 ( counting multiplicity of eigenvalues ) .
further , any other harmonic eigenvector orthogonal to the ones given must be constant on @xmath21 and @xmath22 and so any other harmonic eigenvector which works for @xmath35 also works for @xmath36 .
therefore we can conclude that they have the same set of eigenvalues , i.e. , are cospectral with respect to the normalized laplacian .
the graphs generated by theorem  [ thm : biswitch ] will generally not give cospectral graphs for other matrices .
so it is instructive to examine the proof and try to understand the point at which the fact that we were using the normalized laplacian as compared to some other matrix came into play .
the key is understanding how remains true even when the degrees of vertices in @xmath22 change  if the degree of a vertex in @xmath22 changes , the @xmath55 term on the right side of will change proportionally to the sum on the left so that the equality still holds .    in the case when the degrees of vertices do not change , for example , when @xmath18 and @xmath19 are both @xmath123-biregular , then the proof generalizes .
recall that a bipartite graph is @xmath123-biregular if the vertices can be partitioned into @xmath20 where all edges go between @xmath21 and @xmath22 , the vertices in @xmath21 all have degree @xmath23 , and the vertices in @xmath22 have degree @xmath47 .
we now have the following theorem .
[ thm : biregswitch ] let @xmath18 and @xmath19 be bipartite , cospectral with respect to the adjacency matrix , @xmath123-biregular graphs on the vertex set @xmath20 such that all vertices in @xmath21 have degree @xmath23 and all vertices in @xmath22 have degree @xmath47 and edges go between @xmath21 and @xmath22 .
further , the dimension of the eigenspace associated with @xmath124 for the adjacency matrix intersected with the subspace of vectors which are nonzero only on @xmath21 is the same for both @xmath18 and @xmath19 .
let @xmath24 be a graph on the vertices @xmath25 where @xmath26 $ ] is an arbitrary graph ; @xmath27 $ ] and @xmath28 $ ] have no edges ; there are no edges going between @xmath1 and @xmath21 , between @xmath1 and @xmath22 , or between @xmath29 and @xmath22 ; @xmath30 $ ] is a complete bipartite graph ; and @xmath31=p_1 $ ] .
the graph @xmath32 is defined similarly except that @xmath33=p_2 $ ]
. then @xmath24 and @xmath32 are cospectral with respect to the matrix @xmath0 for @xmath125 arbitrary .
similarly , the graphs @xmath35 and @xmath36 are also cospectral with respect to @xmath0 for @xmath125 arbitrary , where @xmath35 only differs from @xmath24 in that @xmath37 $ ] is the complete graph and similarly @xmath36 only differs from @xmath32 in that @xmath38 $ ] is the complete graph .
note when @xmath126 then @xmath127 , when @xmath128 then @xmath129 , and when @xmath130 then @xmath131 , so that any such pairs of graphs are cospectral with respect to all of the matrices @xmath1 , @xmath12 , and @xmath13 .
further , these are also cospectral for @xmath132 since if @xmath133 , we have @xmath134 since the two graphs have the same degree sequence , the product @xmath135 is the same , so the characteristic polynomials also are the same .    in order to use this theorem , we must find two @xmath123-biregular bipartite graphs which are cospectral and for which the dimensions of the eigenspace of @xmath5 restricted to @xmath21 agree on both graphs .
there is one special case for which this is much easier , namely the theorem does not prohibit the possibility that @xmath18 and @xmath19 are the same graph ( so trivially are cospectral so we only are reduced to checking the dimension of the eigenspaces ) .
this is not so trivial as it might sound since we are distinguishing the two parts of the bipartite graph in the proof of the theorem , so that while @xmath18 and @xmath19 are the same graph we are attaching them in two different ways and so the resulting graphs might not be isomorphic .    as an example of this
consider the graphs shown in figure  [ fig:2bi ] .
these are both the same graph ( just flipped ) and a simple check shows that the dimension of the eigenspace associated with @xmath5 and restricted to @xmath21 is @xmath136 in both graphs .
therefore theorem  [ thm : biregswitch ] applies , and any two graphs to which these two graphs are attached are cospectral . but note that in @xmath19 there are two vertices in @xmath22 with the same neighbors ( marked ) , but there are no such pairs in @xmath18 .
so the resulting graphs are _ not _ isomorphic as long as @xmath137 .
so for example if we consider the graphs @xmath35 and @xmath36 where @xmath138 and we induce a complete graph on the vertices of @xmath21 , then theorem  [ thm : biregswitch ] shows that the resulting graphs are cospectral , and again since there are two vertices with degree @xmath139 sharing common neighbors in one graph but not the other they are non - isomorphic .
these graphs are shown in figure  [ fig : bicom ] .
these graphs also happen to be complements of one another , and so give an example of a graph which is cospectral with its complement but not self - complementary .
( note the graph shown in figure  [ fig:1c ] is another example of a graph which is cospectral with its complement , but in that case it is only cospectral with respect to @xmath13 and not with respect to @xmath1 , @xmath12 or @xmath9 . )
the proof for theorem  [ thm : biregswitch ] will proceed similar to the proof for theorem  [ thm : biswitch ] , so we will skip some of the routine computations and provide an outline of the proof . note that if @xmath60 is an eigenvalue associated with eigenvalue @xmath50 , i.e. , @xmath140 , then at a vertex @xmath4 we have @xmath141 we first start with some simple properties of bipartite graphs that will be helpful in the proof .
the proofs are similar to lemma  [ lem : help ] and we will omit them .    [
lem : bihelp ] let @xmath59 be a bipartite @xmath123-biregular graph on the vertices @xmath20 where all edges go between @xmath21 and @xmath22 .
further , if @xmath60 is an eigenvector for the adjacency matrix , then we can write it as @xmath61 where @xmath62 is the vector @xmath60 restricted to the vertices of @xmath21 , and similarly for @xmath63 . then the following hold :    * if @xmath142 , i.e. , @xmath60 is an eigenvector for the eigenvalue @xmath5 , then @xmath143 and @xmath144 .
in other words , we can divide the eigenspace for the eigenvalue @xmath5 between those which are nonzero only on the vertices of @xmath21 and nonzero only on the vertices of @xmath22 . *
if @xmath61 is an eigenvector associated with @xmath50 , then @xmath67 is an eigenvector associated with @xmath145 .
* if @xmath146 are a set of orthogonal eigenvectors of @xmath59 , then for @xmath70 we have @xmath71 is orthogonal to @xmath72 and @xmath73 is orthogonal to @xmath74 . in particular , for @xmath70 the sum of the entries of @xmath71 is @xmath5
.    we will first consider the case when the vertices on @xmath21 induce an empty graph .
let @xmath146 be a full set of orthogonal eigenvectors of @xmath18 for the adjacency matrix associated with the eigenvalues @xmath147 respectively .
let @xmath79 denote the degree of a vertex in @xmath21 in @xmath24 , and note by construction that @xmath47 will be the degree of a vertex in @xmath22 in @xmath24 .
we now show how to extend each of these eigenvectors to an eigenvector for @xmath24 .
we have the following :    * for an eigenvector of the form @xmath148 ( i.e. , only nonzero on the vertices of @xmath21 ) and @xmath124 then @xmath149 is an eigenvector for @xmath0 on @xmath24 associated with the eigenvalue of @xmath150 . * for an eigenvector of the form @xmath151 ( i.e. , only nonzero on the vertices of @xmath22 ) and @xmath124 then @xmath152 is an eigenvector for @xmath0 on @xmath24 associated with the eigenvalue of @xmath153 . * for an eigenvector @xmath61 associated with an eigenvalue @xmath154 , we can construct two new eigenvectors , namely @xmath155 these are associated with the eigenvalues of @xmath0 for @xmath24 @xmath156 respectively .
+ ( note , as in the proof of theorem  [ thm : biswitch ] , that we are not creating too many new eigenvectors and eigenvalues , i.e. , by lemma  [ lem : bihelp](b ) we can pair up the eigenvalues and eigenvectors .
in particular , @xmath61 and @xmath50 will create the same new eigenvectors and eigenvalues as @xmath157 and @xmath145 . )    the verification that each of these is an eigenvector / eigenvalue pair for @xmath24 reduces to verifying for each vertex , similarly as was done in theorem  [ thm : biswitch ] .
we will skip these routine computations .
next we note by lemma  [ lem : bihelp](c ) that @xmath72 and @xmath74 are orthogonal to each of the new eigenvectors for @xmath70 .
since the dimension of the space orthogonal to the new eigenvectors restricted to @xmath18 is two , any vector orthogonal to all of these new vectors _ must _ be a linear combination of @xmath72 and @xmath74 when restricted to @xmath18 .
in particular , if we let @xmath92 be an eigenvector of @xmath24 that is orthogonal to all of the new eigenvectors , then @xmath92 restricted to @xmath21 and @xmath22 must be of the form @xmath93 for some constants @xmath94 and @xmath95 . in other words ,
the eigenvector @xmath92 is constant on the vertices of @xmath21 and constant on the vertices of @xmath22 .
everything that we have done for @xmath24 carries over to @xmath32 .
in particular , since @xmath18 and @xmath19 are cospectral , then the newly found @xmath96 are the same for both graphs , and further , every other eigenvector orthogonal to the ones found must be constant on the vertices of @xmath21 and the vertices of @xmath22 .
to finish off the first case of the theorem , we now only need to observe that _ any _ eigenvector which is orthogonal to all of the new eigenvectors in @xmath24 is also an eigenvector for @xmath32 for the same eigenvalue .
this again is done in the exact same way as in theorem  [ thm : biswitch ] and we will skip the computations here .    in summary , we were able to find @xmath101 eigenvectors of @xmath0 for @xmath24 and @xmath32 that gave the same set of eigenvalues .
for any other eigenvector orthogonal to these , the same eigenvector worked for both graphs and so the remaining set of eigenvalues also agreed .
so we can conclude that @xmath24 and @xmath32 are cospectral with respect to the matrix @xmath0 .
now we turn to the case when the vertices on @xmath21 induce a complete graph .    using the same notation as in the previous case
, we will extend the eigenvectors of @xmath18 for the adjacency matrix to eigenvectors of @xmath35 for @xmath0 .
we have the following :    * for an eigenvector of the form @xmath148 ( i.e. , only nonzero on the vertices of @xmath21 ) and @xmath124 then @xmath149 is an eigenvector for @xmath0 on @xmath35 associated with the eigenvalue of @xmath158 . * for an eigenvector of the form @xmath151 ( i.e. , only nonzero on the vertices of @xmath22 ) and @xmath124 then @xmath152 is an eigenvector for @xmath0 on @xmath35 associated with the eigenvalue of @xmath153 . * for an eigenvector @xmath61 associated with an eigenvalue @xmath154 , we can construct two new eigenvectors , namely @xmath159 these are associated with the eigenvalues of @xmath0 for @xmath35 @xmath160 respectively .
+ ( note , as in the proof of theorem  [ thm : biswitch ] , that we are not creating too many new eigenvectors and eigenvalues , i.e. , by lemma  [ lem : bihelp](b ) we can pair up the eigenvalues and eigenvectors .
in particular , @xmath61 and @xmath50 will create the same new eigenvectors and eigenvalues as @xmath157 and @xmath145 . )    the verification that each of these is an eigenvector / eigenvalue pair for @xmath35 reduces to verifying for each vertex , similarly as was done in theorem  [ thm : biswitch ] .
we will skip these routine computations .
the remainder of this case now proceeds as before .
namely , everything that we did for @xmath35 carries over for @xmath36 ( counting multiplicity of eigenvalues ) .
further , any other eigenvector orthogonal to the ones given must be constant on @xmath21 and @xmath22 , so any other eigenvector for @xmath0 which works for @xmath35 also works for @xmath36 .
therefore we can conclude that they have the same set of eigenvalues , i.e. , are cospectral with respect to @xmath0 .
we can use theorem  [ thm : biswitch ] to construct large families of mutually cospectral graphs for the normalized laplacian that need not be bipartite nor have the same degree sequences .
let @xmath161 , @xmath23 be positive integers such that @xmath162 is a partition of @xmath161 into @xmath23 positive integer parts .
let @xmath163 be the graph on @xmath164 vertices @xmath165 , where the @xmath166 induce a complete graph , each @xmath167 is only adjacent to exactly @xmath168 vertices from @xmath169 , and each @xmath166 is adjacent to exactly one of the @xmath167 .
the collection of all @xmath163 for partitions @xmath170 give a family of cospectral graphs for the normalized laplacian , i.e. , the spectrum is completely determined for the normalized laplacian by @xmath161 and @xmath23 .
we have dubbed these graphs `` fuzzy balls '' since they consist of a complete graph ( the @xmath166 ) with some extra vertices attached ( the fuzz ) . in figure
[ fig : fuzzy ] , we have shown the cospectral family that corresponds to @xmath171 and @xmath172 .     and @xmath172 .
]    to see that these are all cospectral , we observe by theorem  [ thm : biswitch ] that @xmath163 is cospectral with @xmath173 by using the pair of cospectral graphs @xmath174 and @xmath175 where @xmath21 induces a complete graph . therefore , applying this idea @xmath176 times , we can conclude that for each partition of @xmath161 into @xmath23 parts we have @xmath163 is cospectral with @xmath177 .
in particular , they are all mutually cospectral .
let @xmath161 , @xmath23 be positive integers such that @xmath162 is a partition of @xmath161 into @xmath23 positive integer parts .
let @xmath178 be the graph on @xmath179 vertices @xmath180 , @xmath181 , @xmath182 , @xmath183 , @xmath184 , @xmath182 , @xmath185 , where the @xmath166 and @xmath180 induce a star graph with @xmath180 as the central vertex , each @xmath167 is only adjacent to exactly @xmath168 vertices from @xmath169 , and each @xmath166 is adjacent to exactly one of the @xmath167 . the collection of all @xmath178 for partitions @xmath170 give a family of cospectral graphs , i.e.
, the spectrum is completely determined by @xmath161 and @xmath23 .
we have dubbed these graphs `` inflated stars '' . in figure
[ fig : stars ] , we have shown the cospectral family that corresponds to @xmath171 and @xmath172 .     and @xmath172 . ]    again , to see that these are all cospectral , we observe by theorem  [ thm : biswitch ] that the graph @xmath178 is cospectral with @xmath186 .
therefore , applying this idea @xmath176 times , we can conclude that for each partition of @xmath161 into @xmath23 parts , we have @xmath178 is cospectral with @xmath187 . in particular
, they are all mutually cospectral .    in both of these examples
, we are using partitions to form large families .
the number of partitions grows subexponentially with @xmath161 , however it is possible to use theorems  [ thm : biswitch ] or  [ thm : biregswitch ] to construct large mutually cospectral families that grow exponentially with @xmath161 , either with respect to the normalized laplacian , or more generally with respect to @xmath0 .
[ thm : nlfam ] for @xmath161 large , there exists a family of @xmath188 non - isomorphic , mutually cospectral graphs with respect to the normalized laplacian .
find a graph @xmath16 on @xmath189 vertices that has a trivial automorphism group ( for @xmath190 this is easy to do ) .
for each vertex @xmath4 of @xmath16 , attach one of the two inflated stars shown in figure  [ fig : nlfam ] by identifying @xmath4 with the vertex marked @xmath4 in the inflated star .
this constructs @xmath191 graphs on @xmath192 vertices .     in theorem  [ thm : nlfam ] .
]    we note that for any one of these graphs we can easily recover @xmath16 and identify which attachment was done at @xmath4 and so these graphs are all non - isomorphic .
further , by theorem  [ thm : biswitch ] changing our choice of widget at any single vertex does not change the spectrum and so all of the graphs are cospectral .
we have not tried to optimize the construction .
for example , if instead of using the family of inflated stars with @xmath193 and @xmath194 we used the family of inflated stars with @xmath195 and @xmath196 ( of which there are @xmath197 such non - isomorphic members ) , then we get a family with rate of growth @xmath198 which is faster than @xmath199 .
we similarly have the following result for graphs cospectral with respect to @xmath0 .
[ thm : atdfam ] for @xmath161 large , there exists a family of @xmath200 non - isomorphic , mutually cospectral graphs with respect to the matrix @xmath0 .
find a graph @xmath16 on @xmath201 vertices that has a trivial automorphism group ( for @xmath202 this is easy to do ) .
for each vertex @xmath4 of @xmath16 attach one of the two graphs shown in figure  [ fig : atdfam ] , which are taken from figure  [ fig:2bi ] .
this constructs @xmath203 graphs on @xmath204 vertices .     in theorem  [ thm : atdfam ] .
]    we note that for any one of these graphs we can easily recover @xmath16 and identify which attachment was done at @xmath4 , so these graphs are all non - isomorphic .
further , by theorem  [ thm : biregswitch ] changing our choice of widget at any single vertex does not change the spectrum and so all of the graphs are cospectral .
the method of swapping bipartite subgraphs was discovered by examining all cospectral pairs of graphs on at most eight vertices ( all cospectral pairs on at most seven vertices has previously appeared in tan @xcite ) , which included several simple examples of this type .
the set of all cospectral graphs on at most eight vertices can be easily constructed in sage ( see appendix  [ sagecode ] ) .    by looking at the small cases
, we also discovered some new graphs which are regular , but are cospectral for the normalized laplacian with graphs which are not regular .
we have already seen that the four - cycle and the eight - cycle are members of cospectral pairs ( see figures  [ fig:1d ] and [ fig : other1 ] ) . also , @xmath205 is cospectral with @xmath42 for @xmath206 .
another example of a graph which is regular , but is cospectral with a graph which is not regular , is shown in figure  [ fig : regular ] .
these examples show that the normalized laplacian can not , in general , detect whether a graph is regular .
as noted in the introduction , the number of edges in a graph is not uniquely determined by the spectrum of the normalized laplacian . as a consequence ,
it is possible for a graph to be cospectral with one of its subgraphs ; examples of this are shown in figure  [ fig : sub ] .
further examples of graphs which are cospectral with respect to the normalized laplacian but have differing number of edges are shown in figure  [ fig : dif ] ( there is currently no known method to generate cospectral graphs with differing number of edges other than using complete bipartite graphs ) .
this also indicates the difficulty in counting the number of cospectral graphs . while for the adjacency and the combinatorial laplacian
, we could first subdivide the graphs according to the number of edges and work with this coarsening when finding cospectral pairs ( see @xcite ) , this is no longer possible for the normalized laplacian .
it would be interesting to find techniques for constructing cospectral graphs for the normalized laplacian which differ in the number of edges .
similarly , it would be interesting to investigate the cospectral graphs for the signless laplacian and see if there are any constructions unique to forming cospectral graphs for that matrix ( similar to what theorem  [ thm : biswitch ] is for the normalized laplacian ) .
99 andries e. brouwer and edward spence , _ cospectral graphs on 12 vertices _ , electronic journal of combinatorics * 16 * ( 2009 ) , # n20 , 3pp .
steve butler , _ cospectral graphs for both the adjacency and normalized laplacian matrices _ , linear and multilinear algebra * 58 * ( 2010 ) , 387390 .
michael cavers , the normalized laplacian matrix and general randic index of graphs , ph.d .
thesis , university of regina , 2010 .
fan chung , _ spectral graph theory _
, cbms lecture notes , ams , providence , ri , 1997 .
drago m. cvetkovi , _ graphs and their spectra _ , univ .
* 354 - 356 * ( 1971 ) , 150 .
drago m. cvetkovi , michael doob and horst sachs , _ spectra of graphs : theory and application _ , academic press , new york , ny , 1979 .
drago m. cvetkovi and slobodan k. simi , _ towards a spectral theory of graphs based on the signless laplacian i _ , institut mathmatique . publications .
nouvelle srie * 85(99 ) * ( 2009 ) , 19 - 33 .    drago m. cvetkovi and slobodan k. simi , _ towards a spectral thory of graphs based on the signless laplacian ii _ , linear algebra and its applications * 432 * ( 2010 ) , 2257 - 2272 .    drago m. cvetkovi and slobodan k. simi , _ towards a spectral thory of graphs based on the signless laplacian iii _ ,
applicable analysis and discrete mathematics * 4 * ( 2010 ) , 156 - 166 .
edwin r. van dam , _ graphs with few eigenvalues : an interplay between combinatorics and algebra _ , ph.d .
thesis , tilburg university , center for economic research dissertations series , no .
20 , 1996 .
edwin r. van dam and willem h. haemers , _ which graphs are determined by their spectrum ? _ , linear algebra and its applications * 373 * ( 2003 ) , 241272 .
edwin r. van dam , willem h. haemers and jack h. koolen , _ cospectral graphs and the generalized adjacency matrix _ , linear algebra and its applications * 423 * ( 2007 ) , 3341 .
chris d. godsil and brendan mckay , _ products of graphs and their spectra _ , combinatorial mathematics , iv ( proc .
fourth australian conf .
adelaide , adelaide , 1975 ) , lecture notes in math .
560 , springer , berlin , 1976 , 6172 .
chris d. godsil and brendan mckay , _ constructing cospectral graphs _ , aequationes mathematicae * 25 * ( 1982 ) , 257268 .    willem h. haemers and edward spence , _ enumeration of cospectral graphs _ , european journal of combinatorics * 25 * ( 2004 ) , 199211 .
charles r. johnson and morris newman , _ a note on cospectral graphs _ ,
journal of combinatorial theory , ser .
b * 28 * ( 1980 ) , 96103 .
russell merriss , _ large families of laplacian isospectral graphs _ , linear and multilinear algebra * 43 * ( 1997 ) , 201205 .
johan jacob seidel , _ graphs and two - graphs _ , proceedings of the fifth southeastern conference on combinatorics , graph theory and computing ( florida atlantic univ . , boca raton , fla . , 1974 ) , congressus numerantium , no .
x , utilitas math . ,
winnipeg , man .
, 1974 , 125143 .
william a. stein , et  al .
, _ sage mathematics software ( version 4.6.1 ) _ , the sage development team , 2011 , http://www.sagemath.org .
jinsong tan , _ on isospectral graphs _ , interdisciplinary information sciences * 4 * ( 1998 ) , 117124 .
richard wilson and ping zhu , _ a study of graph spectra for comparing graphs and trees _
, journal of pattern recognition , * 41 * ( 2008 ) , 28332841 .
the code listing below shows how to use sage version 4.5.2 @xcite to generate and check cospectral graphs with 8 vertices . in the listing below , we specifically generate cospectral graphs with respect to the normalized laplacian .
the lines starting with the number sign `` # '' are comment lines .


the term _ cardinal interpolation _ refers to interpolation of data given at the integer lattice ( or multi - integer lattice in higher dimensions ) .
it was i. j. schoenberg s work on cardinal spline interpolation that brought about an intense study of the subject .
many avenues of study have been explored , including forming interpolation operators from translates of certain radial basis functions ( rbfs ) .
works by buhmann , baxter , riemenschneider , and sivakumar @xcite ( see also @xcite and references therein ) have explored many cardinal interpolation operators of this type .
some of the radial basis functions that have been considered are the gaussian kernel , the thin plate spline , the hardy multiquadric , and the inverse multiquadric .
radial basis cardinal interpolation also enjoys a strong connection with classical sampling theory , as evidenced by much of the aforementioned literature .
this connection continues to be explored in recent developments by the second author @xcite , and by parts of this article . as this is one of the most interesting aspects of cardinal interpolation
, we provide some of the motivation .
recall the one - dimensional whittaker  kotelnikov
shannon sampling theorem , which states that a bandlimited function , @xmath10 , ( say with band size @xmath11 ) can be recovered uniformly by the series @xmath12 , where the @xmath13 function is suitably defined so that it takes value 1 at the origin , and 0 at all the other integers .
one observation about this series is that the @xmath13 function decays slowly ( as @xmath14 ) , and so to approximate the series by truncation for example , one would have to use quite a lot of data points of @xmath10 to get a reasonable degree of accuracy .
however , there is a way of approximating the @xmath13 series above : we seek to replace the @xmath13 function with a so - called _ fundamental function _ ( or _ cardinal function _ ) , @xmath15 , that preserves the property that @xmath15 takes value 1 at the origin and 0 at all other integers .
we then form a function @xmath16 the trade - off here is that while @xmath17 is not pointwise equal to the function @xmath10 , it does interpolate @xmath10 at the integer lattice , and moreover , the fundamental function @xmath15 may be constructed so that it decays much more rapidly than the @xmath13 function .
precisely , one may construct a fundamental function from a given radial basis function , @xmath18 , which has the form @xmath19    in the case of the gaussian kernel , @xmath20 , the fundamental function decays exponentially , whereas the fundamental function for the hardy multiquadric , @xmath21 , decays as @xmath22 .
so we give up the pointwise equality of the wks sampling theorem in exchange for a series that converges more rapidly , while also ensuring that @xmath17 is close to @xmath10 in the @xmath23 norm .
this paper primarily considers the fundamental functions and cardinal interpolation operators associated with general multiquadrics , @xmath0 , which have thus far only been considered for certain exponents @xmath24 .
interpolation with fundamental functions has too long a history to recount here ; however , @xcite offers a good introduction using radial basis functions .
using below as a starting point is especially popular since it allows one to solve problems in the fourier transform domain .
many authors have used similar techniques for various radial basis functions . in @xcite , riemenschneider and
sivakumar proved several results pertaining to the gaussian .
multiquadric cardinal interpolation has been studied in a similar way by buhmann and micchelli @xcite , baxter @xcite , baxter and sivakumar @xcite , riemenschneider and sivakumar @xcite , among others .
compactly supported radial basis functions have been studied by buhmann @xcite and wendland @xcite .
the rest of the paper is laid out as follows .
section [ secbasic ] provides the necessary preliminaries and a discussion of applications and calculations of the fundamental functions ; section [ secbaxteranalogue ] shows recovery results for cardinal interpolation of bandlimited functions in any dimension via interpolants of the form discussed above .
section [ secpropfundamental ] contains decay rates and information about the univariate fundamental functions associated with the general multiquadrics for a broad range of exponents . in section [ seceastjournal ] , we consider the cardinal interpolation operators acting on data in traditional sequence spaces and calculate decay rates , bounds on the operator norms , and also explore some convergence properties in terms of the parameter @xmath25 .
section [ secexamples ] provides some interesting concrete examples based on the theoretical results from the previous section .
finally , section [ secproofs ] provides the technical proofs of the statements in section [ secpropfundamental ] .
if @xmath26 is an interval , then let @xmath27 , @xmath28 , be the usual lebesgue space over @xmath29 with its usual norm .
if no set is specified , we mean @xmath30 .
similarly , denote by @xmath31 the usual sequence spaces indexed by the set @xmath32 ; if no index set is given , we refer to @xmath33
. we will use @xmath34 to denote the natural numbers including 0 .
let @xmath35 be the space of schwartz functions on @xmath36 , that is the collection of infinitely differentiable functions @xmath18 such that for all multi - indices @xmath24 and @xmath37 , @xmath38    the fourier transform of a schwartz function @xmath18 is given by @xmath39 thus the inversion formula is @xmath40 in the event that these formulas do not hold , then the fourier transform should be interpreted in the sense of tempered distributions .
recall that if @xmath10 is a tempered distribution , then its fourier transform is the tempered distribution defined by @xmath41 , @xmath42 .
let @xmath43 and @xmath44 be fixed ; then define the @xmath45-dimensional _ general multiquadric _ by @xmath46 where @xmath47 denotes the euclidean distance on @xmath36 .
if @xmath48 , the generalized fourier transform of @xmath4 is given by the following ( see , for example , ( * ? ? ?
* theorem 8.15 ) ) : @xmath49 where @xcite @xmath50 @xmath51 is called the modified bessel function of the second kind .
we note that both @xmath4 and its fourier transform are radial .
it is also clear from the definition that @xmath52 is symmetric in its order ; that is , @xmath53 for any @xmath54 .
if @xmath55 , then the generalized fourier transform of @xmath4 involves a measure and so can not be expressed as a function .
now suppose that @xmath48 is fixed .
to define the _ fundamental function _ associated with the general multiquadric , we first define the following function @xmath56 we will see that @xmath57 , and so the function @xmath58 is well - defined and continuous .
furthermore , we will show that @xmath3 is a _ fundamental function _ , also called a _
cardinal function _ , which means that @xmath59 where @xmath60 is the kronecker delta .
additionally , @xmath3 has the form @xmath61    throughout the paper , we will use @xmath62 to denote an absolute constant due to the use of @xmath25 as the shape parameter of the multiquadric .
the value of the particular constant may change from line to line , and we use subscripts to denote dependence on certain parameters when needed .
interpolation schemes involving fundamental functions as in have been studied for quite some time , and there are many aspects to the theory .
for example , such methods enjoy applications to geoscience @xcite and sampling theory @xcite .
recently , investigations have considered interpolation via radial kernels on manifolds @xcite . for a galerkin type method for solving pdes using meshless interpolation on the sphere ,
see @xcite .
given the widespread applications of radial basis function approximation , it is of import to the computational community to determine stable ways of evaluating the approximants .
consequently , there is a substantial literature dealing with accuracy and stability of different computational methodologies for radial basis function approximation .
we do not claim to list all of these methods , but at least a sampling is in order .
we note that approximating is typically very difficult , especially if one dilates the lattice .
one way around this is the use of indirect computational methods to approximate the rbf interpolant @xcite ( for a discussion specifically related to multiquadrics , see @xcite ) .
another technique involves a change of basis method @xcite , while work by fasshauer and mccourt @xcite uses an eigenfunction decomposition to provide stable reconstruction using gaussians .
another quite promising method has recently been considered in which so - called _ local lagrange functions _ are used to approximate rather than the global ones @xcite .
many of these results revolve around the situation of interpolation at finitely many data sites , which is of a somewhat different nature than we are considering here . for the interested reader
, we also mention that these methods of cardinal interpolation have , at their core , a deep connection to the classical spline theory instigated by schoenberg ( see @xcite and references therein ) and continued by many followers .
the underlying principle is that many of the results in spline interpolation theory have natural analogues via rbfs , and the problem at hand may determine which method is more useful .      here , we provide some brief illustrations of the fundamental functions we have mentioned above .
figure [ figsinc1d ] shows the univariate fundamental function associated with the inverse multiquadric and the sinc function for comparison .
@xmath63{acom_sinc_imq - eps - converted - to.pdf}\ ] ]    @xmath64{multivariatefundamental - eps - converted - to.pdf}\ ] ]    the fundamental function was calculated by truncating the series in and using a fast fourier transform ( fft ) method to approximate @xmath3 .
the same method may be applied for the multivariate version , as figure [ figsinc ] shows .
when an interpolation scheme depends on a parameter , questions of convergence naturally arise . this question has been addressed by several authors . in @xcite , baxter examines the hardy multiquadric , while the gaussian is studied in @xcite by riemenschneider and sivakumar and in @xcite by hangelbroek , narcowich , madych , and ward . in a more general context , ` increasingly flat ' radial basis functions are the focus of driscoll and fornberg in @xcite , while the second author worked with ` regular families ' of cardinal interpolants in @xcite .    in this section , we show that the result obtained by baxter @xcite holds not only for the traditional hardy multiquadric ( corresponding to @xmath65 ) but rather for all @xmath48 .
we consider interpolation of bandlimited ( or paley - wiener ) functions in any dimension , and show that the cardinal interpolant converges to the function as the shape parameter @xmath25 tends to infinity .
general multiquadrics were not considered for quite some time in this setting , but in @xcite convergence results for cardinal interpolation of bandlimited functions were shown for a restricted range of exponents .
however , the analysis there was of a more general nature , so here we show that a more specific analysis yields convergence results for the full range @xmath66    we note that the results of this section have very close ties to classical sampling theory , which studies the reconstruction of certain classes of signals from their discrete samples . as mentioned above
, these considerations lead to alternative methods for approximating the sampling series given by the wks sampling theorem for bandlimited signals .
let @xmath45 be the dimension , and @xmath48 and @xmath44 be fixed .
it is evident from that @xmath67 does not change sign . therefore , @xmath68 for all @xmath69 . from , it is also evident that @xmath70 . to show that @xmath71 we begin with the following lemma .
[ lembaxterphihatdecay ] let @xmath72 , @xmath44 , and @xmath73 . then @xmath74    defining @xmath75 , equations and yield the following series of estimates : @xmath76    we note that if @xmath77 , then @xmath78 , and so we have a purely exponential upper bound .    [ proplchatintegrable ] let @xmath48 and @xmath44 .
then @xmath71 .    first , choose an @xmath79 large . then since @xmath80 for all @xmath81
, we have that @xmath82^d}| { \widehat{l_{\alpha , c}}}(\xi)|d\xi\leq ( 2m)^d.\ ] ] now we need to estimate @xmath83^d}| { \widehat{l_{\alpha , c}}}(\xi)|d\xi.\ ] ]    to do this , we establish a pointwise estimate for @xmath84 .
let @xmath85^d$ ] be fixed .
since @xmath86 is large , there exists some @xmath87 such that @xmath88 .
additionally , there is some constant @xmath89 which depends on @xmath24 and @xmath45 such that if @xmath90 , @xmath91 ( see , for example , ( * ? ? ?
* corollary 5.12 ) ) .
therefore , choose @xmath86 large enough so that for @xmath85^d$ ] , we have @xmath92 .
then if @xmath93 is the constant from lemma [ lembaxterphihatdecay ] , @xmath94 now depending on the sign of @xmath95 , the above expression is minimized by plugging in @xmath96 or @xmath97 for @xmath98 in the appropriate places . consequently , there is a positive constant @xmath99 such that @xmath100 we also find from ( * ? ? ?
* lemma 5.13 ) that for every @xmath101 , @xmath102 consequently , by adjusting @xmath86 if need be so that @xmath103 for @xmath85^d$ ] , we find that there is a positive constant @xmath37 such that @xmath104 .
we conclude that @xmath105^d}| { \widehat{\phi_{\alpha , c}}}(\xi)|d\xi\\ & \leq\beta d^{-1}|\lambda|e^{4\pi c } { \displaystyle\int } _ { { \mathbb{r}}^d\setminus[-m , m]^d}\|\xi\|^{-\alpha-\frac{d}{2}}e^{-c\|\xi\|}d\xi.\end{aligned}\ ] ] the integral on the right hand side above is convergent , and the constants outside are all finite , so @xmath71 .
now we turn our attention to the function @xmath3 .
[ proplfundamental ] let @xmath48 and @xmath44 .
then the function @xmath106 is continuous , square - integrable , and satisfies the interpolatory condition @xmath107 , for every @xmath108 .
proposition [ proplchatintegrable ] implies that @xmath3 is continuous and square - integrable , and indeed that @xmath109 is its fourier transform . to see the interpolatory condition , first define @xmath110^d$ ] .
then we have via the substitution @xmath111 that @xmath112 the interchange of sum and integral in the third line is justified by the dominated convergence theorem , for example .
it is an important observation that much of the cardinal interpolation theory for bandlimited functions revolves around the fact that the fundamental functions converge to the function @xmath113 , which is equivalent to the statement that the fourier transform of the fundamental function converges almost everywhere to the indicator function of the cube @xmath114^d$ ] .
the story is no different here . defining @xmath115 to be the function that takes value 1 whenever @xmath116^d$ ] , and 0 elsewhere
, the following holds .
[ propconvergencetosinc ] let @xmath48 .
then @xmath117 for all @xmath69 such that @xmath118 .
first suppose that @xmath119^d$ ] .
then there exists some @xmath120 such that @xmath121 .
therefore by lemma [ lembaxterphihatdecay ] , @xmath122 consequently , since @xmath67 is of one sign , @xmath123 since the exponent is negative , the limit of the right hand side as @xmath6 is 0 .
therefore , for @xmath119^d$ ] , @xmath124 .
now suppose that @xmath125 .
then for all @xmath126 , @xmath127 .
due to , we may write @xmath128 and therefore it suffices to show that @xmath129 by lemma [ lembaxterphihatdecay ] , @xmath130 the series on the right hand side is convergent and dominated by the convergent series where @xmath25 is replaced by 1 , so @xmath131 as desired . convergence of the series stems from the fact that the series is majorized by @xmath132 , which converges .
we now consider interpolation of bandlimited functions at the lattice @xmath133 by translates of the function @xmath134 .
define the @xmath45-dimensional paley - wiener space by @xmath135^d\}.\ ] ] we begin our analysis with an @xmath23 version of the poisson summation formula :    [ lempsf ] if @xmath136 , then @xmath137 where the second series is convergent in @xmath138 .
[ lemcauchysequence ] let @xmath136 . for @xmath139 , define @xmath140 where @xmath141 for @xmath108 .
then @xmath142 forms a cauchy sequence in @xmath138 .
define @xmath143 via @xmath144 thus , @xmath145 . from lemma [ lempsf ] ,
it is clear that @xmath146 is a cauchy sequence in @xmath147^d$ ] .
so @xmath148^d}|q_m(\xi+2\pi k)-q_\ell(\xi+2\pi k)|^2 \\     & & \hfill\times\left ( { \widehat{l_{\alpha , c}}}(\xi+2\pi k)\right)^2 d\xi\\     \\     & = & { \displaystyle\int}_{[-\pi,\pi]^d}|q_m(\xi)-q_\ell(\xi)|^2 { \underset{k\in { \mathbb{z}}^d } { \displaystyle\sum}}\left ( { \widehat{l_{\alpha , c}}}(\xi+2\pi
k)\right)^2 d\xi\\     \\     & \leq & { \displaystyle\int}_{[-\pi,\pi]^d}|q_m(\xi)-q_\ell(\xi)|^2d\xi\;. \\    \end{array}\ ] ] the interchange of sum and integral is valid by tonelli s theorem , and the last inequality follows from the fact that @xmath149 we also used the fact that for @xmath108 , @xmath150 . we conclude that @xmath151 is a cauchy sequence in @xmath138 because @xmath152^d}$ ] , and the latter is cauchy .
lemmas [ lempsf ] and [ lemcauchysequence ] allow us to define @xmath153 where the series is convergent in @xmath138 . by a periodization argument similar to that in the proof of lemma [ lemcauchysequence ]
, one can show that @xmath154 . thus applying the fourier inversion formula
we see that @xmath155    [ thmbandlimitedconvergence ] let @xmath48 . if @xmath136 , then @xmath156 and @xmath157 uniformly on @xmath36
.    we will first prove uniform convergence .
the proof is the same as in @xcite .
again let @xmath115 be the characteristic function of the cube .
then we see by the inversion formula and the oft - exploited periodization argument , that @xmath158^d}\widehat{f}(\xi ) { \underset{k\in { \mathbb{z}}^d } { \displaystyle\sum}}\left ( { \widehat{l_{\alpha , c}}}(\xi+2\pi k)-i(\xi+2\pi k)\right)\\ & & \hfill\times e^{-i { \left\langlex,\xi+2\pi k\right\rangle}}d\xi.\\    \end{array}\ ] ] therefore , we find that @xmath159^d}|\widehat{f}(\xi)| { \underset{k\in { \mathbb{z}}^d } { \displaystyle\sum}}\left| { \widehat{l_{\alpha , c}}}(\xi+2\pi k)-i(\xi+2\pi k)\right|d\xi\\     \\     & = & \dfrac{1}{(2\pi)^d } { \displaystyle\int}_{[-\pi,\pi]^d}|\widehat{f}(\xi)|\left(1- { \widehat{l_{\alpha , c}}}(\xi)+ { \underset{k\neq0 } { \displaystyle\sum } } { \widehat{l_{\alpha , c}}}(\xi+2\pi k)\right)d\xi.\\    \end{array}\ ] ]    but then by definition , @xmath160 therefore , @xmath161^d}|\widehat{f}(\xi)|(1- { \widehat{l_{\alpha , c}}}(\xi))d\xi.\ ] ] as the integrand is non - negative and bounded by @xmath162^d$ ] , and @xmath163 , the dominated convergence theorem implies that @xmath164 the upper bound is independent of @xmath165 , hence the convergence is uniform .
we now turn to the proof of @xmath23 convergence . by parseval
s identity , it suffices to show that @xmath166 .
this breaks up into two estimates .
we first show this for the cube @xmath114^d$ ] . recall that since @xmath167 is an orthonormal basis for @xmath147^d$ ] , we may write @xmath168 . moreover , @xmath169^d } = \|f(k)\|_{\ell_2 ( { \mathbb{z}}^d)}.\ ] ] thus using , @xmath170^d}^2 & =   { \displaystyle\int}_{[-\pi,\pi]^d}\left| { \underset{k\in { \mathbb{z}}^d } { \displaystyle\sum}}f(k ) ( { \widehat{l_{\alpha , c}}}(\xi)-1)e^{-i { \left\langlek,\xi\right\rangle}}\right|^2d\xi \\ & = { \displaystyle\int}_{[-\pi,\pi]^d}| { \widehat{l_{\alpha , c}}}(\xi)-1|^2\left| { \underset{k\in { \mathbb{z}}^d } { \displaystyle\sum}}f(k)e^{-i { \left\langlek,\xi\right\rangle}}\right|^2d\xi.\end{aligned}\ ] ]    the right hand side is bounded by @xmath171 , and so by the dominated convergence theorem and proposition [ propconvergencetosinc ] , @xmath172^d}=0 $ ] .
now for the rest of the space , for @xmath173 define @xmath174\times\dots\times[-\pi-2\pi l_d,\pi-2\pi l_d]$ ] .
then we see that since @xmath10 is bandlimited , @xmath175^d ) } = \| { \widehat { { \mathscr{i}}_{\alpha , c}f}}\|^2_{l_2 ( { \mathbb{r}}^d\setminus[-\pi,\pi]^d ) } =   { \underset{l\neq 0 } { \displaystyle\sum}}\| { \widehat { { \mathscr{i}}_{\alpha , c}f}}\|_{l_2(q_l)}^2.\ ] ] consequently , @xmath176^d}| { \widehat { { \mathscr{i}}_{\alpha , c}f}}(\xi)|^2d\xi &   =   & { \underset{l\neq 0 } { \displaystyle\sum } } { \displaystyle\int}_{q_l}\left| { \widehat{l_{\alpha , c}}}(\xi ) { \underset{k\in { \mathbb{z}}^d } { \displaystyle\sum}}f(k)e^{-i { \left\langlek,\xi\right\rangle}}\right|^2d\xi\\ \\ & = & { \displaystyle\int}_{[-\pi,\pi]^d } { \underset{l\neq 0 } { \displaystyle\sum}}| { \widehat{l_{\alpha , c}}}(\xi+2\pi l)|^2\left| { \underset{k\in { \mathbb{z}}^d } { \displaystyle\sum}}f(k)e^{-i { \left\langlek,\xi\right\rangle}}\right|^2d\xi,\\ \\   \end{array}\ ] ] by the monotone convergence theorem .    recall that @xmath70 , so @xmath177 , and as calculated above , @xmath178 .
consequently , the integrand is bounded by @xmath179 therefore , the dominated convergence theorem and proposition [ propconvergencetosinc ] imply that @xmath180^d)}=0 $ ] , and the proof is complete .    to illustrate the convergence given by theorem [ thmbandlimitedconvergence ] above , the following figure shows the inverse multiquadric interpolant of the function whose fourier transform is @xmath181 in dimension 1 .
[ figinterp ] @xmath182{interp_mq - eps - converted - to.pdf}\ ] ]
for the rest of the paper , we turn our attentions to the one - dimensional cardinal interpolation operator associated with the general multiquadric .
this section is devoted to the one - dimensional fundamental function @xmath3 , whose fourier transform can be rewritten as @xmath183^{-1}.\ ] ] the proofs of the results in this section are quite technical , so we postpone them until section [ secproofs ] and simply state our conclusions here .    to determine decay rates for @xmath3 , we determine how many derivatives @xmath9 has in @xmath184 , which we accomplish by establishing pointwise estimates .
we begin by fixing @xmath185 , so that our estimates fall into three ranges : @xmath186 , @xmath187 , and the @xmath96-length blocks @xmath188 $ ] for @xmath189 . due to the differing behavior of @xmath190 for positive and negative values of @xmath24
, we must make corresponding distinctions in our calculations .    following the insightful techniques of riemenschneider and sivakumar found in @xcite , we begin by defining some auxiliary functions to aid in the analysis of the fundamental function .
we abbreviate as @xmath191 , where @xmath192 and study the properties of @xmath193 .
[ prop l1 ] suppose that @xmath194 , @xmath195 , @xmath196 , and @xmath197 . if @xmath198 and @xmath199 , then there exists a constant @xmath200 such that @xmath201 as @xmath202 .
this estimate leads to the following bounds on @xmath9 and its derivatives .
[ prop l2 ] suppose that @xmath194 , @xmath195 , @xmath196 , and @xmath197 . if @xmath203 , then there exist constants @xmath204 such that    1 .
@xmath205 whenever @xmath198 , 2 .
@xmath206 whenever @xmath207 $ ] , and 3 .
@xmath208 whenever @xmath209 $ ] and @xmath189 ,    where @xmath210 as @xmath202 .
these pointwise estimates yield the following result .
[ thm l1 ] suppose that @xmath194 , @xmath196 , and @xmath197 .
if @xmath203 , then there exists a constant @xmath211 such that @xmath212    using standard arguments , we have the following estimate for the growth rate of @xmath3 .
[ cor l3 ] if @xmath194 and @xmath213 , then @xmath214 as @xmath215 .
analogous estimates can be made for the case that @xmath216 .
of interest to us are the following results .    [ thm l2 ]
suppose that @xmath216 , @xmath196 , and @xmath197 .
if @xmath217 , then there exists a constant @xmath211 such that @xmath218    [ cor l4 ] if @xmath216 and @xmath213 , then @xmath219 as @xmath215 .
it turns out that the poisson kernel , which is the case @xmath220 , is a special case , and exhibits much better decay because @xmath221 is purely an exponential function .
[ thm kpoissonftdecay ] suppose that @xmath222
. then for every @xmath223 , there exists a constant @xmath224 such that @xmath225    [ cor kpoissondecay ] suppose that @xmath222 .
then @xmath226 as @xmath215 for every @xmath223 .
we may refine the above estimates in the case @xmath227 to find a uniform bound on the @xmath228norm of @xmath229 .
[ thm l3 ] suppose that @xmath230\cup ( 0,\infty)\right)\setminus\mathbb{n } $ ] .
there exists a constant @xmath231 such that for all @xmath213 , @xmath232 .
so far , the upper bounds on @xmath3 may depend on the parameters @xmath24 and @xmath25 ; however , the following still holds .
[ lem h3 ] if @xmath233\cup(0,\infty))\setminus { \mathbb{n}}$ ] and @xmath222 , then @xmath234 for all @xmath235 .
we end the section with a statement on the zeros of @xmath9 .
[ thm h3 ] if @xmath233\cup(0,\infty))\setminus { \mathbb{n}}$ ] and @xmath222 , then @xmath236 for every @xmath237 .
in this section , we show that the results of riemenschneider and sivakumar @xcite have analogues for general multiquadrics . in section [ secpropfundamental ] , the decay of @xmath109 is discussed , and we use that information to uncover growth conditions on data that are suitable to cardinal interpolation . recall from corollary [ cor l3 ] and theorem [ thm h3 ] that for @xmath238 and @xmath222 , @xmath239 this decay is not the best one can get for individual @xmath24 , in fact for @xmath65 , buhmann @xcite proves a decay rate of @xmath22 . according to further work by buhmann and micchelli @xcite
, it appears that the multiquadrics with exponents @xmath240 for @xmath241 are exceptional cases . for these values
, @xmath109 can be shown to have more derivatives than what we have shown for general @xmath24 due to some special symmetry involving the bessel functions in the fourier transforms . moreover , in these cases , decay of the fundamental function is given by @xmath242    for negative exponents , the so - called inverse multiquadrics , the fundamental functions have slightly slower decay ( corollary [ cor l4 ] and theorem [ thm h3 ] ) : @xmath243    as a consequence of the decay of the fundamental functions , we have the following .    [ proprsprop33 ] for @xmath244 , the function @xmath245 is a well - defined , 1-periodic , bounded , continuous function .
periodicity is apparent , so we need only consider @xmath246 $ ] .
if @xmath247 , then gives @xmath248 which yields the result since @xmath3 is continuous on @xmath249 and @xmath250 .    on the other hand ,
if @xmath251 , then @xmath252 as the exponent is less than -1 , the series converges , whence the result .
we next define the cardinal interpolation operator acting on sequences , and show that it is well - defined for sequences that grow at a sufficiently slower rate than the decay of the fundamental function .
[ thmrsthm34 ] if @xmath253 , and @xmath254 for some fixed positive constants @xmath255 and @xmath62 , then the function @xmath256 is well - defined , continuous on @xmath249 , and satisfies @xmath257 @xmath215 .
if @xmath258 , and @xmath259 then @xmath7 as defined in is well - defined , continuous on @xmath249 , and satisfies @xmath260    consider the case @xmath253 .
we first show that @xmath261 is continuous on every interval of the form @xmath262 $ ] , for @xmath263 .
let @xmath264 $ ] .
then @xmath265 the first term on the right hand side is a finite sum of continuous functions , and so it suffices to show that the second sum converges uniformly on @xmath262 $ ] .
indeed , the decays of @xmath266 and @xmath3 yield the estimate @xmath267 which converges because @xmath268 and @xmath269 .    now to consider the decay of @xmath7 ,
let @xmath270 be fixed .
then let @xmath271 be the unique integer such that @xmath272 .
since @xmath273 , @xmath274 , and @xmath275 for every @xmath276 , we have @xmath277    the proof for negative @xmath24 values follows similar reasoning .
now we explore the properties of the cardinal interpolation operator @xmath7 acting on traditional sequence spaces ; we first show boundedness , and follow with estimates on its norm .    [ thmrsthm35 ] for a fixed @xmath244 and @xmath44 , the cardinal interpolation operator @xmath7 is a bounded linear operator from @xmath8 to @xmath278 for @xmath28 .
the proof is the same as in @xcite , which follows the techniques of @xcite .
linearity is evident , and the cases @xmath279 and @xmath280 follow from proposition [ proprsprop33 ] and lemma [ lem h3 ] , respectively .
therefore , let @xmath281 , and @xmath235 be fixed .
as before , let @xmath282 be the unique integer such that @xmath272 .
if @xmath283 , then theorem [ thmrsthm34 ] implies that @xmath261 is continuous , and we write @xmath284 to estimate the @xmath278 norm of the first term , notice that since @xmath234 , @xmath285    for the second term , first assume that @xmath253 . then @xmath286 the quantity in the final line above corresponds to the @xmath8 norm of the discrete convolution @xmath287 , where @xmath288 , and the entries of @xmath289 are given by @xmath290 . from (
* , theorem 7.6 ) , we find that @xmath291,\ ] ] where @xmath292 is a non - increasing rearrangement of @xmath289 .
it is easily checked that the supremum on the right hand side of is finite , and therefore the second term estimated above is @xmath293 , and the theorem follows .
the proof for the case @xmath251 is much the same , except that the power on the elements of @xmath289 will be @xmath294 , which nevertheless results in the supremum of @xmath295 being finite .
[ thmrsthm36 ] let @xmath244 and @xmath281 be fixed .
then @xmath296    consider a fixed @xmath281 , and @xmath297 such that @xmath298 .
let @xmath299 , @xmath300 , and @xmath282 be the unique integer such that @xmath272
. then @xmath301 by hlder s inequality and the fact that @xmath234 , @xmath302 to estimate @xmath303 , we represent @xmath3 by its fourier integral ( see proposition [ proplfundamental ] ) and integrate by parts .
indeed , @xmath304 where the final step follows from fubini s theorem .    from ( * ? ? ?
* proposition 1.3 ) , the _ mixed hilbert transform _ defined by @xmath305 is a bounded linear operator from @xmath8 to @xmath278 , and @xmath306 , where @xmath307 is a constant depending only on @xmath308 .
consequently , @xmath309 the final inequality comes from boundedness of the mixed hilbert transform and theorem [ thm l3 ] .
the conclusion of the theorem follows from the estimates on @xmath310 and @xmath303 .
we now estimate the norms in the cases @xmath280 and @xmath279 .
[ proprs38 ]    suppose @xmath244 and @xmath222 are fixed .
the following hold :    \(i ) let @xmath311 be defined via .
then @xmath312 where @xmath313 is a constant .
\(ii ) @xmath314    \(iii ) @xmath315    we supply the proof for the case @xmath316 , the estimates for the negative range of @xmath24 being wholly similar .
\(i ) by theorems [ thm l1 ] and [ thm l3 ] , @xmath317 since @xmath311 is 1-periodic , it suffices to check the inequality for @xmath246 $ ] .
let @xmath222 and @xmath318 .
then by @xmath319\\ & \leq   a_\alpha\left[1+\sum_{1\leq|j|\leq n-1}|j|^{-1}+\sum_{|j|\geq n}c^{5\alpha+2}|j|^{-2}\right]\\ & \leq   a_\alpha\left[1+\ln(n)+c^{5\alpha+2}n^{-1}\right],\\\end{aligned}\ ] ] whence the result .
\(ii ) simply note that if @xmath320 , then @xmath321 , @xmath235 , and apply ( i ) .
\(iii ) similarly , if @xmath322 , then @xmath323 whereby ( i ) provides the desired bound .
we make note that the case @xmath324 provides an interesting special case .
[ thmrsthm311 ] for any @xmath244 and @xmath222 , @xmath325    first , note that by plancherel s identity and a standard periodization argument , @xmath326 consequently , @xmath327}\max { \underset{k\in { \mathbb{z } } } { \displaystyle\sum}}\left(\widehat{l_{\alpha , c}}(\xi+2\pi k)\right)^2.\ ] ] taking the supremum over @xmath328 in the unit ball of @xmath329 yields @xmath330}\max { \underset{k\in { \mathbb{z } } } { \displaystyle\sum}}\left(\widehat{l_{\alpha , c}}(\xi+2\pi k)\right)^2.\ ] ] that the maximum on the right hand side is at most 1 is the content of , but the maximum is attained at @xmath331 by theorem [ thm h3 ] .
having established some information about the interpolation operators and their norms for different values of @xmath308 , we now shift our gaze to some convergence properties when the shape parameter @xmath6 for a fixed @xmath24 .
as one might expect , for large ( in absolute value ) @xmath24 , we obtain better convergence results owing to the more rapid decay of the fundamental functions .    to begin our discussion , consider the _
whittaker operator _ defined via @xmath332    this operator is bounded from @xmath8 to @xmath278 for @xmath281 , and the following holds .
[ thmconvergencetowhittaker ] if @xmath333 , @xmath281 , then for a fixed @xmath244 , @xmath334    boundedness of @xmath335 , theorem [ thmrsthm36 ] , and the uniform boundedness principle imply that it is sufficient to check convergence on the coordinate basis for @xmath8 , @xmath336 .
that is , if suffices that @xmath337 as @xmath338 by theorem [ thmbandlimitedconvergence ] , @xmath339 converges to 0 uniformly as @xmath6 ; moreover , proposition [ thm l3 ] implies that @xmath340 , whence an application of the dominated convergence theorem yields the statement of the theorem .
suppose a function @xmath10 has sufficiently slow growth so that @xmath261 , with @xmath341 , is well - defined ( see theorem [ thmrsthm34 ] ) .
in this case we let @xmath342 , and call this the cardinal interpolant of @xmath10 due to the identity @xmath343    we will consider pointwise and uniform convergence of @xmath344 to @xmath10 as @xmath6 , but first we must make some preliminary arrangements . consider @xmath24 to be fixed .
then define @xmath345 and @xmath346    [ lemrslem42 ] if @xmath253 and @xmath347 , then @xmath348 , is well - defined , continuous in each of its variables , and @xmath96-periodic in the second variable .    if @xmath258 and @xmath349 , then @xmath348 is well - defined , continuous in each of its variables , and @xmath96-periodic in the second variable .
continuity of @xmath350 is provided by corollaries [ cor l3 ] and [ cor l4 ] , while proposition [ prop l2 ] shows that the series @xmath351 is uniformly convergent on @xmath114 $ ] .
thus @xmath352 is well - defined and uniformly continuous for @xmath353 $ ] , and moreover , we may differentiate the series term by term .
finally , @xmath96-periodicity is evident .
now let @xmath354 be the space of continuous , @xmath96-periodic functions on the real line , and let @xmath355 denote its dual , which is the space of all @xmath96-periodic complex borel measures on the real line , with the total variation norm given by @xmath356 .
following @xcite , given @xmath357 , define the @xmath358-th fourier - stieltjes coefficient of @xmath359 by @xmath360 and the @xmath361-th fejr mean of the fourier - stieltjes series of @xmath359 by @xmath362    we begin by showing that for certain functions , the interpolant exhibits a special form .
[ thmrsthm43 ] let @xmath24 be fixed , and suppose @xmath10 is given by @xmath363 for some @xmath357 and some @xmath364 in the case @xmath253 or @xmath365 in the case @xmath258 .
let @xmath348 be defined by .
then @xmath366    by definition , @xmath367 , thus theorem [ thmrsthm34 ] implies that @xmath344 is well - defined and continuous . by property of
the fejr means , @xmath368 therefore , by the inversion formula and a standard periodization argument , @xmath369 { \widehat{l_{\alpha , c}}}(t)e^{-ixt}dt\nonumber\\ & = \dfrac{1}{2\pi } { \displaystyle\int}_{-\pi}^{\pi}\left [ { \underset{j =- n}{\overset{n}\sum}}\left(1-\frac{|j|}{n+1}\right)f(j)e^{ijt}\right]\phi_{\alpha , c}(x , t)dt.\end{aligned}\ ] ]    by definition , @xmath370 for every @xmath371 .
so if @xmath372 , @xmath373 , and if @xmath374 , @xmath375 .
consequently , if @xmath372 , then    @xmath376\phi_{\alpha , c}(x , t)dt = { \dfrac{1}{2\pi } } { \displaystyle\int}_{-\pi}^\pi\sigma_n(\mu , t)\phi_{\alpha , c}(x , t)dt.\ ] ]    if @xmath374 , then we integrate by parts @xmath377 times .
the boundary terms cancel due to periodicity , so @xmath378\phi_{\alpha , c}(x , t)dt \nonumber\\   & = ( -1)^k { \dfrac{1}{2\pi } } { \displaystyle\int}_{-\pi}^\pi\left [ { \displaystyle\sum}_{1\leq|j|\leq n}\left(1-\dfrac{|j|}{n+1}\right)\dfrac{f(j)}{(ij)^k}e^{ijt}\right]\phi^{(k)}_{\alpha , c}(x , t)dt \nonumber\\   & = { \dfrac{1}{2\pi } } { \displaystyle\int}_{-\pi}^\pi\left [ { \displaystyle\sum}_{1\leq|j|\leq n}\left(1-\dfrac{|j|}{n+1}\right)\widehat{\mu}(j)e^{ijt}\right]\phi_{\alpha , c}^{(k)}(x , t)dt\nonumber\\ & = { \dfrac{1}{2\pi } } { \displaystyle\int}_{-\pi}^\pi\left [ { \underset{j =- n}{\overset{n}\sum}}\left(1-\dfrac{|j|}{n+1}\right)\widehat{\mu}(j)e^{ijt}\right]\phi_{\alpha , c}^{(k)}(x , t)dt ,   \end{aligned}\ ] ] where the final equality comes from the fact that @xmath379 for @xmath380 .    combining , , and , we see that @xmath381    using theorem [ thmrsthm43 ] , we show the following result on uniform convergence .
[ thmrsthm44 ] suppose @xmath24 is as above , and @xmath10 is given by @xmath382 for some @xmath377 as in theorem [ thmrsthm43 ] and @xmath357 such that @xmath383,\ ] ] for some fixed @xmath384
. then @xmath385 with convergence being uniform on compact subsets of @xmath249 . in the case
@xmath372 , convergence is uniform on @xmath249 .
suppose first that @xmath380 .
let @xmath79 be fixed .
we will show that @xmath386 uniformly for @xmath387 . by theorem [ thmrsthm43 ] and the definitions of @xmath10 and @xmath348 , @xmath388    by applying the leibniz rule ,
@xmath389 the first term above converges to 0 uniformly for @xmath390 , which can be seen from the last portion of the proof of proposition [ propconvergencetosinc ] .
additionally , by proposition [ prop l2](i ) the second term above is bounded above by @xmath391 and hence converges uniformly to 0 as @xmath6 .
similarly , write @xmath392.\ ] ]    here , if @xmath390 , then whenever @xmath393 and @xmath394 , @xmath395 falls into the respective ranges for proposition [ prop l2](ii ) and ( iii ) .
applying the estimates of that proposition demonstrates that @xmath396 uniformly as @xmath6 for @xmath387 .    if @xmath372 , then we again split the integral into two pieces , and analyze the corresponding integrands @xmath397 and @xmath398
however , notice that there are no terms of the form @xmath399 , and so the inequalities from proposition [ prop l2 ] give upper bounds on the integrands that do not depend on @xmath165 at all , and so the convergence of @xmath7 to @xmath10 is uniform on @xmath249 .    the condition on the support of the measure @xmath359 in the previous theorem
was essential because of the use of proposition [ prop l2 ] in the proof .
heuristically , the condition on the support should be of no surprise due to the fact that we have no uniform control ( in @xmath25 ) of the derivatives of @xmath109 at the boundary points @xmath400 .
in fact , it is likely that the derivatives get much larger near the boundary as @xmath25 grows , especially given the fact that @xmath109 converges to the characteristic function of @xmath401 .
nevertheless , we may make a weaker assumption on @xmath359 which still yields a convergence result .
[ thmrsthm46 ] suppose @xmath24 is as above , and @xmath10 is given by @xmath402 for some @xmath357 which is absolutely continuous with respect to the lebesgue measure
. then @xmath385 with convergence being uniform on @xmath249 .
the proof is quite similar to that of theorem [ thmrsthm44 ] .
let @xmath403 be arbitrary , and choose @xmath404 such that @xmath405+|\mu|[\pi(1- {   \varepsilon}),\pi]<\gamma$ ] since @xmath406 is absolutely continuous with respect to the lebesgue measure . then , as before , @xmath407    split each integral into one over @xmath408 $ ] and one near the endpoints .
the integral over the interior segment can be handled exactly as in the proof of theorem [ thmrsthm44 ] .
for the integrals near the endpoints , notice that @xmath409 , and so the first integral is at most @xmath410+|\mu|[\pi(1- {   \varepsilon}),\pi])$ ] , which is at most @xmath411 by the choice of @xmath412 .
meanwhile , by proposition [ prop l2](iii ) , the second integrand is at most @xmath413 which can be bounded by a constant depending only on @xmath24 .
thus @xmath414 for some constant @xmath415 independent of @xmath25 , and so the conclusion follows .
we conclude this section with the special consideration of the case @xmath220 , where @xmath4 is called the poisson kernel .
as mentioned above , the fundamental function for the poisson kernel , @xmath416 , decays faster than any polynomial .
thus any results in this section that depend on the decay rate of the fundamental function hold true for the poisson kernel . in particular , propositions [ proprsprop33 ] and [ proprs38 ] , and theorems [ thmrsthm35 ] , [ thmrsthm36 ] , [ thmrsthm311 ] ,
[ thmconvergencetowhittaker ] , and [ thmrsthm46 ] are all valid for the poisson kernel , since existence of the interpolant primarily depends on the decay of the fundamental function .    moreover , because the fundamental function for the poisson kernel decays so rapidly , we find much stronger versions of the remaining results in this section .
we begin by stating the analogue of theorem [ thmrsthm34 ] .
[ thmrsthm34poisson ] suppose that @xmath417 for some @xmath223 and constant @xmath62 .
then the _ cardinal poisson interpolant _ ,
@xmath418 is well - defined , continuous on @xmath249 , and satisfies @xmath419 for any @xmath420 .    as a consequence of the preceding theorem , the function @xmath421 defined by has well - defined derivatives of all orders .
in particular , lemma [ lemrslem42 ] holds , and consequently we find the following analogue of theorem [ thmrsthm43 ] .
suppose @xmath10 is given by @xmath422 for some @xmath357 and some @xmath223 .
then @xmath423    finally , theorem [ thmrsthm44 ] holds for any @xmath223 for the poisson kernel .
in this section , we illustrate the convergence phenomena discussed in the previous section .
the examples are of a similar flavor to those found in @xcite .
[ example1 ] let @xmath253 and @xmath347 , or @xmath258 and @xmath424 , or @xmath220 and @xmath223 .
let @xmath425 be the @xmath96-periodic extension of the measure @xmath426 , where @xmath427 is the usual dirac measure at 0 .
if @xmath428 then theorem [ thmrsthm44 ] implies that @xmath429 uniformly on compact subsets of @xmath249 .
however , theorem [ thm h3 ] allows us to say more given this information .
if @xmath372 , then we find the following identity on account of theorem [ thmrsthm43 ] :    @xmath430    for higher order polynomials , we may use theorems [ thmrsthm43 ] and [ thmrsthm44 ] to show that @xmath431 whereby one can obtain an error bound in terms of @xmath25 via proposition [ prop l2 ] .
this also demonstrates that @xmath432 uniformly on compact subsets of @xmath249 .
let @xmath433 be fixed , and let @xmath359 be the @xmath96-periodic extension of @xmath434}dt$ ] , where @xmath435}$ ] takes value 1 on @xmath436 $ ] and 0 elsewhere .
define @xmath437 since @xmath359 is absolutely continuous with respect to the lebesgue measure , theorem [ thmrsthm46 ] implies that @xmath438 uniformly on @xmath249 .
note that this fact also follows from theorem [ thmbandlimitedconvergence ] , albeit from substantially different reasoning .
however , the hypothesis of theorem [ thmbandlimitedconvergence ] requires the function being interpolated to be bandlimited , and so the following result can not be obtained simply by appealing to that theorem .
let @xmath24 and @xmath377 be as in example [ example1 ] .
let @xmath425 be the periodic extension of @xmath439}dt$ ] for some @xmath440 .
then if @xmath441 theorem [ thmrsthm44 ] implies that @xmath442 , uniformly on compact subsets of @xmath249 .
in this section we prove the various results listed in section [ secpropfundamental ] .
our methods closely resemble those found in @xcite . to reduce the clutter in our calculations , we will henceforth drop explicit dependence upon @xmath24 in our calculations that follow .
we begin by rewriting in terms of a laplace transform to exhibit the singularity at the origin . to do this , we require an integral representation for the bessel function , which we find from @xcite : @xmath443        where @xmath415 is a constant .
we relabel the product of the exponential and the integral in the above expression @xmath447 .
that is , @xmath448 hence may be abbreviated as @xmath449 for @xmath450 , we have @xmath451    we turn our attention to the estimates involving @xmath447 and its derivatives , and begin by noting that @xmath452(c),\ ] ] where @xmath453 denotes the usual laplace transform , @xmath454(s)=\int_{0}^{\infty}f(t)e^{-st}dt$ ] . our estimates for @xmath447 will rely on estimates for the laplace transform . the first result in this direction
is a lower bound .
we obtain this bound by combining the inequalities @xmath457(c ) \geq \int_{0}^{\infty } e^{-ct } ( 2|\xi|)^\alpha t^{\alpha } dt = \dfrac{(2|\xi|)^\alpha \gamma(\alpha+1)}{c^{\alpha+1}}\ ] ] and @xmath457(c ) \geq \int_{0}^{\infty } e^{-ct }   t^{2\alpha } dt = \dfrac{\gamma(2\alpha+1)}{c^{2\alpha+1}}.\ ] ]        1 .   if @xmath461 , @xmath462(c)\right|\leq \dfrac{(4|\xi|)^{\alpha - l}\gamma(\alpha+1)}{c^{\alpha+1 } } + \dfrac{2^{\alpha - l}\gamma(2\alpha+1-l)}{c^{2\alpha+1-l}};\ ] ] 2 .   if @xmath463,@xmath464(c)\right|\leq \dfrac{2^{3\alpha-2l}|\xi|^{2\alpha - l}}{c } + \dfrac{2^{\alpha - l}\gamma(2\alpha+1-l)}{c^{2\alpha+1-l}};\ ] ] 3 .   if @xmath465 ( thus @xmath24 is an odd half integer ) , @xmath464(c)\right|\leq \ln\left ( 1+|\xi|^{-1 } \right ) + \dfrac{e^{-c}}{c(1+|\xi|)}.\ ] ]    for ( i ) , we have the following @xmath466(c )   = & \int_{0}^{2|\xi|}e^{-ct } t^\alpha(t+2|\xi|)^{\alpha - l }   dt    +   \int_{2|\xi|}^{\infty}e^{-ct } t^\alpha ( t+2|\xi|)^{\alpha - l }    dt \\ \leq &    \int_{0}^{2|\xi|}e^{-ct } t^\alpha(4|\xi|)^{\alpha - l }   dt    +   \int_{2|\xi|}^{\infty}e^{-ct } t^\alpha(2t)^{\alpha - l }    dt \\ \leq & \dfrac { ( 4|\xi|)^{\alpha - l}\gamma(\alpha+1)}{c^{\alpha+1 } }   + \dfrac{2^{\alpha - l}\gamma(2\alpha+1-l)}{c^{2\alpha+1-l}}.\end{aligned}\ ] ] inequality ( ii ) follows from a similar calculation , replacing @xmath467 with @xmath468 in the first integral below : @xmath466(c )   = & \int_{0}^{2|\xi| } e^{-ct}t^\alpha(t+2|\xi|)^{\alpha - l }   dt    +   \int_{2|\xi|}^{\infty}e^{-ct } t^\alpha ( t+2|\xi|)^{\alpha - l }    dt \\ \leq &    \int_{0}^{2|\xi|}e^{-ct}2^{3\alpha -2l}|\xi|^{2\alpha -l }   dt    +   \int_{2|\xi|}^{\infty } e^{-ct}2^{\alpha - l}t^{2\alpha - l }    dt \\ \leq & \dfrac{2^{3\alpha -2l}|\xi|^{2\alpha -l }   } { c }   +
\dfrac{2^{\alpha - l } \gamma(2\alpha+1-l ) }   { c^{2\alpha+1-l}}.\end{aligned}\ ] ] finally , for inequality ( iii ) , we have @xmath469(c ) = & \int_{0}^{\infty}e^{-ct}\dfrac{t^\alpha}{(t+2|\xi|)^{\alpha+1}}dt = \int_{0}^{\infty}e^{-c|\xi|t}\dfrac{t^\alpha}{(t+2)^{\alpha+1 } }   dt\\ = & \int_{0}^{|\xi|^{-1}}e^{-c|\xi|t}\dfrac{t^\alpha}{(t+2)^{\alpha+1}}dt     + \int_{|\xi|^{-1}}^{\infty}e^{-c|\xi|t}\dfrac{t^\alpha}{(t+2)^{\alpha+1}}dt \\ \leq &   \int_{0}^{|\xi|^{-1 } }   ( 1+t)^{-1}dt    +    \int_{|\xi|^{-1}}^{\infty } e^{-c|\xi|t } ( 1+t)^{-1}dt\\ \leq & \ln\left(1+|\xi|^{-1}\right ) + \dfrac{e^{-c}}{c(1+|\xi| ) } \leq \ln\left(1+|\xi|^{-1}\right ) +   \dfrac{e^{-c}}{c}.\end{aligned}\ ] ]        applying the leibniz rule to , we see that @xmath473(c ) \right)(\xi).\end{aligned}\ ] ] differentiating under the integral sign , which is justified by the exponential decay of the integrand , and using the triangle inequality yield .              to see ( i ) ,
apply the leibniz rule , triangle inequality , and lemma [ lem jl3 ] to . as a prelude to ( ii ) , we remark that combining ( i ) with lemma [ lem jl2 ] ( i ) and ( ii ) reveals that if @xmath485 , latexmath:[\[\label{eq jl8 }    if @xmath476 , we pick up an extra singularity from lemma [ lem jl2 ] ( iii ) , in which case we split the double sum in ( i ) into three parts , the first of which corresponds to @xmath487 , the second to @xmath488 , and the third to @xmath489 : @xmath490(c)\\ \nonumber = &   a_\alpha c^{2\alpha+1}e^{-c|\xi|}|\xi|^{-4\alpha-2 } \left ( |\xi|^{2\alpha+1}\ln(1+|\xi|^{-1 } ) \right.\\ \nonumber + & \left .
\sum_{l'=0}^{k-1}a_{k , l'}c^{2\alpha+1-l'}|\xi|^{2\alpha+1}\mathcal{l}[t^\alpha(t+2|\xi|)^{\alpha - l'}](c)\right .
\\ \nonumber + & \left.\sum_{l=0}^{k-1}\sum_{l'=0}^{l}a_{k , l , l'}c^{l - l'}|\xi|^{l}\mathcal{l}[t^\alpha(t+2|\xi|)^{\alpha - l'}](c ) \right)\\ \label{eq jl9}\leq & a_\alpha c^{2\alpha+1}e^{-c|\xi|}|\xi|^{-2\alpha-1}\ln(1+|\xi|^{-1})\nonumber\\ + & a_\alpha c^{4\alpha+1-\lfloor\alpha\rfloor}e^{-c|\xi|}|\xi|^{-4\alpha-2}[1+o(1 ) ] , \quad |\xi|\to 0.\end{aligned}\ ] ] to track the largest powers of @xmath25 , we note that there is nothing to do in the logarithmic term ; in the second term we use lemma [ lem jl2 ] and take @xmath491 to obtain @xmath492 ; for the third term we let @xmath493 and @xmath491 to obtain @xmath494 .    to prove ( ii ) , we again apply the leibniz rule to obtain @xmath495 where @xmath496 is the set of increasing non - negative integer partitions of @xmath377 , that is , @xmath497 . for @xmath485 , we may plug , , and into to obtain : @xmath498 , \quad|\xi|\to 0.\end{aligned}\ ] ]                recall that for @xmath503 , @xmath504 . from the leibniz rule
, we have @xmath505 thus we may use and to obtain @xmath506 where @xmath507 $ ] . since @xmath198 , @xmath508 , and @xmath509 , we have @xmath510 this is the desired estimate for @xmath511 . if @xmath476 , we must use and to handle the logarithmic term .
this changes @xmath512 ; however , it is still true that @xmath513 as @xmath202 .
an immediate consequence of proposition [ prop l1 ] is that the function @xmath514 defined in converges for @xmath198 and we can differentiate the series term by term .
in fact , by applying the estimates from , we have latexmath:[\[\label{eq jl14 }            analogous to the argument given in @xcite , we find that if @xmath189 and @xmath521 $ ] , then @xmath522 where @xmath523 .
this means that @xmath524 $ ] , so we may use the leibniz rule together with proposition [ prop l1 ] and lemma [ lem jl5 ] ( letting @xmath525 ) to obtain our next result .              in the special case of @xmath220
, we get the poisson kernel , whose fourier transform is given by @xmath531 using this much simpler formula allow us to simplify our earlier work and get stronger results .
in fact , we see that the analogue of proposition 4.1 is given by @xmath532 where there is no longer a restriction on @xmath197 .
the proof of theorem [ thm kpoissonftdecay ] now follows the same line of reasoning as that used to prove theorem [ thm l1 ] .
we begin by noting that @xmath533 is even , so we need only consider the integral @xmath534 for @xmath32 , we use lemma [ lem jl5 ] with @xmath535 and find that @xmath536 .
the quantity @xmath537 may be estimated similarly using lemmas [ lem jl6 ] and [ lem jl7 ] ( and again letting @xmath538 ) .
we have @xmath539 the last inequality follows from the fact that both terms in parentheses have a global maximum which depends only on @xmath24 .
finally , we show that @xmath533 is monotone on @xmath540 $ ] , which implies that @xmath541 . to that end
, we use the quotient rule to write @xmath542 we recall that @xmath543 and by the formula found in @xcite , @xmath544 . these allow us to rewrite the numerator , @xmath545 where @xmath546 from these expressions , we see that @xmath547 for @xmath548 $ ] and @xmath549 on @xmath540 $ ] .
the latter is true since @xmath550 is decreasing .
this shows that @xmath533 is decreasing on @xmath540 $ ] as desired .
recall that @xmath533 is non - negative as commented in section [ secbaxteranalogue ] , and theorems [ thm l1 ] and [ thm l2 ] show that @xmath551 is continuous and integrable for the given range of @xmath24 .
consequently , @xmath551 is positive definite , and thus for every @xmath235 , @xmath552 the final equality coming from .
we first consider the case @xmath372 .
as in section 3 , write @xmath553 .
proposition [ prop l1 ] with @xmath554 implies that the series converges uniformly for @xmath555 .
consequently , @xmath556 . to see that this limit is 1 requires estimating @xmath557 in a slightly different manner than we have so far . combining and
, we estimate @xmath558(c)}{\mathcal{l}[t^\alpha(t+2|\xi|)^\alpha](c)}\right|,\ ] ] which by lemma [ lem jl1 ] and proposition [ lem jl2](i ) is at most @xmath559 the series is summable , and uniformly bounded for @xmath555 , and it follows that @xmath560 .    now for @xmath561 , let @xmath562 with @xmath555
. then we may write @xmath563}\\   & = \dfrac{|\xi|^{2\alpha+1}\widehat{\phi_c}(\xi+2\pi k)}{a_\alpha c^{2\alpha+1}f_\alpha(\xi)\left[1+\dfrac{c^{-2\alpha-1}|\xi|^{2\alpha+1}}{f_\alpha(\xi ) } { \underset{j\neq0 } { \displaystyle\sum}}\widehat{\phi_c}(\xi+2\pi j)\right]},\end{aligned}\ ] ] the second equality following from . by definition , @xmath564 , and a similar argument to the case @xmath372 above shows that we may allow @xmath81 to tend to 0 , and conclude that @xmath565 .
s. d. riemenschneider and n. sivakumar , on cardinal interpolation by gaussian radial - basis functions : properties of fundamental functions and estimates for lebesgue constants , _ j. anal .
_ * 79 * ( 1999 ) , 33 - 61 .