matlab快速傅里叶变换讲义.pptVIP

资料仅可参考，如有不妥，请联系本人改正或者删除。 * 2020年08月13日 matlab快速傅里叶变换 资料仅可参考，如有不妥，请联系本人改正或者删除。 * 2020年08月13日 matlab快速傅里叶变换 资料仅可参考，如有不妥，请联系本人改正或者删除。 * 2020年08月13日 matlab快速傅里叶变换 资料仅可参考，如有不妥，请联系本人改正或者删除。 * 2020年08月13日 matlab快速傅里叶变换 资料仅可参考，如有不妥，请联系本人改正或者删除。 * 2020年08月13日 matlab快速傅里叶变换 资料仅可参考，如有不妥，请联系本人改正或者删除。 * 2020年08月13日 matlab快速傅里叶变换 资料仅可参考，如有不妥，请联系本人改正或者删除。 * 2020年08月13日 matlab快速傅里叶变换 资料仅可参考，如有不妥，请联系本人改正或者删除。 * 2020年08月13日 matlab快速傅里叶变换 资料仅可参考，如有不妥，请联系本人改正或者删除。 * 2020年08月13日 matlab快速傅里叶变换 * matlab快速傅里叶变换 * 说明 * matlab快速傅里叶变换 * §4.5.2 傅里叶级数与傅里叶积分（阅读） 书上这些内容仅供有兴趣的读者参考。 傅里叶级数：将变量范围为无限的连续的周期函数变换成无限的离散的傅里叶频谱序列。 傅里叶积分：对变量范围为无限的连续的非周期函数变换成无限的连续的傅里叶频谱序列。 离散时间傅里叶变换：将无限长数据序列变换为周期函数。 * matlab快速傅里叶变换 * * matlab快速傅里叶变换 * 快速傅里叶变换(FFT)是离散傅里叶变换(DFT)的快速算法。 §4.5.3 快速傅里叶变换的指令 * matlab快速傅里叶变换 * 指令及语句格式: ??????(??)；??????(??,??) 1.一维的FFT及其逆变换 * matlab快速傅里叶变换 * 指令及其语句格式：????????(??) 2. 二维的FFT及其逆变换 * matlab快速傅里叶变换 * X = [ 4, 0, 3, 6, 2, 9, 6, 5] ; Y=fft(X) Y = 35 %零频分量 -5.0711 +8.6569i -3 + 2i 9.0711 +2.6569i -5 %尼奎斯特频率分量 9.0711 - 2.6569i -3 - 2i -5.0711 - 8.6569i 例 共轭分量 一维傅里叶变换 * matlab快速傅里叶变换 * 取样频率1000Hz, 取样间隔1/1000s=1ms 取样数250 总取样时间0.25S 建立时间轴：从0到0.25s 构造正弦波的信号 包含频率50Hz与100Hz t = 0:0.001:.25; x = sin(2*pi*50*t) + sin(2*pi*120*t); y = x + 2*randn(size(t)); %噪声信号 plot(y(1:50)) title(Noisy time domain signal) 例 * matlab快速傅里叶变换 * Y = fft(y,256); Pyy = Y.*conj(Y)/256; f = 1000/256*(0:127); plot(f,Pyy(1:128)) title(Power spectral density) xlabel(Frequency (Hz)) * matlab快速傅里叶变换 * 例 t = 0:1/100:10-1/100; x = sin(2*pi*15*t) +sin(2*pi*40*t); y = fft(x); m = abs(y); f =(0:numel(y)-1)*100/numel(y); plot(f,m) grid on ylabel(Abs. Magnitude) xlabel(Frequency [Hertz]) 指令numel：计算列阵中的元素总数。 * matlab快速傅里叶变换 * 通过分析300年左右的太阳黑子活动的记录资料可以发现每过11年就会达到一个极大值。 Wolfer数图： 太阳黑子的数目与面积 例 放大的前50年的数据 分析太阳黑子活动的周期性 * matlab快速傅里叶变换 * 图片取自维基百科 * matlab快速傅里叶变换 * 傅里叶系数Y在复平面上的分布 Y = fft(wolfer); Y(1)=[]; p