FW交通均衡求解算法

Q=1500;                                                  %OD需求
W=[30 33 25 30 30 35 20 23 20];                          %路段零流出行时间
C=[600 500 400 500 500 600 300 700 300];
Mxf=[0 1 0 0 0 0 1 1 0;0 1 1 0 0 1 0 0 0;0 1 1 0 0 0 0 1 1;1 0 0 0 1 0 0 1 1;1 0 0 0 1 1 0 0 0;0 0 0 1 0 0 0 1 1;0 0 0 1 0 1 0 0 0];
syms lambda real
for i=1:numx
    syms x（i） real;
end
cont=0;
e=inf;
x=x（1:7）;
X0=zeros（17）;
t=zeros（17）;
tt=t;
while e>1e-3                           %精度判断
    cont=cont+1;                       %迭代次数更新
    D（11）=W（11）*[1+0.15*（（x1（14）+x1（15））/C（11））^4]; 
D（12）=W（12）*[1+0.15*（（x1（11）+x1（12）+x1（13））/C（12））^4]; 
D（13）=W（13）*[1+0.15*（（x1（13）+x1（12））/C（13））^4]; 
D（14）=W（13）*[1+0.15*（（x1（16）+x1（17））/C（14））^4]; 
D（15）=W（13）*[1+0.15*（（x1（14）+x1（15））/C（15））^4]; 
D（16）=W（13）*[1+0.15*（（x1（12）+x1（15）+x1（17））/C（16））^4]; 
D（17）=W（13）*[1+0.15*（x1（11）/C（17））^4]; 
D（18）=W（13）*[1+0.15*（（x1（11）+x1（13）+x1（14）+x1（16））/C（18））^4]; 
D（19）=W（13）*[1+0.15*（（x1（13）+x1（14）+x1（16））/C（19））^4]; 
R（11）=D（12）+D（17）+D（18）;
R（12）=D（12）+D（13）+D（16）;
R（13）=D（12）+D（13）+D（19）+D（18）;
R（14）=D（11）+D（15）+D（19）+D（18）;
R（15）=D（11）+D（15）+D（16）;
R（16）=D（14）+D（19）+D（18）;
R（17）=D（14）+D（16）;
    Y1=Mxf（index:）.*Q;                %全有全无法求辅助流量
    S=Y1-X1;                           %搜索方向
    X2=X1+lambda*S;                    %先将X2用X1和lambda进行表达
    t=（W）.*（1+0.15*（X2./Cmax）.^4）;     %含lambd