利用修正单纯形法解线性规划问题.doc

《利用修正单纯形法解线性规划问题.doc》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、利用修正单纯形法解线性规划问题一软件示意：二代码说明：DimA(1To3,1To6)AsDouble'矩阵ADima1(1To3)AsDouble'矩阵A的第一列向量Dima2(1To3)AsDouble'矩阵A的第二列向量Dima3(1To3)AsDouble'矩阵A的第三列向量Dima4(1To3)AsDouble'矩阵A的第四列向量Dima5(1To3)AsDouble'矩阵A的第五列向量Dima6(1To3)AsDouble'矩阵A的第六列向量DimB_(1To3,1To3)AsDouble'基矩阵B的逆矩阵DimXB(1To3)AsDoub

2、le'基本可行解Dimb(1To3)AsDouble'右端向量bDimC(1To6)AsDouble'检验数DimCB(1To3)AsDouble'基本可行解对应的检验数Dimπ(1To3)AsDouble'单纯形乘子矢量Dimr(1To6)AsDouble'检验矢量rDimr_minAsDouble'检验矢量最小值Dimk_signAsInteger'检验矢量最小值对应的位置DimY(1To3,0To6)AsDouble'矩阵yDimjust_vector(1To3)AsDoubleDimliji_minAsDouble'用于判断离基变量所用值Di

3、mr_signAsInteger'用于记录离基变量对应的位置Dimmain_yuanAsDouble'用于存放主元DimErk(1To3,1To3)AsDoubleDimExchange_B(1To3,1To3)'在矩阵Erk与矩阵B_进行乘法运算时，作为矩阵B_的替换矩阵DimExchange_XB(1To3)'在矩阵Erk与XB_进行乘法运算时，作为XB_的替换矩阵Dimiterative_timeAsInteger'定义迭代的次数DimXB_optimization(1To6)AsDouble'最优解PrivateSubCommand1_Cli

4、ck()'窗口1Foriterative_time=1To1000'开始了迭代循环SelectCasek_signCase1If(r_sign-3)=3ThenCB(1)=C(2)'付值给基本可行解对应的检验数CB(2)=C(3)CB(3)=C(6)XB_optimization(1)=0XB_optimization(2)=XB(1)XB_optimization(3)=XB(2)XB_optimization(4)=0XB_optimization(5)=0XB_optimization(6)=XB(3)EndIfCase2If(r_sign-3

5、)=1ThenCB(1)=C(2)'付值给基本可行解对应的检验数CB(2)=C(5)CB(3)=C(6)XB_optimization(1)=0XB_optimization(2)=XB(1)XB_optimization(3)=0XB_optimization(4)=0XB_optimization(5)=XB(2)XB_optimization(6)=XB(3)EndIfCase3If(r_sign-3)=2ThenCB(1)=C(2)'付值给基本可行解对应的检验数CB(2)=C(3)CB(3)=C(6)XB_optimization(1)=0X

6、B_optimization(2)=XB(1)XB_optimization(3)=XB(2)XB_optimization(4)=0XB_optimization(5)=0XB_optimization(6)=XB(3)EndIfEndSelectFori=1To3'计算单纯形乘子矢量π(i)=CB(1)*B_(1,i)+CB(2)*B_(2,i)+CB(3)*B_(3,i)NextiForj=1To6'计算检验矢量rr(j)=C(j)-(π(1)*A(1,j)+π(2)*A(2,j)+π(3)*A(3,j))Nextjr_min=r(1)'预先给

7、定一个值k_sign=1Fori=1To6'找出最小检验值If(r(i)0Thenjust_vector(i)=Y(i,0)/Y(i,k_sign

8、)EndIfNextiFori=1To3If(just_vector(i))>0Thenliji_min=j