# residu

残差

## 呼び出し手順

    [V]=residu(P,Q1,Q2)

## パラメータ

P, Q1, Q2

  実数または複素数の係数を有する多項式または多項式の行列.

## 説明

`V=residu(P,Q1,Q2)` は, `V(i,j)`が `Q1(i,j)`のゼロで計算される 有理数`P(i,j)/(Q1(i,j)*Q2(i,j))` の残差の合計となるような行列`V`を返します.

`Q1(i,j)` および `Q2(i,j)` は, 共通の根を有さない必要があります.

## 例

    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))  //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    //正確
    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   //正確
    a=a+0*%i;b=b+0*%i;
    real(residu(b*btild,z*a,atild)-2.9488038) //複素数の場合

## 参照

- pfss
- bdiag
- roots
- poly
- gtild