無標度網絡matlab建模-資料下載頁

　　

【正文】 mp。(r2cum_degree(i+1)) choose(2) = i。 break end end while choose(2) == choose(1) %第一個點和第二個點相同的話，重新擇優 r2= rand(1)*total_degree。 for i= 1:iter1 if (r2=cum_degree(i))amp。amp。(r2cum_degree(i+1)) choose(2) = i。 break end end end % 選出第三個和新點相連接的頂點 r3= rand(1)*total_degree。 for i= 1:iter1 if (r3=cum_degree(i))amp。amp。(r3cum_degree(i+1)) choose(3) = i。 break end end while (choose(3)==choose(1))||(choose(3)==choose(2)) r3= rand(1)*total_degree。 for i=1:iter1 if (r3=cum_degree(i))amp。amp。(r3cum_degree(i+1)) choose(3) = i。 break end end end %新點加入網絡后, 對鄰接矩陣進行更新 for k = 1:m adjacent_matrix(iter,choose(k)) = 1。 adjacent_matrix(choose(k),iter) = 1。 end node_degree=zeros(1,iter+1)。 node_degree(2:iter+1) = sum(adjacent_matrix)。endmatrix = adjacent_matrix。function tu_plot(rel,control)%由鄰接矩陣畫連接圖，輸入為鄰接矩陣rel，必須為方陣； %control為控制量，0表示畫出的圖為無向圖，1表示有向圖。默認值為0r_size=size(rel)。 %a=size(x)返回的是一個行向量，該行向量第一個元素是 %x的行數，第2個元素是x的列數if nargin2 %nargin是用來判斷輸入變量個數的函數 control=0。 %輸入變量小于2，即只有一個，就默認control為0endif r_size(1)~=r_size(2) %行數和列數不相等，不是方陣，不予處理 disp(39。Wrong Input! The input must be a square matrix!39。)。 return。endlen=r_size(1)。rho=10。 %限制圖尺寸的大小r=5/^len。 %點的半徑theta=0:(2*pi/len):2*pi*(11/len)。[pointx,pointy]=pol2cart(theta39。,rho)。theta=0:pi/36:2*pi。[tempx,tempy]=pol2cart(theta39。,r)。point=[pointx,pointy]。hold onfor i=1:len temp=[tempx,tempy]+[point(i,1)*ones(length(tempx),1),point(i,2)*ones(length(tempx),1)]。 plot(temp(:,1),temp(:,2),39。b39。)。 text(point(i,1),point(i,2),num2str(i))。 %畫點endfor i=1:len for j=1:len if rel(i,j) link_plot(point(i,:),point(j,:),r,control)。 %連接有關系的點 end endendset(gca,39。XLim39。,[rhor,rho+r],39。YLim39。,[rhor,rho+r])。axis offfunction link_plot(point1,point2,r,control) %連接兩點temp=point2point1。if (~temp(1))amp。amp。(~temp(2)) return。 %不畫子回路endtheta=cart2pol(temp(1),temp(2))。[point1_x,point1_y]=pol2cart(theta,r)。point_1=[point1_x,point1_y]+point1。[point2_x,point2_y]=pol2cart(theta+(2*(thetapi)1)*pi,r)。point_2=[point2_x,point2_y]+point2。if control arrow(point_1,point_2)。else plot([point_1(1),point_2(1)],[point_1(2),point_2(2)],39。m39。)。end