Starting from:

$30

CSE6643 Homework 3 Solved

Notice: For the computational problems, you must design your own ways, such as using tables and plots etc, to present the results you obtain. And you must show the procedure you take, and explain the results you get, and draw conclusions if you have any. Computer code is not considered as an explanation. Please upload your code in T-square, and return your HW in class on or before the Due date.

•    Problem 1: Use inverse power iteration to compute the eigenvector of matrix

 

that corresponds to the eigenvalue close to µ = 1.268. You can use calculators, but not computers to help you. Compute at least three iterations and comment on the accuracy of your results.

•    Problem 2: Use power iteration to design an algorithm that computes the largest, in the absolute value sense, root of a polynomial p(z) = zn+a1zn−1+···+an−1z+an. Justify your answer.

•    Problem 3: Suppose A ∈ Rm×n with m > n and y ∈ Rm. Construct an augmented matrix A¯ = [A,y]. Show that the condition number, measured in the induced matrix 2-norm, of A¯ is larger than that of A.

•    Problem 4: Let A be a matrix in the following form

 a1

 c2



 c3 

 ...



 ...





cn
b2

a2

0

...

...

0
b3 0

a3

...

···
···

···

...

...

... ···
···

···

...

an−1

0
bn 

0 

... 

... .

 0 

an
Give an detailed algorithm using Givens rotation to compute the QR factorization of A. Write your own code to Implement your algorithm. Apply the code to matrices with randomly generated ai, bi and ci values. Test your code with n = 100,200,400,800,··· . Comment on the computational complexity of your tests.

•    Problem 5: Use the matrix A, with λ = 0 and λ = 2 respectively, from of Problem 5 in HW1 to finish the following work. You should run your code for matrix size m = 256,512,1024,···. You need comment on your results as well as computation cost.

(a)    Write your code for the power iteration or the inverse iteration to identify thelargest and smallest eigenvalues of A. How do you speed up the computation?

What is the condition number of A?

(b)    Write your code for QR algorithm to compute the eigenvalues and eigenvectors.

2

More products