main1.m

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

clear;clc;
[sx,fs,nibits]=wavread('d:\aaa\ttt.wav');
x=sx(:,1);
x=x';
x1=delay(x,2,0.3);
x2=delay(x,2,0.4);
snr=60;
t1=awgn(x1,snr);%加入SNR的正态分布白噪声随机序列
t2=awgn(x2,snr);
fs=16000;
M1 = blockFrames(x1, fs, 128, 256);
M2 = blockFrames(x2, fs, 128, 256);
for i