傅里叶变换

                     

贡献者: addis

预备知识 1 傅里叶级数(指数),傅里叶变换(三角)

   用三角傅里叶变换中同样的方法可把指数傅里叶级数的区间长度 l 取极限后拓展为指数傅里叶变换

(1)g(k)=12π+f(x)eikxdx ,
(2)f(x)=12π+g(k)eikxdk .
f(x) 为实函数时,g(k) 的实部是偶函数,虚部是奇函数。

定理 1 实函数、奇函数、偶函数的傅里叶变换

  

  1. f(x) 是实函数,当且仅当 g(k)=g(k)
  2. g(x) 是实函数,当且仅当 f(x)=f(x)
  3. f(x) 是奇(偶)函数,当且仅当 g(k) 是奇(偶)函数。
  4. f(x) 同时是实函数和偶函数,当且仅当 g(k) 同时是实函数和偶函数。
  5. f(x) 同时是实函数和奇函数,当且仅当 g(k) 同时是虚函数和奇函数。

   证明第 1 条:

   f(x) 为实函数的充要条件是 f(x)=f(x),代入式 1

(3)g(k)=12π+f(x)eikxdx .
两边取复共轭,得
(4)g(k)=g(k) .
注意 g(k) 的实部是偶函数,虚部是奇函数,因此往往只需要 k 的正半轴即可得到所有信息。上式就是 f(x) 为实函数的充要条件,要证明充分性,将式 4 其代入式 2 可得 f(x)=f(x)。证毕。(其余证明留做习题:提示,偶函数的傅里叶变换相当于 cos 变换)

例 1 高斯分布的傅里叶变换

   要计算高斯函数

(5)f(x)=eax2(a>0) 
的傅里叶变换,代入式 1 并使用例 2
(6)g(k)=12aexp(k24a) ,
一个方便的记忆法是 x2 前的系数乘以 k2 前的系数相乘等于 1/4

例 2 

  

(7)f(x)={exp(ik0x)cos2(ax)(|x|<π2a) .0(其他)
则傅里叶变换为
(8)g(k)=2πa4a2(kk0)2sinc[π(kk0)2a] ,
其中 sinc 函数见相关页面

例 3 方波

   区间 [l,l],高为 1 的单个方波,

(9)g(k)=2πsin(lk)k .

1. 证明

预备知识 2 狄拉克 delta 函数

  

未完成:以下证明可能存在问题,需要专业人士审核
以下的证明可以用矢量空间和基底的概念得到更深刻的理解,详见 “傅里叶变换与连续正交归一基底”。

   我们把式 2 看作定义,用狄拉克 δ 函数来证明式 1 ,反之同理。把式 2 代入式 1

(10)g(k)=12π+[+g(k)eikxdk]eikxdx .
这就是我们需要证明的。我们把无穷的积分上下限改写为极限,即
(11)g(k)=12πlimnnn[limmmmg(k)eikxdk]eikxdx .
如果内极限可以移动到两个积分外,有
未完成:极限和积分交换的条件是什么?是否满足?
(12)g(k)=12πlimnlimmnnmmg(k)ei(kk)xdkdx .
我们假设 g(k)[m,m] 内绝对值可积,那么有限区间的重积分可以交换顺序,变为
(13)12πlimnlimmmmg(k)nnei(kk)xdxdk=limnlimmmmg(k)δn(kk)dk=limng(k)δn(kk)dk=g(k) .
第二个等号中,δn例 2 中的 delta 函数列,最后一步使用了其性质式 9 。证毕。

   注意这里要求对每个 n=1,2, 积分 g(k)δn(kk)dk 都收敛。只有满足该要求的函数才适用该证明。

2. 性质

   为了书写方便我们用算符 FF1 表示傅里叶变换和反变换,即 Ff=g 以及 F1g=f。算符在这里可以看作 “函数的函数”,即自变量和函数值都是函数。

   平移

(14)F[f(x)eik0x]=g(kk0) .
(15)F[f(xx0)]=g(k)eikx0 .
也就是说,给函数乘以 eik0x 因子再做傅里叶变换,等于先对函数做傅里叶变换,再向右平移 k0;给函数再向右平移 x0 再做反傅里叶变换,等于先对函数做傅里叶变换,再乘以 eix0k。证明留做习题。

   模长不变性

(16)+g(k)g(k)dk=+f(x)f(x)dx .

   拉伸

(17)F[f(ax)]=1ag(ka) .
也就是说把函数在 x 方向压缩 a 倍后,各个频率都变大 a 倍,所以傅里叶变换会在 k 方向拉伸 a 倍,另外归一化不变性易得系数 1/a

   导数

(18)F[f(x)]=ikg(k) ,
同理
(19)F1[g(k)]=ixf(x) .

   作为式 16 的拓展,有

(20)+f1(x)f2(x)dx=+g1(k)g2(k)dk .
这可以理解为傅里叶变换不改变内积,所以是一个无穷维空间中的幺正变换

   如果 f1(x) 可以在 x=0 泰勒展开,有

(21)F[f1(x)f2(x)]=f1(ik)g2(k) .
如果 g1(k) 可以在 k=0 泰勒展开,有
(22)F1[g1(k)g2(k)]=g1(ix)f2(x) .
注意式 18 式 19 是该性质的特殊情况(令式 21 f1(x)=x式 22 g1(k)=k)。记忆方法:在傅里叶变换外面的 i/k 相当于傅里叶变换里面的 x,反傅里叶变换外面的 i/x 相当于反傅里叶变换里面的 k

   平均值

(23)k=+k|g(k)|2dx=+f(x)(iddx)f(x)dx .
(24)x=+x|f(x)|2dx=+g(k)(iddk)g(k)dk .
推导参考 “平均值(量子力学)”。

   不确定性原理

(25)σxσk12 .
其中 σx,σk 分别是 f(x),g(k) 的标准差。推导参考量子力学的 “不确定性原理”。

3. 性质的证明

   证明式 16 :把傅里叶变换看成傅里叶级数在 l 时的极限,使用式 14 ,右边的求和在极限下变为积分即可证明。详细过程留做习题。

   证明式 18 式 19 同理):对式 2 关于 x 求导得

(26)f(x)=12π+[ikg(k)]eikxdk .
把方括号看作一整个 k 的函数,那么上式对应的反变换为
(27)ikg(k)=12π+f(x)eikxdx=F[f(x)] .
其中 g(k)=F[f(x)],证毕。

   证明式 21 式 22 同理):过程和式 18 类似,等式右边为

(28)12πf1(ik)+f2(x)eikxdx=12π+f2(x)f1(ik)eikxdx=12π+f1(x)f2(x)eikxdx .

                     

© 小时科技 保留一切权利