Bài giảng Toán kinh tế - Chương 1: Bài toán quy hoạch tuyến tính phương pháp đơn hình

47 trang huongle 6020

Download

Bạn đang xem 20 trang mẫu của tài liệu "Bài giảng Toán kinh tế - Chương 1: Bài toán quy hoạch tuyến tính phương pháp đơn hình", để tải tài liệu gốc về máy bạn click vào nút DOWNLOAD ở trên

Tài liệu đính kèm:

bai_giang_toan_kinh_te_chuong_1_bai_toan_quy_hoach_tuyen_tin.pdf

Nội dung text: Bài giảng Toán kinh tế - Chương 1: Bài toán quy hoạch tuyến tính phương pháp đơn hình

Toán Kinh tế - Trường Cao Đẳng Công Nghiệp Nam Định Chương 1. BÀI TOÁN QUY HOẠCH TUYẾN TÍNH. PHƯƠNG PHÁP ĐƠN HÌNH 1.1. Bài toán quy hoạch tuyến tính 1.1.1. Một số mô hình thực tế A. Bài toán lập kế hoạch sản xuất Một cơ sở có thể sản xuất hai loại sản phẩm A và B, từ các nguyên liệu I, II, III. Chi phí từng loại nguyên liệu và tiền lãi của một đơn vị sản phẩm, cũng như dự trữ nguyên liệu cho trong bảng sau đây: Nguyên liệu I II III Lãi Sản phẩm A 2 0 1 3 B 1 1 0 5 Dự trữ 8 4 3 Hãy lập bài toán thể hiện kế hoạch sản xuất sao cho có tổng số lãi lớn nhất, trên cơ sở dự trữ nguyên liệu đã có. Lập bài toán: Gọi x, y lần lượt là số sản phẩm A và B được sản xuất ( x,y ≥ 0, đơn vị sản phẩm). Khi đó ta cần tìm x,y ≥ 0 sao cho đạt lãi lớn nhất. f ()3Xxy=+→ 5 max với điều kiện nguyên liệu: 28x +≤y ; 1.y ≤ 4; 1.x ≤ 3; Tức là cần giải bài toán: f (Xxy )=+→ 3 5 max ⎧28xy+=; ⎪ ⎪y ≤ 4; với điều kiện:⎨ ⎪x ≤ 3; ⎩⎪xy,0;≥ - 1 - Nguyễn Hải Đăng - Khoa KHCB&KTCS
Toán Kinh tế - Trường Cao Đẳng Công Nghiệp Nam Định B. Bài toán phân công lao động: Một lớp học cần tổ chức lao động với hai loại công việc: xúc đất và chuyển đất. Lao động của lớp được chia làm 3 loại A, B, C, với số lượng lần lượt là 10, 20, 12. Năng suất của từng loại lao động trên từng công việc cho trong bảng dưới đây: Lao động A(10) B(20) C(12) Công việc Xúc đất 6 5 4 Chuyển đất 4 3 2 Hãy tổ chức lao động sao cho có tổng năng suất lớn nhất. Lập bài toán: Gọi xij là số lao động loại j làm công việc i(j=1,2;xij ≥ 0 , nguyên). Khi đó, năng suất lao động của công việc đào đất sẽ là: 654x11++xx 12 13 ; còn chuyển đất sẽ là : 432x21+ xx 22+ 23; Ta thấy rằng để có năng suất lớn nhất thì không thể có lao động dư thừa, tức là phải cân bằng giữa hai công việc. Vì vậy ta có bài toán sau: 654xxx11++ 12 13 →max; ⎧654432xxxxx11++−++= 12 13 21 22 x 23 0; ⎪ ⎪xx11+= 21 10; với điều kiện ⎨ xx+=20; ⎪ 12 22 ⎪ ⎩xx13+= 23 12; C. Bài toán khẩu phần thức ăn: Một khẩu phần thức ăn có khối lượng P, có thể cấu tạo từ n loại thức ăn. Gía mua một đơn vị thức ăn loại j là cj. Để đảm bảo cơ thể phát triển bình thường thì khẩu phần cần m loại chất dinh dưỡng. Chất dinh dưỡng thứ i cần tối thiểu cho khẩu phần là bi và có trong một đơn vị thức ăn loại j là aij. Hỏi nên cấu tạo một khẩu phần thức ăn như thế nào để ăn đủ no, đủ chất dinh dưỡng mà có giá thành rẻ nhất. Lập bài toán: Gọi xj (xj ≥ 0 ) là số đơn vị thức ăn loại j được cấu tạo trong khẩu phần. Khi đó, giá thành của khẩu phần là: - 2 - Nguyễn Hải Đăng - Khoa KHCB&KTCS
Toán Kinh tế - Trường Cao Đẳng Công Nghiệp Nam Định n f ()Xc= ∑ j xj ; j=1 Vì phải đảm bảo thoả mãn điều kiện đủ no và đủ chất, tức là: nn ∑∑x ji=≥=Paxbim,,jjj1,. jj==11 n Ta có bài toán sau: fX()=→∑ cxjj min j=1 ⎧ n ∑ xPj = ; ⎪ j =1 ⎪ ⎪ n với điều kiện ⎨ ∑ axij j≥= b i ,1, i m; ⎪ j =1 ⎪ xjnj ≥=0, 1, ; ⎩⎪ Ta thấy rằng ba bài toán trên đều thuộc bài toán tổng quát. 1.1.2. Bài toán quy hoạch tuyến tính tổng quát Bài toán tổng quát của QHTT có dạng : n fX()=→∑ cxjj min(ax) m j=1 ⎧ n ∑ axij j== b i ,1, i k ⎪ j =1 ⎪ ⎪ n với điều kiện ⎨∑ axij j≥=+ b i ,1 i k, m ⎪ j =1 ⎪ x j ≥=0,jrrn 1, , ≤ ⎩⎪ Để phân biệt tính chất của các ràng buộc đối với một phương án, ta làm quen với hai khái niệm : ràng buộc chặt và ràng buộc lỏng. Định nghĩa 1: Nếu đối với phương án x mà ràng buộc i thỏa mãn với dấu đẳng n thức, nghĩa là ∑axij j= b i thì ta nói phương án x thỏa mãn chặt ràng buộc i j=1 Nếu đối với phương án x mà ràng buộc i thỏa mãn với dấu bất đẳng thức thực sự, n nghĩa là ∑axij j> b i thì ta nói phương án x thỏa mãn lỏng ràng buộc i j=1 Định nghĩa 2 : Ta gọi một phương án thỏa mãn chặt n ràng buộc độc lập tuyến tính là phương án cực biên. Một phương án cực biên thỏa mãn chặt đúng n ràng buộc - 3 - Nguyễn Hải Đăng - Khoa KHCB&KTCS
Toán Kinh tế - Trường Cao Đẳng Công Nghiệp Nam Định gọi là phương án cực biên không suy biến, thỏa mãn chặt hơn n ràng buộc gọi là phương án cực biên suy biến. Định nghĩa 3: Một phương án mà tại đó hàm mục tiêu đạt cực tiểu ( cực đại ) gọi là phương án tối ưu. Bài toán có ít nhất một phương án tối ưu gọi là bài toán giải được, bài toán không có phương án hoặc có phương án nhưng hàm mục tiêu không bị chặn dưới ( trên ) trên tập phương án gọi là không giải được. Để nhất quán trong lập luận, ta xét bài toán tìm cực tiểu, sau đó ta xét cách đưa bài toán tìm cực đại về bài toán tìm cực tiểu. * Chuyển bài toán tìm cực đại về bài toán tìm cực tiểu : Nếu gặp bài toán tìm max, tức là : n fX()=→∑ cxjj max j=1 XM∈ thì giữ nguyên ràng buộc, ta đưa nó về dạng bài toán tìm min : n gX()=− f () X =−∑ cxjj → min j=1 XM∈ Chứng minh : Nếu bài toán tìm min có phương án tối ưu là X* thì bài toán tìm max cũng có phương án tối ưu là X* và g(X)= - f(X). Thật vậy, X* là phương án tối ưu của bài toán tìm min, tức là nn f ()XcxcxX=≤∀∈∑∑jj jj, M jj==11 nn * ⇒−∑∑cxjj ≥ cx jj, ∀ X ∈ M jj==11 hay − f ()XgXgXX =≥∀∈ () (), M Vậy X* là phương án tối ưu của bài toán max và n * fmax =−∑cxjj =− gmin j=1 1.1.3. Dạng chính tắc của bài toán quy hoạch tuyến tính - 4 - Nguyễn Hải Đăng - Khoa KHCB&KTCS
Toán Kinh tế - Trường Cao Đẳng Công Nghiệp Nam Định Người ta thường xét bài toán QHTT dưới dạng sau: n fX()=→∑ cxjj min (1.1) j=1 n ⎧ ax== b,1, i m (1.2) ⎪ ∑ ij j i với điều kiện ⎨ j =1 ⎪ xjn≥=0, 1, (1.3) ⎩ j Bài toán (1.1), (1.2), (1.3) được gọi là Bài toán Quy hoạch tuyến tính dạng chính tắc. Kí hiệu ma trận hàng cccc= (,12 , ,)nn 1× và các ma trận : ⎛⎞x1 ⎛⎞b1 ⎛⎞a1j ⎜⎟ ⎜⎟x ⎜⎟b a ⎜⎟2 ⎜⎟2 ⎜⎟2j x = , b = , Aa= () , A ==jn1, ⎜⎟ ⎜⎟ ij mn× j ⎜⎟ ⎜⎟ ⎜⎟ ⎜⎟ ⎜⎟a ⎝⎠xn ⎝⎠bm ⎝⎠mj Ta có bài toán ở dạng ma trận như sau: f(x) = cx → min ⎧Axb= Với điều kiện ⎨ ⎩x ≥ 0 và bài toán ở dạng véc tơ như sau: f(x) = cx → min ⎧A11xAx+ 2 2++ Axnn =b Với điều kiện ⎨ ⎩xx12, , , xn ≥ 0 Đối với bài toán dạng chính tắc ta có một số tính chất và khái niệm quan trọng của phương án cực biên như sau : Tính chất 1( Nhận dạng phương án cực biên ) : Phương án x của bài toán dạng chính tắc là cực biên khi và chỉ khi hệ thống các véc tơ {Axjj:0> } độc lập tuyến tính. Với giả thiết hạng[A] = m thì một phương án cực biên không suy biến có đúng m thành phần dương, suy biến có ít hơn m thành phần dương. Chứng minh Lấy một phương án x bất kì, giả sử p thành phần đầu của x là dương, tức là x j >=0(jp 1, ), suy ra xk ==+0(kp 1, n ). Thành lập ma trận tương ứng với các ràng - 5 - Nguyễn Hải Đăng - Khoa KHCB&KTCS
Toán Kinh tế - Trường Cao Đẳng Công Nghiệp Nam Định buộc chặt của phương án x ( bao gồm m ràng buộc đẳng thức và n - p ràng buộc chặt về dấu ), kí hiệu là C, ta có : ⎡⎤aa11 12 aa1pp 1++ 1 a 1 p 2 a1 n ⎢⎥ aa21 22 aa2pp 2++ 1 a 2 p 2 a2 n ⎢⎥m ⎢⎥ P ⎢⎥ ⎢aam12 m aamp mp++1 a mp 2 amn ⎥ C = ⎢ 00 01 0 0⎥ ⎢ ⎥ ⎢ 00 00 1 0⎥ ⎢ 00 00 0 0⎥ n-p ⎢ ⎥ ⎢ ⎥ ⎢ ⎥ ⎣⎦00 00 0 1 Theo cách tính hạng của ma trận thì hạng[C] = n – p + hạng[P] Từ định nghĩa suy ra, x là phương án cực biên khi và chỉ khi hạng[C] = n, nghĩa là khi và chỉ khi hạng[P] = p, nói cách khác thì {Aj :1,jp= } hay {Axjj:> 0}độc lập tuyến tính. Từ đây đối với bài toán dạng chính tắc, không mất tính tổng quát , ta giả thiết : • Hệ (1.2) có đúng m phương trình độc lập tuyến tính. • ∀bi ≥0, i = 1, m • m < n (trong trường hợp m≥ n thì tập phương án có nhiều nhất một điểm, do vậy việc tìm phương án tối ưu là tầm thường). Vì hạng của ma trận A bằng m nên số véc tơ điều kiện Aj độc lập tuyến tính cực đại là m, do đó bất kì phương án cực biên nào cũng tương ứng với ít nhất một hệ véc tơ độc lập tuyến tính cực đại, từ đó ta có định nghĩa sau : Định nghĩa 4: Ta gọi một hệ m véc tơ {Aj} độc lập tuyến tính bao hàm hệ thống các véc tơ tương ứng với các thành phần dương của phương án cực biên x là cơ sở của phương án cực biên ấy, kí hiệu một cách quy ước là J, trong đó J = { j : Aj nằm trong cơ sở } Một phương án cực biên không suy biến có đúng m thành phần dương, m véc tơ Aj tương ứng độc lập tuyến tính nên có một cơ sở duy nhất, đó chính là các véc tơ tương ứng với các thành phần dương ; còn đối với phương án cực biên suy biến thì có - 6 - Nguyễn Hải Đăng - Khoa KHCB&KTCS
Toán Kinh tế - Trường Cao Đẳng Công Nghiệp Nam Định nhiều cơ sở khác nhau, phần chung của chúng là các véc tơ tương ứng với các thành phần dương. Như vậy khi nói phương án cực biên có cơ sở J, cần hiểu J có 3 nội dung sau: • Số phần tử của J: Jm= • {Aj : jJ∈ }độc lập tuyến tính • {}{AjJjj::∈⊃ Axj >0} Tính chất 2 ( Tính hữu hạn của số phương án cực biên ): Số phương án cực biên của mọi bài toán quy hoạch tuyến tính đều hữu hạn Thật vậy: Vì mỗi phương án cực biên đều tương ứng với một hệ n ràng buộc chặt độc lập tuyến tính, số hệ gồm n phương trình độc lập tuyến tính là hữu hạn, do đó số phương án cực biên cũng là hữu hạn. Tính chất 3 ( Sự tồn tại phương án tối ưu ): Nếu bài toán dạng chính tắc có phương án và hàm mục tiêu bị chặn dưới trên tập phương án thì bài toán có phương án cực biên tối ưu. Thật vậy, giả sử bài toán có các phương án cực biên là x12,xx , , k . Đặt f ()xMinfxikli=={ ():1,}. Do tập nghiệm ( hay tập phương án của bài toán quy hoạch tuyến tính) của hệ ràng buộc là một đa diện lồi, cho nên mọi nghiệm (phương án) x đều có thể phân tích thành: kk i xx==∑∑λλii,1 ii==11 kk k ii ll Suy ra f ()xf==≥ (∑∑λλii x ) fx ()( ∑ λ i )() fxfxxM =∀ () ∈ ii==11 i = 1 Vậy xl chính là phương án cực biên tối ưu. Một lớp quan trọng của các bài toán quy hoạch tuyến tính dạng chính tắc là bài toán dưới đây gọi là bài toán dạng chuẩn : - 7 - Nguyễn Hải Đăng - Khoa KHCB&KTCS
Toán Kinh tế - Trường Cao Đẳng Công Nghiệp Nam Định n fx()=→∑ cxjj min j=1 ⎧ x1 ++++ax1m+1 m+1 ax1m+2 m+2 ax1n n=b 1 ⎪ x ++++ax ax ax=b ⎪ 2 2m+1 m+1 2m+2 m+2 2n n 2 ⎪ ⎨ ⎪ x ++++ax ax ax=b ⎪ m mm+1 m+1 mm+2 m+2 mn n m xjn≥=0, 1, ⎩⎪ j trong đó bimi ≥=0( 1, ) , nghĩa là bài toán dạng chính tắc có vế phải không âm và mỗi phương trình đều có một biến số với hệ số bằng 1 đồng thời không có trong các phương trình khác ( gọi là biến cô lập với hệ số bằng 1). Từ hệ phương trình ràng buộc của bài toán dễ dàng suy ra một phương án: 0 xbbb= (12 , , ,m ,0,0, ,0) Đây chính là một phương án cực biên có hệ cơ sở tương ứng là ⎛⎞10 ⎛⎞ ⎛⎞0 ⎜⎟01 ⎜⎟ ⎜⎟0 AA==⎜⎟, ⎜⎟ , , A = ⎜⎟ 12⎜⎟ ⎜⎟ m ⎜⎟ ⎜⎟ ⎜⎟ ⎜⎟ ⎝⎠00 ⎝⎠ ⎝⎠1 1.1.4. Đưa bài toán quy hoạch tuyến tính về dạng chính tắc Phương pháp: Ta có thể đưa bài toán tuyến tính tổng quát về bài toán tuyến tính dạng chính tắc tương đương nhờ các quy tắc sau: • Nếu có max f(X) thì đổ thành {min− f (X )} . n n • Nếu có bất đẳng thức ∑axij j≥ b i hoặc ∑axij j≤ b i thì ta đưa thêm ẩn j=1 j=1 phụ xn+i ≥ 0 ,với hệ số hàm mục tiêu cn+i = 0 để có: n n ∑axij j−= x n+ i b i hoặc ∑axij j+ x n+ i= b i ; j=1 j=1 • Nếu có ẩn xk chưa ràng buộc về dấu, thì ta có thể thay nó bởi hai biến mới x và x không âm, theo công thức: k' k" x = x - x . k k' k" Ví dụ 1.1 Đưa bài toán sau về dạng chính tắc: min {x123−−xx}; - 8 - Nguyễn Hải Đăng - Khoa KHCB&KTCS
Toán Kinh tế - Trường Cao Đẳng Công Nghiệp Nam Định ⎧654432xxxxx11++−++= 12 13 21 22 x 23 0; ⎪ xxx123++=5; ⎪ với điều kiện ⎨xxx123−+≤23; ⎪xxx+−≥4; ⎪ 123 ⎩⎪xx13,0;≥ Giải: Ta thấy có bất đẳng thức xxx123− 2+≤3 nên ta đưa thêm ẩn phụ xx45,0≥ '" Mặt khác, có ẩn x2 chưa ràng buộc về dấu, do đó ta thay x2 bởi x2− x2. Khi đó, bài toán ban đầu được chuyển về dạng sau: '" ()xxxx1223−+− →min '" ⎧ xxxx1223+−+=5 ⎪ '" ⎪ xxxxx1− 22 2+++= 2343 với điều kiện ⎨ '" ⎪ xxxxx12235+−−−=4 ⎪ '" ⎩ xxxxxx122345,,,,,≥ 0 Ví dụ 1.2 Đưa bài toán QHTT sau về dạng chính tắc: 222xx12−+ xx 34 +− x 5 →min ⎧xxxxx12345−+++≤227(1) ⎪xxx++≥−21(2) ⎪ 234 ⎪2310(3)xx34++ x 5 ≥ với điều kiện ⎨ xx+−220 xx += ⎪ 12 34 ⎪xx,0≥ ⎪ 15 ⎩⎪x4 ≤ 0 Giải: '"'" Vì x2, x3 chưa ràng buộc về dấu nên ta thay x2 bởi xxxx2222−(,≥ 0), x3 bởi '"'" '' xxxx3333−≥(, 0), x4 ≤ 0 nên thay x4 bởi −xx44(0≥ ). Vì có các bất đẳng thức (1), (2), (3) nên ta thêm các ẩn phụ x6, x7, x8. Từ đó, ta được bài toán sau: '" '" ' 2(xxx122−−+ )2() xxx 334 −−−→ 2 x 5 min - 9 - Nguyễn Hải Đăng - Khoa KHCB&KTCS
Toán Kinh tế - Trường Cao Đẳng Công Nghiệp Nam Định '" '" ' ⎧xxxxxxxx1−−+−−++=2( 22 ) 33 2 4567 ⎪ ()2()xx'"−+ xx '" −−−=− xx ' 1 ⎪ 22 33 47 ⎪ '" ' Với điều kiện ⎨2(xx22−−+−= ) x 4 3 xx 58 10 ⎪xxx+−−()2()20'" xxx '" −−= ' ⎪ 122 334 ⎪ '"'"' ⎩xxxxxxxxxx1567822334,,,,,,,,,≥ 0 1.2. Phương pháp đơn hình 1.2.1. Tư tưởng của phương pháp đơn hình Xét bài toán quy hoạch tuyến tính dạng chính tắc: n fx()=→∑ cxjj min j=1 ⎧ n ⎪∑ axij ji== b,1, i m ⎨ j=1 ⎪ ⎩x j ≥ 0 Dạng véctơ của bài toán: f ()xcx=→ min ⎧ n Ax b Với điều kiện: ⎪∑ jj= ⎨ j=1 ⎪ ⎩x j ≥=0,jn 1, Ta đã biết rằng : - Nếu bài toán có phương án thì có phương án cực biên - Nếu bài toán có phương án tối ưu thì cũng có phương án cực biên tối ưu. Số phương án cực biên là hữu hạn. Do đó, ta có thể tìm một phương án tối ưu(hay một lời giải của bài toán) trong tập hợp các phương án cực biên. Tập hợp này là hữu hạn. Vì vậy Dantzig đề xuất thuật toán đơn hình như sau: Xuất phát từ một phương án cực biên x0. Kiểm tra xem x0 có là phương án tối ưu hay chưa. Nếu x0 chưa phải là phương án tối ưu thì tìm cách cải tiến nó để được phương pháp khác là x1 tốt hơn x0, tức là f(x1) < f( x0). Qúa trình này lặp lại nhiều lần. Vì số phương án cực biên là hữu hạn nên sau một số hữu hạn lần lặp ta sẽ tìm thấy phương án cực biên tối ưu. Để thực hiện thuật toán đề ra ở trên, ta cần làm rõ hai vấn đề sau: - 10 - Nguyễn Hải Đăng - Khoa KHCB&KTCS
Toán Kinh tế - Trường Cao Đẳng Công Nghiệp Nam Định 1. Làm thế nào để biết một phương án cực biên đã cho là tối ưu hay chưa, tức là cần tìm « dấu hiệu tối ưu ». 2. Làm thế nào để từ một phương án cực biên chưa tối ưu tìm được một phương án cực biên tốt hơn nó. 1.2.2. Biểu diễn qua cơ sở. Dấu hiệu tối ưu 0 Gỉa sử có phương án cực biên x với cơ sở J0(tức là hệ véctơ cột độc lập tuyến tính {Ajj , ∈ J0} và xj > 0). Ta có: n 0 00 AX= b hay bxAx==∑∑jj jjA (1) jj=∈1 J0 0 (vì x j =∀∉0 jJ0 ). Với mỗi kJ∉ 0 , ta biểu diễn véc tơ Ak qua các véc tơ cơ sở A =∀zA() k∉J {Ajj , ∈ J0} : kj∑ kj 0 jJ∉ 0 Gỉa sử x ∈ M là một phương án bất kỳ. Ta có : n bxAxAxAxAxzAx==+=+∑∑∑∑∑∑∑∑jj jj kk jj k jkj =+() j zx jkkA j (2) jj=∈∉∈1 JkJjJk00000 ∉∈∈∉JjJjJkJ Vì {Ajj , ∈ J0} độc lập tuyến tính nên từ (1) và (2) suy ra: 0 x jj=+xzxjJ∑ jkk() ∀∈0 kJ∉ 0 Hay: x jj=−xzxj∑ jkk ∀∈J0 (3) kJ∉ Khi đó, ta có: n f ()xcxcxc==+∑∑∑jj jj kkx (4) jj=∈1 Jk00 ∉J Thay (3) vào ta được : nn 00 f ()xxzxccxcxzcc=−∑∑ (jjkkjkkjjj ) + ∑ = ∑ − ∑∑()kjk −xk jkJkJj=∉1100 ∉ = kJjJ ∉∈00 Ký hiệu: Δ=kj∑ zckjk − c( k ∉ J0 ) gọi là ước lượng của véc tơ cột Ak theo cơ sở J và: jJ∈ 0 0 f ()xfx=−Δ ( ) ∑ kk x kJ∉ 0 - 11 - Nguyễn Hải Đăng - Khoa KHCB&KTCS
Toán Kinh tế - Trường Cao Đẳng Công Nghiệp Nam Định 0 Nhận xét: Do x ≥ 0 nên nếu ∀∉kJ0 ,0 Δ≤k thì f ()xfxxM≥∀∈ ( ), do đó ta có dấu hiệu tối ưu sau đây: 0 Phương án cực biên x với cơ sở J0 là phương án tối ưu khi và chỉ khi Δk0≤ 0,∀∉ k J . 1.2.3. Công thức biến đổi. Bảng đơn hình Giả sử x0 với cơ sở J là một phương án cực biên nhưng chưa phải là phương án tối ưu, khi đó ∃∉kJ0 sao cho Δ>k 0. Giả sử s là một chỉ số trong các chỉ số nói trên: sJ∉Δ>0 ,0s . Theo thuật toán trên ta cần cải tiến x0 để nhận được một phương án cực biên mới tốt hơn. Ta sẽ tìm một phương án cực biên mới x1 ứng với cơ sở J1 chỉ khác J một véc 1 tơ : JJr=(\)0 ∪s, tức là ta đưa véc tơ As vào cơ sở thay cho véc tơ Ar bị loại ra khỏi cơ sở J. Vì mọi thành phần ngoài cơ sở của một phương án cực biên đều bằng 0 nên phương án cực biên mới x1 có cơ sở JJr1 = (\)∪ s là: 1 1 ⎧0)víi k∉≠ J0 , k s(tøc lμ kJ ∉ xk = ⎨ (5) ⎩θ víi k = s Trong đó θ là một số dương sẽ được xác định sau sao cho x1 là một phương án cực biên. * Tìm điều kiện của θ để x1 là một phương án: Thay (5) vào (3) ta được: 10 10 x jj=−xzxxzjJ∑ jkkj.( =−θ js ∀∈0 ) kJ∉ 0 ⎧ 1 0)víi j∉≠J,0 j s(tøc lμ j ∉J 1 ⎪ Hay x j = ⎨θ víi j = s (6) ⎪ 0 ⎩xvíij0−∈θ z js jJ Để x1 là một phương án thì nó phải thỏa mãn các điều kiện buộc Ax=b và x ≥ 0. + Ta thấy rằng với mọiθ , ta có: 11 1 0 0 Ax =+=−+=−∑∑∑x jjAxAxzAAxAzAAbAA jj() jθθ jsj s ∑ jj θ ∑ jsj +=−+ θθθ s s s=b jJ∈∉ jJ jJ ∈ jJ∈00jJ ∈ Vậy với mọi θ thì ràng buộc thứ nhất thỏa mãn. Ta chỉ cần tìm θ sao cho x1 ≥ 0. Có hai trường hợp xảy ra: - 12 - Nguyễn Hải Đăng - Khoa KHCB&KTCS
Toán Kinh tế - Trường Cao Đẳng Công Nghiệp Nam Định 1 1 + Trường hợp 1: Nếu zjs ≤ 0 với mọi j ∈J thì x ≥ 0 với mọi θ > 0 và x là phương án của bài toán. Nhưng do Δ>s 0 : 10 10 fx()=−Δ=−Δ→−∞→ fx ()∑ kk x fx() sθθkhi +∞ kJ∉ 0 Như vậy hàm mục tiêu không bị chặn dưới trên miền ràng buộc. Khi đó bài toán không có lời giải hữu hạn. + Trường hợp 2: Nếu ∃∈j J0 để zjs > 0, khi đó theo (6) số θ không thể lớn tùy ý, nó 0 phải thỏa mãn xzjj−θ s≥∀0, j∈ Jm0jsà z > 0 . Giá trị θ lớn nhất chỉ có thể bằng 0 ⎪⎪⎧⎫x j min⎨⎬ : jJ∈>0jsmà z 0 ⎩⎭⎪⎪z js 0 1 Nếu vượt qua miền đó sẽ có một trong các ( xzjjs−θ 0. zrs Rõ ràng, phương án x1 được xác định theo công thức (7) là phương án cực biên với 1 cơ sở JJr=∪(\)0 s. 1 + 11 + Thật vậy, theo (7) suy ra xr = 0 nên JX()⊆ J với (J (X) = { jx:0j > } ). Ta cần 1 chứng minh hệ véc tơ {Ajj , ∈ J} độc lập tuyến tính. - 13 - Nguyễn Hải Đăng - Khoa KHCB&KTCS
Toán Kinh tế - Trường Cao Đẳng Công Nghiệp Nam Định α x = 0 1 + Thật vậy, giả sử ∑ jj ta cần chứng minh α j = 0∀∈j J . Ta có: jJ∈ 1 0(==∑∑αααααjjx jjxx +=+ ss ∑jj x s ∑ zA jsj = ∑ αααj −+ sjsj zAz) srsrA 1 jJ∈ jJjr∈≠,,jJ∈000 jr ≠ jJ∈ jJ∈ , jr ≠ Vì hệ véc tơ {Aj , jJ∈ 0} độc lập tuyến tính, nên: ⎪⎧ααjsjs−=∀∈≠zjJj0,0 r ⎨ ⎩⎪α srsz = 0 1 Vì zrs > 0 nên suy ra αs = 0 , do đó α j = 0(j ∈≠Jj , r ) . Vậy α j =∀∈0(\)j JJr =0 ∪s. + 11 + 1 1 Vì JX()⊆ J nên hệ véc tơ {Ajj ,(∈ J X)} độc lập tuyến tính. Do đó x là phương án cực biên và J1 là cơ sở của phương án cực biên x1. Như vậy x1 là một phương án cực biên của bài toán. ở trên ta đã tìm các thành phần của phương án cực biên mới x1 cùng với giá trị 1 hàm mục tiêu f(x ) thông qua các hệ số khai triển zjk và các ước lượng Δk trong cơ sở J. 11 1 Như vậy chúng ta cần xác định các đại lượng zvjk à Δk trong cơ sở J . Theo định nghĩa zjk 1 1 và z jk là các hệ số khai triển của véc tơ Ak tương ứng với cơ sở J, J . 1 Akj==∑∑zAkjj zAkj(*) 1 jJ∈ 0 jJ∈ ∑∑z jkAzAz j =+jk j rkA r ( ) jJ∈∈00jJ, jr≠ Từ 1 ⎛⎞ As ==zAjs j Azj js + zAv rs r ì zrs > 0 tac ó A r =−⎜⎟ As zAjs j . ∑∑ z ∑ jJ∈∈00jJ, jr≠ rs ⎝⎠jJ∈0 , jr≠ Thay vào ( ) ta được: zrk ∑∑z jkAzAAz j =+−jk j ()s ∑jsA j jJ∈∈00jJ,, jr≠zrs jJ∈0 jr≠ zzrk rk ( ) =−+∑ (.).zzAAjk js j s jJ∈≠0 , jr zzrs rs 1 Vì {Aj , jJ∈ } độc lập tuyến tính nên từ (*) và ( ) suy ra: - 14 - Nguyễn Hải Đăng - Khoa KHCB&KTCS
Toán Kinh tế - Trường Cao Đẳng Công Nghiệp Nam Định ⎧zrk ⎪ víi j = s 1 ⎪ zrs x jk = ⎨ z ⎪zz−∈rk víi j J,j ≠r ⎪ jk js 0 ⎩ zrs 11 ⎛⎞zzrk rk Δ=kj∑∑zckjkj − c =⎜⎟ zkj − zsjsk c + c − c 1 zz jJ∈ jJ∈≠0 , jr⎝⎠rs rs ⎛⎞zzrk rk ⎛⎞zrk =−⎜⎟zzcccvjk js j +− s k ( ì j=r thì ⎜⎟ zzjk −=js 0) ∑ zz z jJ∈ 0 ⎝⎠rs rs ⎝⎠rs zzrk rk =−−−=Δ−()∑∑zcjk j c k zcjs j c s ) k Δs . jJ∈∈00zzrs jJ rs 1 zrk Vậy: Δ=Δ−kk Δs. zrs Để thuận tiện cho việc tính toán, người ta sắp xếp các số liệu thành một bảng gọi là bảng đơn hình như dưới đây: Hệ Cơ Phương c1 c2 ck cs cn số sở án A1 A2 Ak As An c j A x z j 1 z j 2 z j k z sj z j n 1 j1 j1 1 1 1 1 1 c z z z j2 A x j21 j2 2 z j k z j s j2 n j2 j2 2 2 c j z j 1 z j 2 z j s r A x r r z j k z j s r jr jr r r c z z z jm jm1 jm 2 z z jm s A x jm k jm s jm jm f(x) Δ1 Δ 2 Δ k Δs Δn Các cột ứng với j∈J sẽ là các véc tơ đơn vị với số 1 trên dòng với chỉ số j. Chú ý: Đối với bài toán dạng chuẩn, ta có ngay một phương án cực biên 0 xbbb= (12 , , ,m ,0,0, ,0) với hệ véc tơ cơ sở tương ứng là A12,A , , Am . Đây là các véc tơ đơn vị nên trong bảng đơn hình xuất phát ta có zajkj= k - 15 - Nguyễn Hải Đăng - Khoa KHCB&KTCS
Toán Kinh tế - Trường Cao Đẳng Công Nghiệp Nam Định 1.2.4. Thuật toán đơn hình 0 Bước xuất phát: Tìm một phương án cực biên x và cơ sở J0 tương ứng . Tìm các hệ số khai triển zjk và các ước lượng Δ k . Bước 1: Kiểm tra dấu hiệu tối ưu 0 - Nếu Δ≤∀∉k 0 kJ0 thì x là phương án tối ưu. Thuật toán kết thúc. - Nếu ∃Δk > 0 thì chuyển sang bước 2. Bước 2: Kiểm tra dấu hiệu hàm mục tiêu giảm vô hạn: Với mỗi kJ∉ 0 mà Δ k > 0 thì kiểm tra các hệ số khai triênjk của cột Ak tương ứng: a) Nếu có một Δ k > 0 mà tất cả zjk ≤ 0 ∀j ∈ J0 thì kết luận hàm mục tiêu giảm vô hạn trên miền ràng buộc. Bài toán không có lời giải hữu hạn. Thuật toán kết thúc. b) ∀∉kJ0 mà Δk > 0 đều tồn tại ít nhất một hệ số zjk>0 thì chuyển sang bước 3 Bước 3: Xác định cột xoay, dòng xoay, phần tử trục - Chọn chỉ số s∉ J0 :Δ=skmax{ Δ> 0,kJ ∉0}, đánh dấu cột s là cột xoay 0 ⎧ 0 ⎫ xr ⎪ x j ⎪ - Tìm chỉ số r đạt min:θ ==min⎨ ,z js > 0⎬ , đánh dấu hàng r là hàng xoay. zzrs ⎪⎩⎭js ⎪ 1111 1 Bước 4: Tính các x j ,(),,f xΔkjkz trong cơ sở mới JJr= (\)0 ∪ stheo các công thức trên. Ghi nhận các kết quả trong một bảng mới. Quay trở lại bước 1. - Để nhận được bảng đơn hình mới từ bảng đơn hình cũ ta làm như sau: + Thay Ar bằng As, cr bằng cs. + Chia các phần tử trên hàng xoay ( hàng r) cho phần tử trục zrs ta được hàng r mới gọi là hàng chuẩn. + Mỗi phần tử khác ngoài hàng xoay trừ đi tích của phần tử cùng hàng với nó trên cột xoay với phần tử cùng cột với nó trên hàng chuẩn được phần tử cùng vị trí trong bảng đơn hình mới. a b ' bd a= a − c d c Dòng xoay Cột xoay - 16 - Nguyễn Hải Đăng - Khoa KHCB&KTCS
Toán Kinh tế - Trường Cao Đẳng Công Nghiệp Nam Định 1.2.5. Tính hữu hạn của thuật toán đơn hình Nếu bài toán quy hoạch tuyến tính có phương án và không suy biến thì sau hữu hạn bước lặp theo thủ tục đơn hình ta sẽ tìm thấy phương án tối ưu hoặc phát hiện ra bài toán có hàm mục tiêu giảm vô hạn hay bài toán không có lời giải hữu hạn. 0 0 xr 1 0 Thật vậy, vì bài toán không suy biến nên x j >∀∈0 jJnên θ =>0 suy ra x ≠ x, zrs f(x1) >0), bài toán “M” có dạng: - 17 - Nguyễn Hải Đăng - Khoa KHCB&KTCS
Toán Kinh tế - Trường Cao Đẳng Công Nghiệp Nam Định nm fx()=+∑∑ cxjj M xni+ →min ji==11 ⎧ n ⎪∑aij x j+ x n+ i = b i ( i = 1, m ) Với điều kiện: ⎨ j=1 ⎪ ⎩xj ≥0( j = 1, n ) Bài tóan “M” có ngay phương án cực biên xuất phát: 0 x = (0;0; ;0;bb12 ; ; bm ) với cơ sở J0 là: Em =(An+1; An+2; ;An+m); J0 = {n+1; n+2; ;n+m} Do vậy áp dụng được thuật toán đơn hình để giải bài toán “M”. Từ cách xây dựng bài toán “M” như trên ta thấy: 0 Nếu xxxx= (12 ; ; n ;0;0; ;0) là phương án của bài toán “M” thì x = (xx12 ; ; xn ) là phương án của bài toán ban đầu và ngược lại, đồng thời f(x) = f ()x . 1.3.3. Mối quan hệ giữa bài toán “M” và bài toán ban đầu * * • Nếu bài toán “M” có: xxxx= ( 12, , ,n ,0,0, ,0) là phương án tối ưu thì bài * * toán ban đầu có x = ( xx12, , , xn ) là phương án tối ưu và f ()xfx= (). * * • Nếu bài toán “M” có xxxx= ( 12, , , nm+ ) trong đó tồn tại ít nhất một * xni+ >=0(im 1, ) thì bài toán ban đầu không có phương án tối ưu (bài toán không giải được. • Nếu bài toán “M” vô nghiệm thì bài toán ban đầu cũng vô nghiệm. Chú ý: • Nếu bài toán ban đầu có nghiệm x* thì nghiệm này chỉ có thể nhận được sau ít nhất m + 1 bảng đơn hình. • Nếu ma trận hệ số a đã chứa m véc tơ đơn vị (mm< ) thì khi xây dựng bài ( ij )mn× 1 1 toán “M” chỉ cần thêm mm− 1 biến giả. • Nếu trong quá trình tính toán một biến giả xni+ (1,i= m) được đưa ra khỏi cơ sở thì sẽ không thể vào cơ sở được nữa. Do vậy việc tính toán trong bảng đơn hình tiếp theo đối với cột Ani+ (1,im= )là không cần thiết (chỉ đối với biến giả). - 18 - Nguyễn Hải Đăng - Khoa KHCB&KTCS
Toán Kinh tế - Trường Cao Đẳng Công Nghiệp Nam Định 1.4. Bài tập mẫu Bài 1: F(x) = 2x1 - x2 + 3x3 + 2x4 + x5 + 4x6 => MIN 2x2 + x3 + 2x4 + 3x5 = 6 x2 - x4 + x5 + x6 = 3 x1 -4x4 + 2x5 = 5 x1 >=0, x2 >=0, x3 >=0, x4 >=0, x5 >=0, x6 >=0 Ñaây laø baøi toaùn quy hoaïch tuyeán tính daïng chuaån, ta coù ngay phöông aùn cöïc bieân 0 x = (5,0,6,0,0,3) vôùi cô sôû A3, A6, A1. Töø ñoù ta coù baûng ñôn hình xuaát phaùt nhö sau: 2 -1 3 2 1 4 Heä soá Cô sôû P.a θ A1 A2 A3 A4 A5 A6 3 A3 6 0 2 1 2 3 0 2 4 A6 3 0 1 0 -1 1 1 3 2 A1 5 1 0 0 -4 2 0 5/2 F(x) 40 0 11 0 -8 16 0 Do coøn toàn taïi giaù trò Delta lôùn hôn 0 neân chöa coù phöông aùn toái öu ta caàn tìm bieán ñöa vaøo Coät coù giaù trò lôùn nhaát öùng vôùi A5.Vaäy bieán ñöa vaøo laø A5 Haøng coù giaù trò θ nhoû nhaát öùng vôùi coät ñoù laø haøng 1 2 -1 3 2 1 4 Heä soá Cô sôû P.a θ A1 A2 A3 A4 A5 A6 1 A5 2 0 2/3 1/3 2/3 1 0 3 4 A6 1 0 1/3 -1/3 -5/3 0 1 3 2 A1 1 1 -4/3 -2/3 -16/3 0 0 - F(x) 8 0 1/3 -16/3 -56/3 0 0 Do coøn toàn taïi giaù trò Delta lôùn hôn 0 neân chöa coù phöông aùn toái öu ta caàn tìm bieán ñöa vaøo Coät coù giaù trò lôùn nhaát öùng vôùi A2.Vaäy bieán ñöa vaøo laø A2 Haøng coù giaù trò θ nhoû nhaát öùng vôùi coät ñoù laø haøng 1 2 -1 3 2 1 4 Heä soá Cô sôû P.a θ A1 A2 A3 A4 A5 A6 -1 A2 3 0 1 1/2 1 3/2 0 - - 19 - Nguyễn Hải Đăng - Khoa KHCB&KTCS
Toán Kinh tế - Trường Cao Đẳng Công Nghiệp Nam Định 4 A6 0 0 0 -1/2 -2 -1/2 1 - 2 A1 5 1 0 0 -4 2 0 - F(x) 7 0 0 -11/2 -19 -1/2 0 Phöông aùn toái öu cuûa baøi toaùn laø: (5,3,0,0,0,0) Giaù trò haøm muïc tieâu ñaït ñöôïc laø: F(x) = 7 Bài 2: F(x) = 2x1 - x2 + 3x3 + 2x4 + x5 + 4x6 => MIN 2x1 + x3 + 2x4 + 3x5 + x6 = 6 x2 - x4 + x5 = 3 - x1 + 3x4 + 2x5 ≤ 5 x1 >=0, x2 >=0, x3 >=0, x4 >=0, x5 >=0, x6 >=0 Chuyeån veà daïng chính taéc: F(x) = 2x1 - x2 + 3x3 + 2x4 + x5 + 4x6 => MIN 2x1 + x3 + 2x4 + 3x5 + x6 = 6 x2 - x4 + x5 = 3 - x1 + 3x4 + 2x5 + x7 = 5 x7 laø bieán phuï x1 >=0, x2 >=0, x3 >=0, x4 >=0, x5 >=0, x6 >=0, x7 >=0 Ñaây laø baøi toaùn quy hoaïch tuyeán tính daïng chuaån, ta coù ngay phöông aùn cöïc bieân 0 x = (0,3,0,0,0,6,5) vôùi cô sôû A6, A2, A7. Töø ñoù ta coù baûng ñôn hình xuaát phaùt nhö sau: Heä Cô 2 -1 3 2 1 4 0 P.a θ soá sôû A1 A2 A3 A4 A5 A6 A7 4 A6 6 2 0 1 2 3 1 0 2 -1 A2 3 0 1 0 -1 1 0 0 3 0 A7 5 -1 0 0 3 2 0 1 5/2 F(x) 21 6 0 1 7 10 0 0 Do coøn toàn taïi giaù trò Delta lôùn hôn 0 neân chöa coù phöông aùn toái öu ta caàn tìm bieán ñöa vaøo Coät coù giaù trò lôùn nhaát öùng vôùi A5.Vaäy bieán ñöa vaøo laø A5 Haøng coù giaù trò θ nhoû nhaát öùng vôùi coät ñoù laø haøng 1 Heä Cô 2 -1 3 2 1 4 0 P.a θ soá sôû A1 A2 A3 A4 A5 A6 A7 1 A5 2 2/3 0 1/3 2/3 1 1/3 0 3 -1 A2 1 -2/3 1 -1/3 -5/3 0 -1/3 0 - 0 A7 1 -7/3 0 -2/3 5/3 0 -2/3 1 3/5 F(x) 1 -2/3 0 -7/3 1/3 0 -10/3 0 Do coøn toàn taïi giaù trò Delta lôùn hôn 0 neân chöa coù phöông aùn toái öu ta caàn tìm bieán ñöa vaøo Coät coù giaù trò lôùn nhaát öùng vôùi A4.Vaäy bieán ñöa vaøo laø A4 Haøng coù giaù trò θ nhoû nhaát öùng vôùi coät ñoù laø haøng 3 Heä Cô P.a 2 -1 3 2 1 4 0 θ - 20 - Nguyễn Hải Đăng - Khoa KHCB&KTCS
Toán Kinh tế - Trường Cao Đẳng Công Nghiệp Nam Định soá sôû A1 A2 A3 A4 A5 A6 A7 1 A5 8/5 8/5 0 3/5 0 1 3/5 -2/5 - -1 A2 2 -3 1 -1 0 0 -1 1 - 2 A4 3/5 -7/5 0 -2/5 1 0 -2/5 3/5 - F(x) 4/5 -1/5 0 -11/5 0 0 -16/5 -1/5 Phöông aùn toái öu cuûa baøi toaùn laø: (0,2,0,3/5,8/5,0,0) Giaù trò haøm muïc tieâu ñaït ñöôïc laø: F(x) = 4/5 Bài 3: F(x) = x1 - x2 + 3x3 + 2x4 + x5 + 2x6 => MIN 2x1 + x3 + 2x4 + 3x5 + x6 = 6 2x1 + x2 -2x3 - x4 + x5 = 4 x1 + 3x2 + 3x4 + 2x5 = 7 x1 >=0, x2 >=0, x3 >=0, x4 >=0, x5 >=0, x6 >=0 Vì baøi toaùn chöa ôû daïng chuaån neân ta ñöa vaøo hai bieán giaû x7, x8. Khi ñoùù baøi toaùn “M” coù daïng: F(x) = x1 - x2 + 3x3 + 2x4 + x5 + 2x6 + Mx7 + Mx8 => MIN 2x1 + x3 + 2x4 + 3x5 + x6 = 6 2x1 + x2 -2x3 - x4 + x5+ x7 = 4 x1 + 3x2 + 3x4 + 2x5 + x8 = 7 x7, x8 laø bieán giaû x1 >=0, x2 >=0, x3 >=0, x4 >=0, x5 >=0, x6 >=0 x7>=0, x8>=0 0 Baøi toaùn coù ngay phöông aùn cöïc bieân x = (0,0,0,0,0,6,4,7) vôùi cô sôû A6, A7, A8. Töø ñoù ta coù baûng ñôn hình xuaát phaùt nhö sau: Heä Cô 1 -1 3 2 1 2 M M P.a θ soá sôû A1 A2 A3 A4 A5 A6 A7 A8 2 A6 6 2 0 1 2 3 1 0 0 - M A7 4 2 1 -2 -1 1 0 1 0 4 M A8 7 1 3 0 3 2 0 0 1 7/3 F(x) 12 3 1 -1 2 5 0 0 0 +11M +3M +4M -2M +2M +3M Do coøn toàn taïi giaù trò Delta lôùn hôn 0 neân chöa coù phöông aùn toái öu ta caàn tìm bieán ñöa vaøo Coät coù giaù trò lôùn nhaát öùng vôùi A2.Vaäy bieán ñöa vaøo laø A2 Haøng coù giaù trò θ nhoû nhaát öùng vôùi coät ñoù laø haøng 3 Cô 1 -1 3 2 1 2 M M Heä soá P.a θ sôû A1 A2 A3 A4 A5 A6 A7 A8 2 A6 6 2 0 1 2 3 1 0 3 M A7 5/3 5/3 0 -2 -2 1/3 0 1 1 -1 A2 7/3 1/3 1 0 1 2/3 0 0 7 F(x) 29/3 8/3 0 -1 1 13/3 0 0 +5/3M +5/3M -2M -2M +1/3M Do coøn toàn taïi giaù trò Delta lôùn hôn 0 neân chöa coù phöông aùn toái öu ta caàn tìm bieán ñöa vaøo Coät coù giaù trò lôùn nhaát öùng vôùi A1.Vaäy bieán ñöa vaøo laø A1 Haøng coù giaù trò θ nhoû nhaát öùng vôùi coät ñoù laø haøng 2 - 21 - Nguyễn Hải Đăng - Khoa KHCB&KTCS
Toán Kinh tế - Trường Cao Đẳng Công Nghiệp Nam Định Heä Cô 1 -1 3 2 1 2 M M P.a θ soá sôû A1 A2 A3 A4 A5 A6 A7 A8 2 A6 4 0 0 17/5 22/5 13/5 1 10/11 1 A1 1 1 0 -6/5 -6/5 1/5 0 - -1 A2 2 0 1 2/5 7/5 3/5 0 10/7 F(x) 7 0 0 11/5 21/5 19/5 0 Do coøn toàn taïi giaù trò Delta lôùn hôn 0 neân chöa coù phöông aùn toái öu ta caàn tìm bieán ñöa vaøo Coät coù giaù trò lôùn nhaát öùng vôùi A4.Vaäy bieán ñöa vaøo laø A4 Haøng coù giaù trò θ nhoû nhaát öùng vôùi coät ñoù laø haøng 1 Heä Cô 1 -1 3 2 1 2 M M P.a θ soá sôû A1 A2 A3 A4 A5 A6 A7 A8 2 A4 10/11 0 0 17/22 1 13/22 5/22 20/13 1 A1 23/11 1 0 -3/11 0 10/11 3/11 23/10 -1 A2 8/11 0 1 -15/22 0 -5/22 -7/22 - F(x) 35/11 0 0 -23/22 0 29/22 -21/22 Do coøn toàn taïi giaù trò Delta lôùn hôn 0 neân chöa coù phöông aùn toái öu ta caàn tìm bieán ñöa vaøo Coät coù giaù trò lôùn nhaát öùng vôùi A5.Vaäy bieán ñöa vaøo laø A5 Haøng coù giaù trò θ nhoû nhaát öùng vôùi coät ñoù laø haøng 1 1 -1 3 2 1 2 M M Heä soá Cô sôû P.a θ A1 A2 A3 A4 A5 A6 A7 A8 1 A5 20/13 0 0 17/13 22/13 1 5/13 - 1 A1 9/13 1 0 -19/13 -20/13 0 -1/13 - -1 A2 14/13 0 1 -5/13 5/13 0 -3/13 - F(x) 15/13 0 0 -36/13 -29/13 0 -19/13 Phöông aùn toái öu cuûa baøi toaùn môû roäng laø: (9/13,14/13,0,0,20/13,0,0,0) Do ñoù phöông aùn toái öu cuûa baøi toaùn ban ñaàu laø: (9/13,14/13,0,0,20/13,0) Giaù trò haøm muïc tieâu ñaït ñöôïc laø: F(x) = 15/13 - 22 - Nguyễn Hải Đăng - Khoa KHCB&KTCS
Toán Kinh tế - Trường Cao Đẳng Công Nghiệp Nam Định BÀI TẬP CHƯƠNG I 1. Lập bài toán QHTT Bài 1. Xí nghiệp sản xuất giấy có 3 phân xưởng. Do trang bị kỹ thuật khác nhau nên mức hao phí tre gỗ, axít để sản xuất một tấn giấy thành phẩm cũng khác nhau. Mức hao phí được cho trong bảng dưới đây: Mức hao phí nguyên liệu cho một tấn giấy Nguyên liệu P.Xưởng I P.Xưởng II P.Xưởng III Tre gỗ 1,4 tấn 1,3 1,2 Axít 0,1 0,12 0,15 Số lượng tre gỗ có trong năm là 1.500.000 tấn, Axít là 100.000 tấn. Hãy lập kế hoạch sản xuất sao cho tổng số giấy sản xuất trong năm của xí nghiệp là nhiều nhất? Bài 2. Một xí nghiệp có thể sản xuất bốn loại mặt hàng xuất khẩu H1, H2, H3, H4. Để sản xuất 4 loại mặt hàng này, xí nghiệp sử dụng hai loại nguyên liệu N1, N2. Số nguyên liệu tối đa mà xí nghiệp huy động được tương ứng là 600 kg và 800 kg. Mức tiêu hao mỗi loại nguyên liệu để sản xuất một mặt hàng và lợi nhuận thu được cho trong bảng sau: Định mức tiêu hao H1 H2 H3 H4 nguyên liệu và lợi nhuận N1 0,5 0,2 0,3 0,4 N2 0,1 0,4 0,2 0,5 Lợi nhuận 0,8 0,3 0,5 0,4 Lập mô hình sao cho xí nghiệp sản xuất đạt lợi nhuận cao nhất? Bài 3. Một cơ sở dự định sản xuất tối đa trong một ngày 500 ổ bánh mì dài và 500 ổ bánh mì tròn, muốn đạt lợi nhuận nhiều nhất, với những điều kiện như sau: - Gía bán một ổ bánh mì dài làm từ 400g bột là 325 đồng, một ổ bánh mì tròn làm từ 250g bột là 220 đồng. - Số lượng bột được cung cấp tối đa trong ngày là 225 kg với giá mỗi kg là 300 đồng. - 23 - Nguyễn Hải Đăng - Khoa KHCB&KTCS
Toán Kinh tế - Trường Cao Đẳng Công Nghiệp Nam Định - Lò nướng bánh cho phép nướng 75 ổ bánh mì dài hay 100 ổ bánh mì tròn trong một giờ nhưng không thể nướng hai loại cùng một lúc. Lò nướng hoạt động tối đa 8 giờ trong một ngày. Hãy lập mô hình cho bài toán nêu trên? Bài 4. Ba xí nghiệp A, B, C cùng có thể sản xuất áo vét và quần. Khả năng sản xuất của môic xí nghiệp như sau: Khi đầu tư 1000 USD vào xí nghiệp A thì thu được 35 áo vét và 45 quần; vào xí nghiệp B thì thu được 40 áo vét và 42 quần; vào xí nghiệp C thì thu được 43 áo vét và 30 quần. Lượng vải và giờ công sản xuất được cho trong bảng sau: A B C Xí nghiệp vải/giờ vải/giờ vải/giờ 1 áo vét 3,5m 20giờ 4m 16giờ 3,8m 18giờ 1 quần 2,8m 10giờ 2,6m 12giờ 2,5m 15giờ Tổng số vải huy động được là 10000m. Tổng số giờ công huy động được là 52000 giờ. Theo hợp đồng thì tối thiểu phải có 1500 bộ quần áo, nếu lẻ bộ thì quần là dễ bán Hãy lập kế hoạch đầu tư vào mỗi xí nghiệp bao nhiêu vốn để: - Hoàn thành hợp đồng. - Không khó khăn về tiêu thụ. - Không bị động về vải và giờ công. - Tổng số vốn đầu tư là nhỏ nhất. Bài 5. Một nhà máy cán thép có thể sản xuất hai loại sản phẩm: thép tấm và thép cuộn. Nếu chỉ sản xuất một loại sản phẩm thì nhà máy chỉ có thể sản xuất 200 tấn thép tấm hoặc 140 tấn thép cuộn trong một giờ. Lợi nhuận thu được khi bán một tấn thép tấm là 25 USD, một tấn thép cuộn là 30 USD. Nhà máy làm việc 40 giờ trong một tuần và thị trường tiêu thụ tối đa là 6000 tấn thép tấm và 4000 tấn thép cuộn. Nhà máy cần sản xuất mỗi loại sản phẩm là bao nhiêu trong một tuần để lợi nhuận thu được là cao nhất? 2. Chuyển bài toán về dạng chính tắc Bài 1. −−543xxx123 − → min ⎧23xxx123++≤ 5 ⎪ ⎪42xx12++ x 3 ≤11 với điều kiện ⎨ ⎪34xxx123++≤ 2 8 ⎪ ⎩xxx123,,≥ 0 - 24 - Nguyễn Hải Đăng - Khoa KHCB&KTCS
Toán Kinh tế - Trường Cao Đẳng Công Nghiệp Nam Định Bài 2. 2mxx12−→ax Bài 4. xxx123− 2m−→in ⎧−−xxx12322 + ≤ ⎧xx12+−22 x 3 ≤ ⎪ ⎪ ⎪22xxx123−−≥ xxx123−+≤1 với điều kiện⎨ ⎪ xxx++=4 xx5 ⎪ 123 với điều kiện⎨ 23+≤ ⎪xx,0≥ ⎪22xx−≤ ⎩ 13 ⎪ 12 ⎩⎪xx13,0≥ Bài 3. 24xxxx1+++→ 234 max Bài 5. −5210/3mixx12−− x 3 →n ⎧xxx124++≤34 ⎧2434xxx123++≥6 ⎪ ⎪ ⎪23xx12+≤ ⎪4233xxx135++≤8 với điều kiện⎨ với điều kiện⎨ xxx++≤43 321xx+≤ ⎪ 234 ⎪ 13 ⎪ ⎪ ⎩xxx124,,≥ 0 ⎩xx13,0≥ 3. Giải bài toán bằng phương pháp đơn hình Bμi 1: f ()xxxxxxx=−+++−→123456 Min Bμi 4: ⎧xx14++69 x 6 = F(x) = x1 - x2 + x3 + x4 + x5 => MIN ⎪ ⎪342xx12+− x 3 + x 6 =2 ⎨ 2x1 + x3 - x4 + 3x5 = 6 xxxx1356+++=226 ⎪ -x1 - x4 + x6 = -3 ⎪xj≥∀=01,6 ⎩ j 4x1 + x2 - x4 + 2x5 = 5 x1 >=0, x2 >=0, x3 >=0, x4 >=0, x5 >=0, x6 >=0 Bμi 2: f ()xxxxxx=−+++→12345 Min Bμi 5: f (xM )= 4x1234 +3x +2x +3x → in ⎧26xxx134++= ⎧x+x-x134= 2 ⎪25xxx++= 3 ⎪ ⎪ 145 ⎪-x124 +2x +x≤ 11 ⎨ ⎨ 45xxx124++= 2x +x -3x≥ 8 ⎪ ⎪ 12 4 ⎪ ⎩xjj ≥∀=01,5 ⎪ ⎩x0(1,4)j ≥=j Bμi 3: f ()xx=+12345 3 xxx −++ 4 x → Min Bμi 6: F(x) = 2x1 - x2 + 3x3 + 2x4 + x5 + 4x6 => MIN ⎧xxxx1345+++=278 ⎪ x3 + 2x4 + 3x5 + x6 = 6 ⎪xxx345++≤57 x - x + x = 3 ⎨xxxx++−=423 2 4 5 ⎪ 2345 x -4x -4x + 2x = 5 ⎪ 1 3 4 5 xjj ≥∀=01,5 ⎩ x1 >=0, x2 >=0, x3 >=0, x4 >=0, x5 >=0, x6 >=0 §S: (1,1,0,0,2,0) - 25 - Nguyễn Hải Đăng - Khoa KHCB&KTCS
Toán Kinh tế - Trường Cao Đẳng Công Nghiệp Nam Định Bμi 7: F(x) = 4x1 - x2 + 2x3 + 2x4 + 3x6 => MIN 2x1 + x2 + x3 + 3x5 = 8 x1 + 3x3 - x5 + x6 = 9 x1 + 2x3 + x4 -2x5 = 3 x1 >=0, x2 >=0, x3 >=0, x4 >=0, x5 >=0, x6 >=0 §S: (0,13/2,3/2,0,0,9/2) Bμi 8: F(x) = x1 - x2 + x3 + x4 + 3x5 + 2x6 => MIN 2x1 -2x3 + x5 + 4x6 = 7 x1 + x2 + 3x3 + x6 = 9 x1 + 2x3 + x4 -2x6 = 5 x1 >=0, x2 >=0, x3 >=0, x4 >=0, x5 >=0, x6 >=0 §S: (7/2,11/2,0,3/2,0,0) Bμi 9: F(x) = x1 - x2 - x3 + x5 + 3x6 => MIN 2x1 + x2 + x5 + 3x6 = 3 2x1 - 3x2 - x4 + x6 =-8 x1 + x3 + 3x6 = 6 x1 >=0, x2 >=0, x3 >=0, x4 >=0, x5 >=0, x6 >=0 Bμi 10: F(x) = 2x1 - x2 + 3x3 + 2x4 => MIN 2x1 + x2 + x3 + 2x4 + x5 = 5 2x1 -2x3 - x4 + x5 + x6 = 3 x1 + x3 + 3x4 + 2x5 + 2x6 = 8 x1 >=0, x2 >=0, x3 >=0, x4 >=0, x5 >=0, x6 >=0 (0,23/5,2/5,0,0,19/5,0,0) Bμi 11: F(x) = 2x1 - x2 + 3x3 + 2x4 => MIN 2x1 + x2 + 3x3 + 2x4 = 6 5x1 + x3 - x4 = 5 x1 + 2x2 + x3 + x4 = 4 x1 >=0, x2 >=0, x3 >=0, x4 >=0 (17/22,23/22,25/22,0) Bμi 12: F(x) = x1 - x2 + x3 + x4 => MIN x1 + x2 + x3 + 2x4 = 5 5x1 + 2x2 + x3 + x4 = 3 x1 + 2x2 + x3 + x4 = 1 x1 >=0, x2 >=0, x3 >=0, x4 >=0 (VN) - 26 - Nguyễn Hải Đăng - Khoa KHCB&KTCS
Toán Kinh tế - Trường Cao Đẳng Công Nghiệp Nam Định Chương 2. BÀI TOÁN ĐỐI NGẪU 2.1. Bài toán gốc và cách thành lập bài toán đối ngẫu 2.1.1. Định nghĩa Định nghĩa 1: Cho bài toán quy hoạch tuyến tính dạng chính tắc: n f(x)=®åcxjj min j=1 ìn ïåaijxji==b,im1, với điều kiện íj=1 ()I ï îxj³=0,jn1, Ta xây dựng bài toán quy hoạch tuyến tính khác có dạng, kí hiệu là ()I% như sau: m °f(y)=®åbyii max i=1 ìn ïåaijyij£=c,jn1, với điều kiện íj=1 ()I% ï îyiÎ=¡,im1, Bài toán ()I% được gọi là bài toán quy hoạch tuyến tính đối ngẫu của bài toán ()I . Ta cũng chứng minh được bài toán ()I là bài toán đối ngẫu của bài toán ()I% , do vậy cặp bài toán ()I và ()I% được gọi là cặp bài toán đối ngẫu. Định nghĩa 2 : Ta gọi hai ràng buộc bất dẳng thức (kể cả ràng buộc về dấu) trong hai bài toán cùng tương ứng với một chỉ số là một cặp ràng buộc đối ngẫu. Như vậy, ta có n + m cặp ràng buộc đối ngẫu: n · xj³0Ûåaijyij£=c,jn1, j=1 n · åaijxj=biiÛyÎ=¡,im1, j=1 Từ quan hệ giữa hai bài toán, ta có những nguyên tắc thành lập bài toán đối ngẫu như sau : +) Nếu fx()® min thì °f(yM)® ax và nếu f(xM)® ax thì °f()y® Min . +) Số ràng buộc trong bài toán này bằng số biến trong bài toán kia +) Hệ số hàm mục tiêu trong bài toán này là vế phải hệ ràng buộc trong bài toán kia. +) Ma trận hệ số trong hai bài toán là chuyển vị của nhau. Nguyễn Hải Đăng - Khoa KHCB&KTCS - 30 -
Toán Kinh tế - Trường Cao Đẳng Công Nghiệp Nam Định Trên đây là phương pháp xác định bài toán đối ngẫu của bài toán quy hoạch tuyến tính dạng chính tắc, còn đối với bài toán quy hoạch tuyến tính dạng tổng quát thì ta làm thế nào? Đối với bài toán quy hoạch tuyến tính tổng quát, ta đưa bài toán về dạng chính tắc, xây dựng bài toán đối ngẫu của bài toán này và gọi nó là bài toán đối ngẫu của bài toán đã cho. Từ đó ta có các cặp bài toán đối ngẫu như sau: Bài toán gốc Bài toán đối ngẫu n m f(x)=®cx min åjj °f(y)=®åbyii max j=1 i=1 ìn ìn ax³=b,im1, ay£=c,jn1, ïåijji ()II ïåijij íj=1 íj=1 ()II° ïx³=0,jn1, ï îj îyi³=0,im1, n m f(x)=®cx min åjj °f(y)=®åbyii max j=1 i=1 n ì ' ìn aijxji£=b,im1, ay³=c,jn1, ' ïå ()II ïåijij ± íj=1 íj=1 ()II ïx£=0,jn1, ï îj îyi£=0,im1, Các cặp ràng buộc đối ngẫu của ()II và ()II° là : n · xj³0Ûåaijyij£=c,jn1, j=1 n · åaijxj³biiÛy³=0,im1, j=1 của ()II ' và ()II±' là : n · xj£0Ûåaijyij³=c,jn1, j=1 n · åaijxj£biiÛy£=0,im1, j=1 Từ đó ta có phương pháp để xây dựng bài toán đối ngẫu trong trường hợp tổng quát như sau : Nguyễn Hải Đăng - Khoa KHCB&KTCS - 31 -
Toán Kinh tế - Trường Cao Đẳng Công Nghiệp Nam Định 2.1.2. Cách thành lập bài toán đối ngẫu Để xây dựng bài toán đối ngẫu, ta có thể tuân thủ theo các quy tắc được cho trong bảng sau : Bài toán gốc Bài toán đối ngẫu n m f(x)=®åcxjj min °f(y)=®åbyii max j=1 i=1 n ax=Îb,iI å ijji 1 yi ÎÎ¡,iI1 j=1 n å aijxji³Îb,iI2 yi ³Î0, iI2 j=1 n å aijxji£Îb,iI3 yi £Î0, iI3 j=1 n yj ³Î0, jJ1 å aijyij£Îc, jJ1 j=1 n yj £Î0, jJ2 å aijyij³Îc, jJ2 j=1 n yj ÎÎ¡, jJ3 å aijyij=Îc, jJ3 j=1 2.2. Các tính chất và ứng dụng của cặp bài toán đối ngẫu 2.2.1. Tính chất của cặp bài toán đối ngẫu Định lý 1: Với mọi cặp phương án x và y của cặp bài toán đối ngẫu, ta có : · Nếu fx()® min và °f(yM)® ax thì f(x)³ °fy() · Nếu f(xM)® ax và °f()y® Min thì f(x)£ °fy() Định lý 2: Nếu x*, y* lần lượt là phương án của một cặp bài toán đối ngẫu, thoả mãn CX* = Y*b thì x*, y* lần lượt là phương án tối ưu của mỗi bài toán. Như vậy đối với một cặp bài toán đối ngẫu, bao giờ cũng chỉ xảy ra một trong ba trường hợp sau: +) Nêú hai bài toán cùng không có phương án thì hiển nhiên cả hai bài toán đều không giải được. +) Nếu cả hai bài toán đều có phương án thì cả hai bài toán đều giải được. Khi đó mọi cặp phương án tối ưu x*, y*, ta luôn có +) Nếu một trong hai bài toán không có phương án thì bài toán còn lại nếu có phương án thì cũng không có phương án tối ưu. Nguyễn Hải Đăng - Khoa KHCB&KTCS - 32 -
Toán Kinh tế - Trường Cao Đẳng Công Nghiệp Nam Định Định lý 3 (Định lý độ lệch bù): Điều kiện cần và đủ để hai phương án x, y của một cặp bài toán đối ngẫu tối ưu là trong các cặp ràng buộc đối ngẫu nếu một ràng buộc thỏa mãn với dấu bất đẳng thức thực sự (lỏng) thì ràng buộc kia phải thỏa mãn với dấu bằng (chặt) và ngược lại. Hệ quả: Nếu một ràng buộc là lỏng đối với một phương án tối ưu của bài toán này thì ràng buộc đối ngẫu của nó phải là chặt đối với phương án tối ưu của bài toán kia. 2.2.2. Ứng dụng Kiểm tra một phương án hay một cặp phương án có tối ưu hay không ? · Nếu biết cặp phương án x* và y*, thì ta chỉ cần kiểm tra điều kiện f(x )= °fy(). · Nếu chỉ biết phương án x* thì áp dụng định lý độ lệch bù. Ví dụ 1. Cho bài toán: f(x)=76x1+x2-12x34+®xMax ì2x1-23x2-xx34+=28 ï ï 3x2+2xx34-21£- với điều kiện í2x-+=3xx10 ï134 ï îxjj ³=0,1,4 a) Lập bài toán đối ngẫu của bài toán trên và xác định các cặp ràng buộc đối ngẫu. b) Chứng tỏ xy00=(0,6,0,10);=-(3,0,7) là phương án tối ưu của cặp bài toán đối ngẫu. Ví dụ 2. Cho bài toán : f(x)=-x1-5x2+42x34-®xMin ì x1-4x2+xx34-£613 ï ï x1+2xx24+£39 í-3x-x-xx+£28 ï 1234 ï îxjj ³=0,1,4 Dùng tính chất của bài toán đối ngẫu, chứng tỏ bài toán trên giải được. Ví dụ 3. Cho bài toán : f(x)=5x1-9x2+15xx3++®7645xMin ìx1+31x2-x3-xx45+£ ï ï4x1+x3+24xx45-= với điều kiện và véc tơ x = (0,1,0,2,0) í-x-x+xx-21³- ï 1235 ï îxj ³0,jx=Î2,5, 1¡ a) Viết bài toán đối ngẫu b) Phương án x có là phương án tối ưu không ? Ví dụ 4. Cho bài toán quy hoạch tuyến tính: f(x)=-2x1-6x2+5x3-xx45-®4max Nguyễn Hải Đăng - Khoa KHCB&KTCS - 33 -
Toán Kinh tế - Trường Cao Đẳng Công Nghiệp Nam Định ìx1-+4x22x3-5xx45+=93 ï ïx2-3x3+4xx45-=56 với điều kiện íx-x+xx-=1 ï2345 ï îxjj ³=0,1,5 Chứng minh x* = (0,0,16,31,14) là phương án tối ưu ? 2.3. Phương pháp đơn hình đối ngẫu 2.3.1. Cơ sở lý luận Xét bài toán chính tắc f()x=®cxMin °f()y=®ybMax (I) Axb= ta có bài toán đối ngẫu ()I m ì % ìy Î¡ í í îx³ 0 îyAc£ 0 Giả sử x là phương án cực biên ứng với hệ véc tơ cơ sở {Aj } , Aj là các véc tơ đơn jJÎ 0 * *1- vị, Jm0 = , B==( AEj ) , cc= ()jjJÎ . Xét y= cB 0 jJÎ 0 *1- * Mệnh đề 1: y= cB là phương án của bài toán ()I% khi và chỉ khi Dk=c.0zkk-ck£" *-1* Thật vậy, y là phương án yA£cÛyAk£ck"kÛyAk=cBAk=czk£ckkkÛD£"0 Mệnh đề 2: Nếu tại phương án y= cB*1- có được x=³Bb-1 0 thì x là phương án tối ưu ( và y cũng là phương án tối ưu). Nếu x không thỏa mãn x=³Bb-1 0 thì x được gọi là một giả phương án tối ưu ( Vì có Dk £"0 k ) -1 -1 Rõ ràng Ax==ABbb và x ³ 0 nên x= Bb là phương án. Kết hợp Dk £"0 k nên x là phương án tối ưu. -1 , Mệnh đề 3: Kí hiệu ωi là hàng thứ i của ma trận B . Với mọi β ³ 0 , ta có yy=-βωi là phương án của ()I% khi và chỉ khi Dj£"β zjij , ,, Thật vậy, yy=-βωi là phương án của ()I% ÛyA=yA-βωiAÛyAjj=-yAAβωij , * 1*1 ÛyAj=-yAAjβωijÛcBAj-βωiAj£cjÛcBAj-ββzij£cjÛDj+cj-£zcijj ÛD£jβ zij *1- -1 Mệnh đề 4: Nếu tại phương án y= cB tồn tại xs < 0 ( trong x= Bb) và zjsj ³"0 thì hàm mục tiêu của bài toán ()I% không bị chặn trên tập phương án, do đó ()I% không có phương án tối ưu và (I) cũng vô nghiệm. Nguyễn Hải Đăng - Khoa KHCB&KTCS - 34 -
Toán Kinh tế - Trường Cao Đẳng Công Nghiệp Nam Định , , Thật vậy, Với mọi β ³ 0 , ta có yy=-βωi thỏa mãn Dj£β zij"Þjy là phương án. ,, Þ°f(y)=yb=yb-βωib=°f(y)-βωβiib=-°f()yx Lấy i = r và cho β®+¥Þ°fy(), ®+¥ . Vậy hàm mục tiêu không bị chặn nên bài toán không giải được. *1- -1 Mệnh đề 5: Nếu tại phương án y= cB tồn tại xr < 0 ( trong x= Bb) và tồn tại zrj < 0 thì ìü , ïïDj D s xây dựng được phương án mới yy=-βω0 r (trong đó β0 =Minzíý:0rj <=) tốt zz îþïïrjrs hơn y. -1 , Thật vậy, Giả sử trong x= Bb có xr < 0 . Ta có y là phương án Û³β 0 và ìü D j ïïDj D s Dj£"β zjij Û0£β£"j=<1,nz;0rj . Chọn β0 =Minzíý:0rj <=, hiển nhiên z zz rj îþïïrjrs ,,, β0 ³ 0 và °f(y)=yb=yb-β0 xr³yb=Þ°f()yy tốt hơn y. Đpc/m Từ các mệnh đề trên ta có nhận xét: · Nếu có nhiều xr < 0 thì ta có thể chọn xrr=<Minx{ 0} . khi đó véc tơ Ar được đưa ra khỏi cơ sở. ìü ïïDj D s · Từ việc chọn β0 =Minzíý:0rj <= nên ta sẽ đưa As vào cơ sở của zz îþïïrjrs phương án mới. · Việc xây dựng các số liệu khác được biến đổi trên bảng đơn hình như cách thông thường. · Cấu trúc bảng đơn hình đối ngẫu: Giả c12c ccmn Hệ số Cơ sở p.a A12A AAmn c1 A1 x1 1 0 0 c2 A2 x2 0 1 0 cm Am xm 0 0 1 D1 D2 Dm Dn Nguyễn Hải Đăng - Khoa KHCB&KTCS - 35 -
Toán Kinh tế - Trường Cao Đẳng Công Nghiệp Nam Định 2.3 2. Thuật toán đơn hình đối ngẫu 1. Bước xuất phát: Xuất phát từ hệ m véc tơ độc lập tuyến tính, có B==( AEj ) , jJÎ 0 * 1- Jm0 = sao cho Dj £"0 j . Tìm xx=(,0), trong đó x=cB;0xj="ÏjJ0 Lập bảng đơn hình đối ngẫu xuất phát. 2. Bước 1: Kiểm tra dấu hiệu tối ưu · Nếu xj ³0 "ÎjJ0 thì x tối ưu · Ngược lại, chuyển sang bước 2 3. Bước 2: Kiểm tra dấu hiệu bài toán vô nghiệm · Nếu $ MIN Các ràng buộc 2x1 + x3 - x4 + 3x5 = 6 -x1 - x4 + x6 = -3 4x1 + x2 - x4 + 2x5 = 5 x1 >=0, x2 >=0, x3 >=0, x4 >=0, x5 >=0, x6 >=0 Giải: Xét cơ sở A1 , A5 , A2 gồm các véc tơ đơn vị. Khi đó ma trận B là: æö100 ç÷ B = ç÷010 , ta được giả phương án x = (0, 5, 6, 0, 0, -3) ç÷ èø001 Bảng đơn hình đối ngẫu là: Nguyễn Hải Đăng - Khoa KHCB&KTCS - 36 -
Toán Kinh tế - Trường Cao Đẳng Công Nghiệp Nam Định Giả 1 -1 1 1 1 0 Hệ số Cơ sở p.a A1 A2 A3 A4 A5 A6 1 A3 6 2 0 1 -1 3 0 0 A6 -3 -1 0 0 -1 0 1 -1 A2 5 4 1 0 -1 2 0 -3 0 0 -1 0 0 1 A3 9 3 0 1 0 3 -1 1 A4 3 1 0 0 1 0 -1 -1 A2 8 5 1 0 0 2 -1 4 -2 0 0 -1 0 0 Tại đây ta có x³ 0 nên nó là phương án tối ưu, với x = (0, 8, 9, 3, 0, 0) và fMin = 4 Bài 2: f(x)=®4x1+3x234+2x+3x Min ìx1+x34-x2= ï ï-x1+2x24+x£11 í2x+x-3x8³ ï 124 ï îxj³=0(j 1,4) Đưa về dạng chính tắc: f(x)=®4x1+3x234+2x+3x Min ìx1+x34-x2= ìx1+x34-x2= ï ï ï-x1+2x2+x45+=x 11 ï-x1+2x2+x45+=x 11 Û í2x+x-3x8-=x í-2x-x+3x8+x =- ï 1246 ï 1246 ï ï îxj³=0(j 1,6) îxj³=0(j 1,6) 0 Như vậy ta có hệ véc tơ cơ sở {A3,A56,Ax} Þ=-(0,0,2,0,11,8) (luôn có Dk £"0 k ). Ta có bảng đơn hình đối ngẫu: Giả 1 -1 1 1 1 0 Hệ số Cơ sở p.a A1 A2 A3 A4 A5 A6 2 A3 2 1 0 1 -1 0 0 0 A5 11 -1 2 0 1 1 0 0 A6 -8 -2 -1 0 3 0 1 -2 -3 0 -5 0 0 Nguyễn Hải Đăng - Khoa KHCB&KTCS - 37 -
Toán Kinh tế - Trường Cao Đẳng Công Nghiệp Nam Định 1 A3 -2 0 -1/2 1 1/2 0 1/2 0 A5 15 0 5/2 0 -1/2 1 -1/2 4 A1 4 1 1/2 0 -3/2 0 -1/2 0 -2 0 -10 0 -1 3 A2 4 0 1 -2 -1 0 -1 0 A5 5 0 0 5 2 1 2 4 A1 2 1 0 1 -1 0 0 20 0 0 -4 -8 0 -3 Tại đây ta có x³ 0 nên nó là phương án tối ưu, với x = (2, 4, 0, 0, 5, 0) và fMin = 20 Nguyễn Hải Đăng - Khoa KHCB&KTCS - 38 -
Toán Kinh tế - Trường Cao Đẳng Công Nghiệp Nam Định BÀI TẬP CHƯƠNG 2 I. Bài toán đối ngẫu. Tính chất của bài toán đối ngẫu 1. Viết bài toán đối ngẫu và chỉ ra cặp ràng buộc đối ngẫu của bài toán sau: a. f(X) = x1+3x2-x3+2x4 → min b. f(X) = 2x1-x2+5x3→ max ⎧−+xxx53+ =3 123 ⎧xx12−++32x3 x 4 ≤8 ⎪−+xxx22 − =− ⎪ ⎪ 234 ⎪−+−324xx12 x 3 ≤− ⎨ ⎨ 32xxx124++= 5 −+512xx ≤ ⎪ ⎪ 24 ⎪ xjj ≥=014(,) ⎪ ⎩ ⎩xjj ≥=014(,) c. f(x) = x1+2x2+3x3+4x4 + 5x5→ min d. f(x) = -4x1+2x2-5x3-x4+x5→ max ⎧ −7x1 + 3 x 2 − x 3 + 2 x 4 ≤ 1 ⎪ ⎧−+231xx13−+x4 x 5 ≥−2 ⎪ ⎪4x1+ 2 x 3 − x 4 + 3 x 5 ≥ − 2 xxxx+−+=323 ⎪ ⎪ 1235 ⎨ ⎨−2x2 + 3 x 3 + 5 x 4 − x 5 ≤ −4 43xx−+ xx − ≤ 4 ⎪ ⎪ 12 45 x1+3 x 2 − x 3 + 2 x 5 ≥ 3 ⎪ ⎪ ⎩xjj ≥=015(,) ⎪ ⎩xj ≥0( j = 1,5) e. f(x) = 2x1-3x2+4x3+5x4 -x5→ min ⎧ −x +3 x − 5 x + x − x = 4 ⎪ 1 2 3 5 6 ⎪2 x1− x 2 + x 4 +2 x 6 ≤ 8 ⎪ 3x1 + 7 x3 − 2 x 5 ≥ − 6 ⎪ ⎨5 x2− x 3 +3 x 4 = 12 ⎪−2 x + x +3 x ≥ 10 ⎪ 1 4 6 ⎪ x ≤ 1 ⎪ 6 ⎩⎪ x2, x 4 ≥ 0; x3 , x 5 ≤ 0; x1 , x 6 ∈ R 2. Dùng lý thuyết đối ngẫu chứng tỏ các bài toán sau giải được: a. f(X) = 2x1 +5x3 +3x4 → max ⎧xx134−+22x =−5 ⎪ ⎪34xxx23+−4 =− 9 ⎨ ⎪532xxx345++ = 14 ⎪ ⎩xxx123,,≤∈0 ; xx 25 , R b. f(X) = 3x1 +2x3 +4x5→ min ⎧23xx13−+x5 ≤16 ⎪ ⎪43xxx345−+ ≥ 9 ⎨ ⎪−+xx23 −211 x 5 =− ⎪xxx,,≥∈0 ; xx , R ⎩ 135 24 Nguyễn Hải Đăng - Khoa KHCB&KTCS - 39 -
Toán Kinh tế - Trường Cao Đẳng Công Nghiệp Nam Định 3. Dùng lý thuyết đối ngẫu, chứng minh rằng các bài toán sau không giải được: a. f(X) = -2x1 +x2+4x3+3x4 → min ⎧32xx12+−x3+530 x 4 = ⎪ ⎪2220xx23+− x 4 = ⎨22xxxx++− 2 = 12 ⎪ 1234 ⎪ ⎩xjj ≥=014(,) b. f(X) = 2x1 +3x3 -4x3+x4+2x5→ max ⎧−+25xx23−−xx45 3 ≥− 7 ⎪ xxx235+−211 ≤ ⎪ ⎨−+xxxxx1234546 − + − ≤ 0 ⎪32xx+= 24 ⎪ 34 ⎩⎪xx45,;,,≥∈0 xxx 123 R 4 . Cho bài toán: f(X) = -2x -6x +5x -x -4x →max 1 2 3 4 5 ⎧xx12−+−+42xxx34559 = 3 ⎪ ⎪xxxx2345−+−345 = 6 ⎨xxxx−+− =1 ⎪ 2345 ⎪ ⎩xjj ≥=015(,) Áp dụng lý thuyết đối ngẫu, chứng minh rằng X* = (0,0,16,31,14) là phương án tối ưu củ toán đã cho. HD: Viết bài toán đối ngẫu và sử dụng định lý độ lệch bù. 5 . Biết rằng X* = (0,5,0,3) là phương án tối ưu của bài toán: f(X) = 10x1+5x2+13x3+16x4 →min ⎧23xx12++xx34+ ≥16 ⎪ ⎪xxxx1234+++234 ≥ 22 ⎨34520xx++ x + x ≥ ⎪ 12 3 4 ⎪ ⎩xjj ≥=014(,) Tìm phương án tối ưu của bài toán đối ngẫu: 3 ĐS: Yopt = Y = (0, ,2). 2 6 . Cho bài toán: f(X) = x1+3x2+2x3+3x4+5x5 → min ⎧xx12+−+−=23xxx345 ⎪ ⎪−−xx23 +24 x 4 + x 5 ≥18 ⎨−+xxx32 + ≤ 10 ⎪ 235 ⎪ ⎩xjj ≥=015(,) 1 4 Tìm phương án tối ưu của bài toán đã cho biết rằng Y* = ( , , 0) là phương án tối ưu 3 3 của bài toán đối ngẫu. Nguyễn Hải Đăng - Khoa KHCB&KTCS - 40 -
Toán Kinh tế - Trường Cao Đẳng Công Nghiệp Nam Định II. Giải bằng phương pháp đơn hình đối ngẫu các bài toán quy hoạch tuyến tính sau: 1) f(x)=®x1+3x2+4x34+x Min ìx1-2x24+2x8³ ï ï-3x2+x34-4x£18 ĐS: x = (0, 0, 12, 4, 0, 22, 0) í-3x+x+2x-x³ 20 ï 1234 ï îxj³=0(j 1,4) 2) f(x)=®3x1+5x2-7x34+x Min ì2x2-x34+3x³10 ï ïx1+x2-4x34-x3= ĐS: Vô nghiệm í3x+2x+4x7£ ï 234 ï îxj³=0(j 1,4) 3) f(x)=-2x1+4x2-3x3-x545+®xMin ìx1+x2+3x3+x645-=x ï ï-2x2+x35-4x£15 ĐS: x = (0, 3, 0, 3, 0, 21, 0) í5x+2x+x-³x 18 ï2345 ï îxj³=0(j 1,5) 4) f(x)=®4x1+5x234+3x+2x Min ì2x1+x23+3x£ 21 ï ïx1+2x2+3x34-x³ 27 ĐS: x = (0, 6, 5, 0, 0, 0, 26) í-x+4x+2x+2x8³ ï 1234 ï îxj³=0(j 1,4) 5) f(x)=5x1+x2-3x3+x4+52x56+®xMin ì5x1-4x2+2x56-³3x 18 ï ï-9x1+x2-x56+2x=11 ĐS: x = (2, 33, 0, 0, 4, 0, 0) í6x-5x+x+3x-=4x 24 ï 13456 ï îxj³=0(j 1,6) Nguyễn Hải Đăng - Khoa KHCB&KTCS - 41 -
Toán Kinh tế - Trường Cao Đẳng Công Nghiệp Nam Định Chương 3. BÀI TOÁN VẬN TẢI 5.1. Nội dung và đặc điểm của bài toán 5.1.1. Nội dung bài toán a) Bài toán Có m địa điểm A1, A2, ,Am cùng sản xuất một loại hàng hóa với các lượng hàng tương ứng là a1, a2, ,am. Có n nơi tiêu thụ loại hang đó B1,B2, ,Bn với các yêu cầu tương ứng là b1, b2, ,bn. Để đơn giản ta sẽ gọi Ai: điểm phát i, i=1, m Bj: điểm thu j, j=1, n Hàng có thể chở từ một điểm phát bất kỳ (i) tới một điểm thu bát kỳ (j). Ký hiệu: cij- chi phí vận chuyển chở một đơn vị hàng từ điểm phát (i) đến điểm thu (j) xịj- lượng hàng chuyên chở từ i tới j Bài toán đặt ra là: xác định những đại lượng xij cho mọi con đường (i;j)sao cho tổng chi phí chuyên chở là nhỏ nhất với giả thiết: m n ∑ai = ∑ b j i=1 j=1 Tức là lượng hàng phát ra bằng đúng lượng hàng yêu cầu ở các điểm thu. b) Mô hình toán học Dạng toán học của bài toán vận tải mn ∑∑cxij ij → min ij==11 ⎧ n ⎪∑ xaimij== i ,1, ⎪ j=1 ⎪ m ⎪∑ xbjnij== j ,1, ⎨ i=1 ⎪ximjn≥=0, 1, ; = 1, ⎪ ij ⎪ mn ⎪abij,0;>=∑∑ a i b j ⎩ ij==11 Bài toán trên được gọi là bài toán cân bằng thu phát. Nguyễn Hải Đăng - Khoa KHCB&KTCS - 1 -
Toán Kinh tế - Trường Cao Đẳng Công Nghiệp Nam Định m n - Trường hợp ∑ ai > ∑ b j . Ta đưa về bài toán cân bằng thu phát bằng cách i=1 j=1 m n thêm vào một trạm thu giả Bn+1 với yêu cầu bn+1= ∑ ai − ∑ b j , đồng thời ci i=1 j=1 n+1 = 0 (∀i). Lượng hàng lấy từ trạm phát Ai cung cấp cho trạm thu giả B n+1, nghĩa là lượng hàng được giữ lại ở trạm phát Ai. m n - Trường hợp ∑ ai 0, còn gọi là ô sử dụng + ô loại: là ô không có hàng, tức là xij=0. + Dây chuyền: là một đoạn thẳng hay một dãy liên tiếp các đoạn thẳng gấp khúc mà hai đầu mút là hai ô chỉ nằm trên cùng một hàng hoặc một cột với một ô chọn khác thuộc dây chuyền của bảng vận tải. + Chu trình: là một dây chuyền khép kín Nguyễn Hải Đăng - Khoa KHCB&KTCS - 2 -
Toán Kinh tế - Trường Cao Đẳng Công Nghiệp Nam Định Như vậy một hàng hoặc một cột mà chu trình đi qua thì chỉ đi qua hai ô và do đó, số ô ít nhất của một chu trình là 4. + Ma trận X=(xij)m×n thỏa mãn hệ điều kiện ràng buộc được gọi là một phương án của bài toán. + Phương án X=(xij)m×n được gọi là phương án cực biên của bài toán vận tải nếu tập hợp các ô tương ứng với các thành phần dương của nó không tạo thành chu trình. + Phương án X=(xij)m×n được gọi là phương án cực biên không suy biến nếu nó có đúng m + n – 1 ô chọn. + Phương án X=(xij)m×n được gọi là phương án cực biên suy biến nếu nó có ít hơn m + n – 1 ô chọn. + Phương án X=(xij)m×n được gọi là phương án tối ưu (hay là nghiệm) của bài toán nếu nó thỏa mãn điều kiện (5.1), ký hiệu là X*. 5.1.2. Tính chất chung của bài toán - Một phương án cực biên có tối đa m + n – 1 thành phần dương. - Các véc tơ Aj tương ứng với biến xij có thành phần i và thành phần m+j bằng 1 còn các thành phần còn lại đều bằng 0. - Bài toán luôn luôn có lời giải. 5.2. Phương pháp thế vị để giải bài toán 5.2.1. Phương pháp tìm phương án cực biên xuất phát Để xây dựng môt phương án cực biên xuất phát người ta dùng một trong 3 phương pháp sau: * Phương pháp góc Tây_Bắc: Bước 1: Chọn ô ở dòng 1, cột 1 của bảng vận tải. Bước 2: Phân lượng hàng h = {a1,b1} vào ô (1;1). Bước 3: Đánh dấu hàng (cột), theo đó lượng hàng ở trạm phát (trạm thu) đã hết (đã đủ). Bước 4: Quay trở về bước 1 thực hiện công việc ở những ô còn lại. * Phương pháp chi phí nhỏ nhất:: Bước 1: Chọn ô có cước phí thấp nhất, giả sử là ô (i;j). Bước 2: Phân lượng hàng h = {ai,bj} vào ô (i;j). Bước 3: Đánh dấu các ô thuộc hàng i nếu trạm phát Ai đã hết hàng hoặc cột j nếu trạm thu Bj đã nhận đủ hàng. Bước 4: Quay trở lại bước 1 thực hiện công việc ở những ô còn lại. Nguyễn Hải Đăng - Khoa KHCB&KTCS - 3 -
Toán Kinh tế - Trường Cao Đẳng Công Nghiệp Nam Định * Phương pháp xấp xỉ Fogels: - Định nghĩa độ lệch của hàng(cột) là hiệu số giữa ô có cước phí thấp thứ nhì trừ đi ô có cước phí thấp thứ nhất ở hàng (cột) đó. Bước 1: Chọn hàng hoặc cột có độ lệch lớn nhất Bước 2: Chọn ô có cước phí thấp nhất thuộc hàng hoặc cột đó, giả sử là ô (i;j). Bước 3: Phân lượng hàng h = {ai,bj} vào ô (i;j). -Bước 4: Đánh dấu các ô thuộc hàng i nếu trạm phát Ai đã hết hàng hoặc cột j nếu trạm thu Bj đã nhận đủ hàng.Quay trở về bước 1 tiếp tục thực hiện thuật toán. 5.2.1. Phương pháp thế vị giải bài toán vận tải a) Tiêu chuẩn tối ưu Phương án cực biên không suy biến X=(xij)m×n được gọi là phương án tối ưu khi và chỉ khi tồn tại các số ui (i=1, ,m) cho các hàng và các số vj (j= 1, ,n)cho các cột của bảng vận tải sao cho: ⎪⎧ui+ v j = c ij ( i ; j ) : xij > 0(*) ⎨ ⎩⎪ui+ v j ≤ c ij ( i ; j ) : xij = 0( ) Phương trình (*) ứng với ô (i;j) là ô chọn. Phương trình ( ) ứng với ô (i;j) là ô loại. Các số ui và vj được gọi là hệ thống thế vị, trong đó ui được gọi là thế vị hàng, vj được gọi là thế vị cột. b) Thuật toán thế vị 0 Bước 1: Tìm phương án cực biên xuất phát X = =(xij)m×n Sử dụng một trong 3 phương pháp ở trên. Nếu phương án tìm được là phương án suy biến thì ta bổ sung ô chọn không để được phương án cực biên không suy biến, ô chọn này có vai trò như các ô chọn khác. Bước 2: Kiểm tra tính tối ưu của phương án. + Xây dựng hệ thống thế vị. Cho ui một giá trị tùy ý nào đó thì mọi giá trị khác đều xác định được một cách duy nhất do (*) có (n+m) ẩn và và (m+n-1) phương trình độc lập tuyến tính. + Tính các số kiểm tra Δij Đặt Δij =u i + v j − c ij ( i = 1, m , j = 1, n ) . Tính các Δ ij ứng với các ô loại. • Nếu Δij ≤0 (i = 1, m , j = 1, n ) thì phương án đang xét là phương án tối ưu. • Nếu ∃Δij >0 (i = 1, m , j = 1, n ) thì phương án đang xét chưa tối ưu, chuyển sang bước 3. Bước 3: Xây dựng phương án mới + Chọn ô điều chỉnh: Nguyễn Hải Đăng - Khoa KHCB&KTCS - 4 -
Toán Kinh tế - Trường Cao Đẳng Công Nghiệp Nam Định Δ =max Δ > 0 (i = 1, m , j = 1, n ) Ô (r,s) gọi là ô điều chỉnh nếu : rs { ij } + Tìm chu trình điều chỉnh: Là chu trình với ô xuất phát là ô điều chỉnh, các ô còn lại là ô chọn. Gọi V là tập hợp các ô thuộc chu trình điều chỉnh. + Đánh dấu các ô của chu trình, bắt đầu từ ô điều chỉnh đánh dấu (+) rồi xen kẽ nhau đánh dấu (-), (+) . cho đến hết chu trình. Gọi V+ là tập các ô có dấu (+), V- là tập các ô có dấu (-). Khi đó, VVV=+ ∪ − . + Xác định lượng hàng điều chỉnh: − q= min{ xij : ( i , j )∈ V} , q > 0. 1 + Điều chỉnh sang phương án mới: X= () xij m× n với: ⎧xij ,(,) i j∉ V ⎪ ⎪ + ⎨xij + q,(,) i j∈ V ⎪ − ⎩⎪ xij − q,(,) i j∈ V Gọi X1 đóng vai trò như X0 rồi quay lại bước 2 và lặp cho tới khi tìm được phương án tối ưu. * Chú ý: - Nếu ô điều chỉnh không duy nhất thì ta xét theo hàng từ trên xuống dưới trong bảng vận tải gặp ô nào trước thì ta chọn ô đó làm ô điều chỉnh. -Nếu có thể thì ta chọn ô (i0, j0) thỏa mãn: q.Δ =max{q . Δ} > 0 thì giá trị hàm mục tiêu giảm nhanh hơn. i0 j 0 ij c) Ví dụ Ví dụ 1. Giải bài toán vận tải với số liệu cho trong bảng sau: Thu 46 45 76 20 52 Phát 79 10 1 5 13 8 50 5 6 10 8 13 60 3 2 8 9 6 50 13 5 7 10 13 Giải: Bước 1: Tìm phương án cực biên xuất phát: Sử dụng phương pháp chi phí nhỏ nhất xây dựng phương án cực biên ta được phương án cực biên không suy biến cho trong bảng sau: Nguyễn Hải Đăng - Khoa KHCB&KTCS - 5 -
Toán Kinh tế - Trường Cao Đẳng Công Nghiệp Nam Định Thu 46 45 76 20 52 Phát 10 1 5 13 8 79 x 45 34 x x 5 6 10 8 13 50 x x x 20 30 3 2 8 9 6 60 46 x x x 14 13 5 7 10 13 50 x x 42 x 30 Bước 2: Kiểm tra tính tối ưu của phương án: + Xây dựng hệ thống thế vị: Cho u4 = 0, các ô (4,3) ,(4,5) là các ô cơ sở nên tính được v3 = 7 và v5 = 13. Xét cột 3 chẳng hạn ta thấy (1,3) là ô cơ sở nên tính được u1 = 7- 5 =2, Tiếp tục tính tương tự sẽ xác định được toàn bộ các thế vị hàng và cột như trong bảng dưới: + Tính các số kiểm tra: Tính các Δij (i = 1,4, j = 1,5) ta thấy các Δ15 = 3> 0,Δ 21 = 5> 0 nên phương án đang xét chưa tối ưu, ta chuyển sang bước 3. Thu u 46 45 76 20 52 i Phát 10 1 5 13 8 -2 79 -2 45 34 -7 3 5 6 10 8 13 50 0 5 -3 -3 20 30 3 2 8 9 6 60 -7 46 -6 -8 -8 14 13 5 7 10 13 50 0 -3 -2 42 -2 8 vj 10 3 7 8 13 Bước 3: Xây dựng phương án mới: + Chọn ô điều chỉnh: Nguyễn Hải Đăng - Khoa KHCB&KTCS - 6 -
Toán Kinh tế - Trường Cao Đẳng Công Nghiệp Nam Định Δ21 = max{} 5;3= 5 nên chọn ô (2,1) làm ô điều chỉnh.đánh dấu (+) vào ô điều chỉnh. + Chọn chu trình điều chỉnh: Chu trình điều chỉnh như trong bảng sau: Thu u 46 45 76 20 52 i Phát 10 1 5 13 8 79 -2 -2 45 34 -7 3 5 6 10 8 13 50 0 (+) 5 -3 -3 20 (-) 30 3 2 8 9 6 60 -7 (-) 46 -6 -8 -8 (+) 14 13 5 7 10 13 50 0 -3 -2 42 -2 8 vj 10 3 7 8 13 +Xác định lượng hàng điều chỉnh: q = min{} 46 ; 30= 30 + Điều chỉnh sang phương án mới: 1 1 X = (xij ) theo công thức cho trong bước 3 của thuật toán ta được bảng sau: Thu u 46 45 76 20 52 i Phát 10 1 5 13 8 79 -2 -2 45 34 -2 3 5 6 10 8 13 50 -5 (+) 30 -8 -8 (-) 20 -5 3 2 8 9 6 60 -7 (-) 16 -6 -8 -3 (+) 44 13 5 7 10 13 50 0 -3 -2 42 (+) 3 (-) 8 vj 10 3 7 13 13 + Lặp lại quá trình ta được chu trình như bảng trên. Ta được: q = min{} 16; 20;8= 8 + Điều chỉnh sang phương án mới X2 cho trong bảng dưới: Nguyễn Hải Đăng - Khoa KHCB&KTCS - 7 -
Toán Kinh tế - Trường Cao Đẳng Công Nghiệp Nam Định Thu u 46 45 76 20 52 i Phát 10 1 5 13 8 79 -2 -5 45 34 -5 0 5 6 10 8 13 50 -2 38 -5 -5 12 -5 3 2 8 9 6 60 -4 8 -3 -5 -3 52 13 5 7 10 13 50 0 -6 -2 42 8 -3 vj 7 3 7 10 10 Ta thấy mọi Δij ≤0 (i = 1,4, j = 1,5) nên phương án tương ứng ở bảng trên là tối ưu với giá trị hàm mục tiêu là f* = 45 * 1 + 34 * 5 + 38 * 5 + 12 * 8 + 8*3 + 52 * 6 + 42 * 7 + 8 * 10 = 1211. Nguyễn Hải Đăng - Khoa KHCB&KTCS - 8 -
Toán Kinh tế - Trường Cao Đẳng Công Nghiệp Nam Định Bμi tËp ch−¬ng 3 Giải các bài toán vận tải sau, sử dụng cả 3 phương pháp tìm phương án cực biên xuất phát: Bài 1 Thu 60 70 40 30 ⎡60 10 0 30⎤ Phát * ⎢ ⎥ * Đáp số: X =⎢ 0 60 20 0⎥ ;()f X =460 1 4 3 100 2 ⎣⎢ 0 0 20 0⎦⎥ 80 5 3 2 6 20 6 2 1 5 Bài 2 Thu Đáp số: 10 10 10 20 20 Phát ⎡0 5 0 0 0⎤ ⎢ ⎥ 1 4 6 10 0 0 0 5 5 5 7 X* =⎢ ⎥ ;()f X* =435 ⎢0 5 10 5 0⎥ 15 3 4 2 7 8 ⎢ ⎥ ⎣0 0 0 15 15⎦ 4 3 1 7 20 9 30 6 5 4 9 11 Bài 3 Thu 20 100 145 30 150 Phát Đáp số: ⎡ 0 0 120 0 0 ⎤ 120 6 3 1 4 5 ⎢ 0 0 0 0 150⎥ X* =⎢ ⎥ ;()f X* =940 150 1 2 5 4 3 ⎢20 75 25 30 0 ⎥ ⎢ ⎥ 150 2 4 3 1 6 ⎣ 0 25 0 0 0 ⎦ 25 3 1 4 2 7 Nguyễn Hải Đăng - Khoa KHCB&KTCS - 9 -