zerocrossrate
讨论二阶导数过零率
语法
描述
zerocrossrate (___）没有输出参数的情况下率沿着y-轴和相应的窗口号沿x设在。如果窗口长度等于完整的信号长度，则该函数将窗口长度沿x-轴和窗口中间的交叉速率。
例子
在信号中计算零交叉
考虑一个符号交替的向量。绘制数据图。
X = [1 -1 1 -1 1 -1 1 -1 1 -1 1 -1];情节(x)
计算的过零率x．
R =零交叉率(x)
R = 0.9500
使用第三个输出参数查找交叉发生的位置。情节x和过零点的位置。该函数返回交叉后下一个样本的索引，不一定是准确的交叉位置。第一个样本被标记为交叉点，因为函数考虑的初始状态x默认为零。
[~，~，indices] = zerocrossrate(x);情节(x)在情节(x(指标),‘*’)举行从
计算的过零率x使用比较法。该比率与使用差值法计算的值不同。
rC = zerocrossrate(x,Method=“比较”）
rC = 0.9000
计算的过零率x再次使用差分方法，并指定0为正数。速率等于使用比较方法计算的值。
rZ =零交叉率(x，零正=1)
rZ = 0.9000
现在指定的初始状态x作为1．速率等于之前的结果。
rI = 0 (x,InitialState=1)
rI = 0.9000
计算温度数据中的水平交叉点
载入2011年1月在波士顿洛根机场每小时拍摄的一组摄氏温度读数。创建一个时间表和使用调整时间将数据汇总成日平均值。
负载bostempT =小时(1:24*31)';TT =时间表(t,tempC);rTT = retime(TT，“每天”，“的意思是”）;
计算温度超过月平均值的天数。绘制数据图，并在月平均温度处画一条水平线，以显示交叉点发生的位置。
avg = mean(TT.tempC)
Avg = -1.3007
[~，count] = 0 (rTT,Level=avg)
计数= 9
plot(hours(rTT.t/24)，rTT.tempC) yline(avg) xlabel(自2011年1月1日起已经过的时间(天)) ylabel (每日平均气温(\circC))轴紧
使用零交叉来识别语音和非语音
语音可分为有声和无声两种。表达了发声，比如元音，是在声带震动时发出的。在无声的说话时，如大多数辅音，声带不振动。您可以使用零交叉来分类音频信号中的浊音和浊音区域。
加载一个音频信号到MATLAB®工作空间。那个声音说:“橡树很坚固，也能遮荫。”
[y,fs] = audioread(“oak.m4a”）;输入soundsc(y,fs)
信号在44.1千赫采样。使用比较方法计算10毫秒窗口的过零率。
Win = fs*0.01;(y,WindowLength=win,Method= .“比较”）;
情节率为了可视化每一段的交叉率。语音语音的交叉率较低，而无声语音的交叉率较高。
情节(率)
使用阈值0．1区分浊音段和浊音段。创建一个signalMask对象，有两个类别(“unvocal”和“vocal”)，并绘制感兴趣的区域(roi)。将发声和不发声的区域与每个发声单词的位置进行比较。
您可以使用Audio Toolbox™语音转文本功能从音频文件中提取单词。加载标记的信号集ls从SpeechTranscription.mat进入工作区。标记的信号集包含音频信号、ROI限制和每个口语单词的标签。有关标记音频信号的详细信息，请参见在语音信号中标注口语．在情节中展示所讲的单词。
H = 0.1;Idu = find(rate > h);Idu (1:2) = [];Vi = [(idu-1) idu]*win;M = sigroi2binmask(vi,length(y));mask = signalMask([m ~m]，Categories=[“无声”“表示”), SampleRate = fs);plotsigroi(面具,y)负载SpeechTranscriptionln = getLabelNames(ls);v = getLabelValues(ls,1,ln);v.Value = category (v.Value,v.Value);RL = v.ROILimits;VL = v.值;持有在文本(意思是(RL, 2), -0.7 *(大小(重要的)),六世,HorizontalAlignment =“中心”，.．.字形大小= 11,FontWeight =“大胆”)举行从
流数据过零率
加载一个包含15秒原声吉他音乐的音频文件。采样率为44.1 kHz。要播放音乐，取消最后一行代码的注释。
Fs = 44100;Y = audieread (“guitartune.wav”）;%的声音(y, Fs)
将音频信号缓冲为重叠帧，每个帧有4096个样本。使用512个样本的重叠。
winLength = 4096;重叠= 512;[yB，~] = buffer(y,winLength,overlap，“nodelay”）;
获得每一帧的过零率。要考虑帧重叠，请将初始状态指定为每帧中第一个样本的先前值。
init = winLength - overlap;状态= 0;ZCR = [];为i = 1:size(yB,2) zcr = [zcr;zerocrossrate(yB(:，i)，InitialState=state)];state = yB(init,i);结束
绘制音频信号并覆盖每一帧的过零率。
图yyaxis左x = 0:1/Fs:(数字(y)-1)/Fs;情节(x, y)包含(“秒”) ylabel (“振幅”) yyaxis正确的xx = (1:size(yB,2))*((winLength-overlap)/Fs);情节(xx”,zcr) ylabel (“零点交叉率”）
计算未缓冲信号的过零率。若要获得等效于缓冲信号过零率的结果，请设置WindowLength来4096而且OverlapLength来512．判断两个结果是否相等。
zcr_batch = zerocrossrate(y,WindowLength=winLength,OverlapLength=overlap);isequal (zcr_batch zcr)
ans =逻辑1
输入参数
输出参数
扩展功能
版本历史
R2021b中引入
您也可以从以下列表中选择一个网站: