Chirp-Z 变换的理解_声音感知

个人资料

微博

正文字体大小：大中小

Chirp-Z 变换的理解

(2012-09-11 16:45:35)

标签：

分类：信号处理

Chirp-Z变换并不是一种新的变换，而是利用了傅立叶变换和卷积的性质。它也不能减小运算量，但可以计算某个频率段的傅立叶频谱。

公式如下：

http://s12/middle/6b6cc2bagfde1ba6e67cb&690变换的理解" TITLE="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 ')

参考文献：

阅读┊ 收藏 ┊ 喜欢 ▼ ┊打印┊举报/Report