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<i
c1(g) =s1(p).*2^6+s1(p+1).*2^5+s1(p+2).*2^4+s1(p+3).*2^3+s1(p+4).*2^2+s1(p+5).*2^1+s1(p+6).*2^0;
p=p+7;
g=g+1;
end
c1=c1+1;
c1=fix(abs(c1)) ;
x1=abs(sin(randn));
xn1=x1;yn1=0;
i1=128.*128.*7;
s1=zeros(1,i1);
for j1=1:i1
yn1=1-2.*xn1^2;
s1(j1)=yn1;
xn1=yn1;
end
c11=zeros(1,i1/7);
s11=(sign(s1)+1)/2;
> s11=round(abs(s1))
p1=1;g1=1;
while p1<i1
c11(g1) =s11(p1).*2^6+s11(p1+1).*2^5+s11(p1+2).*2^4+s11(p1+3).*2^3+s11(p1+4).*2^2+s11(p1+5).*2^1+s11(p1+6).*2^0;
p1=p1+7;
g1=g1+1;
end
c11=c11+1;
c11=fix(abs(c11)) ;
I=imread('lena.bmp');
[M,N]=size(I);
T=zeros(M,N);
k=1;l=1;
for d=1:N
for f=1:M
T(c1(k),c11(l))=I(f,d);
k=k+1;l=l+1;
end
end
imshow(T,[]);
> [m,n]=size(T);
w=zeros(M,N);
k=1;l=1;
for d=1:N
for f=1:M
w(f,d)=T(c1(k),c11(l));
k=k+1;l=l+1;
end
end
figure,imshow(w,[])