历史
定义
上图展示了蓝色无周期函数
s
(
x
)
{\displaystyle s(x)}
,它只定义在红色区间
[
0
,
P
]
{\displaystyle [0,P]}
之上。可以认为傅里叶级数,是对下图中最初函数的“周期延拓”的分析,傅里叶级数总是周期函数,即使最初函数
s
(
x
)
{\displaystyle s(x)}
不是周期函数。
傅里叶级数可以用不同的形式来表达,下面将周期 为
P
{\displaystyle P}
的一个周期函数
s
(
x
)
,
x
∈
R
{\textstyle s(x),\ x\in \mathbb {R} }
表达为不同形式的傅里叶级数。
正弦-余弦形式
人们常用
sin
(
x
)
{\displaystyle \sin(x)}
与
cos
(
x
)
{\displaystyle \cos(x)}
的三角级数 来表示
s
(
x
)
{\textstyle s(x)}
,就是将所有
n
{\displaystyle n}
阶谐波
sin
(
2
π
n
x
P
)
{\textstyle \sin({\frac {2\pi nx}{P}})}
与
cos
(
2
π
n
x
P
)
{\textstyle \cos({\frac {2\pi nx}{P}})}
,乘以其各自在
s
(
x
)
{\textstyle s(x)}
中的权重 ,求得它们的总和 ;这些
n
{\displaystyle n}
阶谐波的权重称为傅里叶级数系数,它们可以借由如下积分来获得:
傅里叶级数系数
A
0
=
1
P
∫
P
s
(
x
)
d
x
A
n
=
2
P
∫
P
s
(
x
)
cos
(
2
π
n
x
P
)
d
x
for
n
≥
1
B
n
=
2
P
∫
P
s
(
x
)
sin
(
2
π
n
x
P
)
d
x
for
n
≥
1
{\displaystyle {\begin{aligned}A_{0}&={\frac {1}{P}}\int _{P}s(x)\,dx\\A_{n}&={\frac {2}{P}}\int _{P}s(x)\cos \left({\frac {2\pi nx}{P}}\right)\,dx\qquad {\text{for }}n\geq 1\qquad \\B_{n}&={\frac {2}{P}}\int _{P}s(x)\sin \left({\frac {2\pi nx}{P}}\right)\,dx\qquad {\text{for }}n\geq 1\end{aligned}}}
Eq.1
符号
∫
P
{\textstyle \int _{P}}
表示在选定区间 上的积分,典型的选择为
[
−
P
/
2
,
P
/
2
]
{\displaystyle [-P/2,P/2]}
或者
[
0
,
P
]
{\displaystyle [0,P]}
。注意
A
0
{\displaystyle A_{0}}
是函数
s
(
x
)
{\displaystyle s(x)}
的平均值 [ A] ,这个性质扩展到了类似的变换比如傅里叶变换 。
通过这些系数定义傅里叶级数为:
傅里叶级数,正弦-余弦形式
s
(
x
)
∼
A
0
+
∑
n
=
1
∞
(
A
n
cos
(
2
π
n
x
P
)
+
B
n
sin
(
2
π
n
x
P
)
)
{\displaystyle s(x)\sim A_{0}+\sum _{n=1}^{\infty }\left(A_{n}\cos \left({\frac {2\pi nx}{P}}\right)+B_{n}\sin \left({\frac {2\pi nx}{P}}\right)\right)}
Eq.2
这里使用符号
∼
{\displaystyle \sim }
,表示傅里叶级数的求和不一定总是等于
s
(
x
)
{\displaystyle s(x)}
。普遍来说
n
{\displaystyle n}
是理论上趋近于无限大的,但是就算趋近于无限大,对所有的
x
{\displaystyle x}
(例如在某一点上不连续),傅里叶级数也不一定收敛到
s
(
x
)
{\textstyle s(x)}
。尽管不收敛的可能性始终存在,在科学和工程领域中经常将Eq. 2 中的
∼
{\displaystyle \sim }
直接替代为
=
{\displaystyle =}
。
在傅里叶级数系数中的整数索引
n
{\displaystyle n}
,是级数中相应的
cos
{\displaystyle \cos }
或
sin
{\displaystyle \sin }
,在这个函数的周期
P
{\displaystyle P}
中,形成的圆周 (cycle)的数目。因此对应于
A
n
{\displaystyle A_{n}}
和
B
n
{\displaystyle B_{n}}
的项有着:
波长 等于
P
n
{\displaystyle {\tfrac {P}{n}}}
,并且有着同于
x
{\displaystyle x}
的单位。
频率 等于
n
P
{\displaystyle {\tfrac {n}{P}}}
,并且有着
x
{\displaystyle x}
的倒数单位。
指数形式
下面借由欧拉公式
e
i
x
=
cos
x
+
i
sin
x
{\displaystyle \ e^{ix}=\cos x+i\sin x\ }
,将傅里叶级数系数简化成复数 指数 形式。
根据定义,我们可以得到:
复数傅里叶级数系数
c
0
=
A
0
c
n
=
(
A
n
−
i
B
n
)
/
2
for
n
>
0
c
n
=
(
A
−
n
+
i
B
−
n
)
/
2
for
n
<
0
{\displaystyle {\begin{aligned}c_{0}&=A_{0}&\\c_{n}&=(A_{n}-iB_{n})/2\qquad &{\text{for }}n>0\\c_{n}&=(A_{-n}+iB_{-n})/2\qquad &{\text{for }}n<0\end{aligned}}}
Eq. 3
通过将等式Eq. 1 代入Eq. 3 ,可以证实[ 4] :
复数傅里叶级数系数
c
n
=
1
P
∫
P
s
(
x
)
e
−
2
π
i
n
x
P
d
x
for
n
∈
Z
{\displaystyle c_{n}={\frac {1}{P}}\int _{P}s(x)e^{-{\frac {2\pi inx}{P}}}\,dx\qquad {\text{for}}\ n\in \mathbb {Z} }
给定复数傅里叶级数系数,可以用公式复原出
A
n
{\displaystyle A_{n}}
和
B
n
{\displaystyle B_{n}}
:
复数傅里叶级数系数
A
0
=
c
0
A
n
=
c
n
+
c
−
n
for
n
>
0
B
n
=
i
(
c
n
−
c
−
n
)
for
n
>
0
{\displaystyle {\begin{aligned}A_{0}&=c_{0}&\\A_{n}&=c_{n}+c_{-n}\qquad &{\textrm {for}}~n>0\\B_{n}&=i(c_{n}-c_{-n})\qquad &{\textrm {for}}~n>0\end{aligned}}}
通过这些定义,傅里叶级数可以写为:
傅里叶级数,指数形式
s
(
x
)
∼
∑
n
=
−
∞
∞
c
n
⋅
e
2
π
i
n
x
P
{\displaystyle s(x)\sim \sum _{n=-\infty }^{\infty }c_{n}\cdot e^{\frac {2\pi inx}{P}}}
Eq. 4
这是可推广到复数值域 函数的惯用形式。
n
{\displaystyle n}
的负数值对应于负频率 。
复数值函数
人们习惯将
s
(
x
)
{\textstyle s(x)}
的值域 普遍化到复数 上,设
s
(
x
)
{\textstyle s(x)}
是一个复数值函数,它的实部和虚部,都是实数值函数:
s
(
x
)
=
Re
(
s
(
x
)
)
+
i
⋅
Im
(
s
(
x
)
)
,
x
∈
R
{\displaystyle s(x)=\operatorname {Re} (s(x))+i\cdot \operatorname {Im} (s(x)),\quad x\in \mathbb {R} }
定义
c
n
≜
c
R
n
+
i
⋅
c
I
n
{\displaystyle c_{n}\triangleq c_{_{Rn}}+i\cdot c_{_{In}}}
则:
c
n
=
1
P
∫
P
s
(
x
)
⋅
e
−
i
2
π
p
n
x
d
x
=
1
P
∫
P
Re
(
s
(
x
)
)
⋅
e
−
i
2
π
p
n
x
d
x
+
i
⋅
1
P
∫
P
Im
(
s
(
x
)
)
⋅
e
−
i
2
π
p
n
x
d
x
{\displaystyle c_{n}={\frac {1}{P}}\int _{P}s(x)\cdot e^{-i{\tfrac {2\pi }{p}}nx}\ dx={\frac {1}{P}}\int _{P}\operatorname {Re} (s(x))\cdot e^{-i{\tfrac {2\pi }{p}}nx}\ dx+i\cdot {\frac {1}{P}}\int _{P}\operatorname {Im} (s(x))\cdot e^{-i{\tfrac {2\pi }{p}}nx}\ dx}
c
R
n
=
1
P
∫
P
Re
(
s
(
x
)
)
⋅
e
−
i
2
π
p
n
x
d
x
{\displaystyle c_{_{Rn}}={\frac {1}{P}}\int _{P}\operatorname {Re} (s(x))\cdot e^{-i{\tfrac {2\pi }{p}}nx}\ dx}
c
I
n
=
1
P
∫
P
Im
(
s
(
x
)
)
⋅
e
−
i
2
π
p
n
x
d
x
{\displaystyle c_{_{In}}={\frac {1}{P}}\int _{P}\operatorname {Im} (s(x))\cdot e^{-i{\tfrac {2\pi }{p}}nx}\ dx}
对于这个复数值函数,它的傅里叶级数的实部,是它的实部的傅里叶级数;它的傅里叶级数的虚部,是它的虚部的傅里叶级数:
s
(
x
)
∼
∑
n
=
−
∞
∞
(
c
R
n
+
i
⋅
c
I
n
)
⋅
e
i
2
π
p
n
x
=
∑
n
=
∞
∞
c
R
n
⋅
e
i
2
π
p
n
x
+
i
⋅
∑
n
=
−
∞
∞
c
I
n
⋅
e
i
2
π
p
n
x
{\displaystyle s(x)\sim \sum _{n=-\infty }^{\infty }\left(c_{_{Rn}}+i\cdot c_{_{In}}\right)\cdot e^{i{\tfrac {2\pi }{p}}nx}=\sum _{n=\infty }^{\infty }c_{_{Rn}}\cdot e^{i{\tfrac {2\pi }{p}}nx}+i\cdot \sum _{n=-\infty }^{\infty }c_{_{In}}\cdot e^{i{\tfrac {2\pi }{p}}nx}}
振幅-相位形式
还可以利用三角恒等式
cos
(
α
−
β
)
=
cos
α
cos
β
+
sin
α
sin
β
{\displaystyle \ \cos(\alpha -\beta )=\cos \alpha \cos \beta +\sin \alpha \sin \beta \,}
,把正弦-余弦形式中后面的正弦函数跟余弦函数合并起来:
A
n
⋅
cos
(
2
π
n
P
x
−
φ
n
)
≡
A
n
cos
(
φ
n
)
⏟
a
n
⋅
cos
(
2
π
n
P
x
)
+
A
n
sin
(
φ
n
)
⏟
b
n
⋅
sin
(
2
π
n
P
x
)
{\displaystyle A_{n}\cdot \cos \left(2\pi {\tfrac {n}{P}}x-\varphi _{n}\right)\ \equiv \ \underbrace {A_{n}\cos(\varphi _{n})} _{a_{n}}\cdot \cos \left(2\pi {\tfrac {n}{P}}x\right)+\underbrace {A_{n}\sin(\varphi _{n})} _{b_{n}}\cdot \sin \left(2\pi {\tfrac {n}{P}}x\right)}
然后定义振幅
A
n
≜
a
n
2
+
b
n
2
{\textstyle A_{n}\triangleq {\sqrt {a_{n}^{2}+b_{n}^{2}}}}
,相位
φ
n
≜
arctan2
(
b
n
,
a
n
)
{\textstyle \varphi _{n}\triangleq \operatorname {arctan2} (b_{n},a_{n})}
,这里的
a
n
{\displaystyle a_{n}}
和
b
n
{\displaystyle b_{n}}
对应正弦-余弦形式中
A
n
{\displaystyle A_{n}}
和
B
n
{\displaystyle B_{n}}
。
A
0
2
{\displaystyle {\tfrac {A_{0}}{2}}}
是
s
(
x
)
{\displaystyle s(x)}
的平均值
1
P
∫
P
s
(
x
)
d
x
{\textstyle {\frac {1}{P}}\int _{P}s(x)\,dx}
。
傅里叶级数,振幅-相位形式
s
(
x
)
∼
A
0
2
+
∑
n
=
1
∞
A
n
⋅
cos
(
2
π
n
P
x
−
φ
n
)
{\displaystyle s(x)\sim {\frac {A_{0}}{2}}+\sum _{n=1}^{\infty }A_{n}\cdot \cos \left(2\pi {\tfrac {n}{P}}x-\varphi _{n}\right)}
Eq. 5
部分求和算子
在描述傅里叶级数行为的时候,经常会为一个函数
f
(
x
)
{\displaystyle f(x)}
介入部分求和 算子
S
N
{\displaystyle S_{N}}
[ 5] :
S
N
(
f
)
=
∑
n
=
−
N
N
c
n
e
2
π
i
n
x
P
{\displaystyle S_{N}(f)=\sum _{n=-N}^{N}c_{n}e^{\frac {2\pi inx}{P}}}
Eq. 6
这里的
c
n
{\displaystyle c_{n}}
是
f
{\displaystyle f}
的傅里叶系数。不同于微积分中的级数,傅里叶级数的部分求和必须采用对称形式,否则收敛结果可能不成立。
假设
f
(
x
)
{\displaystyle f(x)}
与
g
(
x
)
{\displaystyle g(x)}
是在
R
{\textstyle \mathbb {R} }
上的可积函数,
f
(
x
)
{\displaystyle f(x)}
与
g
(
x
)
{\displaystyle g(x)}
在
[
−
π
,
π
]
{\displaystyle [-\pi ,\pi ]}
的卷积
(
f
∗
g
)
(
x
)
{\displaystyle (f*g)(x)}
为:
(
f
∗
g
)
(
x
)
=
∫
−
π
π
f
(
τ
)
g
(
x
−
τ
)
d
τ
{\displaystyle (f*g)(x)=\int _{-\pi }^{\pi }f(\tau )g(x-\tau )d\tau }
周期为
2
π
{\displaystyle 2\pi }
的函数
f
(
x
)
{\displaystyle f(x)}
的傅里叶级数的部分求和 ,可以经由
f
(
x
)
{\displaystyle f(x)}
与狄利克雷核
D
n
(
x
)
=
∑
k
=
−
n
n
e
i
k
x
{\textstyle D_{n}(x)=\sum _{k=-n}^{n}e^{ikx}}
的卷积 来表示:
S
N
(
f
)
(
x
)
=
∑
n
=
−
N
N
c
n
e
i
n
x
=
∑
n
=
−
N
N
(
1
2
π
∫
−
π
π
f
(
τ
)
e
−
i
n
τ
d
τ
)
⋅
e
i
n
x
=
1
2
π
∫
−
π
π
f
(
τ
)
(
∑
n
=
−
N
N
e
i
n
(
x
−
τ
)
)
d
τ
=
1
2
π
(
f
∗
D
N
)
(
x
)
{\displaystyle {\begin{aligned}S_{N}(f)(x)&=\sum _{n=-N}^{N}c_{n}e^{inx}\\&=\sum _{n=-N}^{N}\left({\frac {1}{2\pi }}\int _{-\pi }^{\pi }f(\tau )e^{-in\tau }d\tau \right)\cdot e^{inx}={\frac {1}{2\pi }}\int _{-\pi }^{\pi }f(\tau )\left(\sum _{n=-N}^{N}e^{in(x-\tau )}\right)d\tau \\&={\frac {1}{2\pi }}(f*D_{N})(x)\end{aligned}}}
收敛性概要
s
N
(
x
)
{\displaystyle s_{N}(x)}
在
[
x
0
,
x
0
+
P
]
{\displaystyle [x_{0},\ x_{0}+P]}
近似了
s
(
x
)
{\displaystyle s(x)}
,该近似程度会随着
N
→
∞
{\displaystyle N\rightarrow \infty }
逐渐改善。这个无穷和
s
∞
(
x
)
{\displaystyle s_{\infty }(x)}
叫做
s
{\displaystyle s}
的傅里叶级数表示。傅里叶级数的收敛性取决于函数有限数量的极大值和极小值,这就是通常称为傅里叶级数的狄利克雷条件 。参见傅里叶级数的收敛性 之一。对于广义函数或分布也可以用范数或弱收敛 定义傅里叶系数。在
s
(
x
)
{\displaystyle s(x)}
的不可导 点上,如果我们只取无穷级数中的有限项求和,那么在这些点上会有幅度不随
N
{\displaystyle N}
增大而持续变小的起伏,这叫做吉布斯现象 ,一个简单的例子是方波信号 。
在工程 应用中,一般假定傅里叶级数除了在不连续点以外处处收敛,原因是工程上遇到的函数比数学家提供的这个假定的反例表现更加良好。特别地,傅里叶级数绝对收敛 且一致收敛 于
s
(
x
)
{\displaystyle s(x)}
,只要在
s
(
x
)
{\displaystyle s(x)}
的导数(或许不会处处存在)是平方可积的[ 6] 。如果一个函数在区间
[
x
0
,
x
0
+
P
]
{\displaystyle [x_{0},x_{0}+P]}
上是平方可积 的,那么此傅里叶级数在几乎处处 的点都收敛于该函数。
一个相同幅度和频率的锯齿波的近似的可视化
另一个分别采用傅里叶级数的前 1, 2, 3, 4 项近似方波的可视化。(可以在这里
[ 7] 看到一个交互式的动画)
其他常用表示法
符号
c
n
{\displaystyle c_{n}}
在讨论多个不同函数的傅里叶系数时是不够用的。因此习惯上将其替代为函数(这里是函数
s
{\displaystyle s}
)的某种修改形式,即采用函数式符号比如
s
^
[
n
]
{\displaystyle {\hat {s}}[n]}
或
S
[
n
]
{\displaystyle S[n]}
,来替代下标式符号:
s
(
x
)
=
∑
n
=
−
∞
∞
s
^
(
n
)
⋅
e
2
π
i
n
x
/
P
{\displaystyle s(x)=\sum _{n=-\infty }^{\infty }{\hat {s}}(n)\cdot e^{2\pi inx/P}\quad }
常用的数学符号
s
(
x
)
=
∑
n
=
−
∞
∞
S
[
n
]
⋅
e
i
2
π
n
x
/
P
{\displaystyle s(x)=\sum _{n=-\infty }^{\infty }S[n]\cdot e^{i\,2\pi nx/P}\quad }
常用的工程符号
在工程上,特别是在变量
x
{\displaystyle x}
表示时间的时候,系数序列叫做频域 表示。经常使用方括号来强调这个函数的定义域是频率 的离散集合。
另一个常用频域表示,使用傅里叶级数系数,调制 像梳子一样的狄拉克采样函数 :
S
(
f
)
≜
∑
n
=
−
∞
∞
S
[
n
]
⋅
δ
(
f
−
n
P
)
{\displaystyle S(f)\ \triangleq \ \sum _{n=-\infty }^{\infty }S[n]\cdot \delta \left(f-{\frac {n}{P}}\right)}
这里的
f
{\displaystyle f}
表示连续频域。在变量
x
{\displaystyle x}
以秒为单位的时候,
f
{\displaystyle f}
以赫兹 为单位。采样的间隔为基本频率
1
P
{\displaystyle {\tfrac {1}{P}}}
的
n
{\displaystyle n}
倍(即为谐波 )。
s
∞
(
x
)
{\displaystyle s_{\infty }(x)}
可以通过逆傅里叶变换 从这种表示恢复出来:
F
−
1
{
S
(
f
)
}
=
∫
−
∞
∞
(
∑
n
=
−
∞
∞
S
[
n
]
⋅
δ
(
f
−
n
P
)
)
e
i
2
π
f
x
d
f
=
∑
n
=
−
∞
∞
S
[
n
]
⋅
∫
−
∞
∞
δ
(
f
−
n
P
)
e
i
2
π
f
x
d
f
=
∑
n
=
−
∞
∞
S
[
n
]
⋅
e
i
2
π
n
x
/
P
≜
s
∞
(
x
)
{\displaystyle {\begin{aligned}{\mathcal {F}}^{-1}\{S(f)\}&=\int _{-\infty }^{\infty }\left(\sum _{n=-\infty }^{\infty }S[n]\cdot \delta \left(f-{\frac {n}{P}}\right)\right)e^{i2\pi fx}\,df\\&=\sum _{n=-\infty }^{\infty }S[n]\cdot \int _{-\infty }^{\infty }\delta \left(f-{\frac {n}{P}}\right)e^{i2\pi fx}\,df\\&=\sum _{n=-\infty }^{\infty }S[n]\cdot e^{i\,2\pi nx/P}\ \ \triangleq \ s_{\infty }(x)\end{aligned}}}
构造出的函数
S
(
f
)
{\displaystyle S(f)}
,因而通常称为“傅里叶变换 ”,即使一个周期函数的傅里叶积分在这个谐波频率上不收敛[ B] 。
常用的傅里叶级数
基本性质
对称性质
范例
一个简单的傅里叶级数
锯齿波 周期函数的图
前五个部分傅里叶级数的动态图
我们现在用上面的公式给出一个简单函数的傅里叶级数展开式。考虑一个锯齿波:
s
(
x
)
=
x
π
,
f
o
r
−
π
<
x
<
π
{\displaystyle s(x)={\frac {x}{\pi }},\quad \mathrm {for} -\pi <x<\pi }
s
(
x
+
2
π
k
)
=
s
(
x
)
,
f
o
r
−
∞
<
x
<
∞
and
k
∈
Z
{\displaystyle s(x+2\pi k)=s(x),\quad \mathrm {for} -\infty <x<\infty {\text{ and }}k\in \mathbb {Z} }
在这种情况下,傅里叶级数为:
A
n
=
1
π
∫
−
π
π
s
(
x
)
cos
(
n
x
)
d
x
=
0
,
n
≥
0
B
n
=
1
π
∫
−
π
π
s
(
x
)
sin
(
n
x
)
d
x
=
−
2
π
n
cos
(
n
π
)
+
2
π
2
n
2
sin
(
n
π
)
=
2
(
−
1
)
n
+
1
π
n
,
n
≥
1
{\displaystyle {\begin{aligned}A_{n}&={\frac {1}{\pi }}\int _{-\pi }^{\pi }s(x)\cos(nx)\,dx=0,\quad n\geq 0\\[4pt]B_{n}&={\frac {1}{\pi }}\int _{-\pi }^{\pi }s(x)\sin(nx)\,dx\\[4pt]&=-{\frac {2}{\pi n}}\cos(n\pi )+{\frac {2}{\pi ^{2}n^{2}}}\sin(n\pi )\\[4pt]&={\frac {2\,(-1)^{n+1}}{\pi n}},\quad n\geq 1\end{aligned}}}
可以证明,当
s
{\displaystyle s}
可微时,傅立叶级数在每个点
x
{\displaystyle x}
都收敛于
s
(
x
)
{\displaystyle s(x)}
,于是:
s
(
x
)
=
A
0
+
∑
n
=
1
∞
(
A
n
cos
(
n
x
)
+
B
n
sin
(
n
x
)
)
=
2
π
∑
n
=
1
∞
(
−
1
)
n
+
1
n
sin
(
n
x
)
,
f
o
r
x
−
π
∉
2
π
Z
{\displaystyle {\begin{aligned}s(x)&=A_{0}+\sum _{n=1}^{\infty }\left(A_{n}\cos \left(nx\right)+B_{n}\sin \left(nx\right)\right)\\[4pt]&={\frac {2}{\pi }}\sum _{n=1}^{\infty }{\frac {(-1)^{n+1}}{n}}\sin(nx),\quad \mathrm {for} \quad x-\pi \notin 2\pi \mathbb {Z} \end{aligned}}}
Eq.1
当
x
=
π
{\displaystyle x=\pi }
时,傅里叶级数收敛于
0
{\displaystyle 0}
,为在
x
=
π
{\displaystyle x=\pi }
处
s
{\displaystyle s}
的左极限和右极限之和的一半。这是傅里叶级数的狄利克雷定理 的特例。
这个例子为我们引出了巴塞尔问题 的一种解法。
傅里叶诱导
金属板内的热分布,使用傅里叶方法求解
在上例中我们的函数的傅里叶级数展开式看起来不比
s
(
x
)
=
x
π
{\displaystyle s(x)={\tfrac {x}{\pi }}}
简单,因此人们需要傅里叶级数的原因也就不会立即显现出来。但还有很多应用,我们举用傅里叶诱导解热方程 的例子。考虑边长为
π
{\displaystyle \pi }
米的方形金属版,坐标为
(
x
,
y
)
∈
[
0
,
π
]
×
[
0
,
π
]
{\displaystyle (x,y)\in [0,\pi ]\times [0,\pi ]}
。如果板内没有热源,并且四个边中三个都保持在
0
{\displaystyle 0}
摄氏度,而第四条边
y
=
π
{\displaystyle y=\pi }
,对于
x
∈
(
0
,
π
)
{\displaystyle x\in (0,\pi )}
,保持在温度梯度
T
(
x
,
π
)
=
x
{\displaystyle T(x,\pi )=x}
摄氏度。在这种情况下,稳态(或者说很长时间过后的)热分布函数
T
(
x
,
y
)
{\displaystyle T(x,y)}
不能得出解析解 ,但却可以证明:
T
(
x
,
y
)
=
2
∑
n
=
1
∞
(
−
1
)
n
+
1
n
sin
(
n
x
)
sinh
(
n
y
)
sinh
(
n
π
)
{\displaystyle T(x,y)=2\sum _{n=1}^{\infty }{\frac {(-1)^{n+1}}{n}}\sin(nx){\sinh(ny) \over \sinh(n\pi )}}
这里的
sinh
{\displaystyle \sinh }
是双曲正弦 函数。热方程的这个解是通过将
π
s
(
x
)
{\displaystyle \pi s(x)}
的傅里叶级数的每一项乘以
sinh
(
n
y
)
sinh
(
n
π
)
{\displaystyle {\tfrac {\sinh(ny)}{\sinh(n\pi )}}}
得到的。尽管示例的函数
s
(
x
)
{\displaystyle s(x)}
的傅里叶级数似乎很复杂,用傅里叶的方法却可以求解这个热分布问题。
其他例子
我们也可以应用傅里叶级数去证明等周不等式 ,或是构造处处连续处处不可微的函数。
收敛性
至今还没有判断傅里叶级数的收敛性充分必要条件,但是对于实际问题中出现的函数,有很多种判别条件可用于判断收敛性。比如
x
(
t
)
{\displaystyle x(t)}
的可微性或级数的一致收敛性 。在闭区间上满足狄利克雷条件 的函数表示成的傅里叶级数都收敛。狄利克雷条件如下:
在定义区间上,
x
(
t
)
{\displaystyle x(t)}
须绝对可积 ;
在任一有限区间中,
x
(
t
)
{\displaystyle x(t)}
只能取有限个极值点;
在任何有限区间上,
x
(
t
)
{\displaystyle x(t)}
只能有有限个第一类间断点 。
满足以上条件的
x
(
t
)
{\displaystyle x(t)}
傅里叶级数都收敛,且:
1.当
t
{\displaystyle t}
是
x
(
t
)
{\displaystyle x(t)}
的连续点时,级数收敛于
x
(
t
)
{\displaystyle x(t)}
;
2.当
t
{\displaystyle t}
是
x
(
t
)
{\displaystyle x(t)}
的间断点时,级数收敛于
1
2
[
x
(
t
−
)
+
x
(
t
+
)
]
{\displaystyle {\frac {1}{2}}[x(t^{-})+x(t^{+})]}
。
1966年,里纳特·卡尔松 证明了勒贝格二次可积 函数的傅立叶级数一定是几乎处处收敛的,即级数在除了一个勒贝格零测集外均收敛。
傅里叶级数收敛证明
假设一个函数在
f
(
x
)
{\displaystyle f(x)}
在
[
0
,
2
π
]
{\displaystyle [0,2\pi ]}
上是平方可积,则会有:
1
2
π
∫
0
2
π
|
f
(
x
)
−
S
N
(
f
)
(
x
)
|
2
d
x
→
0
{\displaystyle {\frac {1}{2\pi }}\int _{0}^{2\pi }|f(x)-S_{N}(f)(x)|^{2}dx\rightarrow 0}
当
N
→
∞
{\displaystyle N\rightarrow \infty }
证明的第一步:
考虑一系列正交基底,
{
e
n
}
n
∈
Z
{\displaystyle \{e_{n}\}_{n\in \mathbb {Z} }}
,其中
e
n
(
x
)
=
e
−
i
n
x
{\displaystyle e_{n}(x)=e^{-inx}}
,且有
(
e
n
,
e
m
)
=
{
1
,
if
n
=
m
0
,
if
n
≠
m
{\displaystyle (e_{n},e_{m})={\begin{cases}1,&{\text{if }}n=m\\0,&{\text{if }}n\neq m\end{cases}}}
然后有
(
f
,
e
n
)
=
1
2
π
∫
0
2
π
f
(
x
)
e
−
i
n
x
d
x
=
f
^
(
n
)
{\displaystyle (f,e_{n})={\frac {1}{2\pi }}\int _{0}^{2\pi }f(x)e^{-inx}dx={\hat {f}}(n)}
特别的有,
f
(
x
)
{\displaystyle f(x)}
的傅里叶级数的部分和
S
N
(
f
)
(
x
)
=
∑
|
n
|
≤
N
f
^
(
n
)
e
n
{\displaystyle S_{N}(f)(x)=\sum _{|n|\leq N}{\hat {f}}(n)e_{n}}
然后根据
f
=
f
−
∑
|
n
|
≤
N
f
^
(
n
)
e
n
+
∑
|
n
|
≤
N
f
^
(
n
)
e
n
{\displaystyle f=f-\sum _{|n|\leq N}{\hat {f}}(n)e_{n}+\sum _{|n|\leq N}{\hat {f}}(n)e_{n}}
以及勾股定理,可以有:
|
|
f
|
|
2
=
|
|
f
−
∑
|
n
|
≤
N
f
^
(
n
)
e
n
|
|
2
+
|
|
∑
|
n
|
≤
N
f
^
(
n
)
e
n
|
|
2
{\displaystyle ||f||^{2}=||f-\sum _{|n|\leq N}{\hat {f}}(n)e_{n}||^{2}+||\sum _{|n|\leq N}{\hat {f}}(n)e_{n}||^{2}}
替换一下后有
|
|
f
|
|
2
=
|
|
f
−
S
N
(
f
)
(
x
)
|
|
2
+
|
|
∑
|
n
|
≤
N
f
^
(
n
)
e
n
|
|
2
{\displaystyle ||f||^{2}=||f-S_{N}(f)(x)||^{2}+||\sum _{|n|\leq N}{\hat {f}}(n)e_{n}||^{2}}
如果右边第一项收敛到0,再根据正交的性质,可以看出上述式子中的右手边第二项:
|
|
∑
|
n
|
≤
N
f
^
(
n
)
e
n
|
|
2
=
∑
|
n
|
≤
N
|
f
^
(
n
)
|
2
{\displaystyle ||\sum _{|n|\leq N}{\hat {f}}(n)e_{n}||^{2}=\sum _{|n|\leq N}|{\hat {f}}(n)|^{2}}
,这就证明了帕塞瓦尔定理 。
证明的第二步:
回到证明右边第一项,因为函数
f
(
x
)
{\displaystyle f(x)}
可积,找到一个连续函数
g
(
x
)
{\displaystyle g(x)}
,然后根据最佳逼近引理,可以找到一个三角多项式p(x),使得
|
f
−
S
N
(
f
)
(
x
)
|
≤
|
f
(
x
)
−
g
(
x
)
|
+
|
g
(
x
)
−
S
N
(
f
)
(
x
)
|
{\displaystyle |f-S_{N}(f)(x)|\leq |f(x)-g(x)|+|g(x)-S_{N}(f)(x)|}
故当
N
→
∞
{\displaystyle N\rightarrow \infty }
,函数
f
(
x
)
{\displaystyle f(x)}
跟
S
N
(
f
)
(
x
)
{\displaystyle S_{N}(f)(x)}
的差为0。
其他性质
傅里叶级数的唯一性
如果有一个定义在
[
−
π
,
π
]
{\displaystyle [-\pi ,\pi ]}
的函数
f
(
x
)
{\displaystyle f(x)}
和
g
(
x
)
{\displaystyle g(x)}
,其中函数
f
(
x
)
{\displaystyle f(x)}
和
g
(
x
)
{\displaystyle g(x)}
的傅里叶系数
f
^
(
n
)
{\displaystyle {\hat {f}}(n)}
还有
g
^
(
n
)
{\displaystyle {\hat {g}}(n)}
相同,且傅里叶级数都收敛到函数本身,那么可以证明此傅里叶级数具有唯一性,也就是
f
(
x
)
=
g
(
x
)
{\displaystyle f(x)=g(x)}
。换句话说,如果函数
f
(
x
)
{\displaystyle f(x)}
在
[
−
π
,
π
]
{\displaystyle [-\pi ,\pi ]}
上可积,傅里叶系数
f
^
(
n
)
{\displaystyle {\hat {f}}(n)}
为0,对所有的
n
∈
N
{\displaystyle n\in \mathbb {N} }
,那么函数
f
(
x
)
=
0
{\displaystyle f(x)=0}
卷积定理
给定周期为
P
{\displaystyle P}
的函数
s
P
{\displaystyle s_{_{P}}}
和
r
P
{\displaystyle r_{_{P}}}
,它们具有傅里叶级数系数
S
[
n
]
{\displaystyle S[n]}
和
R
[
n
]
{\displaystyle R[n]}
,这里的
n
∈
Z
{\displaystyle n\in \mathbb {Z} }
。
逐点乘积
h
P
(
x
)
≜
s
P
(
x
)
⋅
r
P
(
x
)
{\displaystyle h_{_{P}}(x)\triangleq s_{_{P}}(x)\cdot r_{_{P}}(x)}
,也是周期为
P
{\displaystyle P}
,并且它的傅里叶级数系数是序列
S
{\displaystyle S}
和
R
{\displaystyle R}
的离散卷积 :
H
[
n
]
=
(
S
∗
R
)
[
n
]
{\displaystyle H[n]=(S*R)[n]}
。
周期卷积
h
P
(
x
)
≜
(
s
P
∗
r
)
(
x
)
=
(
s
∗
r
P
)
(
x
)
=
∫
P
s
P
(
τ
)
⋅
r
P
(
x
−
τ
)
d
τ
{\textstyle h_{_{P}}(x)\triangleq (s_{_{P}}*r)(x)=(s*r_{_{P}})(x)=\int _{P}s_{_{P}}(\tau )\cdot r_{_{P}}(x-\tau )\,d\tau }
,也是周期为
P
{\displaystyle P}
,它具有傅里叶级数系数:
H
[
n
]
=
P
⋅
S
[
n
]
⋅
R
[
n
]
{\displaystyle H[n]=P\cdot S[n]\cdot R[n]}
。
在
c
0
(
Z
)
{\displaystyle c_{0}(\mathbb {Z} )}
中的双无限序列
{
c
n
}
n
∈
Z
{\displaystyle \left\{c_{n}\right\}_{n\in Z}}
,是在
L
1
(
[
0
,
2
π
]
)
{\displaystyle L^{1}([0,2\pi ])}
中的傅里叶系数的序列,当且仅当它是在
ℓ
2
(
Z
)
{\displaystyle \ell ^{2}(\mathbb {Z} )}
中的两个序列的卷积 [ 11] 。
微分性质
我们说
f
(
x
)
{\displaystyle f(x)}
属于在
C
k
(
T
)
{\displaystyle C^{k}(\mathbb {T} )}
⇒
{\displaystyle \Rightarrow }
如果
f
(
x
)
{\displaystyle f(x)}
是一个在实数上以
2
π
{\displaystyle 2\pi }
为周期的函数,且
k
{\displaystyle k}
次可微而且
k
{\displaystyle k}
阶连续。
如果
f
(
x
)
{\displaystyle f(x)}
属于在
C
1
(
T
)
{\displaystyle C^{1}(\mathbb {T} )}
,那么
f
′
(
x
)
{\displaystyle f'(x)}
傅里叶系数
f
′
^
(
n
)
{\displaystyle {\hat {f'}}(n)}
可以被用
f
(
x
)
{\displaystyle f(x)}
傅里叶系数
f
^
(
n
)
{\displaystyle {\hat {f}}(n)}
的表示,借由公式
f
′
^
(
n
)
=
i
n
f
^
(
n
)
{\displaystyle {\hat {f'}}(n)=in{\hat {f}}(n)}
如果
f
(
x
)
{\displaystyle f(x)}
属于在
C
k
(
T
)
{\displaystyle C^{k}(\mathbb {T} )}
,
f
(
k
)
^
(
n
)
=
(
i
n
)
k
f
^
(
n
)
{\displaystyle {\hat {f^{(k)}}}(n)=(in)^{k}{\hat {f}}(n)}
。特别的,当固定
k
≥
1
{\displaystyle k\geq 1}
,我们有
f
(
k
)
^
(
n
)
{\displaystyle {\hat {f^{(k)}}}(n)}
趋近于0当
n
→
∞
{\displaystyle n\rightarrow \infty }
,且有
f
(
k
)
^
(
n
)
=
O
(
1
/
n
k
)
{\displaystyle {\hat {f^{(k)}}}(n)=O(1/n^{k})}
。
黎曼-勒贝格定理
如果
S
{\displaystyle S}
是可积函数 ,则
lim
|
n
|
→
∞
S
[
n
]
=
0
{\textstyle \lim _{|n|\to \infty }S[n]=0}
,
lim
n
→
+
∞
a
n
=
0
{\textstyle \lim _{n\to +\infty }a_{n}=0}
而
lim
n
→
+
∞
b
n
=
0
{\textstyle \lim _{n\to +\infty }b_{n}=0}
。
帕塞瓦尔定理
如果函数
f
(
x
)
{\displaystyle f(x)}
属于在
L
2
(
[
−
π
,
π
]
)
{\displaystyle L^{2}([-\pi ,\pi ])}
之中,那么便有
∑
−
∞
∞
|
f
^
(
n
)
|
2
=
1
2
π
∫
−
π
π
|
f
(
x
)
|
2
d
x
=
|
|
f
|
|
{\textstyle \sum _{-\infty }^{\infty }|{\hat {f}}(n)|^{2}={\frac {1}{2\pi }}\int _{-\pi }^{\pi }|f(x)|^{2}dx=||f||}
。
普朗歇尔定理
如果
c
0
,
c
±
1
,
c
±
2
,
…
{\displaystyle c_{0},\,c_{\pm 1},\,c_{\pm 2},\ldots }
是系数,并且
∑
n
=
−
∞
∞
|
c
n
|
2
<
∞
{\textstyle \sum _{n=-\infty }^{\infty }|c_{n}|^{2}<\infty }
,则有一个唯一的函数
s
∈
L
2
(
P
)
{\displaystyle s\in L^{2}(P)}
使得对于所有
n
{\displaystyle n}
有着
S
[
n
]
=
c
n
{\displaystyle S[n]=c_{n}}
。
延伸
希尔伯特空间的解读
正弦和余弦形成了正交集合。正弦、余弦及其乘积的积分,当m 与n 不同或二函数不同时是0(绿色和红色区域相等抵消),仅当m 和n 相等并且函数相同时为π。
所谓的两个不同向量 正交是指它们的内积 为0,这也就意味着这两个向量之间没有任何相关性,例如,在三维欧氏空间中,互相垂直的向量之间是正交的。事实上,正交是垂直在数学上的一种抽象化和一般化。一组n个互相正交的向量必然是线性无关 的,所以必然可以张成一个n维空间,也就是说,空间中的任何一个向量可以用它们来线性表出。
在希尔伯特空间 释义下,函数的集合{en = einx ; n ∈ Z }是[−π, π]平方可积函数L 2 ([−π, π])的正交基 。这个空间实际上是一个希尔伯特空间,有着针对任何两个的元素f 和g 的如下内积:
⟨
f
,
g
⟩
=
d
e
f
1
2
π
∫
−
π
π
f
(
x
)
g
(
x
)
¯
d
x
.
{\displaystyle \langle f,\,g\rangle \;{\stackrel {\mathrm {def} }{=}}\;{\frac {1}{2\pi }}\int _{-\pi }^{\pi }f(x){\overline {g(x)}}\,dx.}
三角函数族的正交性用公式表示出来就是:
∫
−
π
π
cos
(
m
x
)
cos
(
n
x
)
d
x
=
π
δ
m
n
,
m
,
n
≥
1
,
{\displaystyle \int _{-\pi }^{\pi }\cos(mx)\,\cos(nx)\,dx=\pi \delta _{mn},\quad m,n\geq 1,\,}
∫
−
π
π
sin
(
m
x
)
sin
(
n
x
)
d
x
=
π
δ
m
n
,
m
,
n
≥
1
{\displaystyle \int _{-\pi }^{\pi }\sin(mx)\,\sin(nx)\,dx=\pi \delta _{mn},\quad m,n\geq 1}
(这里的δmn 是克罗内克函数 ),而
∫
−
π
π
cos
(
m
x
)
sin
(
n
x
)
d
x
=
0
;
{\displaystyle \int _{-\pi }^{\pi }\cos(mx)\,\sin(nx)\,dx=0;\,}
参阅
注释
^ 一些作者定义了与此不同的
A
0
{\displaystyle A_{0}}
,使得可以用相同的积分定义
A
0
{\displaystyle A_{0}}
和
A
n
{\displaystyle A_{n}}
。这改变了Eq. 2 使得第一项需要除以
2
{\displaystyle 2}
。
^ 因为周期函数的傅里叶变换的积分定义不是收敛的,需要将周期函数和它的变换视为分布 。在这种意义上,
F
{
e
i
2
π
n
x
P
}
{\displaystyle {\mathcal {F}}\{e^{i{\frac {2\pi nx}{P}}}\}}
是一个狄拉克δ函数 ,它是分布的是例子。
引用
^ 详见莫里斯·克莱因 《古今数学思想》,第20章无穷级数,第5节三角级数;第28章十九世纪的偏微分方程,第5节热方程与傅里叶级数。see here, pg.s 209 & 210, (页面存档备份 ,存于互联网档案馆 )
^ 李狗嗨. 如何给文科生解释傅里叶变换? . 知乎专栏. 2019-07-25 [2020-02-07 ] . (原始内容存档 于2020-10-24) (中文) .
^ Théorie analytique de la chaleur . Paris: Firmin Didot Père et Fils. 1822. OCLC 2688081 (法语) .
^ Pinkus, Allan; Zafrany, Samy. Fourier Series and Integral Transforms 1st. Cambridge, UK: Cambridge University Press. 1997: 42 –44. ISBN 0-521-59771-4 .
^ Katznelson, Yitzhak. An introduction to Harmonic Analysis 2nd corrected. New York, NY: Dover Publications, Inc. 1976: 46 . ISBN 0-486-63331-4 .
^ Georgi P. Tolstov. Fourier Series . Courier-Dover. 1976. ISBN 0-486-63317-9 .
^ 这里 (页面存档备份 ,存于互联网档案馆 )
^ 8.0 8.1 8.2 8.3 8.4 Papula, Lothar. Mathematische Formelsammlung: für Ingenieure und Naturwissenschaftler [Mathematical Functions for Engineers and Physicists]. Vieweg+Teubner Verlag. 2009. ISBN 978-3834807571 (德语) .
^ 9.0 9.1 9.2 9.3 Shmaliy, Y.S. Continuous-Time Signals. Springer. 2007. ISBN 978-1402062711 .
^ Proakis, John G.; Manolakis, Dimitris G. Digital Signal Processing: Principles, Algorithms, and Applications 3rd. Prentice Hall. 1996: 291 . ISBN 978-0-13-373762-2 .
^ Characterizations of a linear subspace associated with Fourier series . MathOverflow. 2010-11-19 [2014-08-08 ] .
延伸阅读
William E. Boyce; Richard C. DiPrima. Elementary Differential Equations and Boundary Value Problems 8th. New Jersey: John Wiley & Sons, Inc. 2005. ISBN 0-471-43338-1 .
Joseph Fourier, translated by Alexander Freeman. The Analytical Theory of Heat. Dover Publications. 2003. ISBN 0-486-49531-0 . 2003 unabridged republication of the 1878 English translation by Alexander Freeman of Fourier's work Théorie Analytique de la Chaleur , originally published in 1822.
Enrique A. Gonzalez-Velasco. Connections in Mathematical Analysis: The Case of Fourier Series . American Mathematical Monthly. 1992, 99 (5): 427–441. JSTOR 2325087 . doi:10.2307/2325087 .
Fetter, Alexander L.; Walecka, John Dirk. Theoretical Mechanics of Particles and Continua . Courier. 2003. ISBN 978-0-486-43261-8 .
Felix Klein , Development of mathematics in the 19th century . Mathsci Press Brookline, Mass, 1979. Translated by M. Ackerman from Vorlesungen über die Entwicklung der Mathematik im 19 Jahrhundert , Springer, Berlin, 1928.
Walter Rudin . Principles of mathematical analysis 3rd. New York: McGraw-Hill, Inc. 1976. ISBN 0-07-054235-X .
A. Zygmund . Trigonometric Series third. Cambridge: Cambridge University Press. 2002. ISBN 0-521-89053-5 . The first edition was published in 1935.
电机电子类科《工程数学》,ISBN 978-957-584-377-9 ,作者 陈锡冠、曾致煌,高立出版社。
外部链接