找到一个信号测量
你得到一些数据,并想知道如果它匹配再流测量。互相关允许你做出决定,即使噪声数据。
加载到工作区桌面环锭纺纱的记录。作物一秒钟的片段和听它。
负载(“Ring.mat”)时间= 0:1 / Fs:(长度(y) 1) / Fs;m = min (y);M = max (y);Full_sig =双(y);赈灾发言人timeA = 7;timeB = 8;剪断=赈灾发言人timeA * Fs: timeB * Fs;片段= Full_sig(剪);%,类型soundsc(片段,Fs)
信号和片段的阴谋。突出了片段端点供参考。
情节(时间、Full_sig[赈灾发言人timeA timeB;赈灾发言人timeA timeB], [m m; m m],“r——”)包含(“时间(s)”)ylabel (“清洁”)轴紧
情节(剪/ Fs,片段)包含(“时间(s)”)ylabel (“清洁”)标题(“片段”)轴紧
计算和情节完整的信号的互相关和片段。
[xCorr,滞后]= xCorr (Full_sig片段);情节(滞后/ Fs, xCorr)网格包含(“滞后(s)”)ylabel (“清洁”)轴紧
的滞后互关联最大的是信号之间的时间延迟的起点。改建信号叠加的片段。
[~,我]= max (abs (xCorr));maxt =滞后(我);审判=南(大小(Full_sig));审判(maxt + 1: maxt +长度(片段))=片段;情节(时间、Full_sig时间,试验)包含(“时间(s)”)ylabel (“清洁”)轴紧
重复这个过程,但单独添加噪声信号和片段。声音不能挑出的噪音。
NoiseAmp = 0.2 *马克斯(abs(片段));片段=片段+ NoiseAmp * randn(大小(片段));Full_sig = Full_sig + NoiseAmp * randn(大小(Full_sig));%,类型soundsc(片段,Fs)情节(时间、Full_sig[赈灾发言人timeA timeB;赈灾发言人timeA timeB], [m m; m m],“r——”)包含(“时间(s)”)ylabel (“吵”)轴紧
程序发现丢失的片段,尽管高噪声水平。
[xCorr,滞后]= xCorr (Full_sig片段);情节(滞后/ Fs, xCorr)网格包含(“滞后(s)”)ylabel (“吵”)轴紧
[~,我]= max (abs (xCorr));maxt =滞后(我);审判=南(大小(Full_sig));审判(maxt + 1: maxt +长度(片段))=片段;图绘制(时间、Full_sig时间、试验)包含(“时间(s)”)ylabel (“吵”)轴紧
