- 6
- 0
- 约9.62千字
- 约 62页
- 2018-02-27 发布于天津
- 举报
教学课件课件PPT医学培训课件教育资源教材讲义
或 将(1.20)式乘上一个正数θ>0得 (1.19)+(1.21)并整理得: 由(1.22)式找到满足约束方程组 的另一个点X(1),有 X(1)=(x10-qa1j,x20-qa2j,…,xm0-qamj,0,…,q,…,0)T 其中q是X(1)的第j个坐标的值。 要使X(1)是一个基本可行解,必须使 xi0-qaij≥0 1.23) 令这m个不等式中至少有一个等号成立。故可令 由式(1.24)知 因alj>0,故由矩阵元素组成的行列式不为零,P1,P2,…,Pl-1,Pl,Pl+1…,Pm是一个基。 在上述增广矩阵中作初等变换,将第l行乘上(1/alj),再分别乘以(-aij)(i=1,2,…,l-1,l+1,…,m)加到各行去,则增广矩阵左半部分变成单位矩阵。 所以,X(1)是一个可行解。 与变量xl1,…,x1l-1,,x1l+1,…,xm,xj对应的向量经重新排列后得 又因bl/alj=q,所以 b=(b1-qa1j,…,bl-1-qal-1,j,bl+1-qal+1,j,…,bm-qamj)T 由此X(1)是同X(0)相邻的基可行解,且由基向量组成的矩阵仍为单位矩阵。 3. 最优性检验和解的判别 将基本可行解X(0)和X(1)分别代入目标函数得 式中,因为q>0,所以只要 就有z(1)>z(0)。 最优性检验和解的判别准则 (1)当所有sj≤0时,当前基可行解是线性规划问题的最优解; (2)当所有sj≤0,若对某个非基变量xj有cj-zj=0,则该线性规划问题有无穷多个最优解;若对所有非基变量有sj<0,线性规划问题有唯一最优解。 (3) 若存在sj>0,又Pj≤0,则表明线性规划问题有无界解。 §1-4 单纯形法计算步骤 第一步:求初始基可行解,列出初始单纯形表。 cj→ c1 … cm … cj … cn CB 基 b x1 … xm … xj … xn c1 x1 b1 1 … 0 … a1j … a1n c2 x2 b2 0 … 0 … a2j … a2n ┆ ┆ ┆ ┆ ┆ ┆ … ┆ ┆ ┆ cm xm bm 0 … 1 … amj … amn cj-zj 0 … 0 … … 基变量及其值 问题中所有变量 单位矩阵 非基变量系数向量Pj表示为基向量线性组合时的系数,因基向量是单位向量,故有Pj=a1jP1+a2jP2+…+amj。 各变量在目标函数中的系数值 各基变量在目标函数中对应的系数 检验数sj=cj-zj=cj-(c1a1j+c2a2j+…+cmamj) 返回第一章目录 第二步:最优性检验 若单纯形表中所有检验数cj-zj≤0,且基变量中不含有人工变量,则得到线性规划问题的最优解,计算结束;若存在cj-zj>0,而Pj≤0,则问题为无界解,计算结束;否则转下一步。 第三步:从一个基可行解转换到相邻的目标函数值更大的基可行解,列出新的单纯形表。 1. 确定换入基的变量。只要有检验数sj>0,其对应的xj就可以作为换入基的变量,当有一个以上检验数大于零时,从中找出最大的一个sk,其对应的变量xk为换入基的变量(简称换入变量)。 2. 确定换出基的变量。用Pk列的系数分别去除常数项,找出最小比值 确定xl为换出基的变量,元素alk 决定了从一个基可行解到相邻基可行解的转移去向,所以,称alk为主元。 3. 用换入变量xk替换换出变量xl,作初等变换,得到一个新的单纯形表。 初等变换的方法是:将主元化为1,主元所在列的其它元素化为0。 第四步:重复第二、三两步,直到得出计算结果。 例5 用单纯形法求解线性规划问题 解:在约束方程中加松弛变量,将该线性规划问题化为标准形式 其约束条件系数矩阵的增广矩阵为: P1 P2 P3 P4 b P5 基变量为x3、x4、x5 ;非基变量为x1、x2。 单位矩阵 构成一个基 令非基变量x1、x2等于零,的初始基本可行解为:X(0)=(0,0,15,24,5)T 初始单纯形表 cj → 2 1 0 0 0 CB 基 b x1 x2 x3 x4 x5 0 x3 15 0 5 1 0 0 0 x4 24 6 2 0 1 0 0 x5 5 1 1 0 0 1 cj-zj 2 1 0 0 0 因为存在s1=2>0,s2=1>0。 所以初始基可行解不是最优解。 选择最大检验数对应的非基变量作为换入变量。 求最小比值,确定换入变量。 主元列 主元行 将主元化为1,主元所在列的其它元素化为0。 迭代运算 初始单纯形表 cj → 2 1 0 0 0 CB 基 b x1 x2 x3 x4 x5 0 x3 15 0 5 1 0 0 0 x4 24 [6] 2 0 1 0 0 x5 5 1 1 0 0 1 cj-zj 2 1 0 0 0 第一次迭代的
原创力文档

文档评论(0)