主要内容GydF4y2Ba

cqtGydF4y2Ba

恒定Q非营养性Gabor变换GydF4y2Ba

描述GydF4y2Ba

例子GydF4y2Ba

CFS.GydF4y2Ba= CQT(GydF4y2BaXGydF4y2Ba)GydF4y2Ba返回常数Q变换(CQT),GydF4y2BaCFS.GydF4y2Ba,输入信号GydF4y2BaXGydF4y2Ba.输入信号必须具有至少四个样本。GydF4y2Ba

  • 如果GydF4y2BaXGydF4y2Ba是向量吗GydF4y2BacqtGydF4y2Ba返回与CQT对应的矩阵。GydF4y2Ba

  • 如果GydF4y2BaXGydF4y2Ba是一个矩阵GydF4y2BacqtGydF4y2Ba获取每个列(独立频道)的CQTGydF4y2BaXGydF4y2Ba.该函数返回一个与CQT的最大冗余版本对应的多维数组。GydF4y2Ba

例子GydF4y2Ba

[GydF4y2BaCFS.GydF4y2Ba那GydF4y2BaFGydF4y2Ba) = cqt (GydF4y2BaXGydF4y2Ba)GydF4y2Ba返回近似的带通中心频率,GydF4y2BaFGydF4y2Ba,对应于行的行GydF4y2BaCFS.GydF4y2Ba.频率从0到1排序,并在周期/样本。GydF4y2Ba

[GydF4y2BaCFS.GydF4y2Ba那GydF4y2BaFGydF4y2Ba那GydF4y2BaGGydF4y2Ba那GydF4y2BafshiftsGydF4y2Ba) = cqt (GydF4y2BaXGydF4y2Ba)GydF4y2Ba返回Gabor框架,GydF4y2BaGGydF4y2Ba,用于分析GydF4y2BaXGydF4y2Ba频率会转移,GydF4y2BafshiftsGydF4y2Ba,以离散的傅里叶变换(DFT)频带之间的通带中的行GydF4y2BaCFS.GydF4y2Ba.GydF4y2Ba

CFS.GydF4y2Ba那GydF4y2BaGGydF4y2Ba,GydF4y2BafshiftsGydF4y2Ba需要输入CQT的反转GydF4y2BaICQT.GydF4y2Ba.GydF4y2Ba

[GydF4y2BaCFS.GydF4y2Ba那GydF4y2BaFGydF4y2Ba那GydF4y2BaGGydF4y2Ba那GydF4y2BafshiftsGydF4y2Ba那GydF4y2BafintervalsGydF4y2Ba) = cqt (GydF4y2BaXGydF4y2Ba)GydF4y2Ba返回频率间隔,GydF4y2BafintervalsGydF4y2Ba,对应的行GydF4y2BaCFS.GydF4y2Ba.这GydF4y2BaK.GydF4y2BaTh元素GydF4y2BafshiftsGydF4y2BaDFT箱的频移是GydF4y2Ba((k-1)mod n)GydF4y2Ba和GydF4y2Ba(k mod N)GydF4y2Ba元素GydF4y2BafintervalsGydF4y2Ba和GydF4y2Bak = 0,1,2,...,n-1GydF4y2Ba在哪里GydF4y2BaNGydF4y2Ba是频率偏移的数量。因为matlab.GydF4y2Ba®GydF4y2Ba索引1,GydF4y2Bafshifts (1)GydF4y2Ba包含之间的频移GydF4y2Bafintervals {结束}GydF4y2Ba和GydF4y2Bafintervals {1}GydF4y2Ba那GydF4y2Bafshifts (2)GydF4y2Ba包含之间的频移GydF4y2Bafintervals {1}GydF4y2Ba和GydF4y2Bafintervals {2}GydF4y2Ba, 等等。GydF4y2Ba

[GydF4y2BaCFS.GydF4y2Ba那GydF4y2BaFGydF4y2Ba那GydF4y2BaGGydF4y2Ba那GydF4y2BafshiftsGydF4y2Ba那GydF4y2BafintervalsGydF4y2Ba那GydF4y2BaBW.GydF4y2Ba) = cqt (GydF4y2BaXGydF4y2Ba)GydF4y2Ba返回带宽,GydF4y2BaBW.GydF4y2Ba,在频率间隔的DFT区间,GydF4y2BafintervalsGydF4y2Ba.GydF4y2Ba

例子GydF4y2Ba

