バックナンバーはこちら。
https://www.simulationroom999.com/blog/compare-matlabpythonscilabjulia5-backnumber/
はじめに
フーリエ係数に至る道。
今回は三角関数の直交性をScilabで確認してみる。
登場人物
博識フクロウのフクさん
イラストACにて公開の「kino_k」さんのイラストを使用しています。
https://www.ac-illust.com/main/profile.php?id=iKciwKA9&area=1
エンジニア歴8年の太郎くん
イラストACにて公開の「しのみ」さんのイラストを使用しています。
https://www.ac-illust.com/main/profile.php?id=uCKphAW2&area=1
【再掲】フーリエ係数に至る道
太郎くん
まずは、フーリエ係数に至る道を再掲。
- 偶関数
- 奇関数
- 関数の内積
- 三角関数の加法定理
- 三角関数の積和公式
- 重要な極限値
- 三角関数の直交性
- フーリエ係数
フクさん
今回は三角関数の直交性をScilabで確認してみる。
【再掲】プログラムで確認する内容
太郎くん
以下をプログラムで確認。
- \(\sin(x)\cdot\cos(x)=0\)
- \(\sin(x)\cdot\cos(2x)=0\)
- \(\sin(x)\cdot\sin(x)=\pi\)
- \(\cos(2x)\cdot\cos(2x)=\pi\)
- \(\cos(x)\cdot\sin(2x)=0\)
- \(\sin(x)\cdot\cos(2x)=0\)
フクさん
畳み込みを意識してみると良いだろう。
Scilabコード
フクさん
Scilabコードは以下となる。
N = 1000000; // 要素数
L = %pi; // 0を中心とした±幅
x = linspace(-L,L,N); // x軸
dx = 2*L/N; // Δx
y=sin(x)*cos(x)'*dx;
printf('sin(x)・cos(x)=%.5f\n', y)
y=sin(x)*cos(2*x)'*dx;
printf('sin(x)・cos(2x)=%.5f\n', y)
y=sin(x)*sin(x)'*dx;
printf('sin(x)・sin(x)=%.5f\n', y)
y=cos(2*x)*cos(2*x)'*dx;
printf('cos(2x)・cos(2x)=%.5f\n', y)
y=cos(x)*sin(2*x)'*dx;
printf('cos(x)・sin(2x)=%.5f\n', y)
y=sin(x)*cos(2*x)'*dx;
printf('sin(x)・cos(2x)=%.5f\n', y)処理結果
フクさん
処理結果は以下。
sin(x)・cos(x)=0.00000
sin(x)・cos(2x)=0.00000
sin(x)・sin(x)=3.14159
cos(2x)・cos(2x)=3.14160
cos(x)・sin(2x)=-0.00000
sin(x)・cos(2x)=0.00000考察
太郎くん
これもOKそう。
フクさん
コード自体もMATLABとほぼ一緒だな。
円周率が「pi」ではなく「%pi」になっている点に注意だ。
まとめ
フクさん
まとめだよ。
- 三角関数の直交性をScilabで確認してみた。
- 同一の関数及び角周波数の場合はπになり、それ以外は0になる。
バックナンバーはこちら。
マンガでわかるフーリエ解析
https://amzn.to/4azyPDj
amzn.to
手を動かしてまなぶ フーリエ解析・ラプラス変換
https://amzn.to/4bDH3vu
amzn.to
物理数学 量子力学のためのフーリエ解析・特殊関数
https://amzn.to/4bsygfU
amzn.to
単位が取れるフーリエ解析ノート
https://amzn.to/3V83fIl
amzn.to
今日から使えるフーリエ変換 普及版 式の意味を理解し、使いこなす
https://amzn.to/3ysbfvf
amzn.to
コメント