untitled.m

本次的算法是相位周期展开聚类

       subplot(2,1,1)
        plot(s2,'r');
        hold on;
        plot(s4,'b');
        title('两段16K采样的音乐信号');
        legend('信号x2','信号x1');
        xlabel('t');
        ylabel('y');
        
        subplot(2,1,2)
        x=xcorr(s2,s4);
        x1=xcorr(s2,s2);
        x2=xcorr(s4,s4);
        x11=fft(x1);
        x22=fft(x2);
        y=fft(x);
        r=abs(y)*abs(y)./(x11.*x22);
        q=r./(abs(y).*(1-r));
        y=y./q;
        R=ifft(y);
        plot(R);
        title('GCC-PHAT SNR=10');
        xlabel('t');
        ylabel('y');
        grid on