t=0:99;
x=cos（4*pi*t/15）-0.8*sin（pi*t/30）;
subplot（221）;
stem（tx）;

x1=conj（x）;%共軛
r1=0;
for?n=0:99
????for?t=0:99-n
????????x2=cos（4*pi*（t+n）/15）-0.8*sin（pi*（t+n）/30）;
????????r1=r1+x1.*x2;
????end
????r2=r1/100;
end
for?n=1:99
????r3=conj（r2）;
end
r4=fliplr