目录



Fourier Series (傅里叶级数)


三角函数的正交性


  • 可以定义一个内积空间,内积空间中的每个向量均为周期为




    2


    π



    2\pi


    2π
    的实值函数
    ,函数之间的内积操作定义如下:








    f


    ,


    g





    =






    0



    2


    π




    f


    (


    x


    )


    g


    (


    x


    )


    d


    x



    \langle f, g\rangle=\int_{0}^{2 \pi} f(x) g(x) d x


    f,g=02πf(x)g(x)dx
    其中积分上下限只需满足积分范围为




    2


    π



    2\pi


    2π
    即可
    ,例如








    0



    2


    π





    \int_0^{2\pi}


    02π













    π



    π




    \int_{-\pi}^{\pi}


    ππ
    都可以。该内积操作可以看作是





    R


    n




    \R^n


    Rn
    上向量内积的推广,即把函数看作具有无数个元素的向量。如果函数定义域为




    [


    0


    ,


    T


    ]



    [0,T]


    [0,T]
    而非




    [


    0


    ,


    2


    π


    ]



    [0,2\pi]


    [0,2π]
    ,则内积操作的积分上下限也要相应地变换
  • 三角函数系




    {


    1


    ,


    sin





    x


    ,


    cos





    x


    ,


    sin





    2


    x


    ,


    cos





    2


    x


    ,


    .


    .


    .


    ,


    sin





    n


    x


    ,


    cos





    n


    x


    ,


    .


    .


    .


    }



    {1,\sin x,\cos x,\sin 2x,\cos 2x,…,\sin nx,\cos nx,…}


    {1,sinx,cosx,sin2x,cos2x,,sinnx,cosnx,}
    (




    n


    =


    0


    ,


    1


    ,


    2


    ,


    .


    .


    .



    n=0,1,2,…


    n=0,1,2,
    ) 是上述内积空间的基
    ,即它们之间两两正交 (注意到,(1)




    2


    π



    2\pi


    2π





    cos





    n


    x


    ,


    sin





    n


    x



    \cos nx,\sin nx


    cosnx,sinnx
    的周期
    ,并且 (2) 内积操作的积分区间长度为




    2


    π



    2\pi


    2π
    (不一定是




    [


    0


    ,


    2


    π


    ]



    [0,2\pi]


    [0,2π]
    ,只要长度为




    2


    π



    2\pi


    2π
    即可)),例如












    cos





    m


    x


    ,


    cos





    n


    x












    =






    0



    2


    π




    cos





    m


    x


    cos





    n


    x


    d


    x
















    =



    1


    2







    0



    2


    π




    [


    cos





    (


    m


    x


    +


    n


    x


    )


    +


    cos





    (


    m


    x





    n


    x


    )


    ]


    d


    x
















    =





    1


    2




    [




    sin





    (


    m


    x


    +


    n


    x


    )




    m


    +


    n




    +




    sin





    (


    m


    x





    n


    x


    )




    m





    n




    ]







    0



    2


    π




    =


    0













    cos





    m


    x


    ,


    cos





    m


    x












    =






    0



    2


    π






    cos






    2



    (


    m


    x


    )


    d


    t
















    =



    1


    2







    0



    2


    π




    [


    1


    +


    cos





    (


    2


    m


    x


    )


    ]


    d


    t
















    =





    1


    2




    [


    x


    +




    sin





    (


    2


    m


    x


    )




    2


    m




    ]







    0



    2


    π




    =


    π







    cosmx,cosnx=02πcosmxcosnxdx=1202π[cos(mx+nx)+cos(mxnx)]dx=12[sin(mx+nx)m+n+sin(mxnx)mn]|02π=0cosmx,cosmx=02πcos2(mx)dt=1202π[1+cos(2mx)]dt=12[x+sin(2mx)2m]|02π=π



    cosmx,cosnxcosmx,cosmx=02πcosmxcosnxdx=2102π[cos(mx+nx)+cos(mxnx)]dx=21[m+nsin(mx+nx)+mnsin(mxnx)]


    02π
    =0
    =02πcos2(mx)dt=2102π[1+cos(2mx)]dt=21[x+2msin(2mx)]


    02π
    =π

周期为




2


π



2\pi


2π
的函数展开为傅里叶级数







f


(


t


)


=


f


(


t


+


2


π


)



f(t)=f(t+2\pi)


