zhiluan.asv

来自「一些程序」· ASV 代码 · 共 14 行

%将x1置乱，miyao
function x=zhiluan(x1,miyao)
c=[0 0 0 0 0 1 0 1 0 0 1 1];
x=x1;
a0=miyao;
a=a0;
b=(ones(1,12)*2).^(0:11);
for i=1:4095   
    xulie(i)=sum(a.*b);
    x(i)=x1(xulie(i));
    d=rem(sum(a.*c),2);
    a(2:12)=a(1:11);
    a(1)=d;
end