模仿频谱细化matlab程序

zhangnan3509

2相位补偿细化
相位补偿细化，可以对全部频率范围内的频谱进行细化，这就克服了频移法的缺点。当然，对于只需要在窄带范围内细化的情况，用相位补偿法有点浪费。
设要求的细化因子为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 编辑 ]

zhangnan3509

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与复包络解调法其实在本质上是类似的。

songzy41

原帖由 微尼 于 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.

toney

在重采样的过程中，采样的总点数也自然随之降代低了，fs=fs/N,n=n/n,请教
频率分辨率是如何提高的

yangzj

从这个方向去理解的话,那ZFFT的方法并不能提高分辨率.
看来很多人对ZFFT有个误区,事实上,ZFFT所谓的提高分辨率是指针对作同样点数的离散傅立叶变换而言的,即作细化D倍的N点细化谱实际上原始数据必须采DN点数据,这时候它比用N点原始数据直接做FFT而言频率分辨率提高了D倍,但如果把DN点原始数据全部做FFT,那它和细化谱的分辨率是一样的.
一句话,ZFFT本质上可以说是一种快速算法,它通过滤波重采样来降低采样频率,这样就可以用较少点数的FFT来实现较高的频率分辨率,当然,提高速度的代价就是只能对局部频带进行细化(而如果将ZFFT利用的所有原始数据全部直接做FFT的话,它做出的是整个频域的,而且频率分辨率和细化后的一样,甚至如果考虑细化时滤波所需去掉的点,直接FFT的频率分辨率可以更高).

貌似类似的我说过很多遍了.:lol

shenqingwo

太好了，讲解的太好了，:victory: :victory:

200810100659

各位大虾的水平太牛了

whutzhq

不错，收获颇丰！

chngxu

回复 6 # songzy41 的帖子

你好，请问互功率谱函数程序有没有？现实条件是知道时间序列，我想求这个时间序列和其延时序列的互功率谱，怎么办？谢谢了

zhoulid2010

学习中，对自己很有用