main.m

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

clear;clc;
[sx,fs,nibits]=wavread('d:\aaa\qqq.wav');
x=sx(:,1);
x=x';
%x=x(1000:length(x));
%x=randn(1,40000);
t=length(x);
x1=delay(x,25,0.4);
x2=delay(x,55,0.9);
snr=60;
t1=awgn(x1,snr);%加入SNR的正态分布白噪声随机序列
t2=awgn(x2,snr);
ht1=t1(1:t);
ht2=t2(1:t);
fid=fopen('111.bin','w');
count1=fwrite(fid,ht1,'double');
fclose(fid);
fid=fopen('222.bin','w');
count2=fwrite(fid,ht2,'double');
fclose(fid);
plot(ht1,'r');
hold on;
plot(ht2,'b');
grid on