prim.m

global line_info dot_array;
setglobal2();
arr = zeros(12,12);
for i=1:12
    for j=i+1:12
        arr(i,j)=-line_info(dot_array(i,j),3);
    end
end
A=arr
A = A+A';
result=prim1(A);
arr1 = zeros(1,11);
for i =1:11
    arr1(1,i)=dot_array(result(1,i),result(2,i));
end
value = get_value(arr1);
show_graph(arr1,value);






function result=prim1(a)
% 输入：a—邻接矩阵，a（i，j)是指i到j之间的权值
% 输出：result—第一、二、三行分别代表最小生成树边的起点、终点、权集合
    a(a==0)=inf;
    result=[];
    p=1;
    tb=2:length(a);
    while size(result,2)~=length(a)-1
       temp=a(p,tb);
       temp=temp(:);
       d=min(temp);
       [jb,kb]=find(a(p,tb)==d);
       j=p(jb(1));k=tb(kb(1));
       result=[result,[j;k;d]];
       p=[p,k];
       tb(find(tb==k))=[];
    end
end
function setglobal2()
    global line_info;
    global CityLoc CityPop Clist dot_array CityPop2;
    CityLoc=[39.91667 116.41667;%北京
        45.75000 126.63333;%哈尔滨
        43.45 87.36;%乌鲁木齐
        34.26667,108.95000;%西安
        34.76667,113.65000;%郑州
        31.14 121.29;%上海
        30.35 114.17;%武汉
        29.35 106.33;%重庆
        30.40 104.04;%成都
        29.39 91.08;%拉萨
        25.04 102.42;%昆明
        23.16667,113.23333];%广东
    CityPop = [1961.24 1063.60 311.03 846.78 862.65 2301.91 978.54 2884.62 1404.76...
        55.94 643.20 1035.79]/100;
    CityPop2 = [12397 10910 2398 7691 23161 21235 17299 3392 8262 331 13164 15790]/100;
    axesm utm   %设置投影方式，这是MATLAT自带的Universal Transverse Mercator （UTM）方式
    Z=utmzone(CityLoc);%utmzone根据latlon20里面的数据选择他认为合适的投影区域，可以是一个台站的经纬度，也可以是所有台站的经纬度（此时是平均）
    setm(gca,'zone',Z)
    h = getm(gca);
    R=zeros(size(CityLoc));
    for i=1:length(CityLoc)
        [x,y]= mfwdtran(h,CityLoc(i,1),CityLoc(i,2));
        Clist(i,:)=[x;y]/1e3;
    end
    k=1;
    for i =1:12
        for j = i+1:12
            dot_array(i,j)=k;
            dot_array(j,i)=k;
            k=k+1;
        end
    end
    CiteNum = 12;
    BitNum=(CiteNum-1)*CiteNum/2; %you chan choose 10, 30, 50, 75
    line_value = zeros(12,12);
    SnrReq = [Inf 3000 1200 600];
    RateTable = [0 8 16 32];
    city_num = 12;
    pop = zeros(12,12);
    flow = zeros(12,12);
    for ii =1:city_num
        for jj=ii+1:city_num
               Dist = distance(CityLoc(ii,1),CityLoc(ii,2),CityLoc(jj,1),CityLoc(jj,2))*pi/180*6371;
               %Dist= norm(CityLoc(ii,:)-CityLoc(jj,:));
               Ind = find(SnrReq>Dist);
               ModRate = RateTable(Ind(end));
               line_value(ii,jj)=sqrt(CityPop2(ii)*CityPop2(jj))*ModRate;
               pop(ii,jj) = sqrt(CityPop2(ii)*CityPop2(jj));
               flow(ii,jj) = ModRate;
        end
    end
    line_info = zeros(66,11);%两点序号 最大价值 平均人口 最大流量 实际流量 是否实际联通 是否虚拟联通 调整次数
    k=1;
    for i=1:city_num
        for j=i+1:city_num
            line_info(k,1:2)=[i j];
            line_info(k,3)=line_value(i,j);
            line_info(k,4)=pop(i,j);
            line_info(k,5)=flow(i,j);
            k=k+1;
        end
    end

    
end