Chương1: BÀI TOÁN QUY HOẠCH TUYẾN TÍNH

Chia sẻ: Minh Nhí | Ngày: | Loại File: PDF | Số trang:13

Thêm vào BST

Báo xấu

293
lượt xem 37
download

Download Vui lòng tải xuống để xem tài liệu đầy đủ

Bài toán lập kế hoạch sản xuất Một công ty sản xuất n loại sản phẩm Sj (j=1,2, ..,n) sử dụng m loại nguyên liệu N i (i = 1,2, ..,m). Biết: + Lượng nguyên liệu Ni cần thiết dùng để sản xuất một đơn vị sản phẩm Sj là: aij + Trữ lượng nguyên liệu loại Ni là: bi + Tiền lãi một đơn vị sản phẩm Sj là: cj Hãy xây dựng kế hoạch sản xuất cho công ty để có lợi nhuận nhiều nhất.

Chủ đề:

Bình luận(0) Đăng nhập để gửi bình luận!

Lưu

Nội dung Text: Chương1: BÀI TOÁN QUY HOẠCH TUYẾN TÍNH

QUY HOẠCH TUYẾN TÍNH 02/09/2012 C hương1: BÀI TOÁN QUY HOẠCH TUYẾN TÍNH C hương1: BÀI TOÁN QUY HOẠCH TUYẾN TÍNH Bài toán lập kế hoạch sản xuất Một công ty sản xuất n loại sản phẩm Sj (j=1,2, ..,n) sử dụng m loại nguyên liệu N i (i = 1,2, ..,m). Biết: + Lượng nguyên liệu Ni cần thiết dùng để sản xuất MỘT SỐ BÀI TOÁN MỞ ĐẦU một đơn vị sản phẩm Sj là: aij + Trữ lượng nguyên liệu loại Ni là: bi + Tiền lãi một đơn vị sản phẩm Sj là: cj Hãy xây dựng kế hoạch sản xuất cho công ty để có lợi nhuận nhiều nhất. C hương1: BÀI TOÁN QUY HOẠCH TUYẾN TÍNH C hương1: BÀI TOÁN QUY HOẠCH TUYẾN TÍNH Minh hoạ dữ liệu bài toán: Giải: Gọi xj là lượng Sj (j=1,2, ..,n) cần sản xuất. xj ≥ 0. Sản phẩm Số nguyên S1 S2 … Sn Lượng nguyên liệu N1 dùng cho sản xuất: Nguyên liệu liệu tối đa a11x1 + a12x2 + … + a1nxn ≤ b1 N1 a11 a12 … a1n b1 Lượng nguyên liệu N2 dùng cho sản xuất: a21x1 + a22x2 + … + a2nxn ≤ b2 N2 a21 a22 … a2n b2 ………………………………………………………….. … … … … … … Lượng nguyên liệu Nm dùng cho sản xuất: Nm am1 am2 … amn bm am1x1 + am2x2 + … + amnxn ≤ bm Số tiền lãi mà công ty thu được là: Tiền lãi/đơn vị c1 c2 … cn f = f(x1, x2, … , xn) = c1x1 + c2x2 + … + cnxn sản phẩm C hương1: BÀI TOÁN QUY HOẠCH TUYẾN TÍNH C hương1: BÀI TOÁN QUY HOẠCH TUYẾN TÍNH Bài toán vốn đầu tư Mô hình toán học của bài toán Tìm x = (x1, x2, …, xn) sao cho Một xí nghiệp xử lý giấy có n phân xưởng Sj f(x) = c1x1 + c2x2 + … + cnxn ⟶ max (j=1,2,…,n) xử lý m loại giấy Ni (i=1,2,…,m). với các ràng buộc Biết: a11x1 + a12x2 + … + a1nxn ≤ b1 . + Lượng giấy loại Ni mà cuối năm phân xưởng Sj có a21x1 + a22x2 + … + a2nxn ≤ b2 .. thể xử lý (nếu cùng đầu tư một đơn vị vốn vào các ……........................................................... phân xưởng) là: aij am1x1 + am2x2 + … + amnxn ≤ bm + Lượng giấy Ni tối thiểu cần phải xử lý theo hợp x1, x2, …, xn ≥ 0 . đồng lao động là: bi Đây là một bài toán quy hoạch tuyến tính. Lập kế hoạch đầu tư để xí nghiệp hoàn thành hợp đồng lao động với tổng vốn đầu tư là nhỏ nhất. Chuongnn-hui.blogspot.com 1
QUY HOẠCH TUYẾN TÍNH 02/09/2012 C hương1: BÀI TOÁN QUY HOẠCH TUYẾN TÍNH C hương1: BÀI TOÁN QUY HOẠCH TUYẾN TÍNH Minh hoạ dữ liệu bài toán: Giải: Gọi xj là lượng vốn đầu tư cho Sj (j=1,2, ..,n). xj ≥ 0. Phân xưởng Lượng rác xử S1 S2 … Sn Lượng giấy N1 xử lý được: Loại rác thải lý tối thiểu a11x1 + a12x2 + … + a1nxn ≥ b1 N1 a11 a12 … a1n b1 Lượng giấy N2 xử lý được: a21x1 + a22x2 + … + a2nxn ≥ b2 N2 a21 a22 … a2n b2 ……………………………………… … … ……… … Lượng giấy Nm xử lý được: Nm am1 am2 … amn bm am1x1 + am2x2 + … + amnxn ≥ bm Số vốn mà xí nghiệp cần đầu tư là: f = f(x1, x2, … , xn) = x1 + x2 + … + xn C hương1: BÀI TOÁN QUY HOẠCH TUYẾN TÍNH C hương1: BÀI TOÁN QUY HOẠCH TUYẾN TÍNH Bài toán vận tải Mô hình toán học của bài toán Tìm x = (x1, x2, …, xn) sao cho Xí nghiệp cần vận chuyển hàng hoá từ m điểm phát f(x) = x1 + x2 + … + xn ⟶ min Pi (i=1,2,…,m) đến n điểm thu Tj (j=1,2,…,n). Biết với các ràng buộc + Lượng hàng có ở điểm phát Pi là: ai a11x1 + a12x2 + … + a1nxn ≥ b1 . + Lượng hàng cần ở mỗi điểm thu Tj là: bj a21x1 + a22x2 + … + a2nxn ≥ b2 . + Phí chuyển một đơn vị hàng từ Pi đến Tj là: cij …………………………………………….. Giả sử tổng lượng hàng có ở các điểm phát bằng am1x1 + am2x2 + … + amnxn ≥ bm tổng lượng hàng cần ở các điểm thu. x1, x2, …, xn ≥ 0 . Hãy lập kế hoạch vận chuyển hàng hoá với tổng chi Đây là một bài toán quy hoạch tuyến tính. phí là nhỏ nhất và đảm bảo yêu cầu thu phát. C hương1: BÀI TOÁN QUY HOẠCH TUYẾN TÍNH C hương1: BÀI TOÁN QUY HOẠCH TUYẾN TÍNH Minh hoạ dữ liệu bài toán: Giải: Gọi xij là lượng hàng cần chuyển từ Pi (i=1,2, ..,m) Điểm thu T1:b1 T2:b2 … Tn:bn đến Tj (i=1,2, ..,m). xij ≥ 0. Điểm phát Lượng hàng chuyển đi từ Pi là: c11 c12 c1n P1:a1 … xi1 + xi2 + … + xin = ai …………………………………………… c21 c22 c2n Lượng hàng chuyển đến Tj là: P2:a2 … x1j + x2j + … + xmj = bj … … … … … …………………………………………… Tổng chi phí vận chuyển là: cm1 cm2 cmn Pm:am … f = f(x11, x12, …, xmn) = c11x11 + c12x12 + … + cmnxmn Chuongnn-hui.blogspot.com 2
QUY HOẠCH TUYẾN TÍNH 02/09/2012 C hương1: BÀI TOÁN QUY HOẠCH TUYẾN TÍNH C hương1: BÀI TOÁN QUY HOẠCH TUYẾN TÍNH Bài tập: Mô hình toán học của bài toán Tìm x = (x11, x12, …, xmn) sao cho Lập mô hình toán học của các bài toán sau : f(x) = c11x11 + c12x12 + … + cmnxmn ⟶ min 1. Một xí nghiệp sản xuất ba loại sản phẩm A, B và C chế tạo từ ba loại nguyên liệu I, II và III . Lượng với các ràng buộc nguyên liệu I, II và III mà xí nghiệp có lần lượt là x11 + x12 + … + x1n = a1 20, 40, 30. Lượng nguyên liệu I, II, III cần cho một ……………………………………. đơn vị sản phẩm loại A lần lượt là 3, 3, 2, loại B lần xm1 + xm2 + … + xmn= am lượt là 3, 2, 1, loại C lần lượt là 2, 1, 1 đơn vị. x11 + x21 + … + xm1 = b1 Hãy lập kế hoạch sản xuất để xí nghiệp thu tiền lãi ……………………………………. nhiều nhất, biết tiền lãi một đơn vị sản phẩm loại A x1n + x2n + … + xmn = bn lãi 2 triệu đồng, loại B lãi 4 triệu đồng và loại C lãi x11, x12, …, xmn ≥ 0 . 3 triệu đồng. C hương1: BÀI TOÁN QUY HOẠCH TUYẾN TÍNH C hương1: BÀI TOÁN QUY HOẠCH TUYẾN TÍNH 2. Một Xí nghiệp chăn nuôi cần mua thức ăn tổng 3. Một Xí nghiệp xử lý giấy có ba phân xưởng I, II, hợp T1, T2, T3 cho gia súc với tỷ lệ chất dinh III cùng xử lý ba loại giấy A, B, C. Nếu đầu tư 10 dưỡng là: 1 kg T1 chứa 3 đơn vị D1, 3 đơn vị D2 và triệu đồng vào mỗi phân xưởng thì cuối năm phân 2 đơn vị D3; 1 kg T2 chứa 3 đơn vị D1, 2 đơn vị D2 xưởng I xử lý được 3 tạ giấy A, 3 tạ giấy B, 2 tạ giấy và 1 đơn vị D3; 1 kg T3 chứa 2 đơn vị D1, 1 đơn vị C, phân xưởng II xử lý được 3 tạ giấy A, 2 tạ giấy B, D2 và 1 đơn vị D3. Mỗi bữa ăn, gia súc cần tối thiểu 1 tạ giấy C, phân xưởng III xử lý được 2 tạ giấy A, 1 20 đơn vị D1, 40 đơn vị D2 và 30 đơn vị D3. Biết tạ giấy B, 1 tạ giấy C. rằng 1 kg T1 có giá là 2 ngàn đồng, 1 kg T2 có giá là Theo yêu cầu Xí nghiệp phải xử lý ít nhất 2 tấn giấy 4 ngàn đồng, 1 kg T3 có giá là 3 ngàn đồng. loại A, 4 tấn giấy loại B, 3 tấn giấy loại C. Hỏi cần Hỏi Xí nghiệp phải mua bao nhiêu kg T1, T2, T3 đầu tư vào mỗi phân xưởng bao nhiêu tiền để xí cho một bữa ăn để bảo đảm tốt về chất dinh dưỡng nghiệp hoàn thành công việc với giá tiền đầu tư là và tổng số tiền mua là nhỏ nhất? nhỏ nhất. C hương1: BÀI TOÁN QUY HOẠCH TUYẾN TÍNH C hương1: BÀI TOÁN QUY HOẠCH TUYẾN TÍNH 4. Một xí nghiệp chế biến đồ gỗ dùng 3000 đơn vị 5. Một xí nghiệp có thể sử dụng tối đa 510 giờ máy gỗ nguyên liệu nhóm I, 5000 đơn vị gỗ nguyên liệu cán, 360 giờ máy tiện và 150 giờ máy mài để chế nhóm II và 2000 đơn vị gỗ nguyên liệu nhóm III để tạo 3 sản phẩm A, B và C. Để chế tạo một sản phẩm sản xuất tủ trang trí, bàn ghế và giường cao cấp. A cần 9 giờ máy cán, 5 giờ máy tiện và 3 giờ máy Một bộ tủ trang trí bán giá 500000đ dùng 30 đơn mài; một sản phẩm B cần 3 giờ máy cán, 4 giờ máy vị nhóm I, 10 đơn vị nhóm II và 10 đơn vị nhóm III. tiện; một sản phẩm C cần 5 giờ máy cán, 3 giờ máy Một bộ bàn ghế bán giá 800000đ dùng 40 đơn vị tiện và 2 giờ máy mài. Mỗi sản phẩm A trị giá 48 nhóm I, 20 đơn vị nhóm II và 50 đơn vị nhóm III. ngàn đồng, mỗi sản phẩm B trị giá 16 ngàn đồng và Một bộ giường bán giá 400000đ dùng 10 đơn vị nhóm I, 50 đơn vị nhóm II và 80 đơn vị nhóm III. mỗi sản phẩm C trị giá 27 ngàn đồng. Hãy cho biết Hãy xác định số lượng các sản phẩm cần sản xuất xí nghiệp cần chế tạo mỗi loại bao nhiêu sản phẩm sao cho xí nghiệp đạt lợi nhuận nhiều nhất? để có tổng giá trị sản phẩm lớn nhất . Chuongnn-hui.blogspot.com 3
QUY HOẠCH TUYẾN TÍNH 02/09/2012 C hương1: BÀI TOÁN QUY HOẠCH TUYẾN TÍNH C hương1: BÀI TOÁN QUY HOẠCH TUYẾN TÍNH 6. Một người nông dân muốn sản xuất lúa gạo và 7. Một xưởng làm cửa sắt có những thanh thép dài lúa mì trên 50 ha đất của ông ta với nguồn nước 12 m, cần cắt thành 8 đoạn dài 4m, 5 đoạn dài 5m dự trữ là 90 ngàn m3 nước và nguồn nhân lực gổm và 3 đoạn dài 7m. Có 5 mẫu cắt sau. Mẫu 1: 3 đoạn 250 công. Biết rằng để sản xuất 1 tấn lúa gạo cần dài 4m, không thừa. Mẫu 2: 1 đoạn 4m, 1 đoạn 5m, sử dụng 3 ha đất, 5 ngàn m3 nước và 15 công với thừa 3m. Mẫu 3: 1 đoạn 4m, 1 đoạn 7m, thừa 1m. lợi nhuận là 18 USD, để sản xuất một tấn lúa mì Mẫu 4: 2 đoạn 5m, thừa 2m. Mẫu 5: 1 đoạn 5m, 1 cần sử dụng 3 ha đất, 4 ngàn m3 nước và 12 công đoạn 7m, không thừa. với lợi nhuận là 21 USD. Hỏi cần dùng những mẫu cắt nào để tiết kiệm nhất. Hỏi người nông dân phải sản xuất bao nhiêu tấn lúa mỗi loại để đạt lợi nhuận cao nhất . C hương1: BÀI TOÁN QUY HOẠCH TUYẾN TÍNH C hương1: BÀI TOÁN QUY HOẠCH TUYẾN TÍNH 8. Một xưởng sản xuất dự định mua hai loại máy để in hình vẽ trên vải. Máy A có thể in 100m/phút và chiếm 50 mét vuông diện tích sàn, còn máy B có thể in 200m/phút và chiếm 140 mét vuông diện tích sàn. Xưởng cần in ít nhất 600m/phút và có BÀI TOÁN QUY HOẠCH TUYẾN TÍNH diện tích sàn để đặt máy in tối đa là 350 mét vuông. Mỗi máy A có giá là 22 triệu đồng và mỗi máy B có giá 42 là triệu đồng. Hỏi cần mua bao nhiêu máy in mỗi loại sao cho tốn ít chi phí nhất . C hương1: BÀI TOÁN QUY HOẠCH TUYẾN TÍNH C hương1: BÀI TOÁN QUY HOẠCH TUYẾN TÍNH Bài toán quy hoạch tuyến tính tổng quát V í dụ: + f(x) = x1 – 2x3 – x4 – 2x6 ⟶ max Tìm x = (x1, x2, … , xn) sao cho x1 + 2 x 2 – x4 + x5 ≤ 8 f(x) = c1x1 + c2x2 + … + cnxn ⟶ max(min) (1) x1 – 3 x 3 – x4 + 3 x 6 ≥ 2 với các ràng buộc 2 x1 + x2 + x3 – x4 = 5 ai1x1 + ai2x2 + … + ainxn ≤ bi, i ∊ I1 (2) – x1 + x3 + 2 x 4 – 3 x 5 ≥ 6 ai1x1 + ai2x2 + … + ainxn ≥ bi, i ∊ I2 (3) 2 x2 – 2 x3 – x4 – 2 x6 ≤ 1 1 ai1x1 + ai2x2 + … + ainxn = bi, i ∊ I3 (4) x1, x4 ≥ 0 xj ≥ 0, j ∊ J1 (5) Với I1∪I2∪I3={1,2,…,m}, . x2, x3, x6 ≤ 0 xj ≤ 0, j ∊ J2 (6) . I1∩I2∩I3=∅ x5 ∊ ℝ xj ∊ ℝ, j ∊ J3 (7) và J1∪J2∪J3={1,2,…,n},. J1∩J2∩J3=∅ Chuongnn-hui.blogspot.com 4
QUY HOẠCH TUYẾN TÍNH 02/09/2012 C hương1: BÀI TOÁN QUY HOẠCH TUYẾN TÍNH C hương1: BÀI TOÁN QUY HOẠCH TUYẾN TÍNH Một số khái niệm Dạng chính tắc của bài toán quy hoạch tuyến tính Tìm x = (x1, x2, … , xn) sao cho Hàm mục tiêu: hàm f(x) = c1x1 + c2x2 + … + cnxn f(x) = c1x1 + c2x2 + … + cnxn ⟶ max(min) Phương án : vectơ x = (x1, x2, … , xn) thoả mãn các ai1x1 + ai2x2 + … + ainxn = bi, i = 1, 2, …, m điều kiện ràng buộc từ (2) đến (7). xj ≥ 0, j = 1, 2, …, n . Phương án tối ưu: phương án làm cho hàm mục hay tiêu đạt max hay min nghĩa là thoả mãn (1). f(x) = ∑cjxj ⟶ max(min) Tập phương án tối ưu : tập hợp tất cả các phương ∑aijxj = bi, i = 1, 2, …, m án của bài toán. xj ≥ 0, j = 1, 2, …, n . Giải bài toán quy hoạch tuyến tính : là tìm phương Lưu ý: Mọi bài toán quy hoạch tuyến tính đều có án tối ưu cho bài toán. thể biến đổi đưa về dạng chính tắc. C hương1: BÀI TOÁN QUY HOẠCH TUYẾN TÍNH C hương1: BÀI TOÁN QUY HOẠCH TUYẾN TÍNH V í dụ: Biến đổi về dạng chính tắc Để biến đổi bài toán về dạng chính tắc tìm min ta + f(x) = x1 – 3x2 + 2x3 – x4 ⟶ max thực hiện như sau: 2 x1 + x 2 + x3 – x4 ≤ 8 + Biến hàm mục tiêu f ⟶ max thành g ⟶ min, với x1 – x3 + x4 = 5 g =  f. x1 – 2 x 2 – x3 + 3 x4 ≥ 7 + Biến ràng buộc ∑aijxj ≤ bi thành ∑aijxj + xn+i = bi x1, x3 ≥ 0 bằng cách bổ sung thêm biến xn+i ≥ 0. x2 ≤ 0 g(x) = – x – 2x – 3x + x – x ⟶ min + Biến ràng buộc ∑aijxj ≥ bi thành ∑aijxj – xn+i = bi 1 3 7 8 9 x4 ∊ ℝ bằng cách bổ sung thêm xn+i ≥ 0. 2 x1 + x3 + x5 – x 7 – x8 + x9 = 8 + Đổi biến xj ≤ 0 thành biến x′j ≥ 0 với x′j =  xj x1 – x3 + x8 – x9 = 5 + Đổi biến xj  ℝ thành hiệu x′j  x″j với x′j , x″j ≥ 0 x1 – x3– x6 + 2x7 + 3x8 – 3x9= 7 và x′j  x″j = xj xj ≥ 0, j=1,2,…,9 C hương1: BÀI TOÁN QUY HOẠCH TUYẾN TÍNH C hương1: BÀI TOÁN QUY HOẠCH TUYẾN TÍNH Phương án cực biên + f(x) = 2x1 – x2 + x3 – 3x4 ⟶ min x1 + 2 x 2 + 2 x3 – x4 ≥ 6 Xét bài toán 2 x 1 – x2 – 2 x3 + x 4 = 9 f(x) = ∑cjxj ⟶ min 3x1 –x2 + x3 + 2x4 ≤ 12 ∑aijxj = bi, i = 1, 2, …, m x2, x4 ≥ 0 xj ≥ 0, j = 1, 2, …, n . x1 ≤ 0 f(x) = – x2 – 3x4 – 2x7 + x8 – x9 ⟶ min Vectơ liên kết với biến xj là vectơ cột Aj = [aij]m1 2x2– x4 – x5 – x7 + 2x8 – 2x9 = 6 x3 ∊ ℝ có các thành phần là hệ số tương ứng của biến xj. – x2 + x4 – 2 x7 – 2 x 8 + 2 x9 = 9 Phương án cực biên: phương án mà hệ vectơ liên –x2 + 2x4 + x6 – 3x7 + x8 – x9 = 12 kết với các biến xj > 0 tạo thành một hệ vectơ độc xj ≥ 0, j=1,2,…,9 lập tuyến tính. Lưu ý: Từ đây ta chỉ xét các bài toán dạng chính tắc tìm min. Chuongnn-hui.blogspot.com 5
QUY HOẠCH TUYẾN TÍNH 02/09/2012 C hương1: BÀI TOÁN QUY HOẠCH TUYẾN TÍNH C hương1: BÀI TOÁN QUY HOẠCH TUYẾN TÍNH Biến cơ sở là biến xj > 0, vectơ Aj tương ứng với V í d ụ: biến cơ sở gọi là vectơ cơ sở, biến phi cơ sở là biến + Xét bài toán xj = 0 . f(x) = 4x1 + x2 + x3 ⟶ min x1 +2x2  x3 = 5 Phương án cực biên không suy biến là phương án x1  x 2 + 2 x3 = 5 có đúng m biến cơ sở, nếu số biến cơ sở bé hơn m ta có phương án cực biên suy biến. xj ≥ 0, j=1,2,3 Số các phương án cực biên của một bài toán quy Vectơ nào sau là phương án cực biên không suy hoạch tuyến tính là hữu hạn. biến: x0 = (1, 4, 4), x1 = (5, 0, 0), x2 = (0, 5, 5)? ĐS: x2 = (0, 5, 5) C hương1: BÀI TOÁN QUY HOẠCH TUYẾN TÍNH C hương1: BÀI TOÁN QUY HOẠCH TUYẾN TÍNH + Xét bài toán Lưu ý: Cách thành lập một phương án cực biên của f(x) = x1 + 2x2  x3 ⟶ min bài toán quy hoạch tuyến tính dạng chính tắc: x 1 + x2 + x3 = 4 + Xác định hệ vectơ liên kết với n biến. x 1  x2 = 0 + Tìm các hệ con {Ai}, i = 1, 2, …,m gồm m véctơ xj ≥ 0, j=1,2,3 độc lập tuyến tính của hệ vectơ Số hệ con này hữu Vectơ nào sau là phương án cực biên không suy hạn. biến: x0 = (1, 1, 2), x1 = (2, 2, 0), x2 = (0, 0, 4)? + Biểu diễn vectơ b theo hệ con {Ai} ở trên, ta ĐS: x1 = (2, 2, 0) được các hệ số biểu diễn. Thành lập vectơ x có các thành phần là hệ số biểu diễn. + Loại đi những vectơ x có thành phần âm, các véctơ còn lại là các phương án cực biên. C hương1: BÀI TOÁN QUY HOẠCH TUYẾN TÍNH C hương1: BÀI TOÁN QUY HOẠCH TUYẾN TÍNH V í dụ: Tìm các phương án cực biên của bài toán Lưu ý: + f(x) = 2x1 + x3 + 5x4 ⟶ min Do số thành phần dương của một phương án cực biên không suy biến là bằng m nên suy ra số thành x 1 + x3 + x4 = 5 phần bằng 0 sẽ là n – m. 9 1 ĐS: (5, 1, 0, 0), ( , 0, 0, ), x 2  x3 + 2 x4 = 1 2 2 Từ đó suy ra muốn tìm phương án cực biên ta có (0, 6, 5, 0), (0, 0, 3, 2) xj ≥ 0, j=1,2,3,4 thể cho n – m thành phần bằng 0 rồi tính giá trị của + f(x) = 2x1  x3 + 2x4 ⟶ min m thành phần còn lại bằng cách giải hệ m phương x 1 + x2 + x3 + x 4 = 1 0 trình m ẩn. 16 14 2 x 2 + x3  x4 = 6 ĐS: (0, 0, 8, 2), (0, , 0, ), 3 3 xj ≥ 0, j=1,2,3,4 (4, 0, 6, 0), (7, 3, 0, 0) Chuongnn-hui.blogspot.com 6
QUY HOẠCH TUYẾN TÍNH 02/09/2012 C hương1: BÀI TOÁN QUY HOẠCH TUYẾN TÍNH C hương1: BÀI TOÁN QUY HOẠCH TUYẾN TÍNH Các tính chất chung của bài toán quy hoạch V í dụ: Tìm các phương án cực biên của bài toán tuyến tính + f(x) = x1 + 6x3  5x4 ⟶ min + Tập các phương án của bài toán quy hoạch tuyến x 1 + 2 x3 + 3 x4 = 5 8 75 ĐS: (5, , 0, 0), (0, , , 0), tính là một tập lồi nghĩa là nếu x, y là hai phương 3 x 2  x3 + 2 x4 = 8 3 22 án bất kỳ của bài toán thì mọi tổ hợp x + (1 )y, 14 5 (0, , 0, ) xj ≥ 0, j=1,2,3,4  ∊ ℝ: 0 ≤  ≤ 1 cũng là một phương án của bài 9 3 + f(x) = x1  6x3 + x4  x5 ⟶ max toán. x 1 + 2 x4 + x5 = 8 + Tập các phương án tối ưu của bài toán quy hoạch tuyến tính cũng là một tập lồi. x 2 + x4  x5 = 4 ĐS: (8,4,6,0,0), (2,10,0,0,6), + Nếu bài toán quy hoạch tuyến tính dạng chính (0,0,2,4,0), (0,6,0,2,4) x 3 + x4 + x5 = 6 tắc có tập phương án khác rỗng thì nó có ít nhất xj ≥ 0, j=1,2,…,5 một phương án cực biên. C hương1: BÀI TOÁN QUY HOẠCH TUYẾN TÍNH C hương1: BÀI TOÁN QUY HOẠCH TUYẾN TÍNH V í dụ: Giải các bài toán quy hoạch tuyến tính + Điều kiện cần và đủ để bài toán quy hoạch tuyến + f(x) = x1 + 6x3  5x4 ⟶ min tính dạng chính tắc có phương án tối ưu là nó có tập phương án khác rỗng và hàm mục tiêu bị chặn. x1 + 2 x3 + 3 x4 = 5 14 5 ĐS: x* = (0, , 0, ), + Nếu bài toán quy hoạch tuyến tính dạng chính 3 x2  x3 + 2 x4 = 8 9 3 tắc có phương án tối ưu thì nó có ít nhất một 25 fmin =  xj ≥ 0, j=1,2,3,4 phương án cực biên là phương án tối ưu. 3 + f(x) = x1  6x3 + x4  x5 ⟶ max Lưu ý: Từ đây, ta có thể giải bài toán quy hoạch x1 + 2 x4 + x5 = 8 tuyến tính dạng chính tắc nếu biết nó có phương x2 + x 4  x5 = 4 án tối ưu bằng cách tìm tất cả các phương án cực ĐS: x* = (0,6,0,2,4), biên của bài toán, phương án tối ưu là phương án fmax = 2 x3 + x 4 + x5 = 6 mà giá trị hàm mục tiêu lớn nhất (hay nhỏ nhất). xj ≥ 0, j=1,2,…,5 C hương1: BÀI TOÁN QUY HOẠCH TUYẾN TÍNH C hương1: BÀI TOÁN QUY HOẠCH TUYẾN TÍNH Giải bài toán quy hoạch tuyến tính hai + f(x) = x + 3y ⟶ max biến bằng phương pháp hình học x + 2y ≤ 6 2x + y ≤ 10 V í dụ: ‒x + y ≤ 4 + f(x) = 2x +3y ⟶ min x ≥ 0, y ≥ 0 3x + y ≥ 3 ĐS: fmax = 9 tại A (0;3) x  4y ≤ 6 x + 2y ≤ 6 x ≥ 0, y ≥ 0 ĐS: fmin = 2 tại C(1; 0) Chuongnn-hui.blogspot.com 7
QUY HOẠCH TUYẾN TÍNH 02/09/2012 C hương1: BÀI TOÁN QUY HOẠCH TUYẾN TÍNH C hương1: BÀI TOÁN QUY HOẠCH TUYẾN TÍNH Cơ sở của phương pháp đơn hình Dựa vào phương án hiện có, ta tìm cách đánh giá phương án đó có là phương án tối ưu hay chưa, nếu phương án đang xét đã là phương án tối ưu thì mục đích của ta đã đạt, được nếu chưa là phương PHƯƠNG PHÁP ĐƠN HÌNH án tối ưu thì ta thay thế bởi nó một phương án mới GIẢI BÀI TOÁN QUY HOẠCH TUYẾN TÍNH tốt hơn. Xét bài toán dạng chính tắc (P) f(x) = ∑cjxj ⟶ min ∑aijxj = bi, i = 1, 2, …, m xj ≥ 0, j = 1, 2, …, n C hương1: BÀI TOÁN QUY HOẠCH TUYẾN TÍNH C hương1: BÀI TOÁN QUY HOẠCH TUYẾN TÍNH 00 x0 0 Định lý: ( Dấu hiệu tối ưu) Giả sử ta có phương án = x1 , x 2 , … , x m , 0 , … , 0 với các biến cơ sở tương ứng x1, x2, …, xm và hệ 00 0 Nếu x = x1 , x2 , … , xm , 0, … , 0 là một phương án vectơ cơ sở liên kết là {Ai}, i = 1,2,…,m. cực biên của (P) sao cho ∆j ≤ 0 với mọi j = 1,2,…,n 0 0 Giá trị mục tiêu f x 0 = c1 x1 + c2 x2 + ⋯ + cm xm. 0 thì x là phương án tối ưu. Định lý: ( Dấu hiệu không có phương án tối ưu) Đặt ∆j = c1 x1j + c2 x2j + ⋯ + cm xmj − cj Nếu ngoài cơ sở liên kết của phương án cực biên với xj = (x1j, x2j, …, xmj) là toạ độ của các vectơ Aj tồn tại j sao cho ∆j > 0 và xj ≤ 0 nghĩa là xij≤ 0 với đối với hệ cơ sở {Ai}, gọi là ước lượng của biến xj. mọi i = 1,2,…,m thì (P) không có phương án tối ưu. Lưu ý: Nếu ma trận cột tạo bởi hệ vectơ cơ sở {Ai} Rõ hơn là hàm mục tiêu không bị chặn dưới trên là ma trận đơn vị thì xij = aij. tập phương án. C hương1: BÀI TOÁN QUY HOẠCH TUYẾN TÍNH C hương1: BÀI TOÁN QUY HOẠCH TUYẾN TÍNH V í dụ: + Xét bài toán + Xét bài toán f(x) = 7x1 ‒ 26x2 + 9x3 ⟶ min f(x) = x1 + 6x2 + 9x3 ⟶ min x1 ‒ 2 x 2 = 5 x 1 + 2 x3 = 6 ‒x2 + x3 = 7 x 2 + x3 = 8 xj ≥ 0, j = 1, 2, 3. xj ≥ 0, j = 1, 2, 3. Vectơ x = (5, 0, 7) có phải là phương án tối ưu hay không? Vectơ x = (6, 8, 0) có phải là phương án tối ưu hay không? ĐS: x không phải là phương án tối ưu. ĐS: x là phương án tối ưu. fmin = 54. Chuongnn-hui.blogspot.com 8
QUY HOẠCH TUYẾN TÍNH 02/09/2012 C hương1: BÀI TOÁN QUY HOẠCH TUYẾN TÍNH C hương1: BÀI TOÁN QUY HOẠCH TUYẾN TÍNH Thuật toán đơn hình Định lý: Xét bài toán Nếu ngoài cơ sở liên kết của phương án cực biên tồn tại j sao cho ∆j > 0 và có ít nhất một xij > 0 với i f(x) = ∑cjxj ⟶ min nào đó thì ta luôn có thể tìm được một phương án ∑aijxj = bi, i = 1, 2, …, m cực biên mới tốt hơn x, nghĩa là phương án này xj ≥ 0, j = 1, 2, …, n làm cho hàm mục tiêu nhỏ hơn phương án x. bi ≥ 0, i=1,2,…,m và có sẵn ma trận đơn vị. + Tìm phương án cực biên ban đầu: Giả sử ma trận đơn vị tạo bởi các vectơ cột Ai, i = 1,2,…,m thì ta có phương án cực biên ban đầu 00 x 0 = x1 , x 2 , … , x m , 0 , … , 0 0 C hương1: BÀI TOÁN QUY HOẠCH TUYẾN TÍNH C hương1: BÀI TOÁN QUY HOẠCH TUYẾN TÍNH + Đánh giá phương án hiện có: lập bảng đơn hình + Xây dựng phương án mới – Phương pháp phần như hình bên dưới và căn cứ vào đó để đánh giá tử trục quay: Xác định cột quay s là cột mà s > 0 lớn nhất, biến Biến Hệ Phương x1 x2 … xn  cơ sở tương ứng xs là biến cơ sở đưa vào. cơ sở số án c1 c2 … cn x0 0 x1 c1 x11 x12 … x1n x1 Dòng quay r là dòng mà r nhỏ nhất với θr = r , xks 0 x2 c2 x21 x22 … x2n x2 biến cơ sở tương ứng xr là biến cơ sở đưa ra. … … … … … … … Phần tử quay xrs là phần tử giao giữa dòng quay r 0 xm cm xm xm1 xm2 … xmn và cột quay s. 1 2 n f(x0) … C hương1: BÀI TOÁN QUY HOẠCH TUYẾN TÍNH C hương1: BÀI TOÁN QUY HOẠCH TUYẾN TÍNH Xác định các thành phần phương án x1 mới: V í dụ: Giải các bài toán + f(x) = x1 ‒ x2 ‒ 2x4 + 2x5 ‒ 3x6 ⟶ min PA s↓ x0 r nếu i = r x1 + x 4 + x5 ‒ x6 = 2 i xrs 0 x1 xi xis = i x0 xrs − x0 xis x2 + x 4 + x6 = 1 2 r nếu i ≠ r i r 0 xr xrs xrs x3 + 2 x4 + 4 x5 + 3 x6 = 9 Xác định các hệ số biểu diễn xij mới xj ≥ 0, j = 1, 2, …,6. ĐS: Phương án tối ưu x* = (0, 8, 0, 3, 0, 1). Giá trị j↓ s↓ xrj nếu i = r tối ưu fmin = ‒17. i xrs ′ xij xis xij = xij xrs −xrj xis nếu i ≠ r r xrs xrj xrs Chuongnn-hui.blogspot.com 9
QUY HOẠCH TUYẾN TÍNH 02/09/2012 C hương1: BÀI TOÁN QUY HOẠCH TUYẾN TÍNH C hương1: BÀI TOÁN QUY HOẠCH TUYẾN TÍNH Lưu ý: + f(x) = ‒2x1 ‒ 4x2 + x3 ‒ x4 ⟶ min + Nếu bài toán có bi ≥ 0, i=1,2,…,m và chưa có sẵn x 1 + 3 x2 + x5 = 4 ma trận đơn vị ta có thể thực hiện biến đổi ma trận 2 x 1 + x2 ‒ x3 + x6 = 3 hệ số mở rộng của hệ điều kiện ràng buộc để có ma x 2 + 4 x3 + x4 = 3 trận đơn vị. xj ≥ 0, j = 1, 2, …,6. V í dụ: Giải bài toán + f(x) = ‒ 3x1 + x2 + 3x3 ‒ x4 ⟶ min ĐS: Phương án tối ưu x* = (1, 1, 0, 2, 0, 0). Giá trị tối ưu fmin = ‒8. x1 + 2 x2 ‒ x3 + x4 = 2 2x1‒ 6x2 + 3x3 + 3x4 = 9 x1 ‒ x 2 + x3 ‒ x4 = 6 xj ≥ 0, j = 1,2,3,4. C hương1: BÀI TOÁN QUY HOẠCH TUYẾN TÍNH C hương1: BÀI TOÁN QUY HOẠCH TUYẾN TÍNH ĐS: Bài toán trở thành + f(x) = 3x1 ‒ 4x2 ‒ 5x3 + 6x4 ⟶ min f(x) = ‒ 3x1 + x2 + 3x3 ‒ x4 ⟶ min x1 + x 2 + x3 + 1 3 x4 = 1 4 6 2x1+ x2 + 14x4 = 11 x 1 + x4 = 3 5 3x2 + x3 + 14x4 = 16 12 x2 ‒ x= 2 xj ≥ 0, j = 1,2,3,4. 54 23 x3 ‒ x= 5 54 xj ≥ 0, j = 1, 2, 3, 4. C hương1: BÀI TOÁN QUY HOẠCH TUYẾN TÍNH C hương1: BÀI TOÁN QUY HOẠCH TUYẾN TÍNH Thuật toán đơn hình mở rộng ĐS: Bài toán trở thành Xét bài toán f(x) = 3x1 ‒ 4x2 ‒ 5x3 + 6x4 ⟶ min f(x) = ∑cjxj ⟶ min 27 x1 + x= 4 54 ∑aijxj = bi, i = 1, 2, …, m 16 xj ≥ 0, j = 1, 2, …, n x2 + x= 3 54 với bi ≥ 0, i=1,2,…,m và chưa có sẵn ma trận đơn vị 22 x3 + x= 7 (giả sử bài toán còn thiếu m vectơ đơn vị). 54 Ta bổ sung thêm m biến xn+i, i=1,2,…,m để có được xj ≥ 0, j = 1, 2, 3, 4. ma trận đơn vị và đưa về xét hàm mục tiêu g = ∑cjxj + Mxn+1 + Mxn+2 + … + Mxn+m với M là một số rất lớn, lớn hơn bất cứ số nào khác. Chuongnn-hui.blogspot.com 10
QUY HOẠCH TUYẾN TÍNH 02/09/2012 C hương1: BÀI TOÁN QUY HOẠCH TUYẾN TÍNH C hương1: BÀI TOÁN QUY HOẠCH TUYẾN TÍNH Bài toán trở thành bài toán M sau đây Ta dùng thuật toán đơn hình để giải bài toán M với g(x) = ∑cjxj + Mxn+1 + Mxn+2 + … + Mxn+m ⟶ min một vài lưu ý: ∑aijxj + xn+i = bi, i = 1, 2, …, m + Nếu đã có k vectơ đơn vị thì ta bổ sung vào m – k xj ≥ 0, j = 1, 2, …, n+m biến để có ma trận đơn vị cấp m. Lưu ý: + Hàm mục tiêu g, ước lượng j của các biến xj có Nếu bài toán M có phương án tối ưu (x, t) với t > 0 thể viết dưới dạng g = AM + B, j = jM+ j. t = (xn+1, …, xn+m) thì bài toán chính không có + Để so sánh các ước lượng ta có phương án. Nếu bài toán M có phương án tối ưu αj > αk αj > 0 (x, 0) thì x là phương án tối ưu của bài toán chính ∆j ≥ 0 ⇔ , ∆j ≥ ∆k ⇔ và nếu bài toán M không có phương án tối ưu thì αj = αk , βj ≥ βk αj = 0, βj ≥ 0 bài toán chính cũng không có phương án tối ưu . và các hệ số A, B, ,  được trình bày trên 2 dòng. C hương1: BÀI TOÁN QUY HOẠCH TUYẾN TÍNH C hương1: BÀI TOÁN QUY HOẠCH TUYẾN TÍNH V í d ụ: Giải bài toán + f(x) = 3x1 ‒ 4x2 ‒ 5x3 + 6x4 ⟶ min + f(x) = ‒ 3x1 + x2 + 3x3 ‒ x4 ⟶ min x1 + x 2 + x3 + 1 3 x4 = 1 4 x 1 + 2 x2 ‒ x3 + x4 = 2 2x1+ x2 + 14x4 = 11 2x1‒ 6x2 + 3x3 + 3x4 = 9 3x2 + x3 + 14x4 = 16 x 1 ‒ x2 + x3 ‒ x 4 = 6 xj ≥ 0, j = 1,2,3,4. xj ≥ 0, j = 1,2,3,4. ĐS: Phương án tối ưu x* = (4, 3, 7, 0). Giá trị tối ưu fmin = ‒35. ĐS: Phương án tối ưu x* = (3, 2, 5, 0). Giá trị tối ưu fmin = 8. C hương1: BÀI TOÁN QUY HOẠCH TUYẾN TÍNH C hương1: BÀI TOÁN QUY HOẠCH TUYẾN TÍNH Lưu ý: Với bài toán dạng chính tắc tìm max ta thực 2 17 39 ĐS: Phương án tối ưu x* = ( , , , 0, 0). Giá trị hiện tương tự, chỉ khác là nếu ∆j ≥ 0, ∀j = 1,2,…,n 36 2 88 thì phương án đang xét là tối ưu, nếu có ∆j < 0 mà tối ưu fmax = . 3 xj ≤ 0 thì bài toán không có lời giải, cột quay s là + f(x) = 3x1 ‒ x2 ‒ 2x3 ⟶ max cột tương ứng với ∆s < 0 nhỏ nhất. ‒x1 + 3x2 + x3 + x4 = 7 V í dụ: Giải bài toán 3x1‒ 4x2 + 8x3 + x5 = 10 + f(x) = 2x1 + 3x2 + x3 ⟶ max x 1 ‒ 5 x2 + x3 = 6 4x1‒ 2x2 + x6 = 12 2x1+ 2x2 + x4 = 7 xj ≥ 0, j = 1,2, …,6. ‒x1 + 2x2 + x5 = 5 ĐS: Phương án tối ưu x* = (5, 4, 0, 0, 11, 0). Giá trị xj ≥ 0, j = 1,2,3,4,5. tối ưu fmax = 11. Chuongnn-hui.blogspot.com 11
QUY HOẠCH TUYẾN TÍNH 02/09/2012 C hương1: BÀI TOÁN QUY HOẠCH TUYẾN TÍNH C hương1: BÀI TOÁN QUY HOẠCH TUYẾN TÍNH Giải bài toán quy hoạch tuyến tính với EXCEL trong đó ở dòng phương án, ta gán các giá trị ban đầu là 1 (hay 0) cho các biến xj ≥ 0. Cài thêm công cụ Add-ins Solver + Options ⟶ Add-Ins ⟶ Excel Add-Ins. + Trong hộp thoại Add-Ins ta click chọn Solver Add-In. Tổ chức dữ liệu trên bảng tính Tiến hành nhập liệu trên bảng tính Excel như sau C hương1: BÀI TOÁN QUY HOẠCH TUYẾN TÍNH C hương1: BÀI TOÁN QUY HOẠCH TUYẾN TÍNH Tiến hành giải bài toán + Data ⟶ Solver + Trong hộp thoại Solver Parameters: Set Objective: click chọn ô chứa giá trị mục tiêu B8. Max/Min/Equal To: chọn loại bài toán để giải. By Changing Cells: click chọn địa chỉ tuyệt đối của các ô ghi các giá trị ban đầu của biến. Subject to the Constraints: nhập các ràng buộc. C hương1: BÀI TOÁN QUY HOẠCH TUYẾN TÍNH C hương1: BÀI TOÁN QUY HOẠCH TUYẾN TÍNH Ta có thể sửa đổi, thêm bớt các ràng buộc. Muốn thêm vào ràng buộc ta click chọn nút Add. trong hộp thoại Add Constraint: Cell Reference: click chọn địa chỉ chứa công thức. Ô dấu: lựa chọn dấu của các ràng buộc tương ứng. Constraint: Ô chứa giá trị vế phải của các ràng buộc tương ứng (ta cũng có thể nhập trực tiếp giá trị vế phải của ràng buộc tương ứng). Chuongnn-hui.blogspot.com 12
QUY HOẠCH TUYẾN TÍNH 02/09/2012 C hương1: BÀI TOÁN QUY HOẠCH TUYẾN TÍNH C hương1: BÀI TOÁN QUY HOẠCH TUYẾN TÍNH Select a solving method: click chọn phương pháp Hãy lập kế hoạch sản xuất sao cho lãi thu được giải, ở đây là Simplex LP phương pháp đơn hình. nhiều nhất? Trong hộp thoại Solver Result: Minh hoạ dữ liệu bài toán: Keep Solver Solution: Lấy kết quả, in ra bảng tính. Sản phẩm Số nguyên A B Restore Original Values: Huỷ kết quả vừa tìm được Nguyên liệu liệu tối đa và trả các biến về tình trạng ban đầu. I 2 1 8 Save Scenario: Lưu kết quả vừa tìm được thành II 0 6 24 một tình huống để có thể xem lại sau này. III 4 0 12 Lưu ý: Có 3 loại báo cáo là Answer, Sensitivity và Tiền lãi/đơn vị Limits. 3 5 sản phẩm C hương1: BÀI TOÁN QUY HOẠCH TUYẾN TÍNH C hương1: BÀI TOÁN QUY HOẠCH TUYẾN TÍNH Tương tự hãy lập mô hình toán học cho bài toán : Mô hình toán học của bài toán Một công ty sản xuất 2 loại sơn nội thất và sơn Tìm x = (x1, x2) sao cho ngoài trời, sử dụng 2 loại nguyên liệu A, B với trữ f(x) = 3x1 + 5x2 ⟶ max lượng là 6 tấn và 8 tấn. Một tấn sơn nội thất cần 2 tấn nguyên liệu A và 1 tấn nguyên liệu B, một tấn với các ràng buộc sơn ngoài trời cần 1 tấn nguyên liệu A và 2 tấn 2 x 1 + x2 ≤ 8 nguyên liệu B. Lượng sơn nội thất không hơn sơn 6 x2 ≤ 2 4 ngoài trời quá 1 tấn, lượng sơn nội thất tối đa là 2 tấn. Giá bán một tấn sơn nội thất là 2000USD, giá 4 x1 ≤ 1 2 bán một tấn sơn ngoài trời là 3000USD. x1, x2 ≥ 0 Hỏi cần sản xuất mỗi loại sơn bao nhiêu tấn để có doanh thu lớn nhất ? C hương1: BÀI TOÁN QUY HOẠCH TUYẾN TÍNH Bài toán tổng quát: Một công ty sản xuất n loại sản phẩm Sj (j=1,2, ..,n) sử dụng m loại nguyên liệu N i (i = 1,2, ..,m). Biết: + Lượng nguyên liệu Ni cần thiết dùng để sản xuất một đơn vị sản phẩm Sj là: aij + Trữ lượng nguyên liệu Ni là: bi + Tiền lãi một đơn vị sản phẩm Sj là: cj Hãy xây dựng kế hoạch sản xuất cho công ty để có lợi nhuận nhiều nhất. Chuongnn-hui.blogspot.com 13