timedelay 时延估计源代码

时延估计源代码，互相关函数和自适应滤波法，可用于声源定位

clear all;              
close all;
N=1024;                   
k=128;                    
u=0.0001;
n=1:N;
a=5;
s1=a*sin（0.05*pi*n）;     
s2=a*sin（0.05*pi*（n-25））;
s1=awgn（s120‘measured‘）;
s2=awgn（s220‘measured‘）;
figure（1）;
subplot（211）
plot（nreal（s1））;        %   信号s1时域波形
title（‘信号s1时域波形‘）;
xlabel（‘n‘）;
ylabel（‘s1‘）;
subplot（212）;
plot（nreal（s2））;
title（‘信号s2时域波形‘）; % 信号s2时域波形
xlabel（‘n‘）;
ylabel（‘s2‘）;
% 设置初值
y=zeros（N1）;            %   输出信号y
y（1:k）=s1（1:k）;          %   将输入信号xn的前k个值作为输出y的前k个值
w=zeros（1N）;            %   设置抽头加权初值
e=zeros（N1）;            %   误差信号
x=zeros（k1）;
% 用LMS算法迭代滤波
for i=k+1:N ;
x=s1（i-k:i-1）;
y（i）=w（11:k）*x‘ ;   
e（i）=s2（i）-y（i）;
w=w+u*e（i）*s1;
end
[CI]=max（w）;
I=mod（I40）
figure（2）;
subplot（211）;
plot（e）
title（

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     文件         716  2014-03-27 17:06  时延估计\Untitled2.m
     文件         689  2014-04-03 17:22  时延估计\Untitled3.m
     文件         935  2014-03-27 16:56  时延估计\final.m
     文件         912  2013-12-16 09:58  时延估计\gcc.m
     目录           0  2014-04-02 14:47  时延估计\