主要内容

光谱衰减点

音频信号和听觉谱图的谱滚点

描述

实例

rolloffpoint.= spectralRolloffPoint (x,F)返回信号的光谱衰减点,x,随着时间的推移。函数如何解释x取决于形状F

实例

rolloffpoint.= spectralRolloffPoint (x,F,名称、值)指定使用一个或多个选项的选项名称、值对参数。

例子

全部折叠

读入一个音频文件,使用默认参数计算rolloff点,然后绘制结果。

[audioIn, fs] = audioread (“计数-16-44p1-mono-15秒波形”);rolloffPoint = spectralRolloffPoint (audioIn, fs);t = linspace(0,大小(audioIn 1) / fs,大小(rolloffPoint, 1));情节(t, rolloffPoint)包含(“时间(s)”) ylabel (“滚边点(Hz)”)

图中包含一个坐标轴。轴包含一个线型对象。

读入音频文件,然后使用光谱图作用计算mel谱图随时间的衰减点。绘制结果。

[audioIn, fs] = audioread (“计数-16-44p1-mono-15秒波形”);(年代,cf, t) = melSpectrogram (audioIn fs);rolloffPoint = spectralRolloffPoint(年代,cf);情节(t, rolloffPoint)包含(“时间(s)”) ylabel (“滚边点(Hz)”)

图中包含一个坐标轴。轴包含一个线型对象。

读入音频文件。

[audioIn, fs] = audioread (“计数-16-44p1-mono-15秒波形”);

计算功率谱随时间的衰减点。计算重叠25 ms的50 ms汉明窗数据的衰减点。使用62.5 Hz到财政司司长/2为滚动点计算。策划的结果。

rolloffPoint=频谱rolloffPoint(音频输入,fs,...“窗口”,汉明(圆形(0.05*fs)),...“重叠长度”,圆形(0.025*fs),...“范围”, 62.5, f / 2);t = linspace(0,大小(audioIn 1) / fs,大小(rolloffPoint, 1));情节(t, rolloffPoint)包含(“时间(s)”) ylabel (“滚边点(Hz)”)

图中包含一个坐标轴。轴包含一个线型对象。

创建一个音频文件读取器对象以逐帧读取音频数据。创建一个dsp。SignalSink记录光谱衰减点计算。

fileReader=dsp.AudioFileReader(“计数-16-44p1-mono-15秒波形”);记录器=dsp.SignalSink;

在音频流循环中:

  1. 读入一帧音频数据。

  2. 计算音频帧的频谱滚动点。

  3. 记录谱的滚动点,以便以后绘图。

为了仅计算给定的输入框架的频谱卷口点,请指定具有与输入相同数量的样本数量的窗口,并将重叠长度设置为零。绘制记录的数据。

win=hamming(fileReader.sampleperframe);~isDone(fileReader) audioIn = fileReader();rolloffPoint = spectralRolloffPoint (audioIn fileReader。SampleRate,...“窗口”,赢了,...“重叠长度”,0); 记录器(滚动偏移点)终止情节(logger.Buffer) ylabel (“滚边点(Hz)”)

图中包含一个坐标轴。轴包含一个线型对象。

使用异步缓冲器如果

  • 音频流循环的输入具有可变的每帧采样数。

  • 对音频流循环的输入与分析窗口的每个帧具有不一致的样本光谱衰减点

  • 你想要计算重叠数据的谱滚转点。

创建一个异步缓冲器对象,重置日志记录器,并释放文件读取器。

buff=dsp.AsyncBuffer;重置(记录器)释放(文件读取器)

指定频谱滚动点是为50毫秒帧和25毫秒重叠计算的。

