为什么使用 fftshift(fft(fftshift(x))) 而不是 fft(x)

详细请查阅《Why use fftshift(fft(fftshift(x))) instead of fft(x) in Matlab?》

使用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))) * dt

IFFT: fftshift(ifft(fftshift(X))) * fs


其次:

MATLAB中实现的IFFT已经除以序列长度L,FFTW中没有除以序列长度,需要自己补上

  • 10
    点赞
  • 33
    收藏
    觉得还不错? 一键收藏
  • 2
    评论
评论 2
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值