f(t)=f(t+2π)


  • 由于三角函数系是




    C


    [


    0


    ,


    2


    π


    ]



    C[0,2\pi]


    C[0,2π]
    的内积空间上的一组基,因此任何定义域在




    [


    0


    ,


    2


    π


    ]



    [0,2\pi]


    [0,2π]
    上的连续函数均可以表示为三角函数的线性组合
    ,下式即为傅里叶级数





    f


    (


    x


    )


    =




    a


    0



    2



    +







    n


    =


    1







    (



    a


    n



    cos





    n


    x


    +



    b


    n



    sin





    n


    x


    )



    f(x)=\frac{a_0}{2}+\sum_{n=1}^\infty(a_n\cos nx+b_n\sin nx)


    f(x)=2a0+n=1(ancosnx+bnsinnx)
    其中,





    a


    k



    ,



    b


    k




    a_k,b_k


    ak,bk





    f



    f


    f
    Fourier coefficients. 傅里叶系数相当于是




    f



    f


    f
    在各个基向量上的投影,因此由正交投影可知,






    a


    k



    =







    f


    ,


    cos





    k


    x










    cos





    k


    x


    ,


    cos





    k


    x







    ,




    b


    k



    =







    f


    ,


    sin





    k


    x










    sin





    k


    x


    ,


    sin





    k


    x







    ,



    k





    1



    a_{k}=\frac{\langle f, \cos k x\rangle}{\langle\cos k x, \cos k x\rangle}, \quad b_{k}=\frac{\langle f, \sin k x\rangle}{\langle\sin k x, \sin k x\rangle}, \quad k \geq 1


    ak=coskx,coskxf,coskx,bk=sinkx,sinkxf,sinkx,k1
    由于







    cos





    k


    x


    ,


    cos





    k


    x





    =


    π



    \langle \cos kx, \cos kx\rangle =\pi


    coskx,coskx=π








    sin





    k


    x


    ,


    sin





    k


    x





    =


    π



    \langle \sin kx, \sin kx\rangle=\pi


    sinkx,sinkx=π
    ,因此有






    a


    k



    =



    1


    π







    0



    2


    π




    f


    (


    x


    )


    cos





    k


    x


    d


    x


    ,




    b


    k



    =



    1


    π







    0



    2


    π




    f


    (


    x


    )


    sin





    k


    x


    d


    x


    ,



    k





    1



    a_{k}=\frac{1}{\pi} \int_{0}^{2 \pi} f(x) \cos k xd x, \quad b_{k}=\frac{1}{\pi} \int_{0}^{2 \pi} f(x) \sin k xd x, \quad k \geq 1


    ak=π102πf(x)coskxdx,bk=π102πf(x)sinkxdx,k1
    而 (constant) function




    1



    1


    1
    的傅里叶系数为










    f


    ,


    1










    1


    ,


    1







    =



    1



    2


    π








    0



    2


    π




    f


    (


    x


    )





    1


    d


    x


    =



    1


    2




    [



    1


    π







    0



    2


    π




    f


    (


    x


    )


    cos





    (


    0





    x


    )


    d


    x


    ]



    =




    a


    0



    2




    \frac{\langle f, 1\rangle}{\langle 1,1\rangle}=\frac{1}{2 \pi} \int_{0}^{2 \pi} f(x) \cdot 1 d x=\frac{1}{2}\left[\frac{1}{\pi} \int_{0}^{2 \pi} f(x) \cos (0 \cdot x) d x\right]=\frac{a_{0}}{2}


    1,1f,1=2π102πf(x)1dx=21[π102πf(x)cos(0x)dx]=2a0


除了用投影的思想以外,也可以用如下方式求出傅里叶系数。例如,要求





a


k




a_k


ak
时,有












f


(


x


)


,


cos





k


x












=






0



2


π






a


0



2



cos





0


x


cos





k


x


d


x


+







n


=


1







(






0



2


π





a


n



cos





n


x


cos





k


x


d


x


+






0



2


π





b


n



sin





n


x


cos





k


x


d


x


)







f(x),coskx=02πa02cos0xcoskxdx+n=1(02πancosnxcoskxdx+02πbnsinnxcoskxdx)



f(x),coskx=02π2a0cos0xcoskxdx+n=1(02πancosnxcoskxdx+02πbnsinnxcoskxdx)由三角函数的正交性可知,












f


(


x


)


,


cos





k


x












=






0



2


π





a


