主要内容

光谱衰减点

音频信号和听觉频谱图的频谱衰减点

描述

实例

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

实例

rolloffpoint.=频谱衰减点(x,F,名称、值)指定使用一个或多个选项的选项名称、值配对参数。

例子

全部崩溃

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

[audioIn,fs]=音频读取(“计数-16-44p1-mono-15秒波形”); rolloffPoint=频谱rolloffPoint(音频输入,fs);t=linspace(0,大小(audioIn,1)/fs,大小(rolloffPoint,1));绘图(t,滚动偏移点)xlabel(‘时间’)伊拉贝尔(“滚动点(Hz)”)

图中包含一个axes对象。axes对象包含一个line类型的对象。

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

[audioIn,fs]=音频读取(“计数-16-44p1-mono-15秒波形”);[s,cf,t]=光谱图(音频输入,fs);rolloffPoint=光谱rolloffPoint(s,cf);绘图(t,rolloffPoint)xlabel(‘时间’)伊拉贝尔(“滚动点(Hz)”)

图中包含一个axes对象。axes对象包含一个line类型的对象。

读入音频文件。

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

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

rolloffPoint=频谱rolloffPoint(音频输入,fs,...“窗口”,汉明(圆形(0.05*fs)),...“重叠长度”,圆形(0.025*fs),...“范围”,[62.5,fs/2]);t=linspace(0,大小(audioIn,1)/fs,大小(rolloffPoint,1));绘图(t,滚动偏移点)xlabel(‘时间’)伊拉贝尔(“滚动点(Hz)”)

图中包含一个axes对象。axes对象包含一个line类型的对象。

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

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); 记录器(滚动偏移点)终止plot(logger.Buffer)ylabel(“滚动点(Hz)”)

图中包含一个axes对象。axes对象包含一个line类型的对象。

使用异步缓冲器如果

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

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

  • 要计算重叠数据的光谱衰减点。

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

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

指定为重叠25毫秒的50毫秒帧计算光谱衰减点。

fs=fileReader.SampleRate;samplesPerFrame=round(fs*0.05);samplesOverlap=round(fs*0.025);samplesPerHop=samplesPerFrame-samplesOverlap;win=hamming(samplesPerFrame);虽然~isDone(fileReader)audioIn=fileReader();写入(buff,audioIn);虽然buff.numReadSamples>=samplesPerHop audioBuffered=read(buff,samplesPerFrame,samplesOverlap);rolloffPoint=spectralRolloffPoint(audioBuffered,fs,...“窗口”...“重叠长度”,0); 记录器(滚动偏移点)终止终止发布(文件阅读器)

绘制记录的数据。

plot(logger.Buffer)ylabel(“滚动点(Hz)”)

图中包含一个axes对象。axes对象包含一个line类型的对象。

输入参数

全部崩溃

输入信号,指定为向量、矩阵或三维数组。函数如何解释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.

数据类型:仅有一个的|双重的

名称值参数

指定可选的逗号分隔的字符对名称、值论据。名称是参数名和价值是对应的值。名称必须出现在引号内。您可以按任意顺序指定多个名称和值对参数,如下所示:名称1,值1,…,名称,值.

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

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

数据类型:仅有一个的|双重的

笔记

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

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

数据类型:仅有一个的|双重的

相邻窗口之间重叠的样本数,指定为逗号分隔对,由“重叠长度”和范围为[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®编码器生成C和C++代码™.

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

在R2019a中引入