(1)求h(n)的主值序列
(2)用FFT求 的付里叶变换 H(k)=FFT[ ] L点
(3)对x(n)加权并补零
g(n)=
(4)G(k)=FFT[g(n)]L点
(5)Y(k)=G(k)H(k) L点
(6)y(n)=IFFT[Y(k)] L点
(7) 0≤k≤M-1
乘法计算量估算:
(1)、(2)两步可以事先计算,不必实时计算。(3)、(7)两步两次加权,共计N+M次复乘。形成Y(k),需L次复乘。一个FFT与一个IFFT需Llog2L次乘法。所以,总乘法数为L+N+M+Llog2L,而直接计算乘法数为NM。可见,N及M较大时,用FFT实现Chirp-Z变换,速度上有很大的改进。
Chirp-z变换的特点:(与FFT相比)
(1)输入序列的长度N与输出序列的长度M不需要相等;
(2)N及M不必是合成数,二者均可为素数;
(3)zk点的角间隔φ0是任意的,因此频率分辨率也是任意的;
(4)周线不必是z平面上的圆;
(5)起始点z0可任意选定,因此可从任意频率上开始对输入数据进行窄带高分辨率力分析;
(6)若A=1,M=N,则可用Chirp-z变换计算DFT,即使N为素数。
需要处理成什么?
chirp-z是频谱分析工具, 能直接得出频谱.
10 kHz载频
50 Hz 调幅信号
最上面是调制波形局部.
中间是幅度频谱, 可以看到,在9950Hz和10050Hz位置的都有一个峰值, 这是调幅的特点.
最下面是相位频谱.
欢迎分享,转载请注明来源:内存溢出
评论列表(0条)