Fourier Series

最新推荐文章于 2021-07-29 17:54:27 发布

兔龟

最新推荐文章于 2021-07-29 17:54:27 发布

阅读量271

点赞数

分类专栏：信号处理文章标签：数学

本文链接：https://blog.csdn.net/qq_17213753/article/details/78691065

版权

信号处理专栏收录该内容

1 篇文章 0 订阅

订阅专栏

introduction

本文对傅里叶级数的基本思想和推导过程进行介绍。

其基本思想是周期函数可以表示为一系列正(余)弦函数的叠加——对它进行正交分解。因为正(余)弦函数也是周期函数，并且拥有无限可导、表示简单(只需要三个参数)等性质，所以对其进行数学处理比较方便。

lemma

满足<script type="math/tex" id="MathJax-Element-3"><\alpha,\beta>=0</script>的非零向量 $\alpha,\beta$ 正交的，内积空间中两两正交的一组非零向量称为正交组。
在 $\mathbb{R}[a,b]$ 上定义的 <f(x),g(x)>=∫baf(x)g(x)dx <script type="math/tex" id="MathJax-Element-6"> = \int_a^b f(x)g(x)dx</script>满足内积条件，从而构成一个内积空间。

正交组
1. 三角函数系 $\{1,\cos x,\sin x,\cdots,\cos nx,\sin nx\},n\in N_+$ 在任意长度为 $2\pi$ 的区间上两两正交
2. 复指数函数集 $\{e^{jn\omega t}\},n \in Z$ 任意长度为 $2\pi$ 的区间上两两正交

内积是求相似性的运算，而正交意味着两个事物(向量或函数)不相似，因此正交分解会有比较好的数学性质。

construction

考虑对某周期函数进行正弦分解，容易写出其表达式

x (t) = \sum i = 1 \infty A i sin (ω i t + φ i)

$x(t) = \sum\limits_{i=1}^\infty A_i\sin(\omega_i t + \varphi_i)$
将

x(t+T)=x(t) $x(t+T) = x(t)$ 代入式子，易得

ωiT=2kπ,k∈Z $\omega_i T = 2k\pi,k \in Z$ ，即

ωi=kω,ω=2πT $\omega_i = k\omega,\omega = \frac{2\pi}{T}$ 。由于正弦函数为奇函数，所以

ωi<0 $\omega_i < 0$ 可以转化为

ωi>0 $\omega_i > 0$ 的情况；而当

ωk=0 $\omega_k = 0$ 时第

k $k$ 项变为常数，因此可以上式写成

x (t) = c 0 + \sum n = 1 \infty c n sin (n ω t + φ n)

$x(t) = c_0 + \sum\limits_{n=1}^\infty c_n \sin(n\omega t + \varphi_n)$
此即傅里叶级数的三角函数形式，类似地有余弦分解后的函数表达式，由于

cosθ $\cos\theta$ 和

sinθ $\sin\theta$ 只相差

π2 $\frac{\pi}{2}$ 个相位，所以它们实际上是等价的。

derivation

为了得到完整的表达式，需要求出各项系数的值，首先通过二角和公式展开为一般形式

x (t) = c 0 + \sum n = 1 \infty c n [sin (n ω t) cos (φ n) + cos (n ω t) sin (φ n)] = a 0 + \sum n = 1 \infty [a n sin (n ω t) + b n cos (n ω t)]

$\begin{align} x(t) &= c_0 + \sum\limits_{n=1}^\infty c_n[\sin(n\omega t)\cos(\varphi_n) + \cos(n\omega t )\sin(\varphi_n)]\\ &=a_0 + \sum\limits_{n=1}^\infty[a_n\sin(n\omega t) + b_n\cos(n\omega t )] \end{align}$
其中

a0=c0,an=cncos(φn),bn=cnsin(φn) $a_0 = c_0,a_n = c_n\cos(\varphi_n),b_n = c_n\sin(\varphi_n)$ ，然后利用三角函数系的正交性，令

u=ωt $u = \omega t$ 构造

\int π - π x (u ω) sin (n u) d (u ω) = 1 ω \int π - π a n sin 2 (n u) d u = a n π ω

$\int_{-\pi}^{\pi}x(\frac{u}{\omega})\sin(nu)d(\frac{u}{\omega}) = \frac{1}{\omega}\int_{-\pi}^{\pi}a_n\sin^2(nu)du = \frac{a_n\pi}{\omega}$
由于

u∈[−π,π],T=2πω $u \in [-\pi,\pi],T= \frac{2\pi}{\omega}$ ，那么

t∈[−T2,T2] $t \in [-\frac{T}{2},\frac{T}{2}]$ 代回原式求得

a n = 2 T \int T 2 - T 2 x (t) sin (n ω t) d t

$a_n = \frac{2}{T}\int_{-\frac{T}{2}}^{\frac{T}{2}}x(t)\sin(n\omega t)dt$
同理有

bn=2T∫T2−T2x(t)cos(nωt)dt,c0=1T∫T2−T2x(t)dt $b_n = \frac{2}{T} \int_{-\frac{T}{2}}^{\frac{T}{2}}x(t)\cos(n\omega t)dt,c_0 = \frac{1}{T}\int_{-\frac{T}{2}}^{\frac{T}{2}}x(t)dt$ ，另外使

c0=a0 $c_0 = a_0$ ，它们合称三角傅里叶系数。由于三角函数系在任意长度为

2π $2\pi$ 的区间上的正交性，容易证明三角傅里叶系数可以表示为

a 0 = 1 T \int T x (t) d t, a n = 2 T \int T x (t) sin (n ω t) d t, b n = 2 T \int T x (t) cos (n ω t) d t

$a_0 = \frac{1}{T}\int_Tx(t)dt,a_n = \frac{2}{T}\int_Tx(t)\sin(n\omega t)dt,b_n = \frac{2}{T} \int_T x(t)\cos(n\omega t)dt$

至此，我们验证了上述分解的可行性，同时也知道了如果想求各项系数，只需取某一周期段的原函数与三角函数系相乘再进行积分即可。

在进行分解之后，可以通过频谱和相位谱对原始函数进行观察和分析，当然它们都只反映了一部分原始函数的信息。

通过引入复指数 $e^{j\theta} = \cos\theta + j\sin\theta$ ，有 $\sin\theta = \frac{e^{j\theta}-e^{-j\theta}}{2j}$ ，可以得到更优雅的傅里叶级数的形式

x (t) = c 0 + \sum n = 1 \infty c n 2 j [e j (n ω t + φ n) - e - j (n ω t + φ n)] = c 0 + \sum n = - \infty, n \neq 0 \infty c n 2 j e j n ω t e j φ n = \sum n = - \infty \infty F n e j n ω t

$\begin{align} x(t) &= c_0 + \sum\limits_{n=1}^\infty \frac{c_n}{2j} [e^{j(n\omega t + \varphi_n)}-e^{-j(n\omega t + \varphi_n)}] \\ &= c_0 + \sum\limits_{n=-\infty,n \neq 0}^\infty \frac{c_n}{2j} e^{jn\omega t} e^{j\varphi_n} \\ &= \sum\limits_{n=-\infty}^\infty F_ne^{jn\omega t}\end{align}$
此即傅里叶级数的指数形式，其中