Salama网络拓扑随机生成算法Matlab源码

Salama网络拓扑随机生成算法Matlab源码(2007-03-0610:19:00)
%Salama博士的网络拓扑随机生成算法
%Border_length----正方形区域的边长，单位：km
%Node_amount------网络节点的个数
%Alpha------------网络特征参数，Alpha越大，短边相对长边的比例越大
%Beta-------------网络特征参数，Beta越大，边的密度越大
%Sxy--------------用于存储节点的序号，横坐标，纵坐标的矩阵
%Cost-------------用于存储边的费用的邻接矩阵，费用在[2,10]之间随机选取，无边的取无穷大
%Delay------------用于存储边的时延的邻接矩阵，时延等于边的距离除以三分之二光速，无边的取无穷大
function[Sxy,Cost,Delay]=Net_Create(Border_length,Node_amount,Alpha,Beta)
%参数初始化
Sxy=zeros(3,Node_amount);
Cost=zeros(Node_amount,Node_amount);
Delay=Cost;
%在正方形区域内随机均匀选取Node_amount个节点
for i=1:Node_amount
Sxy(2,i)=Border_length*rand;
Sxy(3,i)=Border_length*rand;
end
%按横坐标由小到大的顺序重新为每一个节点编号
temp=Sxy;
Sxy2=Sxy(2,:);
Sxy2_sort=sort(Sxy2);
for i=1:Node_amount
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;
%在节点间随机产生边，并构造延时矩阵和费用矩阵
for i=1:(Node_amount-1)
for j=(i+1):Node_amount
Distance=((Sxy(2,i)-Sxy(2,j))^2+(Sxy(3,i)-Sxy(3,j))^2)^0.5;
P=Beta*exp(-Distance/(Alpha*Border_length));
if P>rand
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);
else
Delay(i,j)=inf;
Delay(j,i)=inf;
Cost(i,j)=inf;
Cost(j,i)=inf;
end
end
end
Net_plot(Border_length,Node_amount,Sxy,Cost,Delay)
%用于绘制网络拓扑的函数
function Net_plot(Border_length,Node_amount,Sxy,Cost,Delay)
%画节点
plot(Sxy(2,:),Sxy(3,:),'k.')
%设置图形显示范围
xlim([0,Border_length]);
ylim([0,Border_length]);
hold on
%为节点标序号
for i=1:Node_amount
Str=int2str(i);
text(Sxy(2,i),Sxy(3,i),Str);
hold on
end
%画边，并给边标注费用和延时
for i=1:(Node_amount-1)
for j=(i+1):Node_amount
if isinf(Cost(i,j))==0
plot([Sxy(2,i),Sxy(2,j)],[Sxy(3,i),Sxy(3,j)]);
xx=0.5*(Sxy(2,i)+Sxy(3,i));
yy=0.5*(Sxy(2,j)+Sxy(3,j));
Str1=num2str(Cost(i,j));
Str2=num2str(Delay(i,j));
Str1=Str1(1:3);
Str2=Str2(1:3);
%text(xx,yy,[Str1,',',Str2])
hold on
end
end
end。