Linear Algebra and its Applications

Similar documents
Linear Algebra and its Applications

So reconstruction requires inverting the frame operator which is often difficult or impossible in practice. It follows that for all ϕ H we have


University of Missouri Columbia, MO USA

arxiv:math/ v1 [math.fa] 14 Sep 2003

Applied and Computational Harmonic Analysis

A NEW IDENTITY FOR PARSEVAL FRAMES

Real, Tight Frames with Maximal Robustness to Erasures

Preconditioning of Frames

2 PETER G. CASAZZA, MANUEL T. LEON In this paper we generalize these results to the case where I is replaced by any positive selfadjoint invertible op

THE KADISON-SINGER PROBLEM AND THE UNCERTAINTY PRINCIPLE

BANACH FRAMES GENERATED BY COMPACT OPERATORS ASSOCIATED WITH A BOUNDARY VALUE PROBLEM

j jf, S K cf = j K c j jf, f H.

A DECOMPOSITION THEOREM FOR FRAMES AND THE FEICHTINGER CONJECTURE

Approximating scalable frames

G-frames in Hilbert Modules Over Pro-C*-algebras

Minimizing Fusion Frame Potential

FRAME DUALITY PROPERTIES FOR PROJECTIVE UNITARY REPRESENTATIONS

A FRAME THEORY PRIMER FOR THE KADISON-SINGER PROBLEM

A primer on the theory of frames

C -Algebra B H (I) Consisting of Bessel Sequences in a Hilbert Space

Frame Diagonalization of Matrices

g-frame Sequence Operators, cg-riesz Bases and Sum of cg-frames

CONSTRUCTIVE PROOF OF THE CARPENTER S THEOREM

Decompositions of frames and a new frame identity

Dual and Similar Frames in Krein Spaces

Real Equiangular Frames

Constructive Proof of the Carpenter s Theorem

Duals of g-frames and g-frame Sequences

On the Feichtinger conjecture

LOCAL AND GLOBAL STABILITY OF FUSION FRAMES

Operators with Closed Range, Pseudo-Inverses, and Perturbation of Frames for a Subspace

Density results for frames of exponentials

Journal of Mathematical Analysis and Applications. Properties of oblique dual frames in shift-invariant systems

Shift Invariant Spaces and Shift Generated Dual Frames for Local Fields

MULTIPLEXING AND DEMULTIPLEXING FRAME PAIRS

Spanning and Independence Properties of Finite Frames

On the Equality of Fusion Frames 1

Subsequences of frames

Yimin Wei a,b,,1, Xiezhang Li c,2, Fanbin Bu d, Fuzhen Zhang e. Abstract

Optimal dual fusion frames for probabilistic erasures

Affine and Quasi-Affine Frames on Positive Half Line

Ole Christensen 3. October 20, Abstract. We point out some connections between the existing theories for

On Riesz-Fischer sequences and lower frame bounds

The Drazin inverses of products and differences of orthogonal projections

PROJECTIONS AND THE KADISON-SINGER PROBLEM

Spectrally Uniform Frames And Spectrally Optimal Dual Frames

Fast Angular Synchronization for Phase Retrieval via Incomplete Information

Scalable Frames and Convex Geometry

Applied Mathematics Letters

Introduction to Hilbert Space Frames

Parseval Frame Construction

DS-GA 1002 Lecture notes 0 Fall Linear Algebra. These notes provide a review of basic concepts in linear algebra.

Means of unitaries, conjugations, and the Friedrichs operator

Linear Algebra and its Applications

Decomposition of Riesz frames and wavelets into a finite union of linearly independent sets

Intrinsic products and factorizations of matrices

Robustness of Fusion Frames under Erasures of Subspaces and of Local Frame Vectors

Constructing tight fusion frames

Frame Wavelet Sets in R d

A CONSTRUCTIVE APPROACH TO THE FINITE WAVELET FRAMES OVER PRIME FIELDS

Lecture 7: Positive Semidefinite Matrices

Uniqueness of the Solutions of Some Completion Problems

Linear Algebra and its Applications

1 Last time: inverses

Loosen Up! An introduction to frames.

RANK AND PERIMETER PRESERVER OF RANK-1 MATRICES OVER MAX ALGEBRA

PRECONDITIONING TECHNIQUES IN FRAME THEORY AND PROBABILISTIC FRAMES

