matlab 广义互相关,广义互相关求时延的问题

clear all

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%获得语音文件%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

[x,fs,nbits]=wavread('bluesky3.wav');     %%%%%%打开wav文件

%fs=8000

%nbits=16

x=x/max(abs(x));                     %%%%%幅度归一化

N=length(x);                         %%%%%设x的采样长度为wav_length

tt=[0:N-1]/fs;                       %%%%%采样时间序列tt

wlen=200;                                  %%%%%%%帧长200，帧移80

inc=80;

fn=fix((N-wlen)/inc)+1;       %%%%%%%求出一共分了多少帧

frametime=frame2time(fn,wlen,inc,fs);     %%%%%%%%计算每帧对应的时间

x1=x';                               %%%%%构造延时800个采样点，即0.1S的信号

x2=[zeros(1,800),x1];

x2=x2(1,[1:end-800]); %%%%为保持x2与x1有相同的采样长度，将x2的最后800个采样点删去

figure(1)

subplot(211)

plot(tt,x1);

subplot(212)

plot(tt,x2);