主要内容

spectralSlope

光谱斜率为音频信号和听觉谱图

自从R2019a

所有的页面崩溃

语法

斜率= spectralSlope (x, f)
斜率= spectralSlope (x, f, Name =值)
spectralSlope (___)

描述

例子

坡= spectralSlope (x,f)返回信号的光谱斜率,x随着时间的推移。函数如何解释x取决于的形状f。

例子

坡= spectralSlope (x,f,名称=值)使用一个或多个名称参数指定选项。

例子

spectralSlope (___)没有输出参数的光谱斜率的阴谋。您可以指定一个输入组合的任何以前的语法。

  • 如果输入是在时域,光谱斜率是绘制与时间。

  • 如果输入是在频域,光谱斜率是绘制帧数。

例子

全部折叠

打开生活的脚本

读入一个音频文件,计算边坡使用默认参数。

[audioIn, fs] = audioread (“Counting-16-44p1-mono-15secs.wav”);斜率= spectralSlope (audioIn, fs);

绘制光谱斜率与时间。

spectralSlope (audioIn fs)

图包含一个坐标轴对象。坐标轴对象包含时间(s), ylabel斜率包含一个类型的对象。

打开生活的脚本

读入一个音频文件,然后计算梅尔声谱图使用melSpectrogram函数。随着时间的推移计算梅尔声谱图的斜率。策划的结果。

[audioIn, fs] = audioread (“Counting-16-44p1-mono-15secs.wav”);(年代,cf, t) = melSpectrogram (audioIn fs);斜率= spectralSlope(年代,cf);情节(t,斜率)包含(“时间(s)”)ylabel (“坡”)

图包含一个坐标轴对象。坐标轴对象包含时间(s), ylabel斜率包含一个类型的对象。

打开生活的脚本

读入一个音频文件。

[audioIn, fs] = audioread (“Counting-16-44p1-mono-15secs.wav”);

计算级频谱随时间变化的斜率。计算斜率为50毫秒汉明窗的数据与25 ms重叠。使用范围从62.5赫兹fs/ 2的斜率计算。

斜率= spectralSlope (audioIn fs,…窗口=汉明(圆(0.05 * fs)),…OverlapLength = (0.025 * fs),…范围= [62.5,f / 2]);

绘制光谱斜率与时间。

spectralSlope (audioIn fs,…窗口=汉明(圆(0.05 * fs)),…OverlapLength = (0.025 * fs),…范围= [62.5,f / 2]);

图包含一个坐标轴对象。坐标轴对象包含时间(s), ylabel斜率包含一个类型的对象。

打开生活的脚本

创建一个dsp.AudioFileReader对象读取音频数据帧。创建一个dsp.SignalSink记录光谱斜率计算。

fileReader = dsp.AudioFileReader (“Counting-16-44p1-mono-15secs.wav”);记录器= dsp.SignalSink;

在一个音频流循环:

  1. 读取音频数据的帧。

  2. 计算光谱斜率的音频帧。

  3. 记录光谱斜率为以后策划。

计算光谱斜率只有一个给定的输入框,指定一个窗口与相同数量的样本作为输入,并将重叠长度设置为零。情节记录的数据。

赢得=汉明(fileReader.SamplesPerFrame);而~结束(fileReader) audioIn = fileReader ();斜率= spectralSlope (audioIn fileReader.SampleRate,…“窗口”,赢了,…“OverlapLength”,0);记录器(坡)结束情节(logger.Buffer) ylabel (“坡”)

图包含一个坐标轴对象。ylabel斜率的坐标轴对象包含一个类型的对象。

使用dsp.AsyncBuffer如果

  • 输入你的音频流循环变量samples-per-frame。

  • 输入你的音频流循环不一致samples-per-frame分析窗口spectralSlope。

  • 你想计算重叠光谱斜率数据。

创建一个dsp.AsyncBuffer对象,重置记录器和发布的文件阅读器。

浅黄色= dsp.AsyncBuffer;重置(logger)释放(fileReader)

为50毫秒指定的光谱斜率计算帧25 ms重叠。

fs = fileReader.SampleRate;samplesPerFrame =圆(fs * 0.05);samplesOverlap =圆(fs * 0.025);samplesPerHop = samplesPerFrame - samplesOverlap;赢得=汉明(samplesPerFrame);而~结束(fileReader) audioIn = fileReader ();写(浅黄色,audioIn);而迷。NumUnreadSamples >= samplesPerHop audioBuffered = read(buff,samplesPerFrame,samplesOverlap); slope = spectralSlope(audioBuffered,fs,…“窗口”,赢了,…“OverlapLength”,0);记录器(坡)结束结束发行版(fileReader)

情节记录的数据。

情节(logger.Buffer) ylabel (“坡”)

图包含一个坐标轴对象。ylabel斜率的坐标轴对象包含一个类型的对象。

输入参数

全部折叠

输入信号,指定为一个向量,矩阵,或三维数组。函数如何解释x取决于的形状f。

数据类型:单|双

采样率和频率在赫兹矢量,分别指定为一个标量或矢量。函数如何解释x取决于的形状f:

  • 如果f是一个标量,x是解释为一个时域信号,f解释为采样率。在这种情况下,x必须是一个真正的向量或矩阵。如果x被指定为一个矩阵,列是解释为个人渠道。

  • 如果f是一个矢量,x是解释为频域信号,f被解释为频率,在赫兹,对应的行吗x。在这种情况下,x必须是一个真正的l——- - - - - -米——- - - - - -N数组,l是在给定频率的谱值的数量吗f,米是个体数量的光谱,N通道的数量。

  • 的行数x,l,必须等于要素的数量f。

数据类型:单|双

名称-值参数

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

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

例子:窗口=汉明(256)

请注意

适用以下名称参数x是一种时域信号。如果x是一个频域信号,名称参数被忽略。

应用于时域窗口,指定为一个真正的向量。向量中的元素的数量必须在[1,大小(x,1)]。向量中的元素的数量也必须大于OverlapLength。

数据类型:单|双

数量的样本之间的重叠相邻窗口,指定为一个整数范围在[0,大小(窗口,1))。

数据类型:单|双

使用的箱子数量计算DFT的窗口的输入样本,指定为一个积极的标量整数。如果未指定的,FFTLength默认为元素的数量窗口。

数据类型:单|双

在赫兹频率范围,指定为一个双元素的行向量增加实际价值范围在[0,f/ 2)。

数据类型:单|双

光谱类型,指定为“权力”或“级”:

  • “权力”——光谱斜率计算单边功率谱。

  • “级”——片面的光谱斜率计算级频谱。

数据类型:字符|字符串

输出参数

全部折叠

光谱斜率在赫兹,作为一个标量,返回向量或矩阵。每一行的坡对应于光谱斜率的窗口x。每一列的坡对应于一个独立的通道。

算法

中描述的光谱斜率计算[1]:

坡 = ∑ k = b 1 b 2 ( f k − μ f ) ( 年代 k − μ 年代 ) ∑ k = b 1 b 2 ( f k − μ f ) 2

在哪里

  • fk是赫兹的频率对应于本吗k。

  • μf是平均频率。

  • 年代k在本是光谱值吗k。

  • μ年代是均值谱值。

  • b1和b2乐队边缘,在垃圾箱,计算光谱斜率。

引用

[1]我们,亚历山大。介绍音频内容分析应用于信号处理和音乐信息。皮斯卡塔韦,新泽西:IEEE出版社,2012年。

扩展功能

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

GPU数组
加速代码运行在一个图形处理单元(GPU)使用并行计算工具箱™。

版本历史

介绍了R2019a

另请参阅

|

主题