【原创】改进的Salama网络拓扑随机生成算法通用Matlab程序
现有的几个网络拓扑随机发生器，其实很难生成理想的网络拓扑结构，其主要原因在于很难控制节点的疏密和间距。我们提出来的这个改进算法，在随机抛撒节点的时候使用了K均值聚类，由本算法作为网络拓扑发生器，网络节点分布均匀且疏密得当，边的分布也比较均衡，GreenSim团队免费推荐您使用。function =NetCreate(BorderLength,NodeAmount,Alpha,Beta,PlotIf,FlagIf)
%% 改进的Salama博士的网络拓扑随机生成算法
%% GreenSim团队原创作品，转载请注明
%% 欲与原作者进行技术交流，请发邮件至aihuacheng@gmail.com
%% 输入参数列表
% BorderLength-----正方形区域的边长，单位：km
% NodeAmount-------网络节点的个数
% Alpha------------网络特征参数，Alpha越大，短边相对长边的比例越大
% Beta-------------网络特征参数，Beta越大，边的密度越大
% PlotIf-----------是否画网络拓扑图，如果为1则画图，否则不画
% FlagIf-----------是否标注参数，如果为1则将标注边的参数，否则不标注
%% 输出参数列表
% Sxy--------------用于存储节点的序号，横坐标，纵坐标的矩阵
% AM---------------01存储的邻接矩阵，AM(i,j)=1表示存在由i到j的有向边
% Cost-------------用于存储边的费用的邻接矩阵，费用在之间随机选取，无边的取无穷大
% Delay------------用于存储边的时延的邻接矩阵，时延等于边的距离除以三分之二光速，无边的取无穷大
% DelayJitter------用于存储边的延时抖动的邻接矩阵，在1～3微秒之间随机选取，无边的取无穷大
% PacketLoss-------用于存储边的丢包率，在0～0.01之间随机选取，无边的取无穷大
%% 参考参数设置
% =NetCreate(10,25,10,20,1,1)
%%
%参数初始化
NN=10*NodeAmount;
SSxy=zeros(NN,2);
%在正方形区域内随机均匀选取NN个节点
for i=1:NN
SSxy(i,1)=BorderLength*rand;
SSxy(i,2)=BorderLength*rand;
end
=kmeans(SSxy,NodeAmount);
Sxy=[',C]';
%按横坐标由小到大的顺序重新为每一个节点编号
temp=Sxy;
Sxy2=Sxy(2,:);
Sxy2_sort=sort(Sxy2);
for i=1:NodeAmount
pos=find(Sxy2==Sxy2_sort(i));
if length(pos)>1
error('仿真故障，请重试！');
end
temp(1,i)=i;
temp(2,i)=Sxy(2,pos);
temp(3,i)=Sxy(3,pos);
end
Sxy=temp;
%在节点间随机产生边，并构造延时矩阵和费用矩阵
AM=zeros(NodeAmount,NodeAmount);
Cost=zeros(NodeAmount,NodeAmount);
Delay=zeros(NodeAmount,NodeAmount);
DelayJitter=zeros(NodeAmount,NodeAmount);
PacketLoss=zeros(NodeAmount,NodeAmount);
for i=1:(NodeAmount-1)
for j=(i+1):NodeAmount
Distance=((Sxy(2,i)-Sxy(2,j))^2+(Sxy(3,i)-Sxy(3,j))^2)^0.5;
P=Beta*exp(-Distance^5/(Alpha*BorderLength));
if P>rand
AM(i,j)=1;
AM(j,i)=1;
Delay(i,j)=0.5*Distance/100000;
Delay(j,i)=Delay(i,j);
Cost(i,j)=2+8*rand;
Cost(j,i)=Cost(i,j);
DelayJitter(i,j)=0.000001*(1+2*rand);
DelayJitter(j,i)=DelayJitter(i,j);
PacketLoss(i,j)=0.01*rand;
PacketLoss(j,i)=PacketLoss(i,j);
else
Delay(i,j)=inf;
Delay(j,i)=inf;
Cost(i,j)=inf;
Cost(j,i)=inf;
DelayJitter(i,j)=inf;
DelayJitter(j,i)=inf;
PacketLoss(i,j)=inf;
PacketLoss(j,i)=inf;
end
end
end
Net_plot(BorderLength,NodeAmount,Sxy,Cost,Delay,DelayJitter,PacketLoss,PlotIf,FlagIf)
%用于绘制网络拓扑的函数
function Net_plot(BorderLength,NodeAmount,Sxy,Cost,Delay,DelayJitter,PacketLoss,PlotIf,FlagIf)
%画节点
if PlotIf==1
plot(Sxy(2,:),Sxy(3,:),'k.')
%设置图形显示范围
xlim();
ylim();
hold on
%为节点标序号
for i=1:NodeAmount
Str=int2str(i);
text(Sxy(2,i)+BorderLength/100,Sxy(3,i)+BorderLength/100,Str,'FontSize',12);
hold on
end
end
%画边，并给边标注费用和延时
if PlotIf==1
for i=1:(NodeAmount-1)
for j=(i+1):NodeAmount
if isinf(Cost(i,j))==0
plot(,);
if FlagIf==1
xx=0.5*(Sxy(2,i)+Sxy(2,j));
yy=0.5*(Sxy(3,i)+Sxy(3,j));
Str1=num2str(Cost(i,j));
Str2=num2str(1000000*Delay(i,j));
Str3=num2str(1000000*DelayJitter(i,j));
Str4=num2str(100*PacketLoss(i,j));
Str1=Str1(1:3);
Str2=Str2(1:3);
Str3=Str3(1:3);
Str4=Str4(1:3);
text(xx,yy,['(',Str1,',',Str2,',',Str3,',',Str4,')'],'FontSize',6);
end
hold on
end
end
end
end 先收藏了。
需要用到了还请楼主多指教。
页:
[1]