最近突然有想法在搞完手头的事情后,做些基于傅里叶的应用研究,加之最近上这方面的课而且前不久看到一篇好文,感觉可以先写篇博客理理自己对于傅里叶的理解。(有错望揪)
正交基
先我们来看两个向量内机(inner product):
内机近似表达的含义是两个向量的相关程度(还没有归一化)。两个完全不相关的向量我们则称其为正交向量。
正交向量一般作为基使用。:
两个相互正交向量可以组成一个平面。三个可以组成三维空间,等等。
函数其实是一个无限维的向量:
定义一个有界函数
对
于
a
<
x
<
b
,
可
表
示
为
[
f
(
a
)
,
f
(
a
+
Δ
x
)
,
f
(
a
+
2
Δ
x
)
,
f
(
a
+
3
Δ
x
)
,
.
.
.
,
f
(
b
)
]
,
i
f
Δ
x
→
0
对于a<x<b,可表示为[f(a),f(a+ \Delta x),f(a+2\Delta x),f(a+3\Delta x),...,f(b)], if\ \Delta x \rightarrow 0
对于a<x<b,可表示为[f(a),f(a+Δx),f(a+2Δx),f(a+3Δx),...,f(b)],if Δx→0
函数的内积可以表示为:
我们再来看傅里叶的那组基
h
p
(
x
)
=
e
i
2
π
p
x
/
T
h_p(x)=e^{i2\pi p x/T}
hp(x)=ei2πpx/T
对于每个可以得到(复数内积其中之一要取共轭;三角函数正交性需要在一个公共周期内取得)
即我们找到无穷个相互正交的函数,由此得到的就是无穷维空间。
我们常见的傅里叶变换基:
h
w
1
(
x
)
=
e
i
w
1
x
h_{w1}(x)=e^{iw_1 x}
hw1(x)=eiw1x
我们也能找到一个公共大周期使得周期函数正交:
T
=
2
Π
/
w
1
∗
2
Π
/
w
2
∗
2
Π
/
w
3
∗
.
.
.
T=2\Pi/w_1*2\Pi/w_2*2\Pi/w_3*...
T=2Π/w1∗2Π/w2∗2Π/w3∗...
这里不管这个T了。
上面我们说的概括来就是找到一组拥有无穷个互相正交的基,又因为函数本身就是一个无穷维的向量,那么我们就能用那组无穷个的正交基来表示任意的无穷维向量(任意函数)
频谱
傅里叶反变换:
f
(
t
)
=
∫
F
(
i
w
)
e
i
w
t
d
w
f(t)=\int F(iw)e^{iwt} dw
f(t)=∫F(iw)eiwtdw
可以把
e
i
w
t
e^{iwt}
eiwt
看作基,那么频谱
就可以看作正交基的权值图,傅里叶变换就是描述这种关系。
无误差拟合
要想无误差拟合,条件是频谱是带限的,采样够密:
但是一般带限的频谱对应无界的时谱,带限的时谱对应无界的频谱(傅里叶变换公式可以发现)。所以因为无法无穷采样,一般无误差拟合也不现实。
(拟合其实是很多内插函数的加和,而一个带限的频谱对应的是一个无界的时谱,则需要无穷多的内插函数相加和才能无误差的拟合。这样显然不现实!)
内插函数部分见数字信号书本。
emmm,知道越写越水,懒了,emmm,其实个人印象傅里叶最牛逼的是那组数量无穷的正交基。
参考:http://songcy.net/posts/story-of-basis-and-kernel-part-1/