This preview has intentionally blurred sections. Sign up to view the full version.View Full Document
Unformatted text preview: Program 1 Foundations of Computational Math 1 Fall 2010 Due date: via email by 11:59PM on Monday, 4 October General Task Consider the matrix A ∈ R n × k with linearly independent columns and the transformation H k H k- 1 ··· H 1 A = R where H i , 1 ≤ i ≤ k are Householder reflectors and R ∈ R k × k is a nonsingular upper triangular matrix. Implement a code that is capable of performing this transformation and storing the information about the H i and R efficiently in-place in a 2-dimensional n × k array and a small number of additional 1-dimensional arrays of length n or k . Also implement the necessary additional routines to solve the linear least squares problem min x ∈ R k k b- Ax k 2 given b ∈ R n . You must demonstrate your code on multiple examples with a wide range of values of n , k , and b for three situations: 1. n = k , i.e., a square nonsingular matrix A where x min = A- 1 b . 2. n > k and Ax = b for b ∈ R n and b ∈ R ( A ) i.e., a rectangular matrix A with full column rank and a vector b that define a consistent set of overdetermined equations. 3. n > k and b ∈ R n and b 6∈ R ( A ) i.e., a rectangular matrix A with full column rank and a vector b that define a linear least squares problem with a nonzero residual....
View Full Document
This note was uploaded on 07/25/2011 for the course MAD 5403 taught by Professor Gallivan during the Spring '11 term at University of Florida.
- Spring '11