DTN.rar_DTN matlab_dtn_matlab DTN_matlab mobile_user location

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% % DYNAMIC TRIANGULATION ALGORITHM % %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %program for generating received signal strength(rss) tic clear all clc r=input('enter the cell radius in meters='); N=input('enter the no of cells='); %number of sensor nods S=3; %finding the coordinates for S sensor nodes x(1)=1; y(1)=1; p(1,:)=[x(1),y(1)]; for i=2:N x(i)=(i-1)*1.5*r; h(i)=(i-1)*sqrt(3)*r; y(i)=h(i)^2-x(i)^2; p(i,:)=[x(i),y(i)]; end for i=1:S k=randint(1,1,[1 N]); P(i,:)=p(k,:); end disp(P) %finding the received signal strength captured by all sensor nodes theta=2*pi*rand; X=r*cos(theta); Y=r*sin(theta); M=[X,Y]; sum=0; for i=1:S for j=1:2 a=sum+(M(1,j)-P(i,j))^2; end d(i)=sqrt(a); end d0=10; n=4; MRSS=-10*n*log(d0); for i=1:S RSS(i)=(MRSS-10*n*log(d(i)/d0))/5; end %finding stongest rssi and generating mapping circle for i=1:S d(i)=d0*10^((MRSS-RSS(i))/(10*n)); end for i=1:S for j=1:2 t(i,j)=(d(i))^2+M(1,j); end end [v p]=min(d); xth =x(p)+d(p)*cos(theta); yth=y(p)+d(p)*sin(theta); %generate the cost function cost %finding minimum cost function and angle for k=1:359 if(p==1) cf(k)=sqrt(abs(e3(k)-e2(k))); elseif(p==2) cf(k)=sqrt(abs(e3(k)-e1(k))); else cf(k)=sqrt(abs(e2(k)-e1(k))); end subplot(3,1,3); plot(k,cf(k),'r.'); hold on grid on xlabel('angle in degrees') ylabel('cost function') title('cost function with respect to change in angle') end [CF ang]=min(cf); plot(ang,CF,'go'); legend('cost function','minimum cost function') theta=ang; D(1,:)=[x(p)+d(p)*cos(theta),y(p)+d(p)*sin(theta)]; figure, plot(t(1:3,1),t(1:3,2),'rs','MarkerEdgeColor','k',... 'MarkerFaceColor','g',... 'MarkerSize',4) grid on hold on xlabel('x co-ordinate of sensor nodes') ylabel('y co-ordinates of sensor nodes') legend plot(D(1,1),D(1,2),'bs','MarkerEdgeColor','k',... 'MarkerFaceColor','r',... 'MarkerSize',4) %results fprintf('RSSIs collected from 3 sensor nodes=') disp(RSS) fprintf('distance of mobile user from all 3 sensor nodes=') disp(d) fprintf('minimum cost function and angle=') disp([CF theta]) fprintf('co-ordinates of the new location of the mobile user are=') disp(D) toc