快速傅立叶变换(Matlab实现)-林志立.pptVIP

快速傅立叶变换(Matlab实现) 林志立 华侨大学信息学院 zllin@ Lord Kelvin on Fourier’s theorem Fourier’s theorem is not only one of the most beautiful results of modern analysis, but it may be said to furnish an indispensable instrument in the treatment of nearly every recondite question in modern physics. Lord Kelvin (1824-1907) Joseph Fourier（1768 – 1830） Fourier was obsessed with the physics of heat and developed the Fourier series and transform to model heat-flow problems. 今天的主角 为什么要进行傅立叶变换? We desire a measure of the frequencies present in a wave. This will lead to a definition of the term, the spectrum. Plane waves have only one frequency, w. This light wave has many frequencies. And the frequency increases in time (from red to blue). It will be nice if our measure also tells us when each frequency occurs. Light electric field Time 傅立叶变换－针对对非周期函数 复数的F(w)的物理意义到底是什么？ 若f(t)为实数，一般情况下F(w)是复数。 这公式什么来的？ 从傅立叶级数说起 假设有一个以 T为周期的函数 ，满足狄利克雷（Dirichlet）条件, 其中： 根据欧拉（Euler）公式，可以得到傅立叶级数的复数形式 Euler公式 从傅立叶级数说起 可以写成单纯的正弦函数形式 所以周期函数可以表示为一系列正弦函数之和。 从傅立叶级数说起 也可以写成单纯的余弦函数形式 所以周期函数可以表示为一系列余弦函数之和。 傅立叶变换 频率间隔： T趋于无穷大时，周期函数变为非周期。 则有： 傅立叶变换及其存在性 傅立叶积分公式 傅立叶变换对 实际上， 其中， 为无数频率连续的正弦函数之和。 离散傅立叶变换 g(t) 连续函数： t 离散化： 即使采样点数不多，不过近似的精度很高。 离散傅立叶变换(采样) 离散傅立叶变换(算法) 1965年,第一篇关于fft算法论文 离散傅立叶变换 根据计算机语言系统的特点，为了能进行快速傅立叶变换，需要对数据进行重新编排，使得n和m取值为： (1)取非负整数（C语言）(0,1,2,3,…,N-1) (2)正整数（Matlab语言) (1,2,3,…,N) 对于偶数N, 对于奇数N, Matlab中的FFT函数 Matlab Labview N为偶数 Matlab Labview N为奇数 注意实际相位问题 t t=0 实际时间轴 fft时间轴 t0 则： 注意实际相位问题 t t=0 实际时间轴 fft时间轴 t0 则： 注意要从g(0)开始采样。t=0 fftshift与ifftshift N为偶数时，它们没有区别。 a=[1,2,3,4]; fftshift(a) ifftshift(a) ans = 3 4 1 2 ans = 3 4 1 2 N为奇数时，它们有区别。 a=[1,2,3,4,5]; fftshift(a) ifftshift(a) ans = 4 5 1 2 3 ans = 3 4 5 1 2 Matlab实现离散傅立叶变换的指令 N为偶数时， gl=fftshift(g); Gl=delta*fft(gl); G=fftshift(Gl); G=fftshift(fft(fftshift(g)))*delta; N为奇数时， gl=ifftshift(g); Gl=delta*fft(gl); G=fftshift(Gl); G=fft