توزیع نودها در فضاهای نامتقارن(Distribution of nodes in asymmetric spaces)

اغلب نواحی که برای شبیه سازی شبکه‌‌های بی‌سیم مورد استفاده قرار می‌گیرند، نواحی متقارن و عموما مربعی شکل می‌باشند ، حال آنکه به شکل طبیعی امکان توزیع متقارن نودها وجود نداره، در این پست سعی کردیم که کد توزیع نودها در فضاهای غیرمتقارن را منتشر کنیم، امیدواریم که براتون مفید واقع شده باشه.

در توزیع نودها از روش انتشار نقاط منحنی استفاده شده است.

کد متلب برای توضیحات بالا

%%%%%%%%%%%%%%%%%%%%%%%%%%-----------software details-----------%%%%%%%%%%%%%%%%%%%%%%
%devloped in Iran(ardebil)                                       
%Copyright (C) 2020(1398) by moussa hasanzadeh as the head devloper %%
%All rights reserved.                                            
%programed in mathlab 2016                                 
%gmail:mo30no@gmail.com
%phone 09147082079

%% clear commands
delete(allchild(0));% delete all figures and graphical objects
clear;% Remove items from workspace, freeing up system memory
clc;% clear the Console

%% path commands
file=matlab.desktop.editor.getActive;% get current script address
try % if dont occur any error this block will run
      fileDetail=dir(file.Filename);% get current script details
      folderName=fileDetail.folder;% get current script folder name
catch% else if an  error occurs then
      index=strfind( file.Filename,'\');% find index of back space(\) in path text
      folderName=file.Filename(1:index(end)-1);% select all path text from 1 to last Occurrence of \ as script folder name
end
paths=genpath(folderName);% make current path and all sub paths
addpath(paths);% add all paths in known paths of matlab
cd(folderName);% go to path that current script runned from that


%% net building  commands
net=[];% create an empty var as net object
net.para.width=1000;% width param of net
net.para.length=1000;% width param of net
net.para.height=5;% width param of net
net.para.side=6;% side of Polygon
net.para.nodNum=1000 ;% number of nodes of net work
net.para.simPer=1000;% number of simulation period
net.para.sInd=0;%period counter



 
ph=pi/net.para.side;% used for shifting of area to be Attractive!
t=linspace(0,1,net.para.side+1);% Generate linearly spaced vector( used for making polygon)
x1=sin(2*pi*t+ph); x1=(((x1-min(x1))/(max(x1)-min(x1)))-0.5)* net.para.length;% x coordinate of area(after normalising)
y1=-cos(2*pi*t+ph);y1=(((y1-min(y1))/(max(y1)-min(y1)))-0.5)*net.para.width;% y coordinate of area(after normalising)


  [x1,y1]=divider(x1,y1,10,1);% adding points between each pair of vertices

rs1=randsample(normalVect(rand(1,length(x1)),1,.5),length(x1));
figure;
x1=x1.*rs1;
y1=y1.*rs1;
[x1,y1]=divider(x1,y1,500);% adding points between each pair of vertices
x1(end+1)=x1(1);
y1(end+1)=y1(1);

plot(x1,y1,'r');hold on
rs=sqrt(randsample(normalVect(rand(1,length(x1)),.0,1),length(x1)));

plot(x1.*rs,y1.*rs,'.')






% give a vector or matrix for normalizing between tow number or 0 and 1
% vect=normalVect([1,-2,4,1,11,4,5])
% vect=normalVect([1,-2,4,1,11,4,5],-2,1)
% vect=normalVect([1,-2,4;11,4,5])

function vect=normalVect(varargin)
if nargin==0
    msgbox('give a 1 argoment atleast');
    vect=[];
elseif ~isnumeric(varargin{1})
    msgbox('abnormal argoment');
    vect=[];
else
   vect=varargin{1}; 
if nargin==1
    vect=(vect-min(vect(:)))./(max(vect(:))-min(vect(:)));
elseif nargin==3 ||  nargin>3
    if varargin{2}~=varargin{3}
        vect=(vect-min(vect(:)))./(max(vect(:))-min(vect(:)));
        vMin=min([varargin{2},varargin{3}]);
        vMax=max([varargin{2},varargin{3}]);
        vDif=vMax-vMin;
        vect=(vect*vDif)+vMin;
        if nargin>3
            msgbox('we used 3 argoment only');
        end
    else
        vect=(vect-min(vect(:)))./(max(vect(:))-min(vect(:)));
        msgbox('abnormal limits: normaled 0: 1');
    end
end

end
end





%version 1.1
%isFair=1 adding points is based on distance isFair=2 adding of points based on edges
function   [X,Y]=divider(X,Y,n,isFair)



try
    if isFair==1
    else
        isFair=0 ;
    end
catch
    
isFair=0;
end
switch isFair
    case 0
T1=0;
T2=0;

for i=1:length(X)-1
    temp1=0;
    temp1=linspace(X(i),X(i+1),n+2)  ;
    
    T1=[T1,temp1(1:end-1)];
end
T1=[T1,temp1(end)];
X=T1(2:end);


for i=1:length(Y)-1
    temp2=0;
    temp2=linspace(Y(i),Y(i+1),n+2)  ;
    
    T2=[T2,temp2(1:end-1)];
end
T2=[T2,temp2(end)];
Y=T2(2:end);
    case 1

dists= sqrt( (X(2:end)-X(1:end-1)).^2 + (Y(2:end)-Y(1:end-1)).^2) ;


n1=round((dists./min(dists))*n);
T1=0;
T2=0;

for i=1:length(X)-1
  
    temp1=0;
    temp1=linspace(X(i),X(i+1),n1(i)+2)  ;
    
    T1=[T1,temp1(1:end-1)];
end
T1=[T1,temp1(end)];
X=T1(2:end);


for i=1:length(Y)-1
    temp2=0;
    temp2=linspace(Y(i),Y(i+1),n1(i)+2)  ;
    
    T2=[T2,temp2(1:end-1)];
end
T2=[T2,temp2(end)];
Y=T2(2:end);
end


end

نمونه خروجی برای کد بالا(Output sample for the above code):