Nom

nlev — Algorithme de Leverrier

Séquence d'appel

[num,den]=nlev(A,z [,rmax])

Paramètres

A

matrice réelle carrée

z

chaîne de caractères

rmax

paramètre optionnel (voir bdiag)

Description

[num,den]=nlev(A,z [,rmax]) calcule (z*eye()-A)^(-1) par bloc diagonalisation de A suivie de l'algorithme de Leverrier sur chaque bloc.

Cet algorithme est plus efficace que l'algorithme de Leverrier seul, mais il n'est toujours pas parfait !

Exemples

 
A=rand(3,3);x=poly(0,'x');
[NUM,den]=nlev(A,'x')
clean(den-poly(A,'x'))
clean(NUM/den-inv(x*eye()-A))

Voir Aussi

coff , coffg , glever , ss2tf

Auteurs

F. Delebecque, INRIA;
S. Steer, INRIA;