Numerical Methods and Modeling in Science数值方法建模这门课最近进入到了一个 我完全不会的区域,傅立叶变换。
傅立叶变换的历史就是让人胆颤的历史。学渣如我,闻之色变。
最近的三周课应该都是讲傅立叶变换,写一遍博客也许就记住了----美好的愿望。
引例
Orthonormal标准正交
Orthonormal=Orthogonal+Unit length
在三维空间中,我们有x,y,z满足
x⋅x=y⋅y=z⋅z=0
x⋅y=y⋅z=z⋅x=1
那么任意向量V=ax+by+cz
在区间[−π,π]中函数g和h的点乘可以表示为
⟨g,h⟩=∫−ππg(x)h(x)dx
一些函数
2π1
π1cosx, π1cos2x ,π1cos3x
π1sinx ,π1sin2x ,π1sin3x
可以被证明是orthonormal的
⟨2π1,2π1⟩=1
⟨π1cos2x,π1cos2x⟩=1
⟨π1sin2x,π1sin2x⟩=1
⟨2π1,π1cos2x⟩=0
…
那么类似于V
函数f(x)=c02π1+c1π1cosx+c2π1cos2x+...+d1π1sinx+d2π1sin2x+...
我们可以利用上面的orthonormal关系计算得到系数
c0=⟨f(x),2π1⟩=∫−ππ2π1f(x)dx
cn=⟨f(x),π1cosnx⟩=∫−πππ1f(x)cosnxdx
dn=⟨f(x),π1sinnx⟩=∫−πππ1f(x)sinnxdx
因此
f(x)=(∫−πππ1f(x)cosnxdx)2π1+∑n=1∞(∫−πππ1f(x)cosnxdx)π1cosnx+∑n=1∞(∫−πππ1f(x)sinnxdx)π1sinnx
化简
f(x)=(π1∫−ππf(x)cosnxdx)21+∑n=1∞(π1∫−ππf(x)cosnxdx)cosnx+∑n=1∞(π1∫−ππf(x)sinnxdx)sinnx
此处有狄利克雷定理,但是我看了好几个版本,看了好久也没太懂,但是仿佛也不重要

然后就是一个小例子

下面的演示是随着使用的term的增多我们会看到多个三角函数叠加在一起的效果,趋近于原函数
1个

2个

3个

12个

个人认为,学习傅立叶变换的思路就是我要用三角函数去模拟原来的函数。为什么要用三角函数呢?我猜是为了后面转换为复变函数。
目前这一部分内容还不是很懂。
留下王一教授给我的朋友圈评论作为激励。
“量子力学告诉我们,位置空间越不apart,Fourier空间越apart。”