arnoldperiod.m

来自「图像置乱」· M 代码 · 共 27 行

function Period=ArnoldPeriod(N)% 求周期

if ( N<2 )
    Period=0;
    return;
end

%初始位置

n=1;
x=1;
y=1;

%通过循环寻找周期
while (n~=0)
    xn=x+y;
    yn=x+2*y;
    
    if ( mod(xn,N)==1 & mod(yn,N)==1 )%再次的回到原来的位置，完成一次的周期
        Period=n;
        return;
    end
    
    x=mod(xn,N);
    y=mod(yn,N);
    n=n+1;
end