最新北航最优化方法大作业参考.docxVIP

PAGE PAGE # X=((X X=((X12,X13,…,X75)1 X13)。 in =■ 1:［亠 2. d\sA 瞒 m 3: 3- 丄山 1流量工程问题 问题重述 定义一个有向网络G=(N,E)，其中N是节点集，E是弧集。令A是网络G的点弧关 联矩阵，即N X E阶矩阵，且第I列与弧里(l,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 X13)T， 4 -4 0 0 ,如一 4 ■ 0 —4 0 ,ba = 0 —电 4 D \切一 4 ■ -0 0 n 0 0 0 o 0 0 0 0 0 0 0 —4 根据上述四个约束条件，分别求得四个情况下的最优决策变量 图1网络拓扑和流量需求 7节点算例求解 11-1 1 1 -1 -1 1 1 -1 -1 -1 1 1 -1 -1 1 -1 算例 1 (b 算例 1 (b仁[4;-4;0;0;0;0;0] T) 转化为线性规划问题： Mini mize cTx1 Subject to Ax 仁b1 x1=0 利用Matlab编写对偶单纯形法程序，可求得： 最优解为 x1*=[4 0 0 0 0 0 0 0 0 0 0 0 0] T 对应的最优值cTx仁20 算例 2( b2=[4;0;-4;0;0;0;0] T) Mini mize cTx2 Subject toAx2=b2 Subject to Ax2=b2 X2=0 利用Matlab编写对偶单纯形法程序，可求得 最优解为 x2*=[0 4 0 0 0 0 0 0 0 0 0 0 0] T 对应的最优值cTx2=20 算例 3( b3=[0;-4;4;0;0;0;0] T) Mi ni mize Jx3 Subject to Ax3=b3 X3=0 利用Matlab编写对偶单纯形法程序，可求得： 最优解为 x3*=[4 0 0 0 4 0 0 0 0 0 0 0 0] T 对应的最优值cTx3=40 算例 4 (b4=[4;0;0;0;0;0;-4] T) Mi ni mize Jx4 Subject to Ax4=b4 X4=0 利用Matlab编写对偶单纯形法程序，可求得： 最优解为 x4*=[4 0 0 4 0 0 0 0 0 4 0 0 0] T 对应的最优值Jx4=60 1.3计算结果及结果说明 算例 1( b1=[4;-4;0;0;0;0;0] T) 算例1中，由b1可知，节点2为需求节点，节点1为供给节点，由节点1将信息传 输至节点2的最短路径为弧1。 图2算例1最优传输示意图 求得的最优解为x1*=[4 0 0 0 0 0 0 0 0 0 0 0 0] T,即只经过弧1运输4个单位流量， 其余弧无流量。又因为，每条弧的费用均为 5，所以最小费用为20。 经分析，计算结果合理可信。 算例 2( b2=[4;0;-4;0;0;0;0] T) 算例2中，由b2可知，节点3为需求节点，节点1为供给节点，由节点1将信息传 输至节点2的最短路径为弧2。 流量4 流量4 r 1 2 L 图3算例2最优传输示意图 求得的最优解为x2*=[0 4 0 0 0 0 0 0 0 0 0 0 0] T,即只经过弧2运输4个单位流量, 其余弧无流量。又因为，每条弧的费用均为 5，所以最小费用为20。 经分析，计算结果合理可信。 算例 3 (b3=[0;-4;4;0;0;0;0] T) 算例3中，由b3可知，节点2为需求节点，节点3为供给节点，由节点3将信息传 输至节点2的最短路径为弧5-＞弧1。 图4算例3最优传输示意图 求得的最优解为x3*=[4 0 0 0 4 0 0 0 0 0 0 0 0] T,即经过弧5运输4个单位流量至节 点1，再经弧1运输4个单位流量至节点2,其余弧无流量。又因为，每条弧的费用均为5， 所以最小费用为40。 经分析，计算结果合理可信。 算例 4( b4=[4;0;0;0;0;0;-4] T) 算例4中，由b4可知，节点7为需求节点，节点1为供给节点，由节点1将信息传 输至节点7的最短路径为弧1-＞弧4-＞弧10。 弧4 涂量4 图5算例3最优传输示意图 求得的最优解为x4*=[4 0 0 4 0 0 0 0 0 4 0 0 0] T,即经过弧1运输4个单位流量至节 点2，再经弧4运输4个单位流量至节点5,最后经弧