zf2.rar_Zero_detector

close all; clear all; clc Q=8; % 8 psk K=2; % no of source M=6; % no of relays %%% ML detector Pbud1=[]; Symbol1=[]; IIT=10000; for iter=1:IIT Pbud=[]; Symbol=[]; % ML SER analysis for Ps= 2:2:30 Pr= Ps; zigma_sr= 0.8; zigma_rd= 0.5;zigma_sd= 0; a=Ps*zigma_rd; si= randint(1,1,Q); sj=randint(1,1,Q); b=1/4*Ps*zigma_sr*(norm(si-sj)^2); c=K*Ps*zigma_sr+1; SER =(1/12)*(((1/(b+1))+((b*c)/(a*(b+1)^2)*exp(c/(a*(b+1))*expint(c/(a*(b+1))))))^M)+(1/6)*(((3/((4*b)+3))+((12*b*c)/(a*((4*b)+3))^2)*exp((3*c)/(a*((4*b)+3)^2)*expint((3*c)/(a*((4*b)+3)))))^M); Pbud= [Pbud,Ps]; Symbol= [Symbol,SER]; end Pbud1=[Pbud1;Pbud]; Symbol1=[Symbol1;Symbol]; end semilogy(sum(Pbud1)./IIT,sum(Symbol1)./IIT,'r');hold on; Q=8; % 8-psk K=2; % no of source M=8; % no of relays % ZF SER analysis zigma_sr= 0.4; zigma_rd= 0.31;zigma_sd=0; c1= 1-(1/sqrt(Q)); c2= 3/(2*(Q-1)); c3= 2/(Q-1); Pbud= []; Symbol= []; for Ps= 2:2:30 Pr= Ps; gamma_est=(Ps*zigma_sr*Pr*zigma_rd)/(K*Ps*zigma_sr + Pr*zigma_rd+1); AZF= (1/3)*((1/(c2*gamma_est))*exp(1/(c2*gamma_est))*expint(1/(c2*gamma_est)))^(M-K+1)+(2/3)*((1/(c3*gamma_est))*exp(1/(c3*gamma_est))*expint(1/(c3*gamma_est)))^(M-K+1); SER= c1 *AZF; Pbud= [Pbud,Ps]; Symbol= [Symbol,SER]; end semilogy(Pbud,Symbol); xlabel('Power in(db)'); ylabel('SER'); title('ERROR PERFORMANCE OF ML and ZF DETECTOR'); legend( '8 psk ML','8 psk ZF'); grid on;