利用互相关对齐两个信号

11视图(30天)

显示旧的评论

露易丝·威尔逊 2021年8月12日

0
链接

这个问题直接联系

//www.tianjin-qmedu.com/matlabcentral/answers/896787-align-two-signals-using-cross-correlation

编辑: 图像分析 2021年8月12日

我有两个四通道记录仪,共有六个麦克风。由于两者之间的时间差录音机~ 5秒,我需要调整信号。我以前这样做成功的使用代码函数AlignWave文件交换,在下面的代码我对齐两个渠道第二记录器与第一频道第一个录音机。然而在我47个记录的数据集,有一个文件,这并不工作,我不知道为什么。在这种情况下,只有一两个文件的第二个录音机是正确对齐。

我很乐意分享数据,但是超过5 mb的文件是压缩格式。

                          % %对齐。wav从第二与第一记录仪记录
                         
                          为i = 1:高度(uPa_out.recorder_671105061)
                         
                          wav1 = uPa_out.recorder_671105061.pos1{我,1};%的datetime文件在第一次记录器
                         
                          wav1_dt = strsplit (uPa_out.recorder_671105061.wavFile{1},我,“。”);
                         
                          wav1_dt = datetime (wav1_dt (2),“InputFormat”,“yyMMddHHmmss”);
                         
                          start_range = wav1_dt-seconds (5);
                         
                          rec2_files = {uPa_out.recorder_335577141.pos5{1},我uPa_out.recorder_335577141.pos2{我1}};% datetime文件的第二个录音机
                         
                          rec2_files_dt = strsplit (uPa_out.recorder_335577141.wavFile{1},我,“。”);
                         
                          rec2_files_dt = datetime (rec2_files_dt (2),“InputFormat”,“yyMMddHHmmss”);
                         
                          %只有对齐频道不同记录如果有< = 5秒
                         
                          %之间
                         
                          如果(rec2_files_dt > = start_range) & (rec2_files_dt < = wav1_dt)
                         
                          为b = 1:长度(rec2_files)
                         
                          wav2 = rec2_files {1, b};
                         
                          x = xcorr (wav1 wav2);%利用互相关对齐数据
                         
                          [~ d] = max (x);% d指数最大的x的值
                         
                          延迟= d-max(长度(wav1), (wav2));
                         
                          如果延迟> 0
                         
                          yfrag =[0(延迟,1);wav2];
                         
                          其他的
                         
                          yfrag = wav2(延迟+ 1:结束);
                         
                          结束
                         
                          wav2_files_new {1, b} = yfrag;
                         
                          wav2_files_new{2,} =延迟+ 1;
                         
                          结束
                         
                          Data.five_new = wav2_files_new {1};
                         
                          Data.two_new = wav2_files_new {1,2};
                         
                          %的正确长度渠道长度都是一样的
                         
                          如果延迟> 0
                         
                          其他的%如果第二数据记录器是左移到与数据一致
                         
                          %的第一个录音机,那么你需要删除在年底“额外的”数据
                         
                          %频道第一个录音机
                         
                          长度(1,- 1)=长度(wav2_files_new {1 1});%每个新调整通道的长度会有所不同
                         
                          长度(2,1)=长度(wav2_files_new {1,2});
                         
                          min_length = min(长度);%找到最短的重塑通道,减少所有渠道的长度
                         
                          uPa_out.recorder_671105061.pos1{我1}= uPa_out.recorder_671105061.pos1{1}我(1:min_length);
                         
                          uPa_out.recorder_671105061.pos6{我1}= uPa_out.recorder_671105061.pos6{1}我(1:min_length);
                         
                          uPa_out.recorder_671105061.pos3{我1}= uPa_out.recorder_671105061.pos3{1}我(1:min_length);
                         
                          uPa_out.recorder_671105061.pos4{我1}= uPa_out.recorder_671105061.pos4{1}我(1:min_length);
                         
                          %也需要削减第二长的通道记录仪
                         
                          %就减少两个,一个将保持不变,一个将会改变
                         
                          %根据信道的时间越长
                         
                          uPa_out.recorder_335577141.pos5_new{我1}= Data.five_new (1: min_length);
                         
                          uPa_out.recorder_335577141.pos2_new{我1}= Data.two_new (1: min_length);
                         
                          结束
                         
                          其他的
                         
                          %不对齐文件如果他们有不同的开始时间
                         
                          警告(的文件从不同的录音机有不同的开始时间);
                         
                          结束
                         
                          结束
                         
                          %看看每个通道的信号
                         
                          图(我);
                         
                          次要情节(611);情节(uPa_out.recorder_671105061.pos1{1}我)
                         
                          次要情节(612);情节(uPa_out.recorder_671105061.pos6{1}我)
                         
                          次要情节(613);情节(uPa_out.recorder_671105061.pos3{1}我)
                         
                          次要情节(614);情节(uPa_out.recorder_671105061.pos4{我1});
                         
                          次要情节(615);情节(uPa_out.recorder_335577141.pos5{1}我)
                         
                          次要情节(616);情节(uPa_out.recorder_335577141.pos2{我1});
                         
                          次要情节(615);情节(Data.five_new);
                         
                          次要情节(616);情节(Data.two_new);
                         
                          次要情节(615);情节(uPa_out.recorder_335577141.pos5_new{我1});
                         
                          次要情节(616);情节(uPa_out.recorder_335577141.pos2_new{我1});