【Matlab】快速傅⾥叶变换FFTfftshiftfftshift(fft(fftshi。。。
【⾃我理解】
fft:可以指定点数的快速傅⾥叶变换
fftshift:将零频点移到频谱的中间⽤法:
Y=fftshift(X)
Y=fftshift(X,dim)
描述:fftshift移动零频点到频谱中间,重新排列fft,fft2和fftn的输出结果。将零频点放到频谱的中间对于观察傅⽴叶变换是有⽤的。fftshift(fft(fftshift(x)))
先将s搬到中⼼,然后fft变换,再将变换后的移到中⼼。使⽤fftshift(fft(fftshift(x)))后的效果:1.不改变频谱的幅度和相位
2.使得频谱的范围为-Fs/2到Fs/2,中⼼频率为0准确的频率是:
int n=N%2==0?N/2:((N-1)/2);最左边频率 -Fs/N*n最右边频率 Fs/N*(N-n-1)其中N是系列长度,Fs是采样频率3.频谱的实部全为正数,震荡性消失使⽤注意:
为使得输出结果与理论分析相同,需要乘以系数FFT: fftshift(fft(fftshift(x))) * dtIFFT: fftshift(ifft(fftshift(X))) * fs其次:
MATLAB中实现的IFFT已经除以序列长度L,FFTW中没有除以序列长度,需要⾃⼰补上