博客
关于我
CDS最小支配集的MATLAB仿真
阅读量：247 次
发布时间：2019-03-01

本文共 5034 字，大约阅读时间需要 16 分钟。

clc;

clear;
close all;

Node_Num=200;%size of nodes: 100

Posxy=zeros(Node_Num,2);%node_coordinate.the(:,1)for the x, and (:,2) for y
Radius=10;
Scale=100;
Area_Size_W=100;%area: square 300*300

d1=zeros(Node_Num,Node_Num+1);%store one_hop neighbors information

%%

%*************************Intialize the network***************************

for i=1:Node_Num

    Posxy(i,:)=Scale*rand(1,2);
end

%%

%***************node get the one_hop neighbors information****************

for i=1:Node_Num

    num1=0;
    for j=1:Node_Num
        if i~=j
            dist=(Posxy(i,1)-Posxy(j,1))^2+...
                (Posxy(i,2)-Posxy(j,2))^2;
            dist=sqrt(dist);
            if dist<Radius&&dist>0
                num1=num1+1;
                d1(i,num1+1)=j;
            end
        end
    end
    d1(i,1)=num1;
end

%%

%*********************show the connected picture**************************

hold on

for i=1:Node_Num%draw the line between two nodes
    for j=1:d1(i,1)
        plot([Posxy(i,1),Posxy(d1(i,j+1),1)],...
            [Posxy(i,2),Posxy(d1(i,j+1),2)],'-ok');
        hold on;
    end
end

%%

%************node get the two_hop neighbors information*******************

Max_Degree=max(d1(:,1));

Neighbor_hop2=zeros(Max_Degree,Max_Degree+1,Node_Num);%store two_hop neighbors information
for i=1:Node_Num
   for j=1:d1(i,1)
       Neighbor_hop2(j,1,i)=d1(i,j+1);
       for m=1:d1(d1(i,j+1),1)
           Neighbor_hop2(j,m+1,i)=d1(d1(i,j+1),m+1);
       end
   end
end

%

%**********color the node with gray and white based on the rule**********
%0-white;1-gray;

Color_N=zeros(Node_Num,1);%intial the color of nodes

for i=1:Node_Num

    if d1(i,1)<2
        Color_N(i,1)=0;
    else
        for j=1:d1(i,1)
            color=0;
            for n=j+1:d1(i,1)
                state=0;
                for m=2:d1(Neighbor_hop2(j,1,i),1)+1
                    if Neighbor_hop2(n,1,i)==Neighbor_hop2(j,m,i)
                        state=1;
                        break;
                    end
                end
                if state==0
                    color=1;
                    break;
                end
            end
            if color==1
                break;
            end
        end
        Color_N(i,1)=color;
    end
end

%%

%**********color Black based on the rule**********
%0-white;1-gray;2-balck;3-gray';

%Balck

for i=1:Node_Num
    if Color_N(i,1)==0&&d1(i,1)>0
        temp_degree=0;
        temp_id=0;
        for j=1:d1(i,1)
            if Color_N(d1(i,j+1),1)~=0
                if Color_N(d1(i,j+1),1)~=2
                    Color_N(d1(i,j+1),1)=3;
                end
                if d1(d1(i,j+1),1)>temp_degree
                    temp_degree=d1(d1(i,j+1),1);
                    temp_id=d1(i,j+1);
                elseif d1(d1(i,j+1),1)==temp_degree
                    if temp_id>d1(i,j+1)
                        temp_degree=d1(d1(i,j+1),1);
                        temp_id=d1(i,j+1);
                    end
                end
            end
        end
        if temp_id~=0
            Color_N(temp_id,1)=2;
        end
    end
end
%end Balck

%%

%%**********color Real_Violet based on the rule**********
%0-white;1-gray;2-balck;3-gray';
%remove the Violet node form to get a dominating set
%1 and 3 组成DS

Count2=0;

for m=1:Node_Num
    if Color_N(m,1)==1
        Count2=Count2+1;
    end
