Solving a System of Multiplicative Matrices

  • #1

Main Question or Discussion Point

Please do not be offended by my literary style. I find thinking about mathematical problems in such a way helps me learn better.

A is a 2x2 matrix of complex numbers, call this "apple"
B is a 2x2 matrix of complex numbers, call this "banana"

Let a "Fruit Salad" be defined as follows:
S = AB

Suppose that I had two apples and two bananas in my fruit bowl
A1 A2 B1 B2

However, before I was able to measure my apples/bananas, my wife cut them up and mixed them to make four different servings of fruit salad:
S11 = A1B1
S12 = A1B2
S21 = A2B1
S22 = A2B2

So now I have a plethora of fruit salad {S11, S12, S21 and S22}, but I really want to know about the fruit that was put into it {A1, A2, B1 and B2}. How can I solve for A1, A2, B1 and B2 given S11, S12, S21 and S22?

That's the end of my question. The remainder of this post is to demonstrate my thoughts.

If my apples, bananas and fruit salads were scalar quantities I could set up the following simultaneous equations to solve via Gaussian elimination:
ln(S11) = ln(A1) + ln(B1)
ln(S12) = ln(A1) + ln(B2)
ln(S21) = ln(A2) + ln(B1)
ln(S22) = ln(A2) + ln(B2)

I am aware of the concept of a "matrix exponential" and its inverse, the "matrix logarithm". However, heading down this path looks like it is going to be messy. I thought I would first check to see if anyone knows of a nicer way to solve this problem.
 
Last edited:

Answers and Replies

  • #2
UltrafastPED
Science Advisor
Gold Member
1,912
216
Shouldn't S21 = A1B2 actually be A2B1?
 
  • #3
You are correct UltrafastPED. I have modified the original post. My apologies.
 
  • #4
Office_Shredder
Staff Emeritus
Science Advisor
Gold Member
3,750
99
I don't think this will work in general.... eA+B = eAeB only if A and B commute, so I assume logarithm has the same property - that log(AB) = log(A) + log(B) only if A and B commute. You'll have to use the more complicated

http://en.wikipedia.org/wiki/Baker–Campbell–Hausdorff_formula

which doesn't really help since you don't get a set of linear equations
 
  • #5
You raise a good point Office_Shredder.

I don't think I that my original question contains enough information to yield a unique solution. In the particular problem I am looking to solve A and B do not always commute. However, suppose that the following is known:
T11 = B1A1
T12 = B1A2
T21 = B2A1
T22 = B2A2

(That is, my wife made eight different types of fruit salad with the available fruit.)
 
  • #6
SteamKing
Staff Emeritus
Science Advisor
Homework Helper
12,798
1,666
Your 'fruit bowel' is actually a 'fruit bowl'. The 'bowel' is another name for the lower digestive tract.
 
  • #7
Thank you Office_Shredder. Your comment about "eA+B = eAeB only if A and B commute" has highlighted something to me. My original question does not contain enough information to yield a unique solution. This is why an assumption (such as commutativity) is required for the mathematics to work out.

Suppose that
A = [ a11 a12 ; a21 a22 ]
B = [ b11 b12 ; b21 b22 ]
S = [ s11 s12 ; s21 s22 ]

In my original question I had assumed that I had provided 16 equations and 16 unknowns.

The unknowns were:
a11, a12, ...
b11, ...

The equations were of the the form:
a11b11 + a12b21 = s11
a11b12 + a12b22 = s12
etc.

I was wrong in assuming this. What I had really constructed was a set of 16 linear equations with 32 unknowns, because the unknowns were actually of the following form:
a11b11, a12b21, a11b12, a12b22, ...

If I were to introduce some extra information into the question - For example, suppose my wife had made the following fruit salads in addition to the ones mentioned earlier:
T11 = B1A1
T12 = B1A2
T21 = B2A1
T22 = B2A2

(Note that assuming commutativity is the same as assuming Txy = Sxy where x,y = {1,2})

With this information I can construct a system of 32 linear equations with 32 unknowns. Gaussian elimination will enable me to solve for these 32 unknowns, which are of the form:
a11b11, a12b21, a11b12, a12b22, ...

I can then take the logarithm of each of these to set up another system of equations:
ln(a11) + ln(b11) = ln(a11b11)
ln(a12) + ln(b21) = ln(a12b21)
...

Following from this, I can use Gaussian elimination to solve for my unknowns:
ln(a11)
ln(b11)
ln(a12)
ln(b21)

From these terms I can calculate the following:
a11, a12, ...
b11, ...

Hence solving for A1, A2, B1 and B2. (The components of my fruit salads.)

These are my thoughts so far. There may be an error somewhere in my line of reasoning, so I won't close the thread just yet. I shall go away to see if the extra information (i.e. T11, T12, T21 and T22) was what I needed to know to solve the problem.
 
Last edited:
  • #8
Your 'fruit bowel' is actually a 'fruit bowl'. The 'bowel' is another name for the lower digestive tract.
How embarrassing. The typographical error has been corrected.
 
  • #9
SteamKing
Staff Emeritus
Science Advisor
Homework Helper
12,798
1,666
Seems like a lot of work, finding matrix exponentials and logarithms. Why don't you just multiply the coefficient matrices together?

It's not clear what the entries of A, B, and S are supposed to represent. Slices of apples, bananas and the number mixed together ...?
 
  • #10
These are my thoughts so far. There may be an error somewhere in my line of reasoning, so I won't close the thread just yet. I shall go away to see if the extra information (i.e. T11, T12, T21 and T22) was what I needed to know to solve the problem.
Alas, no luck. I am not sure how to solve this one. For now I think I am going to resort to coding up an algorithm to find a numeric approximation.
 

Related Threads on Solving a System of Multiplicative Matrices

Replies
4
Views
500
Replies
6
Views
3K
  • Last Post
Replies
8
Views
4K
  • Last Post
Replies
4
Views
2K
  • Last Post
Replies
2
Views
393
  • Last Post
Replies
2
Views
5K
Replies
4
Views
10K
  • Last Post
Replies
3
Views
1K
Replies
2
Views
1K
  • Last Post
Replies
1
Views
1K
Top