mrrc1xnr_cor.m

MIMO多天线收发 各种性能比较（BER曲线） 包括MMRC EGC SC Alamouti V-BLAST。信道为瑞利信道

% M file for MRRC 1xNr and Alamouti 2xNr

function ber_m1_MRRC12=MRRC1xNr_Cor(snr,Nr,p)

Nt=1;
Nr=Nr;
p=p;
no_bit_sym=1;
no_it_x_SNR=10000;
%iter=0;
tot_err_m1=0;
% while tot_err_m1<soglia
%     iter=iter+1;
    for i=1:no_it_x_SNR
        %DATA
        X=(2*round(rand(Nt,1))-1);
        %Channel
        H=Cor_H(Nr,Nt,p);
        %Noise
        sig=sqrt(0.5/(10^(snr/10)));
        n=sig*(randn(Nr,Nt)+j*randn(Nr,Nt));
        %Received Signal
        R=H*X+n;
        %Combiner
        s0=0;
        for n=1:Nr
            s0=conj(H(n,1))*R(n,1)+s0;
        end
        if(real(s0)>0)
            decoded=1;
        else
            decoded=-1;
        end
%         S=[s0];
%         dh=[1 -1];
%         d11=((dh(1)-real(S(1)))^2+(imag(S(1)))^2);
%         d12=((dh(2)-real(S(1)))^2+(imag(S(1)))^2);
%         D1=[d11 d12];
%         for k=1:2
%             X1_dec(k)=((abs(dh(k)))^2)*sum(sum((abs(H)).^2)-1)+D1(k);
%         end
%         
%         Decision
%         [sceltal,posizione1]=min(X1_dec);
%         decoded=[dh(posizione1)];
         err_m1=sum(round(X')~=round(decoded));
         tot_err_m1=err_m1+tot_err_m1;
         end
         %end
ber_m1_MRRC12=tot_err_m1/(no_it_x_SNR)