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

clear all

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

[x,fs,nbits]=wavread(luesky3.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);