www.pudn.com > zishiyingjiami.rar > RUTLDJ.m, change:2011-12-08,size:352b


function A1=RUTLDJ(A) 
[m,n]=size(A); 
A1=zeros(m,n); 
%ϵ½ 
for i=1:m 
    for j=n:-1:1 
        if i==m | j==1 
           A1(i,j)=A(i,j); 
        else 
               tp1=bitxor(A(i,j),A(i,j-1)); 
               tp2=bitxor(tp1,A(i+1,j-1)); 
               tp3=bitxor(tp2,A(i+1,j)); 
               A1(i,j)=tp3; 
        end 
    end 
end