15.11.2014 Views

MATLAB Mathematics - SERC - Index of

MATLAB Mathematics - SERC - Index of

MATLAB Mathematics - SERC - Index of

SHOW MORE
SHOW LESS

Create successful ePaper yourself

Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.

1 Matrices and Linear Algebra<br />

The elements <strong>of</strong> A are binomial coefficients. Each element is the sum <strong>of</strong> its<br />

north and west neighbors. The Cholesky factorization is<br />

R = chol(A)<br />

R =<br />

1 1 1 1 1 1<br />

0 1 2 3 4 5<br />

0 0 1 3 6 10<br />

0 0 0 1 4 10<br />

0 0 0 0 1 5<br />

0 0 0 0 0 1<br />

The elements are again binomial coefficients. The fact that R'*R is equal to A<br />

demonstrates an identity involving sums <strong>of</strong> products <strong>of</strong> binomial coefficients.<br />

Note The Cholesky factorization also applies to complex matrices. Any<br />

complex matrix which has a Cholesky factorization satisfies A' = A and is said<br />

to be Hermitian positive definite.<br />

The Cholesky factorization allows the linear system<br />

Ax<br />

=<br />

b<br />

to be replaced by<br />

R′Rx = b<br />

Because the backslash operator recognizes triangular systems, this can be<br />

solved in <strong>MATLAB</strong> quickly with<br />

x = R\(R'\b)<br />

If A is n-by-n, the computational complexity <strong>of</strong> chol(A) is O(n 3 ), but the<br />

complexity <strong>of</strong> the subsequent backslash solutions is only O(n 2 ).<br />

1-28

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!