文档

CZT

chirp z-transform

句法

y = czt(x,m,w,a)
y = czt(x)

描述

y = czt(x,m,w,a)返回信号的chirp z转换X。chirp z-transform是z转换X沿着螺旋轮廓定义w一个m是指定转换长度的标量w是点之间的比率z- 感兴趣的平面螺旋轮廓和标量一个是该轮廓上的复杂起点。轮廓,螺旋形或“ chirp”z- 平面,由

z = a*(w。^ - (0:m-1))

y = czt(x)使用以下默认值:

  • m=长度(x)

  • w=EXP(-j*2*pi/m)

  • 一个=1

使用这些默认值CZT返回Z变换Xm单位圆周围的分布点。这等同于离散的傅立叶变换X, 或者FFT(x)。空矩阵[]指定参数的默认值。

如果X是矩阵,czt(x,m,w,a)转换列X

例子

全部收缩

创建一个随机向量,X,长度1013。使用CZT

RNG默认x = randn(1013,1);y = czt(x);

利用CZTto zoom in on a narrow-band section of a filter's frequency response.

使用窗口方法设计30阶低通fir滤波器。指定1 kHz的样本率和125 Hz的截止频率。使用矩形窗口。找到过滤器的传输函数。

FS = 1000;d = distionfilt(“ Lowpassfir”,,,,“ FilterOrder”,30,“ cutofffrequency',125,...“ DesignMethod”,,,,'窗户',,,,'窗户',@rectwin,'采样率',fs);h = tf(d);

计算过滤器的DFT和CZT。将CZT的频率范围限制在100至150 Hz之间。在每种情况下生成1024个样品。

M = 1024;y = fft(h,m);F1 = 100;F2 = 150;w = exp(-j*2*pi*(f2-f1)/(m*fs));a = exp(j*2*pi*f1/fs);Z = CZT(H,M,W,A);

绘制转换。放大感兴趣的领域。

fn =(0:m-1)'/m;fy = fs*fn;fz =(f2-f1)*fn + f1;子图(2,1,1)图(FY,ABS(Y))轴([F1 F2 0 1.2])标题('FFT')子图(2,1,2)图(FZ,ABS(Z))轴([F1 F2 0 1.2])标题('CZT')Xlabel(“频率(Hz)”

诊断

如果m,,,,w, 或者一个不是标量CZT给出以下错误消息:

输入M,W,A必须是标量。

算法

CZT在计算时,请使用下一个2的长度FFT进行快速卷积z- 在指定的chirp轮廓上转换[1]

References

[1] Rabiner,Lawrence R.和Bernard Gold。数字信号处理的理论和应用。恩格尔伍德悬崖,台北:新世纪,1975,pp。393–399.

扩展功能

也可以看看

|

在R2006a之前引入