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

《利用修正单纯形法解线性规划问题(精).doc》由会员分享，可在线阅读，更多相关《利用修正单纯形法解线性规划问题(精).doc（7页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、精品文档，仅供学习与交流，如有侵权请联系网站删除利用修正单纯形法解线性规划问题一 软件示意：二 代码说明：Dim A(1 To 3, 1 To 6) As Double 矩阵ADim a1(1 To 3) As Double 矩阵A的第一列向量Dim a2(1 To 3) As Double 矩阵A的第二列向量Dim a3(1 To 3) As Double 矩阵A的第三列向量Dim a4(1 To 3) As Double 矩阵A的第四列向量Dim a5(1 To 3) As Double 矩阵A的第五列向量Dim a6(1 To 3) As Double 矩阵A的第六列向量Dim B_(1

2、 To 3, 1 To 3) As Double 基矩阵B的逆矩阵Dim XB(1 To 3) As Double 基本可行解Dim b(1 To 3) As Double 右端向量bDim C(1 To 6) As Double 检验数Dim CB(1 To 3) As Double 基本可行解对应的检验数Dim (1 To 3) As Double 单纯形乘子矢量Dim r(1 To 6) As Double 检验矢量rDim r_min As Double 检验矢量最小值Dim k_sign As Integer 检验矢量最小值对应的位置Dim Y(1 To 3, 0 To 6) As

3、Double 矩阵yDim just_vector(1 To 3) As DoubleDim liji_min As Double 用于判断离基变量所用值Dim r_sign As Integer 用于记录离基变量对应的位置Dim main_yuan As Double 用于存放主元Dim Erk(1 To 3, 1 To 3) As DoubleDim Exchange_B(1 To 3, 1 To 3) 在矩阵Erk与矩阵B_进行乘法运算时，作为矩阵B_的替换矩阵Dim Exchange_XB(1 To 3) 在矩阵Erk与XB_进行乘法运算时，作为XB_的替换矩阵Dim iterativ

4、e_time As Integer 定义迭代的次数Dim XB_optimization(1 To 6) As Double 最优解Private Sub Command1_Click()窗口1For iterative_time = 1 To 1000 开始了迭代循环 Select Case k_sign Case 1 If (r_sign - 3) = 3 Then CB(1) = C(2) 付值给基本可行解对应的检验数 CB(2) = C(3) CB(3) = C(6) XB_optimization(1) = 0 XB_optimization(2) = XB(1) XB_optimi

5、zation(3) = XB(2) XB_optimization(4) = 0 XB_optimization(5) = 0 XB_optimization(6) = XB(3) End If Case 2 If (r_sign - 3) = 1 Then CB(1) = C(2) 付值给基本可行解对应的检验数 CB(2) = C(5) CB(3) = C(6) XB_optimization(1) = 0 XB_optimization(2) = XB(1) XB_optimization(3) = 0 XB_optimization(4) = 0 XB_optimization(5) =

6、 XB(2) XB_optimization(6) = XB(3) End If Case 3 If (r_sign - 3) = 2 Then CB(1) = C(2) 付值给基本可行解对应的检验数 CB(2) = C(3) CB(3) = C(6) XB_optimization(1) = 0 XB_optimization(2) = XB(1) XB_optimization(3) = XB(2) XB_optimization(4) = 0 XB_optimization(5) = 0 XB_optimization(6) = XB(3) End If End Select For i

7、 = 1 To 3 计算单纯形乘子矢量 (i) = CB(1) * B_(1, i) + CB(2) * B_(2, i) + CB(3) * B_(3, i) Next i For j = 1 To 6 计算检验矢量r r(j) = C(j) - (1) * A(1, j) + (2) * A(2, j) + (3) * A(3, j) Next j r_min = r(1) 预先给定一个值 k_sign = 1 For i = 1 To 6 找出最小检验值 If (r(i) r_min) Then r_min = r(i) 最小检验值 k_sign = i 最小检验值所在位置,对应进基矢量

