Code covered by the BSD License

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

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

# Fractional Matrix Powers with Frechet Derivatives and Condition Number Estimate

### Lijing Lin (view profile)

03 May 2013 (Updated )

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

File Information
Description

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 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/

Required Products MATLAB
MATLAB release MATLAB 8.0 (R2012b)
16 Sep 2013 Puppet Master Ivaschenko

### Puppet Master Ivaschenko (view profile)

Great stuff! Thanx a lot!