rx_diversity_proc2.m

MIMO-OFDM 信道的自适应仿真程序

function output = rx_diversity_proc2(rx_signal,Modulation,H,ll,s, K)
%rx_signal
r1=rx_signal(1,1:ll);
r2=rx_signal(2,1:ll);
r3=rx_signal(3,1:ll);
r4=rx_signal(4,1:ll);

H=H.';

H_sum=0;

for i=1:4
    for j=1:4
        H_sum=H_sum+abs(H(i,j))^2;
    end
end

for  i=1:8:8*floor(ll/8);
    for j=1:K
        b11=r1(i)*conj(H(1,1))+r1(i+1)*conj(H(2,1))+r1(i+2)*conj(H(3,1))+r1(i+3)*conj(H(4,1))+conj(r1(i+4))*H(1,1)+conj(r1(i+5))*H(2,1)+conj(r1(i+6))*H(3,1)+conj(r1(i+7))*H(4,1);
        b12=r2(i)*conj(H(1,2))+r2(i+1)*conj(H(2,2))+r2(i+2)*conj(H(3,2))+r2(i+3)*conj(H(4,2))+conj(r2(i+4))*H(1,2)+conj(r2(i+5))*H(2,2)+conj(r2(i+6))*H(3,2)+conj(r2(i+7))*H(4,2);
        b13=r3(i)*conj(H(1,3))+r3(i+1)*conj(H(2,3))+r3(i+2)*conj(H(3,3))+r3(i+3)*conj(H(4,3))+conj(r3(i+4))*H(1,3)+conj(r3(i+5))*H(2,3)+conj(r3(i+6))*H(3,3)+conj(r3(i+7))*H(4,3);
        b14=r4(i)*conj(H(1,4))+r4(i+1)*conj(H(2,4))+r4(i+2)*conj(H(3,4))+r4(i+3)*conj(H(4,4))+conj(r4(i+4))*H(1,4)+conj(r4(i+5))*H(2,4)+conj(r4(i+6))*H(3,4)+conj(r4(i+7))*H(4,4);
        b1_temp=b11+b12+b13+b14;
        b1(j)=abs(b1_temp-s(j))^2+abs(2*H_sum-1)*abs(s(j))^2;
        
        
        b21=r1(i)*conj(H(2,1))-r1(i+1)*conj(H(1,1))-r1(i+2)*conj(H(4,1))+r1(i+3)*conj(H(3,1))+conj(r1(i+4))*H(2,1)-conj(r1(i+5))*H(1,1)-conj(r1(i+6))*H(4,1)+conj(r1(i+7))*H(3,1);
        b22=r2(i)*conj(H(2,2))-r2(i+1)*conj(H(1,2))-r2(i+2)*conj(H(4,2))+r2(i+3)*conj(H(3,2))+conj(r2(i+4))*H(2,2)-conj(r2(i+5))*H(1,2)-conj(r2(i+6))*H(4,2)+conj(r2(i+7))*H(3,2);
        b23=r3(i)*conj(H(2,3))-r3(i+1)*conj(H(1,3))-r3(i+2)*conj(H(4,3))+r3(i+3)*conj(H(3,3))+conj(r3(i+4))*H(2,3)-conj(r3(i+5))*H(1,3)-conj(r3(i+6))*H(4,3)+conj(r3(i+7))*H(3,3);
        b24=r4(i)*conj(H(2,4))-r4(i+1)*conj(H(1,4))-r4(i+2)*conj(H(4,4))+r4(i+3)*conj(H(3,4))+conj(r4(i+4))*H(2,4)-conj(r4(i+5))*H(1,4)-conj(r4(i+6))*H(4,4)+conj(r4(i+7))*H(3,4);
        b2_temp=b21+b22+b23+b24;
        b2(j)=abs(b2_temp-s(j))^2+abs(2*H_sum-1)*abs(s(j))^2;
        
        b31=r1(i)*conj(H(3,1))+r1(i+1)*conj(H(4,1))-r1(i+2)*conj(H(1,1))-r1(i+3)*conj(H(2,1))+conj(r1(i+4))*H(3,1)+conj(r1(i+5))*H(4,1)-conj(r1(i+6))*H(1,1)-conj(r1(i+7))*H(2,1);
        b32=r2(i)*conj(H(3,2))+r2(i+1)*conj(H(4,2))-r2(i+2)*conj(H(1,2))-r2(i+3)*conj(H(2,2))+conj(r2(i+4))*H(3,2)+conj(r2(i+5))*H(4,2)-conj(r2(i+6))*H(1,2)-conj(r2(i+7))*H(2,2);
        b33=r3(i)*conj(H(3,3))+r3(i+1)*conj(H(4,3))-r3(i+2)*conj(H(1,3))-r3(i+3)*conj(H(2,3))+conj(r3(i+4))*H(3,3)+conj(r3(i+5))*H(4,3)-conj(r3(i+6))*H(1,3)-conj(r3(i+7))*H(2,3);
        b34=r4(i)*conj(H(3,4))+r4(i+1)*conj(H(4,4))-r4(i+2)*conj(H(1,4))-r4(i+3)*conj(H(2,4))+conj(r4(i+4))*H(3,4)+conj(r4(i+5))*H(4,4)-conj(r4(i+6))*H(1,4)-conj(r4(i+7))*H(2,4);
        b3_temp=b31+b32+b33+b34;
        b3(j)=abs(b3_temp-s(j))^2+abs(2*H_sum-1)*abs(s(j))^2;

        b41=r1(i)*conj(H(4,1))-r1(i+1)*conj(H(3,1))+r1(i+2)*conj(H(2,1))-r1(i+3)*conj(H(1,1))+conj(r1(i+4))*H(4,1)-conj(r1(i+5))*H(3,1)+conj(r1(i+6))*H(2,1)-conj(r1(i+7))*H(1,1);
        b42=r2(i)*conj(H(4,2))-r2(i+1)*conj(H(3,2))+r2(i+2)*conj(H(2,2))-r2(i+3)*conj(H(1,2))+conj(r2(i+4))*H(4,2)-conj(r2(i+5))*H(3,2)+conj(r2(i+6))*H(2,2)-conj(r2(i+7))*H(1,2);
        b43=r3(i)*conj(H(4,3))-r3(i+1)*conj(H(3,3))+r3(i+2)*conj(H(2,3))-r3(i+3)*conj(H(1,3))+conj(r3(i+4))*H(4,3)-conj(r3(i+5))*H(3,3)+conj(r3(i+6))*H(2,3)-conj(r3(i+7))*H(1,3);
        b44=r4(i)*conj(H(4,4))-r4(i+1)*conj(H(3,4))+r4(i+2)*conj(H(2,4))-r4(i+3)*conj(H(1,4))+conj(r4(i+4))*H(4,4)-conj(r4(i+5))*H(3,4)+conj(r4(i+6))*H(2,4)-conj(r4(i+7))*H(1,4);
        b4_temp=b41+b42+b43+b44;
        b4(j)=abs(b4_temp-s(j))^2+abs(2*H_sum-1)*abs(s(j))^2;
    end
    
    [x1,o1]=min(b1);
    [y1,p1]=min(x1);
    
    [x2,o2]=min(b2);
    [y2,p2]=min(x1);
  
    [x3,o3]=min(b3);
    [y3,p3]=min(x1);

    [x4,o4]=min(b4);
    [y4,p4]=min(x1);
    
    out(i:i+7)=[s(o1) s(p1) s(o2) s(p2) s(o3) s(p3) s(o4) s(p4)];

end

    output=out(1:2:8*floor(length(r1)/8));
        
    clear b11 b12 b13 b14 b21 b22 b23 b24
    clear b31 b32 b33 b34 b41 b42 b43 b44
    clear b1 b2 b3 b4 x1 x2 x3 x4 y1 y2 y3 y4
    clear o1 o2 o3 o4 p1 p2 p3 p4