Scribd
Upload
38 views2 pages

MATH 3650 Homework 8

Uploaded by

corymoody5
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as PDF, TXT or read online on Scribd
38 views2 pages

MATH 3650 Homework 8

Uploaded by

corymoody5
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as PDF, TXT or read online on Scribd
You are on page 1/ 2

Homework 08

MATH 3650 Spring 2024 ©Timo Heister, [email protected]


Note: Several of these questions require you to write code, let it produce output, and include
interpretation. My suggestion for this: 1. submit all your code as .m files. 2. Create a Word
document where you copy and paste the output and add your interpretation, then print to pdf
and attach.

1. (a) Compute (on paper) the eigenvalues, the corresponding set of all eigenvectors for each
eigenvalue, and one example eigenvector for each eigenvalue (pick it to have k · k∞ = 1
of the matrices
 
  0 8 0  
1 3 −3 0
A= , B = 3 −2 −3 ,
  C= .
0 2 −1 0
0 0 4

Note: please include all steps of the computation!


(b) Compute the characteristic polynomial for the following two matrices:
 
  1 0 1
0 −1
D= , E = 1 2 3
1 0
0 1 3

2. (Submit hw08q2.m , your modified power method, and a pdf with output and your interpre-
tations/answers) You are given the matrices:
   
0 2 3 2 3 2
A = 4 5 6  , B= 1 0 −2 .
7 8 −9 −1 −3 −1

(a) Test PowerMethod.m on the matrix A (pick a tolerance of 1e-8). You can choose the
starting vector. Use the eig function to verify that your answer is correct.
(b) Now test the matrix B with starting value [1;1;1] . Why does the power method fail
(check with eig )? Create a modified power method that prints the vector x in each
iteration (and aborts after 15 iterations). What happens and why does this happen?
(c) With the same B, try running with the starting vector x=[2 ; -2; 2] . What happens
now? Can you explain why?

3. Implement the function [lambda,x] = InvPowerMethod(A, x0, tolerance) that finds the
smallest eigenvalue of a matrix A using the Inverse Power Method without computing the
inverse of A (as described in class). Include an example of using InvPowerMethod by checking
with the matrix A in question 2.
submit on canvas: InvPowerMethod.m and hw08q3.m (that prints the smallest eigenvalue of
A).

1
4. Consider the population example in chapter 7.3 (also see canvas). Suppose you have birth
rates b1 = 0.3, b2 = 0.3, b3 = 0.3, b4 = 0.1 and death rates d1 = 0.1, d2 = 0.2, d3 = 0.5,
d4 = 0.9.

(a) What is the biggest eigenvalue of the resulting matrix? What do you expect to happen
with the population over a long time period?
(b) Suppose you have P1 = 100, P2 = 200, P3 = 150, P4 = 75 in year 0. What will the
population be in year 3000? Does this agree with a)?
(c) Suppose we change the problem by reducing the death rate of P4 to d4 = 0.02 (from
0.9). What is the largest eigenvalue now? What do you expect to happen after a long
time period? What is the population after 1000 years?

Submit your code hw08q4.m , the output, and your interpretations/answers.

You might also like