Đang tải... (xem toàn văn)
27316Nghiên cứu một số bài toàn mở rộng: vận tải ba chỉ số không hạn chế và có hạn chế khả năng thông qua, bài toán vận tải ba chỉ số khoảng và giới thiệu một số bài toán vận tải đa mục tiêu
Website: http://www.docs.vn Email : lienhe@docs.vn Tel (: 0918.775.368 lêi më đầu Cùng với phát triển mạnh mẽ khoa học kỹ thuật, toán tối u xuất ngày nhiều tính phức tạp chúng ngày lớn Phạm vi khả ứng dụng toán tối u ngày đa dạng phong phú Lớp toán tối u quan trọng đợc nghiên cứu đợc ứng dụng nhiều toán quy hoạch tuyến tính (linear programming) Đó mô hình toán học lớp rộng lớn toán ứng dụng kinh tế kỹ thuật Do cấu trúc lớp toán quy ho¹ch tun tÝnh cã nhiỊu tÝnh chÊt rÊt tèt mặt toán học, ngời ta đà tìm đợc thuật giải hữu hiệu cho toán Năm 1947 nhà toán học Mỹ G.B Dantzig đà nghiên cứu đề xuất thuật toán đơn hình (simplex method) để giải toán quy hoạch tuyến tính Thuật toán đơn hình đợc phát triển mạnh mẽ năm sau đợc xem phơng pháp kinh điển để giải toán quy hoạch tuyến tính Đây phơng pháp đợc sử dụng nãi lµ réng r·i nhÊt Cã ba lý chÝnh: Một là: Rất nhiều vấn đề thực tế, nhiều lĩnh vực khác đa toán quy hoạch tuyến tính Hai là: Trong nhiều phơng pháp giải toán phi tuyến, toán tuyến tính xuất nh toán phụ cần phải giải nhiều bớc lặp Ba là: Phơng pháp đơn hình phơng pháp hiệu để giải toán quy hoạch tuyến tính Ngày nay, thuật toán đơn hình dạng cải biên chúng, ngời ta giải nhanh toán QHTT cỡ lớn Lớp toán vận tải trờng hợp đặc biệt quy hoạch tuyến tính, dùng phơng pháp quy hoạch tuyến tính để giải Tuy nhiên, tính chất đặc thù riêng nó, ngời ta xây dựng phơng pháp giải riêng Thông thờng nói đến toán vận tải ta thờng liên hệ đến toán vận tải hai số, toán vận tải kinh điển có phơng pháp giải hay Bên cạnh đó, ngời ta xét số toán vận tải mở rộng nh toán vận tải ba số, toán vận tải khoảng, toán vận tải đa mục tiêu nhiều toán khác, biến thể toán vận tải kinh điển Trong khuôn khổ khoá luận này, em xem xét nghiên cứu số toán mở rộng lớp toán vận tải mở rộng Đó toán: Bài toán vận tải ba số (solid transport problem) không hạn chế có hạn chế khả thông qua, Bài toán vận tải ba số khoảng (interval solid transport problem) giới thiệu số Bài toán vận tải đa mục tiªu Trang: Website: http://www.docs.vn Email : lienhe@docs.vn Tel (: 0918.775.368 Cuối cùng, em xin bày tỏ lòng biết ơn sâu sắc thày giáo hớng dẫn Thạc sỹ Vũ Tiến Việt, ngời đà tận tình bảo, giúp đỡ em trình hoàn thành khoá luận Em xin chân thành cảm ơn giúp đỡ nhiệt tình thầy cô khoa Toán - Tin, Häc viƯn An ninh Nh©n d©n Trang: Website: http://www.docs.vn Email : lienhe@docs.vn Tel (: 0918.775.368 Ch¬ng I Bài toán quy hoạch tuyến tính Trong việc nghiên cứu toán tối u nói chung, giải tích lồi giữ vai trò quan trọng Nó đợc sử dụng làm sở toán học việc xây dựng thuật toán Quy hoạch tuyến tính lớp toán tối u đợc nghiên cứu trọng vẹn phơng diện lý thuyết lẫn thực hành, Bài toán vận tải dạng đặc biệt QHTT Do chơng nhằm giới thiệu số khái niệm kiến thức giải tích lồi QHTT 1.1 Một số khái niệm giải tÝch låi 1.1.1 Kh«ng gian Euclude Mét vector n chiỊu trờng số thực đợc thứ tù gåm n sè thùc x=(x1, x2, , xn) C¸c xi, i =1, , n gọi thành phần hay toạ độ vector Ví dụ x=(4,5,10,20) Hai vectơ x vµ y gäi lµ b»ng x=y, nÕu xi=yi, i =1, , n Xét hai phép toán vector: PhÐp céng: x+y=(x1+y1, x2+y2, , xn+yn) PhÐp nh©n: x=(x1, x2, , xn), R Khi tập hợp tất vector n chiều xác định phép cộng vector, nhân số thực với vector nh tạo thành không gian tuyến tính n chiều trờng số thực R, ký hiệu Rn Các vector x(i) Rn, i =1, , m đợc gọi độc lËp tuyÕn tÝnh nÕu: m x i (i ) i 0, i 1, m i 1 m NÕu: x i x (i ) víi Ýt nhÊt mét i th× x gäi tổ hợp tuyến tính i m x , i =1, , m Hơn i 0, i =1, , m vµ (i) i x gọi tổ hợp i lồi cđa c¸c x(i), i =1, , m Trong Rn cã n vector độc lập tuyến tính lập thành sở Giả sử e(1), e(2), , e(n) sở Rn vector x Rn tổ hợp tuyến tính vector e(1), e(2), , e(n) Ta gäi tÝch v« híng cđa hai vector x=(x1, x2, , xn) vµ y=(y1, y2, , yn), ký hiƯu, , lµ mét sè b»ng Trang: Website: http://www.docs.vn Email : lienhe@docs.vn Tel (: 0918.775.368 n x, y xi y i i 1 Tích vô hớng dạng song tuyến tính, đối xứng, không âm, tức là: = x,y Rn =< x(1), y >+< x(2), y> x(1), x(2), y Rn = x,y Rn 0, x Rn dÊu b»ng xÈy vµ chØ x= Độ dài vector x=(x1, x2, , xn) số xác định x x, x n x i i 1 Kho¶ng cách hai vector x y số xác định bởi: x, y x y x y , x y n x i yi i Không gian vector có tích vô hớng khoảng cách nh gọi không gian Euclude 1.1.2 TËp compact D·y {x(k) }Rn k=1, 2, đợc gọi có giới hạn x(0) k vµ viÕt lim x(k) = x(0), nÕu k lim x ( k ) , x ( ) k Hình cầu tâm a bán kính tập S=xRn :x-a Hình cầu tạo nên - lân cận điểm a, hay gọi lân cận a * Nếu tập ARn chứa với điểm x lân cận x gọi điểm A Nếu lân cận x A có điểm A điểm không thuộc A x gọi điểm biên tập hợp A * Một tập ARn gọi giới nội đợc chứa hình cầu tâm O đó, tức tồn số đủ lớn cho với mäi xA,x Mét d·y {x(k)} héi tơ th× giới nội * Một tập hợp GRn đợc gọi mở với xG tồn hình cầu tâm x nằm gọn G Một tập FRn đợc gọi đóng với dÃy héi x (k ) F tơ{x(k)} F ta ®Ịu có: klim Một tập chứa điểm biên tập đóng Trang: Website: http://www.docs.vn Email : lienhe@docs.vn Tel (: 0918.775.368 * Tập C đợc gọi tập Compact từ dÃy vô hạn {x (k)} thc C ®Ịu cã thĨ trÝch mét d·y {x(ki)} héi tơ tíi phÇn tư thc C TËp C Compact C đóng giới néi TËp Compact M cđa tËp ®ãng C cịng ®ãng C TËp ®ãng M cđa tËp Compact cịng Compact Hàm f(x) liên tục tập Compact C đạt cực trị tập 1.1.3 Tập lồi Cho hai điểm a, b Rn Ta gọi đờng thẳng qua a, b tập điểm có dạng xRn : x = a + (1-)b, R Đoạn thẳng nối hai điểm a, b tập lồi điểm cã d¹ng xRn :x = x + (1-)y, * Một tập MRn đợc gọi đa tạp affine với hai điểm x, y M toàn đờng thẳng qua hai điểm thuộc M Tức x + (1-)y M : x,y M, R * Mét siªu phẳng không gian Rn tập hợp tất điểm x=(x1, x2, , xn) Rn thỏa mÃn phơng tr×nh tuyÕn tÝnh a1x1+ a2x2+ + anxn = ®ã a1, a2, , an , R * Tập hợp điểm x=(x1, x2, , xn) Rn thoản mÃn bất phơng trình tuyến tính a1x1+ a2x2+ + anxn đợc gọi nửa không gian đóng * Nửa không gian đợc cho a1x1+ a2x2+ + anxn < đợc gọi nửa không gian mở * Tập XRn đợc gọi tập lồi cïng víi viƯc chøa hai ®iĨm x, y nã chøa đoạn thẳng chứa hai điểm ấy, tức chứa tất điểm có dạng: x + (1-)y, VÝ dơ vỊ c¸c tËp låi: Không gian Euclide, nửa không gian, mặt phẳng, nửa mặt phẳng, hình chữ nhật, hình vuông, hình elip, hình hộp, hình cầu * Một tập hợp giao số hữu hạn nửa không gian đóng đợc gọi tập lồi đa diện Mệnh đề: Giao hai tập lồi tập lồi Hệ Giao cđa mét sè bÊt kú tËp hỵp låi tập lồi Hệ Miền chứa nghiệm hệ bất phơng trình tuyến tính dạng tËp låi (®a diƯn låi) Mét tËp låi ®a diƯn giới nội gọi đa diện a1 x1 a x1 a x n1 a1 x a 22 x a n x a1n x n a n x n a nn x n b b2 bn Giao tất tËp låi chøa tËp X gäi lµ bao låi cđa nã, ký hiƯu [X] 1.1.4 Hµm låi * Mét hµm số f(x) xác định tập lồi C Rn đợc gọi hàm lồi C, với x, y C vµ ta cã f(x + (1-)y) f(x) + (1-)f(y) Trang: Website: http://www.docs.vn Email : lienhe@docs.vn Tel (: 0918.775.368 * Hàm f(x) đợc gọi hàm lồi chặt víi mäi x, y C vµ ta cã f(x + (1-)y) f(x) + (1-)f(y) * Hàm f(x) đợc gọi hàm lõm (lõm chặt) - f(x) hàm lồi (lồi chặt) * Hàm f(x) xác định C đạt cực tiểu tuyệt đối x* C f(x*) f(x): xC * Hàm f(x) đạt cực tiểu địa phơng x* C tồn lân cận mở U x* cho f(x*) f(x): xC U MƯnh ®Ị 1: BÊt kỳ điểm cực tiểu địa phơng hàm lồi tập lồi điểm cực tiểu tuyệt đối Hệ quả: Bất kỳ điểm cực đại địa phơng hàm lõm cực đại tuyệt đối Mệnh đề 2: Cực đại hàm lồi (nếu có) tập lồi có điểm cực biên đạt điểm cực biên 1.2 Bài toán Quy hoạch tuyến tính QHTT bắt nguồn từ nghiên cứu nhà toán học Nga tiếng, Viện sỹ L.V Kantorovich loạt công trình toán kế hoạch hoá sản xuất, công bố năm 1938 Năm 1947 nhà toán học Mỹ G.B Dantzig đà nghiên cứu đề xuất phơng pháp đơn hình (Simplex method) để giải toán QHTT Năm 1952 phơng pháp đơn hình đà đợc chạy máy tính điện tử Mỹ 1.2.1 Bài toán quy hoạch tuyến tính Bài toán tổng quát Để quán lập luận ta xét toán tìm cực đại, sau ta xét cách chuyển toán tìm cực tiểu sang tìm cực đại Bài toán tổng quát QHTT có dạng: Ký hiệu: A=(aij)mxn ma trận với phần tử aij n 1.1 c x max j j j 1 n a x , , b , i 1, , m ij j 1.2 i j 1 1.3 (1.1) gọi hàm mục tiêu, (1.2) buộc Nếu gặp toán Min, tức x j 0, j 1, , n n f x c j x j j 1 D Thì giữ nguyên ràng buộcxvà đa toán Max b»ng c¸ch n f x c j x j max j 1 xD Trang: Website: http://www.docs.vn Email : lienhe@docs.vn Tel (: 0918.775.368 NÕu toán Max có phơng án tối u x* toán có phơng án x* fmin=-fmax Thật vậy, x* phơng án tối u toán Max nên ta có: n n f max c j x *j c j x j , j 1 x D hay j 1 n n c j 1 * j x j c j x j , x D j Chứng tỏ x* phơng án tối u toán Min n f c j x *j f max j 1 Dạng chuẩn dạng tắc Ngời ta thờng xét toán quy hoạch tuyến tính dới hai dạng sau: -D¹ng chuÈn: n c j x j max j 1 n a x j bi , ij i 1, , m j 1 x j 0, j 1, , n -Dạng tắc: n c j x j max j 1 n a ij x j bi , i 1, m j 1 x j 0, j 1, , n Đa toán QHTT dạng chuẩn dạng tắc Bất kỳ QHTT đa hai dạng chuẩn tắc nhờ phép biến đổi tuyến tính sau: i) Mét rµng buéc n a ij x j bi j Có thể đa ràng buộc lại n j 1 n j 1 a ij x j bi cách nhân hai vế với (-1) viÕt a ' ij x j b' i ii) Một ràng buộc đẳng thức n a ij x j bi j 1 cã thĨ thay b»ng hai rµng bc bất đẳng thức: n aij x j bi , j 1 n a ij x j bi j 1 Trang: Website: http://www.docs.vn Email : lienhe@docs.vn Tel (: 0918.775.368 iii) Một biến xj không bị ràng buộc dÊu cã thĨ thay thÕ bëi hiƯu cđa hai biÕn không âm cách đặt: iv) Một ràng buộc bất ®¼ng thøc n a ij x j bi j 1 x j x j x j , víi x j 0, x j Có thể đa ràng buộc đẳng thức cách đa vào biến phụ yi 0: n a ij x j y i bi j Về nguyên tắc, áp dụng nhiều lần phép biến đổi (i), (ii) (iii) ta đa toán QHTT dạng chuẩn, sau áp dụng nhiều lần phép biến đổi (iv) ta đa dạng tắc Giải toán QHTT phơng pháp hình học Xét toán QHTT dới dạng chuẩn với hai biến số: c1 x1 c x max a i1 x1 a i1 x bi , D j 1,2 x j 0, i 1, , m Từ ý nghĩa hình học ta biết bất phơng trình tuyến tính ai1x1+ai2x2 bi xác định nửa mặt phẳng Nh miền ràng buộc D đợc xác định nh giao nửa mặt phẳng đa giác lồi mặt phẳng Phơng trình c1x1+c2x2= thay đổi xác định mặt phẳng đờng thẳng song song với mà ta gọi đờng mức (với giá trị mức ) Mỗi điểm D nằm mét ®êng møc víi x x1 , x2 mức Bài toán đặt thể phát biểu theo ngôn ngữ hình học nh sau: số đ ccó x1 c2 x2 êng møc c¾t tËp D, hÃy tìm đờng mức với gía trị lớn Nếu dịch chuyển song song đờng mức theo hớng vector ph¸p tun cđa chóng n c , c giá trị mức tăng, dịch chuyển theo hớng ngợc lại giá trị mức giảm Vì để giải toán đặt ra, ta tiến hành nh sau Bắt đầu từ đờng mức cắt D, ta dịch chuyển song song đờng mức theo hớng vector pháp tuyến (c1,c2) việc dịch chuyển làm cho đờng mức không cắt D dừng Điểm D (có thể nhiều điểm) nằm đờng mức cuối lời giải tối u cần tìm, giá trị hàm mục tiêu giá trị tối u toán Ví dụ: Xét to¸n: f(x)= 4x1+5x2max 8 x1 x 7 x1 x x2 3 x1 0, Trang: x2 0 Website: http://www.docs.vn Email : lienhe@docs.vn Tel (: 0918.775.368 XÐt ®êng mức: 4x1+5x2=10 Đờng mức qua hai điểm (0,2) vµ (2.5,0) Ta cã x*=(3,2), fmax=22 y vµ x* lµ đỉnh D Qua phơng pháp hình học ta thÊy r»ng: n - NÕu quy ho¹ch tuyÕn tÝnh cã phơng án tối u có đỉnh tối u Sở dĩ nói có trờng hợp đờng mức vị trí giới hạn trùng với cạnh * D tất điểm cạnh xphơng án tối u, ®ã cã hai ®Ønh - NÕu miỊn rµng bc D giới nội khác rỗng chắcx chắn có phơng án tối u - Nếu miền ràng buộc không giới nội nhng hàm mục tiêu bị chặn trên miền ràng buộc chắn có phơng án tèi u 1.2.2 Mét sè tÝnh chÊt chung MƯnh ®Ị 1: Tập hợp tất phơng án toán QHTT tập lồi Tập lồi D phơng án toán QHTT xác định toàn ràng buộc (1.2) (1.3) Tập D rỗng, đa diện lồi tập lồi đa diện không giới nội Nếu D đa diện lồi toán có phơng án, giá trị tối u hàm mục tiêu đa diện lồi hữu hạn việc tìm phơng án tối u đa đến việc chọn điểm đa diện D có số đỉnh (điểm cực biên hay phơng án cực biên) hữu hạn Mệnh đề 2: Hàm mục tiêu toán QHTT đạt Max điểm cực biên tập D Nếu hàm mục tiêu không nhận Max điểm cực biên tập lồi D mà nhiều điểm cực biên đạt giá trị cực đại điểm tổ hợp lồi điểm Ký hiệu Aj, j=1, , n vector cét cđa ma trËn A Khi Êy hƯ rµng buéc Ax =b cã thÓ viÕt: x1A1 + x2A2 + + xnAn = b (1.4) MƯnh ®Ị 3: NÕu vector A1, A2, , Ak độc lập tuyến tính thoả mÃn x1A1+x2A2+ +xnAn=b xj 0, j=1, k điểm x=(x1,x2, ,xk,0, ,0) Trang: Website: http://www.docs.vn Email : lienhe@docs.vn Tel (: 0918.775.368 điểm cực biên tập lồi đa diện D Mệnh đề 4: Nếu x =(x1,x2, ,xn) điểm cực biên tập lồi đa diện D vector Aj biểu diễn (1.4) ứng với thành phần x j lập thành hệ độc lập tuyến tính Vì ma trận A có m dòng nên từ suy điểm cực biên m thành phần dơng Các mệnh đề mệnh đề gộp lại thành mệnh đề sau: Mệnh đề 5: Để x =(x1,x2 ,xn) phơng án cực biên QHTT dới dạng tắc cần đủ vector cột Aj ma trận A ứng với thành phần x j độc lập tuyến tính 1.2.3 Phơng pháp đơn hình giải toán QHTT Cơ sở phơng pháp đơc G.B Dantzig công bố năm 1947 có tên gọi phơng pháp đơn hình Sở dĩ có tên gọi nh toán đợc giải phơng pháp có ràng buộc dạng: n 1.5ràng x j 1, x j 0,điểm j 1, , nxRn thoả mÃn Mà tậphợp buộc đơn hình j không gian n chiều Đờng lối chung sở thuật toán i) Đờng lối chung Phơng pháp đơn hình dựa hai nhận xét sau: toán QHTT có phơng án tối u, đa diện lồi D có số hữu hạn đỉnh Nh phải tồn thuật toán hữu hạn Thuật toán gồm hai giai đoạn: - Giai đoạn 1: Tìm phơng án cực biên (một ®Ønh) - Giai ®o¹n 2: KiĨm tra ®iỊu kiƯn tèi u đối phơng án tìm đợc giai đoạn Nếu điểu kiện tối u đợc thoả mÃn phơng án tối u Nếu không, ta chuyển sang phơng án cực biên cho cải tiến giá trị hàm mục tiêu Tiếp theo lại kiểm tra điều kiện tối u phơng án Ngời ta thùc hiƯn mét d·y c¸c thđ tơc nh vËy nhận đợc phơng án tối u, đến tình toán phơng án tối u ii) Cơ sở thuật toán Xét toán QHTT dới dạng tắc: c, x max Ax = b (1.6) (1.7) x0 (1.8) Trong ®ã A ma trận kích thớc mxn giả sử hạng ma trận A m (điều không làm tính tổng quát) x vector Trang: 10 ... Chơng Bài toán vận tải toán vận tải mở rộng 2.1 Bài toán vận tải hai số 2.1.1 Phát biểu toán tính chất Có m địa điểm A1, A2, , An sản xuất loại hàng hóa với lợng hàng tơng ứng a1, a2, , an Có n... tập lồi đa diện không giới nội Nếu D đa diện lồi toán có phơng án, giá trị tối u hàm mục tiêu đa diện lồi hữu hạn việc tìm phơng án tối u đa đến việc chọn điểm đa diện D có số đỉnh (điểm cực biên... cuối (không kể f ) có số âm xem thử có cột cắt dòng cuối số âm mà số cột âm hay không Nếu có cột nh toán phơng án tối u Nếu trái lại chọn cột cho 0 s k k Råi chän ( số dòng cắt cột s số
