You can create your script for the sensor. Only adding noise (not necessarily gaussian) to each sensor in the network for the dynamic model.

You can define random r and theta, for example N users in a Circular area with radius R and a uniform distribution for users:

r=R*rand(1,N);

theta=2*pi*rand(1,N);

i=1:N %user number

%location of each user:

x(1,i)=r(1,i)*cos(theta(1,i));

y(1,i)=r(1,i)*sin(theta(1,i));

%D, coverage range

If |r(1,i)-r(1,j)|