加载中…[转载]Chirp-Z 变换的理解
(2018-03-28 20:49:40)
标签：
转载 |
分类： 技术点滴 |
!!!
参考文献：
原文地址：Chirp-Z 变换的理解 作者：声音感知
Chirp-Z变换并不是一种新的变换，而是利用了傅立叶变换和卷积的性质。它也不能减小运算量，但可以计算某个频率段的傅立叶频谱。
公式如下：
其出发点是将z写成A和W指数次乘积。
例如（MATLAB[2]）：
Establish frequency and CZT parameters:
fs = 1000; f1 = 100; f2 = 150; % in hertz
m = 1024;
w = exp(-j*2*pi*(f2-f1)/(m*fs));% 角度变换量 2*pi*(f2-f1)/(m*fs)
a = exp(j*2*pi*f1/fs);% 初始角度2*pi*f1/fs
Compute the frequency response of the filter using fft and czt:
y = fft(h,1000);
fy = (0:length(y)-1)'*1000/length(y);
subplot(211);
plot(fy(1:500),abs(y(1:500))); axis([1 500 0 1.2])
xlabel('Hz'); ylabel('Magnitude');
title('Magnitude Response using FFT')
z = czt(h,m,w,a);
fz = ((0:length(z)-1)'*(f2-f1)/length(z)) + f1;
subplot(212);
plot(fz,abs(z)); axis([f1 f2 0 1.2])
xlabel('Hz'); ylabel('Magnitude');
title('Magnitude Response using CZT ')
喜欢
0
赠金笔