Frames and operator representations of frames

Linear Algebra and its Applications

IN AN ALGEBRA OF OPERATORS

arxiv: v2 [math.fa] 11 Nov 2016

DUALITY PRINCIPLE IN g-frames

Fourier and Wavelet Signal Processing

Semi-orthogonal wavelet frames on positive half-line using the Walsh Fourier transform

On Sum and Restriction of Hypo-EP Operators

SPARSE NEAR-EQUIANGULAR TIGHT FRAMES WITH APPLICATIONS IN FULL DUPLEX WIRELESS COMMUNICATION

The minimum rank of matrices and the equivalence class graph

1. General Vector Spaces

Invariances of Frame Sequences under Perturbations

PRECONDITIONING TECHNIQUES IN FRAME THEORY AND PROBABILISTIC FRAMES

Preconditioning techniques in frame theory and probabilistic frames

MATH 20F: LINEAR ALGEBRA LECTURE B00 (T. KEMP)

Linear Algebra and its Applications

INVARIANCE OF A SHIFT-INVARIANT SPACE

Math 3191 Applied Linear Algebra

Where is matrix multiplication locally open?

IMA Preprint Series # 2066

MATH 2331 Linear Algebra. Section 2.1 Matrix Operations. Definition: A : m n, B : n p. Example: Compute AB, if possible.

PROBLEMS ON PAVING AND THE KADISON-SINGER PROBLEM. 1. Notation

Math 407: Linear Optimization

Boolean Inner-Product Spaces and Boolean Matrices

Approximately dual frames in Hilbert spaces and applications to Gabor frames

Houston Journal of Mathematics. University of Houston Volume, No.,

PERTURBATION OF FRAMES FOR A SUBSPACE OF A HILBERT SPACE

Spectrally arbitrary star sign patterns

The Residual Spectrum and the Continuous Spectrum of Upper Triangular Operator Matrices

Linear Algebra and its Applications

On Euclidean distance matrices

arxiv:math/ v1 [math.fa] 5 Aug 2005

Transcription:

Linear Algebra and its Applications 439 (2013) 1330 1339 Contents lists available at SciVerse ScienceDirect Linear Algebra and its Applications journal homepage: www.elsevier.com/locate/laa Maximum robustness and surgery of frames in finite dimensions < Martin S. Copenhaver a, Yeon Hyang Kim b,cortneylogan c, Kyanne Mayfield d, Sivaram K. Narayan b,, Jonathan Sheperd e a School of Mathematics, Georgia Institute of Technology, United States b Department of Mathematics, Central Michigan University, United States c Department of Mathematics, Stonehill College, United States d Department of Mathematics, University of Portland, United States e Department of Mathematics, University of Notre Dame, United States ARTICLE INFO Article history: Received 7 February 2013 Accepted 10 April 2013 Availableonline13May2013 Submitted by R.A. Brualdi AMS classification: Primary 42C15 05B20 15A03 ABSTRACT We consider frames in a finite-dimensional Hilbert space H n where frames are exactly the spanning sets of the vector space. We present a method to determine the maximum robustness of a frame. We present results on tight subframes and surgery of frames. We also answer the question of when length surgery resulting in a tight frame set for H n is possible. 2013 Elsevier Inc. All rights reserved. Keywords: Frames Tight frames Diagram vectors Robustness Erasures Redundancy Gramian operator Subframes Surgery Length surgery < Research supported by NSF-REU Grant DMS 08-51321. This work was done as a part of the REU program in Summer 2011. Corresponding author. E-mail addresses: copenhaver@gatech.edu (M.S. Copenhaver), kim4y@cmich.edu (Y.H. Kim), logan.cort@gmail.com (C. Logan), mayfield13@up.edu (K. Mayfield), sivaram.narayan@cmich.edu (S.K. Narayan), jsheperd@nd.edu (J. Sheperd). 0024-3795/$ - see front matter 2013 Elsevier Inc. All rights reserved. http://dx.doi.org/10.1016/j.laa.2013.04.016

