function [dk]=pric(taomin,taomax,K,t)
n=2;
bq0=0.01;
tao0=(taomax-taomin)/K;
N=length(t);
for i=1:K
tao(i)=(i-0.5)*tao0+taomin;
num(i)=0;
dk(i)=0;
star(i)=0;
end
ks=0.005;
while (n<=N)
m=n-1;
while(m>=1)
ta=t(n)-t(m);
if ta<=(1-ks)*taomin
m=m-1;
elseif ta>(1+ks)*taomax
break;
else
k1=ceil((ta/(1+ks)-taomin)/tao0+1);%ceil(x) 不小于x的最小整数
k2=floor((ta/(1-ks)-taomin)/tao0+1);%floor（x） 不超过x的最大整数
for k=k1:k2
if num(k)==0
star(k)=t(n);
end
grad0=(t(n)-star(k))/tao(k);
v=floor(grad0+0.5-(10e-10));
if v==0
else
bq=grad0/v-1;
if v==1&&t(m)==star(k)
star(k)=t(n);
elseif abs(bq)<=bq0
star(k)=t(n);
end
end
sta=(t(n)-star(k))/tao(k);
dk(k)=abs(dk(k)+exp(2*pi*j*sta));
num(k)=num(k)+1;
end
end
m=m-1;
end
n=n+1;
end
plot(tao,dk)
