Giải thuật đơn hình

Giải thuật đơn hình

CHƯƠNG II

GIẢI THUẬT ĐƠN HÌNH

Chương này trình bày một cách chi tiết nội dung của giải thuật đơn hình Sau phần cơ sở lý thuyết của giải thuật là các ví dụ tương ứng Các ví dụ được trình bày đúng theo các bước của giải thuật Kiến thức trong chương này cần thiết cho việc lập trình giải quy hoạch tuyến tính trên máy tính

Nội dung chi tiết của chương bao gồm :

I- GIẢI THUẬT ĐƠN HÌNH CƠ BẢN

1- Cơ sở xây dựng giải thuật đơn hình cơ bản

2- Định lý về sự hội tụ

3- Giải thuật đơn hình cơ bản

4- Chú ý trong trường hợp suy biến

II- GIẢI THUẬT ĐƠN HÌNH CẢI TIẾN

1- Một cách tính ma trận nghịch đảo

2- Quy hoạch tuyến tính dạng chuẩn

3- Giải thuật đơn hình cải tiến

4- Phép tính trên dòng - Bảng đơn hình

III- PHƯƠNG PHÁP BIẾN GIẢ CẢI BIÊN

1- Bài toán cải biên

a- Cải biên bài toán quy hoạch tuyến tính

b- Quan hệ giữa bài toán xuất phát và bài toán cải biên

2- Phương pháp hai pha

3- Phương pháp M vô cùng lớn

IV- QUY HOẠCH TUYẾN TÍNH SUY BIẾN

1- Các ví dụ về quy hoạch tuyến tính suy biến

2- Xử lý quy hoạch tuyến tính suy biến

CHƯƠNG II: GIẢI THUẬT ĐƠN HÌNH

I- GIẢI THUẬT ĐƠN HÌNH CƠ BẢN

Chương này trình bày một phương pháp để giải bài toán quy hoạch tuyến tính

đó là phương pháp đơn hình Phương pháp đơn hình được George Bernard Dantzig

đưa ra năm 1947 cùng lúc với việc ông khai sinh ra quy hoạch tuyến tính Đây là một

phương pháp thực sự có hiệu quả để giải những bài toán quy hoạch tuyến tính cở lớn

trong thực tế Với cách nhìn hiện đại ý tưởng của phương pháp đơn hình rất đơn giản

Có nhiều cách tiếp cận phương pháp đơn hình, chương này trình bày một trong các

cách đó

1- Cơ sở xây dựng giải thuật đơn hình cơ bản

Xét bài toán quy hoạch tuyến tính chính tắc :

bAx

xcz(x)

Giả sử rằng B0 là một cơ sở khả thi xuất phát của bài toán ( không nhất thiết là

m cột đầu tiên của ma trận A ) Thuật toán đơn hình cơ bản được xây dựng dựa trên

bBxxN

1

B : phương án cơ sở khả thi tương ứng

bB

b= − 1

NBcc

N

T N

i s

a

b0a ,a

bmin

Về mặt hình học, giải thuật này được hiểu như là một quá trình duyệt qua các điểm cực biên của đa diện lồi S các phương án khả thi của bài toán

Về mặt đại số, giải thuật này được hiểu như là một quá trình xác định một chuỗi các ma trận cơ sở kề B0 B1 B2 mà các phương án cơ sở tương ứng x0 x1

3- Giải thuật đơn hình cơ bản

Xét bài toán quy hoạch tuyến tính chính tắc

0x

bAx

x

