[原创]zoomfft

miao7mijao · 发表于 2006-10-16 13:40

原帖由 yangzj 于 2006-9-23 08:11 发表

输入序列x要足够长，至少要D*N+2*M点数据.

错误提示为：Attempted to access x(383989); index must be a positive integer or logical.
是不是说程序中的wN要大于x的值呢？？

miao7mijao · 发表于 2006-10-23 09:33

大家好
我想问大家一个问题，就是说对一个信号时频变换以后，频率-幅值谱和频率-对数幅值谱不一样么？
我用x=10*sin（2*pi*50*t）；而在频率-对数幅值谱的52.4处也会出现幅值，那是为什么呢？

songzy41 · 发表于 2006-10-23 14:39

你把程序传上来看看。

realhappy · 发表于 2006-11-9 16:04

请问f=fl:df:fl+(N/2-1)*df;这个语句是否正确，fl是下限频率，fh才是上限频率，是否应该为f=fl:df:fh+(N/2-1)*df;呢

yangzj · 发表于 2006-11-9 17:28

原帖由 realhappy 于 2006-11-9 16:04 发表
请问f=fl:df:fl+(N/2-1)*df;这个语句是否正确，fl是下限频率，fh才是上限频率，是否应该为f=fl:df:fh+(N/2-1)*df;呢

没错的,应该来说
fl+(N/2-1)*df=fh
细化后的频谱就是fl到fh间的

realhappy · 发表于 2006-11-9 23:09

原帖由 yangzj 于 2006-11-9 17:28 发表

没错的,应该来说
fl+(N/2-1)*df=fh
细化后的频谱就是fl到fh间的

呵呵,明白了,事实上就是f=fl:df:fh;

cloudofsky · 发表于 2006-11-17 22:23

采样率fs=10240，采样点数102400，分辨率＝0.1，信号中两频率间距＝0.5Hz，直接做FFT即可区分求解两根谱线。感觉细化除了简化了计算以外，并无其它优势啊，点解？

cloudofsky · 发表于 2006-11-17 23:04

分辨率提高了，但是是以增加采样时间和存储空间为代价，瓶颈依然存在，很难适用于实时信号处理，也许只能在实验室使用吧。

Aprilgogo · 发表于 2007-1-18 17:03

有没有考虑过实现论文中说的两级降采样？这样可以达到更高的精度。

有一点疑惑的是经过第一级滤波后实信号已经变成有理信号(实部+虚部)了，那怎样再通过第二级的复解析滤波器呢？

Aprilgogo · 发表于 2007-1-18 17:05

似乎在滤波器的长度M的设置方面有点随意。应该符合M=4D/a最好吧^_^。

yangzj · 发表于 2007-1-19 08:54

原帖由 Aprilgogo 于 2007-1-18 17:03 发表
有没有考虑过实现论文中说的两级降采样？这样可以达到更高的精度。

有一点疑惑的是经过第一级滤波后实信号已经变成有理信号(实部+虚部)了，那怎样再通过第二级的复解析滤波器呢？

两级细化我也做过了,有兴趣的也可以试着做下

Aprilgogo · 发表于 2007-1-20 17:59

原帖由 yangzj 于 2007-1-19 08:54 发表

两级细化我也做过了,有兴趣的也可以试着做下

不知效果怎么样？不知能不能达到文中说的0.001Hz分辨率？

最近正在研究这份代码，等我也把两级细化的代码写一写再讨论吧^_^

yangzj · 发表于 2007-1-20 18:29

细化的效果还是很好的,不过正如大家所说的,时域长度还是要很长

zhulw · 发表于 2007-1-23 19:36

ZoomFFT程序调通了吗？我也需要该程序，跟我联系，
如果可以，我们可交换程序如何？

我的Email是vibsys@vip.sina.com朱先生

yixiaofan · 发表于 2007-4-12 22:21

function [f xz]=ZoomFFT(x,fs,N,fe,D,L,M)

k=1:M;
w=0.5+0.5*cos(pi*k/M); %Hanning窗 ???

我也看了丁康老师的文章了
Hanning窗不是w=0.5-0.5*cos(2*pi*k/M);吗?

账号		自动登录	找回密码
密码			我要加入

[FFT] [原创]zoomfft