主要内容

CZT.

Chirp Z-Transform

页面上倒塌

句法

y = czt（x，m，w，a）

描述

例子

y= czt（X那m那W.那一种）返回长度 -mchirp z-transform（czt）的X沿着螺旋形轮廓Z.-plane定义W.和一种通过z =一种*W.。^ - （0：m-1）。

默认值m那W.， 和一种那CZT.返回z变换X在m单位圈周围等间隔点，结果相当于离散的傅里叶变换（DFT）X据给出FFT.（X）。

例子

全部收缩

打开直播脚本

创建一个随机向量，X，长度为1013.计算其DFT使用CZT.。

RNG.默认x = RANDN（1013,1）;y = czt（x）;
打开直播脚本

采用CZT.要放大滤波器频率响应的窄带部分。

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

FS = 1000;d = designfilt（'低通道'那'筛选道'，30，'cutfffrequency'，125，......'DesignMethod'那'窗户'那'窗户'，@ entrwin，'采样率'，fs）;h = tf（d）;

计算过滤器的DFT和CZT。将CZT的频率范围限制为75到175 Hz之间的带。在每种情况下生成1024个样本。

m = 1024;Y = FFT（H，M）;F1 = 75;F2 = 175;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;图（FY，ABS（Y），FZ，ABS（Z））XLIM（[50 200]）传奇（'fft'那'czt'）xlabel（'频率（Hz）'）

图包含轴。轴包含2个类型的型号。这些物体代表FFT，CZT。

输入参数

全部收缩

输入信号，指定为向量，矩阵或3-D阵列。如果X是一个矩阵，函数转换列的列X。如果X是3-D阵列，该功能沿着第一个阵列尺寸运行，大小大于1。

例子：SIN（PI ./ [4; 2] *（0：159））'指定双通道正弦曲线。

数据类型：单身的|双倍的
复数支持：金宝app是的

变换长度，指定为正整数标量。

数据类型：单身的|双倍的

螺旋轮廓点之间的比率，指定为复标量。

数据类型：单身的|双倍的
复数支持：金宝app是的

螺旋形轮廓初始点，指定为复标量。

例子：EXP（1J * PI / 4）沿着单位圈Z.-Plane并使真实轴与45度的角度。

数据类型：单身的|双倍的
复数支持：金宝app是的

输出参数

全部收缩

Chirp z变换，作为向量或矩阵返回。

算法

CZT.在指定的Chirp轮廓上计算z变换时，使用下一个电源-2长度FFT来执行快速卷积[1]。

参考

[1] Rabiner，Lawrence R.和Bernard Gold。数字信号处理的理论与应用。Englewood Cliffs，NJ：Prentice-Hall，1975年。

扩展能力

也可以看看

|

话题

在R2006A之前介绍

信号处理工具箱文档

金宝app

利用MATLAB的信号处理深入学习

利用MATLAB的信号处理深入学习

下载白皮书