c)x(z

Giả sử rằng sau khi hoán vị các cột trong A ta chọn được ma trận cơ sở B thoả

sự phân hoạch sau đây :

A = [ B N ]

]c c[

cT = B N

]x x[

xT = B N Giải thuật đơn hình cơ bản được thực hiện như sau :

bbBxx

N

1 B

Giá trị hàm mục tiêu z(x) =cBTxB Ma trận = B N -1N

c- Xét dấu hiệu tối ưu :

T B

T N 1 T B

T N

bbBxx

N

1 B

và giá trị hàm mục tiêu là :

B

T

Bxc)x(

- Nếu tồn tại cs ∈cN mà cs > 0thì sang bước d

d- Xác định chỉ số của phần tử pivot trong ma trận N

Xác định chỉ số cột s của pivot

cs =max {ck >0∈cN}

Nếu Nis ≤0 thì giải thuật dừng, bài toán không có phương án tối ưu Ngược lại thì tiếp tục

Xác định chỉ số dòng r của pivot

m)1,2, ,(i

N

b0N ,N

bmin

rs

r is

Phần tử Nrs trong ma trận được gọi là phần tử pivot N

Trong trường hợp bài toán min

c- Xét dấu hiệu tối ưu :

T B

T N 1 T B

T N

bbBxx

N

1 B

B

T

Bxc)x(

- Nếu tồn tại cs ∈cN mà cs <0thì sang bước d

d- Xác định chỉ số của phần tử pivot trong ma trận N

N

b0N ,N

bmin

rs

r is

Biến xs trong xNT với biến xr trong xBT

f- Quay về (a)

Ví dụ : Tìm phương án tối ưu cho bài toán quy hoạch tuyến tính chính tắc sau đây bằng giải thuật đơn hình cơ bản

−

=++

=+

−

+

=

1,2,3,4,5)(j

0x

2xx2x

6xx2x

3xxx

xx2)x(z max

j

5 2 1

4 2 1

3 2 1

2 1

Ta có :

T B

T N T

T B

T N

5 4 3 2

1 T

c c

0 0 0

| 1 2 c

x x

xxx

|xxx

B

N

26

3b 1

00

|2 1

0 10

|2 1

0 0 1

|11A

010

001B

xx

b263

263

1 0 0

0 1 0

0 0 1bBxx

xxx

2

1 N

1

5 4

3 B

Giá trị hàm mục tiêu :

[ ] 0

26

3000xc)x(

2 1

11

21

2 1

11

100

010

001NB

21

1100012Ncc

B

T N

1

1

N1

Xác định chỉ số dòng pivot r :

11

1 21

2 11

1 is

i

N

b31

6,1

3minN

b,N

bminN

Cột thứ s=1 trong ma trận N và cột thứ r=1 trong ma trận B Phần tử thứ s=1 trong cNT với phần tử thứ r=1 trong T

Bc

Biến thứ s=1 trong xNT với biến thứ r=1 trong T

Bx

|20

011

|20

001

|11A100

|21

010

|21

001

|11A

5 4 3 2

1

001B

1

01

011

001

xx

b533

263

1 0 1

0 1 1

0 0 1bBxx

xxx

2

3 N

1

5 4

1 B

Giá trị hàm mục tiêu :

53

3002xc)x(

31-

11

20

2 0

11

101

011-

001NB

31-

110 0 210Ncc

B

T N

=

Vậy s=2

3

1-

N2

Xác định chỉ số dòng pivot r :

22

2 23

3 22

2 is

i

N

b11

5,3

3minN

b,N

bminN

Cột thứ s=2 trong ma trận N và cột thứ r=2 trong ma trận B Phần tử thứ s=2 trong cNT với phần tử thứ r=2 trong T

Bc

Biến thứ s=2 trong xNT với biến thứ r=2 trong T

Bx

121

|00

021

|10

011

|01A101

|20

011

|20

001

|11A

1- 3

4

0 3

1 31

0 3

1 3

2

B 1

21

021

01-1

b- Tính các tham số

Phương án cơ sở khả thi tốt hơn :

b414

263

1 3

1- 3

4

0 3

1 31

0 3

1 3

2

bBxx

xxx

4

3 N

1

5 2

1 B

Giá trị hàm mục tiêu :

41

4012xc)x(

1 313

1 32

00

10

01

1 3

1- 3

4

0 3

1 31

0 3

1 3

2

NB

1 313

1 32

0 1 200Ncc

B

T N

xx

414x

x

xx

4

3 N

5 2

1 B

Giá trị hàm mục tiêu là z(x) = 9 với x1 = 4 và x2 = 1

4- Chú ý trong trường hợp suy biến

Trong trường hợp bài toán suy biến, nghĩa là br =0, ta có :

0a

bxrs

z ∧ = + s ∧s =

Vậy thì, có thể sau một số lần thay đổi cơ sở lại quay trở về cơ sở đã gặp và lặp như vậy một cách vô hạn Người ta có nhiều cách để khắc phục hiện tượng này bằng cách xáo trộn một chút các dữ liệu của bài toán, sử dụng thủ tục từ vựng, quy tắc

B

−

∧

-1 Để làm điều đó chỉ cần nhân (bên trái) B-1 với một ma trận đổi cơ sở được xác định như sau :

rcôt

r dòng

1

a

a

00

a

1

00

a

a

10

0

a

a

01

rs ms rs

rs 2s rs 1s

Khi đó :

1 1

isa

a

− : đối với thành phần i ≠ r

rsa

1 : đối với thành phần r Khi mà ma trận cở sở xuất phát là ma trận đơn vị, sau một số bước đổi cơ sở

B0 B1 B2 Bq tương ứng với các ma trận đổi cơ sở µ0 µ1 µ2 … µq-1 người ta có cách tính ma trận nghịch đảo như sau :

[ ]Bq −1 =µ0.µ1 µq − 1

2- Quy hoạch tuyến tính dạng chuẩn

Quy hoạch tuyến tính dạng chuẩn là quy hoạch tuyến tính chính tắc mà trong

đó có thể rút ra một ma trận cơ sở là ma trận đơn vị Quy hoạch tuyến tính chuẩn có dạng :

bxN]

