Đang tải... (xem toàn văn)
Xem xét & nghiên cứu Một số bài toán mở rộng trong lớp các bài toán vận tải mở rộng (Toán - Tin)
lời mở đầuCùng với sự phát triển mạnh mẽ của khoa học kỹ thuật, các bài toán tối u xuất hiện ngày càng nhiều và tính phức tạp của chúng ngày càng lớn. Phạm vi và khả năng ứng dụng của các bài toán tối u cũng ngày càng đa dạng và phong phú.Lớp bài toán tối u quan trọng đợc nghiên cứu đầu tiên và đợc ứng dụng nhiều nhất là bài toán quy hoạch tuyến tính (linear programming). Đó là mô hình toán học của một lớp rộng lớn các bài toán ứng dụng trong kinh tế và kỹ thuật. Do đó cấu trúc của lớp bài toán quy hoạch tuyến tính có nhiều tính chất rất tốt về mặt toán học, ngời ta đã tìm đợc các thuật giải rất hữu hiệu cho bài toán này. Năm 1947 nhà toán học Mỹ G.B. Dantzig đã nghiên cứu và đề xuất ra thuật toán đơn hình (simplex method) để giải bà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ẽ trong những năm sau đó và đợc xem là một phơng pháp kinh điển để giải các bài toán quy hoạch tuyến tính. Đây là một phơng pháp đợc sử dụng có thể nói là rộng rãi nhất. Có ba lý do chính:Một là: Rất nhiều vấn đề thực tế, trong nhiều lĩnh vực khác nhau có thể đa về bài toán quy hoạch tuyến tính.Hai là: Trong nhiều phơng pháp giải các bài toán phi tuyến, bài toán tuyến tính xuất hiện nh là một bài toán phụ cần phải giải trong nhiều bớc lặp.Ba là: Phơng pháp đơn hình là phơng pháp hiệu quả để giải bài toán quy hoạch tuyến tính. Ngày nay, bằng thuật toán đơn hình và các dạng cải biên của chúng, ngời ta có thể giải rất nhanh các bài toán QHTT cỡ lớn.Lớp các bài toán vận tải là trờng hợp đặc biệt của quy hoạch tuyến tính, bởi vậy có thể dùng các phơng pháp của quy hoạch tuyến tính để giải. Tuy nhiên, do tính chất đặc thù riêng của nó, ngời ta xây dựng các phơng pháp giải riêng.Thông thờng khi nói đến bài toán vận tải ta thờng liên hệ ngay đến bài toán vận tải hai chỉ số, bởi đây là bài toán vận tải kinh điển có những phơng pháp giải hay. Bên cạnh đó, ngời ta còn xét một số các bài toán vận tải mở rộng nh bài toán vận tải Trang: 1 ba chỉ số, bài toán vận tải khoảng, bài toán vận tải đa mục tiêu và rất nhiều bài toán khác, đó là các biến thể của bài toán vận tải kinh điển trên.Trong khuôn khổ khoá luận này, em xem xét và nghiên cứu một số bài toán mở rộng trong lớp các bài toán vận tải mở rộng đó. Đó là các bài toán: Bài toán vận tải ba chỉ số (solid transport problem) 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 (interval solid transport problem) và giới thiệu một số Bài toán vận tải đa mục tiêu. Cuối cùng, em xin bày tỏ lòng biết ơn sâu sắc nhất đối với thày giáo hớng dẫn Thạc sỹ Vũ Tiến Việt, ngời đã tận tình chỉ bảo, giúp đỡ em trong quá trình hoàn thành khoá luận này. Em cũng xin chân thành cảm ơn sự giúp đỡ nhiệt tình của các thầy cô trong khoa Toán - Tin, Học viện An ninh Nhân dân. Trang: 2 Chơng I. Bài toán quy hoạch tuyến tínhTrong việc nghiên cứu các bài toán tối u nói chung, giải tích lồi giữ một vai trò rất quan trọng. Nó đợc sử dụng làm cơ sở toán học trong việc xây dựng các thuật toán. Quy hoạch tuyến tính là một trong những lớp bài toán tối u đợc nghiên cứu trọng vẹn cả về phơng diện lý thuyết lẫn thực hành, Bài toán vận tải là một dạng đặc biệt của QHTT. Do đó chơng này nhằm giới thiệu một số khái niệm và kiến thức cơ bản về giải tích lồi và QHTT.1.1 Một số khái niệm về giải tích lồi1.1.1 Không gian Euclude Một vector n chiều trên trờng số thực là một bộ đợc sắp thứ tự gồm n số thực x=(x1, x2, ., xn). Các xi, i =1, ., n gọi là các thành phần hay toạ độ của vector. Ví dụ x=(4,5,10,20).Hai vectơ x và y gọi là bằng nhau x=y, nếu xi=yi, i =1, ., n.Xét hai phép toán trên các vector:Phép cộng: x+y=(x1+y1, x2+y2, ., xn+yn)Phép nhân: x=(x1, x2, ., xn), RKhi đó tập hợp tất cả các vector n chiều trong đó xác định phép cộng các vector, nhân một số thực với vector nh trên tạo thành không gian tuyến tính n chiều trên tr-ờng số thực R, ký hiệu Rn.Các vector x(i) Rn, i =1, ., m đợc gọi là độc lập tuyến tính nếu:Trang: 3mixiimii,1,00)(1==== Nếu: ==miiixx1)( với ít nhất một i 0 thì x gọi là tổ hợp tuyến tính của các x(i), i =1, ., m. Hơn nữa nếu i > 0, i =1, ., m và ==mii11 thì x gọi là tổ hợp 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 cơ sở của nó. Giả sử e(1), e(2), ., e(n) là một cơ sở của Rn thì bất kỳ một vector x Rn đều là tổ hợp tuyến tính của các 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, <x,y>, là một số bằng.Tích vô hớng là một dạng song tuyến tính, đối xứng, không âm, tức là:1. <x,y> = <y, x>. x,y Rn2. <x(1) + x(2), y >=< x(1), y >+< x(2), y>. x(1), x(2), y Rn3. <x,y> = <x,y>. x,y Rn4. <x,x> 0, x Rn dấu bằng xẩy ra khi và chỉ khi x= 0.Độ dài của vector x=(x1, x2, ., xn) là một số xác định bởi.Khoảng cách giữa hai vector x và y là một số xác định bởi: Không gian vector trong đó có tích vô hớng và khoảng cách nh trên gọi là không gian Euclude.1.1.2 Tập compactDãy {x(k) }Rn k=1, 2, . đợc gọi là có giới hạn x(0) khi k và viết lim x(k) = x(0), nếuTrang: 4==><=niixxxx12,( )0,lim)0()(=xxkk==><niiiyxyx1,k ( ) ( )==><==niiiyxyxyxyxyx12,, Hình cầu tâm a bán kính là tập S={xRn :x-a }. Hình cầu này tạo nên - lân cận của điểm a, hay gọi là lân cận của a.* Nếu tập ARn chứa cùng với điểm x một lân cận của nó thì x gọi là điểm trong của A. Nếu trong lân cận bất kỳ của x A có các điểm của A và các điểm không thuộc A thì x gọi là điểm biên của tập hợp A.* Một tập ARn gọi là giới nội nếu nó đợc chứa trong một hình cầu tâm O nào đó, tức là tồn tại số đủ lớn sao cho với mọi xA,x . Một dãy {x(k)} hội tụ thì bao giờ cũng giới nội.* Một tập hợp GRn đợc gọi là mở nếu với mọi xG đều tồn tại một hình cầu tâm x nằm gọn trong G. Một tập FRn đợc gọi là đóng nếu với mọi dãy hội tụ{x(k)} F ta đều có: Fxkk)(lim Một tập chứa mọi điểm biên của nó là tập đóng.* Tập C đợc gọi là tập Compact nếu từ mọi dãy vô hạn {x(k)} thuộc C đều có thể trích ra một dãy con {x(ki)} hội tụ tới phần tử thuộc C. Tập C là Compact khi và chỉ khi C đóng và giới nội. Tập Compact M của tập đóng C cũng đóng trong C. Tập con đóng M của tập Compact cũng Compact.Hàm f(x) liên tục trên tập Compact C thì sẽ đạt cực trị trên tập ấy.1.1.3 Tập lồiCho hai điểm a, b Rn. Ta gọi đờng thẳng qua a, b là tập điểm có dạng xRn : x = a + (1-)b, R.Đoạn thẳng nối hai điểm a, b là tập lồi các điểm có dạngxRn :x = x + (1-)y, 0 1* Một tập MRn đợc gọi là một đa tạp affine nếu với hai điểm bất kỳ x, y M thì toàn bộ đờng thẳng đi qua hai điểm đó cũng thuộc M. Tức là x + (1-)y M : x,y M, R.* Một siêu phẳng trong không gian Rn là tập hợp tất cả các điểm x=(x1, x2, ., xn) Rn thỏa mãn phơng trình tuyến tính Trang: 5 a1x1+ a2x2+ . + anxn = trong đó a1, a2, ., an , R* Tập hợp các đ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 là nửa không gian đóng. * Nửa không gian đợc cho bởi a1x1+ a2x2+ . + anxn < đợc gọi là nửa không gian mở.* Tập XRn đợc gọi là tập lồi nếu cùng với việc chứa hai điểm x, y nó chứa cả đoạn thẳng chứa hai điểm ấy, tức là chứa tất cả các điểm có dạng:x + (1-)y, 0 1Ví dụ về các tập lồi: Không gian Euclide, các 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 là giao của một số hữu hạn các nửa không gian đóng đợc gọi là tập lồi đa diện. Mệnh đề: Giao của hai tập lồi là một tập lồi.Hệ quả 1. Giao của một số bất kỳ tập hợp lồi là tập lồi.Hệ quả 2. Miền chứa nghiệm của một hệ bất phơng trình tuyến tính dạng. là một tập lồi (đa diện lồi). Một tập lồi đa diện giới nội gọi là một đa diện. Giao của tất cả các 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 trên tập lồi C Rn đợc gọi là hàm lồi trên C, nếu với mọi x, y C và 0 1 ta có f(x + (1-)y) f(x) + (1-)f(y).* Hàm f(x) đợc gọi là hàm lồi chặt nếu với mọi x, y C và 0 1 ta có. f(x + (1-)y) < f(x) + (1-)f(y).* Hàm f(x) đợc gọi là hàm lõm (lõm chặt) nếu - f(x) là hàm lồi (lồi chặt)Trang: 6 22112222212111212111+++++++++nnnnnnnnnnbxaxaxabxaxaxabxaxaxa * Hàm f(x) xác định trên C đạt cực tiểu tuyệt đối tại x* C nếu f(x*) f(x): xC* Hàm f(x) đạt cực tiểu địa phơng tại x* C nếu tồn tại lân cận mở U của x* sao cho f(x*) f(x): xC UMệnh đề 1: Bất kỳ điểm cực tiểu địa phơng nào của hàm lồi trên tập lồi cũng là điểm cực tiểu tuyệt đối.Hệ quả: Bất kỳ điểm cực đại địa phơng nào của hàm lõm cũng là cực đại tuyệt đối.Mệnh đề 2: Cực đại của một hàm lồi (nếu có) trên một tập lồi có điểm cực biên bao giờ cũng đạt tại một điểm cực biên.1.2 Bài toán Quy hoạch tuyến tínhQHTT bắt nguồn từ những nghiên cứu của nhà toán học Nga nổi tiếng, Viện sỹ L.V. Kantorovich trong một loạt các công trình về bài 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 và đề xuất phơng pháp đơn hình (Simplex method) để giải bài toán QHTT. Năm 1952 ph-ơng pháp đơn hình đã đợc chạy trên máy tính điện tử của Mỹ.1.2.1 Bài toán quy hoạch tuyến tính Bài toán tổng quát.Để nhất quán lập luận ta xét bài toán tìm cực đại, sau đó ta xét cách chuyển bài toán tìm cực tiểu sang tìm cực đại. Bài toán tổng quát của QHTT có dạng:Ký hiệu: A=(aij)mxn là ma trận với các phần tử aij(1.1) gọi là hàm mục tiêu, (1.2) là các rằng buộc.Nếu gặp bài toán Min, tức làTrang: 7( )Dxxcxfjnjj==min1( )1.1max1=jnjjxc( ) ( )( )3.1, .,1,02.1 .,,1,,,1njxmibxajinjjij==== Thì giữ nguyên ràng buộc và đa về bài toán Max bằng cáchNếu bài toán Max có phơng án tối u là x* thì bài toán min cũng có phơng án là x* và fmin=-fmax Thật vậy, vì x* là phơng án tối u của bài toán Max nên ta có:Chứng tỏ x* là phơng án tối u của bài toán Min và Dạng chuẩn và dạng chính tắc.Ngời ta thờng xét bài toán quy hoạch tuyến tính dới hai dạng sau:-Dạng chuẩn:-Dạng chính tắc: Đa bài toán QHTT về dạng chuẩn hoặc dạng chính tắc.Bất kỳ QHTT nào cũng có thể đa về một trong hai dạng chuẩn hoặc chính tắc nhờ các phép biến đổi tuyến tính sau:Trang: 8( )Dxxcxfjnjj==max1DxxcxchayDxxcxcfnjjjnjjjjnjjnjjj=====,,11*11*max===njjjfxcf1max*minnjxmibxaxcjinjjijnjjj, .,1,0, .,1,max11====njxmibxaxcjinjjijjnjj, .,1,0, .1,max11===== i) Một ràng buộcCó thể đa về ràng buộc bằng cách nhân hai vế với (-1) và viết lạiii) Một ràng buộc đẳng thứccó thể thay bằng hai ràng buộc bất đẳng thức: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 bằng cách đặt:iv) Một ràng buộc bất đẳng thứcCó thể đa về ràng buộc đẳng thức bằng cách đa vào biến phụ yi 0:Về nguyên tắc, áp dụng nhiều lần các phép biến đổi (i), (ii) và (iii) ta có thể đa một bài toán QHTT bất kỳ về dạng chuẩn, sau đó áp dụng nhiều lần phép biến đổi (iv) ta sẽ đa nó về dạng chính tắc. Giải bài toán QHTT bằng phơng pháp hình học.Xét bài toán QHTT dới dạng chuẩn với hai biến số:Trang: 9=njijijbxa1injjijbxa==1injjijinjjijbxabxa==11,0,0,=++jjjjjxxxxx vớiijnjijbxa=1injijijbyxa=+=1ijnjijbxa=1.''1ijnjijbxa===+=+2,1,0, .,1,max21112211jxmibxaxaDxcxcjiii Từ ý nghĩa hình học ta biết rằng mỗi bất phơng trình tuyến tính ai1x1+ai2x2 bi xác định một nửa mặt phẳng.Nh vậy miền ràng buộc D đợc xác định nh là giao của một nửa mặt phẳng và sẽ là một đa giác lồi trên mặt phẳng. Phơng trình c1x1+c2x2= khi thay đổi sẽ xác định trên mặt phẳng các đờng thẳng song song với nhau mà ta sẽ gọi là các đờng mức (với giá trị mức ). Mỗi điểm D sẽ nằm trên một đờng mức với mức Bài toán đặt ra có thể phát biểu theo ngôn ngữ hình học nh sau: trong số các đ-ờng mức cắt tập D, hãy tìm đờng mức với gía trị lớn nhất.Nếu dịch chuyển song song các đờng mức theo hớng vector pháp tuyến của chúng thì giá trị mức sẽ tăng, nếu dịch chuyển theo hớng ngợc lại thì giá trị mức sẽ giảm. Vì vậy để giải bài toán đặt ra, ta có thể tiến hành nh sau.Bắt đầu từ một đờng mức cắt D, ta dịch chuyển song song các đờng mức theo h-ớng vector pháp tuyến (c1,c2) cho đến khi việc dịch chuyển tiếp theo làm cho đờng mức không còn cắt D nữa thì dừng. Điểm của D (có thể nhiều điểm) nằm trên đờng mức cuối cùng này sẽ là lời giải tối u cần tìm, còn giá trị của hàm mục tiêu tại đó chính là giá trị tối u của bài toán. Ví dụ: Xét bài toán: f(x)= 4x1+5x2maxXét đờng mức: 4x1+5x2=10. Đờng mức này đi qua hai điểm (0,2) và (2.5,0). Ta có x*=(3,2), fmax=22 Trang: 100,0372822122121++xxxxxxx( )21, xxx=.2211xcxc+=( )21,ccn=yn*xx [...]... -2 3 -2 4 -1 2 11 2 -1 0 -1 0 -1 8 -1 1 0 -1 7 14 -1 -8 12 12 0 7 0 0 21 28* 16 -4 -8 -1 2 -1 5 0 -1 7 5 7 0 ai 11 16 10 Ta có h = 4 và x111=6, x112= 1, x122= 4, x222= 0,x223= 9, x242= 7, x332= 4, x343= 6 *Lần lặp 3: u = (0, -6 , -2 ), v = (0, 3, -1 , 3), w= (3, 4, 9) Bảng delta 3: bj 7 4 13 13 0 0 -1 5 -2 3 -2 4 2 -3 -1 2 -1 0 -1 0 -4 -1 1 0 -3 +0 -1 5 -2 3 -2 -2 -0 -7 -1 4 0 -7 0 2 -4 -8 2* -5 0 -3 -9 -7 0 ai 11 16 10... Tìm các thế vị: u1 = 0 v2 = c12 - u1 = 2 - 0 = 2 u3 = c32 - v2 = 3 - 2 = 1 v3 = c33 - u3 = 9 - 1 = 8 u2 = c23 - v3 = 8 - 8 = 0 v1 = c21 - u2 = 1 - 0 = 1 v4 = c34 - u3 = 7 - 1 = 6 Bớc 3 Tính các ớc lợng 11 = u1 + v1 - c11 = 0 + 1 - 4 = -3 < 0 Trang: 30 13 = u1 + v3 - c13 = 0 + 8 - 10 = -2 < 0 14 = u1 + v4 - c14 = 0 + 6 - 6 = 0 22 = u2 + v2 - c22 = 0 + 2 - 3 = -1 < 0 24 = u2 + v4 - c24 = 0 + 6 - 12 = -6 ... (0, -6 , -5 ), v= (0, 3, 13, 20), w= (3, 4, -5 ) Bảng delta là: bj 7 4 13 13 0 0 -1 5 -2 3 -2 4 -1 2 -6 -1 5 -2 7 -1 0 -1 8 -1 1 0 -1 7 -3 -1 8 -5 12 12 0 7 0 0 4 11 -1 13 9 5 12 17 0 5 7 0 ai 11 16 10 Ta có h=3 và các giá trị của x tơng ứng là: x111=6, x112=4, x222=0, x232=8, x233=5, x242=5, x343=10 Trang: 35 *Lần lặp 2: u= (0, -6 , 12), v= (0, 3, 13, 3), w=(3, 4, -5 ) Bảng delta 2 : bj 7 4 13 13 0 0 -1 5 -2 3 -2 4 -1 2... nhận xét mô hình của bài toán STP là dạng tổng quát của mô hình bài toán vận tải hai chỉ số thông thờng TP (Transport Problem) nếu chúng ta chỉ nghiên cứu trong một phơng thức vận chuyển duy nhất (l=1) 2.2.2 Một số định nghĩa cơ bản i) Một bộ giá trị {xijk} i=1, , m, j=1, , n, k=1, , l thoả mãn các ràng buộc đợc gọi là một phơng án của bài toán ii) Một phơng án mà hệ thống các vector hệ số ứng với các. .. . còn xét một số các bài toán vận tải mở rộng nh bài toán vận tải Trang: 1 ba chỉ số, bài toán vận tải khoảng, bài toán vận tải đa mục tiêu và rất nhiều bài. bài toán khác, đó là các biến thể của bài toán vận tải kinh điển trên .Trong khuôn khổ khoá luận này, em xem xét và nghiên cứu một số bài toán mở rộng trong