M.S. Copenhaver et al. / Linear Algebra and its Applications 439 (2013) 1330 1339 1331 1. Introduction Abasis{f i } k i=1 of vectors in a finite-dimensional inner product space H n can be used to represent every vector f as a linear combination of the elements in {f i } k i=1 : f = i a i f i. This representation gives us characteristics of f in terms of the coefficients a i. However, uniqueness of this representation is not always an advantage. For example, in applications such as image and signal processing, the loss of a single coefficient during data transmission will prevent the recovery of the original signal. A frame is a generalization of a basis that includes redundancy. That is, a frame in finite dimensions is a redundant set of vectors that span a finite-dimensional vector space. This redundancy yields robustness, which makes frame representations less sensitive to transmission errors. The study of frames began in 1952 with their introduction by Duffin and Schaeffer [10] and has since been expanded by Daubechies [8] and others [2,6,5,17]. In [14], the authors characterize frames in R n that are robust to k erasures and a necessary and sufficient condition is given for performing (r, k)-surgery on unit-norm tight frames in R 2.Inthis paper, we present a method to determine the maximum robustness of any given frame, and generalize results on surgery from [14]. We also answer the question of when length surgery resulting in a tight frame set for H n is possible. We begin by defining various notions that are mentioned above. A good introduction to frames in finite dimensions can be found in [9,15]. Let I be a subset of N and H n be an n-dimensional real or complex Hilbert space. A frame in H n is a sequence of vectors {f i } i I for which there exist constants 0 < A B < such that for all f H n, A f 2 f, f i 2 B f 2. i I When A = B = λ, {f i } i I is called a λ tight frame. Whenλ = 1, the frame is called a Parseval frame. A unit-norm frame is a frame such that each vector in the frame has norm one. In a finite-dimensional Hilbert space H n, a sequence of vectors is a frame if and only if it spans H n. Given a sequence of vectors {f i } k i=1 in H n,wedefinetheanalysis operator to be the linear map θ : H n l 2 ({1,, k}) defined by (θf )(i) = f, f i. The adjoint θ such that θ : l 2 ({1,, k}) H n is called the synthesis operator. InH n, the analysis operator associated with a sequence of vectors {f i } k i=1 can be written with respect to a basis as the k n matrix f 1 θ =., f k and the synthesis operator as the n k matrix θ = f 1 f k. The frame operator S of a sequence of vectors {f i } k i=1 (not necessarily a frame) is defined as θ θ.for all f H n, Sf = θ θf = k f, f i f i. i=1

1332 M.S. Copenhaver et al. / Linear Algebra and its Applications 439 (2013) 1330 1339 For {f i } k i=1 H n,thegramian operator G is the k k matrix defined by G = θθ = ( f j, f i ) k. i,j=1 Given a sequence of vectors {f i } k i=1 in H n, it is known that the frame operator S of the sequence has rank n if and only if the sequence is a frame. The frame operator S = λi n if and only if {f i } k i=1 is a λ tight frame. Moreover, S = I n if and only if {f i } k i=1 is a Parseval frame [9]. 2. Robustness In this section, we present a method to determine the maximum robustness of any given frame. Definition 2.1. AframeF := {f i } k i=1 for H n is said to be robust to r erasures if {f i } i I C is still a frame for any index set I {1, 2,...,k} with I =r. Themaximum robustness of a frame F for H n, which we denote rob(f), is defined to be the maximum number r such that F is robust to r erasures. Observation 2.2. From the definition it follows that rob(f) k n. In the case when F is a full spark frame, i.e., every set of n vectors in the frame F form a basis for H n [1], we have rob(f) = k n. The goal is to develop methods to easily compute the maximum robustness of any given frame. There are some known results that allow us to check whether a frame is robust to a particular number of erasures. One obvious method to determine maximum robustness is to use these results to find the number r such that the frame is robust to r, butnotr + 1, erasures. We first state two known results. Theorem 2.1 ([5]). Let F ={f i } k i=1 be a frame for Rn. The following are equivalent: (1) F is a frame robust to one erasure. (2) There are scalars c i = 0, for1 i k, such that k c i f i = 0. i=1 The following is a generalization of Theorem 2.1. Theorem 2.2 ([14]). Let F ={f i } k i=1 be a frame for Rn with synthesis operator θ. The following are equivalent: (1) F is a frame robust to r erasures. (2) For all index sets I {1, 2,...,k} with I = r 1, I C { supp(f ) : f null(θ ) }, where supp(f ) is the set of indices where the vector f has nonzero components. We note that condition (2) can be checked in practice for a given frame and given I C as follows: Let A I be the submatrix of θ consisting of the I C columns and let v 1,, v p be a spanning set of null(a I ) (for example, using SVD). Then p i=1 supp(v i) ={1,, I c } if and only if there exists v null(a I ) with supp(v) ={1,, I c }, which is equivalent to the existence f null(θ ) with supp(f ) = I C. This result leads directly to an exponential algorithm for determining robustness of a frame. For more details about erasure we refer the reader to [12,16,11]. The following lemma states that for finding maximum robustness of a frame we can ignore zero vectors in the frame.

