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 = H.num; D = H.den; 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) // 複素数の場合
参照
| Report an issue | ||
| << polyint | Polynomials | roots >> |