MathcotさんのHomePageです


フーリエ級数展開

Author: Mathcot

初版:2007.7.1

Update:2013.05.12


フー リエ級数(Fourier series)の定義式
一般の場合
FS-001a.PNG
 f(t)が周期Tの周期関数とする時、次の周期関数にフーリエ級数展開できる。
 wo=2π/t とおく。

 f(t)= a0/2 + Σ[n=1,∞] ( an cos(nwot) + bn sin(nwot) )

   a0= 2/T ∫[-T/2,T/2] f(t) dt

   an= 2/T ∫
[-T/2,T/2] f(t) cos(nwot) dt    ( n ≧1 )

   bn= 2/T ∫[-T/2,T/2] f(t) sin(nwot) dt    ( n ≧1 )

 または

   a0= 2/T ∫[0,T] f(t) dt

   an= 2/T ∫
[0,T] f(t) cos(nwot) dt    ( n ≧1 )

   bn= 2/T ∫[0,T] f(t) sin(nwot) dt    ( n ≧1 )

 
 複素指数によるフーリエ級数展開は次式で与えられる。

 f(t)= Σ[n=-∞,∞] ( cn exp(jnwot) )

   c0  = a0/2 = 1/T ∫[-T/2,T/2] f(t) dt

   cn  = ( an - i bn )/2 = 1/T ∫
[-T/2,T/2] f(t) exp( i nwot) dt    ( n ≧1 )

   c-n = ( an + i bn )/2     ( n ≧1 )

 または

   c0  = a0/2 = 1/T ∫[0,T] f(t) dt

   cn  = ( an - i bn )/2 = 1/T ∫
[0,T] f(t) exp( i nwot) dt    ( n ≧1 )

   c-n = ( an + i bn )/2     ( n ≧1 )

偶関数の場合
FS-001b.PNG

奇関数の場合
FS-001c.PNG


参考URL
[1]フー リエ級数(Wikipedia)
[2]
[3]
[4]
[5]

時間関数(周期T):1周期分
フーリエ級数 (周期 T=2π)
矩形波 (奇関数)
f(x) = 0 (-π≦ x < 0 )
    = 1 ( 0 ≦ x ≦π )
T = 2π(-π〜π)
FS-002-1a.PNG

積分区間:[-π→π]
a0 = 1/π∫f(x)dx  = 1
an = 1/π ∫f(x)cos(nx)dx = 0 (n≧1)
bn = 1/π ∫f(x)sin(nx)dx = 2/(nπ) (n:odd),
   = 0 (n:even)

f(x) = 1/2+(2/π){sin(x)+(1/3)sin(3x)+(1/5)sin(5x)+…}
矩形波(奇関数)
f(x) = -1 (-π≦ x < 0 )
    = 1 ( 0 ≦ x ≦π )
T = 2π(-π〜π)
FS-001-K02.PNG

T=2π,ωo=2π/T=1
an=0(n=0,1,2,…)
bn=(2/T)∫[-T/2,T/2] f(x)sin(2nπt/T)dt
 =(2/π)∫[0,π] sin(nt) dt
 =(2/π)[-cos(nx)/n][0,π] =(2/π){1-cos(nπ)}
 =(4/(nπ))(n=2m-1), =0(n=2m)
f(t)=(4/π)Σ[m=1,n] sin((2m-1)x)/(2m-1) (n→∞)


矩形波(偶関数)
f(x) = 0 (π/2≦ | x | ≦ π)
    = 1 ( | x | ≦π/2 )
T = 2π(-π〜π)
FS-002-1b.PNG
積分区間:[-π→π]
a0 = 1/π∫f(x)dx  = 1
an = 1/π ∫f(x)cos(nx)dx = -{(-1)^n}*2/(nπ) (n:odd),
   = 0 (n:even)
bn = 1/π ∫f(x)sin(nx)dx = 0

f(x) = 1/2+(2/π) { sin(x)-(1/3)sin(3x)+(1/5)sin(5x)-…}
鋸波
f(x) = x (-π≦x ≦π)
T = 2π(-π〜π)
FS-002-2a.PNG
積分区間:[-π→π]
an = (1/π) ∫x cos(nx)dx = 0
bn = (1/π) ∫x sin(nx)dx
   = (2/π) ∫x sin(nx)dx   (積分区間:[ 0 →π] )
   = (-1)^(n+1) 2/n
 
f(x) = 2{ sin(x)-(1/2)sin(2x)+(1/3)sin(3x)-・・・)


FS-002-N02.PNG
f(t)=t (t=-π〜π)
フーリエ級数展開係数はf(t)が奇関数であるから
 an=0 (n=0,1,2,…)
 bn=(4/T)∫(0〜T/2) t sin(nωot) (n=1,2,3,…)
  =(2/π)∫(0〜π) t sin(nt)dt
部分積分して
 bn=(2/π){[t(-1/n)cos(nt)](t:0〜π)+(1/n)∫(0〜π) cos(nt)dt}
  =(2/π){-(π/n)cos(nπ)+(1/n)[(1/n)sin(nt)(t:0〜π)}
  =(2/π){-(π/n)*(-1)^n+((1/n)^2)sin(nπ)}
  =(2/n)*(-1)^(n+1) (n=1,2,3,…)
フーリエ級数展開は
 f(t)=Σ(n=1,∞) (2/n)*((-1)^(n+1))sin(nt)
となる。
放物線y=(π^2)- (x^2) (-π≦x ≦π)
T = 2π(-π〜π)
FS-002-3c.PNG


FS-002-3b.PNG
 y=(2/3)π^2 +4cos(x)-cos(2x)+(4/9)cos(3x)-(1/4)cos(4x)+(4/25)cos(5x)-…
FS-002-3a.PNG

放物線y = x^2 (-π≦x ≦π)
T = 2π(-π〜π)
FS-002-4a.PNG


FS-002-4b.PNG

y=(1/3)(π^2)-4cos(x)+cos(2x)-(4/9)cos(3x)+(1/4)cos(4x)-(4/25)cos(5x)+…
y=cos (x/2) (-π≦x≦π)
 T=2π(-π〜π)
FS-002-5a.PNG


[Maple10]
> 
> an=
  
> 
  = (-1)^(n+1) *(4/π)(1/((4n^2)-1))
> 
  
  

フーリエ級数 展開 積分区間[0〜1]
周期 T=1
y=f(x)=x2 (0≦x≦1), f(x)=f(x+1)
To=1,ωo=2π
f(x)=a0/2 +Σ[n=1,∞] {an cos(nωo x)+bn sin(nωo x)}
a0=1/To)∫[0,To] f(x)dx
an=(1/To)∫[0,To] f(x)cos(nωo x)dx
bn=(1/To)∫[0,To] f(x)sin(nωo x)dx
FS-003-1a.PNG
[Maple10]
a0=∫[0,1] x2dx = x3/3| [0,1]=1/3
an=∫[0,1] x2cos (2nπx) dx =1/(2π2 n2 )
bn=∫[0,1] x2sin (2nπx) dx=-1/(2nπ)
f(x)=(1/6)
  +Σ[n=1,∞] [{1/(2π2n2)} cos(2nπx)-{1/(2nπ)}sin(2nπx)]




フーリエ級数 展開 積分区間 [-1〜1] 周期 T=2
y=f(x)=0 (-1≦x<0), x (0≦x≦1

[Maple10]

参考URL
[1]
[2]
[3]
[4]
[5]
 

Copyright(C) 2007-2013 Mathcot HI. All rights reserved.
初版:2007.7.1
Update:2007.12.16
Update:2013.05.11




inserted
      by FC2 system inserted by FC2 system