MathGroup Archive 2000

[Date Index] [Thread Index] [Author Index]

Search the Archive

General matrix inverses

  • To: mathgroup at
  • Subject: [mg23516] General matrix inverses
  • From: bt585 at (Michael Chang)
  • Date: Tue, 16 May 2000 02:45:02 -0400 (EDT)
  • Organization: The National Capital FreeNet
  • Sender: owner-wri-mathgroup at

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!


  • Prev by Date: Re: plot discrete spectrum
  • Next by Date: Re: Problem with partial differential equations
  • Previous by thread: Responses for Box-Cox transformations
  • Next by thread: [Question] How to load(get) Packages