ドキュメンテーション

最新のリリースでは、このページがまだ翻訳されていません。 このページの最新版は英語でご覧になれます。

chebyshevT

第 1 種チェビシェフ多項式

構文

chebyshevT(n,x)

説明

chebyshevT(n,x) は、点 x における n 次の第 1 種チェビシェフ多項式を表します。

最初の 5 つの第 1 種チェビシェフ多項式

変数 x の最初の 5 つの第 1 種チェビシェフ多項式を求めます。

syms x
chebyshevT([0, 1, 2, 3, 4], x)
ans =
[ 1, x, 2*x^2 - 1, 4*x^3 - 3*x, 8*x^4 - 8*x^2 + 1]

数値引数およびシンボリック引数のチェビシェフ多項式

引数に応じて、chebyshevT は浮動小数点解またはシンボリック厳密解の結果を返します。

これらの点で 5 次の第 1 種チェビシェフ多項式の値を求めます。これらの数値はシンボリック オブジェクトではないため、chebyshevT は浮動小数点の結果を返します。

chebyshevT(5, [1/6, 1/4, 1/3, 1/2, 2/3, 3/4])
ans =
    0.7428    0.9531    0.9918    0.5000   -0.4856   -0.8906

シンボリック オブジェクトに変換された同じ数値に対する 5 次の第 1 種チェビシェフ多項式の値を求めます。シンボリック数に対し、chebyshevT はシンボリック厳密解の結果を返します。

chebyshevT(5, sym([1/6, 1/4, 1/3, 1/2, 2/3, 3/4]))
ans =
[ 361/486, 61/64, 241/243, 1/2, -118/243, -57/64]

浮動小数点数によるチェビシェフ多項式の求解

chebyshevT の直接呼び出しによるチェビシェフ多項式の浮動小数点求解は数値的に安定しています。しかし、最初にシンボリック変数を使用して多項式を計算し、それからこの式に可変精度の値を代入すると数値的に不安定になる場合があります。

1/3 および vpa(1/3) における 500 次の第 1 種チェビシェフ多項式の値を求めます。浮動小数点求解は数値的に安定しています。

chebyshevT(500, 1/3)
chebyshevT(500, vpa(1/3))
ans =
    0.9631
 
ans =
0.963114126817085233778571286718

続いて、シンボリック多項式 T500 = chebyshevT(500, x) を求め、結果に x = vpa(1/3) を代入します。この方法は数値的に不安定です。

syms x
T500 = chebyshevT(500, x);
subs(T500, x, vpa(1/3))
ans =
-3293905791337500897482813472768.0

vpa を使用して多項式の係数を近似し、結果に x = sym(1/3) を代入します。この方法も数値的に不安定です。

subs(vpa(T500), x, sym(1/3))
ans =
1202292431349342132757038366720.0

第 1 種チェビシェフ多項式のプロット

最初の 5 つの第 1 種チェビシェフ多項式をプロットします。R2016a より前の場合は、fplot の代わりに ezplot を使用します。

syms x y
fplot(chebyshevT(0:4, x))
axis([-1.5 1.5 -2 2])
grid on

ylabel('T_n(x)')
legend('T_0(x)', 'T_1(x)', 'T_2(x)', 'T_3(x)', 'T_4(x)', 'Location', 'Best')
title('Chebyshev polynomials of the first kind')

入力引数

すべて折りたたむ

多項式の次数。非負の整数、シンボリックな変数、式または関数として、あるいは数値、シンボリックな数値、変数、式または関数のベクトルまたは行列として指定します。

求解点。数値、シンボリックな数値、変数、式または関数、あるいは数値、シンボリックな数値、変数、式または関数のベクトルまたは行列として指定します。

詳細

すべて折りたたむ

第 1 種チェビシェフ多項式

第 1 種チェビシェフ多項式は Tn(x) = cos(n*arccos(x)) として定義されます。

これらの多項式は次の再帰式を満たします。

T(0,x)=1,T(1,x)=x,T(n,x)=2xT(n1,x)T(n2,x)

第 1 種チェビシェフ多項式は区間 -1 ≤ x ≤ 1 において次の重み関数と直交します。

w(x)=11x2

第 1 種チェビシェフ多項式は次のヤコビ多項式、

T(n,x)=22n(n!)2(2n)!P(n,12,12,x)

および次のゲーゲンバウアー多項式の特殊なケースです。

T(n,x)=n2G(n,0,x)

ヒント

  • chebyshevT は、シンボリック オブジェクトではない数値引数に対し浮動小数点の結果を返します。

  • chebyshevT は非スカラーの入力の要素ごとに働きます。

  • 少なくとも 1 つの入力引数はスカラーであるか、両方の引数は同じサイズのベクトルまたは行列でなければなりません。一方の入力引数がスカラーであり、もう一方の入力引数がベクトルまたは行列である場合、chebyshevT によってスカラーは、すべての要素がそのスカラーと等しい、もう一方の引数と同じサイズのベクトルまたは行列に拡張されます。

参照

[1] Hochstrasser, U. W. “Orthogonal Polynomials.” Handbook of Mathematical Functions with Formulas, Graphs, and Mathematical Tables. (M. Abramowitz and I. A. Stegun, eds.). New York: Dover, 1972.

R2014b で導入