【运筹学】线性规划 单纯形法 ( 基矩阵 | 基变量 | 非基矩阵 | 非基变量 | 矩阵分块形式 | 逆矩阵 | 基解 | 基可行解 )
文章目录
- I . 基矩阵 B
- II . 基向量 PjP_jPj
- III . 基变量
- IV . 非基矩阵 NNN
- V . 系数矩阵分块形式 A=(BN)A = ( B N )A=(BN)
- VI . 基变量向量 XBX_BXB 非基变量向量 XNX_NXN 及 分块形式
- VII . 分块形式的计算公式
- VIII . 逆矩阵
- IX . 解基变量
- X . 基解
- XI . 基可行解
I . 基矩阵 B
线性规划标准形式 , 约束方程的系数矩阵是 AAA , 如下 :
A=[a11a12⋯a1na21a22⋯a2n⋮⋮⋱⋮am1am2⋯amn]A = \begin{bmatrix}\\\\ & a_{11} & a_{12} & \cdots & a_{1n} &\\\\ & a_{21} & a_{22} & \cdots & a_{2n} &\\\\ & \vdots & \vdots & \ddots & \vdots &\\\\ & a_{m1} & a_{m2} & \cdots & a_{mn} &\\\\ \end{bmatrix}A=⎣⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢⎡a11a21⋮am1a12a22⋮am2⋯⋯⋱⋯a1na2n⋮amn⎦⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎤
该矩阵 AAA 是 m×nm \times nm×n 阶矩阵 , 有 mmm 行 , nnn 列 , 代表 mmm 个约束方程 , nnn 个变量 , 并且 n>mn > mn>m ;
基矩阵 BBB :
- ① 满秩子矩阵 : 矩阵 AAA 的 满秩子矩阵 BBB , 矩阵 BBB 的秩是 mmm ;
- ② 列向量线性无关 : 该矩阵中的 列向量 线性无关 , 即 每一列不能通过 乘以系数 加减的方式得到另外一列列向量 ,
- ③ 基矩阵 BBB : 这样的 系数矩阵 AAA 的 m×mm \times mm×m 阶满秩矩阵 BBB 就是基矩阵 ;
B=[a11a12⋯a1ma21a22⋯a2m⋮⋮⋱⋮am1am2⋯amm]=[P1P2⋯Pm]B= \begin{bmatrix}\\\\ & a_{11} & a_{12} & \cdots & a_{1m} &\\\\ & a_{21} & a_{22} & \cdots & a_{2m} &\\\\ & \vdots & \vdots & \ddots & \vdots &\\\\ & a_{m1} & a_{m2} & \cdots & a_{mm} &\\\\ \end{bmatrix} = \begin{bmatrix} & P_1 & P_2 & \cdots & P_m & \end{bmatrix}B=⎣⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢⎡a11a21⋮am1a12a22⋮am2⋯⋯⋱⋯a1ma2m⋮amm⎦⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎤=[P1P2⋯Pm]
II . 基向量 PjP_jPj
基向量 :
- ① 概念 : 基矩阵 BBB 中的每个列向量 , 都是一个 基向量 , 记作 PjP_jPj , 其中 j=1,2,⋯,mj = 1 , 2 , \cdots , mj=1,2,⋯,m ;
- ② 基向量个数 : 每个基矩阵中有 mmm 个列向量 ;
III . 基变量
基变量 : 每个基向量都对应一个变量 , 基向量是列向量 , 该列向量是 xjx_jxj 变量的系数组成 , 这个对应的 xjx_jxj 变量就是基变量 ;
IV . 非基矩阵 NNN
非基矩阵 NNN : 确定一个基矩阵 , 剩下的列向量就是 非基向量 , 这些非基向量 组成 非基矩阵 NNN ;
N=[a1m+1a1m+2⋯a1na2m+1a2m+2⋯a2n⋮⋮⋱⋮amm+1amm+2⋯amn]=[Pm+1Pm+2⋯Pn]N= \begin{bmatrix}\\\\ & a_{1m+1} & a_{1m+2} & \cdots & a_{1n} &\\\\ & a_{2m+1} & a_{2m+2} & \cdots & a_{2n} &\\\\ & \vdots & \vdots & \ddots & \vdots &\\\\ & a_{mm+1} & a_{mm+2} & \cdots & a_{mn} &\\\\ \end{bmatrix} = \begin{bmatrix} & P_{m+1} & P_{m+2} & \cdots & P_{n} & \end{bmatrix}N=⎣⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢⎡a1m+1a2m+1⋮amm+1a1m+2a2m+2⋮amm+2⋯⋯⋱⋯a1na2n⋮amn⎦⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎤=[Pm+1Pm+2⋯Pn]
V . 系数矩阵分块形式 A=(BN)A = ( B N )A=(BN)
系数矩阵 AAA , 可以写成分块形式 :
A=[a11a12⋯a1na21a22⋯a2n⋮⋮⋱⋮am1am2⋯amn]=[BN]A = \begin{bmatrix}\\\\ & a_{11} & a_{12} & \cdots & a_{1n} &\\\\ & a_{21} & a_{22} & \cdots & a_{2n} &\\\\ & \vdots & \vdots & \ddots & \vdots &\\\\ & a_{m1} & a_{m2} & \cdots & a_{mn} &\\\\ \end{bmatrix} = \begin{bmatrix} & B & N & \end{bmatrix}A=⎣⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢⎡a11a21⋮am1a12a22⋮am2⋯⋯⋱⋯a1na2n⋮amn⎦⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎤=[BN]
VI . 基变量向量 XBX_BXB 非基变量向量 XNX_NXN 及 分块形式
基变量向量 XBX_BXB :
XB=[x1x2⋮xm]X_B = \begin{bmatrix}\\\\ & x_1 &\\\\ &x_2&\\\\ &\vdots& \\\\ &x_m&\\\\ \end{bmatrix}XB=⎣⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢⎡x1x2⋮xm⎦⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎤
非基变量向量 XNX_NXN :
XB=[xm+1xm+2⋮xn]X_B = \begin{bmatrix}\\\\ & x_{m + 1} &\\\\ &x_{m + 2}&\\\\ &\vdots& \\\\ &x_n&\\\\ \end{bmatrix}XB=⎣⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢⎡xm+1xm+2⋮xn⎦⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎤
向量 XXX 可以写成 XBX_BXB 和 XNX_NXN 分块形式 :
X=[x1x2⋮xn]=[xBxN]X = \begin{bmatrix}\\\\ & x_1 &\\\\ &x_2&\\\\ &\vdots& \\\\ &x_n&\\\\ \end{bmatrix} = \begin{bmatrix}\\\\ & x_B &\\\\ &x_N &\\\\ \end{bmatrix}X=⎣⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢⎡x1x2⋮xn⎦⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎤=⎣⎢⎢⎢⎢⎡xBxN⎦⎥⎥⎥⎥⎤
VII . 分块形式的计算公式
矩阵分块形式方程代入 : 约束方程组 AX=bAX = bAX=b ;
bbb 是大于 000 的常数组成的向量 ;
将上述分块形式的 矩阵 AAA 和 矩阵 XXX 代入 上述 AX=bAX = bAX=b 公式 ;
A=[BN]A = \begin{bmatrix} & B & N & \end{bmatrix}A=[BN]
X=[XBXN]X = \begin{bmatrix}\\\\ & X_B &\\\\ &X_N &\\\\ \end{bmatrix}X=⎣⎢⎢⎢⎢⎡XBXN⎦⎥⎥⎥⎥⎤
得到
[BN]×[XBXN]=b\begin{bmatrix} & B & N & \end{bmatrix} \times \begin{bmatrix}\\\\ & X_B &\\\\ & X_N &\\\\ \end{bmatrix} = b[BN]×⎣⎢⎢⎢⎢⎡XBXN⎦⎥⎥⎥⎥⎤=b
BXB+NXN=bBX_B + NX_N = bBXB+NXN=b
VIII . 逆矩阵
逆矩阵 : 其中矩阵 BBB 是满秩的 m×mm \times mm×m 阶矩阵 , 该矩阵是可逆的 ( 非奇异矩阵 ) , 必定存在一个 B−1B^{-1}B−1 , 使得
B×B−1=EB \times B^{-1} = EB×B−1=E
单位矩阵 : 这里的 矩阵 EEE 是单位矩阵 , 即 左上角到右下角 对角线 上 的元素 为 111 , 其它元素为 000 ;
主对角线 : 左上角 到 右下角 的对角线称为 主对角线 ;
单位矩阵 示例 如下 :
E=[100010001]E=\begin{bmatrix} & 1 & 0 & 0 & \\\\ & 0 & 1 & 0 &\\\\ & 0 & 0 & 1 & \end{bmatrix}E=⎣⎢⎢⎢⎢⎡100010001⎦⎥⎥⎥⎥⎤
IX . 解基变量
解基变量 :
BXB+NXN=bBX_B + NX_N = bBXB+NXN=b
将 NXNNX_NNXN 提到公式右边 :
BXB=b−NXNBX_B = b - NX_NBXB=b−NXN
公式两边乘以 B−1B^{-1}B−1 :
BXB×B−1=(b−NXN)×B−1BX_B \times B^{-1} = ( b - NX_N ) \times B^{-1}BXB×B−1=(b−NXN)×B−1
XB=B−1b−B−1NXNX_B = B^{-1}b - B^{-1}NX_NXB=B−1b−B−1NXN
X . 基解
引入基解 : 令非基变量 XNX_NXN 中所有变量为 000 , 此时上述公式为 :
XB=B−1bX_B = B^{-1}bXB=B−1b
约束方程的解为
X=[XB0]=[B−1b0]X = \begin{bmatrix} & X_B & \\\\ & 0 & \end{bmatrix}=\begin{bmatrix} & B^{-1}b & \\\\ & 0 & \end{bmatrix}X=⎣⎡XB0⎦⎤=⎣⎡B−1b0⎦⎤
上述解为基解 , 矩阵 BBB 是满秩的 , 其秩为 mmm , 将非基变量赋值 000 , 剩余的 mmm 个变量 , mmm 个等式 , 必能解出一组唯一解 ; 即
∑j=1mpjxj=b\sum_{j = 1}^{m}p_j x_j = bj=1∑mpjxj=b
方程组有唯一解
XB=[x1x2⋮xm]X_B = \begin{bmatrix} & x_1 & \\\\ & x_2 &\\\\ & \vdots &\\\\ & x_m & \end{bmatrix}XB=⎣⎢⎢⎢⎢⎢⎢⎢⎢⎢⎡x1x2⋮xm⎦⎥⎥⎥⎥⎥⎥⎥⎥⎥⎤
该解 XBX_BXB 是线性规划的一个基解 ;
XI . 基可行解
基可行解 : 如果上述解出的基解 XBX_BXB , 满足线性规划数学模型 标准形式 的变量非负约束 , 即所有的变量都大于等于 000 , 该解称为基可行解 ;
并不是所有的基解都是基可行解 , 只有部分基解是基可行解 ;
总结
以上是生活随笔为你收集整理的【运筹学】线性规划 单纯形法 ( 基矩阵 | 基变量 | 非基矩阵 | 非基变量 | 矩阵分块形式 | 逆矩阵 | 基解 | 基可行解 )的全部内容,希望文章能够帮你解决所遇到的问题。
- 上一篇: 【运筹学】线性规划 单纯形法 ( 原理
- 下一篇: 【Kotlin】IntelliJ IDE