k



cos





k


x


cos





k


x


d


x


=


π



a


k








f(x),coskx=02πakcoskxcoskxdx=πak



f(x),coskx=02πakcoskxcoskxdx=πak







  • m



    m


    m
    阶傅里叶近似
    即为





    f


    (


    x


    )







    a


    0



    2



    +







    n


    =


    1



    m



    (



    a


    n



    cos





    n


    x


    +



    b


    n



    sin





    n


    x


    )



    f(x)\approx\frac{a_0}{2}+\sum_{n=1}^m(a_n\cos nx+b_n\sin nx)


    f(x)2a0+n=1m(ancosnx+bnsinnx)
    可以证明,当




    m









    m\rightarrow\infty


    m
    时,




    f



    f


    f
    和傅里叶近似之间的均方误差也趋近于 0.




    m



    m


    m
    阶傅里叶近似相当于是




    f



    f


    f





    {


    1


    ,


    sin





    x


    ,


    cos





    x


    ,


    sin





    2


    x


    ,


    cos





    2


    x


    ,


    .


    .


    .


    ,


    sin





    m


    x


    ,


    cos





    m


    x


    }



    {1,\sin x,\cos x,\sin 2x,\cos 2x,…,\sin mx,\cos mx}


    {1,sinx,cosx,sin2x,cos2x,,sinmx,cosmx}
    张成的子空间




    W



    W


    W
    上的正交投影,并且由正交投影的线性性质可知







    proj






    W



    (


    f


    +


    g


    )


    =




    proj






    W



    f


    +




    proj






    W



    g



    \operatorname{proj}_{W}(f+g)=\operatorname{proj}_{W} f+\operatorname{proj}_{W} g


    projW(f+g)=projWf+projWg





    f


    +


    g



    f+g


    f+g





    m



    m


    m
    阶傅里叶近似等于




    f



    f


    f





    m



    m


    m
    阶傅里叶近似和




    g



    g


    g





    m



    m


    m
    阶傅里叶近似之和



EXAMPLE


Find the




n



n


n
th-order Fourier approximation to the function




f


(


t


)


=


t



f(t)= t


f(t)=t
on the interval




[


0


,


2


π


]



[0, 2\pi]


[0,2π]
.


SOLUTION


  • Compute
    and for




    k


    >


    0



    k > 0


    k>0
    , using integration by parts,
  • Thus the




    n



    n


    n
    th-order Fourier approximation of




    f


    (


    t


    )


    =


    t



    f(t)= t


    f(t)=t
    is
    可以看到,三角函数频率越高系数越小,




    f


    (


    t


    )



    f(t)


    f(t)
    低频分量占了主要部分
  • Figure 3 shows the third- and fourth-order Fourier approximations of




    f



    f


    f
    .

周期为




2


L



2L


2L
的函数展开为傅里叶级数







f


(


t


)


=


f


(


t


+


2


L


)



f(t)=f(t+2L)


