QR factorization of a n x 1 matrix

Join the discussion
Ask a follow-up here, or get your own question answered by working scientists, mathematicians and engineers — people, not an autocomplete.
Real named experts · corrections over time · the nuance an AI answer skips
2 replies · 3K views
abajaj2280
Messages
2
Reaction score
0

Homework Statement


Consider the vector a as an n × 1 matrix.

A) Write out its reduced QR factorization, showing the matrices [itex]\hat{Q}[/itex] and [itex]\hat{R}[/itex] explicitly.

B) What is the solution to the linear least squares problem ax ≃ b where b is a given n-vector.


Homework Equations


I was using the equation from 1.1 (http://www.math.ucla.edu/~yanovsky/Teaching/Math151B/handouts/GramSchmidt.pdf) to help me solve this problem.

The Attempt at a Solution


I haven't taken linear algebra for about 2 years and this is kind of hazy. I'm really confused here, and I really don't know where to start. I know that I'm supposed to come to this website with some sort of progress, but I'm really confused here.
 
Physics news on Phys.org
So, the first step in the Gram-Schmidt process is to think of the matrix as being a row of column vectors. Since your matrix is n x 1, it's like having a matrix of only one column vector. Thus [itex]A = [\mathbf{a_1}][/itex] in this case. So, going by the pdf you provided, let [itex]\mathbf{u_1} = \mathbf{a_1}[/itex] and then
[itex]\mathbf{e_1} = \frac{\mathbf{u_1}}{\left\| \mathbf{u_1} \right\|}[/itex] [itex]= \frac{\mathbf{a_1}}{\sqrt{(a_{11})^2+(a_{21})^2+...+(a_{n1}^2)}}[/itex].

In this case, [itex]Q = [\frac{ \mathbf{a_1} }{\left\| \mathbf{a_1} \right\|}][/itex] and R is the 1x1 matrix
[itex][ \mathbf{a_1} \bullet (\frac{ \mathbf{a_1} }{\left\| \mathbf{a_1} \right\|}) ][/itex] = [itex][ \frac{ \mathbf{a_1} \bullet \mathbf{a_1} }{ \left\| \mathbf{a_1} \right\| } ][/itex] = [itex][\frac{(\left\| \mathbf{a_1} \right\|)^2}{\left\| \mathbf{a_1} \right\|}] = [\left\|\mathbf{a_1}\right\|][/itex].

Then, [itex]QR = [ \frac{\mathbf{a_1}}{\left\|\mathbf{a_1}\right\|} \left\|\mathbf{a_1}\right\| ] = [\mathbf{a_1}][/itex].

In this case, the result isn't very interesting, because it's an nx1 matrix. But I think that's also to (supposedly) make it easier for you. I can see how in this case it made it even more confusing.
 
Thank you so much for your help. I have an exam in this class in one week, so I will be referring back to this when studying.