交通分配之用户均衡分配模型之三(matlab源码)

例
总流量为100,走行函数为：
⎪⎭⎫ ⎝⎛+=40)(6.04)(111t x x c ⎪⎭
⎫ ⎝⎛+=40)(9.06)(222t x x c ⎪⎭
⎫ ⎝⎛+=60)(3.02)(333t x x c ⎪⎭⎫ ⎝⎛+=40)(75.05)(444t x x c ⎪⎭
⎫ ⎝⎛+=40)(45.03)(555t x x c 模型求解的Matlab 源码：
syms lambda ;
numf = 3; %路径总数
numx = 5;%路段总数
Q=100;%总流量
fid=fopen('D:\Program Files\MATLAB\R2011b\bin\我的matlab\traffic\UECOM.txt','w'); %设置运行结果输出文件
T = [4 6 2 5 3 ]; %路段走行时间函数参数
cap = [(0.6/40) (0.9/40) (0.3/60) (0.75/40) (0.45/40) ]; %路段走行时间函数参数 Mxf = [1 0 0 1 0;
0 1 0 0 1 ;
1 0 1 0 1]; % 路段转路径矩阵
% Mfx = Mxf'; % 路径转路段矩阵
%========================================================== %以上为程序需要输入的变量
xx= zeros(1,numx);
t = zeros(1,numx);
t = T + cap .* xx ;%路段走行时间函数
ft = (Mxf * t')'; %三条路径的走行时间初值。

路径1为路段1,4 ，路径2为路段2,5 ，路径3为路段1,3,5
N= 15; %最大迭代次数，也可使用其他收敛条件
[Min,index] = min(ft) ;
xx = Mxf(index,:).*Q ; % 全有全无法为最短路径上的路段分配流量
for i =1 :N
y = zeros(numx); %辅助路段流
t = T + cap .* xx ;%路段时间
ft = (Mxf * t')' ; %路径时间
fprintf(fid,'\n========================================================== =================\n' );
fprintf(fid,'\n\n第%2d 次迭代的路径时间值：\n' , i );
for (ii = 1: numf)
fprintf(fid,'%11.4f ：' ,ft(1,ii) );
end
[Min,index] = min(ft) ;
y = Mxf(index,:).*Q ; % 全有全无法为最短路径上的路段分配流量
fprintf(fid,'\n第%d 次迭代的辅助流量值是：\n' , i);
for (ii = 1: numx)
fprintf(fid,'%11.4f: ' ,y(1,ii) );
end
zz = xx + lambda * (y-xx); % 按方向(y-xx)进行一维搜索,步长为lamda
t = T + cap .* zz ;
f = sum( (y -xx).* t ,2);
lambda1 =double( solve(f)) ; %求解方程，确定步长。

k = length(lambda1); % 如步长lambda1的解不唯一，取实数，且大于0 小于1；
for m=1: k
if lambda1(m,1) > 0 && lambda1(m,1) < 1 && isreal(lambda1(m,1))
lambda2 =lambda1(m,1);
end
end
xx = xx + lambda2*(y - xx ); % 得到下一步的流量值，且进行下一次迭代
fprintf(fid,'\n第%2d 次迭代后的路段流量值是：\n' , i);
for (ii = 1: numx)
fprintf(fid,'%11.4f: ' ,xx(1,ii) );
end
end
status = fclose('all')
ft。