lcmdiag

Paramètres

H

matrice de polynômes

N

matrice de polynômes

D

matrice diagonale de polynômes

flag

chaîne de caractères : 'row' ou 'col' (par défaut)

Description

[kN,D]=lcmdiag(H,'row') calcule une factorisation D*H=N, i.e. H=D^(-1)*N où D est une matrice diagonale avec D(k,k)=ppcm de la k-ième ligne de H.den .

[N,D]=lcmdiag(H) ou [N,D]=lcmdiag(H,'col) renvoie H=N*D^(-1) avec D diagonale et D(k,k)=ppcm de ma k-ième colonne de H.den .

Exemples

s=poly(0,'s');
H=[1/s,(s+2)/s/(s+1)^2;1/(s^2*(s+2)),2/(s+2)];
[N,D]=lcmdiag(H);
N/D-H

Voir aussi

• lcm — Plus petit commun multiple (PPCM) de polynômes ou d'entiers
• gcd — calcul de PGCD
• bezout — équation de Bezout pour les polynômes