function y=demodulation(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,1)+x(1,2)*ctranspose(x(2,2)))*(c1(k)+c2(n))+(ctranspose(x(3,1))*x(4,1)+x(3,2)*ctranspose(x(4,2)))*(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(1,1))*x(2,1)+x(1,2)*ctranspose(x(2,2)))*(c1(k)+c2(n))+(ctranspose(x(3,1))*x(4,1)+x(3,2)*ctranspose(x(4,2)))*(c1(k)-c2(n)));i=i+1;end;end;[m,n1]=min(r);p=mod(n1,4);if p==0 p=4;end; q=(n1-p)/4;q=q+1;y1=c1(q);y2=c2(p);y=cat(2,y1,y2);y=transpose(y);