end
while(Count2)%Simulation of distributed computing
    for i=1:Node_Num
        if Color_N(i,1)==1
            temp_degree=d1(i,1);
            state=1;
            for j=1:d1(i,1)
                if Color_N(d1(i,j+1))==1
                    if temp_degree>d1(d1(i,j+1),1)
                        state=1;
                    elseif temp_degree==d1(d1(i,j+1),1)
                        if i<d1(i,j+1)
                            state=1;
                        else
                            state=0;
                            break;
                        end
                    else
                        state=0;
                        break;
                    end
                end
            end
            if state==1
                Color_N(i,1)=4;
                Count2=Count2-1;
                for j=1:d1(i,1)
                    if Color_N(d1(i,j+1))==1
                        Color_N(d1(i,j+1),1)=5;
                        Count2=Count2-1;
                    end
                end
            end
        end
    end
end

%%

%choes some node to form the connect nodes who's color number eqult to 2 or
%4
tempall=1;
temp_self=Neighbor_hop2(:,:,1)*0;%2 hop's information

for i=1:Node_Num

    if Color_N(i,1)==4||Color_N(i,1)==2

        temp_self=temp_self*0;%2 hop's information
        
        %node i formulation the 2_hop connected graph
        for n=1:d1(i,1)
            temp2=d1(d1(i,n+1),1);
            temp_count=1;
            state=1;
            for m=2:temp2+1
                temp_self(n,1)=d1(i,n+1);
                temp=Neighbor_hop2(n,m,i);
                state=1;
                for p=1:d1(i,1)
                    if temp==i
                        state=0;
                        break;
                    elseif temp==d1(i,p+1)
                        state=0;
                        break;
                    end
                end
                if state==1

                    temp_count=temp_count+1;

                    temp_self(n,temp_count)=temp;

                end

            end
        end
        
        
        
        
%%        
        %chose some node in two_hops of node i to connect the2_hop's dominating nodes
        Already_handle=zeros(Max_Degree*Max_Degree,1);
        Already_handle_result=Already_handle;
        handle_count=0;
        for n=1:d1(i,1)
            if temp_self(n,1)==0;
                break;                
            end
            for m=2:Max_Degree+1
                if temp_self(n,m)==0
                    break;
                end
                if Color_N(temp_self(n,m),1)==4||Color_N(temp_self(n,m),1)==2
                    state=0;
                    for p=1:handle_count
                        if Already_handle(p,1)==temp_self(n,m)
                            state=1;
                            if Already_handle_result(p,1)==0
                                break;
                            else
                                if Color_N(temp_self(n,1))==2||Color_N(temp_self(n,1))==4
                                    Already_handle_result(p,1)=0;
                                    break;
                                elseif  d1(temp_self(n,1),1)>d1(Already_handle_result(p,1),1)
                                    Already_handle_result(p,1)=temp_self(n,1);
                                elseif d1(temp_self(n,1),1)==d1(Already_handle_result(p,1),1)
                                    if temp_self(n:1)<Already_handle_result(p,1)
                                        Already_handle_result(p,1)=temp_self(n,1);
                                    end
                                end
                            end
                        end
                    end
                    if state==0;
                        if Color_N(temp_self(n,1))==2||Color_N(temp_self(n,1))==4
                            Already_handle(handle_count+1,1)=temp_self(n,m);
                            Already_handle_result(handle_count+1,1)=0;
                            handle_count=handle_count+1;
                        else
                            Already_handle(handle_count+1,1)=temp_self(n,m);
                            Already_handle_result(handle_count+1,1)=temp_self(n,1);
                            handle_count=handle_count+1;
                        end
                    end
                end
            end
        end
        
        for n=1:handle_count
           if  Already_handle_result(n,1)==0
              
           else
               Color_N(Already_handle_result(n,1))=6;
               
               temp1=Already_handle_result(n,1);
               temp2=Already_handle(n,1);
               plot([Posxy(i,1),Posxy(temp1,1)],[Posxy(i,2),Posxy(temp1,2)],'o-.c')
               plot([Posxy(temp2,1),Posxy(temp1,1)],[Posxy(temp2,2),Posxy(temp1,2)],'o-.c')
               
           end
        end
 

D-38

转载地址：http://ptwx.baihongyu.com/

你可能感兴趣的文章
喝酒易醉，品茶养心，人生如梦，品茶悟道，何以解忧？唯有杜康！-- 愿君每日到此一游！
当前时间: 2025-07-17 14:15:13   当前IP: 10.2.23.86   联系邮箱:javaeecc@qq.com     Copyright © 2020 - 2025 baihongyu.com 京ICP备2021015314号-2