cqt公司GydF4y2Ba

恒-Q非平稳Gabor变换GydF4y2Ba

折叠页面中的所有内容GydF4y2Ba

句法GydF4y2Ba

CFS = CQT（x）的GydF4y2Ba
[cfs，f]=cqt（x）GydF4y2Ba
[CFS，F，G，fshifts] = CQT（x）的GydF4y2Ba
[CFS，F，G，fshifts，fintervals] = CQT（x）的GydF4y2Ba
[cfs，f，g，fshifts，fintervals，bw]=cqt（x）GydF4y2Ba
[GydF4y2Ba___GydF4y2Ba]=cqt(GydF4y2Ba___GydF4y2Ba，名称，值）GydF4y2Ba
CQT（GydF4y2Ba___GydF4y2Ba）GydF4y2Ba

描述GydF4y2Ba

例GydF4y2Ba

CFSGydF4y2Ba= CQT（GydF4y2BaXGydF4y2Ba）GydF4y2Ba返回常量-Q变换（CQT），GydF4y2BaCFSGydF4y2Ba中，输入信号的GydF4y2BaXGydF4y2Ba。输入信号必须具有至少四个样品。GydF4y2Ba

  • 如果GydF4y2BaXGydF4y2Ba是一个向量，那么GydF4y2Bacqt公司GydF4y2Ba返回对应于CQT的矩阵。GydF4y2Ba

  • 如果GydF4y2BaXGydF4y2Ba是一个矩阵，然后GydF4y2Bacqt公司GydF4y2Ba获得的CQT的每一列（独立信道）GydF4y2BaXGydF4y2Ba. 该函数返回与CQT的最大冗余版本相对应的多维数组。GydF4y2Ba

例GydF4y2Ba

[GydF4y2BaCFSGydF4y2Ba，GydF4y2BaFGydF4y2Ba]=cqt(GydF4y2BaXGydF4y2Ba）GydF4y2Ba返回近似的带通中心频率，GydF4y2BaFGydF4y2Ba中，对应于各行GydF4y2BaCFSGydF4y2Ba。频率被排序从0到1并且周期/样本。GydF4y2Ba

[GydF4y2BaCFSGydF4y2Ba，GydF4y2BaFGydF4y2Ba，GydF4y2BaGGydF4y2Ba，GydF4y2BaF换档GydF4y2Ba]=cqt(GydF4y2BaXGydF4y2Ba）GydF4y2Ba返回的Gabor帧GydF4y2BaGGydF4y2Ba，用于分析GydF4y2BaXGydF4y2Ba频率变化，GydF4y2BaF换档GydF4y2Ba在离散傅立叶变换的通带之间（DFT）在仓中的行GydF4y2BaCFSGydF4y2Ba。GydF4y2Ba

CFSGydF4y2Ba，GydF4y2BaGGydF4y2Ba，和GydF4y2BaF换档GydF4y2Ba需要输入的CQT与反转GydF4y2BaicqtGydF4y2Ba。GydF4y2Ba

[GydF4y2BaCFSGydF4y2Ba，GydF4y2BaFGydF4y2Ba，GydF4y2BaGGydF4y2Ba，GydF4y2BaF换档GydF4y2Ba，GydF4y2Ba芬兰GydF4y2Ba]=cqt(GydF4y2BaXGydF4y2Ba）GydF4y2Ba返回频率间隔，GydF4y2Ba芬兰GydF4y2Ba，对应于GydF4y2BaCFSGydF4y2Ba。该GydF4y2BaķGydF4y2Ba的第i个元素GydF4y2BaF换档GydF4y2Ba是DFT容器中GydF4y2Ba（（K-1）MOD N）GydF4y2Ba和GydF4y2Ba（k模式N）GydF4y2Ba元素GydF4y2Ba芬兰GydF4y2Ba同GydF4y2BaK = 0,1,2，...，N-1GydF4y2Ba哪里GydF4y2BañGydF4y2Ba是频移的数量。由于MATLABGydF4y2Ba®GydF4y2Ba从1索引，GydF4y2BaF换档（1）GydF4y2Ba包含GydF4y2Bafintervals {端}GydF4y2Ba和GydF4y2Bafintervals {1}GydF4y2Ba，GydF4y2BaF换档（2）GydF4y2Ba包含GydF4y2Bafintervals {1}GydF4y2Ba和GydF4y2Ba结果{2}GydF4y2Ba， 等等。GydF4y2Ba

[GydF4y2BaCFSGydF4y2Ba，GydF4y2BaFGydF4y2Ba，GydF4y2BaGGydF4y2Ba，GydF4y2BaF换档GydF4y2Ba，GydF4y2Ba芬兰GydF4y2Ba，GydF4y2BaBWGydF4y2Ba]=cqt(GydF4y2BaXGydF4y2Ba）GydF4y2Ba返回带宽，GydF4y2BaBWGydF4y2Ba，在频率间隔的DFT仓，GydF4y2Ba芬兰GydF4y2Ba。GydF4y2Ba

例GydF4y2Ba

[GydF4y2Ba___GydF4y2Ba]=cqt(GydF4y2Ba___GydF4y2Ba，GydF4y2Ba名称，值GydF4y2Ba）GydF4y2Ba返回CQT与由一个或多个指定的附加选项GydF4y2Ba名称，值GydF4y2Ba配对参数，使用任何前述语法。GydF4y2Ba

例GydF4y2Ba

CQT（GydF4y2Ba___GydF4y2Ba）GydF4y2Ba不带输出参数绘制在目前的数字CQT。绘图仅支持矢量输入。万博1manbetx如果输入信号是真实的，GydF4y2BaFSGydF4y2Ba是采样频率，CQT被绘制在范围GydF4y2Ba[0，FS / 2]GydF4y2Ba. 如果信号很复杂，则在该范围内绘制CQT[GydF4y2Ba0，FSGydF4y2Ba）。GydF4y2Ba

注意GydF4y2Ba

为了可视化稀疏CQT，系数必须进行插值。当发生插值，情节可以有显著污渍和难以解释。如果你要绘制的CQT，我们建议使用默认GydF4y2Ba变形金刚GydF4y2Ba值GydF4y2Ba'充分'GydF4y2Ba。GydF4y2Ba

例子GydF4y2Ba

全部收缩GydF4y2Ba

打开实时脚本GydF4y2Ba

加载的信号，并获得恒定-Q变换。GydF4y2Ba

加载GydF4y2BanoisdoppGydF4y2BaCFS = CQT（noisdopp）;GydF4y2Ba
打开实时脚本GydF4y2Ba

加载一个实数值的信号，并获得恒定-Q变换。返回的近似带通中心频率。GydF4y2Ba

加载GydF4y2Ba汉德尔GydF4y2Ba[cfs，f]=cqt（y）；GydF4y2Ba

通过奈奎斯特频率绘制在对数刻度的带通中心频率。GydF4y2Ba

LFREQ =长度（F）;nyquistBin =地板（LFREQ / 2）+1;情节（F（1：nyquistBin））标题（GydF4y2Ba“带通中心频率”GydF4y2Ba）网格GydF4y2Ba在GydF4y2Ba集（GCA，GydF4y2Ba'yscale'GydF4y2Ba，GydF4y2Ba“日志”GydF4y2Ba）GydF4y2Ba

为了确认连续频率对的比率是恒定的，绘制比率图。自从GydF4y2Bacqt公司GydF4y2Ba默认使用每倍频程12个箱，该比值应等于GydF4y2Ba 2GydF4y2Ba 1GydF4y2Ba /GydF4y2Ba 1GydF4y2Ba 2GydF4y2Ba 。由于DC和奈奎斯特频率不是中心频率的几何序列的成员，但被包括在频率向量，从图中排除它们。GydF4y2Ba

图形绘制（f（3:nyquistBin-1）。/f（2:nyquistBin-2）网格GydF4y2Ba在GydF4y2Ba标题（[GydF4y2Ba“比：”GydF4y2Ba，num2str（2 ^（1/12））]）GydF4y2Ba

打开实时脚本GydF4y2Ba

获得音频信号的最小冗余恒Q变换。使用Blackman-Harris窗口作为Gabor框架的原型函数。GydF4y2Ba

加载GydF4y2Ba汉德尔GydF4y2BaDF = FS / numel（Y）;[CFS，F，G，fshifts，fintervals，BW] = CQT（Y，GydF4y2Ba'采样频率'GydF4y2Ba，FS，GydF4y2Ba'TransformType'GydF4y2Ba，GydF4y2Ba“疏”GydF4y2Ba，GydF4y2Ba'窗口'GydF4y2Ba，GydF4y2Ba“黑曼哈利斯”GydF4y2Ba）;GydF4y2Ba

CFSGydF4y2Ba是一个单元阵列，其中所述阵列中的每个元素对应于一个带通中心频率和Gabor帧。绘制与奈奎斯特频率相关联的伽柏帧。GydF4y2Ba

LF =长度（F）;IND =地板（LF / 2）+1;gFrame = fftshift（克{IND}）;fvec = F（IND-1）：DF：F（IND + 1）-df;图（fvec，gFrame）xlabel（GydF4y2Ba'频率（Hz）'GydF4y2Ba）网格GydF4y2Ba在GydF4y2Ba标题（{[GydF4y2Ba'的Gabor框架 - 频率：'GydF4y2Ba，num2str（F（IND）），GydF4y2Ba“赫兹”GydF4y2Ba]; [GydF4y2Ba'带宽'GydF4y2Ba，num2str（BW（IND）* FS / numel（Y）），GydF4y2Ba“赫兹”GydF4y2Ba]}）GydF4y2Ba

在恒Q变换中，将Gabor框架应用于输入信号的离散Fourier变换，并进行离散Fourier逆变换。第k个Gabor帧应用于GydF4y2Ba芬兰GydF4y2Ba。以离散傅立叶变换的信号并绘制其幅度谱。采用GydF4y2Ba芬兰GydF4y2Ba以指示哪些傅里叶系数是与奈奎斯特频率相关联的Gabor帧。GydF4y2Ba

yDFT = FFT（y）基lyDFT =长度（yDFT）;情节（FS *（0：lyDFT-1）/ lyDFT，ABS（yDFT））格GydF4y2Ba在GydF4y2BafIntervalGabor = fintervals {IND};MX = MAX（ABS（yDFT））;保持GydF4y2Ba在GydF4y2Ba情节（[DF * fIntervalGabor（1）DF * fIntervalGabor（1）]，[0 MX]GydF4y2Ba'R-'GydF4y2Ba，GydF4y2Ba'行宽'GydF4y2Ba，2）图（[DF * fIntervalGabor（结束）DF * fIntervalGabor（结束）]，[0 MX]GydF4y2Ba'R-'GydF4y2Ba，GydF4y2Ba'行宽'GydF4y2Ba，2）海峡= sprintf的（GydF4y2Ba'的Gabor帧间隔（赫兹）：[％3.2F，3.2F％]'GydF4y2Ba，DF * fIntervalGabor（1），DF * fIntervalGabor（结束））;标题（STR）GydF4y2Ba

利用Gabor框架对区间内的Fourier系数进行窗口化，并进行逆离散Fourier变换。对计算结果进行归一化处理，并与计算出的常Q系数进行比较，确定它们是相等的。GydF4y2Ba

lGframe =长度（gFrame）;INDX = 1：lGframe;INDX = fftshift（INDX）;winDFT（INDX）= yDFT（fIntervalGabor）* fftshift（gFrame（INDX））;cqCoefs = IFFT（winDFT）;cqCoefs =（2 * lGframe /长度（Y））* cqCoefs;MAX（ABS（cqCoefs（：） -  CFS {} IND（:)））GydF4y2Ba
ANS = 0GydF4y2Ba
打开实时脚本GydF4y2Ba

加载音频信号。使用变换的最大冗余版本和每倍频程12个二进制数绘制常量Q变换（CQT）。GydF4y2Ba

加载GydF4y2Ba汉德尔GydF4y2Bacqt（是，GydF4y2Ba'采样频率'GydF4y2Ba，FS）GydF4y2Ba

执行使用每倍频程48个二进制位相同的信号的CQT。设置在其上CQT具有对数的频率响应是可允许的最低频率到2kHz的频率范围内。GydF4y2Ba

minFreq = FS /长度（Y）;maxFreq = 2000;图CQT（Y，GydF4y2Ba'采样频率'GydF4y2Ba，FS，GydF4y2Ba'BinsPerOctave'GydF4y2Ba，48岁，GydF4y2Ba'FrequencyLimits'GydF4y2Ba，[最小频率最大频率]）GydF4y2Ba

输入参数GydF4y2Ba

全部收缩GydF4y2Ba

输入信号，指定为实数或复数矢量或矩阵。GydF4y2BaXGydF4y2Ba必须至少有四个样本。GydF4y2Ba

数据类型：GydF4y2Ba单一的GydF4y2Ba|GydF4y2Ba双GydF4y2Ba
复数支持：万博1manbetxGydF4y2Ba是GydF4y2Ba

名称 - 值对参数GydF4y2Ba

指定可选的用逗号分隔的对GydF4y2Ba名称，值GydF4y2Ba参数。GydF4y2Ba名称GydF4y2Ba是参数的名称和GydF4y2Ba价值GydF4y2Ba是相应的值。GydF4y2Ba名称GydF4y2Ba必须出现引号内。您可以按照任何顺序指定多个名称和值对参数GydF4y2Ba名1，值1，...，NameN，值NGydF4y2Ba。GydF4y2Ba

例：GydF4y2Ba'采样频率'，20，'采样频率'，15GydF4y2Ba

采样频率，单位为赫兹，指定为逗号分隔对，由GydF4y2Ba'采样频率'GydF4y2Ba一个正标量。GydF4y2Ba

每倍频程频段在CQT使用，指定为正整数从1到96的数。GydF4y2Ba

恒定Q值的类型的变换来执行，指定为逗号分隔的一对组成的GydF4y2Ba'TransformType'GydF4y2Ba和GydF4y2Ba'充分'GydF4y2Ba或GydF4y2Ba'疏'GydF4y2Ba. 稀疏变换是常数Q变换的最小冗余形式。GydF4y2Ba

频率限制在其上CQT与每倍频频率仓的指定数目的对数频率响应，指定为逗号分隔的一对GydF4y2Ba'FrequencyLimits'GydF4y2Ba和两元件实数向量。GydF4y2Ba

  • 第一个元素必须大于或等于GydF4y2BaFS / NGydF4y2Ba，其中GydF4y2BaFSGydF4y2Ba是采样频率和GydF4y2BañGydF4y2Ba是信号的长度。GydF4y2Ba

  • 第二元素必须是严格大于奈奎斯特频率以下。GydF4y2Ba

窗口用作非平稳Gabor框架的原型函数，指定为GydF4y2Ba“翰”GydF4y2Ba，GydF4y2Ba“汉明”GydF4y2Ba，GydF4y2Ba'blackmanharris'GydF4y2Ba，GydF4y2Ba'itersine'GydF4y2Ba，或GydF4y2Ba“巴特利特GydF4y2Ba. 这些紧支撑函数是在频率上定义的。万博1manbetx对于规范化频率，它们在间隔上定义GydF4y2Ba（-1/2,1/2）GydF4y2Ba. 如果指定采样频率，GydF4y2BaFSGydF4y2Ba，它们在所定义的时间间隔GydF4y2Ba（-Fs/2，Fs/2）GydF4y2Ba。GydF4y2Ba

输出参数GydF4y2Ba

全部收缩GydF4y2Ba

恒定-Q变换，返回作为基质，多维数组，单元阵列，或结构阵列。GydF4y2Ba

  • 如果GydF4y2Ba'TransformType'GydF4y2Ba指定为GydF4y2Ba'充分'GydF4y2Ba无GydF4y2Ba'FrequencyLimits'GydF4y2Ba，GydF4y2BaCFSGydF4y2Ba是一个矩阵或多维数组。GydF4y2Ba

    • 如果GydF4y2BaXGydF4y2Ba是一个向量，那么GydF4y2Bacqt公司GydF4y2Ba返回对应于CQT的矩阵。GydF4y2Ba

    • 如果GydF4y2BaXGydF4y2Ba是一个矩阵，然后GydF4y2Bacqt公司GydF4y2Ba获得的CQT的每一列（独立信道）GydF4y2BaXGydF4y2Ba. 该函数返回与CQT的最大冗余版本相对应的多维数组。GydF4y2Ba

    的阵列，GydF4y2BaCFSGydF4y2Ba，对应于CQT的最大冗余版本。每一行GydF4y2BaCFSGydF4y2Ba对应于具有1/2对应于奈奎斯特频率0和1的归一化频率之间对数间隔归一化的中心频率（周期/样品）的通带。列，或啤酒花，对应于最大带宽的中心频率，这通常发生一个频率窗口的下方或奈奎斯特仓以上的数量。GydF4y2Ba

  • 如果GydF4y2Ba'TransformType'GydF4y2Ba指定为GydF4y2Ba'充分'GydF4y2Ba并指定频率限制，GydF4y2BaCFSGydF4y2Ba作为具有以下四个字段的结构数组返回。GydF4y2Ba

    • CGydF4y2Ba- 为对特定频率范围内频率多维数组的系数矩阵。这既包括正和“负”的频率。GydF4y2Ba

    • DCcfsGydF4y2Ba-通带的系数向量或矩阵，从0到频率下限。GydF4y2Ba

    • 奈奎斯特GydF4y2Ba- 用于从所述频率上限奈奎斯特通带系数矢量或矩阵。GydF4y2Ba

    • 奈奎斯特宾GydF4y2Ba- 对应于奈奎斯特频率DFT仓。反转CQT时，该字段被使用。GydF4y2Ba

  • 如果GydF4y2Ba'TransformType'GydF4y2Ba指定为GydF4y2Ba'疏'GydF4y2Ba，GydF4y2BaCFSGydF4y2Ba是用元件的数量的单元阵列等于带通频率的数量。单元阵列中的每个元素，GydF4y2BaCFSGydF4y2Ba是矢量或矩阵的行数等于所述带宽在DFT仓值，GydF4y2BaBWGydF4y2Ba。GydF4y2Ba

CFSGydF4y2Ba，GydF4y2BaGGydF4y2Ba，和GydF4y2BaF换档GydF4y2Ba需要输入的CQT与反转GydF4y2BaicqtGydF4y2Ba。GydF4y2Ba

对应于以下行的近似带通中心频率GydF4y2BaCFSGydF4y2Ba，作为实值向量返回。频率从0到1排列，以周期/采样为单位。如果您指定GydF4y2Ba'采样频率'GydF4y2Ba， 然后GydF4y2BaFGydF4y2Ba在赫兹。GydF4y2Ba

Gabor框架在结构分析中的应用GydF4y2BaXGydF4y2Ba，返回为实值向量的单元阵列。在每个向量GydF4y2BaGGydF4y2Ba对应的行GydF4y2BaCFSGydF4y2Ba。GydF4y2Ba

CFSGydF4y2Ba，GydF4y2BaGGydF4y2Ba，和GydF4y2BaF换档GydF4y2Ba需要输入的CQT与反转GydF4y2BaicqtGydF4y2Ba。GydF4y2Ba

在离散傅里叶频移变换箱，返回为实值向量。所述位移在的行中的通带之间GydF4y2BaCFSGydF4y2Ba。GydF4y2Ba

CFSGydF4y2Ba，GydF4y2BaGGydF4y2Ba，和GydF4y2BaF换档GydF4y2Ba需要输入的CQT与反转GydF4y2BaicqtGydF4y2Ba。GydF4y2Ba

与行对应的频率间隔GydF4y2BaCFSGydF4y2Ba，返回为单元阵列。中的每个元素GydF4y2Ba芬兰GydF4y2Ba是实值向量。该GydF4y2BaķGydF4y2Ba的第i个元素GydF4y2BaF换档GydF4y2Ba是DFT容器中GydF4y2Ba（（K-1）MOD N）GydF4y2Ba和GydF4y2Ba（k模式N）GydF4y2Ba元素GydF4y2Ba芬兰GydF4y2Ba同GydF4y2BaK = 0,1,2，...，N-1GydF4y2Ba哪里GydF4y2BañGydF4y2Ba是频移的数量。因为从1 MATLAB指标，GydF4y2BaF换档（1）GydF4y2Ba包含GydF4y2Bafintervals {端}GydF4y2Ba和GydF4y2Bafintervals {1}GydF4y2Ba，GydF4y2BaF换档（2）GydF4y2Ba包含GydF4y2Bafintervals {1}GydF4y2Ba和GydF4y2Ba结果{2}GydF4y2Ba， 等等。GydF4y2Ba

在频率间隔的DFT仓带宽GydF4y2Ba芬兰GydF4y2Ba，作为实值向量返回。GydF4y2Ba

算法GydF4y2Ba

全部收缩GydF4y2Ba

非平稳Gabor框架GydF4y2Ba

用于频率自适应分析的非平稳Gabor（NSG）帧理论和使用NSG帧进行分析和合成的有效算法是由Dórfler、Holighaus、Grill和Velasco提出的GydF4y2Ba[1]GydF4y2Ba，GydF4y2Ba[2]GydF4y2Ba。在CQT和ICQT使用的算法是由Dörfler，Holighaus，烤肉，和贝拉斯科开发的并且描述于GydF4y2Ba[1]GydF4y2Ba，GydF4y2Ba[2]GydF4y2Ba. 在GydF4y2Ba[3]GydF4y2Ba，Schörkhuber，Klapuri，Holighaus，和Dörfler开发和用于提供算法的经相位校正的CQT变换其将由幼稚卷积来获得CQT系数相匹配。大时频分析工具箱（GydF4y2Bahttps://github.com/ltfatGydF4y2Ba）提供了广泛的算法套件非平稳的Gabor帧GydF4y2Ba[4]GydF4y2Ba。GydF4y2Ba

完美重建GydF4y2Ba

与非平稳的Gabor帧实现的恒定Q值分析的完全重构属性，GydF4y2Bacqt公司GydF4y2Ba在内部预先设置零频率（DC），并将奈奎斯特频率附加到频率间隔。负频率是正中心频率和带宽的镜像版本GydF4y2Ba

参考GydF4y2Ba

[1] Holighaus，N.，M。Dörfler，G. A.贝拉斯科和T.烧烤。“为可逆实时恒Q值变换的框架。”GydF4y2BaIEEE交易对音频，语音和语言处理。GydF4y2Ba卷。21，第4号，2013，第775-785。GydF4y2Ba

[2] Velasco、G.A.，N.Holighaus、M.Dórfler和T.Grill。”用非平稳Gabor框架构造一个可逆的常Q变换GydF4y2Ba第14届国际会议上的数字音效（DAFX-11）的诉讼GydF4y2Ba。法国巴黎：2011。GydF4y2Ba

[3] 舍尔库伯，C.，A.克拉普里，N.霍利加斯和M.德弗勒。”一个Matlab工具箱，用于以对数频率分辨率高效完美地重建时频变换GydF4y2BaAES第53届国际会议上语义音频GydF4y2Ba. 英国伦敦：2014年。GydF4y2Ba

[4] 普罗沙，Z.，P.L.Søndergaard，N.Holighaus，C.Wiesmeyr和P.Balazs。GydF4y2Ba大时频分析工具箱2.0GydF4y2Ba。声音，音乐和运动，在计算机科学2014年，第419-442讲义。GydF4y2Ba

也可以看看GydF4y2Ba

主题GydF4y2Ba

介绍了在R2018aGydF4y2Ba

小波工具箱文档GydF4y2Ba
万博1manbetx

尝试MATLAB、Sim万博1manbetxulink和其他产品s manbetx 845GydF4y2Ba

现在就试用GydF4y2Ba