15.12.2012 Views

scipy tutorial - Baustatik-Info-Server

scipy tutorial - Baustatik-Info-Server

scipy tutorial - Baustatik-Info-Server

SHOW MORE
SHOW LESS

You also want an ePaper? Increase the reach of your titles

YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.

SciPy Reference Guide, Release 0.8.dev<br />

Notes<br />

This is an interface to the LAPACK routines dgeqrf, zgeqrf, dorgqr, and zungqr.<br />

Examples<br />

>>> from <strong>scipy</strong> import random, linalg, dot<br />

>>> a = random.randn(9, 6)<br />

>>> q, r = linalg.qr(a)<br />

>>> allclose(a, dot(q, r))<br />

True<br />

>>> q.shape, r.shape<br />

((9, 9), (9, 6))<br />

>>> r2 = linalg.qr(a, mode=’r’)<br />

>>> allclose(r, r2)<br />

>>> q3, r3 = linalg.qr(a, econ=True)<br />

>>> q3.shape, r3.shape<br />

((9, 6), (6, 6))<br />

schur(a, output=’real’, lwork=None, overwrite_a=0)<br />

Compute Schur decomposition of a matrix.<br />

The Schur decomposition is<br />

A = Z T Z^H<br />

where Z is unitary and T is either upper-triangular, or for real Schur decomposition (output=’real’), quasi-upper<br />

triangular. In the quasi-triangular form, 2x2 blocks describing complex-valued eigenvalue pairs may extrude<br />

from the diagonal.<br />

See Also:<br />

Parameters<br />

a : array, shape (M, M)<br />

Matrix to decompose<br />

output : {‘real’, ‘complex’}<br />

Construct the real or complex Schur decomposition (for real matrices).<br />

lwork : integer<br />

Work array size. If None or -1, it is automatically computed.<br />

overwrite_a : boolean<br />

Whether to overwrite data in a (may improve performance)<br />

Returns<br />

T : array, shape (M, M)<br />

Schur form of A. It is real-valued for the real Schur decomposition.<br />

Z : array, shape (M, M)<br />

An unitary Schur transformation matrix for A. It is real-valued for the real Schur<br />

decomposition.<br />

rsf2csf<br />

Convert real Schur form to complex Schur form<br />

240 Chapter 3. Reference

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

Saved successfully!

Ooh no, something went wrong!