M.S. Copenhaver et al. / Linear Algebra and its Applications 439 (2013) 1330 1339 1333 Lemma 2.3. Let F ={f i } k i=1 be a frame for H n.thenrob(f) = rob(f \{0}). The proof of Lemma 2.3 follows from the observation that for any subset S of frame vectors, span(s) = span(s \{0}). The next proposition gives an upper bound for rob(f). Proposition 2.4. Let F ={f i } k i=1 be a frame for H n and let S ={I {1, 2,...,k} :span({f i } i I ) = H n } be the collection of index sets that correspond to spanning sets. Then rob(f) log 2 S. Proof. Suppose F is robust to m = log 2 S + 1 erasures. Then {f i } k i=m+1 is a frame. Let S ={I {m + 1,...,k} :span({f i } i I ) = H n }. Then for any J {1, 2,, m} and I S, I J S.Thuswehave S 2 m S, which implies S S < S = 1. 2 m 2 log 2 S Consequently, S =. This is a contradiction since {m + 1,, k} S. We note that if F is a frame robust to r erasures, then r rob(f) and r n S i=0 n r + i since each set of at least n r vectors is a frame. In the next proposition, we express the maximum robustness of a frame using the cardinality of maximum nonspanning set. Proposition 2.5. Let F ={f i } k i=1 be a frame for H n,andletn ={I {1, 2,...,k} :span({f i } i I ) = H n } be the collection of index sets that correspond to nonspanning sets. Then rob(f) = k 1 max{ I :I N}. Proof. Let I 0 N be such that I 0 = max{ I : I N}. Suppose any k 1 I 0 vectors are removed from F. Then I 0 +1 vectors remain. Since the remaining set spans H n, F is robust to k 1 I 0 erasures. Suppose the k I 0 vectors whose indices are not in I 0 are removed from F.The remaining vectors corresponding to I 0 do not span H n,sof is not robust to k I 0 erasures. Therefore, rob(f) = k 1 I 0. Using Proposition 2.5 wecanfindavectory H n so that the maximum robustness of a frame F ={f i } k i=1 in H n is equal to one less than the count of nonzero numbers in the set { f i, y } k i=1.thisis stated in the next theorem. Theorem 2.3. If F ={f i } k i=1 is a frame for H n, then there exists a y H n such that rob(f) = rob({ f i, y } k i=1 ). Proof. Let N be a largest non-spanning subset of {f i } k i=1.thenwehaverob({f i} k ) = i=1 k 1 N by Proposition 2.5 and span (N) = span (y) for some y in H n.since f i, y = 0 if and only if f i N, { f i, y } k i=1 contains N zero vectors and k N nonzero vectors, which implies that rob({ f i, y } k ) = i=1 k 1 N. The following is an example of the implementation of this theorem. We note that in general finding the largest non-spanning set is a combinatorially hard problem.

