Mercurial > matrix-functions

diff mftoolbox/funm_ev.m @ 0:8f23314345f4 draft

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Create local repository for matrix toolboxes. Step #0 done.
author Antonio Pino Robles <data.script93@gmail.com>
date Wed, 06 May 2015 14:56:53 +0200
parents
children
line wrap: on
line diff
--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/mftoolbox/funm_ev.m	Wed May 06 14:56:53 2015 +0200
@@ -0,0 +1,9 @@
+function F = funm_ev(A,fun)
+%FUNM_EV   Evaluate general matrix function via eigensystem.
+%   F = FUNM_EV(A,FUN) evaluates the function FUN at the
+%   square matrix A using the eigensystem of A.
+%   This function is intended for diagonalizable matrices only
+%   and can be numerically unstable.
+
+[V,D] = eig(A);
+F = V * diag(feval(fun,diag(D))) / V;