%%%%+++++++++++++++++++++++++      INPUT PARAMETERSNc=127;   %%%%%%%  code lengthMc=3;     %%%%%%%   # of samples per chipNs=Nc*Mc; %%%%%%%   # of samples per information symbolIn=10;     %%%%%%%   no. of information symbols processed in one snapshot Nl=5*Mc;    %%%%%%%  max. delay spread in sample periodsK=20;                    %%%%%%%  number of interferersLs=Nl;    %%%%%%%  step in computing covariance of I+NL=5;       %%%%%%%  number of antenassensors=(0:L-1).'/2;     %%%%%%%  positions of antennassqn=sqrt(0.5);           %%%%%%%  1/2 of st.dev. of noisethetaI=[-K/2:-1 1:K/2]*2/3*pi/(K+1);   %%%%% angular positions of interferersAi=exp(i*2*pi*sensors*sin(thetaI));thetaS=[0 -15]*2*pi/180; %%%%%  angular positions of SOI (one side path)M=2;                     %%%%%  number of multipathsi=sqrt(-1);As=exp(i*2*pi*sensors*sin(thetaS));Em=zeros(M,Nl);Em(1,1)=1;  Em(M,Nl)=1;  %%%%% valid for M=2 only. Side path of SOI is delayed by Nl.%%+++++++++++++++++++++++++++++=++++ COMPUTATION OF R_M, R_L, and R_E%first_row_RL= Nc*[Mc:-1:1 zeros(1,Nl-Mc)];R_L=toeplitz(first_row_RL',first_row_RL);first_row_RM=zeros(1,Nl);k1=1:Mc; k2=Mc+1:2*Mc;first_row_RM(1:2*Mc)=Nc*[3*k1.^3-k1.^2*(6*Mc+9)+6*k1*(2*Mc+1)+4*Mc*(Mc^2-1) ...                         -(k2-2*Mc).*(k2-2*Mc-1).*(k2-2*Mc-2)]/(6*Mc);R_M=toeplitz(first_row_RM',first_row_RM);                        R_E=Em*R_L*inv(R_M)*R_L*Em';%%+++++++++++++++++++++++++++++=++++ COMPUTATION OF QRIN=Ai*Ai'+2*sqn*eye(L);Q=As'*inv(RIN)*As;%%+++++++++++++++++++++++++++++=++++ COMPUTATION OF SINR_outSINR1=real(sum(diag(Q*R_E)));SINR2=SINR1+(sum(diag(Q^2*R_E)))^2/sum(diag(Q*R_E))-sum(diag(Q^3*R_E));SINR2=real(SINR2);[SINR1 SINR2]