28.12.2012 Views

Figure Properties - SERC

Figure Properties - SERC

Figure Properties - SERC

SHOW MORE
SHOW LESS

Create successful ePaper yourself

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

toeppen—Pentadiagonal Toeplitz matrix (sparse)<br />

gallery<br />

P = gallery('toeppen',n,a,b,c,d,e) returns the n-by-n sparse,<br />

pentadiagonal Toeplitz matrix with the diagonals: P(3,1) = a, P(2,1) = b,<br />

P(1,1) = c, P(1,2) = d, and P(1,3) = e, where a, b, c, d, and e are scalars.<br />

By default, (a,b,c,d,e) = (1,-10,0,10,1), yielding a matrix of Rutishauser.<br />

This matrix has eigenvalues lying approximately on the line segment<br />

2*cos(2*t) + 20*i*sin(t).<br />

tridiag—Tridiagonal matrix (sparse)<br />

A = gallery('tridiag',c,d,e) returns the tridiagonal matrix with<br />

subdiagonal c, diagonal d, and superdiagonal e. Vectors c and e must have<br />

length(d)-1.<br />

A = gallery('tridiag',n,c,d,e), where c, d, and e are all scalars, yields the<br />

Toeplitz tridiagonal matrix of order n with subdiagonal elements c, diagonal<br />

elements d, and superdiagonal elements e. This matrix has eigenvalues<br />

d + 2*sqrt(c*e)*cos(k*pi/(n+1))<br />

where k = 1:n. (see [1].)<br />

A = gallery('tridiag',n) is the same as<br />

A = gallery('tridiag',n,-1,2,-1), which is a symmetric positive definite<br />

M-matrix (the negative of the second difference matrix).<br />

triw—Upper triangular matrix discussed by Wilkinson and others<br />

A = gallery('triw',n,alpha,k) returns the upper triangular matrix with<br />

ones on the diagonal and alphas on the first k >= 0 superdiagonals.<br />

Order n may be a 2-element vector, in which case the matrix is n(1)-by-n(2)<br />

and upper trapezoidal.<br />

Ostrowski [“On the Spectrum of a One-parametric Family of Matrices, J. Reine<br />

Angew. Math., 1954] shows that<br />

cond(gallery('triw',n,2)) = cot(pi/(4*n))^2,<br />

2-947

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

Saved successfully!

Ooh no, something went wrong!