2

矩形波・のこぎり波のフーリエ展開の初等的計算

1859
0

前書き

矩形波、のこぎり波とは以下の図のような角ばった波だったり、
まさにのこぎりの刃のような波のことを言う。
矩形波 矩形波 のこぎり波 のこぎり波
名前自体はせいぜい概形を表してそう言うという水準で
具体的な数値を伴うものではないと思われるのでここではそれぞれ
次のような関数をそれとして取り扱う。

矩形波・のこぎり波

矩形波を表す関数
f(x)=1 (0<x<1/2),f(x)=1 (1/2<x<1),f(x+1)=f(x)
のこぎり波を表す関数
g(x)=x (0<x<1),g(x+1)=g(x)

どちらも不連続となる点での値を定義していないが、そこはいろいろと流儀があったりで
今回そこにこだわった話はないので、一旦未定義としておく。
これらは区分的に連続な周期関数なのでフーリエ展開できる。

矩形波・のこぎり波のフーリエ展開

f(x)=4πk=0sin((2k+1)2πx)2k+1

g(x)=122πk=1sin(2kπx)2k

この記事では、これらのフーリエ展開について、極力初等的な範疇で証明することを目的とする。
なお、のこぎり波についての同趣旨の記事として
tyamada氏の Gauss関数のFourier級数表示の初等的証明 があるので
本記事のタイトルで興味をもって開いてくれた方は
こちらもご覧いただくことをお勧めする。
本記事は上記記事とは異なる方針となる。

方針

f(x),g(x)をなす級数の部分和である三角多項式は
そのままでは計算が難しいがそれを項別微分した物の
計算は難しくない。なのでここでは
(1)部分和を微分した物を計算し、
(2)それを積分して、
(3)極限を取る。
という直球勝負の作戦をとる。
(この直球勝負がしたかったというのが本記事の趣旨である)
ということで、和の計算を先に確認する

(1)cosの和の計算

2k=0n1cos((2k+1)x)=sin(2nx)sin(x)

2k=1n1cos(2kx)=sin((2n1)x)sin(x)1
※和の取り方が上下で若干違う(k=0k=1)ことに注意

いわゆる積和公式により
2sin(x)cos((2k+1)x)=sin(x+(2k+1)x)+sin(x(2k+1)x)=sin(2(k+1)x)sin(2kx)
ゆえ
2sin(x)k=0n1sin(x)cos((2k+1)x)=1sin(x)k=0n1(sin(2(k+1)x)sin(2kx))=sin(2nx)sin(x)
同様に
2sin(x)cos(2kx)=sin(x+2kx)+sin(x2kx)=sin((2k+1)x)sin((2k1)x)
ゆえ
2sin(x)k=1n1sin(x)cos(2kx)=1sin(x)k=1n1(sin((2k+1)x)sin((2k1)x))=sin((2n1)x)sin(x)sin(x)=sin((2n1)x)sin(x)1

(2) 部分和の計算

得られた和を積分してやれば部分和の表示を得ることができる。
次のようになる。

0<x<πのとき
2k=0n1sin((2k+1)x)2k+12k=0n1(1)k2k+1=π/2xsin(2nt)sin(t)dt
2k=1n1sin(2kx)2k=π/2xsin((2n1)t)sin(t)dtx+π2

証明はどちらも補題2の式をπ/2からxまでで定積分してやればよい。
(sin(x)=0とならないようxの範囲を制限している)

(3) 極限の計算

ここまでの結果からnの極限を考えるにあたって検討が必要となるのは次の形の極限である。

mを自然数としてmのとき
π/2xsin(mt)sin(t)dt0 (0<x<π)

これは「リーマン・ルベーグの補題」として知られる結果のごく特別な場合となるが、ここではそういったことを踏まえず計算によって示す。

分子を積分する形で部分積分を考えると
π/2xsin(mt)sin(t)dt=[cos(mt)msin(t)]π/2x1mπ/2xcos(mt)cos(t)sin2(t)dt
        =1m[cos(mt)sin(t)]π/2x1mπ/2xcos(mt)cos(t)sin2(t)dt
右辺の第1項および第2項はともに1/mを係数に持ち
0<x<πであることから0<sin(x)であるので
発散しない有限の値を取る。
よってmでともに0に収束する。すなわち
π/2xsin(mt)sin(t)dt0
が言える。

この補題により計算を進めることができる

 定理1の証明

補題4より、まず
2k=1n1sin(2kx)2k=π/2xsin((2n1)t)sin(t)dtx+π2
についてnとすると、積分の項が0となるので
2k=1sin(2kx)2k=x+π2
x=π22k=1sin(2kx)2k
xπxと置き換えると
πx=π22k=1sin(2kπx)2k
x=122πk=1sin(2kπx)2k (0<x<1)
右辺の無限級数の項であるsin(2kπx)はいずれの自然数kについても周期1を持つ周期関数であるので
以上より定理1ののこぎり波g(x)の場合が示せた。

続いて
2k=0n1sin((2k+1)x)2k+12k=0n1(1)k2k+1=π/2xsin(2nt)sin(t)dt
について、nとすると、積分の項が0となり
左辺第2項はいわゆるライプニッツ級数となるので
2k=0sin((2k+1)x)2k+12π4=0
k=0sin((2k+1)x)2k+1=π4 (0<x<π)
x2πxと置き換えると
k=0sin((2k+1)2πx)2k+1=π4 (0<x<1/2)

さらに上記についてxx1/2に置き換えると
k=0sin((2k+1)2π(x1/2))2k+1=π4 (1/2<x<1)
 sin((2k+1)2π(x1/2))
=sin((2k+1)2πx)cos((2k+1)2π1/2)cos((2k+1)2πx)sin((2k+1)2π1/2)
=sin((2k+1)2πx)cos((2k+1)π)cos((2k+1)2πx)sin((2k+1)π)
=sin((2k+1)2πx)(1)cos((2k+1)2πx)0
=sin((2k+1)2πx)
ゆえ
k=0sin((2k+1)2πx)2k+1=π4 (1/2<x<1)
k=0sin((2k+1)2πx)2k+1=π4 (1/2<x<1)
ここまでの結果とsin((2k+1)2πx)が周期1をもつ周期関数であることから
係数π/4を整理することで定理1の矩形波f(x)の場合が示せる。

終わりに

フーリエ級数展開を見たままの形からどういった関数になるのかという計算をできたので個人的な目的な十分に達せられ、ほぼほぼ満足した。
矩形波でライプニッツ級数が出てくるところは、うまいこと処理すれば
π/4を値にとる平易な積分になったりしないかと思わなくもないが
とはいえライプニッツ級数の計算自体もこの記事で現れるのと同程度の道具立てで計算できるので、まぁよかろう。

本記事で扱った関数はどちらもsinの級数であるが、cosの級数であっても
係数が同じであれば計算は同じようにできる。
関数としてはlog(sin(x))のようなやや見慣れない関数になるのだが
これを利用していろいろやることを考えており、これもいずれ記事にしたい。

投稿日:2023113
更新日:2023113
OptHub AI Competition

この記事を高評価した人

高評価したユーザはいません

この記事に送られたバッジ

バッジはありません。
バッチを贈って投稿者を応援しよう

バッチを贈ると投稿者に現金やAmazonのギフトカードが還元されます。

投稿者

TKSS
TKSS
45
5049

コメント

他の人のコメント

コメントはありません。
読み込み中...
読み込み中
  1. 前書き
  2. 方針
  3. (1)cosの和の計算
  4. (2) 部分和の計算
  5. (3) 極限の計算
  6. 終わりに