第9回: 課題ファイル#
課題1#
課題1-1: 三角関数系の直交性をsympyで検証し,手計算の結果と一致することを確認せよ.ただし,すべての組み合わせを確認することはできないので,\(1, \cos x, \cos 2x \cos 3x, \sin x, \sin 2x, \sin 3x\) の中で検証せよ.
課題1-2: 次の関数 \(P_n(x), n=0,1,2,\ldots,\) は \(n\)次のルジャンドルの多項式と呼ばれる.
\[\begin{split}
\begin{align}
P_0(x) &= 1 \\
P_1(x) &= x \\
P_2(x) &= \frac{1}{2}(3 x^{2} - 1) \\
P_3(x) &= \frac{1}{2}(5 x^{3} - 3x) \\
P_4(x) &= \frac{1}{8}(35x^{4}-30x^{2}+3) \\
P_5(x) &= \frac{1}{8}(63x^{5}-70x^{3}+15x) \\
P_6(x) &= \frac{1}{16}(231x^{6}-315x^{4}+105x^{2}-5) \\
P_7(x) &= \frac{1}{16}(429x^{7}-693x^{5}+315x^{3}-35x) \\
P_8(x) &= \frac{1}{128}(6435x^{8}-12012x^{6}+6930x^{4}-1260x^{2}+35) \\
\vdots
\end{align}
\end{split}\]
そして,ルジャンドル多項式は区間 \([-1,1]\) 上で直交関数系である.
\[\begin{split}
\int^{1}_{-1} P_n(x)P_m(x)dx =\left\{\begin{matrix}
\frac{2}{2n+1} \quad m = n \\
0 \quad m \neq n
\end{matrix}\right.
\end{split}\]
ルジャンドル多項式が直交関数系であることをsympyを使って確認せよ(ヒント:区間 \([-1,1]\) で積分して自身以外との内積が0になることを示せば良い).課題1-1と同様にすべての組み合わせを検証することはできないので一部の項のみで良い.
課題2#
課題2-1: 区間 \([-\pi, \pi]\) 上の以下の関数 \(f(t)\) をフーリエ級数展開せよ.
\[
f(t) = t^2
\]
課題2-2: 区間 \([-\pi, \pi]\) 上の以下の関数 \(f(t)\) をフーリエ級数展開せよ.
\[
f(t) = \cos t \sin t
\]
課題3#
課題3-1: ホームページのサンプルコードを参考に,近似対象の関数やフーリエ級数の展開の項数を変更して結果を確認せよ.