[GydF4y2Ba___GydF4y2Ba) = cqt (GydF4y2Ba___GydF4y2Ba那GydF4y2Ba名称,价值GydF4y2Ba)GydF4y2Ba返回带有由一个或多个指定的附加选项的CQTGydF4y2Ba名称,价值GydF4y2Ba配对参数,使用前面的任何语法。GydF4y2Ba

例子GydF4y2Ba

cqt(GydF4y2Ba___GydF4y2Ba)GydF4y2Ba没有输出参数在当前数字中绘制CQT。仅支持矢量输入的绘图。金宝app如果输入信号是真实的GydF4y2BaFS.GydF4y2Ba是采样频率,CQT在范围内绘制GydF4y2Ba[0,FS / 2]GydF4y2Ba.如果信号是复杂的,则CQT在范围[GydF4y2Ba0, FsGydF4y2Ba)。GydF4y2Ba

笔记GydF4y2Ba

为了可视化稀疏CQT,必须内容系数。当发生插值时,曲线可能具有显着的涂抹并且难以解释。如果您想绘制CQT,我们建议使用默认值GydF4y2BaTransformTypeGydF4y2Ba价值GydF4y2Ba“全部”GydF4y2Ba.GydF4y2Ba

例子GydF4y2Ba

全部收缩GydF4y2Ba

加载信号并获得常数Q变换。GydF4y2Ba

加载GydF4y2Banoisdopp.GydF4y2Bacfs = cqt(noisdopp);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

图中包含一个坐标轴。标题为“带通中心频率”的轴包含一个类型为line的对象。GydF4y2Ba

为了确定连续频率对的比率是常数,绘制比率图。自GydF4y2BacqtGydF4y2Ba默认使用每个八度的12个箱,比率应该等于GydF4y2Ba 2GydF4y2Ba 1GydF4y2Ba /GydF4y2Ba 1GydF4y2Ba 2GydF4y2Ba .由于DC和NYQUIST频率不是中心频率的几何序列的成员,而是包括在频率向量中,从绘图中排除它们。GydF4y2Ba

图绘制(f (3: nyquistBin-1)。/ f (2: nyquistBin-2))网格GydF4y2Ba在GydF4y2Ba标题([GydF4y2Ba'比率: 'GydF4y2Ba,num2str(2 ^(1/12))])GydF4y2Ba

图中包含一个坐标轴。具有标题比率的轴:1.0595包含类型线的对象。GydF4y2Ba

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

加载GydF4y2Ba汉德尔GydF4y2Badf = fs / numel(y);[CFS,F,G,FShift,Fintervals,BW] = CQT(Y,GydF4y2Ba'采样频率'GydF4y2Ba,fs,GydF4y2Ba'transformype'GydF4y2Ba那GydF4y2Ba“疏”GydF4y2Ba那GydF4y2Ba'窗户'GydF4y2Ba那GydF4y2Ba“blackmanharris”GydF4y2Ba);GydF4y2Ba

CFS.GydF4y2Ba是一个小区阵列,其中阵列中的每个元素对应于带通中心频率和Gabor帧。绘制与奈奎斯特频率相关的Gabor帧。GydF4y2Ba

lf =长度(f);IND =地板(LF / 2)+1;gframe = fftshift(g {ind});FVEC = F(IND-1):DF:F(IND + 1)-DF;绘图(FVEC,GFRAME)XLABEL(GydF4y2Ba'频率(Hz)'GydF4y2Ba) 网格GydF4y2Ba在GydF4y2Ba标题({[GydF4y2Ba'Gabor Frame  -  Freq:'GydF4y2Ba,num2str(f(ind)),GydF4y2Ba'赫兹'GydF4y2Ba]; [GydF4y2Ba“带宽”GydF4y2Ba,num2str(bw(ind)* fs / numel(y)),GydF4y2Ba'赫兹'GydF4y2Ba]})GydF4y2Ba

图中包含一个坐标轴。标题为Gabor Frame - Freq: 4096 Hz Bandwidth 412.3283 Hz的轴包含一个类型为line的对象。GydF4y2Ba

在常q变换中,将Gabor框架应用于输入信号的离散傅里叶变换,并进行离散傅里叶反变换。将第k个Gabor框架应用于规定的第k个频率区间GydF4y2BafintervalsGydF4y2Ba.采用信号的离散傅立叶变换并绘制其幅度谱。采用GydF4y2BafintervalsGydF4y2Ba以指示与奈奎斯特频率相关联的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(end)df * fintervalgabor(end)],[0 mx],GydF4y2Ba'r-'GydF4y2Ba那GydF4y2Ba'行宽'GydF4y2Ba,2)str = sprintf(GydF4y2Ba'Gabor帧间隔(Hz):[%3.2f,%3.2f]'GydF4y2Ba,df * fintervalgabor(1),df * fintervalgabor(端));标题(str)GydF4y2Ba

