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-.