Halmonicratio.
谐波比例
描述
例子
计算谐波比
在音频文件中读取，使用默认参数计算谐波比，然后绘制结果。
[AudioIn，FS] = audioread（'ranswoscthree-24-96-stereo-13secs.aif'）;张视网膜=均值（AudioIn，2）;hr = houronicratio（isausimmono，fs）;t =（0：长度（isausimmono）-1）/ fs;子图（2,1,1）plot（t，audioinmono）ylabel（'振幅'）t = linspace（0，尺寸（isausimmono，1）/ fs，尺寸（hr，1））;子图（2,1,2）绘图（t，hr）xlabel（'时间''）ylabel（'谐波比'）
指定非默认参数
在音频文件中读取。
[AudioIn，FS] = audioread（'Counting-16-44p1-mono-15secs.wav'）;
使用25毫秒重叠使用50毫秒HANN Windows的音频文件的谐波比。绘制结果。
HR = Houronicratio（AudioIn，FS，......'窗户'，HANN（圆形（FS。* 0.05），'定期'），......'overlaplencth'，圆形（FS。* 0.025））;t = linspace（0，尺寸（audioin，1）/ fs，尺寸（hr，1））;绘制（t，hr）xlabel（'时间''）ylabel（'谐波比'）
谐波比指示音频谐波部分中的能量与音频的总能量的比率。因为这个例子中的音频信号有近沉默的区域，所以总能量非常低，谐波比在言论和沉默地区之间的差异辨别的差。向音频信号添加白噪声，然后计算和绘制谐波比。
AudioIn = AudioIn + 0.1 * Randn（尺寸（isaisin））;HR = Houronicratio（AudioIn，FS，......'窗户'，HANN（圆形（FS。* 0.05），'定期'），......'overlaplencth'，圆形（FS。* 0.025））;t = linspace（0，尺寸（audioin，1）/ fs，尺寸（hr，1））;绘制（t，hr）xlabel（'时间''）ylabel（'谐波比'）
计算流音频的谐波比
创建一个dsp.audiofilereader.对象在立体声音频数据逐帧中读取。创建一个dsp.signalsink.对象来记录谐波比计算。
filereader = dsp.audiofilereader（'ranswoscthree-24-96-stereo-13secs.aif'）;logger = dsp.signalsink;
在一个音频流循环中：
读入音频数据帧。
计算音频帧的每个通道的谐波比。
记录后来绘制的谐波比。
为了仅计算给定的输入框架的谐波比，请指定具有与输入相同数量的样本数量的窗口，并将重叠长度设置为零。绘制记录的数据。
win = hamping（filereader.samplesperframe，'定期'）;尽管〜ISDONE（FILEREADER）AUDION = FILEREADER（）;hr = houronicratio（isainin，filereader.samplerate，......'窗户'，赢，......'overlaplencth'，0）;记录器（HR）结尾plot（logger.buffer）ylabel（'谐波比'） 传奇（'左频道'那'右频道'）
如果对音频流循环的输入有一个可变的样本样本，则使用分析窗口大小的每帧不一致的样本Halmonicratio.，或者如果要计算重叠数据的谐波比，请使用dsp.asyncbuffer.。
创建一个dsp.asyncbuffer.对象，重置记录器，然后释放文件读取器。
buff = dsp.asyncuffer;重置（记录器）释放（Filereader）
使用25 ms重叠的50 ms帧计算谐波比。
fs = filereader.samplerate;SamplesPerframe = ROUNG（FS * 0.05）;Samplesoverlap = ROUNG（FS * 0.025）;Samplesperhop = SamplesPerframe - Samplesoverlap;WIN =汉明（样品框架）;尽管〜ISDONE（FILEREADER）AUDION = FILEREADER（）;写（buff，audioin）;尽管buff.numunReadsamples> = samplesperhop audiobuffered = read（buff，samplesperframe，samplesoverlap）;hr = houronicratio（有obiobuffered，fs，......'窗户'，赢，......'overlaplencth'，0）;记录器（HR）结尾结尾发布（Filereader）
绘制记录的数据。
plot（logger.buffer）ylabel（'谐波比'） 传奇（'左频道'那'右频道'）
色调和白噪声的谐波比
与总信号中的能量相比，谐波比测量信号的色调部分中的能量。
纯音的谐波比
创建纯音，然后使用默认参数计算谐波比。缺省情况下，谐波比率为30毫秒跳跃30毫秒的汉敏窗口。绘制结果。谐波比在1附近，这是理论最大值。
FS = 48E3;OSC =探听器（'频率'，500，......'samplesperframe'，192E3，'采样率'，fs）;sinewave = osc（）;hr = houronicratio（sinewave，fs）;t = linspace（0，尺寸（正弦波，1）/ fs，尺寸（hr，1））;绘制（t，hr）xlabel（'时间''）ylabel（'谐波比'） 标题（'sinusoid - 默认参数'）
窗口所需的短时间分析降低了从理论值的谐波比率。为了减少窗口的效果，您可以增加窗口大小。使用100毫秒汉明窗口和10 ms跳，观察谐波比更接近一个而不是使用默认窗口长度时的谐波比。
Win =汉明（圆形（FS * 0.1），'定期'）;重叠=圆形（FS * 0.099）;hr = houronicratio（sinewave，fs，'窗户'，赢，'overlaplencth'，重叠）;t = linspace（0，尺寸（正弦波，1）/ fs，尺寸（hr，1））;绘制（t，hr）xlabel（'时间''）ylabel（'谐波比'） 标题（'Sinusoid - 100 Ms窗口'）
白噪声的谐波比
创建5秒的白噪声，然后使用默认参数计算谐波比。缺省情况下，谐波比率为30毫秒跳跃30毫秒的汉敏窗口。绘制结果。谐波比为0。
FS = 48E3;噪声=兰特（FS * 5,1）;hr = houronicratio（噪音，fs）;t = linspace（0，尺寸（噪音，1）/ fs，尺寸（hr，1））;绘制（t，hr）xlabel（'时间''）ylabel（'谐波比'） 标题（'噪音 - 默认参数'）
输入参数
输出参数
算法
如上所述计算谐波比[1]。以下算法独立应用于每个音频数据窗口。信号的归一化自相关的确定为：
在哪里
S.是单一的音频数据帧N.元素。
m计算中的最大滞后。最大滞后为40毫秒，其对应于25Hz的最小频率。
在给定范围内确定谐波比的第一估计是归一化自相关的最大值：
在哪里m0.是搜索范围的下边缘，被确定为归一化自相关的第一零交叉。
最后，使用抛物线插值改善谐波比率估计，如[2]。
参考
[1] Kim，Hyoung-Gook，Nicholas Moreau和Thomas Sikora。MPEG-7音频及更大：音频内容索引和检索。John Wiley＆Sons，2005。
[2]光谱峰的二次插值。访问2018年10月11日。https://ccrma.stanford.edu/~jos/sasp/quadratic_interpolation_spectral_peaks.html
扩展能力
也可以看看
您还可以从以下列表中选择一个网站：