图中包含一个坐标轴。具有标题Gabor帧间隔(Hz)的轴:[3889.89,4302.11]包含3个类型的线。GydF4y2Ba

将傅里叶系数在Gabor坐标系的区间内加窗,并进行离散傅里叶反变换。将结果归一化,并与计算出的常q系数进行比较,确认它们相等。GydF4y2Ba

lgframe =长度(gframe);Indx = 1:LGFrame;Indx = FFTShift(Indx);Windft(Indx)= YDFT(FINTERVALGABOR)。* FFTSHIFT(GFRAME(INDX));cqcoefs = ifft(windft);cqcoefs =(2 * lgframe / length(y))* cqcoefs;max(abs(cqcoefs(:)  -  cfs {ind}(:)))GydF4y2Ba
ans = 0.GydF4y2Ba

加载音频信号。使用变换的最大冗余版本并使用每个八度音阶的12个bins绘制常量q变换(CQT)。GydF4y2Ba

加载GydF4y2Ba汉德尔GydF4y2Bacqt (y,GydF4y2Ba'采样频率'GydF4y2Ba,fs)GydF4y2Ba

图中包含一个坐标轴。具有标题常数Q变换的轴包含类型表面的对象。GydF4y2Ba

使用每个八度音程48个箱执行相同信号的CQT。将CQT的频率范围设置为对数频率响应的频率范围,以使最小允许频率为2 kHz。GydF4y2Ba

minfreq = fs / length(y);maxfreq = 2000;图cqt(y,GydF4y2Ba'采样频率'GydF4y2Ba,fs,GydF4y2Ba'binsperoctave'GydF4y2Ba现年48岁的GydF4y2Ba'surformlimits'GydF4y2Ba, (minFreq maxFreq])GydF4y2Ba

图中包含一个坐标轴。具有标题常数Q变换的轴包含类型表面的对象。GydF4y2Ba

输入参数GydF4y2Ba

全部收缩GydF4y2Ba

输入信号,指定为真实或复杂的向量或矩阵。GydF4y2BaXGydF4y2Ba必须至少有四个样本。GydF4y2Ba

数据类型:GydF4y2Ba双倍的GydF4y2Ba
复数支持:金宝appGydF4y2Ba是的GydF4y2Ba

名称值对参数GydF4y2Ba

指定可选的逗号分离对GydF4y2Ba名称,价值GydF4y2Ba论点。GydF4y2Ba名称GydF4y2Ba是参数名称和GydF4y2Ba价值GydF4y2Ba为对应值。GydF4y2Ba名称GydF4y2Ba必须出现在引号内。您可以以任何顺序指定多个名称和值对参数GydF4y2BaName1, Value1,…,的家GydF4y2Ba.GydF4y2Ba

例子:GydF4y2Ba“BinsPerOctave”“SamplingFrequency”,20日,15GydF4y2Ba

采样频率,以Hz为单位,指定为逗号分隔对,由GydF4y2Ba'采样频率'GydF4y2Ba一个正标量。GydF4y2Ba

要在CQT中使用的每个八度音阶的容器数,指定为从1到96的正整数。GydF4y2Ba

常数Q变换类型以执行,指定为包括的逗号分隔对GydF4y2Ba'transformype'GydF4y2Ba和GydF4y2Ba“全部”GydF4y2Ba或GydF4y2Ba'疏'GydF4y2Ba.稀疏变换是常q变换的最小冗余版本。GydF4y2Ba

CQT的频率限制具有对数频率响应,每个八位音程数为指定的频率箱,指定为逗号分隔对GydF4y2Ba'surformlimits'GydF4y2Ba和一个两个元素的真正的矢量。GydF4y2Ba

  • 第一个元素必须大于或等于GydF4y2BaFS / N.GydF4y2Ba, 在哪里GydF4y2BaFS.GydF4y2Ba是采样频率和GydF4y2BaNGydF4y2Ba是信号的长度。GydF4y2Ba

  • 第二个元素必须严格小于奈奎斯特频率。GydF4y2Ba