1334 M.S. Copenhaver et al. / Linear Algebra and its Applications 439 (2013) 1330 1339 Example 2.6. Let F = 3, 2, 1, 1, 0, 1 0 0 0 1 1 1 R2. The largest non-spanning set is N = 3, 2, 1. As in the proof of Theorem 2.3, 0 0 0 we take y = 0. Thus, taking the inner product of vectors in F with y, wehave{0, 0, 0, 1, 1, 1}, 1 which clearly has a maximum robustness of two. Therefore F must have a maximum robustness of two as well. The following theorem tells us that a transformation onto a smaller space can also preserve the maximum robustness of a particular frame. Theorem 2.4. Let F ={f i } k i=1 be a frame for H n with rob(f) = r and let P be an orthogonal projection onto a subspace U. The frame PF for U has maximum robustness r if and only if there exists a subset F of k r 1 vectors which does not span H n and span(f ) U. Proof. ( ) Without loss of generality, we assume that {P(f i )} k r 1 i=1 does not span U. ThenF = {f i } k r 1 i=1 does not span H n.sincerob(f) = r, span(f ) has dimension one. Since PF does not span U, there exists a nonzero vector e U such that e, Pf = e, f = 0, for f F, which implies that span(f ) = span (e) U. ( ) Suppose there is some set F of k r 1vectorsfromF such that F does not span H n and span(f ) U. Then for any nonzero vector e span(f ) and f F,0= e, f = e, Pf since e U. Thus span(f ) has dimension at least one. Therefore P(F ) does not span U, which implies that P(F) is not robust to r + 1 erasures. Since the projected frame is necessarily robust to r erasures, itsmaximumrobustnessisr. Definition 2.7. [3] LetF = {f i } k i=1 be a frame for H n.foreachx S, theredundancy function R F : S R + is defined by R F (x) = k P fi (x) 2, i=1 where S = {x H n : x = 1} is the unit sphere in H n and P fi (x) is the orthogonal projection of x onto the span of f i. In [3], the concepts of upper and lower redundancy are also defined. The upper redundancy R + F of a frame is the maximum of its redundancy function taken over the unit sphere S.Thelower redundancy R F of a frame is the minimum of its redundancy function taken over the unit sphere S. According to [3], R F is the maximum number of disjoint spanning sets in the frame F and R + F is the minimum number of disjoint linearly independent sets in the frame F. If a frame F in R n has a maximum of R F disjoint spanning sets, then by removing a vector from all but one of the disjoint spanning sets the frame is robust to at least R F 1 erasures. This observation together with Proposition 2.4 gives us both a lower and an upper bound for maximum robustness of a frame, namely R F 1 rob(f) min{ log2 S, k n}.

M.S. Copenhaver et al. / Linear Algebra and its Applications 439 (2013) 1330 1339 1335 3. Tight subframes and surgery on frames While there exist constructions that take a frame {f i } k i=1 H n and add vectors to produce a tight frame [6], it is not clear whether we can instead choose some subset {f i } i I {f i } k i=1 that is a tight frame? This is useful to consider for several reasons. We may be interested in having vectors of certain norms. While we may apply some method of construction that could possibly keep the norms of the added vectors within a specified range, removing vectors from a frame to produce a tight frame would leave the norms unaffected. This ultimately relies only on what we have, so it is possible that this trimming method may preserve some special features of the initial frame. Definition 3.1. Given a frame {f i } k i=1 for H n, we say that {f i } i I {f i } k i=1 is a tight subframe if {f i} i I is itself a tight frame for H n. We say that a (p, q)-surgery on {f i } k i=1 is possible if and only if there exist I {1, 2,...,k} with I =pand {g j } q j=1 H n so that {f i } i I C {g j } q j=1 is a tight frame for H n.when the relevant set {f i } k i=1 is clear, we may simply say that (p, q)-surgery is possible or impossible. In the case of surgery on unit-norm frames, we require that the new collection {g j } q j=1 H n contains only unit-norm vectors. Definition 3.2 ([7]). For any vector f R n, we define the diagram vector associated with f,denoted f,by f 2 (1) f 2 (2) f = 1 n 1. f 2 (n 1) f 2 (n) R n(n 1), 2nf (1)f (2).. 2nf (n 1)f (n) where the difference of squares f 2 (i) f 2 (j) and the product f (i)f (j) occur exactly once for i < j, i = 1, 2,, n 1. For any vector f C n, we define the diagram vector associated with f,denoted f,by f (1)f (1) f (2)f (2) f = 1 n 1. f (n 1)f (n 1) f (n)f (n) nf (1)f (2) C 3n(n 1)/2, nf (1)f (2). nf (n 1)f (n) nf (n 1)f (n) where the difference of the form f (i)f (i) f (j)f (j) occurs exactly once for i < j, i = 1, 2,, n 1 and the product of the form f (i)f (j) occurs exactly once for i = j. In order to give a formulation of tight subframes and surgeries, we present the following characterization of tight frames. We use the following remark in the next proposition.

