用户均衡流量分配模型UE的matlab求解程序

clear;
clc;
syms lambda;
numf=4;%路径总数
numx=7;%路段总数
Q1=3000;%OD对1总流量
Q2=3000;%OD对2总流量

t0=[8 9 2 6 3 3 4];%路段走行时间函数参数
Sa=[2000 2000 2000 4000 2000 2500 2500];%路段走行时间函数参数
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];%路段转路径矩阵
%Mfx=Mxf;%路径转路段矩阵
%%以上为程序需要输入的变量
xx=zeros（1numx）;
t=zeros（1numx）;
t=t0+t0.*0.15.*（xx./Sa）.^4;%路段走行时间函数
ft=（Mxf*t‘）‘;%几条路径的走行时间初值，前两个为OD对1的，后两个为OD对2的路径
N=15;%最大迭代次数
[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