main.m

整数时延和分数时延的方法。分数时延是用的基于拉格朗日插值的FIR滤波！然后256个点读取数据 进行GCCPHAT的算法

clear;clc;
[sx,fs,nibits]=wavread('d:\aaa\ttt.wav');
x=sx(:,1);
x=x';
t=length(x);
%x=randn(1,40000);
%t=length(x);
x1=delay(x,3,0.4);
x2=delay(x,20,0.9);
snr=20;
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);
grid on