声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

查看: 1273|回复: 2

[FFT] 关于FFT频谱的细化谱,请高手帮忙

[复制链接]
发表于 2009-5-31 17:59 | 显示全部楼层 |阅读模式

马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。

您需要 登录 才可以下载或查看,没有账号?我要加入

x
本帖最后由 wdhd 于 2016-9-18 13:59 编辑

  N=5120;
  f=5120;
  n=0:(N-1);
  t=n/f;
  x=sin(2*pi*50*t);
  X=fft(x);
  f=f*(0:N-1)/N;
  % 被修改
  subplot(2,1,1);
  plot(f,abs(X)*2/N,'-');
  % 被修改
  axis([0,115,0,1.2]);
  % 被修改
  set(gca, 'XTickMode', 'manual', 'XTick', [0, 20,30,40,45,50,55,60,70,80,90,100,110,115]);
  xlabel('Frequency');
  ylabel('|F(k)|');
  grid on;
  f0=45;%移动的频率
  h=exp(-j*2*pi*f0*t);
  y=x.*h;%数字变频,把w=100点移动到原点
  y1=resample(y,1,10);%重新采样,采样频率为fs/N
  k1=0:1:511;
  h1=exp((j*2*pi*f0*k1)/512);
  y1=y1.*h1;%反移频,使得前后频率一致
  y1=abs(fft(y1));%快速傅里叶变换
  subplot(2,1,2);
  plot(k1,y1*2/512);
  % 被修改
  axis([45,55,0,1.2]);
  % 被修改
  set(gca, 'XTickMode', 'manual', 'XTick', [45,46,47,48,49,50,51,52,53,54,55]);
  title('细化10倍后的频率特性');
  xlabel('频率');
  ylabel('幅度值');
  grid on;
  请解释一下红色部分的程序。我看不懂。谢谢了
回复
分享到:

使用道具 举报

发表于 2009-6-1 08:36 | 显示全部楼层
本帖最后由 wdhd 于 2016-9-18 13:59 编辑
原帖由 wxlwxl08 于 2009-5-31 17:59 发表
N=5120;
f=5120;
n=0:(N-1);
t=n/f;
x=sin(2*pi*50*t);
X=fft(x);
f=f*(0:N-1)/N;
% 被修改
subplot(2,1,1);
plot(f,abs(X)*2/N,'-');
% 被修改
axis([0,115,0,1.2]);
% 被修改
set(gca, 'XTickMode',  ...

LZ的红色部分主要对细化的原理不了解,建议LZ看一下王济和胡晓编的“ MATLAB在振动信号处理中的应用”(中国水利水电出版社)中有关ZOOM-FFT的有关部分,或高怀纲 “一种分析频谱局部特性的快速算法” 火控雷达技术 1999年9月
 楼主| 发表于 2009-6-1 12:04 | 显示全部楼层

回复 沙发 songzy41 的帖子

额,我去图书馆找找看。谢谢了,但我不一定能找到。大侠,你能不能给我解释一下呀?
您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

QQ|小黑屋|Archiver|手机版|联系我们|声振论坛

GMT+8, 2024-11-11 17:08 , Processed in 0.068480 second(s), 18 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

快速回复 返回顶部 返回列表