nzy.m

一种用于噪声测试的程序

function [x,y]=nzy(a,b)
originimage=imread('100.bmp');    %%%%%%%%一定是方阵 
origin=originimage;
    
temp=origin;


n=length(temp) ;
m=0;
flag=0;

%%%%%%%%%%%cat映射要求  下标从0开始， 0～(n-1)

while flag==0
      for i=0:n-1                       %%%%%%%%%%%  置换一次  
           for j=0:n-1
           x=mod(i+a*j,n); 
           y=mod(b*i+(a*b+1)*j,n); 
           lastimage(x+1,y+1)=temp(i+1,j+1);
           end
      end 
      m=m+1;                           %%  置换一次完毕
         
      if m==32
         %figure,imshow(lastimage),title('40')
         a1=lastimage;
      end
         
      if(origin==lastimage)            %%%%置换多少次出现重复
            flag=1;
          %figure,imshow(lastimage),title('出现重复的图像')  
      else
         temp=lastimage;    
      end
end
imwrite(a1,'E:\MATLAB6p5p1\work\p.bmp')