Linear algebra -- compute the following without finding invA

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
4 replies · 2K views
SchrodingersMu
Messages
14
Reaction score
0

Homework Statement



upload_2015-3-6_12-4-18.png

Homework Equations


A=LU, U^-1 * L^-1= A^-1 , U^-1 * L^-1 * U^-1 * L^-1 = A^-2,

The Attempt at a Solution


I used MATLAB and the relations:
U^-1 * L^-1= A^-1 , U^-1 * L^-1 * U^-1 * L^-1 = A^-2,
to find a solution

I found U^-1*L^-1 , let =B

Then, found B^2 and took the inverse of B to get A^-2.

B=A^-1

C=inv(B^2)

so we have B*x+C*y= [2;5;10]

My question is: is there an easier way to do this? We are supposed to do this problem by hand. This would have taken a good amount fo effort to find so many inverse matrices. Any input is appreciated!

Thank you!
 
Physics news on Phys.org
SchrodingersMu said:

Homework Statement



View attachment 80035

Homework Equations


A=LU, U^-1 * L^-1= A^-1 , U^-1 * L^-1 * U^-1 * L^-1 = A^-2,

The Attempt at a Solution


I used MATLAB and the relations:
U^-1 * L^-1= A^-1 , U^-1 * L^-1 * U^-1 * L^-1 = A^-2,
to find a solution

I found U^-1*L^-1 , let =B

Then, found B^2 and took the inverse of B to get A^-2.

B=A^-1

C=inv(B^2)

so we have B*x+C*y= [2;5;10]

My question is: is there an easier way to do this? We are supposed to do this problem by hand. This would have taken a good amount fo effort to find so many inverse matrices.Any input is appreciated!

Thank you!
L and U are triangular matrices, so I don't think it would have been all that hard to compute their inverses by hand, especially since they're only 3 x 3. Also, once you have A-1, just multiply it by itself to get (A-1)2. You don't have to do the product of four matrices you show.
 
SchrodingersMu said:

Homework Statement



View attachment 80035

Homework Equations


A=LU, U^-1 * L^-1= A^-1 , U^-1 * L^-1 * U^-1 * L^-1 = A^-2,

The Attempt at a Solution


I used MATLAB and the relations:
U^-1 * L^-1= A^-1 , U^-1 * L^-1 * U^-1 * L^-1 = A^-2,
to find a solution

I found U^-1*L^-1 , let =B

Then, found B^2 and took the inverse of B to get A^-2.

B=A^-1

C=inv(B^2)

so we have B*x+C*y= [2;5;10]

My question is: is there an easier way to do this? We are supposed to do this problem by hand. This would have taken a good amount fo effort to find so many inverse matrices.

Any input is appreciated!

Thank you!
When you found matrices, B and C, that violated your instruction which was that you should not form A-1 and A-2 explicitly.

Simply find L-1 and U-1 and use those. They are easy to find by hand.
 
As you say, since A= LU, [itex]A^{-1}= L^{-1}U^{-1}[/itex] and [itex]A^{-2}= L^{-1}U^{-1}L^{-1}U^{-1}[/itex] so the first thing I would do is find [itex]L^{-1}[/itex] and [itex]U^{-1}[/itex]. Since L and U are triangular, that should be easy to do "by hand".

(Well, that was three answer in quick succession- and we are all saying the same thing.)
 
SchrodingersMu said:

Homework Statement



View attachment 80035

Homework Equations


A=LU, U^-1 * L^-1= A^-1 , U^-1 * L^-1 * U^-1 * L^-1 = A^-2,

The Attempt at a Solution


I used MATLAB and the relations:
U^-1 * L^-1= A^-1 , U^-1 * L^-1 * U^-1 * L^-1 = A^-2,
to find a solution

I found U^-1*L^-1 , let =B

Then, found B^2 and took the inverse of B to get A^-2.

B=A^-1

C=inv(B^2)

so we have B*x+C*y= [2;5;10]

My question is: is there an easier way to do this? We are supposed to do this problem by hand. This would have taken a good amount fo effort to find so many inverse matrices.Any input is appreciated!

Thank you!

You can do it by solving equations with ##A## on the left, which is doable easily from the LU decomposition. You want to compute ##v+w##, where ##v = A^{-1} x## and ##w = A^{-2} y##. We can get ##v## as the solution of the system ##A v = x## and we can get ##w## as the solution of ##A^2 w = y##, and this, in turn, can be obtained from the systems ##A w_0 = y##, ##A w = w_0##.

I don't know if this is what the questioner really wants.