主要内容

cwtfreqbounds

CWT最大和最小频率或时间

描述

例子

(minfreq,maxfreq)= cwtfreqbounds (N)返回最小和最大的子波带通频率周期/样本长度的一个信号N。最小和最大频率确定为默认莫尔斯(60)小波。最低频率是决定这两个时间默认小波跨度的标准差N分粗尺度的信号。最大频率是最高频率波带通滤波器的峰值大小降到½奈奎斯特频率。

(minfreq,maxfreq)= cwtfreqbounds (N,Fs)返回赫兹的频率带通采样频率Fs

(maxperiod,minperiod)= cwtfreqbounds (N,Ts)返回带通采样周期的时间Tsmaxperiodminperiod标量的持续时间与相同的格式吗Ts。如果标准偏差的数量设置,这样log2 (maxperiod / minperiod) < 1 / NV在哪里NV是声音的数量每倍频程maxperiod调整到minperiod×2 ^ (1 / NV)

例子

(___)= cwtfreqbounds (___,名称=值)返回最小和最大的子波带通频率或时间由一个或多个指定附加选项名称=值参数。例如,[minf, maxf] = cwtfreqbounds(1000年,TimeBandwidth = 30)设置默认的莫尔斯的时间带宽参数小波为30。

例子

全部折叠

获得最小和最大的子波带通频率信号的1000个样本使用默认值。

[minfreq, maxfreq] = cwtfreqbounds (1000)
minfreq = 0.0033
maxfreq = 0.4341

获得最小和最大的子波带通频率为默认莫尔斯信号的小波长度10000和1 kHz的采样频率。设置截止至100%,最高频率波奈奎斯特频率的带通滤波器的峰值500 Hz。

sigLength = 10000;Fs = 1 e3;[minfreq, maxfreq] = cwtfreqbounds (sigLength Fs,截止= 100);

构造一个类使用值滤波器组cwtfreqbounds的回报。银行获得滤波器的频率响应。

fb = cwtfilterbank (SignalLength = sigLength SamplingFrequency = Fs,FrequencyLimits = [minfreq maxfreq]);[psidft f] = freqz (fb);

构造一个第二CWT过滤器银行相同的第一,而是使用默认频率限制。获得第二滤波器的频率响应。

fb2 = cwtfilterbank (SignalLength = sigLength SamplingFrequency = Fs);[psidft2 ~] = freqz (fb2);

对于每一个滤波器组,画出滤波器的频率响应与中心频率最高。从第一个过滤器银行确认频率响应峰值在奈奎斯特频率响应和第二滤波器组是在奈奎斯特50%的峰值大小。

:情节(f, psidft (1)):情节(f, psidft2 (1))标题(“频率响应”)包含(“频率(赫兹)”)ylabel (“级”)传说(“第一个过滤器银行”,“第二过滤器银行”,位置=“西北”)

图包含一个坐标轴对象。坐标轴对象标题频率响应,包含频率(赫兹),ylabel级包含2线类型的对象。这些对象代表第一过滤器银行,第二滤波器组。

获得的最小和最大频率撞击小波的信号长度5000和10 kHz的采样频率。指定的截断值 1 0 0 × 1 0 - - - - - - 8 / 2 的最高频率波带通滤波器衰减 1 0 - - - - - - 8 尼奎斯特。

[minf, maxf] = cwtfreqbounds (5 e3, 1 e4,小波=“撞”,截止= 100 * 1 e-8/2);

使用返回的值构造滤波器组cwtfreqbounds。情节的频率响应。

fb = cwtfilterbank (SignalLength = 5 e3,小波=“撞”,SamplingFrequency = 1 e4, FrequencyLimits = [minf maxf]);freqz(神奇动物)

图包含一个坐标轴对象。坐标轴对象标题CWT滤波器组,包含频率(赫兹),ylabel级包含80行类型的对象。

获得小波带通频率信号的最小和最大长度为4096。指定一个截止为0。显示最小和最大带通频率。

sLength = 4096;有限公司= 0;[minfreq, maxfreq] = cwtfreqbounds (sLength截止=有限公司);流(“最小频率:% f周期/样本\ nMax频率:% f周期/样本”,minfreq maxfreq)
最小频率:0.000805周期/样本最大频率:0.103574周期/样品

创建一个过滤器银行使用频率限制。获得双向小波的频率响应。

fb = cwtfilterbank (SignalLength = sLength FrequencyLimits = [minfreq maxfreq]);[psif f] = freqz (fb, FrequencyRange =“双侧”);

获得最小和最大的子波带通频率信号的长度为4096,但这一次指定的截止 1 0 0 × 1 0 - - - - - - 8 / 2 。创建第二个过滤器银行使用这些新的频率。确认第二频率范围大于第一频率范围。

有限公司= 100 * (1 e-8/2);[minfreq2, maxfreq2] = cwtfreqbounds (sLength截止=有限公司);fb2 = cwtfilterbank (SignalLength = sLength FrequencyLimits = [minfreq2 maxfreq2]);流(“最小频率:% f周期/样本\ nMax频率:% f周期/样本”,minfreq2 maxfreq2);
最小频率:0.000805周期/样本最大频率:0.281770周期/样品

获得第二的双向小波频率响应滤波器组。

[psif2, f2] = freqz (fb2 FrequencyRange =“双侧”);

画出滤波器的频率响应。

