how to do this in matlab by the multiply constant matrix and the inversion matrix

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);
Q4 = reshape(Q, 4, 4, numel(Q)/16);

[Guillaume]

Personally, I'd go back to whoever wrote this assignment and ask them to rewrite it using grammatically correct sentences.

[Jan]

As I have written already in another of your questions: The line A=sort(x) is a waste oif time, because A is overwritten directly. Using the 2nd output of sort for a permutation works, but is inefficient. I've mentioned the better solution already:

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

Step 1 is solved already. Now step 2 is: "Get a new 4x4 matrix... shown if Fig 4". We cannot see Fig 4. Therefore this question cannot be answered yet and I will close it. Add the missing information to re-open it.

[Sultan Mehmood]

[Jan]

After step 1 you have a set of 4x4 matrices. Now the text of step 2 is full of syntactical errors. I'm not sure, what it means. Maybe:

Q4 = reshape(Q, 4, 4, numel(Q)/16);
R  = zeros(size(Q4));
C  = [2,3,1,1; 1,2,3,1; 1,1,2,3; 3,1,1,2];  % You should be able to write this!
for k = 1:size(Q4, 3)
    R(:, :, k) = Q4(:, :, k) * C;
end

The part with the inversion matrix is not clear. Ask the one who has written the question.