I

xc)x(z max

3- Giải thuật đơn hình cải tiến

Từ những kết quả trên người ta xây dựng giải thuật đơn hình cải tiến đối với bài toán qui hoạch tuyến tính (max) dạng chuẩn như sau :

a- Khởi tạo

A

A0 =b

b0 =

b- Thực hiện bước lặp với k = 0,1,2,

Xác định phương án cơ sở khả thi :

bxx

k

k

N

k B k

Tính giá trị hàm mục tiêu :

k

T B B

T B

k) c x c bx

k

k

N

k B

k là phương án tối ưu

k

T B B

T B

k) c x c bx

.Tính bk + 1 = µkbk.Tăng số lần lặp k=k+1

0x

2x2xx

6x2xx

3xxx

x2xz(x)max

j

5 2 1

4 2 1

3 2 1

2 1

−

=++

=+

0 0

B N

26

3b 1

00

|21

010

|21

001

|11AA

T N T

0

0 cc

000

|12

26

3bxx

xx

3 B

0

[ ] 0

26

30 0 0bc)

0 1 0 2 1

0 0 1 1- 1 0 0 00 0 0 1 2Ac

011

001

011

001

0 0 1 1- 1

0 1 1- 3 0

0 0 1 1- 1

011

001

53

3bxx

xx

1 B

1

53

30 0 2bc)

cc

0 1 1- 3 0

0 0 1 1- 1

= [ 0 3 -2 0 0 ]

1-

10

03

1 0

03

1 11

=µ

03

1 0

03

1 1

0 1 1- 3 0

0 0 1 1- 1

1- 3

4 0 0

0 3

1 3

1- 1 0

0 3

1 3

2 0 1

=µ

03

1 0

03

1 1

41

4bxx

xx

1 B

2

41

40 1 2bc)

cc

1- 3

4 0 0

0 3

1 3

1- 1 0

0 3

1 3

2 0 1

= [ 0 0 -1 -1 0 ] : thoả dấu hiệu tối ưu

Vậy kết quả của bài toán là :

Phương án tối ưu x = x2 =

Giá trị hàm mục tiêu z(x) = 9

4- Phép tính trên dòng - Bảng đơn hình

Các bước thực hiện giải thuật đơn hình cải tiến được trình bày lần lượt trong các bảng, gọi là bảng đơn hình Trong thực hành, để cập nhật những giá trị mới ta có thể làm như sau :

Khử các phần tử trên cột chứa pivot

Tính dấu hiệu tối ưu