8、位置 End If Next i If r_min 0 Then For i = 1 To 3 对应进基矢量所得的y值 Y(i, k_sign) = B_(i, 1) * A(1, k_sign) + B_(i, 2) * A(2, k_sign) + B_(i, 3) * A(3, k_sign) Next i For i = 1 To 3 确定离基变量 If Y(i, k_sign) 0 Then just_vector(i) = Y(i, 0) / Y(i, k_sign) End If Next i For i = 1 To 3 If (just_vector(i) 0 Then li

9、ji_min = just_vector(i) 预先给定一个值 Exit For End If Next i r_sign = 4 For i = 1 To 3 找出最小检验值 If (just_vector(i) 0) Then liji_min = just_vector(i) 最小检验值 r_sign = i + 3 最小检验值所在位置,对应进基矢量位置 End If Next i main_yuan = Y(r_sign - 3, k_sign) If main_yuan 0 Then Select Case (r_sign - 3) 计算矩阵Erk Case 1 Erk(1, 1)

10、= 1 / main_yuan Erk(2, 1) = -Y(2, k_sign) / main_yuan Erk(3, 1) = -Y(3, k_sign) / main_yuan For i = 2 To 3 For j = 1 To 3 Erk(j, i) = A(j, i + 3) Next j Next i Case 2 Erk(1, 2) = -Y(1, k_sign) / main_yuan Erk(2, 2) = 1 / main_yuan Erk(3, 2) = -Y(3, k_sign) / main_yuan For j = 1 To 3 Erk(j, 1) = A(j,

11、 4) Next j For j = 1 To 3 Erk(j, 3) = A(j, 6) Next j Case 3 Erk(1, 3) = -Y(1, k_sign) / main_yuan Erk(2, 3) = -Y(2, k_sign) / main_yuan Erk(3, 3) = 1 / main_yuan For i = 1 To 2 For j = 1 To 3 Erk(j, i) = A(j, i + 4) Next j Next i End Select End If For i = 1 To 3 给矩阵Exchange_B(i, j)付值 For j = 1 To 3

12、Exchange_B(i, j) = B_(i, j) Next j Next i For i = 1 To 3 计算B的新可逆矩阵B_ For j = 1 To 3 B_(i, j) = Erk(i, 1) * Exchange_B(1, j) + Erk(i, 2) * Exchange_B(2, j) + Erk(i, 3) * Exchange_B(3, j) Next j Next i For i = 1 To 3 给Exchange_XB(i)付值 Exchange_XB(i) = XB(i) Next i For i = 1 To 3 计算基本可行解 XB(i) = Erk(i,

13、 1) * Exchange_XB(1) + Erk(i, 2) * Exchange_XB(2) + Erk(i, 3) * Exchange_XB(3) Y(i, 0) = XB(i) Next i Else MsgBox (优化完毕!) Exit For End IfNext iterative_timeCommand2.Enabled = True 使能按钮Command3.Enabled = True 使能按钮Form1.Text1(1).Text = XB_optimization(1) 显示最优解x1Form1.Text1(2).Text = XB_optimization(2)

14、 显示最优解x2Form1.Text1(3).Text = XB_optimization(3) 显示最优解x3Form1.Text1(4).Text = XB_optimization(4) 显示最优解x4Form1.Text1(5).Text = XB_optimization(5) 显示最优解x5Form1.Text1(6).Text = XB_optimization(6) 显示最优解x6Form1.Text1(0).Text = C(1) * XB_optimization(1) + C(2) * XB_optimization(2) + C(3) * XB_optimization

15、(3) + C(4) * XB_optimization(4) + C(5) * XB_optimization(5) + C(6) * XB_optimization(6) 显示最优值ZForm1.Text1(7) = iterative_time - 1 显示迭代次数End SubPrivate Sub Command2_Click()End 结束程序End SubPrivate Sub Command3_Click() Form1.Show 调用优化结果显示窗口End SubPrivate Sub Form_Load()Dim i, j As IntegerFor i = 1 To 3

16、付值给系数矩阵A For j = 1 To 6 A(i, j) = Quotiety_A(j + 6 * (i - 1).Text 付值给系数矩阵A Next jNext iFor i = 1 To 6 将矩阵A列分块（a1,a2,a3,a4,a5,a6） For j = 1 To 3 Select Case i Case 1 a1(j) = A(j, i) Case 2 a2(j) = A(j, i) Case 3 a3(j) = A(j, i) Case 4 a4(j) = A(j, i) Case 5 a5(j) = A(j, i) Case 6 a6(j) = A(j, i) End

17、Select Next jNext iFor i = 1 To 3 付值给基矩阵B的逆矩阵B_ For j = 1 To 3 Select Case i Case 1 B_(j, i) = a4(j) Case 2 B_(j, i) = a5(j) Case 3 B_(j, i) = a6(j) End Select Next jNext iFor i = 1 To 6 付值给检验数 Select Case i Case 1 C(1) = C_quotiety(0).Text Case 2 C(2) = C_quotiety(1).Text Case 3 C(3) = C_quotiety(2

18、).Text Case 4 C(4) = C_quotiety(3).Text Case 5 C(5) = C_quotiety(4).Text Case 6 C(6) = C_quotiety(5).Text End SelectNext ib(1) = b_quotiety(0).Text 付值给右端向量bb(2) = b_quotiety(1).Textb(3) = b_quotiety(2).TextFor i = 1 To 3 计算初始基本可行解 XB(i) = B_(i, 1) * b(1) + B_(i, 2) * b(2) + B_(i, 3) * b(3) Y(i, 0) = XB(i)Next iCB(1) = C(4) 付值给初始基本可行解对应的检验数CB(2) = C(5)CB(3) = C(6)End SubPrivate Sub Command1_Click()窗口2 Unload Form1 关闭优化结果显示窗口End Sub【精品文档】第 7 页