Code covered by the BSD License  

Highlights from
Fractional Matrix Powers with Frechet Derivatives and Condition Number Estimate

5.0 | 1 rating Rate this file 11 Downloads (last 30 days) File Size: 14.3 KB File ID: #41621 Version: 1.1

Fractional Matrix Powers with Frechet Derivatives and Condition Number Estimate



03 May 2013 (Updated )

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

| Watch this File

File Information

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):


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

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.

Required Products MATLAB
MATLAB release MATLAB 8.0 (R2012b)
Tags for This File   Please login to tag files.
Please login to add a comment or rating.
Comments and Ratings (1)
16 Sep 2013 Puppet Master Ivaschenko

Great stuff! Thanx a lot!

06 May 2013 1.1

typos, tags

Contact us