residu

Paramètres

P, Q1, Q2

polynômes ou matrice de polynômes à coefficients réels ou complexes.

Description

V=residu(P,Q1,Q2) renvoie la matrice V telle que V(i,j) est la somme des résidus de la fraction rationnelle P(i,j)/(Q1(i,j)*Q2(i,j)) calculés aux zéros de Q1(i,j).

Q1(i,j) et Q2(i,j) ne doivent pas avoir de racines communes.

Exemples

s=poly(0,'s');
H=[s/(s+1)^2,1/(s+2)];N=numer(H);D=denom(H);
w=residu(N.*horner(N,-s),D,horner(D,-s));  //N(s) N(-s) / D(s) D(-s)
sqrt(sum(w))  // La norme H2
h2norm(tf2ss(H))

p=(s-1)*(s+1)*(s+2)*(s+10);a=(s-5)*(s-1)*(s*s)*((s+1/2)**2);
b=(s-3)*(s+2/5)*(s+3);
residu(p,a,b)+531863/4410    // Exact
z=poly(0,'z');a=z^3+0.7*z^2+0.5*z-0.3;b=z^3+0.3*z^2+0.2*z+0.1;
atild=gtild(a,'d');btild=gtild(b,'d');
residu(b*btild,z*a,atild)-2.9488038   // Exact
a=a+0*%i;b=b+0*%i;
real(residu(b*btild,z*a,atild)-2.9488038) // Cas complexe

Voir aussi

• pfss — partial fraction decomposition
• bdiag — bloc-diagonalisation, vecteurs propres généralisés
• roots — racines d'un polynôme
• poly — définition d'un polynôme
• gtild — tilde operation

Auteurs

F Delebecque INRIA