How can I sort data like this?

I have data like this:

          A=  [6.5048 + 4.7841i
                6.5048 - 4.7841i
                1.2539 + 3.7939i
                1.2539 - 3.7939i
               -3.0305 + 3.2324i
               -3.1216 + 5.9919i
               -3.0305 - 3.2324i
               -3.1216 - 5.9919i];

how can I sort to be following:

               A=  [6.5048 + 4.7841i
                  6.5048 - 4.7841i
                  1.2539 + 3.7939i
                  1.2539 - 3.7939i
                 -3.0305 + 3.2324i
                 -3.0305 - 3.2324i
                 -3.1216 + 5.9919i
                 -3.1216 - 5.9919i];

where all data are followed directly by the conjugates?

Best regard

Ali

[Azzi Abdelmalek]

How are they sorted?

[Husni]

I am sorry it is not clear, they sorted based on this: first data is 6.5048 + 4.7841i, it then will be followed by its conjugate 6.5048 - 4.7841i which automatically becomes second data. the third data is 1.2539 + 3.7939i and will be followed by 1.2539 - 3.7939i (third & fourth data). The fifth data is -3.0305 + 3.2324i.but in the original data, the conjugate of fifth does not follow, therefore its should be reordered. I hope it is clear?

Answer 1

[Azzi Abdelmalek]

sort(A)

[Husni]

But if sort(A) the data will be:

   1.2539 - 3.7939i
   1.2539 + 3.7939i
  -3.0305 - 3.2324i
  -3.0305 + 3.2324i
  -3.1216 - 5.9919i
  -3.1216 + 5.9919i
   6.5048 - 4.7841i
   6.5048 + 4.7841i

[Azzi Abdelmalek]

Is there a problem?

[Matt Fig]

Why did you accept the answer if you are not satisfied??

[Husni]

Sorry, I try to find good way to explain. Imagine that A is complex power. I want to sort this. The top element has energy 6.5048 + 4.7841i and its complex conjugate 6.5048 - 4.7841i. the second element has energy 1.2539 + 3.7939i and its conjugate 1.2539 - 3.7939i. the third element has energy -3.0305 + 3.2324i but the following data is not its conjugate. Its conjugate locates at data number 7. so the 7 th data is moved to 5 th data (since it is the conjugate of -3.0305 + 3.2324i ) and so on. hopefully it is quite clear :-)

many thanks anyway :-)

[Matt Fig]

See my solution below.

Answer 2

[Matt Fig]

You need more than to simply call SORT. The SORT function sorts complex arrays by magnitude, which is a problem for an array like this:

A = [2+3i;
     3+2i;
     sqrt(13/2)+sqrt(13/2)*i;
     2-3i;
     3-2i;
     sqrt(13/2)-sqrt(13/2)*i]
% Note all(abs(A)==sqrt(13))

For problems like this you need to do something more involved.

[RA,J] = sort(real(A));
% Now we can put them together
newA = complex(RA,imag(A(J)))

[Husni]

Thanks matt, I got idea from your code

Answer 3

[Azzi Abdelmalek]

try this

 B=A
 while length(A)>0;
    idx=find(A==conjug(A(1)));
    B(ii)=A(1);
    B(ii+1)=A(idx);
    A(idx)=[];
    A(1)=[];
    ii=ii+2;
 end

[Husni]

Thanks Azzi for your help :-)