静态交通分配采用Dial算法

    %%输入初始信息
    O=[1,1,2,2,3,4,4,5,5,6,7,8,9];%网络节点
    D=[2,4,3,5,6,5,7,6,8,9,8,9,9];%网络节点 
    W=[2,2,2,2,2,1,2,1,2,2,2,2,0];%权重
    linkweight=sparse(O,D,W);%构建稀疏矩阵
    
    %起点
    origin=1;
    destination=9;
    
    
    num_nodes=length(O);
    %求起点到终点最短路
    %dis=zeros(9,9);
    for i=1:9  %可改
    for j=1:9
    [dist,path]=graphshortestpath(linkweight,i,j);%最短路工具箱
    dis(i,j)=dist;
    end 
    end 
    r=dis(1,:);%起点到任何一节点的最短路
    s=dis(:,9);%任何一点到终点的最短路
    
    %定义发点
    for i=1:9
       m=find(O==i);
       Oi{i}=[D(m)];
    end 
    %定