Hi everyone,
I have a question about how one might use Mathematica 3.0.x to perform
general matrix inverses for square matrices.
Specifically, suppose that I have a square complex matrix A having
dimensions (m+n)x(m+n). Let A be partitioned as (please bear with me in
the following ascii representation)
A = [ A11 A12
0 A22]
where A11, A12, A22 are complex matrices of dimensions (nxn), (nxm), and
(mxm), respectively. (The zero is a zero matrix having dimensions (mxn).)
Is there a way that I can represent the above generalization in
Mathematica without explicitly defining dimensions n or m? If so, how can
I get Mathematica to produce
A^(-1) = [A11^(-1) -A11^(-1)A12*A22^(-1)
0 A22^(-1)]
assuming, of course, that inverses for A11 and A22 exist?
Moreover, can Mathematica 3.0.x recognize and utilize many well known
(finite dimensional) matrix identities (e.g. matrix inversion lemma) for
general matrix simplification, where, generally, the size of individual
matrices is left unspecified (but all matrix dimensions are assumed to be
compatable/consistent), and the invertability of certain matrices is
assumed? (Since this is all symbolic, I am hoping that this is the case
...) Or is this something that Mathematica 4.x can do?
Many thanks in advance!
Mike