Hello, dear friend, you can consult us at any time if you have any questions, add WeChat: daixieit

MATH20602

Numerical Analysis 1

Problem Sheet 6

Part A

(6.1)    The purpose of this first problem is to review some linear algebra and matrix computation rules.

1. The transpose of a matrix A is denoted AT . Write down AT , where A =  1   2    A =    A = 10(0)0

2. To say that A is an m K n matrix means it has m rows and n columns. If 北 is an n K 1 vector and y is a 1 K m vector, then yA and A are well defined. Compute A北 and yA when

A =     =  .(、)     y = (10 10)

3. When we say that  à Rn , or  is an n-vector, we mean that it is a column vector. Within paragraphs, we often write = (1 3)T to define a column vector without breaking across lines of text. Let

A = 3(1)   4(2)     = (1 3)T

Compute A and 北T A.

4. In the previous example, A北T and 北A are not defined. Let

A =     =  .(、)     y = 1   2    z = 1   2   3   4

Which of the following are well dened?

A北     Az ﹐    zA ﹐    北T A ﹐    yA.

(6.2)    This problem deals with the computational complexity of matrix arithmetic, in the sense of the number of arithmetical operations required.

1. How many additions and multiplications are necessary to compute A北, where A is a 2 K 2 matrix and is a 2-vector?

2. How many if A is n K n, and is an n-vector?

3. How many to compute AB for two n K n matrices A and B ?

Part B

(6.3)    Consider the linear system

2|1 ì |2 + |3  = 1 ﹐

2|1 + 2|2 + 2|3  = 4 ﹐

ì|1 ì |2 + 2|3  =  ì 5.

Formulate the Jacobi and Gauss-Seidel iterations for this linear system.

(6.4)    Formulate the Gauss-Seidel iteration for nding the solution x  à R3  of the linear system

ì(2)1    2(ì)1    ì(0)1.(、) x =  1(0).(、) 

Perform two iterations of the Gauss-Seidel method, with initial guess x0  = (0 1 0)T . (6.5)    Compute |x|o , |x|1 , and |x|2 for x = (1 ì2)T and x = (ì1 2 ì3)T .

(6.6)    Show that

|x| :=   max  |i

is not a vector norm on Rn .