Divide the shuffled sequence Q={Q1，Q2，Q3，···，Qm×n}. into MN/16 matrices which are 4×4. M=256 and N=256

x=0.3;
p=0.343;
for n=2:65536;
if x(n-1)>=0 & x(n-1)<=p
 x(n)=x(n-1)/p;
else
 x(n)=(1-x(n-1))/(1-p);
end
end
A=sort(x);
[A,T]=sort(x);
Q=R(T);

whats next.? thanks in advance.

[Jan]

I've formatted your code today. Please do this by your own using the toolbar on top of the section for editing the message.

Answer 1

[Jan]

Replace

A=sort(x);
[A,T]=sort(x);
Q=R(T);

by

Q = Q(randperm(65536, 65536));

Then:

Q4 = reshape(Q, 4, 4, numel(Q)/16)

Now Q4(:, :, n) is the n.th matrix.