互相关的时延估计算法

function  [myk]=Tde
%[cl]=Tde（kfig）
%用互相关的方法计算time delay estimate。
%产生的两列数据为sin （t） 和 sin（t+k）。
%k为时间的延迟，fig是状态函数，为1则画图，0则不画图，默认为0。
k=input（‘input th value of k:‘）;
fig=input（‘input th value of fig（1为作图）:‘）;


t=linspace（099.999100000）;
s1=sin（t）;
s2=sin（t+k）;
if fig==1
    figure（1）;
    plot（ts1ts2‘r-.