f(t)=f(t+2L)


  • 函数




    f


    (


    t


    )



    f(t)


    f(t)
    周期为




    2


    L



    2L


    2L
    时需要进行换元,令




    x


    =



    π


    L



    t



    x=\frac{\pi}{L}t


    x=Lπt
    ,则




    t


    =



    L


    π



    x



    t=\frac{L}{\pi}x


    t=πLx





    g


    (


    x


    )


    =


    f


    (


    t


    )


    =


    f


    (



    L


    π



    x


    )



    g(x)=f(t)=f(\frac{L}{\pi}x)


    g(x)=f(t)=f(πLx)
    的周期为




    2


    π



    2\pi


    2π
    .




    g


    (


    x


    )



    g(x)


    g(x)
    的傅里叶级数






    g


    (


    x


    )


    =




    a


    0



    2



    +







    n


    =


    1







    (



    a


    n



    cos





    n


    x


    +



    b


    n



    sin





    n


    x


    )



    g(x)=\frac{a_0}{2}+\sum_{n=1}^\infty(a_n\cos nx+b_n\sin nx)


    g(x)=2a0+n=1(ancosnx+bnsinnx)
    其中,






    a


    k



    =



    1


    π







    0



    2


    π




    g


    (


    x


    )


    cos





    k


    x


    d


    x


    ,




    b


    k



    =



    1


    π







    0



    2


    π




    g


    (


    x


    )


    sin





    k


    x


    d


    x


    ,



    k





    0




    a


    0



    =



    1


    π







    0



    2


    π




    g


    (


    x


    )


    d


    x



    a_{k}=\frac{1}{\pi} \int_{0}^{2 \pi} g(x) \cos k xd x, \quad b_{k}=\frac{1}{\pi} \int_{0}^{2 \pi} g(x) \sin k xd x, \quad k \geq 0\ a_{0}=\frac{1}{\pi} \int_{0}^{2 \pi} g(x)d x


    ak=π102πg(x)coskxdx,bk=π102πg(x)sinkxdx,k0a0=π102πg(x)dx





    x


    =



    π


    L



    t



    x=\frac{\pi}{L}t


    x=Lπt






    1


    π







    0



    2


    π




    d


    x


    =



    1


    L







    0



    2


    L




    d


    t



    \frac{1}{\pi}\int_0^{2\pi} dx=\frac{1}{L}\int_0^{2L} dt


    π102πdx=L102Ldt
    代入可知,




    f


    (


    t


    )



    f(t)


    f(t)
    的傅里叶级数






    f


    (


    t


    )


    =




    a


    0



    2



    +







    n


    =


    1







    (



    a


    n



    cos







    n


    π



    L



    t


    +



    b


    n



    sin







    n


    π



    L



    t


    )



    f(t)=\frac{a_0}{2}+\sum_{n=1}^\infty(a_n\cos \frac{n\pi}{L}t+b_n\sin \frac{n\pi}{L}t)


    f(t)=2a0+n=1(ancosLt+bnsinLt)
    其中,






    a


    k



    =



    1


    L







    0



    2


    L




    f


    (


    t


    )


    cos







    k


    π



    L



    t


    d


    t


    ,




    b


    k



    =



    1


    L







    0



    2


    L




    f


    (


    t


    )


    sin







    k


    π



    L



    t


    d


    t


    ,



    k





    0




    a


    0



    =



    1


    L







    0



    2


    L




    f


    (


    t


    )


    d


    t



    a_{k}=\frac{1}{L} \int_{0}^{2 L} f(t) \cos \frac{k\pi}{L}td t, \quad b_{k}=\frac{1}{L} \int_{0}^{2 L} f(t) \sin \frac{k\pi}{L}td t, \quad k \geq 0\ a_{0}=\frac{1}{L} \int_{0}^{2L} f(t)d t


    ak=L102Lf(t)cosLtdt,bk=L102Lf(t)sinLtdt,k0a0=L102Lf(t)dt
  • 基频率




    w


    =




    2


    π



    T



    =



    π


    L




    w=\frac{2\pi}{T}=\frac{\pi}{L}


    w=T2π=Lπ
    ,则周期为




    T


    =


    2


    L



    T=2L


    T=2L
    的函数展开为傅里叶级数
    可以写为如下形式:





    f


    (


    t


    )


    =




    a


    0



    2



    +







    n


    =


    1







    (



    a


    n



    cos





    (


    n


    w


    t


    )


    +



    b


    n



    sin





    (


    n


    w


    t


    )


    )



    f(t)=\frac{a_0}{2}+\sum_{n=1}^\infty(a_n\cos (nwt)+b_n\sin (nwt))


    f(t)=2a0+n=1(ancos(nwt)+bnsin(nwt))
    其中,






    a


    k



    =



    2


    T







    0


    T



    f


    (


    t


    )


    cos





    (


    k


    w


    t


    )


    d


    t


    ,




    b


    k



    =



    2


    T







    0


    T



    f


    (


    t


    )


    sin





    (


    k


    w


    t


    )


    d


    t


    ,



    k





    0




    a


    0



    =



    2


    T







    0


    T



    f


    (


    t


    )


    d


    t



    a_{k}=\frac{2}{T} \int_{0}^{T} f(t) \cos (kwt)d t, \quad b_{k}=\frac{2}{T} \int_{0}^{T} f(t) \sin (kwt)d t, \quad k \geq 0\ a_{0}=\frac{2}{T} \int_{0}^{T} f(t)d t


    ak=T20Tf(t)cos(kwt)dt,bk=T20Tf(t)sin(kwt)dt,k0a0=T20Tf(t)dt
    因此,任何周期函数都可以看作是不同振幅,不同相位正弦波的叠加。在频域,0 频率也被称为直流分量,在傅里叶级数的叠加中,它仅仅影响全部波形相对于数轴整体向上或是向下而不改变波的形状