fs = fileReader.SampleRate;samplesPerFrame =圆(fs * 0.05);samplesOverlap =圆(fs * 0.025);samplesPerHop = samplesPerFrame - samplesOverlap;赢得=汉明(samplesPerFrame);~isDone(fileReader)audioIn=fileReader();写入(buff,audioIn);迷。NumUnreadSamples >= samplesPerHop audioBuffered = read(buff,samplesPerFrame,samplesOverlap); rolloffPoint = spectralRolloffPoint(audioBuffered,fs,...“窗口”,赢了,...“重叠长度”,0); 记录器(滚动偏移点)终止终止发布(文件阅读器)

绘制记录的数据。

情节(logger.Buffer) ylabel (“滚边点(Hz)”)

图中包含一个坐标轴。轴包含一个线型对象。

输入参数

全部折叠

输入信号,指定为矢量、矩阵或三维阵列。函数如何解释x取决于形状F

数据类型:仅有一个的|

以Hz为单位的采样率或频率矢量,分别指定为标量或矢量。函数如何解释x取决于形状F:

  • 如果F是一个标量,x被解释为时域信号,并且F为采样率。在这种情况下,x必须是实向量或矩阵。如果x指定为矩阵时,列被解释为单个通道。

  • 如果F是一个向量,x被解释为频域信号,并且F被解释为频率,单位为Hz,对应于x.在这种情况下,x一定是真的L——- - - - - -M——- - - - - -N数组,在哪里L在给定频率的谱值的数目是多少F,M是个别光谱的数量,和N是通道数。

  • 表的行数x,L,必须等于F

数据类型:仅有一个的|

名称-值对的观点

指定可选的逗号分隔的字符对名称、值参数。名称是参数名和价值为对应值。名称必须出现在引号内。可以以任意顺序指定多个名称和值对参数Name1, Value1,…,的家

例子:“窗口”,哈明(256)

滚动点阈值,指定为逗号分隔对,由“门槛”和零和一个之间的标量,独占。

数据类型:仅有一个的|

笔记

以下名称-值对参数适用于以下情况:x是一个时域信号。如果x是频域信号,忽略名称-值对参数。

在时域中应用的窗口,指定为逗号分隔对组成“窗口”和一个实向量。向量中的元素数必须在[1,大小(x, 1)].向量中的元素数也必须大于OverlapLength

数据类型:仅有一个的|

相邻窗口之间重叠的样本数,指定为由逗号分隔的对组成“重叠长度”和范围为[0]的整数,大小(, 1)).

数据类型:仅有一个的|

用于计算窗口化输入样本DFT的容器数,指定为逗号分隔对,包括“FFTLength”一个正的标量整数。如果未指定的,FFT长度属性中的元素数量默认为

数据类型:仅有一个的|

频率范围(以Hz为单位),指定为逗号分隔对,包括“范围”一个两元素的行向量在[0,F/2].

数据类型:仅有一个的|

谱类型,指定为逗号分隔对,由“光谱类型”“权力”“级”:

  • “权力”–针对单侧功率谱计算频谱衰减点。

  • “级”–为单侧幅值谱计算谱衰减点。

数据类型:字符|一串

输出参数

全部折叠

以Hz为单位的谱滚动点,以标量、矢量或矩阵返回。每一行的rolloffpoint.对应于某窗口的光谱滚动点x. 每列rolloffpoint.对应于一个独立的通道。

算法

光谱衰减点的计算如中所述[1]:

rolloffpoint. =

以致

K = B 1. s K = κ K = B 1. B 2. s K

哪里

  • sK是bin处的光谱值K

  • B1.B2.是波段边缘,在箱子里,在上面计算光谱扩散。

  • κ之间包含的总能量的百分比B1..你可以设定κ使用阈值

参考

[1] Scheirer, E.,和M. Slaney,“鲁棒多特征语音/音乐鉴别器的构建和评估”,IEEE声学、语音和信号处理国际会议.第2卷,1997年,1221-1224页。

扩展能力

C / c++代码生成
使用MATLAB®Coder™生成C和c++代码。

GPU阵列
使用并行计算工具箱™在图形处理单元(GPU)上运行加速代码。

在R2019a中引入