北航最优化方法大作业参考教学幻灯片.docxVIP

1 流量工程问题 问题重述 定义一个有向网络G=(N,E)，其中N是节点集，E是弧集。令A是网络G的点弧关联矩阵，即N×E阶矩阵，且第l列与弧里(I,j)对应，仅第i行元素为1，第j行元素为-1，其余元素为0。再令bm=(bm1,…,bmN)T，fm=(fm1,…,fmE)T，则可将等式约束表示成： Afm=bm 本算例为一经典TE算例。算例网络有7个节点和13条弧，每条弧的容量是5个单位。此外有四个需求量均为4个单位的源一目的对，具体的源节点、目的节点信息如图所示。这里为了简单，省区了未用到的弧。此外，弧上的数字表示弧的编号。此时，c=((5,5…,5)1×13)T， 根据上述四个约束条件，分别求得四个情况下的最优决策变量x=((x12,x13,…,x75)1×13)。 图 SEQ 图 \* ARABIC 1 网络拓扑和流量需求 7节点算例求解 算例1（b1=[4;-4;0;0;0;0;0]T） 转化为线性规划问题： Minimize cTx1 Subject to Ax1=b1 x1=0 利用Matlab编写对偶单纯形法程序，可求得: 最优解为x1*=[4 0 0 0 0 0 0 0 0 0 0 0 0]T 对应的最优值cTx1=20 算例2（b2=[4;0;-4;0;0;0;0]T） Minimize cTx2 Subject to A