horner

Parameters

P

多項式または有理数の行列

x

数値または多項式または有理数の配列

説明

多項式の変数sがxで 置換される時, 多項式または有理数行列P = P(s)を評価します:

horner(P,x)=P(x)

例 (双一次変換): P = P(s) が有理行列と仮定すると, 有理行列P((1+s)/(1-s))は horner(P,(1+s)/(1-s))により得られます.

指定した周波数で有理行列を評価する場合は, freqプリミティブを使用する方が良いでしょう.

例

//数値ベクトルに関する多項式の評価
P=poly(1:3,'x')
horner(P,[1 2 5])
horner(P,[1 2 5]+%i)

//有理数の評価
s=poly(0,'s');M=[s,1/s];
horner(M,1)
horner(M,%i)
horner(M,1/s)

//数値行列に関する多項式の評価
X= [1 2;3 4]
p=poly(1:3,'x','c')
m=horner(p, X)
1*X.^0+2*X.^1+3*X.^2

参照

• freq — 周波数応答
• repfreq — 周波数応答
• evstr — 式の評価