III- PHƯƠNG PHÁP BIẾN GIẢ CẢI BIÊN

1- Bài toán cải biên

a- Cải biên bài toán quy hoạch tuyến tính

Người ta có thể biến đổi một bài toán quy hoạch tuyến tính chính tắc thành dạng chuẩn bằng cách cộng một cách phù hợp vào vế trái của ràng buộc i một biến giả xn+i≥ 0 để làm xuất hiện ma trận đơn vị Vì các biến giả cải biên có ảnh hưởng đến hàm mục tiêu nên cũng sẽ có

sự cải biên hàm mục tiêu

Vậy, người ta có thể biến đổi bài toán quy hoạch tuyến tính tổng quát, gọi là bài toán xuất phát, thành bài toán dạng chuẩn, gọi là bài toán cải biên (mở rộng)

Ví dụ :

Biến đổi bài toán quy hoạch tuyến tính sau đây thành dạng chuẩn

)4,3,2,1j( 0 x

28x

8x3

18x

6xx4

25x

5x5x

xxxx2)x(zmax

j

4 2

4 3 2

4 2 1

4 3 2 1

=+

−

−

=++

−++

=

Bài toán xuất phát có các biến, ma trận ràng buộc và chi phí :

]1- 1 1 2[c

8 0 3 0

6 1- 4- 0

5 0 5 1 A

] x x xx[x

T

4 3 2 1 T

Bằng cách thêm biến giả x5, x6 lần lượt vào ràng buộc 2 và 3 Ta được bài toán cải biên :

)6,5,4,3,2,1j( 0x

28x

x8x3

18xx6xx4

25x

5x5x

)xx(Mxxxx2)x(z max

j

6 4 2

5 4 3 2

4 2 1

6 5 4

3 2 1

=++

−

−

=++

+

−

−++

z′ là hàm mục tiêu cải biên sẽ được giải thích trong phần tiếp theo

Các biến, ma trận ràng buộc các hệ số và chi phí của bài toán cải biên là

]M- M- 1- 1 1 2[c

1 0 8 0 3 0

0 1 6 1- 4- 0

0 0 5 0 5 1 A

] x x x x xx[x

T

6 5 4 3 2 1 T

0 0 x

- Nếu bài toán cải biên (dạng chuẩn) có phương án tối ưu thì cũng sẽ phương

án cơ sở tối ưu

Ví dụ

1- Xét bài toán :

0x

3

2x3

1x3

4x3

2x3

1x

52x5x7xx

09x3x

5xx2xx)x(z min

j

5 4

3 2

1

5 4 3 2

4 3

5 4 2 1

++

0x

75x5x

3

1xx3

1x32

9x7x16xz(x)

min

j

2 1

3 2 1

3 2 1

−

=+

−

−

++

70x

xx

70x

x

x1 2 33- Xét bài toán :

1,2,3)(j

x

18xxx

502x

x2x

27x

2xx

2x4x2xz(x) min

j

3 2 1

3 2 1

3 2 1

3 2 1

=+

−

−+

2- Phương pháp hai pha

Đây là phương pháp thuận lợi cho việc lập trình ứng dụng giải thuật đơn hình cải tiến

Ví dụ : Xét bài toán quy hoạch tuyến tính

1,2,3)(j

0x

3

7x3x2x

3

8x2x2x

xx4x3)x(z max

j

3 2 1

3 2 1

3 2 1

≤++

++

=

Đưa bài toán về dạng chính tắc bằng cách thêm biến phụ x4 , x5 ta được

1,2,3,4,5)(j

0x

3

7xx3x2x

3

8xx2x2x

xx4x3)x(z max

j

5 3 2 1

4 3 2 1

3 2 1

=+++

++

0 1 2 2 1

Áp dụng phương pháp đơn hình cải biên hai pha như sau :

Pha 1

Thêm biến giả (cải biên ) x6 ≥ 0 vào ràng buộc thứ hai để được ma trận đơn vị Khi đó bài toán cải biên có dạng :

6)1,2,3,4,5,(j

0x

3

7xxx3x2x

3

8xx2x2x

x)x(w min

j

6 5 3 2 1

4 3 2 1

−++

=+++

0 0 1 2 2 1AGiải bài toán cải biên bằng giải thuật đơn hình cải tiến

8b 110321

001221

1 6 1 2 3 0 -1 1

37

T 0

37

T 1

Ta được phương án tối ưu Xong pha 1 Chuyển sang pha 2

Loại bỏ biến giả cải biên x6 ≥ 0

10b

3

1013

23

2103

23

1A

0 0

cT =Bước lặp k=0

T 0

T 1

314

T 2

316

Bước lặp k=3

3

B

c i B3 x1 x2 x3 x4 x5 b3

0 5 0 0 -1 1 1

31

38

T 3

Kết quả của bài toán đã cho :

Phương án tối ưu

0x

0x

0x3

8x

5 4 3 2 1

Giá trị hàm mục tiêu z(x)=z(x3)= 8

3- Phương pháp M vô cùng lớn

Phương pháp M vô cùng lớn ( M là số vô cùng lớn ) tương tự như phương pháp hai pha, ngoại trừ ở pha 1 hàm mục tiêu cải biên có dạng sau đây cho bài toán max/min

max [z(x) - M*( tổng các biến giả cải biên) ] min [z(x) + M*( tổng các biến giả cải biên) ] Bằng phương pháp này, trong quá trình tối ưu, các biến giả cải biên sẽ được loại dần ra khỏi ma trận cơ sở : tất cả đều bằng 0 Nếu trong quá trình tìm phương án tối ưu mà không loại bỏ được các biến giả cải biên ra khỏi cơ sở thì bài toán vô nghiệm

So với phương pháp hai pha thì phương pháp này tránh được việc phải cập nhật lại dữ liệu cho bài toán gốc nhưng không tiện lợi bằng trong lập trình ứng dụng

Ví dụ : Xét bài toán tương tự như trên

0x

3

7xx3x2x

3

8xx2x2x

xx4x3)x(z max

j

5 3 2 1

4 3 2 1

3 2 1

=+++

++

=

Thêm biến giả cải biên x6 ≥ 0 vào ràng buộc thứ hai đồng thời cải biên hàm mục tiêu theo như trên ta được :

6)1,2,3,4,5,(j

0x

3

7xxx3x2x

3

8xx2x2x

Mxx

x4x3)x(wmax

j

6 5 3 2 1

4 3 2 1

6 3

2 1

−++

=+++

−++

8b

110321

001221

37

T 0

3

M7

1 3

3

13

3

1

M3

5 −

−

97

Do x6 = 0 (vì ngoài cơ sở) nên bị loại ra khỏi bảng và ta tiếp tục tìm phương

án tối ưu cho bài toán gốc đã cho có phương án cơ sở khả thi được khởi tạo như sau :

T 0

3

1

97

Các bước tiếp theo được thực hiện giống như phương pháp hai pha

IV- QUY HOẠCH TUYẾN TÍNH SUY BIẾN

Khi thực hiện thuật toán đơn hình trường hợp bất thường có thể xảy ra là khi xác định biến ra thì tồn tại tỷ số 0

a

bik

i = , tức là tồn tại bi=0, hay không có tỷ số nào dương thật sự Người ta xem đây là trường hợp suy biến Khi một bảng đơn hình rơi vào tình trạng suy biến thì có thể gây khó khăn mà cũng có thể không khi ta tiếp tục thực hiện thuật toán đơn hình

1- Các ví dụ về quy hoạch tuyến tính suy biến

Ví dụ 1 : xét quy hoạch tuyến tính :

0x,x

0x2

6x3

2x2x

xx7z(x)min

2 1 1 1

2 1

2 1

=

Đưa bài toán về dạng chuẩn :

0x,x

0xx2

6xx3

2xx2x

xx7z(x)min

2 1

5 1

4 1

3 2 1

2 1

−

=+

−

=+

−

−+

=

với ma trận hệ số là :

x1 x2 x3 x4 x5 b

1 -2 1 0 0 2 -3 0 0 1 0 6 -2 0 0 0 1 0

có chứa ma trận đơn vị Áp dụng thuật toán đơn hình cải tiến ta được :

cB iB x1 x2 x3 x4 x5 b

0 3 1 -2 1 0 0 2

0 4 -3 0 0 1 0 6

0 5 -2 0 0 0 1 0 T

=

≥+

=

≥+

=

0x00x

0x06x

0x22x

2 3

2 4

2 3

0x2

2x

2 2 2

Như vậy x2 có thể lớn tùy ý nên hàm mục tiêu không bị giới nội Vậy bài toán không có phương án tối ưu Trường hợp này ở bảng đơn hình không có tỷ số nào dương thật sự để xác định biến ra

Ví dụ 2 : xét quy hoạch tuyến tính :

0x,x

0x2

6x3

2x2x

xx7z(x)min

2 1 1 1

2 1

2 1

−+

=

Đưa bài toán về dạng chuẩn :

0x,x

0xx2

6xx3

2xx2x

xx7z(x)min

2 1

5 1

4 1

3 2 1

2 1

−

=+

−

=++

−+

=

với ma trận hệ số là :

x1 x2 x3 x4 x5 b

1 2 1 0 0 2 -3 0 0 1 0 6 -2 0 0 0 1 0

có chứa ma trận đơn vị Áp dụng thuật toán đơn hình cải tiến ta được :

cB iB x1 x2 x3 x4 x5 b

0 4 -3 0 0 1 0 6

0 5 -2 0 0 0 1 0 T

T

cB iB x1 x2 x3 x4 x5 b -1 2

Đây là bảng đơn hình tối ưu

Ví dụ 3 : xét quy hoạch tuyến tính :

0x,x,x

0xxx

1x2

1x21

x2

3x2x2

1-3w(x) min

3 2 1

3 2 1

3 1

3 2

−

≤+

+

−+

=

Đưa bài toán về dạng chuẩn :

0x,x,x,x,x

0xxxx

1xx2

1x21

x2

3x2x2

1-3w(x) min

5 4 3 2 1

5 3 2 1

4 3 1

3 2

−+

−

=++

+

−+

=

với ma trận hệ số :

có chứa ma trận đơn vị Áp dụng giải thuật đơn hình cải tiến :

Đây là bảng đơn hình tối ưu

Ví dụ 4 : xét quy hoạch tuyến tính

1x

0xx5,0x5,1x5,0

0x9x5,2x5,5x5,0

x24x9x57x10z(x)

min

4 3 2 1 1

4 3 2

1

4 3 2

1

4 3

2 1

−

−

≤+

−

−

+++

−

=

Đưa bài toán về dạng chuẩn

0x,x,x,x,x,x,x

1xx

0xxx5,0x5,1x5,0

0xx9x5,2x5,5x5,0

x24x9x57x10w(x)

min

7 6 5 4 3 2 1

7 1

6 4 3 2

1

5 4 3 2

1

4 3

2 1

=++

−

−

=++

−

−

+++

9 3 2 0 1 -8 -1,5 5,5 0 0

0 7 1 0 0 0 0 0 1 1 T

2- Xử lý trường hợp suy biến

Theo các ví dụ trên, trong trường hợp quy hoạch tuyến tính suy biến thì sau một số lần lặp có thể phương án nhận được vẫn như cũ mà không có sự thay đổi nào,

có thể phương án nhận được tốt hơn, có thể phương án nhận được là một phương án

đã nhận trước đó rồi và từ đó cứ xoay vòng mãi Do đó nếu không có biện pháp phòng ngừa thì thuật toán đơn hình sẽ có thể kéo dài vô tận

Khi thực hiện thuật toán đơn hình thì hiện tượng suy biến xảy ra khi có sự tình

cờ khử lẫn nhau làm cho tồn tại b nào đó bằng 0 Trong trường hợp này có thể có i

nhiều biến thỏa điều kiện của biến ra Gặp trường hợp này cần phải lựa chọn biến ra sao cho tránh được hiện tượng xoay vòng