声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

查看: 1546|回复: 7

[综合讨论] 分辨率加倍,工频漂移

[复制链接]
发表于 2012-7-25 16:36 | 显示全部楼层 |阅读模式

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

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

x
  1. clc;
  2. clear;
  3. wave=load('30002.txt');
  4. N=length(wave);
  5. for i=1:N;
  6. t(i)=wave(i,2);
  7. end
  8. t=detrend(t);
  9. figure;
  10. subplot(211)
  11. plot(t)
  12. title('原始信号');

  13. Fs=10000

  14. Y=fft(t)
  15. Ayy = (abs(Y)); %取模

  16. Ayy=Ayy/(N/2);   %换算成实际的幅度
  17. Ayy(1)=Ayy(1)/2;
  18. F=([1:N]-1)*Fs/N; %换算成实际的频率值
  19. subplot(212)
  20. plot(F(1:N/2),Ayy(1:N/2));   %显示换算后的FFT模值结果
  21. title('幅度-频率曲线图');
复制代码

以上是程序,数据在附件,结果如上图,现在想提高分辨率,分辨率=采样频率/采样点数,减小采样频率从而提高分辨率,这时结果如下图,问题:工频应该是50Hz的,怎么提高分辨率以后,就漂移了呢?
未命名.JPG

30002.txt

93.67 KB, 下载次数: 4

回复
分享到:

使用道具 举报

 楼主| 发表于 2012-7-26 08:54 | 显示全部楼层
频谱图的横坐标变了,原来50Hz的振幅移到了25Hz。为什么?
发表于 2012-7-26 09:39 | 显示全部楼层
http://forum.vibunion.com/thread-23614-1-1.html
看看这个帖子是否有所帮助

点评

赞成: 5.0
赞成: 5
  发表于 2012-7-26 17:03
发表于 2012-7-26 23:55 | 显示全部楼层
1.相同类似帖没必要发新帖
2.减小采样频率并非单纯将Fs降低而已!
3.for i=1:N; t(i)=wave(i,2); end
  这种写法不好, 直接t=wave(:,2);即可
  1. clc; clear;
  2. wave=load('30002.txt');
  3. t=wave(:,2); t=detrend(t); N=length(t);
  4. Fs=10000; Ayy=abs(fft(t))*2/N; Ayy(1)=Ayy(1)/2;
  5. F=([1:N]-1)*Fs/N; %换算成实际的频率值
  6. figure; subplot(211); plot(t); title('原始信号');
  7. subplot(212); plot(F(1:N/2),Ayy(1:N/2)); title('幅度-频率曲线图');  %显示换算后的FFT模值结果

  8. t=t(1:2:end); t=detrend(t); N=length(t);
  9. Fs=5000; Ayy=abs(fft(t))*2/N; Ayy(1)=Ayy(1)/2;
  10. F=([1:N]-1)*Fs/N; %换算成实际的频率值
  11. figure; subplot(211); plot(t); title('原始信号');
  12. subplot(212); plot(F(1:N/2),Ayy(1:N/2)); title('幅度-频率曲线图');  %显示换算后的FFT模值结果
复制代码
 楼主| 发表于 2012-7-27 09:05 | 显示全部楼层
回复 4 # ChaChing 的帖子

你这两段程序,分辨率分别是多少?我的目的是想提高分辨率,并且横坐标随信号中真实频率变化
 楼主| 发表于 2012-7-27 09:21 | 显示全部楼层
回复 4 # ChaChing 的帖子

第一段:采样点数是N=4620点,Fs是10KHz,分辨率=Fs/N=2.164Hz
第二段:采样隔一个采一次,采样点数是=2310点,Fs是5000Hz,分辨率=2.164Hz
没有达到提高分辨率的目的。我是想问为什么分辨率提高以后(因为分辨率=采样频率/采样点数,所以减小采样频率从而提高分辨率),原有的工频50Hz,却到了25Hz处?
发表于 2012-7-27 14:41 | 显示全部楼层
LZ所提供的“30002.txt”中的时间信号序列,对应采样频率为10kHZ已是无法改变的事实,除非按照ChaChing的方法对该序列进行重采样,才可相应变动采样频率。
因此,要想提高频谱图的频率分辨率,只能人为提高FFT的线数,然后计算出需参与FFT运算的序列长度,即FFT点数。

评分

1

查看全部评分

发表于 2012-8-1 10:21 | 显示全部楼层
额。。。。。好深奥。。。。
您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

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

GMT+8, 2024-11-11 07:01 , Processed in 0.062597 second(s), 22 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

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