Local_comm1.rar_matlab plot equation

function [E_local_1]= Local_comm1(dm,H) Nt=1; %General parameters fc = 2.5*10^9; %Carrier frequency in Hz k = 3.5; p=1; N0 = 2*(10^(-174/10))*10^(-3);% W/Hz equivallent to -170 dBm/Hz, =N0/2 B = 10*10^(3);%Hz //B=Rs// b=12; % Optimized value for dm=1 meter M=2^b; %Power elements %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %Power of different blocks Pmix = 30.3*10^(-3);%W Psyn = 50*10^(-3);%W PLNA = 20*10^(-3);%W PIFA = 3*10^(-3);%W Pmax = 250*10^(-3);%W Pfilt = 2.5*10^(-3);% W Pfilr = 2.5*10^(-3);% W % Parameters for PDAC and PADC Beta = 1; Vdd = 3; %V Lmin = 0.5*10^(-6);%m n1 = 10; n2 = 10; fcor = 1*10^(6);% MHz I0 = 10*10^(-6);%A Cp = 1*10^(-12);%F % Calculation of PDAC Ps = (1/2)*Vdd*I0*(2^n1-1); Pd = n1*Cp*(2*B + fcor)*Vdd^2; PDAC = Beta*(Ps + Pd); % Calculation of PADC PADC = (3*(Vdd^2)*Lmin*(2*B+fcor))/10^(-0.1525*n2+4.838); %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %%%%%%%%%%%%%%%%%%%%for b=2:3 %b=2 % Optimized for coded case %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%Ton=L/B*b %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%Gd = G1.*(d.^k).*Ml %Calculation of Rbeff Rs=10*10^3; % Symbol rate in Hz Rb=b*Rs; % Bit rate in Hz v= 0.82;%m/s velocity equivallent to 60 mile/hour 5.82 c=3*10^8; lambda=c/fc; fm=v/lambda; % maximum doppler shift Tc=3/(4*fm*sqrt(pi));% Fading coherence time F=floor(Tc*Rs); Rbeff=(F-p)*Rb/F; %Calculation of total power %Total circuit power Pct=(Pmix+Pfilt+PDAC); Pcr = (Pmix+PLNA+Pfilr+PIFA+PADC); Pc=Pct+Psyn+Pcr; %Total amplifier power %Alpha calculation eta = 0.35 ;% Virtual MIMO based Cooperative communication for .... Jayaweera Page 985 zeta = 3*(M-2*sqrt(M)+1)/(M-1);% same ref (For MQAM) alpha = (zeta/eta)-1; %for Class A power amplifier %Elements for Pout calculation %%%%%%%%%%%%%%%%Ttr = 5*10^(-6)% s Ml = 10000; % 40dB %%%%%%%%%%%%%%%%%%%%T = 100 *10^(-3)%s Pb = 10^(-3); Nf = 10; % 10dB GtGr=3.16; %5dBi for Eb=10^(-21):10^(-22):10^(-15) gamma_b=(H^2/Nt)*(Eb/N0); x=sqrt(((3*b)/(M-1))*gamma_b); pb=(4/b)*(1-(1/(2^(b/2))))*q_fun(x); if pb<10^(-3) break end end %%%%%G1 = 1000 % 30dB %%%%%Gc = 3 % 4.7 dB %Calculation Pout=(((4*pi)^2)*(dm.^k)*Ml*Nf*Eb*Rb)/(GtGr*lambda^2); Pa=(1+alpha)*Pout; %Calculation of total energy Eb %Finding Ebt=(Pa+Pc)./Rbeff; E_local_1=Ebt;