zerocrossrate
讨论二阶导数过零率
语法
描述
zerocrossrate (___)没有输出参数的阴谋率沿着y设在沿着数量和相应的窗口x设在。如果窗口长度等于完整的信号,那么该函数块窗户沿的长度x设在和十字路口中间的窗口。
例子
计数信号的零交叉
考虑一个向量的交替的迹象。图数据。
x = [1 1 1 1 1 1 1 1 1 1];情节(x)
计算的零交点速度x。
r = zerocrossrate (x)
r = 0.9500
使用第三输出参数找到口岸发生的地点。情节x和零交点位置。函数返回一个索引在下次样品在穿越之后,不一定准确的交叉位置。第一个示例是标记为一个交叉点,因为函数认为的初始状态x默认为零。
[~,~,指数]= zerocrossrate (x);情节(x)在情节(x(指标),‘*’)举行从
计算的零交点速度x使用比较法。率不同于使用差分法计算值。
rC = = zerocrossrate (x)方法“比较”)
rC = 0.9000
计算的零交点速度x再次使用差分法和零指定为正数。速度等于使用比较法计算值。
rZ = zerocrossrate (x, ZeroPositive = 1)
rZ = 0.9000
现在指定的初始状态x作为1。速度等于前面的结果。
国际扶轮= zerocrossrate (x, InitialState = 1)
国际扶轮= 0.9000
计数水平过境点的温度数据
加载一组温度在摄氏度每小时在波士顿洛根机场为整个月1月,2011年。创建一个时间表和使用调整时间将数据聚合到日常的意思。
负载bostempt =小时(一24 * 31)';TT =时间表(t, tempC);rTT =调整时间(TT,“每天”,“的意思是”);
数的天数穿过月平均温度。图的数据,包括一条水平线的月平均温度可视化口岸的发生。
avg =意味着(TT.tempC)
avg = -1.3007
[~,计数]= zerocrossrate (rTT,水平= avg)
数= 9
情节(小时(rTT.t / 24), rTT.tempC) yline (avg)包含(2011年1月1日以来的时间(天))ylabel (“每日平均温度(\ circC)”)轴紧
识别语音浊音和清音的零交叉使用
语音特征可以作为有声的或无声的笑。表达了演讲中,如元音、声带振动时发生。在无声的演讲,如大多数辅音,声带不振动。您可以使用零交叉分类地区浊音和清音的音频信号。
一个音频信号加载到MATLAB®工作区。的声音说,“橡树强,也给阴影”。
(y, fs) = audioread (“oak.m4a”);%,类型soundsc (y, fs)
在44.1 kHz信号采样。计算零点交叉率为10 ms windows使用比较法。
赢得= f * 0.01;率= zerocrossrate (= y, WindowLength =赢,方法“比较”);
情节率每段交叉率形象化。言论表示预计将有一个交叉率低,而清音预计将有交叉率高。
情节(率)
使用一个阈值0.1对浊音和清音区分段。创建一个signalMask对象有两个类别(“无声”和“表达”)和情节感兴趣的区域(roi)。比较区域语音浊音和清音的每个口语的位置。
IBM®沃森演讲文本API和音频工具箱™软件可以用来提取文字从一个音频文件。负载Transcription.mat进入工作区。标记信号集包含音频信号,ROI限制,为每个口语和标签。有关详细信息,请参见标签口语单词音频信号。显示情节上的口语词汇。
h = 0.1;idu =找到(率> h);idu (1:2) = [];vi = [(idu-1)诱导器]*赢;m = sigroi2binmask (vi、长度(y));掩码= signalMask (= [[m ~ m],类别“无声”“表示”),SampleRate = fs);plotsigroi(面具,y)负载转录ln = getLabelNames (transcribedAudio);v = getLabelValues (transcribedAudio 1 ln);v。值=分类(v.Value v.Value);RL = v.ROILimits;六世= v.Value;持有在文本(意思是(RL, 2), -0.7 *(大小(重要的)),六世,HorizontalAlignment =“中心”,…字形大小= 11,FontWeight =“大胆”)举行从
讨论二阶导数过零率的流数据
加载一个原声吉他音乐的音频文件包含15秒。采样率是44.1 kHz。玩音乐,取消最后一行代码。
Fs = 44100;y = audioread (“guitartune.wav”);%的声音(y, Fs)
缓冲音频信号重叠的帧,每个国家都有4096个样本。使用一个重叠的512个样本。
winLength = 4096;重叠= 512;(yB ~) =缓冲区(y, winLength,重叠,“nodelay”);
得到每一帧的零交点。占帧重叠,指定初始状态的前一个值在每一帧首样。
init = winLength -重叠;国家= 0;zcr = [];为i = 1:尺寸(yB, 2) zcr = [zcr; zerocrossrate (yB (:, i), InitialState =状态)];我状态= yB (init);结束
情节的音频信号和覆盖每一帧的零交点。
图yyaxis左x = 0:1 / Fs:(元素个数(y) 1) / Fs;情节(x, y)包含(“秒”)ylabel (“振幅”)yyaxis正确的xx =(1:尺寸(yB, 2)) * ((winLength-overlap) / Fs);情节(xx”, zcr) ylabel (“零点交叉率”)
计算的零交点率无缓冲的信号。获得结果相当于零交点的缓冲信号,设置WindowLength来4096年和OverlapLength来512年。确定两个结果是相等的。
zcr_batch = zerocrossrate (y, WindowLength = winLength OverlapLength =重叠);isequal (zcr_batch zcr)
ans =逻辑1
