clear?all
clc
N=1000;
t1=0:333;
t2=0.1:sqrt（2）:（0.1+332*sqrt（2））;
t3=0.2:sqrt（5）:（0.2+332*sqrt（5））;
t=[t1?t2?t3];
clear?t1?t2?t3
t=sort（t）;
K=201;
taumin=0;
taumax=10;
b=（taumax-taumin）/K;
D=zeros（1K）;
for?i=1:K
????tauk（i）=（i-1/2）*（taumax-taumin）/K+taumin;
end
n=2;
while?n<=N
????m=n-1;
????while?m>=1
????????tau=t（n）-t（m）;
????????if?（tau>taumin）&（tau<=taumax）
????????????for?k=1:K
????????????????if?（tau>（tauk（k）-b/2））&（tau<=（tauk（k）+b/2））
????????????????????D（k）=D（k）+1;????%exp（2*pi*t（n）*j/tau）
????????????????end
????????????end
????????elseif?tau>taumax
????????????break
????????else?
????????????;
????????end
????????m=m-1;
????end
????n=n+1;
end
plot（taukD）
axis（[0?10?0?800]）

?屬性????????????大小?????日期????時(shí)間???名稱
-----------?---------??----------?-----??----
?????目錄???????????0??2011-09-09?20:13??PRI_transform\
?????文件?????????696??2010-04-13?20:01??PRI_transform\autocorrelation_function.m
?????文件????????7339??2010-05-12?19:01??PRI_transform\data_generate.m
?????文件????????1549??2010-04-05?23:09??PRI_transform\improved2.m
?????文件????????2700??2010-04-19?12:51??PRI_transform\improved_thresholding.m
?????文件?????????960??2010-04-17?15:14??PRI_transform\pritransform.m
?????文件????????1280??2010-04-27?20:53??PRI_transform\pritransform2.m