clear;
clc;
syms?lambda;
numf=4;%路徑總數(shù)
numx=7;%路段總數(shù)
Q1=3000;%OD對1總流量
Q2=3000;%OD對2總流量

t0=[8?9?2?6?3?3?4];%路段走行時間函數(shù)參數(shù)
Sa=[2000?2000?2000?4000?2000?2500?2500];%路段走行時間函數(shù)參數(shù)
Mxf=[1?0?0?0?0?0?0;0?0?1?1?0?1?0;0?0?0?1?1?0?1;0?1?0?0?0?0?0];%路段轉(zhuǎn)路徑矩陣
%Mfx=Mxf;%路徑轉(zhuǎn)路段矩陣
%%以上為程序需要輸入的變量
xx=zeros（1numx）;
t=zeros（1numx）;
t=t0+t0.*0.15.*（xx./Sa）.^4;%路段走行時間函數(shù)
ft=（Mxf*t‘）‘;%幾條路徑的走行時間初值，前兩個為OD對1的，后兩個為OD對2的路徑
N=15;%最大迭代次數(shù)
[Min1index1]=min（ft（:1:2））;
[Min2index2]=min（ft（:3:4））;index2=index2+2;%Min為最短路徑時間，index為最短路徑的位置
xx1=Mxf（index1:）.*Q1;
xx2=Mxf（index2:）.*Q2;%全有全無法為最短路徑上的路段分配流量
xx=xx1+xx2;
for?i=1:N
????y=zeros（numx）;
????y1=zeros（numx）;%輔助路段流
????y2=zero