作者realtemper (Lunatic Pandora)
看板NTUEE112HW
标题[闲聊] 傅立叶转换的基础概念
时间Mon Dec 15 08:32:54 2008
如果把函数想像成向量
sin(nπx/L) 跟 cos(nπx/L) 想像成基底向量(即 e1, e2 一类的东西)
那麽一言以蔽之,傅立叶转换就是在计算函数的某个分量。
举个直角坐标的例子
e1=(1,1,0)/√2 , e2=(-1,1,0)/√2, e3=(0,0,1),三者两两正交
若我们要计算向量 v 的 v1 分量,就是用 v‧e1,相信大家都知道的。
为什麽这件事情是对的呢?
那是因为在正交坐标系里头,"▁▁‧ei" 这个动作(算子)有一个特性
就是 ej‧ei = δij,其他成分跟自己的内积皆为零,
所以假如我们要取出 v = v1e1 + v2e2 + v3e3 的 v1 成分
那麽只要把等式两边 apply "▁▁‧e1" 算子
就变成
v‧ei = v1e1‧e1 + v2e2‧e1 + v3e3‧e1
^^^^^^^^只有这项留下
那傅立叶转换跟这有什麽关系?
为什麽 sin(nπx/L) 跟 cos(nπx/L) 可以当作基底向量,
并以线性组合的方式来凑出函数f(x)呢?
这是一个数学上的事实(这边只讲大意,精确的版本要参考高微):
1.有限闭区间(长度L)上的连续函数
都可以表示成 sin(nπx/L) 跟 cos(nπx/L) 的线性组合
2.无限区间上的连续函数,
都可以表示成 sinkx 与 coskx 的线性组合(或者等价的,exp(ikx))
但是从物理的角度来看,这两件事情是很显然的:
1.在闭区间上,若把函数f(x)想像成弦的振动
那麽你们都知道,弦上只有一些稳定的振动模式能够存活
就是驻波,sin(nπx/L) 跟 cos(nπx/L)
所以任何稳定的驻波振动,都必然是这两种函数的线性组合,
不可能包含其他的成分。
故以这些驻波模式做为基底,是一种很自然的看法。
2.如果是无限区间,那麽弦振动的方式就可以被随意决定,
所以会包含任意频率的振动模式
(也有一种看法是 L→∞ 时,sin(nπx/L)就变成连续的了)
所以「如何以 sin(nπx/L) 跟 cos(nπx/L) 的线性组合凑出连续函数 f(x)」
这个凑法的存在性有数学定理保证
所以 f(x) = Σan sin(nπx/L) + Σbn cos(nπx/L) +b0/2 这种写法是存在的,
问题就是如何找出 an 跟 bn(即各成分的振幅)。
怎麽找呢?
注意到
L/2 L/2
转换算子 ∫▁▁sin(nπx/L) dx 或 ∫▁▁cos(nπx/L) dx
-L/2 -L/2
其实具有内积的特性:
(1)∫sin(mπx/L)cos(nπx/L) dx = 0
(2)∫sin(mπx/L)sin(nπx/L) dx = L/2 δmn
(3)∫cos(mπx/L)cos(nπx/L) dx = L/2 δmn
而这点也可以用物理的角度来看:
若把∫f(x)g(x)dx 解释成相关系数,fg同号的区间贡献正值,反之则贡献负值,
那麽两个独立的驻波振动模式,同正或同负的区间长度将会均等。
画个第1泛音跟第2泛音,观察其正负便可明白。
如果把 L/2 除掉(normalize),
那麽 (2/L)∫▁▁sin(nπx/L) dx 跟 (2/L)∫▁▁cos(nπx/L) dx这种动作
就相当於直角坐标下的"▁▁‧ei"算子
因此,我们若要计算 下面等式的各项系数
f(x) = Σan sin(nπx/L) + Σbn cos(nπx/L) +b0/2
只要两边 apply (2/L)∫▁▁sin(nπx/L) dx 或 (2/L)∫▁▁cos(nπx/L) dx 即可
所以才会得到下面公式
an = (2/L)∫f(x)sin(nπx/L) dx
bn = (2/L)∫f(x)cos(nπx/L) dx (包含b0)
无限区间的内积牵涉到 delta-function
不过除此之外,概念基本上是类似的。
--
※ 发信站: 批踢踢实业坊(ptt.cc)
◆ From: 203.67.40.87
※ 编辑: realtemper 来自: 203.67.40.87 (12/15 09:37)
1F:推 yunmeow:唔唔......怀念>O< 学弟妹要好好学习喔XDDD 12/15 09:51
2F:推 cookiecake31:学姐好呛 12/15 23:50
3F:推 kickpp:懂了 推助教认真~~! 12/16 00:54
4F:推 cjshine:有读有推 12/16 17:53
5F:推 jl3000x:有转有推 12/17 01:18
6F:推 cokewolf:推认真助教 12/17 20:39
7F:→ cokewolf:不过这算闲聊等级的吗 XDDDDD 12/17 20:41
8F:推 herbert570:Fourier 的影子大二科科都会看到啊... 12/19 12:18
9F:→ realtemper:当然是闲聊啊 我只是讲基本的想法 根本没算任何东西 12/19 21:01
10F:→ realtemper:对了,老师教bra-ket符号了没? 12/19 21:01
11F:→ realtemper:如果提到的话跟我讲一声 我会在助教时间好好解释 12/19 21:03
12F:→ realtemper:这套符号(包括"表象vs物自身",还有dual-space的概念) 12/19 21:04
13F:→ cjshine:还没 12/19 22:01