次要情节(2,1,1)情节(f, psif)标题(“频率响应:零截止滤光片银行”)ylabel (“级”)包含(“归一化频率(周期/样本)”次要情节(2,1,2)情节(f2, psif2)标题(频率响应:非零截过滤器银行”)ylabel (“级”)包含(“归一化频率(周期/样本)”)

图包含2轴对象。轴与标题频率响应对象1:零截止滤波器组,包含归一化频率(周期/样本),ylabel级包含71行类型的对象。轴与标题频率响应对象2:非零截滤波器组,包含归一化频率(周期/样本),ylabel级包含85行类型的对象。

小波滤波器的中心频率最高的在每一个滤波器组,获得在奈奎斯特频率响应的大小。观察者有最小的两个值之间的差别。

流(“零截止/级奈奎斯特:% g”psif(1楼(大小(psif, 2) / 2)))
零截止/级奈奎斯特:2.43333 e - 309
流(“非零截/级奈奎斯特:% g”psif2(1楼(大小(psif2, 2) / 2)))
非零截/级奈奎斯特:1.02265 e-08

输入参数

全部折叠

信号长度指定为一个正整数大于或等于4。

数据类型:

采样频率在赫兹,指定为一个积极的标量。

例子:[minf, maxf] = cwtfreqbounds (2048、100)

数据类型:

采样周期,指定为一个积极的标量持续时间

例子:[minp, maxp] = cwtfreqbounds(2048秒(2))

数据类型:持续时间

名称-值参数

指定可选的双参数作为Name1 = Value1,…,以=家,在那里的名字参数名称和吗价值相应的价值。名称-值参数必须出现在其他参数,但对的顺序无关紧要。

例子:[minf, maxf] = cwtfreqbounds(1000,小波=“撞”,VoicesPerOctave = 10)返回最小和最大使用凹凸小波和带通频率与1000年10每八度的声音信号样本。

R2021a之前,用逗号来分隔每一个名称和值,并附上的名字在报价。

例子:[minf, maxf] = cwtfreqbounds(1000年,“小波”、“撞”、“VoicesPerOctave”, 10)

小波分析用于确定最小和最大频率或时间,指定为“莫尔斯”,“埃莫”,或“撞”。这些字符串指定分析莫尔斯、Morlet和凹凸小波,分别。默认小波分析莫尔斯(60)小波。

莫尔斯小波,你也可以用参数表示使用的小波TimeBandwidthWaveletParameters名称-值参数。

例子:[minp, maxp] = cwtfreqbound(2048秒(1),小波=“撞”)

尼奎斯特比例的峰值大小,指定为0到100之间的一个标量。设置截止0表明小波在奈奎斯特频率响应衰减到0。设置截止One hundred.显示的值小波带通滤波器在奈奎斯特山峰。

cwtfilterbank分析小波滤波器的峰值在2的值。因此,您可以确保最高频率波衰减值α在奈奎斯特频率的设置截止100×α/ 2。在这种情况下,你必须有0≤α≤2

请注意

除非您的应用程序需要一个严格的截止值为0,考虑设置截止小非零值,例如,的108。通过指定一个较小的值,您可以增加频率范围(minfreq,maxfreq]而且还获得小波在奈奎斯特频率响应,有效地衰减到0。看到频率范围严格零和有效截止值为零

数据类型:

数量的时间标准差用来确定最低频率(最长的规模),指定为一个正整数大于或等于2。莫尔斯,分析Morlet,凹凸小波,四个标准差通常确保小波衰减为零的信号支持。金宝app递增StandardDeviations4的倍数,例如4 *,确保整个小波适合信号长度。如果标准偏差的数量设置,这样log2 (minfreq / maxfreq) > 1 / NV,在那里NV是声音的数量每倍频程minfreq调整到maxfreq×2 ^ (1 / NV)

数据类型:

莫尔斯波时间带宽,指定为一个积极的标量。莫尔斯波的对称性(γ)被认为是3。时间带宽参数越大,越分散小波在时间和更窄的小波的频率。莫尔斯的标准差小波在时间大约是√TimeBandwidth / 2)。标准差的频率大约是1/2 *√2 / TimeBandwidth

如果您指定TimeBandwidth,你就不能指定WaveletParameters

数据类型:

莫尔斯小波参数,指定为一个双元素向量。第一个元素是对称参数(γ),必须大于或等于1。第二个元素是时间带宽参数,必须大于或等于γ。时间带宽参数γ的比例不能超过40。

γ= 3时,莫尔斯在频域小波是完全对称的。偏态等于0。γ的值大于3导致正偏态,而γ的值小于3导致负偏态。

如果您指定WaveletParameters,你就不能指定TimeBandwidth

数据类型:

每个八度的声音数使用在决定必要的最小值和最大值之间的分离尺度,指定为一个整数1 - 48。最小和最大尺度是等价的最小和最大频率、最大和最小周期,分别。

数据类型:

输出参数

全部折叠

最小值小波带通频率,作为一个标量返回。minfreq在周期/样本如果SamplingFrequency没有指定。否则,minfreq在赫兹。

数据类型:

最大的小波带通频率,作为一个标量返回。maxfreq在周期/样本如果SamplingFrequency没有指定。否则,maxfreq在赫兹。

数据类型:

最大的小波带通时期,作为一个标量返回持续时间与相同的格式Ts

如果标准偏差的数量设置,这样log2 (maxperiod / minperiod) < 1 / NV,在那里NV是声音的数量每倍频程maxperiod调整到minperiod×2 ^ (1 / NV)

数据类型:持续时间

最小值小波带通时期,作为一个标量返回持续时间与相同的格式Ts

如果标准偏差的数量设置,这样log2 (maxperiod / minperiod) < 1 / NV,在那里NV是声音的数量每倍频程maxperiod调整到minperiod×2 ^ (1 / NV)

数据类型:持续时间

扩展功能

版本历史

介绍了R2018a

另请参阅

|