• 设为首页
  • 点击收藏
  • 手机版
    手机扫一扫访问
    迪恩网络手机版
  • 关注官方公众号
    微信扫一扫关注
    迪恩网络公众号

[转]阿英 Matlab fftshift 详解 - 信号处理基本功 - casio1374633

原作者: [db:作者] 来自: [db:来源] 收藏 邀请

http://anony3721.blog.163.com/blog/static/511974201031395110529/

愿消三障诸烦恼,愿得智慧真明了,普愿罪障悉消除,世世常行菩萨道。

愿以此功德。庄严佛净土。上报四重恩。下济三途苦。若有见闻者。悉发菩提心。尽此一报身。同生极乐国。

 

一. 实信号情况

因为实信号以fs为采样速率的信号在 fs/2 处混叠,所以实信号fft的结果中前半部分对应[0, fs/2],后半部分对应[ -fs/2, 0]

1)实信号fft的结果前半部分对应[0, fs/2]是正频率的结果,后半部分对应[ -fs/2, 0]是负频率的结果。大于fs/2的部分的频谱实际上是实信号的负频率加fs的结果。故要得到正确的结果,只需将视在频率减去fs即可得到频谱对应的真实负频率

2)如果要让实信号fft的结果与[-fs/2, fs/2]对应,则要fft后fftshift一下即可,fftshift的操作是将fft结果以fs/2为中心左右互换

3)如果实信号fft的绘图频率f从[-fs/2, fs/2],并且没有fftshift,则fft正频谱对应f在[0, fs/2]的结果将混叠到(f - fs/2)的位置;

fft负频谱对应f在[-fs/2, 0]的结果混叠到 f + fs - fs/2 的位置,注意这里f为负值,也就是说此种情况下fft负频谱对应的视在频率减去fs/2即可得到频谱对应的真实负频率

 

二. 复信号情况

1)复信号没有负频率,以fs为采样速率的信号,fft的频谱结果是从[0, fs]的。

2)在 f > fs/2 时,对复信号的fft结果进行fftshift会产生频率混叠(将下面的示例2中的频率从f=15改为f=85可以验证f=85的谱线在fftshift后跑到 f = -15 = 85 - fs = 85 - 100的位置了),所以复信号也一般要求 f <= fs/2

3)在对雷达的慢时间维(复信号)进行fft后,由于要用doppler = ((0:LFFT-1)/LFFT  - 0.5)*PRF; 计算多普勒频率,所以对该慢时间信号fft后要fftshift下,以便和正确的频率单元相对应。注意多普勒频率fd < = PRF/2 时才测的准!

 

fftshift

作用:将零频点移到频谱的中间

用法:

Y=fftshift(X)
Y=fftshift(X,dim)

描述:fftshift移动零频点到频谱中间,重新排列fft,fft2和fftn的输出结果。将零频点放到频谱的中间对于观察傅立叶变换是有用的。

示例1 - 实信号的情况:

clf;

fs=100;N=256;   %采样频率和数据点数
n=0:N-1;t=n/fs;   %时间序列
x=0.5*sin(2*pi*15*t)+2*sin(2*pi*40*t); %信号

y1=fft(x,N);    %对信号进行快速Fourier变换
y2=fftshift(y1);

mag1=abs(y1);     %求得Fourier变换后的振幅
mag2=abs(y2);    

f1=n*fs/N;    %频率序列
f2=n*fs/N-fs/2;

subplot(3,1,1),plot(f1,mag1,\'r\');   %绘出随频率变化的振幅
xlabel(\'频率/Hz\');
ylabel(\'振幅\');title(\'图1:usual FFT\',\'color\',\'r\');grid on;

subplot(3,1,2),plot(f2,mag1,\'b\');   %绘出随频率变化的振幅
xlabel(\'频率/Hz\');
ylabel(\'振幅\');title(\'图2:FFT without fftshift\',\'color\',\'b\');grid on;

subplot(3,1,3),plot(f2, mag2,\'c\');   %绘出随频率变化的振幅
xlabel(\'频率/Hz\');
ylabel(\'振幅\');title(\'图3:FFT after fftshift\',\'color\',\'c\');grid on;

结论:

1)如果期望绘制的幅频图的频率范围为0~fs,则无需运行fftshift变换,正频率对应在[0, fs/2],

大于fs/2的频谱的频率值为对应[-fs/2  , 0 ]负频率f + fs,注意f是负频率,是个负数。如图1。

2)如果期望绘制的幅频图的频率范围为-fs/2~fs/2,则需要运行fftshift变换,如图3;

如果不变换,图示的响应频点会发生变换,如图2,分析见顶端。

示例2 - 复信号的情况:

close all; clear; clf;

fs=100;N=256;   %采样频率和数据点数
n=0:N-1;t=n/fs;   %时间序列
x=0.5*exp(j*2*pi*15*t)+2*exp(j*2*pi*40*t); %信号

y1=fft(x,N);    %对信号进行快速Fourier变换
y2=fftshift(y1);

mag1=abs(y1);     %求得Fourier变换后的振幅
mag2=abs(y2);    

f1=n*fs/N;    %频率序列
f2=n*fs/N-fs/2;

subplot(3,1,1),plot(f1,mag1,\'r\');   %绘出随频率变化的振幅
xlabel(\'频率/Hz\');
ylabel(\'振幅\');title(\'图1:usual FFT\',\'color\',\'r\');grid on;

subplot(3,1,2),plot(f2,mag1,\'b\');   %绘出随频率变化的振幅
xlabel(\'频率/Hz\');
ylabel(\'振幅\');title(\'图2:FFT without fftshift\',\'color\',\'b\');grid on;

subplot(3,1,3),plot(f2, mag2,\'c\');   %绘出随频率变化的振幅
xlabel(\'频率/Hz\');
ylabel(\'振幅\');title(\'图3:FFT after fftshift\',\'color\',\'c\');grid on;

--
Anony

专注的力量成就梦想
posted on 2013-10-29 21:56  casio1374633  阅读(943)  评论(0编辑  收藏  举报

鲜花

握手

雷人

路过

鸡蛋
该文章已有0人参与评论

请发表评论

全部评论

专题导读
热门推荐
阅读排行榜

扫描微信二维码

查看手机版网站

随时了解更新最新资讯

139-2527-9053

在线客服(服务时间 9:00~18:00)

在线QQ客服
地址:深圳市南山区西丽大学城创智工业园
电邮:jeky_zhao#qq.com
移动电话:139-2527-9053

Powered by 互联科技 X3.4© 2001-2213 极客世界.|Sitemap