x=0:351
for i=1:352
 if x(i)<=31
                y(i)=1-x(i)/31;
 elseif (x(i)<=126)&(x(i)>=32)
                y(i)=[rand(1)-0.5]*0.2;
 elseif (x(i)<=159)&(x(i)>=127)
                y(i)=(x(i)-127)/32;
 elseif (x(i)<=191)&(x(i)>=160)
                y(i)=1-(x(i)-160)/32;
 elseif (x(i)<=287)&(x(i)>=192)
                y(i)=[rand(1)-0.5]*0.2;
 elseif (x(i)<=319)&(x(i)>=288)
                y(i)=(x(i)-287)/32;
 elseif (x(i)<=351)&(x(i)>=320)
                y(i)=1-(x(i)-320)/32;
         else 
               y(i)=[rand(1)-0.5]*0.2;
 end
end
for i=149:170
    y(i)=0.98
end
for i=
y=y*0.65+[rand(1,352)-0.5]*0.2;
plot(x,y)
axis([0 351 -0.5  1.2]);
xlabel('采樣點');
ylabel('循環前綴的相關值')
title('高斯白噪下的循環前綴相關特性')