Fourier 级数入门

第十一周考了期中,感觉裂开(我tcl;任何周期为 2π2\pi 的函数都可以表示为傅里叶级数(一种三角级数),然后就可以将难以积分、求导的函数变化为易于积分的三角级数。

定义1(三角级数)

akR, k=0,1,, bkR, k=1,2,a_k\in \mathbb R,\ k=0,1,\cdots,\ b_k\in\mathbb R,\ k=1,2,\cdots 称级数

a02+k=1+(akcoskx+bksinkx)\frac{a_0}{2}+\sum_{k=1}^{+\infty}(a_k\cos kx+b_k\sin kx)

为三角级数。


EulerEuler 公式知,eiθ=cosθ+isinθe^{i\theta} = \cos \theta+i\sin\theta,则

a02+k=1+(akcoskx+bksinkx)= a02ei0x+k=1+(akeikx+eikx2+ibkeikxeikx2)= a02ei0x+k=1+akibk2eikx+k=1+ak+ibk2eikx\begin{aligned} \frac{a_0}{2}+\sum_{k=1}^{+\infty}(a_k\cos kx+b_k\sin kx) = &\ \frac{a_0}{2}e^{i\cdot 0x}+\sum_{k=1}^{+\infty}\left(a_k\frac{e^{ikx}+e^{-ikx}}{2}+ib_k\frac{e^{-ikx}-e^{ikx}}{2}\right)\\ = &\ \frac{a_0}{2}e^{i\cdot 0x}+\sum_{k=1}^{+\infty}\frac{a_k-ib_k}{2}e^{ikx}+\sum_{k=1}^{+\infty}\frac{a_k+ib_k}{2}e^{-ikx} \end{aligned}

所以,三角级数还可以表示为幂指数的形式。

定义2(f可展开三角级数)

f:RRf:\mathbb R\rightarrow \mathbb R2π2\pi 周期的函数,如果 ak, k=0,1,, bk, k=1,2,\exists a_k, \ k=0,1,\cdots,\ b_k,\ k = 1, 2, \cdots,使得级数

a02+k=1+(akcoskx+bksinkx)\frac{a_0}{2}+\sum_{k=1}^{+\infty}(a_k\cos kx+b_k\sin kx)

收敛,并且

f=a02+k=1+(akcoskx+bksinkx)f=\frac{a_0}{2}+\sum_{k=1}^{+\infty}(a_k\cos kx+b_k\sin kx)

则称 FF 可展成三角级数。


为了证明所有的以 2π2\pi 为周期的函数都能表示为三角级数的形式,所以先对这些系数进行一个推导,然后在进行详细的证明,确定了这些系数后的三角函数就成为 FourierFourier 级数。

下面先对 ff 进行限制,以助于下面的不严谨推导。

f:[a,b]Rf:[a,b]\rightarrow \mathbb R,如果 ff RiemannRiemann 可积或者 ff 有有限个奇点,ff 绝对收敛(可积),则称 ff 可积或者 ff 绝对可积,记 fL1([a,b])f\in L^1([a,b])(这里 LL 指的是 LebesgueLebesgue 可积)

f:RRf:\mathbb R\rightarrow \mathbb R2π2\pi 周期,f[π,π]f\biggl|_{[-\pi,\pi]} 是可积或绝对可积,设

f=a02+k=1+(akcoskx+bksinkx)f=\frac{a_0}{2}+\sum_{k=1}^{+\infty}(a_k\cos kx+b_k\sin kx)

下面对 a0,ak,bka_0,a_k,b_k 的取值进行推导(形式推导):(先计算几个积分)

ππsinkxdx=ππcoskxdx=0k=1,2,ππsinkxcosmxdx=0ππsin2kxdx=ππcos2kxdx=πk=1,2,ππcoskxcosmxdx=ππsinkxsinmxdx=0k,m0\begin{aligned} &\int_{-\pi}^{\pi}\sin kx\,dx = \int_{-\pi}^{\pi}\cos kx\,dx = 0&k = 1, 2,\cdots\\ &\int_{-\pi}^{\pi}\sin kx\cos mx\,dx = 0\\ &\int_{-\pi}^{\pi}\sin^2kx\,dx = \int_{-\pi}^{\pi}\cos^2kx\,dx = \pi&k = 1, 2,\cdots\\ &\int_{-\pi}^{\pi}\cos kx\cos mx\,dx = \int_{-\pi}^{\pi}\sin kx\sin mx\,dx = 0& k,m\neq 0 \end{aligned}

第四条可以理解为一组基的正交性

ππf=πa0a0=1πππfdxππf(x)cosnxdx=ππancos2nxdx=πanan=ππf(x)cosnxdx同理可得,bn=1πππf(x)sinnxdx\begin{aligned} &\int_{-\pi}^{\pi}f = \pi a_0\Rightarrow a_0 = \frac{1}{\pi}\int_{-\pi}^{\pi}f\,dx\\ &\int_{-\pi}^{\pi}f(x)\cos nx\,dx = \int_{-\pi}^{\pi}a_n\cos^2nx\,dx = \pi a_n\Rightarrow a_n = \int_{-\pi}^{\pi}f(x)\cos nx\,dx\\ &\text{同理可得,} b_n = \frac{1}{\pi}\int_{-\pi}^{\pi}f(x)\sin nx \,dx\\ \end{aligned}

定义3(FourierFourier 级数)

f:RRf:\mathbb R\rightarrow \mathbb R2π2\pi 为周期,f[π,π]L1([a,b])f\biggl|_{[-\pi,\pi]}\in L^1([a,b]),定义

ak=1πππf(x)coskxdxk=0,1,2,bk=1πππf(x)sinkxdxk=1,2,\begin{aligned} &a_k=\frac{1}{\pi}\int_{-\pi}^{\pi}f(x)\cos kx \,dx &k=0,1,2,\cdots\\ &b_k=\frac{1}{\pi}\int_{-\pi}^{\pi}f(x)\sin kx \,dx &k=1,2,\cdots \end{aligned}

称级数

a02+k=1+(akcoskx+bksinkx)\frac{a_0}{2}+\sum_{k=1}^{+\infty}(a_k\cos kx+b_k\sin kx)

FourierFourier 级数,称 ak, k=0,1,2,, bk, k=1,2,a_k,\ k = 0, 1, 2,\cdots,\ b_k,\ k = 1, 2,\cdotsffFourierFourier 系数,并记

fa02+k=1+(akcoskx+bksinkx)f\sim\frac{a_0}{2}+\sum_{k=1}^{+\infty}(a_k\cos kx+b_k\sin kx)


例一

f(x):RRf(x):\mathbb R\rightarrow \mathbb R2π2\pi 周期函数,f(x)=cosαx, x[π,π], αZf(x) = \cos \alpha x,\ x \in [-\pi,\pi],\ \alpha \notin\mathbb Z,则

ak= 2π0πf(x)coskxdx=1π0π{cos(α+k)x+cos(αk)x}dx= 1π{sin(α+k)πα+k+sin(αk)παk}=1π(sinαπα+k+sinαπαk)(1)k= 2α(1)ksinαππ(α2k2)bk= 1πππf(x)sinkxdx=1πππcosαxsinkxdx=0\begin{aligned} a_k =&\ \frac{2}{\pi}\int^{\pi}_0f(x)\cos kx\,dx = \frac{1}{\pi}\int_{0}^{\pi}\{\cos(\alpha+k)x+\cos(\alpha -k)x\}\,dx\\ =&\ \frac{1}{\pi}\left\{\frac{\sin(\alpha+k)\pi}{\alpha + k}+\frac{\sin(\alpha - k)\pi}{\alpha -k}\right\}=\frac{1}{\pi}\left(\frac{\sin\alpha\pi}{\alpha+k}+\frac{\sin\alpha\pi}{\alpha -k}\right)(-1)^k\\ =&\ \frac{2\alpha(-1)^k\sin\alpha\pi}{\pi(\alpha^2-k^2)}\\ b_k =&\ \frac{1}{\pi}\int_{-\pi}^{\pi}f(x)\sin kx\,dx = \frac{1}{\pi}\int_{-\pi}^{\pi}\cos \alpha x\sin kx\,dx=0 \end{aligned}


Fourier 级数入门
https://wty-yy.github.io/posts/35053/
作者
wty
发布于
2021年11月26日
许可协议