function y=demodulation2(x)%DEMODULATION Summary of this function goes here%  Detailed explanation goes herefor i=1:4;    c1(i)=exp(j*(pi*(i-1)/2));    c2(i)=exp(j*(pi*(i-1)/2+pi/6));end;i=1;    for k=1:4;        for n=1:4;r(i)=abs(c1(k)+c2(n))^2*(abs(x(2,1))^2+abs(x(2,2))^2)/x(9,1)+abs(c1(k)-c2(n))^2*(abs(x(4,1)^2)+abs(x(4,2)^2))/x(9,2)-2*real((ctranspose(x(1,1))*x(2,2)-x(1,2)*ctranspose(x(2,1)))*(c1(k)+c2(n))+(ctranspose(x(3,1))*x(4,2)-x(3,2)*ctranspose(x(4,1)))*(c1(k)-c2(n)));s(i)=abs(c1(k)+c2(n))^2*(abs(x(6,1))^2+abs(x(6,2))^2)/x(9,1)+abs(c1(k)-c2(n))^2*(abs(x(8,1)^2)+abs(x(8,2)^2))/x(9,2)-2*real((ctranspose(x(5,1))*x(6,2)-x(5,2)*ctranspose(x(6,1)))*(c1(k)+c2(n))+(ctranspose(x(7,1))*x(8,2)-x(7,2)*ctranspose(x(8,1)))*(c1(k)-c2(n)));t(i)=r(i)+s(i);i=i+1;end;end;[m,n2]=min(t);p=mod(n2,4);if p==0 p=4;end; q=(n2-p)/4;q=q+1;y1=c1(q);y2=c2(p);y=cat(2,y1,y2);y=transpose(y);