窗口作为非平稳Gabor框架的原型函数,指定为GydF4y2Ba'Hann'GydF4y2Ba那GydF4y2Ba'汉明'GydF4y2Ba那GydF4y2Ba'blackmanharris'GydF4y2Ba那GydF4y2Ba“itersine”GydF4y2Ba,或GydF4y2Ba'Bartlett'GydF4y2Ba.这些紧凑的支持函数用频率定义。金宝app对于归一化频率,它们是在区间上定义的GydF4y2Ba(1/2, 1/2)GydF4y2Ba.如果你指定一个采样频率,GydF4y2BaFS.GydF4y2Ba,它们在区间上定义GydF4y2Ba(Fs - f / 2 / 2)GydF4y2Ba.GydF4y2Ba

输出参数GydF4y2Ba

全部收缩GydF4y2Ba

恒定Q变换,返回为矩阵,多维数组,单元阵列或结构阵列。GydF4y2Ba

  • 如果GydF4y2Ba'transformype'GydF4y2Ba被指定为GydF4y2Ba“全部”GydF4y2Ba没有GydF4y2Ba'surformlimits'GydF4y2Ba那GydF4y2BaCFS.GydF4y2Ba是矩阵或多维阵列。GydF4y2Ba

    • 如果GydF4y2BaXGydF4y2Ba是向量吗GydF4y2BacqtGydF4y2Ba返回与CQT对应的矩阵。GydF4y2Ba

    • 如果GydF4y2BaXGydF4y2Ba是一个矩阵GydF4y2BacqtGydF4y2Ba获取每个列(独立频道)的CQTGydF4y2BaXGydF4y2Ba.该函数返回一个与CQT的最大冗余版本对应的多维数组。GydF4y2Ba

    阵列,GydF4y2BaCFS.GydF4y2Ba,对应于CQT的最大冗余版本。每一行的页码GydF4y2BaCFS.GydF4y2Ba对应于具有归一化中心频率(循环/样本)的通带(循环/样本)在0到1之间。归一化频率为1/2对应于奈奎斯特频率。列或跳数对应于最大的带宽中心频率,通常发生在奈奎斯特箱下方或高于奈奎斯特箱上方的一个频率箱。GydF4y2Ba

  • 如果GydF4y2Ba'transformype'GydF4y2Ba被指定为GydF4y2Ba“全部”GydF4y2Ba并且您指定频率限制,GydF4y2BaCFS.GydF4y2Ba作为具有以下四个字段的结构数组返回。GydF4y2Ba

    • CGydF4y2Ba- 用于指定频率限制内的频率的多维阵列的系数矩阵。这包括正面和“负”频率。GydF4y2Ba

    • DCCFS.GydF4y2Ba-从0到低频限的通频带的系数向量或矩阵。GydF4y2Ba

    • NyquistcfsGydF4y2Ba- 从上频率到奈奎斯特的与通带的系数矢量或矩阵。GydF4y2Ba

    • NyquistBinGydF4y2Ba- 与奈奎斯特频率对应的DFT箱。当反相CQT时使用此字段。GydF4y2Ba

  • 如果GydF4y2Ba'transformype'GydF4y2Ba被指定为GydF4y2Ba'疏'GydF4y2Ba那GydF4y2BaCFS.GydF4y2Ba是一个单元阵列,其单元数等于带通频率的数目。单元格数组的每个元素,GydF4y2BaCFS.GydF4y2Ba,是矢量或矩阵,具有等于DFT箱中带宽的值的行数,GydF4y2BaBW.GydF4y2Ba.GydF4y2Ba

CFS.GydF4y2Ba那GydF4y2BaGGydF4y2Ba,GydF4y2BafshiftsGydF4y2Ba需要输入CQT的反转GydF4y2BaICQT.GydF4y2Ba.GydF4y2Ba

对应于行的近似带通中心频率GydF4y2BaCFS.GydF4y2Ba,返回实值向量。频率从0到1排序,并在周期/样本。如果你指定GydF4y2Ba'采样频率'GydF4y2Ba, 然后GydF4y2BaFGydF4y2Ba在赫兹。GydF4y2Ba

在分析中使用的Gabor框架GydF4y2BaXGydF4y2Ba,返回作为实值向量的单元格阵列。每个矢量GydF4y2BaGGydF4y2Ba对应于一行GydF4y2BaCFS.GydF4y2Ba.GydF4y2Ba

CFS.GydF4y2Ba那GydF4y2BaGGydF4y2Ba,GydF4y2BafshiftsGydF4y2Ba需要输入CQT的反转GydF4y2BaICQT.GydF4y2Ba.GydF4y2Ba

