单纯形法的matlab实现极小化问题.docx

..

实验报告

实验题目:单纯形法的matlab实现学生XX:

学号:

实验时间:2021-4-15

一．实验名称:单纯形法的MATLAB实现二．实验目的及要求:

1.了解单纯形算法的原理及其matlab实现.

2.运用MATLAB编辑单纯形法程序解决线性规划的极小化问题,求出最优解及目标函

数值.

三．实验内容:

1.单纯形方法原理:

单纯形方法的根本思想,是从一个根本可行解出发,求一个使目标函数值有所改善的根本可行解;通过不断改进根本可行解,力图到达最优根本可行解.

对问题

其中A是一个m×n矩阵,且秩为m,c为n维行向量,x为n维列向量,b为m维非

负列向量.符号

负列向量.符号〞表示右端的表达式是左端的定义式,即目标函数f的具体形式就是

cx.

记

令A=(B,N),B为基矩阵,N为非基矩阵,设

是根本可行解,在x(0)处的目标函数值

.v.

..

.v.

f0=cx(0)=(cB,cN)Bb=cBB-1b,

其中cB是c中与基变量对应的分量组成的m维行向量;cN是c中与非基变量对应的分量组

成的n-m维行向量.

现由根本可行解x(0)出发求解一个改进的根本可行解.

LxN」设x=「|xB]|是任一可行解,那么由Ax

LxN」

BN,x=B-1b-B-

BN,

在点x处的目标函数值

LxN」j=Rf=cx=(cB,cN)「|xB]|=f0-x(zj-

LxN」j=R

其中R是非基变量下标集,

zj=cBB-1pj.

2.单纯形方法计算步骤:

首先给定一个初始根本可行解,设初始基为B,然后执行以下主要步骤:

(1)解BxB=b,求得xB=B-1b=b_,令xN=0,计算目标函数值f=cBxB.

(2)求单纯形乘子w,解wB=cB,得到w=cBB-1.对于所有非基变量,计算判别数

zj-cj=wjpj-cj.令

zk-ck=x{zj-cj}.

假设zk-ck不0,那么对于所有非基变量zj-cj不0,对应基变量的判别数总是为零,因此停顿计算,现行根本可行解是最优解.否那么,进展下一步.

(3)解Byk=pk,得到yk=B-1pk,假设yk不0,即yk的每个分量均非正数,那么停顿计算,问题不存在有限最优解.否那么进展步骤〔4〕.

(4)确定下标r,使

..

xk=