This preview shows pages 1–3. Sign up to view the full content.
This preview has intentionally blurred sections. Sign up to view the full version.View Full Document
Unformatted text preview: Solutions for Homework 2 Numerical Linear Algebra 1 Fall 2011 Problem 2.1 Let x C n and y C n be two arbitrary vectors. Consider determining a circulant matrix C C n n such that y = Cx 2.1.a . Assume that C exists for a given pair ( x,y ), show how to construct it. 2.1.b . When is C unique for a given pair ( x,y )? 2.1.c . When does C not exist for a given pair ( x,y )? Solution: There are two equivalent lines of reasoning. The first is in the spatial domain. Denote C = C ( c ) to indicate that it is generated by its first row c T . We have Cx = y C ( c ) x = y H ( x ) c = y C ( x r ) c r = y where H ( x ) C n n . H ( x ) is not however circulant in the same fashion as C ( c ). It is a Hankel matrix not a Toeplitz matrix. A Hankel matrix is constant along its antidiagonals. It is therefore circulant with respect to the first row using a circulant shift to the left (not to the right as with C ( c )) and with respect to the first column using a circulant shift up (not down as with C ( c )). The Hankel/left circulant form of the matrix vector product can be transformed into the Toeplitz/right circulant form of the matrix vector product. This is done by using a Toeplitz/right circulant matrix and a righthand side vector x r and c r that are the vectors x and c with their elements in reverse order respectively, i.e., P T = ( e n e n- 1 ... e 1 ) , x r = Px, c r = Pc Since x and y are given, one way to find c is to solve the system with the matrix, H ( x ), defined by x . Therefore, to have a solution we must have y R ( H ( x )) or y R ( C ( x r )). 1 These facts are easily seen by example with n = 4. x = 1 2 3 c = 1 2 3 y = C ( c ) x = + 1 1 + 2 2 + 3 3 3 + 1 + 1 2 + 2 3 2 + 3 1 + 2 + 1 3 1 + 2 1 + 3 2 + 3 = + 1 1 + 2 2 + 3 3 1 + 1 2 + 2 3 + 3 2 + 1 3 + 2 + 3 1 3 + 1 + 2 1 + 3 2 = H ( ) c = 3 3 + 2 2 + 1 1 + 3 + 2 3 + 1 2 + 1 3 1 + 2 + 1 3 + 2 3...
View Full Document
This note was uploaded on 11/10/2011 for the course MAD 5932 taught by Professor Gallivan during the Fall '06 term at FSU.
- Fall '06