1336 M.S. Copenhaver et al. / Linear Algebra and its Applications 439 (2013) 1330 1339 Remark 3.3. ([9], p.121) Let A, B be self-adjoint positive semidefinite matrices such that Ax, x = Bx, x for all x H n.thena = B. Proposition 3.4. Let k n, λ>0 and F = {f i } k i=1 be a sequence of vectors in Rn (or C n ),notallof which are zero. Let G be the Gramian associated with F.Let G be the Gramian associated with the diagram vectors { fi } k i=1. Then the following conditions are equivalent. (1) F is a λ-tight frame. (2) G has rank n and G 2 = λg. (3) σ(g) = (λ,...,λ, 0,...,0). }{{}}{{} ntimes k n times (4) k f i=1 i = 0. (5) (1, 1,, 1) T null( G). Proof. Let θ and S = θ θ be the synthesis and frame operators corresponding to {f i } k i=1, respectively. We recall that S = λi n if and only if F is a λ-tight frame. (1) (2) Since θ θ = λi n,wehaveθθ θθ = θλi n θ, which is equivalent to G 2 = λg. Also, rank(g) = rank(s) = n. (2) (1) Since rank(s) = rank(g) = n and S is an n n matrix, S is invertible. Hence θ is onto. Since G 2 = λg is equivalent to θsθ = θλi n θ,wehave(θ x) S(θ x) = (θ x) λi n (θ x) for all x R n (or C n ). By Remark 3.3, we conclude S = λi n. (1) (3) Since S = θ θ is an n n matrix and G = θθ is a k k matrix, the result now follows from σ(s) (0,...,0) = σ(g). }{{} k n times (1) (4) follows from Proposition 2.4 and 2.8 in [7]. (1) (5) follows from Proposition 4.2 in [7]. The following result gives a necessary condition for the existence of tight subframes. For more details about tight subframes we refer the reader to [13]. Proposition 3.5. Let F ={f i } k i=1 be a tight frame for Rn (or C n ).IfF has a tight subframe, then k 2n. Proof. If there exists a tight subframe {f i } i I F,thensince k i=1 f i = 0 and i I f i = 0, by condition (4) of Proposition 3.4, {f i } i I C is also a tight subframe. This implies that k = I + I C 2n since each subframe must span H n. We also observe that if {f i } k i=1 is a tight frame for Rn (or C n ) with n 2 and k n,thennontrivial (0, r), (r, 0)-surgery is impossible for r = 1, 2,...,n 1. The next proposition follows from this observation and condition (5) of Proposition 3.4. Proposition 3.6. If F ={f i } k i=1 is a unit-norm frame for Rn (or C n ), then the following conditions are equivalent. (1) There exists a unit-norm tight subframe { f j }j I F where I ={i 1,, i m }. (2) An (r, 0)-surgery on F which leaves a unit-norm tight frame is possible for some r {1, 2,...,k n}. (3) Each row sum of the Gramian of { f i1, f i2,..., f im } is zero. Next, we provide a generalization of the necessary condition for (p, q)-surgery presented in [14]. Theorem 3.1. Suppose F ={f i } k i=1 is a unit-norm tight frame for Rn (or C n ).Ifa(p, q)-surgery on F which leaves a unit-norm tight frame is possible then the sum of the entries of the Gramian ( f j, f i ) k p i,j=1 is

M.S. Copenhaver et al. / Linear Algebra and its Applications 439 (2013) 1330 1339 1337 bounded above by q 2,where f 1,..., f k p denote the diagram vectors that remain after removing p vectors from F. Proof. Let W = f 1 + + f k p. If a (p, q)-surgery on F is possible, then W q since f i = f. Note that W 2 = f 1 + + f k p, f 1 + + f k p = k p f i, f j. Thus, if (p, q)-surgery is possible, then the sum of the entries of the Gramian ( f j, f i ) k p i,j=1 is bounded above by q 2. We conclude our results on general surgeries on frames with the following theorem. Theorem 3.2. If (p, q)-surgery on a frame {f i } k i=1 H n is possible, then (r, r p+q)-surgery is possible for r = p + 1, p + 2,...,k. Proof. Suppose that (p, q)-surgery on {f i } k i=1 is possible. Without loss of generality, let {f 1,...,f k p } be the set of vectors remaining after removal of p vectors and let {g 1,...,g q } be the q vectors added to {f 1,...,f k p } so that {f i } k p i=1 {g j} q j=1 is a tight frame. Let r {p+1, p+2,...,k}.then(r, r p+q)- surgery on a frame {f i } k i=1 H n is possible, by excising {f k r+1, f k r+2,...,f k }, and then adding back the r p + q vectors {f k r+1, f k r+2,...,f k p } {g j } q j=1 so that i,j=1 {f 1, f 2,...,f k r } {f k r+1, f k r+2,...,f k p } {g j } q j=1 ={f i} k p i=1 {g j} q j=1, which is tight. 4. Length surgery While frame surgery involves adding vectors to a frame or removing vectors from a frame, length surgery deals only with the norms of vectors in a frame. Definition 4.1 ([14]). Let a 1, a 2,..., a k denote norms of k vectors in H n.a(p, q)-length surgery on {a i } k i=1 removes p numbers and replaces them with q nonnegative numbers. An interesting question while performing length surgery is whether the norms corresponding to the set {a i } k i=1 or any modified set are norms of vectors that form a tight frame. Definition 4.2. Asetofnumbers{a i } k i=1 is a tight frame set for H n if {a i } k i=1 contains the norms of vectors of a tight frame in H n. The following theorem describes exactly when a set of numbers is a tight frame set for H n. Theorem 4.1 ([4]). Given an n-dimensional Hilbert space H n and a sequence of positive scalars {a i } k i=1, there exists a tight frame {f i } k i=1 for H n of lengths f i =a i for all i = 1,...,k if and only if max a 2 i i 1 k a 2. i n i=1 (1) Due to its profound importance, the inequality (1) is often called the fundamental inequality [9]. An immediate consequence of Theorem 4.1 is that {a i } n i=1 is a tight frame set for H n if and only if a 1 = = a n.

1338 M.S. Copenhaver et al. / Linear Algebra and its Applications 439 (2013) 1330 1339 When performing (p, q)-length surgery,we would like to know when length surgery resulting in a tight frame set for H n is possible, and furthermore, specifically what numbers can be added to a set of nonnegative numbers for it to become or remain a tight frame set. Proposition 4.3. Suppose that {a i } k i=1 is a tight frame set for H n and let m q := max i=q+1,,k a 2 i q = 1, 2,...,k 1. Thenwehavethat for m q 1 n q In particular, we have k i=q+1 a 2 i. k nm 1 a 2 i m 1 1 k a 2 n. i 1 i=2 i=2 Proof. By Theorem 4.1,wehavemax i a 2 i 1 ( ) k n i=1 a2 i, which is equivalent to n max i a 2 q i i=1 a2 i ( ) k i=q+1 a2 i.sincen max i a 2 ( ) ( ) q i i=1 a2 i n max i a 2 i q max i a 2 i (n q)mq,wehavem q 1 k n q i=q+1 a2 i. The left inequality for m 1 is equivalent to m 1 1 k n 1 i=2 a2 i. In the following theorem we state what nonnegative numbers b should replace a 1 in a tight frame set {a i } k i=1 so that it remains a tight frame set. Theorem 4.2. Suppose that {a i } k i=1 is a tight frame set for H n and b is a nonnegative scalar. Then {b, a 2,, a k } is a tight frame set for H n if and only if b 2 [ nm k i=2 a2, 1 ] k i n 1 i=2 a2 i,where m := max i=2,,k a 2. i Proof. Let us assume b 2 m.now,b 2 nm ( k i=2 a2 i is equivalent to m 1 b 2 + ) k n i=2 a2 i.thisis equivalent to {b, a 2,, a k } being a tight frame set for H n by Theorem 4.1 since { m = max b 2, a 2,, } 2 a2 k. Suppose b 2 > m. Now,b 2 1 k n 1 i=2 a2 i, which is equivalent to ( b 2 1 b 2 + ) k n i=2 a2 i.thisisequivalentto{b, a 2,, a k } being a tight frame set for H n by Theorem 4.1 since b 2 = max { b 2, a 2 2,, a2 k }. The previous theorem states when (1, 1)-length surgery resulting in a tight frame set is possible. The next theorem describes when (p, q)-length surgery resulting in a tight frame set is possible. This result is a restatement of Theorem 4.9 in [14] which gives a condition for (0, q)-length surgery resulting in a tight frame set; the proof given here is different from that of [14]. Theorem 4.3. Let n 2, 0 p kand{a i } k i=1 be a sequence of positive scalars. Then (p, q)-length surgery on {a i } k i=1 resulting in a tight frame set is possible if and only if there exists a subset { } k p b j j=1 of {a i } k i=1 with max 1 j k p b 2 j 1 k p n q j=1 b2. j Proof. ( ) Thesetafter(p, q)-length surgery on {a i } k i=1 satisfies the desired inequality from Proposition 4.3 by letting {a i } q i=1 be the set of the numbers which are added to the set of k p remaining numbers. ( ) Choose b k p+i := ( ) 1/2 max 1 j k p b 2 j for i = 1,, q. Then max 1 j k p b 2 j 1 k p n q j=1 b2 j implies that (n q) max 1 j k p b 2 j k p j=1 b2 j. By adding q j=1 b2 k p+i to both sides of this inequality, we obtain the fundamental inequality for {b j} k p+q i=1.

M.S. Copenhaver et al. / Linear Algebra and its Applications 439 (2013) 1330 1339 1339 Acknowledgements Copenhaver, Logan, Mayfield, Narayan, and Sheperd were supported by the NSF-REU Grant DMS 08-51321. Kim was supported by the Central Michigan University ORSP Early Career Investigator (ECI) grant #C61373. The authors thank M. Petro and the referee for their valuable comments. References [1] B. Alexeev, J. Cahill, D.G. Mixon, Full spark frames, J. Fourier Anal. Appl. 18 (6) (2012) 1167 1194. [2] J.J. Benedetto, M. Fickus, Finite normalized tight frames, Adv. Comput. Math. 18 (2003) 357 385. [3] B.G. Bodmann, P.G. Casazza, G. Kutyniok, Upper and lower redundancy of finite frames, Proc. Conf. Inf. Sci. Syst. (2010) 1 6. [4] P.G. Casazza, M. Fickus, J. Kovačević, M.T. Leon, J.C. Tremain, A physical interpretation of finite frames, Appl. Numer. Harmon. Anal. 2 3 (2006) 51 76. [5] P.G. Casazza, J. Kovačević, Equal-norm tight frames with erasures, Adv. Comput. Math. 392 (2006) 387 430. [6] P.G. Casazza, N. Leonhard, Classes of finite equal norm Parseval frames, frames and operator theory in analysis and signal processing, in: Contemp. Math., vol. 451, Amer. Math. Soc., Providence, RI, 2008, pp. 11 31, [7] M.S. Copenhaver, Y. Kim, C. Logan, K. Mayfield, S.K. Narayan, M.J. Petro, J. Sheperd, Diagram vectors and tight frame scaling in finite dimensions, preprint. [8] I. Daubechies, Ten Lectures on Wavelets, SIAM, Philadelphia, PA, 1992. [9] H. Deguang, K. Kornelson, D. Larson, E. Weber, Frames for Undergraduates, in: Student Mathematical Library, vol. 40, American Mathematical Society, Providence, RI, 2007 [10] R.J. Duffin, A.C. Shaeffer, A class of nonharmonic Fourier series, Trans. Amer. Math. Soc. 72 (1952) 341 366. [11] M. Fickus, D.G. Mixon, Numerically erasure-robust frames, Linear Algebra Appl. 437 (2012) 1394 1407. [12] R.B. Holmes, V.I. Paulsen, Optimal frames for erasures, Linear Algebra Appl. 377 (2004) 31 51. [13] J. Lemvig, C. Miller, K.A. Okoudjou, Prime tight frames, preprint. [14] S.K. Narayan, E.L. Radzwion, S.P. Rimer, R.L. Tomasino, J.L. Wolfe, A.M. Zimmer, Robustness and surgery of frames, Linear Algebra Appl. 434 (2011) 1893 1901. [15] N.Q. Nguyen, Surgery on frames, Ph.D. thesis, Texas A&M University, August 2008. [16] M. Puschel, J. Kovačević, Real, tight frames with maximal robustness to erasures, Proc. Data Compr. Conf. (2005) 63 72. [17] A. Ron, Z. Shen, Frames and stable bases for shift-invariant subspaces of L 2 (R d ), Canad. J. Math. 47 (5) (1995) 1051 1094.