demodulation.m

来自「一个基于四发二收天线系统的matlab无线传输例程」· M 代码 · 共 34 行

function y=demodulation(x)
%DEMODULATION Summary of this function goes here
%  Detailed explanation goes here
for 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(5,1))*x(6,1)+x(5,2)*ctranspose(x(6,2)))*(c1(k)+c2(n))+(ctranspose(x(7,1))*x(8,1)+x(7,2)*ctranspose(x(8,2)))*(c1(k)-c2(n)));
t(i)=r(i)+s(i);
i=i+1;
end;
end;
[m,n1]=min(t);
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);