ml.m

信号空间分集ssd的误码率分析

function [ output ] = ML( input )
%ML Summary of this function goes here
%   Detailed explanation goes here
r=(sqrt(2))*[input(1,1) input(2,1) input(3,1) input(4,1) input(5,1) input(6,1) input(7,1) input(8,1)];
M=[1 0 1 0 1 0 1 0;
    0 1 0 1 0 1 0 1;
    cos(pi/8) -sin(pi/8) cos((5*pi)/8) -sin((5*pi)/8) cos((9*pi)/8) -sin((9*pi)/8) cos((13*pi)/8) -sin((13*pi)/8);
    sin(pi/8) cos(pi/8) sin((5*pi)/8) cos((5*pi)/8) sin((9*pi)/8) cos((9*pi)/8) sin((13*pi)/8) cos((13*pi)/8);
    cos((2*pi)/8) -sin((2*pi)/8) cos((10*pi)/8) -sin((10*pi)/8) cos((18*pi)/8) -sin((18*pi)/8) cos((26*pi)/8) -sin((26*pi)/8);
    sin((2*pi)/8) cos((2*pi)/8) sin((10*pi)/8) cos((10*pi)/8) sin((18*pi)/8) cos((18*pi)/8) sin((26*pi)/8) cos((26*pi)/8);
    cos((3*pi)/8) -sin((3*pi)/8) cos((15*pi)/8) -sin((15*pi)/8) cos((27*pi)/8) -sin((27*pi)/8) cos((39*pi)/8) -sin((39*pi)/8);
    sin((3*pi)/8) cos((3*pi)/8) sin((15*pi)/8) cos((15*pi)/8) sin((27*pi)/8) cos((27*pi)/8) sin((39*pi)/8) cos((39*pi)/8)];% rotation M
alpha=[input(10,1) input(12,1) input(11,1) input(9,1) input(15,1) input(14,1) input(13,1) input(16,1)];
j=sqrt(-1);
u=(1/sqrt(10))*[1 -1 3 -3];
min=100;
k1=1;
k2=1;
k3=1;
k4=1;
k5=1;
k6=1;
k7=1;
k8=1;
for j1=1:4
    for j2=1:4
        for j3=1:4
            for j4=1:4
                for j5=1:4
                    for j6=1:4
                        for j7=1:4
                            for j8=1:4
                                b=[u(j1) u(j2) u(j3) u(j4) u(j5) u(j6) u(j7) u(j8)];
                                c=sum((r-b*M.*alpha).^2);
                                if c<min;
                                min=c;
                                k1=j1;k2=j2;k3=j3;k4=j4;k5=j5;k6=j6;k7=j7;k8=j8;
                                end
                            end
                        end
                    end
                end
            end
        end
    end
end
u1=[u(k1) u(k2) u(k3) u(k4) u(k5) u(k6) u(k7) u(k8)];
output=[u1(1)+j*u1(5);u1(2)+j*u1(6);u1(3)+j*u1(7);u1(4)+j*u1(8)];