Sparse representations from moments
|
|
- Bruno Claud Wilkins
- 5 years ago
- Views:
Transcription
1 Sparse representations from moments Bernard Mourrain Inria Méditerranée, Sophia Antipolis
2 Sparse representation of signals Given a function or signal f (t): decompose it as f (t) = (a i cos(µ i t) + b i sin(µ i t))e ν i t = r i=1 r ω i e ζ it i=1 B Mourrain Sparse representations from moments 2 / 29
3 Prony s method (1795) For the signal f (t) = r i=1 ω ie ζit, (ω i, ζ i C), Evaluate f at 2 r regularly spaced points: σ 0 := f (0), σ 1 := f (1), Compute a non-zero element p = [p 0,, p r ] in the kernel: σ 0 σ 1 σ r σ 1 σ r+1 σ r 1 σ 2r 1 σ 2r 1 Compute the roots ξ 1 = e ζ 1,, ξ r = e ζr Solve the system 1 1 ξ 1 ξ r ξ r 1 1 ξr r 1 ω 1 ω 2 ω r p 0 p 1 p r = 0 of p(x) := r i=0 p ix i = σ 0 σ 1 σ r 1 B Mourrain Sparse representations from moments 3 / 29
4 Symmetric tensor decomposition and Waring problem (1770) Tensor decomposition problem: given a homogeneous polynomial Ψ of degree d in the variables x = (x 0, x 1,, x n ): Ψ(x) = ( ) d σ α x α, α α =d find a minimal decomposition of Ψ of the form Ψ(x) = r ω i (ξ i,0 x 0 + ξ i,1 x ξ i,n x n ) d i=1 for ξ i = (ξ i,0, ξ i,1,, ξ i,n ) C n+1, ω i C The minimal r in such a decomposition is called the rank of τ B Mourrain Sparse representations from moments 4 / 29
5 Sylvester approach (1851) Theorem The binary form Ψ(x 0, x 1 ) = d i=0 σ d ) i( i x d i sum of r distinct powers of linear forms 0 x i 1 can be decomposed as a Ψ = r ω k (α k x 0 + β k x 1 ) d k=1 iff there exists a polynomial p(x 0, x 1 ) := p 0 x0 r + p 1x0 r 1 x p r x1 r st σ 0 σ 1 σ r p 0 σ 1 σ r+1 p 1 σ d r σ d 1 σ d and of the form p = c r k=1 (β kx 0 α k x 1 ) p r = 0 B Mourrain Sparse representations from moments 5 / 29
6 Sparse interpolation Given a black-box polynomial function f (x) find what are the terms inside from output values Find r N, ω i C, α i N such that f (x) = r i=1 ω i x α i B Mourrain Sparse representations from moments 6 / 29
7 Choose ϕ C Compute the sequence of terms σ 0 = f (1),, σ 2r 1 = f (ϕ 2r 1 ); Construct the matrix H = [σ i+j ] and its kernel p = [p 0,, p r ] st σ 0 σ 1 σ r σ 1 σ r+1 σ r 1 σ 2r 1 σ 2r 1 Compute the roots ξ 1 = ϕ α 1,, ξ r = ϕ αr deduce the exponents α i = log ϕ (ξ i ) p 0 p 1 p r = 0 of p(x) := r i=0 p ix i and Deduce the weights W = [ω i ] by solving V Ξ W = [σ 0,, σ r 1 ] where V Ξ is the Vandermonde system of the roots ξ 1,, ξ r B Mourrain Sparse representations from moments 7 / 29
8 Decoding An algebraic code: E = {c(f ) = [f (ξ 1 ),, f (ξ m )] f K[x]; deg(f ) d} Encoding messages using the dual code: C = E = {c c [f (ξ 1 ),, f (ξ l )] = 0 f V = x a F[x]} Message received: r = m + e for m C where e = [ω 1,, ω m ] is an error with ω j 0 for j = i 1,, i r and ω j = 0 otherwise Find the error e B Mourrain Sparse representations from moments 8 / 29
9 Berlekamp-Massey method (1969) Compute the syndrome σ k = c(x k ) r = c(x k ) e = r j=1 ω i j ξi k j Compute the matrix σ 0 σ 1 σ r σ 1 σ r+1 σ r 1 σ 2r 1 σ 2r 1 and its kernel p = [p 0,, p r ] p 0 p 1 Compute the roots of the error locator polynomial p(x) = r i=0 p i x i = p r r j=1 (x ξ i j ) Deduce the errors ω ij p r = 0 B Mourrain Sparse representations from moments 9 / 29
10 Sequences, series, duality (1D) Sequences: σ = (σ k ) k N K N indexed by k N Formal power series: σ = y k σ k k! K[[y]] (or σ k z k K[[z]]) k=0 k=0 Linear functionals: K[x] = {Λ : K[x] K linear} Example: p coefficient of x i in p = 1 i! i (p)(0) e ζ : p p(ζ) B Mourrain Sparse representations from moments 10 / 29
11 Series as linear functionals: For σ = k=0 σ k y k k! K[[y]], σ : p = k p k x k σ p = k σ k p k ( y k k! ) is the dual basis of the monomial basis (x k ) k N Example: e ζ = k=0 ζk y k k! = eζy K[[y]] Structure of K[x]-module: p(x) Λ : q Λ(p q) x σ = k=1 σ k y k 1 p(x) σ = p( )(σ(y)) (k 1)! = (σ(y)) B Mourrain Sparse representations from moments 11 / 29
12 Sequences, series, duality (nd) Multi-index sequences: σ = (σ α ) α N n K Nn α = (α 1,, α n ) N n Taylor series: indexed by σ(y) = y α σ α α! K[[y 1,, y n ]] (or σ α z α K[[z 1,, z n ]]) α N n α N n where α! = α i! for α = (α 1,, α n ) N Linear functionals: σ R = {σ : R K, linear} σ : p = α p α x α σ p = α σ α p α The coefficients σ x α = σ α K, α N n are called the moments of σ Structure of R-module: p R, σ R, p σ : q σ p q p σ = p( 1,, n )(σ)(y) B Mourrain Sparse representations from moments 12 / 29
13 Hankel operator: For σ K[[y]], H σ : K[x] K[[y]] p p σ = p( 1,, n )(σ)(y) σ is the symbol of H σ Truncated Hankel operator: V, V K[x], H V,V σ : p V p σ V V Example: V = x A, V = x B K[x] Ideal: H A,B σ = [σ α+β ] α A,β B I σ = ker H σ = {p K[x] p σ = 0}, = {p = α p α x α β N n α p α σ α+β = 0} Linear recurrence relations on the sequence σ = (σ α ) α N n Quotient algebra: A σ = K[x 1,, x n ]/I σ Studied case: dim A σ < B Mourrain Sparse representations from moments 13 / 29
14 Structure of an Artinian algebra A Hypothesis: A = K[x]/I is Artinian, ie dim K A < Hilbert nullstellensatz: A = K[x]/I Artinian V K (I ) = {ξ 1,, ξ r } is finite Assuming K = K is algebraically closed, we have I = Q 1 Q r where Q i is m ξi -primary where V K (I ) = {ξ 1,, ξ r } A = K[x]/I = A 1 A r, with Ai = u i A K[x 1,, x n ]/Q i, u 2 i = u i, u i u j = 0 if i j, u u r = 1 dim R/Q i = µ i is the multiplicity of ξ i B Mourrain Sparse representations from moments 14 / 29
15 Structure of the dual A Sparse series: PolExp = { σ(y) = } r ω i (y) e ξi (y) ω i (y) K[y], i=1 where e ξi (y) = e y ξ i = e y 1ξ 1,i + +y nξ n,i with ξ i,j K Inverse system generated by ω 1,, ω r K[y] Theorem For K = K algebraically closed, V K (I ) = {ξ 1,, ξ r } D i = D(ω i,1,, ω i,li ) = Q i D(ω 1,, ω r ) = α y (ω i ), α N n A = r i=1 D i e ξi (y) PolExp with ω i,j K[y] and I = Q 1 Q r µ(ω i,1,, ω i,li ) := dim K (D i ) = µ i multiplicity of ξ i B Mourrain Sparse representations from moments 15 / 29
16 The roots by eigencomputation Hypothesis: V K (I ) = {ξ 1,, ξ r } A = K[x]/I Artinian M a : A A u a u M t a : A A Λ a Λ = Λ M a Theorem The eigenvalues of M a are {a(ξ 1 ),, a(ξ r )} The eigenvectors of all (M t a) a A are (up to a scalar) e ξi : p p(ξ i ) Proposition If the roots are simple, the operators M a are diagonalizable Their common eigenvectors are, up to a scalar, interpolation polynomials u i at the roots and idempotent in A B Mourrain Sparse representations from moments 16 / 29
17 Example Roots of polynomial systems { f1 = x 2 1 x 2 x 2 1 I = (f 1, f 2 ) C[x] f 2 = x 1 x 2 x 2 A = C[x]/I 1, x 1, x 2 I = (x 2 1 x 2, x 1 x 2 x 2, x 2 2 x 2) M 1 = , M 2 = common eigvecs of M t 1, Mt 2 = 1 0 0, I = Q 1 Q 2 where Q 1 = (x 2 1, x 2), Q 2 = m (1,1) = (x 1 1, x 2 1) I = Q 1 Q 2 Q 1 = 1, y 1 = 1, y 1 e (0,0) (y) Q 2 = 1 e (1,1) (y) = e y1+y2 Solution of partial differential equations (with constant coeff) { 2 y1 y2 σ 2 y 1 σ = 0 f 1 σ = 0 σ I = Q 1 Q 2 y1 y2 σ y2 σ = 0 f 2 σ = 0 σ = a + b y 1 + c e y1+y2 a, b, c C B Mourrain Sparse representations from moments 17 / 29
18 Hankel operators and quotient algebra Hankel operator: For σ K[[y]], H σ : K[x] K[[y]] p p σ = p( 1,, n )(σ)(y) Quotient algebra: A σ = K[x 1,, x n ]/I σ A = K[x 1,, x n ]/I Artinian iff dim K A < A Gorenstein iff σ A such that A = A σ is a free A Isomorphism with the dual space A σ: 0 I σ K[x] H σ A σ 0 p p σ Correspondence between sequences σ K Nn with rankh σ < and Artinian Gorenstein algebras A σ := K[x]/I σ B Mourrain Sparse representations from moments 18 / 29
19 Univariate series: Kronecker (1881) Multivariate series: The Hankel operator H σ : (p m ) C N,finite ( m σ m+np m ) n N C N is of finite rank r iff ω 1,, ω r C[y] and ξ 1,, ξ r C distincts st σ(y) = y n r σ n n! = ω i (y)e ξi (y) with (deg(ω i )+1) = r n N i=1 i=1 Theorem (Generalized Kronecker Theorem) H σ is of rank r iff σ = r i=1 ω i(y) e ξi (y) PolExp with r = r i=1 µ(ω i) In this case, we have V C (I σ ) = {ξ 1,, ξ r } I σ = Q 1 Q r with Q i = D(ω i ) e ξi (y) A σ = A σ σ (free A σ -module of rank 1) (a, b) σ ab is non-degenerate in A σ B Mourrain Sparse representations from moments 19 / 29 r
20 The structure of A σ For σ = r i=1 ω i e ξi, with ω i C \ {0} and ξ i C n distinct rank H σ = r and the multiplicity of the points ξ 1,, ξ r in V(I σ ) is 1 For B, B be of size r, H B,B σ A σ = K[x]/I σ invertible iff B and B are bases of The matrix M i of multiplication by x i in the basis B of A σ is such that H B,x i B σ = H B,B x i σ = H B,B σ The common eigenvectors of M i are (up to a scalar) the Lagrange interpolation polynomials { u ξi at the points ξ i, i = 1,, r 1 ifi = j, u ξi (ξ j )= u 0 otherwise 2 ξi u ξi, r i=1 u ξ i 1 The common eigenvectors of Mi t [B(ξ i )], i = 1,, r M i are (up to a scalar) the vectors B Mourrain Sparse representations from moments 20 / 29
21 Decomposition algorithm Input: The first coefficients (σ α ) α A of the series σ = y α r σ α α! = ω i e ξi (y) α N n i=1 1 Compute bases B, B x A st that H B,B invertible and B = B = r = dim A σ ; 2 Deduce the tables of multiplications M i := (H B,B σ ) 1 H B,x i B σ 3 Compute the eigenvectors v 1,, v r of i l im i for a generic l = l 1 x l n x n ; 4 Deduce the points ξ i = (ξ i,1,, ξ i,n ) st M j v i ξ i,j v i = 0 and the weights ω i = 1 v i (ξ i ) σ v i Output: The decomposition σ = r 1 i=1 v i (ξ i ) σ v i e ξi (y) B Mourrain Sparse representations from moments 21 / 29
22 Multivariate Prony method (1) Let h(t 1, t 2 ) = t 1 2 t 2 3 t 1, σ = yα α N2 h(α) α! = 2e (1,1) (y) + 3e (2,2) (y) e (3,1) (y) Take B = {1, x 1, x 2 } and compute h(0, 0) h(1, 0) h(0, 1) H 0 := Hσ B,B = h(1, 0) h(2, 0) h(1, 1) H 1 := H B,x1 B σ = h(0, 1) h(1, 1) h(0, 2) =, H 2 := H B,x2 B σ = , Compute the generalized eigenvectors of (ah 1 + bh 2, H 0 ): U = 1/2 0 1/2 and H 0 U = /2 1 1/ This yields the weights 2, 3, 1 and the roots (1, 1), (2, 2), (3, 1) B Mourrain Sparse representations from moments 22 / 29
23 Multivariate Prony method (2) h(t 1, t 2 ) := r i=1 ω ie f 1t 1 +f 2 t 2 with f := i i i i i ω := i i i i For the sampling [ 1 50, ], B = {1, x 1, x 2, x 1 x 2 }, the SVD of Hσ B,B 92 [ , , , ] and the computed decomposition is f = i i i i i i i i i ω i = i i B Mourrain Sparse representations from moments 23 / 29 is
24 Sparse interpolation f (x) = r i=1 ω i x α i σ = γ f (ϕγ ) yγ γ! = r i=1 ω i e ϕ α i (y) Example: f (x 1, x 2 ) = x 33 1 x x 1x Compute σ α = f (ϕ α 1 1, ϕα 2 2 ) for α 1 + α 2 3 and ϕ 1 = ϕ 2 = e 2iπ Compute the Hankel matrix H 1,2 σ : i i i i i i i i i i i i i i 9 Deduce the decomposition of σ = 3 i=1 ω ie ξi : Ξ = i i e 11 i e 10 i i i ω = and the exponents 50 Ξ 2πi mod 50 of the terms of f : e 7 i e 6 i, e e 8 i e e 8 i e 6 i e 6 i e 7 i e 7 i e 8 i B Mourrain Sparse representations from moments 24 / 29
25 Symmetric tensor and apolarity Apolar product: For f = α =d f α ( d α) x α, g = α =d g α ( d α) x α K[x] d, f, g d = α d ( ) d f α g α α Property: f, (ξ 0 x ξ n x n ) d = f (ξ 0,, ξ n ) Duality: τ = α N n+1, α =d ( ) d τ α x α τ = α where (α 1,, α n ) = (d i α i, α 1,, α n ) N n+1 α N n, α d y α τ α α! τ(x) = r i=1 ω i (1 + ξ i,1 x ξ i,n x n ) d τ (y) = r i=1 ω i e ξi (y) + (y) d+1 B Mourrain Sparse representations from moments 25 / 29
26 Symmetric tensor decomposition τ = (x 0 x x 2 ) 4 + (x 0 + x 1 + x 2 ) 4 3 (x x x 2 ) 4 = x x 3 0 x 1 8 x 3 0 x 2 60 x 2 0 x x 2 0 x 1 x 2 12 x x 2 96 x 0 x x 0 x 2 1 x x 0 x 1 x x 0 x x x 3 1 x x 2 1 x x 1 x x 2 τ = e ( 1,3) (y) + e (1,1) (y) 3e (2,2) (y) (by apolarity) For B = {1, x 2, x 1 }, H 2,2 τ := H B,B τ = , HB,x 1 B τ = , HB,x 2 B τ = B Mourrain Sparse representations from moments 26 / 29
27 The matrix of multiplication by x 1 in B = {1, x 2, x 1 } is M 1 = (H B,B τ ) 1 H B,x 1 B τ = Its eigenvalues are [ 1, 1, 2] and the eigenvectors: U := that is the polynomials U(x) = [ 1 2 x x x x x x 1 We deduce the weights and the frequencies: H [1,x 1,x 2 ],U τ = Weights: 1, 1, 3; frequencies: ( 1, 3), (1, 1), (2, 2) Decomposition: τ (y) = e ( 1,3) (y) + e (1,1) (y) 3 e (2,2) (y) + (y) 4 ] B Mourrain Sparse representations from moments 27 / 29
28 A general framework F the functional space, in which the signal lives S 1,, S n : F F commuting linear operators: S i S j = S j S i : h F [h] C a linear functional on F Generating series associated to h F: Eigenfunctions: Generalized eigenfunctions: σ h (y) = α N n [S α (h)] yα α! = α N n σ α y α α! S j (E) = ξ j E, j = 1,, n σ E = e ξ (y) S j (E k ) = ξ j E k + k <k m j,k E k ) σ Ek = ω i (y)e ξ (y) If h σ h is injective unique decomposition of f as a linear combination of generalized eigenfunctions B Mourrain Sparse representations from moments 28 / 29
29 Other applications Decomposition of measures as sums of spikes from moments (images, spectroscopy, radar, astronomy, ) Decomposition of convolution operators of finite rank Sparse interpolation of PolyLog functions Polynomial optimisation and convex relaxations Vanishing ideal of points: σ = r i=1 e ξ i (y) Change of ordering for Grobner basis or basis representation for zero-dimensional (Gorenstein) ideals: σ α = u, N(x α ), Challenges, open questions Numerical stability, correction of errors, Efficient construction of basis, complexity, Super-resolution, collision of points, Super-extrapolation, Best low rank approximation, Thanks for your attention B Mourrain Sparse representations from moments 29 / 29
arxiv: v3 [math.ag] 20 Oct 2017
Polynomial-exponential decomposition from moments arxiv:160905720v3 [mathag] 20 Oct 2017 Bernard Mourrain Université Côte d Azur, Inria, aromath, France bernardmourrain@inriafr October 23, 2017 Abstract
More informationTensor decomposition and moment matrices
Tensor decomposition and moment matrices J. Brachat (Ph.D.) & B. Mourrain GALAAD, INRIA Méditerranée, Sophia Antipolis Bernard.Mourrain@inria.fr SIAM AAG, Raleigh, NC, Oct. -9, 211 Collaboration with A.
More informationDecomposition of Low Rank Multi-Symmetric Tensor
Decomposition of Low Rank Multi-Symmetric ensor Jouhayna Harmouch, Bernard Mourrain, Houssam Khalil o cite this version: Jouhayna Harmouch, Bernard Mourrain, Houssam Khalil. Decomposition of Low Rank Multi-
More informationReal solving polynomial equations with semidefinite programming
.5cm. Real solving polynomial equations with semidefinite programming Jean Bernard Lasserre - Monique Laurent - Philipp Rostalski LAAS, Toulouse - CWI, Amsterdam - ETH, Zürich LAW 2008 Real solving polynomial
More informationRanks of Real Symmetric Tensors
Ranks of Real Symmetric Tensors Greg Blekherman SIAM AG 2013 Algebraic Geometry of Tensor Decompositions Real Symmetric Tensor Decompositions Let f be a form of degree d in R[x 1,..., x n ]. We would like
More informationA Sparse Flat Extension Theorem for Moment Matrices
A Sparse Flat Extension Theorem for Moment Matrices Monique Laurent and Bernard Mourrain Abstract. In this note we prove a generalization of the flat extension theorem of Curto and Fialkow [4] for truncated
More information1. Let r, s, t, v be the homogeneous relations defined on the set M = {2, 3, 4, 5, 6} by
Seminar 1 1. Which ones of the usual symbols of addition, subtraction, multiplication and division define an operation (composition law) on the numerical sets N, Z, Q, R, C? 2. Let A = {a 1, a 2, a 3 }.
More informationHomework 2 - Math 603 Fall 05 Solutions
Homework 2 - Math 603 Fall 05 Solutions 1. (a): In the notation of Atiyah-Macdonald, Prop. 5.17, we have B n j=1 Av j. Since A is Noetherian, this implies that B is f.g. as an A-module. (b): By Noether
More informationQuadraticity and Koszulity for Graded Twisted Tensor Products
Quadraticity and Koszulity for Graded Twisted Tensor Products Peter Goetz Humboldt State University Arcata, CA 95521 September 9, 2017 Outline Outline I. Preliminaries II. Quadraticity III. Koszulity
More informationRes + X F F + is defined below in (1.3). According to [Je-Ki2, Definition 3.3 and Proposition 3.4], the value of Res + X
Theorem 1.2. For any η HH (N) we have1 (1.1) κ S (η)[n red ] = c η F. Here HH (F) denotes the H-equivariant Euler class of the normal bundle ν(f), c is a non-zero constant 2, and is defined below in (1.3).
More information12. Hilbert Polynomials and Bézout s Theorem
12. Hilbert Polynomials and Bézout s Theorem 95 12. Hilbert Polynomials and Bézout s Theorem After our study of smooth cubic surfaces in the last chapter, let us now come back to the general theory of
More informationADVANCED TOPICS IN ALGEBRAIC GEOMETRY
ADVANCED TOPICS IN ALGEBRAIC GEOMETRY DAVID WHITE Outline of talk: My goal is to introduce a few more advanced topics in algebraic geometry but not to go into too much detail. This will be a survey of
More informationChap. 3. Controlled Systems, Controllability
Chap. 3. Controlled Systems, Controllability 1. Controllability of Linear Systems 1.1. Kalman s Criterion Consider the linear system ẋ = Ax + Bu where x R n : state vector and u R m : input vector. A :
More informationOn deflation and multiplicity structure
On deflation and multiplicity structure Jonathan D. Hauenstein a, Bernard Mourrain b, Agnes Szanto c a Department of Applied and Computational Mathematics and Statistics, University of Notre Dame, Notre
More informationEigenvectors of Tensors and Waring Decomposition
Eigenvectors of Tensors and Waring Decomposition Luke Oeding University of California Berkeley Oeding, Ottaviani (UCB, Firenze) Giorgio Ottaviani Universita degli Studi di Firenze Waring Decomposition
More informationFast Polynomial Multiplication
Fast Polynomial Multiplication Marc Moreno Maza CS 9652, October 4, 2017 Plan Primitive roots of unity The discrete Fourier transform Convolution of polynomials The fast Fourier transform Fast convolution
More informationSparse polynomial interpolation: sparse recovery, super resolution, or Prony?
Sparse polynomial interpolation: sparse recovery, super resolution, or Prony? Cédric Josz, Jean-Bernard Lasserre, Bernard Mourrain To cite this version: Cédric Josz, Jean-Bernard Lasserre, Bernard Mourrain.
More informationAn Interpolation Algorithm for List Decoding of Reed-Solomon Codes
An Interpolation Algorithm for List Decoding of Reed-Solomon Codes Kwankyu Lee Department of Mathematics San Diego State University San Diego, USA Email: kwankyu@sogangackr Michael E O Sullivan Department
More informationALGEBRA QUALIFYING EXAM PROBLEMS LINEAR ALGEBRA
ALGEBRA QUALIFYING EXAM PROBLEMS LINEAR ALGEBRA Kent State University Department of Mathematical Sciences Compiled and Maintained by Donald L. White Version: August 29, 2017 CONTENTS LINEAR ALGEBRA AND
More informationx 3y 2z = 6 1.2) 2x 4y 3z = 8 3x + 6y + 8z = 5 x + 3y 2z + 5t = 4 1.5) 2x + 8y z + 9t = 9 3x + 5y 12z + 17t = 7
Linear Algebra and its Applications-Lab 1 1) Use Gaussian elimination to solve the following systems x 1 + x 2 2x 3 + 4x 4 = 5 1.1) 2x 1 + 2x 2 3x 3 + x 4 = 3 3x 1 + 3x 2 4x 3 2x 4 = 1 x + y + 2z = 4 1.4)
More informationChapter 6 Reed-Solomon Codes. 6.1 Finite Field Algebra 6.2 Reed-Solomon Codes 6.3 Syndrome Based Decoding 6.4 Curve-Fitting Based Decoding
Chapter 6 Reed-Solomon Codes 6. Finite Field Algebra 6. Reed-Solomon Codes 6.3 Syndrome Based Decoding 6.4 Curve-Fitting Based Decoding 6. Finite Field Algebra Nonbinary codes: message and codeword symbols
More informationAPPENDIX A. Background Mathematics. A.1 Linear Algebra. Vector algebra. Let x denote the n-dimensional column vector with components x 1 x 2.
APPENDIX A Background Mathematics A. Linear Algebra A.. Vector algebra Let x denote the n-dimensional column vector with components 0 x x 2 B C @. A x n Definition 6 (scalar product). The scalar product
More informationA Brief Introduction to Functional Analysis
A Brief Introduction to Functional Analysis Sungwook Lee Department of Mathematics University of Southern Mississippi sunglee@usm.edu July 5, 2007 Definition 1. An algebra A is a vector space over C with
More informationMatrix Theory. A.Holst, V.Ufnarovski
Matrix Theory AHolst, VUfnarovski 55 HINTS AND ANSWERS 9 55 Hints and answers There are two different approaches In the first one write A as a block of rows and note that in B = E ij A all rows different
More informationMATH 8253 ALGEBRAIC GEOMETRY WEEK 12
MATH 8253 ALGEBRAIC GEOMETRY WEEK 2 CİHAN BAHRAN 3.2.. Let Y be a Noetherian scheme. Show that any Y -scheme X of finite type is Noetherian. Moreover, if Y is of finite dimension, then so is X. Write f
More informationThe converse is clear, since
14. The minimal polynomial For an example of a matrix which cannot be diagonalised, consider the matrix ( ) 0 1 A =. 0 0 The characteristic polynomial is λ 2 = 0 so that the only eigenvalue is λ = 0. The
More informationMath 113 Winter 2013 Prof. Church Midterm Solutions
Math 113 Winter 2013 Prof. Church Midterm Solutions Name: Student ID: Signature: Question 1 (20 points). Let V be a finite-dimensional vector space, and let T L(V, W ). Assume that v 1,..., v n is a basis
More informationContents. Preface for the Instructor. Preface for the Student. xvii. Acknowledgments. 1 Vector Spaces 1 1.A R n and C n 2
Contents Preface for the Instructor xi Preface for the Student xv Acknowledgments xvii 1 Vector Spaces 1 1.A R n and C n 2 Complex Numbers 2 Lists 5 F n 6 Digression on Fields 10 Exercises 1.A 11 1.B Definition
More informationMATH32062 Notes. 1 Affine algebraic varieties. 1.1 Definition of affine algebraic varieties
MATH32062 Notes 1 Affine algebraic varieties 1.1 Definition of affine algebraic varieties We want to define an algebraic variety as the solution set of a collection of polynomial equations, or equivalently,
More informationCHAPTER 1. AFFINE ALGEBRAIC VARIETIES
CHAPTER 1. AFFINE ALGEBRAIC VARIETIES During this first part of the course, we will establish a correspondence between various geometric notions and algebraic ones. Some references for this part of the
More informationMath Linear Algebra II. 1. Inner Products and Norms
Math 342 - Linear Algebra II Notes 1. Inner Products and Norms One knows from a basic introduction to vectors in R n Math 254 at OSU) that the length of a vector x = x 1 x 2... x n ) T R n, denoted x,
More informationSOS TENSOR DECOMPOSITION: THEORY AND APPLICATIONS
SOS TENSOR DECOMPOSITION: THEORY AND APPLICATIONS HAIBIN CHEN, GUOYIN LI, AND LIQUN QI Abstract. In this paper, we examine structured tensors which have sum-of-squares (SOS) tensor decomposition, and study
More informationAbstract Vector Spaces and Concrete Examples
LECTURE 18 Abstract Vector Spaces and Concrete Examples Our discussion of linear algebra so far has been devoted to discussing the relations between systems of linear equations, matrices, and vectors.
More informationFrobenius maps on injective hulls and their applications Moty Katzman, The University of Sheffield.
Frobenius maps on injective hulls and their applications Moty Katzman, The University of Sheffield. Consider a complete local ring (S, m) of characteristic p > 0 and its injective hull E S = E S (S/m).
More informationNumerical Solutions to Partial Differential Equations
Numerical Solutions to Partial Differential Equations Zhiping Li LMAM and School of Mathematical Sciences Peking University Introduction to Hyperbolic Equations The Hyperbolic Equations n-d 1st Order Linear
More informationLinear Algebra for Computing Gröbner Bases of Linear Recursive Multidimensional Sequences
Linear Algebra for Computing Gröbner Bases of Linear Recursive Multidimensional Sequences Jérémy Berthomieu, Brice Boyer, Jean-Charles Faugère To cite this version: Jérémy Berthomieu, Brice Boyer, Jean-Charles
More informationRECURSIVE RELATIONS FOR THE HILBERT SERIES FOR CERTAIN QUADRATIC IDEALS. 1. Introduction
RECURSIVE RELATIONS FOR THE HILBERT SERIES FOR CERTAIN QUADRATIC IDEALS 1 RECURSIVE RELATIONS FOR THE HILBERT SERIES FOR CERTAIN QUADRATIC IDEALS AUTHOR: YUZHE BAI SUPERVISOR: DR. EUGENE GORSKY Abstract.
More information10. Noether Normalization and Hilbert s Nullstellensatz
10. Noether Normalization and Hilbert s Nullstellensatz 91 10. Noether Normalization and Hilbert s Nullstellensatz In the last chapter we have gained much understanding for integral and finite ring extensions.
More informationKernel Method: Data Analysis with Positive Definite Kernels
Kernel Method: Data Analysis with Positive Definite Kernels 2. Positive Definite Kernel and Reproducing Kernel Hilbert Space Kenji Fukumizu The Institute of Statistical Mathematics. Graduate University
More informationAlgebraic Varieties. Notes by Mateusz Micha lek for the lecture on April 17, 2018, in the IMPRS Ringvorlesung Introduction to Nonlinear Algebra
Algebraic Varieties Notes by Mateusz Micha lek for the lecture on April 17, 2018, in the IMPRS Ringvorlesung Introduction to Nonlinear Algebra Algebraic varieties represent solutions of a system of polynomial
More informationMATH 326: RINGS AND MODULES STEFAN GILLE
MATH 326: RINGS AND MODULES STEFAN GILLE 1 2 STEFAN GILLE 1. Rings We recall first the definition of a group. 1.1. Definition. Let G be a non empty set. The set G is called a group if there is a map called
More informationLinear Algebra- Final Exam Review
Linear Algebra- Final Exam Review. Let A be invertible. Show that, if v, v, v 3 are linearly independent vectors, so are Av, Av, Av 3. NOTE: It should be clear from your answer that you know the definition.
More informationCOHEN-MACAULAY RINGS SELECTED EXERCISES. 1. Problem 1.1.9
COHEN-MACAULAY RINGS SELECTED EXERCISES KELLER VANDEBOGERT 1. Problem 1.1.9 Proceed by induction, and suppose x R is a U and N-regular element for the base case. Suppose now that xm = 0 for some m M. We
More informationON THE FIRST FALL DEGREE OF SUMMATION POLYNOMIALS
ON THE FIRST FALL DEGREE OF SUMMATION POLYNOMIALS STAVROS KOUSIDIS AND ANDREAS WIEMERS Abstract We improve on the first fall degree bound of polynomial systems that arise from a Weil descent along Semaev
More informationCourse 311: Michaelmas Term 2005 Part III: Topics in Commutative Algebra
Course 311: Michaelmas Term 2005 Part III: Topics in Commutative Algebra D. R. Wilkins Contents 3 Topics in Commutative Algebra 2 3.1 Rings and Fields......................... 2 3.2 Ideals...............................
More informationLocal properties of plane algebraic curves
Chapter 7 Local properties of plane algebraic curves Throughout this chapter let K be an algebraically closed field of characteristic zero, and as usual let A (K) be embedded into P (K) by identifying
More informationPolynomial interpolation over finite fields and applications to list decoding of Reed-Solomon codes
Polynomial interpolation over finite fields and applications to list decoding of Reed-Solomon codes Roberta Barbi December 17, 2015 Roberta Barbi List decoding December 17, 2015 1 / 13 Codes Let F q be
More informationhttps://goo.gl/kfxweg KYOTO UNIVERSITY Statistical Machine Learning Theory Sparsity Hisashi Kashima kashima@i.kyoto-u.ac.jp DEPARTMENT OF INTELLIGENCE SCIENCE AND TECHNOLOGY 1 KYOTO UNIVERSITY Topics:
More information5 Quiver Representations
5 Quiver Representations 5. Problems Problem 5.. Field embeddings. Recall that k(y,..., y m ) denotes the field of rational functions of y,..., y m over a field k. Let f : k[x,..., x n ] k(y,..., y m )
More informationLecture 2. (1) Every P L A (M) has a maximal element, (2) Every ascending chain of submodules stabilizes (ACC).
Lecture 2 1. Noetherian and Artinian rings and modules Let A be a commutative ring with identity, A M a module, and φ : M N an A-linear map. Then ker φ = {m M : φ(m) = 0} is a submodule of M and im φ is
More informationHilbert function, Betti numbers. Daniel Gromada
Hilbert function, Betti numbers 1 Daniel Gromada References 2 David Eisenbud: Commutative Algebra with a View Toward Algebraic Geometry 19, 110 David Eisenbud: The Geometry of Syzygies 1A, 1B My own notes
More informationEach is equal to CP 1 minus one point, which is the origin of the other: (C =) U 1 = CP 1 the line λ (1, 0) U 0
Algebraic Curves/Fall 2015 Aaron Bertram 1. Introduction. What is a complex curve? (Geometry) It s a Riemann surface, that is, a compact oriented twodimensional real manifold Σ with a complex structure.
More informationNOTES ON MODULES AND ALGEBRAS
NOTES ON MODULES AND ALGEBRAS WILLIAM SCHMITT 1. Some Remarks about Categories Throughout these notes, we will be using some of the basic terminology and notation from category theory. Roughly speaking,
More information4.5 Hilbert s Nullstellensatz (Zeros Theorem)
4.5 Hilbert s Nullstellensatz (Zeros Theorem) We develop a deep result of Hilbert s, relating solutions of polynomial equations to ideals of polynomial rings in many variables. Notation: Put A = F[x 1,...,x
More informationFinite Elements. Colin Cotter. January 18, Colin Cotter FEM
Finite Elements January 18, 2019 The finite element Given a triangulation T of a domain Ω, finite element spaces are defined according to 1. the form the functions take (usually polynomial) when restricted
More informationGeometric Modeling Summer Semester 2012 Linear Algebra & Function Spaces
Geometric Modeling Summer Semester 2012 Linear Algebra & Function Spaces (Recap) Announcement Room change: On Thursday, April 26th, room 024 is occupied. The lecture will be moved to room 021, E1 4 (the
More informationContradiction. Theorem 1.9. (Artin) Let G be a finite group of automorphisms of E and F = E G the fixed field of G. Then [E : F ] G.
1. Galois Theory 1.1. A homomorphism of fields F F is simply a homomorphism of rings. Such a homomorphism is always injective, because its kernel is a proper ideal (it doesnt contain 1), which must therefore
More informationDIRECT SUM DECOMPOSABILITY OF POLYNOMIALS AND FACTORIZATION OF ASSOCIATED FORMS
DIRECT SUM DECOMPOSABILITY OF POLYNOMIALS AND FACTORIZATION OF ASSOCIATED FORMS MAKSYM FEDORCHUK Abstract. We prove two criteria for direct sum decomposability of homogeneous polynomials. For a homogeneous
More informationDuality of finite-dimensional vector spaces
CHAPTER I Duality of finite-dimensional vector spaces 1 Dual space Let E be a finite-dimensional vector space over a field K The vector space of linear maps E K is denoted by E, so E = L(E, K) This vector
More information3.1. Derivations. Let A be a commutative k-algebra. Let M be a left A-module. A derivation of A in M is a linear map D : A M such that
ALGEBRAIC GROUPS 33 3. Lie algebras Now we introduce the Lie algebra of an algebraic group. First, we need to do some more algebraic geometry to understand the tangent space to an algebraic variety at
More informationMAKSYM FEDORCHUK. n ) = z1 d 1 zn d 1.
DIRECT SUM DECOMPOSABILITY OF SMOOTH POLYNOMIALS AND FACTORIZATION OF ASSOCIATED FORMS MAKSYM FEDORCHUK Abstract. We prove an if-and-only-if criterion for direct sum decomposability of a smooth homogeneous
More informationMATH 101B: ALGEBRA II PART A: HOMOLOGICAL ALGEBRA
MATH 101B: ALGEBRA II PART A: HOMOLOGICAL ALGEBRA These are notes for our first unit on the algebraic side of homological algebra. While this is the last topic (Chap XX) in the book, it makes sense to
More informationChapter 4 Euclid Space
Chapter 4 Euclid Space Inner Product Spaces Definition.. Let V be a real vector space over IR. A real inner product on V is a real valued function on V V, denoted by (, ), which satisfies () (x, y) = (y,
More informationA Little Beyond: Linear Algebra
A Little Beyond: Linear Algebra Akshay Tiwary March 6, 2016 Any suggestions, questions and remarks are welcome! 1 A little extra Linear Algebra 1. Show that any set of non-zero polynomials in [x], no two
More informationVector Spaces and Linear Transformations
Vector Spaces and Linear Transformations Wei Shi, Jinan University 2017.11.1 1 / 18 Definition (Field) A field F = {F, +, } is an algebraic structure formed by a set F, and closed under binary operations
More informationChapter 2 The Approach of Moments for Polynomial Equations
Chapter 2 The Approach of Moments for Polynomial Equations Monique Laurent and Philipp Rostalski 2.1 Introduction Computing all points x K n (K = R or C) at which a given system of polynomials in n variables
More informationHilbert Polynomials. dimension and counting monomials. a Gröbner basis for I reduces to in > (I)
Hilbert Polynomials 1 Monomial Ideals dimension and counting monomials 2 The Dimension of a Variety a Gröbner basis for I reduces to in > (I) 3 The Complexity of Gröbner Bases a bound on the degrees of
More informationBare-bones outline of eigenvalue theory and the Jordan canonical form
Bare-bones outline of eigenvalue theory and the Jordan canonical form April 3, 2007 N.B.: You should also consult the text/class notes for worked examples. Let F be a field, let V be a finite-dimensional
More informationChapter 3. Rings. The basic commutative rings in mathematics are the integers Z, the. Examples
Chapter 3 Rings Rings are additive abelian groups with a second operation called multiplication. The connection between the two operations is provided by the distributive law. Assuming the results of Chapter
More informationReview 1 Math 321: Linear Algebra Spring 2010
Department of Mathematics and Statistics University of New Mexico Review 1 Math 321: Linear Algebra Spring 2010 This is a review for Midterm 1 that will be on Thursday March 11th, 2010. The main topics
More informationDunkl operators and Clifford algebras II
Translation operator for the Clifford Research Group Department of Mathematical Analysis Ghent University Hong Kong, March, 2011 Translation operator for the Hermite polynomials Translation operator for
More informationLiouvillian solutions of third order differential equations
Article Submitted to Journal of Symbolic Computation Liouvillian solutions of third order differential equations Felix Ulmer IRMAR, Université de Rennes, 0 Rennes Cedex, France felix.ulmer@univ-rennes.fr
More informationTracking multiplicities
Tracking multiplicities Jean-Claude Yakoubsohn Institut de Mathématiques de Toulouse Santalo s Summer School Recent Advances in Real Complexity and Computation Work in collaboration with Marc Giusti. Jean-Claude
More informationCounting Zeros over Finite Fields with Gröbner Bases
Counting Zeros over Finite Fields with Gröbner Bases Sicun Gao May 17, 2009 Contents 1 Introduction 2 2 Finite Fields, Nullstellensatz and Gröbner Bases 5 2.1 Ideals, Varieties and Finite Fields........................
More informationInfinite-Dimensional Triangularization
Infinite-Dimensional Triangularization Zachary Mesyan March 11, 2018 Abstract The goal of this paper is to generalize the theory of triangularizing matrices to linear transformations of an arbitrary vector
More informationNUMERICAL MONOIDS (I)
Seminar Series in Mathematics: Algebra 2003, 1 8 NUMERICAL MONOIDS (I) Introduction The numerical monoids are simple to define and naturally appear in various parts of mathematics, e.g. as the values monoids
More informationNovember 20, Interpolation, Extrapolation & Polynomial Approximation
Interpolation, Extrapolation & Polynomial Approximation November 20, 2016 Introduction In many cases we know the values of a function f (x) at a set of points x 1, x 2,..., x N, but we don t have the analytic
More informationON THE JORDAN DECOMPOSITION OF TENSORED MATRICES OF JORDAN CANONICAL FORMS
Math J Okayama Univ 51 (2009), 133 148 ON THE JORDAN DECOMPOSITION OF TENSORED MATRICES OF JORDAN CANONICAL FORMS Kei-ichiro IIMA and Ryo IWAMATSU Abstract Let k be an algebraically closed field of characteristic
More informationSPRING 2006 PRELIMINARY EXAMINATION SOLUTIONS
SPRING 006 PRELIMINARY EXAMINATION SOLUTIONS 1A. Let G be the subgroup of the free abelian group Z 4 consisting of all integer vectors (x, y, z, w) such that x + 3y + 5z + 7w = 0. (a) Determine a linearly
More informationThe Waring rank of the Vandermonde determinant
The Waring rank of the Vandermonde determinant Alexander Woo (U. Idaho) joint work with Zach Teitler(Boise State) SIAM Conference on Applied Algebraic Geometry, August 3, 2014 Waring rank Given a polynomial
More informationINTRODUCTION TO LIE ALGEBRAS. LECTURE 10.
INTRODUCTION TO LIE ALGEBRAS. LECTURE 10. 10. Jordan decomposition: theme with variations 10.1. Recall that f End(V ) is semisimple if f is diagonalizable (over the algebraic closure of the base field).
More informationRings and groups. Ya. Sysak
Rings and groups. Ya. Sysak 1 Noetherian rings Let R be a ring. A (right) R -module M is called noetherian if it satisfies the maximum condition for its submodules. In other words, if M 1... M i M i+1...
More informationOn the Grothendieck Ring of a Hopf Algebra
On the Grothendieck Ring of a Hopf Algebra H is a fin. dim l Hopf algebra over an alg. closed field k of char p 0, with augmentation ε, antipode S,.... (1) Grothendieck ring and character algebra The Grothendieck
More informationJournal of Algebra 226, (2000) doi: /jabr , available online at on. Artin Level Modules.
Journal of Algebra 226, 361 374 (2000) doi:10.1006/jabr.1999.8185, available online at http://www.idealibrary.com on Artin Level Modules Mats Boij Department of Mathematics, KTH, S 100 44 Stockholm, Sweden
More informationne varieties (continued)
Chapter 2 A ne varieties (continued) 2.1 Products For some problems its not very natural to restrict to irreducible varieties. So we broaden the previous story. Given an a ne algebraic set X A n k, we
More informationMath 113 Final Exam: Solutions
Math 113 Final Exam: Solutions Thursday, June 11, 2013, 3.30-6.30pm. 1. (25 points total) Let P 2 (R) denote the real vector space of polynomials of degree 2. Consider the following inner product on P
More informationMath 603, Spring 2003, HW 6, due 4/21/2003
Math 603, Spring 2003, HW 6, due 4/21/2003 Part A AI) If k is a field and f k[t ], suppose f has degree n and has n distinct roots α 1,..., α n in some extension of k. Write Ω = k(α 1,..., α n ) for the
More informationAlgebraic varieties and schemes over any scheme. Non singular varieties
Algebraic varieties and schemes over any scheme. Non singular varieties Trang June 16, 2010 1 Lecture 1 Let k be a field and k[x 1,..., x n ] the polynomial ring with coefficients in k. Then we have two
More informationSymmetric Tensor Decomposition and Algorithms
Symmetric Tensor Decomposition and Algorithms Author: María López Quijorna Directors: M.E.Alonso and A.Díaz-Cano Universidad Complutense de Madrid Academic Year: 2010-2011 Final Project of Master in Mathematic
More informationFormal power series rings, inverse limits, and I-adic completions of rings
Formal power series rings, inverse limits, and I-adic completions of rings Formal semigroup rings and formal power series rings We next want to explore the notion of a (formal) power series ring in finitely
More informationDe Nugis Groebnerialium 5: Noether, Macaulay, Jordan
De Nugis Groebnerialium 5: Noether, Macaulay, Jordan ACA 2018 Santiago de Compostela Re-reading Macaulay References Macaulay F. S., On the Resolution of a given Modular System into Primary Systems including
More informationREAL ANALYSIS II HOMEWORK 3. Conway, Page 49
REAL ANALYSIS II HOMEWORK 3 CİHAN BAHRAN Conway, Page 49 3. Let K and k be as in Proposition 4.7 and suppose that k(x, y) k(y, x). Show that K is self-adjoint and if {µ n } are the eigenvalues of K, each
More informationReduction of linear systems based on Serre s theorem
Reduction of linear systems based on Serre s theorem Mohamed S. Boudellioua and Alban Quadrat Abstract. Within a module-theoretic approach, we study when a (multidimensional) linear system can be defined
More informationLECTURE 7. k=1 (, v k)u k. Moreover r
LECTURE 7 Finite rank operators Definition. T is said to be of rank r (r < ) if dim T(H) = r. The class of operators of rank r is denoted by K r and K := r K r. Theorem 1. T K r iff T K r. Proof. Let T
More informationMaths for Signals and Systems Linear Algebra in Engineering
Maths for Signals and Systems Linear Algebra in Engineering Lectures 13 15, Tuesday 8 th and Friday 11 th November 016 DR TANIA STATHAKI READER (ASSOCIATE PROFFESOR) IN SIGNAL PROCESSING IMPERIAL COLLEGE
More informationGRÖBNER BASES AND POLYNOMIAL EQUATIONS. 1. Introduction and preliminaries on Gróbner bases
GRÖBNER BASES AND POLYNOMIAL EQUATIONS J. K. VERMA 1. Introduction and preliminaries on Gróbner bases Let S = k[x 1, x 2,..., x n ] denote a polynomial ring over a field k where x 1, x 2,..., x n are indeterminates.
More informationPacific Journal of Mathematics
Pacific Journal of Mathematics GROUP ACTIONS ON POLYNOMIAL AND POWER SERIES RINGS Peter Symonds Volume 195 No. 1 September 2000 PACIFIC JOURNAL OF MATHEMATICS Vol. 195, No. 1, 2000 GROUP ACTIONS ON POLYNOMIAL
More informationGlobal polynomial optimization with Moments Matrices and Border Basis
Global polynomial optimization with Moments Matrices and Border Basis Marta Abril Bucero, Bernard Mourrain, Philippe Trebuchet Inria Méditerranée, Galaad team Juillet 15-19, 2013 M. Abril-Bucero, B. Mourrain,
More informationTopics in linear algebra
Chapter 6 Topics in linear algebra 6.1 Change of basis I want to remind you of one of the basic ideas in linear algebra: change of basis. Let F be a field, V and W be finite dimensional vector spaces over
More informationThe Cayley-Hamilton Theorem and the Jordan Decomposition
LECTURE 19 The Cayley-Hamilton Theorem and the Jordan Decomposition Let me begin by summarizing the main results of the last lecture Suppose T is a endomorphism of a vector space V Then T has a minimal
More information