construct.m

一个基于8发1收的多天线系统的准正交差分调制例程

function y = construct( x )

%CONSTRUCT Summary of this function goes here
%  Detailed explanation goes here
d1=x(5,1)+x(5,4);d2=x(5,2)+x(5,5);d3=x(5,3)+x(5,6);d4=x(5,1)-x(5,4);d5=x(5,2)-x(5,5);d6=x(5,3)-x(5,6);
s1=0.5*x(5,7)*(x(1,1)*d1+x(1,3)*d2+x(1,4)*d3)+0.5*x(5,8)*(x(1,5)*d4+x(1,7)*d5+x(1,8)*d6);
s4=0.5*x(5,7)*(x(1,1)*d1+x(1,3)*d2+x(1,4)*d3)-0.5*x(5,8)*(x(1,5)*d4+x(1,7)*d5+x(1,8)*d6);
s2=0.5*x(5,7)*(x(3,1)*d1+x(3,3)*d2+x(3,4)*d3)+0.5*x(5,8)*(x(3,5)*d4+x(3,7)*d5+x(3,8)*d6);
s5=0.5*x(5,7)*(x(3,1)*d1+x(3,3)*d2+x(3,4)*d3)-0.5*x(5,8)*(x(3,5)*d4+x(3,7)*d5+x(3,8)*d6);
s3=0.5*x(5,7)*(x(4,1)*d1+x(4,3)*d2+x(4,4)*d3)+0.5*x(5,8)*(x(4,5)*d4+x(4,7)*d5+x(4,8)*d6);
s6=0.5*x(5,7)*(x(4,1)*d1+x(4,3)*d2+x(4,4)*d3)-0.5*x(5,8)*(x(4,5)*d4+x(4,7)*d5+x(4,8)*d6);
z=zeros(4,4);
a(1,1)=s1+s4;
a(1,2)=0;
a(1,3)=-(conj(s2)+conj(s5));
a(1,4)=-(conj(s3+s6));
a(2,1)=0;
a(2,2)=s1+s4;
a(2,3)=s3+s6;
a(2,4)=-(s2+s5);
a(3,1)=s2+s5;
a(3,2)=-(conj(s3+s6));
a(3,3)=conj(s1+s4);
a(3,4)=0;
a(4,1)=s3+s6;
a(4,2)=conj(s2)+conj(s5);
a(4,3)=0;
a(4,4)=conj(s1+s4);
b(1,1)=s1-s4;
b(1,2)=0;
b(1,3)=-(conj(s2)-conj(s5));
b(1,4)=-(conj(s3-s6));
b(2,1)=0;
b(2,2)=s1-s4;
b(2,3)=s3-s6;
b(2,4)=-(s2-s5);
b(3,1)=s2-s5;
b(3,2)=-(conj(s3-s6));
b(3,3)=conj(s1-s4);
b(3,4)=0;
b(4,1)=s3-s6;
b(4,2)=conj(s2)-conj(s5);
b(4,3)=0;
b(4,4)=conj(s1-s4);
y1=cat(2,a,z);
y2=cat(2,z,b);
y=cat(1,y1,y2);