抱歉,您的浏览器无法访问本站
本页面需要浏览器支持(启用)JavaScript
了解详情 >

傅里叶级数与变换推导

一、三角函数正交性

组成三角级数的函数系

1,cosωt,sinωt,cos2ωt,sin2ωt,……,cosnωt,sinnωt,……1,cos\omega t,sin\omega t,cos2\omega t,sin2\omega t,\ldots\ldots,cosn\omega t,sinn\omega t,\ldots\ldots

在区间[−π,π][-\pi,\pi]或[0,2π][0,2\pi]上正交,即:其中任意两个不同函数之积在区间[−π,π][-\pi,\pi]或[0,2π][0,2\pi]积分为0.

二、周期为2π2\pi的函数展开为傅里叶级数

设周期为2π2\pi的函数f(x)=f(x+2π)f(x)=f(x+2\pi),其可以展开为傅里叶级数,形如:

f(x)=a02+∑n=1∞[ancos⁡(nx)+bnsin⁡(nx)]f(x)=\frac{a_0}{2}+\sum_{n=1}^{\infty}\left[a_n \cos (n x)+b_n \sin (n x)\right]

上式中写作a02\frac{a_0}{2}主要是为了保持表达式的统一性

如何求a0,an,bna_0,a_n,b_n?核心在于利用三角函数的正交性

  • 求a0a_0

    对等式两端在区间[−π,π][-\pi,\pi]上积分

    ∫−ππf(x)dx=∫−ππa02dx+∫−ππ∑n=1πancos⁡nxdx+∫−ππ∑n=1∞bnsin⁡xdx∫−ππf(x)dx=π⋅a0a0=1π∫−ππf(x)dx\begin{aligned} \int_{-\pi}^{\pi}f(x)dx&=\int_{-\pi}^{\pi}\frac{a_{0}}{2}dx+\int_{-\pi}^{\pi}\sum_{n=1}^{\pi}a_{n}\cos nxdx+\int_{-\pi}^{\pi}\sum_{n=1}^{\infty}b_{n}\sin xdx\\ \int_{-\pi}^{\pi}f(x)dx&=\pi \cdot a_0\\ a_0&=\frac{1}{\pi}\int_{-\pi}^{\pi}f(x)dx \end{aligned}

  • 求ana_n

    先同乘cos⁡mx\cos mx,再对等式两端在区间[−π,π][-\pi,\pi]上积分

    ∫−ππf(x)cos⁡mxdx=∫−ππa02cos⁡mxdx+∫−ππ∑n=1∞ancos⁡nx⋅cos⁡mxdx+∫−ππ∑n=1∞bn∫−πsin⁡nx⋅cos⁡mxdx∫−ππf(x)cos⁡mxdx=∫−ππ∑n=1∞ancos⁡nx⋅cos⁡mxdx∫−ππf(x)cos⁡nxdx=∫−ππan(cos⁡nx)2dx=anπan=1π∫−ππf(x)cos⁡nxdx\begin{aligned} \int_{-\pi}^{\pi}f(x)\cos mxdx&=\int_{-\pi}^{\pi}\frac{a_0}{2}\cos mxdx\\ +\int_{-\pi}^{\pi}\sum_{n=1}^{\infty}a_n\cos nx\cdot\cos mxdx & +\int_{-\pi}^{\pi}\sum_{n=1}^{\infty}b_n\int_{-\pi}\sin nx\cdot\cos mxdx \\ \int_{-\pi}^{\pi}f(x)\cos mxdx &=\int_{-\pi}^{\pi}\sum_{n=1}^{\infty}a_n\cos nx\cdot \cos mxdx \\ \int_{-\pi}^{\pi}f(x)\cos nxdx&=\int_{-\pi}^{\pi}a_n(\cos nx)^2dx =a_n\pi\\ a_n&=\frac{1}{\pi}\int_{-\pi}^{\pi}f(x)\cos nxdx \end{aligned}

    上式中第一行到第二行利用了三角函数的正交性,两个积分始终为0。第二行到第三行去掉∑n=1∞\sum_{n=1}^{\infty}则是由于,当且仅当n=mn=m时,积分不为0。

  • 同理可求得bnb_n

    bn=1π∫−ππf(x)sin⁡nxdxb_n=\frac{1}{\pi}\int_{-\pi}^{\pi}f(x)\sin nxdx

综上,周期为2π2\pi的函数f(x)f(x)的傅里叶展开为:

f(x)=a02+∑n=1∞[ancos⁡(nx)+bnsin⁡(nx)]a0=1π∫−ππf(x)dxan=1π∫−ππf(x)cos⁡nxdxbn=1π∫−ππf(x)sin⁡nxdx\begin{aligned} f(x)&=\frac{a_0}{2}+\sum_{n=1}^{\infty}\left[a_n \cos (n x)+b_n \sin (n x)\right]\\ a_0&=\frac{1}{\pi}\int_{-\pi}^{\pi}f(x)dx\\ a_n&=\frac{1}{\pi}\int_{-\pi}^{\pi}f(x)\cos nxdx\\ b_n&=\frac{1}{\pi}\int_{-\pi}^{\pi}f(x)\sin nxdx \end{aligned}

三、周期为2L2L的函数展开为傅里叶级数

设周期为2L2L的函数f(t)=f(t+2L)f(t)=f(t+2L),使用换元法,令x=πLx=\frac{\pi}{L},则t=Lπxt=\frac{L}{\pi}x,此时

x=πLtf(t)=f(Lπx)≜g(x)=g(x+2π)cos⁡nx=cos⁡nπLtsin⁡nx=sin⁡nπLt1π∫−ππdx=1L∫−LLdtx=\frac{\pi}{L}t \\ f(t)=f(\frac{L}{\pi}x)\triangleq g(x)=g(x+2\pi)\\ \cos nx=\cos\frac{n\pi}{L}t \\ \sin nx=\sin\frac{n\pi}{L}t \\ \frac{1}{\pi}\int_{-\pi}^{\pi}dx=\frac{1}{L}\int_{-L}^{L}dt

由此可知,周期为2L2L的函数f(t)f(t)的傅里叶展开为:

f(t)=a02+∑n=1∞ancos⁡nπLt+∑n=1∞bnsin⁡nπLta0=1L∫−LLf(t)dtan=1L∫−LLf(t)cos⁡nπLtdtbn=1L∫−LLf(t)sin⁡nπLtdt\begin{aligned} f(t) & =\frac{a_0}{2}+\sum_{n=1}^\infty a_n\cos\frac{n\pi}{L}t+\sum_{n=1}^\infty b_n\sin\frac{n\pi}{L}t \\ a_0 & =\frac{1}{L}\int_{-L}^Lf(t)dt \\ a_n & =\frac{1}{L}\int_{-L}^Lf(t)\cos\frac{n\pi}{L}tdt \\ b_n & =\frac{1}{L}\int_{-L}^Lf(t)\sin\frac{n\pi}{L}tdt \end{aligned}

实际工程中,t>0t>0,怎么∫−LLf(t)dt\int_{-L}^Lf(t)dt呢?

评论