傅里叶级数的复数形式


  • 欧拉公式





    e



    i


    θ




    =


    cos





    θ


    +


    i


    sin





    θ



    e^{i\theta}=\cos\theta+i\sin \theta


    eiθ=cosθ+isinθ
    可得,





    cos





    θ


    =



    1


    2



    (



    e



    i


    θ




    +



    e






    i


    θ




    )



    sin





    θ


    =






    1


    2



    i


    (



    e



    i


    θ








    e






    i


    θ




    )



    \cos\theta=\frac{1}{2}(e^{i\theta}+e^{-i\theta})\ \sin\theta=-\frac{1}{2}i(e^{i\theta}-e^{-i\theta})


    cosθ=21(eiθ+eiθ)sinθ=21i(eiθeiθ)
    将上述 2 式代入傅里叶级数展开式有









    f


    (


    t


    )









    =




    a


    0



    2



    +







    n


    =


    1







    (



    a


    n



    cos





    (


    n


    w


    t


    )


    +



    b


    n



    sin





    (


    n


    w


    t


    )


    )
















    =




    a


    0



    2



    +







    n


    =


    1







    (



    a


    n




    1


    2



    (



    e



    i


    n


    w


    t




    +



    e






    i


    n


    w


    t




    )






    b


    n




    1


    2



    i


    (



    e



    i


    n


    w


    t








    e






    i


    n


    w


    t




    )


    )
















    =




    a


    0



    2



    +







    n


    =


    1










    a


    n






    i



    b


    n




    2




    e



    i


    n


    w


    t




    +







    n


    =


    1










    a


    n



    +


    i



    b


    n




    2




    e






    i


    n


    w


    t


















    =







    n


    =


    0



    0





    a


    0



    2




    e



    i


    n


    w


    t




    +







    n


    =


    1










    a


    n






    i



    b


    n




    2




    e



    i


    n


    w


    t




    +







    n


    =













    1







    a






    n




    +


    i



    b






    n





    2




    e



    i


    n


    w


    t


















    =







    n


    =














    C


    n




    e



    i


    n


    w


    t









    f(t)=a02+n=1(ancos(nwt)+bnsin(nwt))=a02+n=1(an12(einwt+einwt)bn12i(einwteinwt))=a02+n=1anibn2einwt+n=1an+ibn2einwt=n=00a02einwt+n=1anibn2einwt+n=1an+ibn2einwt=n=Cneinwt



    f(t)=2a0+n=1(ancos(nwt)+bnsin(nwt))=2a0+n=1(an21(einwt+einwt)bn21i(einwteinwt))=2a0+n=12anibneinwt+n=12an+ibneinwt=n=002a0einwt+n=12anibneinwt+n=12an+ibneinwt=n=Cneinwt
    其中,






    C


    n



    =



    {















    a


    0



    2








    n


    =


    0



















    a


    n






    i



    b


    n




    2





      


    n


    >


    0



















    a






    n




    +


    i



    b






    n





    2




      


    n


    <


    0









    C_n=\left{
    a02n=0anibn2  n>0an+ibn2  n<0
    \right.


    Cn=





    2a0n=02anibn  n>02an+ibn  n<0
  • 傅里叶级数系数的表达式如下:






    a


    k



    =



    2


    T







    0


    T



    f


    (


    t


    )


    cos





    (


    k


    w


    t


    )


    d


    t


    ,




    b


    k



    =



    2


    T







    0


    T



    f


    (


    t


    )


    sin





    (


    k


    w


    t


    )


    d


    t


    ,



    k





    0




    a


    0



    =



    2


    T







    0


    T



    f


    (


    t


    )


    d


    t



    a_{k}=\frac{2}{T} \int_{0}^{T} f(t) \cos (kwt)d t, \quad b_{k}=\frac{2}{T} \int_{0}^{T} f(t) \sin (kwt)d t, \quad k \geq 0\ a_{0}=\frac{2}{T} \int_{0}^{T} f(t)d t


    ak=T20Tf(t)cos(kwt)dt,bk=T20Tf(t)sin(kwt)dt,k0a0=T20Tf(t)dt
    将其带入





    C


    n




    C_n


    Cn
    可得






    C


    n



    =



    {














    1


    T







    0


    T



    f


    (


    t


    )


    d


    t


    =



    1


    T




    (






    0


    T



    f


    (


    t


    )


    [


    cos





    (


    0


    )





    i


    sin





    (


    0


    )


    ]


    d


    t


    )





    n


    =


    0

















    1


    T




    (






    0


    T



    f


    (


    t


    )


    [


    cos





    (


    n


    w


    t


    )





    i


    sin





    (


    n


    w


    t


    )


    ]


    d


    t


    )




      


    n


    >


    0

















    1


    T




    (






    0


    T



    f


    (


    t


    )


    [


    cos





    (


    n


    w


    t


    )





    i


    sin





    (


    n


    w


    t


    )


    ]


    d


    t


    )




      


    n


    <


    0









    C_n=\left{
    1T0Tf(t)dt=1T(0Tf(t)[cos(0)isin(0)]dt)n=01T(0Tf(t)[cos(nwt)isin(nwt)]dt)  n>01T(0Tf(t)[cos(nwt)isin(nwt)]dt)  n<0
    \right.


    Cn=





    T10Tf(t)dt=T1(0Tf(t)[cos(0)isin(0)]dt)n=0T1(0Tf(t)[cos(nwt)isin(nwt)]dt)  n>0T1(0Tf(t)[cos(nwt)isin(nwt)]dt)  n<0
    整理可得









    C


    n









    =



    1


    T




    (






    0


    T



    f


    (


    t


    )


    [


    cos





    (


    n


    w


    t


    )





    i


    sin





    (


    n


    w


    t


    )


    ]


    d


    t


    )

















    =



    1


    T




    (






    0


    T



    f


    (


    t


    )



    e






    i


    n


    w


    t




    d


    t


    )








    Cn=1T(0Tf(t)[cos(nwt)isin(nwt)]dt)=1T(0Tf(t)einwtdt)



    Cn=T1(0Tf(t)[cos(nwt)isin(nwt)]dt)=T1(0Tf(t)einwtdt)
  • 总结可得,傅里叶级数的复数形式









    f


    (


    t


    )









    =







    n


    =














    C


    n




    e



    i


    n


    w


    t












    s


    .


    t


    .




    C


    n










    =



    1


    T




    (






    0


    T



    f


    (


    t


    )



    e






    i


    n


    w


    t




    d


    t


    )











    f


    (


    t


    )









    =


    f


    (


    t


    +


    T


    )







    f(t)=n=Cneinwts.t.Cn=1T(0Tf(t)einwtdt)f(t)=f(t+T)



    f(t)s.t.Cnf(t)=n=Cneinwt=T1(0Tf(t)einwtdt)=f(t+T)
    其中,




    w


    =




    2


    π



    T




    w=\frac{2\pi}{T}


    w=T2π
    基频率

复数基函数的性质


  • 下面我们可以先看一下当




    n



    n


    n
    取不同值时




    f


    (


    t


    )


    =



    e



    i


    n


    w


    t





    f(t)=e^{inwt}


    f(t)=einwt
    是怎样的函数,















    n


    =


    0




    f


    (


    t


    )


    =



    e


    0



    =


    1
















    n


    =


    1




    f


    (


    t


    )


    =



    e



    i


    w


    t






    T


    =




    2


    π



    w

















    n


    =


    2




    f


    (


    t


    )


    =



    e



    2


    i


    w


    t






    T


    =



    π


    w

















    .


    .


    .
















    n


    =





    1




    f


    (


    t


    )


    =



    e






    i


    w


    t






    T


    =




    2


    π



    w

















    n


    =





    2




    f


    (


    t


    )


    =



    e






    2


    i


    w


    t






    T


    =



    π


    w

















    .


    .


    .







    n=0f(t)=e0=1n=1f(t)=eiwtT=2πwn=2f(t)=e2iwtT=πw...n=1f(t)=eiwtT=2πwn=2f(t)=e2iwtT=πw...



    n=0f(t)=e0=1n=1f(t)=eiwtT=w2πn=2f(t)=e2iwtT=wπn=1f(t)=eiwtT=w2πn=2f(t)=e2iwtT=wπ
    其中,




    n


    =


    0



    n=0


    n=0
    时,




    f


    (


    t


    )



    f(t)


    f(t)
    是一个常函数。




    n


    >


    0



    n>0


    n>0
    时,




    f


    (


    t


    )


    =



    e



    i


    n


    w


    t




    =


    cos





    n


    w


    t


    +


    i


    sin





    n


    w


    t



    f(t)=e^{inwt}=\cos nwt+i\sin nwt


    f(t)=einwt=cosnwt+isinnwt
    可以看作是复平面上沿着单位圆逆时针旋转的函数,周期为转一圈所需的时间也就是




    n


    w


    T


    =


    2


    π





    T


    =




    2


    π




    n


    w





    nwT=2\pi\Rightarrow T=\frac{2\pi}{nw}


    nwT=2πT=nw2π
    .




    n


    <


    0



    n<0


    n<0
    时,




    f


    (


    t


    )


    =



    e






    i


    n


    w


    t




    =


    cos





    n


    w


    t





    i


    sin





    n


    w


    t



    f(t)=e^{-inwt}=\cos nwt-i\sin nwt


    f(t)=einwt=cosnwtisinnwt
    可以看作是复平面上沿着单位圆顺时针旋转的函数,周期也是




    T


    =




    2


    π




    n


    w





    T=\frac{2\pi}{nw}


    T=nw2π
    . 除了




    n


    =


    0



    n=0


    n=0
    以外,




    T


    =




    2


    π



    w




    T=\frac{2\pi}{w}


    T=w2π
    是每个函数的公共周期,并且对应同一周期的函数有两个
  • 事实上,上述函数




    {


    .


    .


    .


    ,



    e






    3


    i


    w


    t




    ,



    e






    2


    i


    w


    t




    ,



    e






    i


    w


    t




    ,


    1


    ,



    e



    i


    w


    t




    ,



    e



    2


    i


    w


    t




    ,



    e



    3


    i


    w


    t




    ,


    .


    .


    .


    }



    {…,e^{-3iwt},e^{-2iwt},e^{-iwt},1,e^{iwt},e^{2iwt},e^{3iwt},…}


    {,e3iwt,e2iwt,eiwt,1,eiwt,e2iwt,e3iwt,}
    周期为




    T



    T


    T
    的复值函数的一组复数基
    ,它的内积被定义为 (内积需满足







    f


    ,


    f








    0



    \langle f, f\rangle\geq0


    f,f0
    )








    f


    ,


    g





    =






    0


    T



    f


    (


    x


    )




    g


    (


    x


    )







    d


    x



    \langle f, g\rangle=\int_{0}^{T} f(x) \overline{g(x)} d x


    f,g=0Tf(x)g(x)dx
    其中,






    g


    (


    x


    )








    \overline{g(x)}


    g(x)





    g


    (


    x


    )



    g(x)


    g(x)
    的共轭。因此有













    e



    i


    m


    w


    t




    ,



    e



    i


    n


    w


    t














    =






    0


    T




    e



    i


    m


    w


    t





    e






    i


    n


    w


    t




    d


    t
















    =






    0


    T




    e



    i


    (


    m





    n


    )


    w


    t




    d


    t
















    =



    {













    0



     


    m





    n
















    T



    m


    =


    n













    <span class=”MathJax” id=”MathJax-Element-10-Frame” tabindex=”0” style=”position: relative;” data-mathml=”\begin{aligned}&#xA0;\langle&#xA0;e^{imwt},&#xA0;e^{inwt}\rangle&amp;=\int_{0}^{T}&#xA0;e^{imwt}e^{-inwt}&#xA0;d&#xA0;t\&#xA0;&amp;=\int_{0}^{T}&#xA0;e^{i(m-n)wt}&#xA0;d&#xA0;t&#xA0;\&amp;=\left{\begin{aligned}&amp;0\quad&#xA0;\&#xA0;m\neq&#xA0;n\&amp;T\quad&#xA0;m=&#xA0;n\end{aligned}“ role=”presentation”>\begin{aligned} \langle e^{imwt}, e^{inwt}\rangle&=\int_{0}^{T} e^{imwt}e^{-inwt} d t\ &=\int_{0}^{T} e^{i(m-n)wt} d t \&=\left{\begin{aligned}&0\quad \ m\neq n\&T\quad m= n\end{aligned}
    \right. \end{aligned}


    eimwt,einwt=0Teimwteinwtdt=0Tei(mn)wtdt={0 m=nTm=n

Fourier Transform (FT, 傅里叶变换)


从时域到频域


  • 傅里叶级数的复数形式为





    f


    (


    t


    )


    =







    n


    =














    C


    n




    e



    i


    n


    w


    t





    f(t)=\sum_{n=-\infty}^\infty C_ne^{inwt}


    f(t)=n=Cneinwt
    可以看到,区分不同周期函数的仅仅是





    C


    n




    C_n


    Cn
    ,借此可以将一个周期函数从时域转换为频域
    上图将





    C


    n




    C_n


    Cn
    表示在了复平面中,实部和虚部都表示出来了 (振幅谱 + 相位谱)。但有时会将





    C


    n




    C_n


    Cn
    的幅度








    C


    n







    |C_n|


    Cn
    拿出来用以表示




    f


    (


    t


    )



    f(t)


    f(t)
    在不同频率下的强度 (振幅谱)
  • 很多在时域看似不可能做到的数学操作,在频域相反很容易。这就是需要傅里叶变换的地方。尤其是从某条曲线中去除一些特定的频率成分,这在工程上称为滤波,是信号处理最重要的概念之一,只有在频域才能轻松的做到;在图像处理领域,离散傅里叶变换将图像的离散像素数据转化为频谱数据,进行高频分量的舍弃与处理可以实现图像滤波以及图像压缩等目的

傅里叶变换及其逆变换


  • 上面我们展示了周期函数时域到频域的变换,而非周期函数呢?将非周期函数由时域转换为频域的变换即为傅里叶变换。非周期函数可以看作




    T









    T\rightarrow\infty


    T
    ,因此





    Δ


    w


    =


    (


    n


    +


    1


    )



    w


    0






    n



    w


    0



    =



    w


    0



    =




    2


    π



    T






    0



    \Delta w=(n+1)w_0-nw_0=w_0=\frac{2\pi}{T}\rightarrow0


    Δw=(n+1)w0nw0=w0=T2π0
    也就是说,基频率




    w



    w


    w
    由离散变为了连续
    ,周期函数频谱中横坐标离散的




    n



    w


    0




    nw_0


    nw0
    在非周期函数中就变为了连续的




    w



    w


    w
  • 傅里叶级数的复数形式为









    f


    (


    t


    )









    =







    n


    =














    C


    n




    e



    i


    n



    w


    0



    t




    =







    n


    =














    1


    T




    (











    T


    2





    T


    2




    f


    (


    t


    )



    e






    i


    n



    w


    0



    t




    d


    t


    )




    e



    i


    n



    w


    0



    t


















    =







    n


    =















    Δ


    w




    2


    π





    (











    T


    2





    T


    2




    f


    (


    t


    )



    e






    i


    n



    w


    0



    t




    d


    t


    )




    e



    i


    n



    w


    0



    t









    f(t)=n=Cneinw0t=n=1T(T2T2f(t)einw0tdt)einw0t=n=Δw2π(T2T2f(t)einw0tdt)einw0t



    f(t)=n=Cneinw0t=n=T1(2T2Tf(t)einw0tdt)einw0t=n=2πΔw(2T2Tf(t)einw0tdt)einw0t





    T









    T\rightarrow\infty


    T
    时,有














    T


    2





    T


    2




    d


    t


















    +







    d


    t



    n



    w


    0






    w








    n


    =













    Δ


    w


    =















    +







    d


    w



    \int_{-\frac{T}{2}}^{\frac{T}{2}}d t\rightarrow\int_{-\infty}^{+\infty}d t\ nw_0\rightarrow w\ \sum_{n=-\infty}^\infty\Delta w=\int_{-\infty}^{+\infty} d w


    2T2Tdt+dtnw0wn=Δw=+dw
    代入傅里叶级数的复数形式可得









    f


    (


    t


    )









    =



    1



    2


    π

















    +








    (















    +







    f


    (


    t


    )



    e






    i


    w


    t




    d


    t


    )




    e



    i


    w


    t




    d


    w







    f(t)=12π+(+f(t)eiwtdt)eiwtdw



    f(t)=2π1+(+f(t)eiwtdt)eiwtdw
  • 傅里叶变换即为中间括号里的积分项









    F


    (


    w


    )









    =















    +







    f


    (


    t


    )



    e






    i


    w


    t




    d


    t







    F(w)=+f(t)eiwtdt



    F(w)=+f(t)eiwtdt
    傅里叶变换的逆变换即为









    f


    (


    t


    )









    =



    1



    2


    π

















    +







    F


    (


    w


    )



    e



    i


    w


    t




    d


    w







    f(t)=12π+F(w)eiwtdw



    f(t)=2π1+F(w)eiwtdw

参考文献