EAST.rar_east

clc clear all tic %AN ADAPTIVE ENERGY EFFICIENT TRANSMISSION SCHEME USING MOBILE SINK FOR WIRELESS SENSOR NETWORKS. %Let suppose we have 100 nodes wireless sensor network(100*100)m. %CALCULATING ETA ETA=(250000/83500000); %CALCULATING RECEIVER SNR SNR=ETA*6.7608; %We have randomly selected the temperature in range from [-10(C)- 53(C)]. T=zeros(100,1); Nnodes=100; T=ceil(-10+63.*rand(Nnodes,1)); %Then to estimate RSSI-loss(dBm) for given temperature. for j=1:1:100 Rssi(j,1)=0.1996*((T(j,1))-25); end % To estimate required transmitter power level for given RSSI. for k=1:1:100 Plevel(k,1)=power(((Rssi(k,1)+40)/12),2.91); end %Division of network into three regions on basis of RSSI-loss. %Finding minimum value in array. RssiL=min(Rssi); %Finding maximum value in array. RssiH=max(Rssi); %Finding average value in array. RssiA=(RssiL+RssiH)/2; %Finding the number of nodes in A region(RssiH). RssiLimitH=RssiA+2; countH=0; A=zeros(100,1); for l=1:1:100 if (Rssi(l,1)<=RssiH) if (Rssi(l,1)>RssiLimitH) countH=countH+1; A(l,1)=Rssi(l,1); end end end %Finding the number of nodes in B region(RssiM). RssiLimitMH=RssiA+2; RssiLimitML=RssiA-2; countM=0; B=zeros(100,1); for m=1:1:100 if (Rssi(m,1)<=RssiLimitMH) if (Rssi(m,1)>RssiLimitML) countM=countM+1; B(m,1)=Rssi(m,1); end end end %Finding the number of nodes in C region(RssiL). RssiLimitLH=RssiA-2; countL=0; C=zeros(100,1); for n=1:1:100 if (Rssi(n,1)<=RssiLimitLH) if (Rssi(n,1)>=RssiL) countL=countL+1; C(n,1)=Rssi(n,1); end end end % Estimating RSSI and Plevel for each node in region A. count1=1; for o=1:1:length(A) if (A(o,1)~=0) RssiAAA(count1,1)=A(o,1); PlevelAAA(count1,1)=power(((A(o,1)+40)/12),2.91); count1=count1+1; end end % Estimating RSSI and Plevel for each node in region B. count2=1; for p=1:1:length(B) if (B(p,1)~=0) RssiBBB(count2,1)=B(p,1); PlevelBBB(count2,1)=power(((B(p,1)+40)/12),2.91); count2=count2+1; end end % Estimating RSSI and Plevel for each node in region C. count3=1; for q=1:1:length(C) if (C(q,1)~=0) RssiCCC(count3,1)=C(q,1); PlevelCCC(count3,1)=power(((C(q,1)+40)/12),2.91); count3=count3+1; end end %Plotting mean RSSI (Threshold) for each region. RssiAAAmean=mean(RssiAAA(1:count1-1)); RssiBBBmean=mean(RssiBBB(1:count2-1)); RssiCCCmean=mean(RssiCCC(1:count3-1)); %Now we will estimate corresponding power level. PlevelA=power(((RssiAAAmean+40)/12),2.91); PlevelB=power(((RssiBBBmean+40)/12),2.91); PlevelC=power(((RssiCCCmean+40)/12),2.91); %Defining Rounds. for rmax=1:1:100 xm=100; ym=100; n=100; c=n/10; for i=1:1:n S(i).xd=xm*rand(1,1); S(i).yd=ym*rand(1,1); end for b=1:1:n if(b<c) S(n+1).xd=50; S(n+1).yd=50; for i=1:1:n dis(i)=sqrt((S(i).xd-S(n+1).xd)^(2)+(S(i).yd-S(n+1).yd)^(2)); end end if(b>c && b<(2*c)) S1(n+1).xd=0; S1(n+1).yd=0; for i=1:1:n dis1(i)=sqrt((S(i).xd-S1(n+1).xd)^(2)+(S(i).yd-S1(n+1).yd)^(2)); end end if(b>(2*c) && b<(3*c)) S2(n+1).xd=0; S2(n+1).yd=100; for i=1:1:n dis2(i)=sqrt((S(i).xd-S2(n+1).xd)^(2)+(S(i).yd-S2(n+1).yd)^(2)); end end if(b>(3*c) && b<(4*c)) S3(n+1).xd=100; S3(n+1).yd=100; for i=1:1:n dis3(i)=sqrt((S(i).xd-S3(n+1).xd)^(2)+(S(i).yd-S3(n+1).yd)^(2)); end end if(b>(4*c) && b<(5*c)) S4(n+1).xd=100; S4(n+1).yd=0; for i=1:1:n dis4(i)=sqrt((S(i).xd-S4(n+1).xd)^(2)+(S(i).yd-S4(n+1).yd)^(2)); end end end % To estimate required transmitter power for given RSSI. for r=1:1:100 Pt(r,1)=(power(((Rssi(r,1)+40)/12),2.91))+((10*log((SNR)*(3.16)*(1.38*(10^(-23)))*(300)*(83500000)*(power(102.61*dis(r),2))))+5+30); end % To estimate required transmitter power for given RSSI. for r=1:1:100 Pt1(r,1)=(power(((Rssi(r,1)+40)/12),2.91))+((10*log((SNR)*(3.16)*(1.38*(10^(-23)))*(300)*(83500000)*(power(102.61*dis1(r),2))))+5+30); end % To estimate required transmitter power for given RSSI. for r=1:1:100 Pt2(r,1)=(power(((Rssi(r,1)+40)/12),2.91))+((10*log((SNR)*(3.16)*(1.38*(10^(-23)))*(300)*(83500000)*(power(102.61*dis2(r),2))))+5+30); end % To estimate required transmitter power for given RSSI. for r=1:1:100 Pt3(r,1)=(power(((Rssi(r,1)+40)/12),2.91))+((10*log((SNR)*(3.16)*(1.38*(10^(-23)))*(300)*(83500000)*(power(102.61*dis3(r),2))))+5+30); end % To estimate required transmitter power for given RSSI. for r=1:1:100 Pt4(r,1)=(power(((Rssi(r,1)+40)/12),2.91))+((10*log((SNR)*(3.16)*(1.38*(10^(-23)))*(300)*(83500000)*(power(102.61*dis4(i),2))))+5+30); end %Generating scenario for region A. countAT=0; countBT=0; Anumber=ceil(9.*rand(1,1)); NdesiredAmax=countH-5; NcurrentA(rmax,1)=((countH-Anumber)/countH)*100; for t=1:1:(countH-Anumber) if (RssiAAAmean<=RssiAAA(t)) if ((countH-Anumber)>=NdesiredAmax) RssinewA(t,1)=RssiAAAmean; end if ((countH-Anumber)<NdesiredAmax) RssinewA(t,1)=RssiAAA(t,1); end countAT=countAT+1; end if (RssiAAAmean>RssiAAA(t)) % if ((countH-Anumber)<NdesiredAmax) % RssinewA(t,1)=RssiAAAmean; % end % if ((countH-Anumber)>=NdesiredAmax) RssinewA(t,1)=RssiAAA(t,1); % end countBT=countBT+1; end end PlevelnewA=power(((RssinewA+40)/12),2.91); %Generating scenario for region B. countAT1=0; countBT1=0; Bnumber=ceil(9.*rand(1,1)); NdesiredBmax=countM-5; NcurrentB(rmax,1)=((countM-Bnumber)/countM)*100; for v=1:1:(countM-Bnumber) if (RssiBBBmean<=RssiBBB(v)) if ((countM-Bnumber)>=NdesiredBmax) RssinewB(v,1)=RssiBBBmean; end if ((countM-Bnumber)<NdesiredBmax) RssinewB(v,1)=RssiBBB(v,1); end countAT1=countAT1+1; end if (RssiBBBmean>RssiBBB(v)) % if ((countM-Bnumber)<NdesiredBmax) % RssinewB(v,1)=RssiBBBmean; % end % if ((countM-Bnumber)>=NdesiredBmax) RssinewB(v,1)=RssiBBB(v,1); % end countBT1=countBT1+1; end end PlevelnewB=power(((RssinewB+40)/12),2.91); %Generating scenario for region C. countAT2=0; countBT2=0; Cnumber=ceil(9.*rand(1,1)); NdesiredCmax=countL-5; NcurrentC(rmax,1)=((countL-Cnumber)/countL)*100; for x=1:1:(countL-Cnumber) if (RssiCCCmean<=RssiCCC(x)) if ((countL-Cnumber)>=NdesiredCmax) RssinewC(x,1)=RssiCCCmean; end if ((countL-Cnumber)<NdesiredCmax) RssinewC(x,1)=RssiCCC(x,1); end countAT2=countAT2+1; end if (RssiCCCmean>RssiCCC(x)) % if ((countL-Cnumber)<NdesiredCmax) % RssinewC(x,1)=RssiCCCmean; % end % if ((countL-Cnumber)>=NdesiredCmax) RssinewC(x,1)=RssiCCC(x,1); % end countBT2=countBT2+1; end end PlevelnewC=power((