现有的几个网络拓扑随机发生器，其实很难生成理想的网络拓扑结构，其主要原因在于很难控制节点的疏密和间距

function [Sxy,AM,Cost,Delay,DelayJitter,PacketLoss]=NetCreate(BorderLength,NodeAmount,Alpha,Beta,PlotIf,FlagIf) %% 改进的Salama博士的网络拓扑随机生成算法 % BorderLength-----正方形区域的边长，单位：km % NodeAmount-------网络节点的个数 % Alpha------------网络特征参数，Alpha越大，短边相对长边的比例越大 % Beta-------------网络特征参数，Beta越大，边的密度越大 % PlotIf-----------是否画网络拓扑图，如果为1则画图，否则不画 % FlagIf-----------是否标注参数，如果为1则将标注边的参数，否则不标注 %% 输出参数列表 % Sxy--------------用于存储节点的序号，横坐标，纵坐标的矩阵 % AM---------------01存储的邻接矩阵，AM(i,j)=1表示存在由i到j的有向边 % Cost-------------用于存储边的费用的邻接矩阵，费用在[2,10]之间随机选取，无边的取无穷大 % Delay------------用于存储边的时延的邻接矩阵，时延等于边的距离除以三分之二光速，无边的取无穷大 % DelayJitter------用于存储边的延时抖动的邻接矩阵，在1～3微秒之间随机选取，无边的取无穷大 % PacketLoss-------用于存储边的丢包率，在0～0.01之间随机选取，无边的取无穷大 %% 参考参数设置 % [Sxy,AM,Cost,Delay,DelayJitter,PacketLoss]=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 [IDX,C]=kmeans(SSxy,NodeAmount); Sxy=[[1:NodeAmount]',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([0,BorderLength]); ylim([0,BorderLength]); 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([Sxy(2,i),Sxy(2,j)],[Sxy(3,i),Sxy(3,j)]); 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