www.pudn.com > hundunzhiluan.rar > hundunzhiluan.m

x0=abs(sin(randn)); 
xn=x0;yn=0; 
i=128.*128.*7; 
s=zeros(1,i); 
for j=1:i 
    yn=1-2.*xn^2;  
    s(j)=yn; 
    xn=yn; 
end 
c1=zeros(1,i/7); 
s1=(sign(s)+1)/2; 
% s1=round(abs(s)); 
    p=1;g=1; 
    while p