离散傅里叶变换箱中的频率偏移,作为真实值的向量返回。班次在行的通带之间GydF4y2BaCFS.GydF4y2Ba.GydF4y2Ba

CFS.GydF4y2Ba那GydF4y2BaGGydF4y2Ba,GydF4y2BafshiftsGydF4y2Ba需要输入CQT的反转GydF4y2BaICQT.GydF4y2Ba.GydF4y2Ba

的行对应的频率间隔GydF4y2BaCFS.GydF4y2Ba,作为一个单元格阵列返回。每个元素GydF4y2BafintervalsGydF4y2Ba是一个真实值的矢量。这GydF4y2BaK.GydF4y2BaTh元素GydF4y2BafshiftsGydF4y2BaDFT箱的频移是GydF4y2Ba((k-1)mod n)GydF4y2Ba和GydF4y2Ba(k mod N)GydF4y2Ba元素GydF4y2BafintervalsGydF4y2Ba和GydF4y2Bak = 0,1,2,...,n-1GydF4y2Ba在哪里GydF4y2BaNGydF4y2Ba是频率偏移的数量。因为Matlab索引从1,GydF4y2Bafshifts (1)GydF4y2Ba包含之间的频移GydF4y2Bafintervals {结束}GydF4y2Ba和GydF4y2Bafintervals {1}GydF4y2Ba那GydF4y2Bafshifts (2)GydF4y2Ba包含之间的频移GydF4y2Bafintervals {1}GydF4y2Ba和GydF4y2Bafintervals {2}GydF4y2Ba, 等等。GydF4y2Ba

频率区间的DFT箱的带宽,GydF4y2BafintervalsGydF4y2Ba,返回实值向量。GydF4y2Ba

算法GydF4y2Ba

全部收缩GydF4y2Ba

非营养的Gabor框架GydF4y2Ba

用于频率自适应分析的非平稳Gabor (NSG)框架理论以及使用NSG框架进行分析和合成的高效算法都归功于Dörfler、Holighaus、Grill和VelascoGydF4y2Ba[1]GydF4y2Ba那GydF4y2Ba[2]GydF4y2Ba.CQT和ICQT中使用的算法由Dörfler,Holighaus,Grill和Velasco开发,并描述于GydF4y2Ba[1]GydF4y2Ba那GydF4y2Ba[2]GydF4y2Ba.在GydF4y2Ba[3]GydF4y2Ba,Schörkhuber,klapuri,holighaus和dörfler为相位校正的CQT变换开发并提供算法,该算法与天真卷积获得的CQT系数匹配。大型时频分析工具箱(GydF4y2Bahttps://github.com/ltfat.GydF4y2Ba)为非平稳Gabor框架提供了一套广泛的算法GydF4y2Ba[4]GydF4y2Ba.GydF4y2Ba

完美的重建GydF4y2Ba

通过非间断Gabor框架实现恒定Q分析的完美重建性能,GydF4y2BacqtGydF4y2Ba在内部添加零频率(DC)并将Nyquist频率附加到频率间隔。负频率是正中心频率和带宽的镜像版本GydF4y2Ba

参考文献GydF4y2Ba

[1] Holighaus,N.,M.Dörfler,G.A.Velasco和T. Grill。“可逆实时常数Q变换的框架。”GydF4y2Ba《IEEE音频、语音和语言处理汇刊》。GydF4y2Ba卷。21,第4,2013,第77.785。775-785。GydF4y2Ba

[2] Velasco, g.a., N. Holighaus, M. Dörfler,和T. Grill。用非平稳Gabor框架构造可逆常数q变换在GydF4y2Ba第14届数字音响效应国际会议的诉讼程序(DAFX-11)GydF4y2Ba.巴黎,法国:2011年。GydF4y2Ba

[3] Schörkhuber, C., A. Klapuri, N. Holighaus, M. Dörfler。一个Matlab工具箱,用于Log-Frequency分辨率的高效完美重构时频变换。提交给GydF4y2BaAES第53届语义音频国际会议GydF4y2Ba.英国伦敦:2014。GydF4y2Ba

[4] Průša, Z., P. L. Søndergaard, N. Holighaus, C. Wiesmeyr, P. Balazs。GydF4y2Ba大型时频分析工具箱2.0GydF4y2Ba.声音,音乐和运动,计算机科学讲义2014,PP 419-442。GydF4y2Ba

扩展能力GydF4y2Ba

在R2018A介绍GydF4y2Ba