静态交通分配（Dail 算法）

最新推荐文章于 2024-03-18 09:26:47 发布

Yang@

最新推荐文章于 2024-03-18 09:26:47 发布

阅读量3.3k

点赞数 5

分类专栏：交通分配文章标签：交通分配 matlab Dial 算法

本文链接：https://blog.csdn.net/qq_39987070/article/details/89440359

版权

本文介绍了Dial算法在静态交通分配中的使用，包括主函数和其他辅助函数的实现，如路段权重和似然估计函数。通过这些函数，实现了路段流量的计算。文章简述了理论基础，并预告将深入探讨算法原理。

摘要由CSDN通过智能技术生成

理论最后介绍

先附上代码

主函数

%%输入初始信息
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 
%定义收点
for i=1:9
   m=find(D==i);
   Di{i}=[O(m)];
end
   Di{1}=1;

%计算likelihood
theta=1;
for i=1:9
    for j=1:9
        if r(i)<r(j) && s(i)>s(j)
            L(i,j)=exp(theta*(r(j)-r(i)-linkweight(i,j)));
        end 
    end 
end 
L(9,:)=zeros(1,9);%补全矩阵

%setp 2 先前计算权重,将L变为稀疏矩

for i=1:leng