acousticFluctuation
波动强度的声波信号
语法
描述
波动= acousticFluctuation (audioIn,fs,calibrationFactor)
波动= acousticFluctuation (specificLoudnessIn)
(还返回特定的波动力量。波动,specificFluctuation)= acousticFluctuation (___)
(还返回占主导地位的调制频率。波动,specificFluctuation,在你)= acousticFluctuation (___)
acousticFluctuation (___)没有输出参数情节波动强度和特定的波动强度和显示文本的调制频率。如果输入是立体声,三维图显示两个通道的总和。
例子
波动的测量使用校准麦克风
建立一个实验显示的图。
创建一个audioDeviceReader从麦克风和读取一个对象audioDeviceWriter对象编写你的扬声器。
fs = 48 e3;deviceReader = audioDeviceReader (fs,“SamplesPerFrame”,2048);deviceWriter = audioDeviceWriter (fs);
创建一个audioOscillator对象来生成一个1 kHz正弦信号。
osc = audioOscillator (“正弦”1 e3,“SampleRate”fs,“SamplesPerFrame”,2048);
创建一个dsp.AsyncBuffer对象从麦克风获得缓冲数据。
大调的= 5;浅黄色= dsp.AsyncBuffer(大调的* fs);
五秒钟,正弦信号通过你的扬声器和记录使用您的麦克风。音频流的同时,注意报道的响度SPL计。一旦完成,读取缓冲区的内容对象。
numFrames =大调的* (fs / osc.SamplesPerFrame);为2 = 1:numFrames audioOut = osc ();deviceWriter (audioOut);audioIn = deviceReader ();写(浅黄色,audioIn);结束SPLreading = 60.4;micRecording =阅读(浅黄色);
计算校正因子的麦克风,使用calibrateMicrophone函数。
calibrationFactor = calibrateMicrophone (micRecording (fs + 1:,:), deviceReader.SampleRate, SPLreading);
你现在可以使用校准系数决定测量任何声音的波动通过相同的麦克风录音链。
再次执行实验,这一次,添加100%的调幅4赫兹。要创建调制信号,使用audioOscillator并指定振幅0.5和直流偏置0.5摆动在0和1之间。
mod = audioOscillator (“正弦”4“SampleRate”fs,…“振幅”,0.5,“DCOffset”,0.5,“SamplesPerFrame”,2048);大调的= 5;浅黄色= dsp.AsyncBuffer(大调的* fs);numFrames =大调的* (fs / osc.SamplesPerFrame);为2 = 1:numFrames audioOut = osc()。*国防部();deviceWriter (audioOut);audioIn = deviceReader ();写(浅黄色,audioIn);结束micRecording =阅读(浅黄色);
调用acousticFluctuation与麦克风录音,采样率和校正因子。从报告的波动acousticFluctuation使用真实的声音响度测量按照532 - 1。在5秒内显示的平均波动强度。
波动= acousticFluctuation (micRecording deviceReader.SampleRate calibrationFactor);流(的平均波动= % d (vacil) ',意味着(波动(501:最终,:)))
平均波动= 1.413824 e + 00 (vacil)
acousticFluctuation (micRecording deviceReader.SampleRate calibrationFactor)
从特定的响度测量波动
读入一个音频文件。
[audioIn, fs] = audioread (“Engine-16-44p1-stereo-20sec.wav”);
调用acousticLoudness计算具体的响度。
[~,specificLoudness] = acousticLoudness (audioIn fs,“TimeVarying”,真正的);
调用acousticSharpness没有任何输出情节声锐度。
acousticSharpness (specificLoudness“TimeVarying”,真正的)
调用acousticFluctuation没有任何输出绘制声波动。
acousticFluctuation (specificLoudness)
频率调制对声学波动的影响
生成一个纯质的音调与1500赫兹中心频率和大约700赫兹频率偏差在0.25赫兹的调制频率。
fs = 48 e3;你=0.25;大调的=
20.;numSamples = f大调的*;t = (0: numSamples-1) / fs;语气=罪(2 *π* t *作用)';fc =
1500年;excursionRatio =
0.47;游览= 2 *π* (fc * excursionRatio / fs);audioIn =调节(语气,fc、fs、“调频”,游览);
听第一个5秒的音频和绘制光谱图。
声音(audioIn (1:5 * fs), fs)谱图(audioIn (1:5 * fs),损害(512年“周期”fs), 256年,1024年,“桠溪”)
调用acousticFluctuation没有输出参数绘制声波动强度。
acousticFluctuation (audioIn fs);
指定已知调制频率
的acousticFluctuation函数允许您指定一个已知的波动频率。如果你不指定一个已知频率波动,波动的功能检测。
创建一个dsp.AudioFileReader对象读取音频信号帧。创建一个audioOscillator对象创建一个调制波。应用调制波的音频文件。
fileReader = dsp.AudioFileReader (“Engine-16-44p1-stereo-20sec.wav”);你=10.8;振幅=
0.15;osc = audioOscillator (的正弦作用,…“DCOffset”,0.5,…“振幅”振幅,…“SampleRate”fileReader.SampleRate,…“SamplesPerFrame”,fileReader.SamplesPerFrame);testSignal = [];而~结束(fileReader) x = fileReader ();testSignal = [testSignal; osc ()。* fileReader ()];结束
听两秒的测试信号波形和情节。
samplesToView = 1:2 * fileReader.SampleRate;声音(testSignal (samplesToView:), fileReader.SampleRate);情节(samplesToView / fileReader.SampleRate testSignal (samplesToView,:))包含(“时间(s)”)
绘制声波动。的检测频率调制显示文本。
acousticFluctuation (testSignal fileReader.SampleRate);
指定已知的调制频率,然后画出声波波动了。
acousticFluctuation (testSignal fileReader.SampleRate,“ModulationFrequency”,在你)
输入参数
输出参数
算法
声学波动强度是缓慢的知觉测量调节幅度或频率。描述的声音响度算法[1]和实现的acousticLoudness函数。声学波动计算中描述[2]。声学波动的算法概述如下。
在哪里f国防部是检测或调制频率和Δ吗l感知到的调制深度。如果没有指定调制频率当调用acousticFluctuation,它是自动peak-picking的频域表示声音响度。感知到的调制深度,Δl,通过修正计算具体的响度乐队通过围绕½八度过滤器f国防部,紧随其后的是一个低通滤波器来确定信封。
引用
[1]ISO 532 - 1:2017 (E)。“音响——响度的工程计算方法——第1部分:茨威格方法。”国际标准化组织。
[2]茨威格,埃伯哈德和h . Fastl。心理声学:事实和模型。2日更新,施普林格,1999年。
