File Exchange

image thumbnail

Fractional Matrix Powers with Frechet Derivatives and Condition Number Estimate

version 1.1 (14.3 KB) by

Computing matrix power A^p in complex/real arithmetic, with condition number and Frechet derivatives

2 Downloads

Updated

View License

Computes the p'th power A^p of the matrix A for arbitrary real -1<p<1 and A with no nonpositive real eigenvalues, by the Schur-Pade algorithm. It also computes the Frechet derivative of A^p in any direction E and estimates the condition number for computing the matrix power.

This submission contains two functions: powerm_pade_fre.m uses complex arithmetic;
powerm_pade_fre_real.m uses real arithmetic which is intended for the case where both A and E are real.

The codes can be called in the following ways (same for powerm_pade_fre_real.m):

  X = POWERM_PADE_FRE(A,P)
 [X,~,COND] = POWERM_PADE_FRE(A,P)
 [X,F] = POWERM_PADE_FRE(A,P,E)
 [X,F,COND] = POWERM_PADE_FRE(A,P,E)
 [X,F,COND,NSQ,M] = POWERM_PADE_FRE(A,P,E)

where X is A^p, F is the Frechet derivative at A in the direction E, COND is the condition number estimate, NSQ is the number of matrix square roots computed and M is the degree of the Pade approximant used in the algorithm.

Function TEST_GALLERY.M runs a simple test of the codes. Matrix Function Toolbox (MFT) must be installed. Obtain it from http://www.maths.manchester.ac.uk/~higham/mftoolbox

More details can be found in:

N. J. Higham and L. Lin,
An Improved Schur--Pade Algorithm for Fractional Powers of a Matrix and their Frechet Derivatives
MIMS Eprint 2013.1, January 2013, revised May 2013.
http://eprints.ma.man.ac.uk/1972/

Comments and Ratings (1)

Great stuff! Thanx a lot!

Updates

1.1

typos, tags

MATLAB Release
MATLAB 8.0 (R2012b)

Download apps, toolboxes, and other File Exchange content using Add-On Explorer in MATLAB.

» Watch video