このページの内容は最新ではありません。最新版の英語を参照するには、ここをクリックします。

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 種チェビシェフ多項式をプロットします。

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')

入力引数

n — 多項式の次数

非負の整数 | シンボリック変数 | シンボリック式 | シンボリック関数 | ベクトル | 行列

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

x — 評価点

数値 | シンボリック数 | シンボリック変数 | シンボリック式 | シンボリック関数 | ベクトル | 行列

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

詳細

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

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

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

  $$T\left(0,x\right)=1,T\left(1,x\right)=x,T\left(n,x\right)=2xT\left(n-1,x\right)-T\left(n-2,x\right)$$

• 第 1 種チェビシェフ多項式は区間 -1 ≤ x ≤ 1 において重み関数 $$w\left(x\right)=\frac{1}{\sqrt{1-x^2}}$$ と直交します。

  $${\displaystyle \underset{-1}{\overset{1}{\int }}\frac{T(n,x)T(m,x)}{\sqrt{1-x^2}}dx}=\{\begin{array}{cc}0& \text{if }n\ne m\\ \pi & \text{if }n=m=0\\ \frac{\pi }{2}& \text{if }n=m\ne 0.\end{array}$$

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

  $$T\left(n,x\right)=\frac{2^{2n}{\left(n!\right)}^2}{\left(2n\right)!}P\left(n,-\frac{1}{2},-\frac{1}{2},x\right)$$

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

  $$T\left(n,x\right)=\{\begin{array}{cc}\frac{1}{2}\underset{a\to 0}{\lim }\frac{n+a}{a}G\left(n,a,x\right)& \text{if }n\ne 0\\ \underset{a\to 0}{\lim }G\left(0,a,x\right)=1& \text{if }n=0\end{array}$$

ヒント

• 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.

[2] Cohl, Howard S., and Connor MacKenzie. “Generalizations and Specializations of Generating Functions for Jacobi, Gegenbauer, Chebyshev and Legendre Polynomials with Definite Integrals.” Journal of Classical Analysis, no. 1 (2013): 17–33. https://doi.org/10.7153/jca-03-02.

バージョン履歴

R2014b で導入

参考

chebyshevU | gegenbauerC | hermiteH | jacobiP | laguerreL | legendreP