相位补偿细化,可以对全部频率范围内的频谱进行细化,这就克服了频移法的缺点。当然,对于只需要在窄带范围内细化的情况,用相位补偿法有点浪费。
设要求的细化因子为D=2^n,则采样个数为:N’=DN=DT/t,式中T—原分析长度,t—采样间隔。
将相距D个采样间隔的样本抽出来集合为一个子序列,每个子序列有N个样本,共有D个子序列。总频谱是D次N点的FFT结果Xd相应乘以W^k DN后叠加的结果。这里的W^kDN即为谱线Xd的相位补偿量。
[ 本帖最后由 zhangnan3509 于 2007-6-29 21:47 编辑 ] 3级联Zoom-FFT
假定样本数据是x(n+mN),n={0,…,N-1;m=0,…, M-1。其中数据序列被划分为M个区组,互不重叠。每一个区组有N个样本点。N*M点FFT的频率分辨率是:f=1/NMTS,运算量为NM log(NM)/2。
N*M点的DFT可以简化成两次DFT,第一次DFT是对M个区组作N点DFT,而第二次DFT是对所关心的谱线k做M点DFT。总计算量为:(MnlogN+MlogM)/2,比上述计算量减少(N-1)MlogM/2。级联FFT与复包络解调法其实在本质上是类似的。 原帖由 微尼 于 2007-6-29 18:05 发表
请问songzy41,这个频谱细化程序在使用resample函数重新采样之前是不是没有滤波? 在resample函数重新采样之前没有通过滤波器,在上程序中是利用了resample函数自带的滤波功能:
RESAMPLE applies an anti-aliasing (lowpass) FIR filter to X during the
resampling process, and compensates for the filter's delay.The filter
is designed using FIRLS. 在重采样的过程中,采样的总点数也自然随之降代低了,fs=fs/N,n=n/n,请教
频率分辨率是如何提高的
回复 #19 toney 的帖子
从这个方向去理解的话,那ZFFT的方法并不能提高分辨率.看来很多人对ZFFT有个误区,事实上,ZFFT所谓的提高分辨率是指针对作同样点数的离散傅立叶变换而言的,即作细化D倍的N点细化谱实际上原始数据必须采DN点数据,这时候它比用N点原始数据直接做FFT而言频率分辨率提高了D倍,但如果把DN点原始数据全部做FFT,那它和细化谱的分辨率是一样的.
一句话,ZFFT本质上可以说是一种快速算法,它通过滤波重采样来降低采样频率,这样就可以用较少点数的FFT来实现较高的频率分辨率,当然,提高速度的代价就是只能对局部频带进行细化(而如果将ZFFT利用的所有原始数据全部直接做FFT的话,它做出的是整个频域的,而且频率分辨率和细化后的一样,甚至如果考虑细化时滤波所需去掉的点,直接FFT的频率分辨率可以更高).
貌似类似的我说过很多遍了.:lol 太好了,讲解的太好了,:victory: :victory: 各位大虾的水平太牛了{:{39}:} 不错,收获颇丰!{:{03}:} 回复 6 # songzy41 的帖子
你好,请问互功率谱函数程序有没有?现实条件是知道时间序列,我想求这个时间序列和其延时序列的互功率谱,怎么办?谢谢了 学习中,对自己很有用
页:
1
[2]