Phương pháp số trong lý thuyết tối ưu

Chương 2. BÀI TOÁN QUY HOẠCH TUYẾN TÍNH142.1. Thuật toán đơn hình dạng bảng162.2. Thuật toán đơn hình dạng ma trận nghịch đảo272.3. Thuật toán đơn hình đối ngẫu33Chương 3. BÀI TOÁN QUY HOẠCH PHI TUYẾN KHÔNG RÀNG BUỘC383.1. Phương pháp Gradient393.2. Phương pháp Newton423.3. Phương pháp tìm kiếm trực tiếp44Chương 4. BÀI TOÁN QUY HOẠCH PHI TUYẾN CÓ RÀNG BUỘC544.1. Phương pháp tuyến tính hóa544.2. Phương pháp hàm phạt65

BỘ GIÁO DỤC VÀ ĐÀO TẠO TRƯỜNG ĐẠI HỌC SƯ PHẠM TP HỒ CHÍ MINH Lê Đại Phước PHƯƠNG PHÁP SỐ TRONG LÝ THUYẾT TỐI ƯU https://sites.google.com/site/sitevnbookworm/ Thành phố Hồ Chí Minh - 2013 LỜI CẢM ƠN Đầu tiên, xin bày tỏ lòng biết ơn sâu sắc TS TRỊNH CÔNG DIỆU khoa Toán trường Đại học Sư phạm Thành phố Hồ Chí Minh dành thời gian công sức tận tình hướng dẫn giúp hoàn thành luận văn Tôi xin gửi lời cảm ơn đến quý Thầy Cô hội đồng chấm luận văn dành thời gian đọc, chỉnh sửa đóng góp ý kiến giúp cho hoàn thành luận văn cách hoàn chỉnh Bên cạnh đó, xin chân thành cảm ơn Ban giám hiệu trường Đại học Sư phạm Thành phố Hồ Chí Minh, Ban chủ nhiệm khoa Toán quý Thầy Cô giảng dạy, tạo điều kiện cho hoàn thành khóa học Và để có kết ngày hôm nay, nhận lời động viên, đóng góp ý kiến bạn bè người thân Cuối cùng, trình viết luận văn này, khó tránh khỏi thiếu sót, mong nhận ý kiến đóng góp bạn đọc Xin chân thành cảm ơn! Thành phố Hồ Chí Minh, tháng 10 năm 2013 Lê Đại Phước MỤC LỤC trang Trang phụ bìa Lời cảm ơn Mục lục Danh mục ký hiệu, chữ viết tắt Danh mục bảng Danh mục hình vẽ, đồ thị MỞ ĐẦU Chương KIẾN THỨC CHUẨN BỊ 1.1 Một số khái niệm kết từ giải tích lồi 1.2 Bài toán tối ưu 11 Chương BÀI TOÁN QUY HOẠCH TUYẾN TÍNH 14 2.1 Thuật toán đơn hình dạng bảng 16 2.2 Thuật toán đơn hình dạng ma trận nghịch đảo 27 2.3 Thuật toán đơn hình đối ngẫu 33 Chương BÀI TOÁN QUY HOẠCH PHI TUYẾN KHÔNG RÀNG BUỘC 38 3.1 Phương pháp Gradient 39 3.2 Phương pháp Newton 42 3.3 Phương pháp tìm kiếm trực tiếp 44 Chương BÀI TOÁN QUY HOẠCH PHI TUYẾN CÓ RÀNG BUỘC 54 4.1 Phương pháp tuyến tính hóa 54 4.2 Phương pháp hàm phạt 65 KẾT LUẬN VÀ KIẾN NGHỊ 74 TÀI LIỆU THAM KHẢO 76 DANH MỤC CÁC KÝ HIỆU, CÁC CHỮ VIẾT TẮT  tập số thực n không gian Euclid n chiều xD x thuộc tập D xD x không thuộc tập D  tập rỗng C\D hiệu tập C D CD hợp tập C tập D CD giao tập C tập D CD C tập D x, y tích vô hướng x y x chuẩn Euclid x x giá trị tuyệt đối x dom f miền xác định hữu hiệu f epi  f  epigraph hàm f f   x0 , d  đạo hàm theo hướng hàm f theo hướng d x0 f  x  vectơ gradient hàm f điểm x AT ma trận chuyển vị ma trận A A1 ma trận nghịch đảo ma trận khả nghịch A Im ma trận đơn vị cấp m rank  A  hạng ma trận A x với x x tồn x DANH MỤC CÁC BẢNG Bảng 2.1 Bảng đơn hình 16 Bảng 2.2 Bảng đơn hình dạng ma trận nghịch đảo 27 Bảng 2.3 Bảng đơn hình đối ngẫu 33 Bảng 4.1 Kết tính 71 DANH MỤC CÁC HÌNH VẼ, ĐỒ THỊ Hình 3.1 Đồ thị hàm số f  x   sin  x   sin  3x   sin  x   sin 13x  47 MỞ ĐẦU Lý chọn đề tài Có nhiều tình xã hội, từ sống đời thường đến hoạt động kinh tế, kỹ thuật, công nghệ quản lý đại… người ta phải quan tâm tới toán tìm phương án tốt để đạt mục tiêu mong muốn điều kiện ràng buộc định Đó toán tối ưu Chính nỗ lực nhằm giải toán tối ưu góp phần kích thích phát triển Giải tích Toán học kỷ XVII – XVIII với đóng góp to lớn nhà toán học lỗi lạc thời đại: Fermat, Leibniz, Euler…Tuy nhiên, phải đến năm 30, 40 kỷ XX, Quy hoạch toán học (Mathematical Programming) hay gọi Toán Tối Ưu (Optimization) hình thành với tư cách lý thuyết độc lập với nhiều hướng nghiên cứu khác nhau: Quy hoạch tuyến tính (Linear Programming), tiếp Quy hoạch lồi (Convex Programming), Quy hoạch toàn cục (Global Programming), Lý thuyết điều khiển tối ưu (Optimization Control) Danh sách tài liệu chủ đề có đến hàng trăm Nó nói lên vai trò quan trọng việc tìm cực trị toán ứng dụng Có nhiều tình thực tế việc tìm kết xác phức tạp cần kết gần sử dụng được, đồng thời cần có thuật toán hiệu cho máy tính theo xử lý Học viên quan tâm đến phương pháp chuyển đổi công thức lý thuyết đẹp đẽ thành kết số cụ thể Mục đích đề tài Vì lý trên, học viên thực luận văn khảo sát phương pháp số tiêu biểu có nhiều ứng dụng để giải vài toán quy hoạch tuyến tính phi tuyến Đối tượng phạm vi nghiên cứu Xét toán quy hoạch tuyến tính phi tuyến Với toán trình bày thuật toán, sở lý luận thuật toán, cho ví dụ minh họa Ý nghĩa khoa học thực tiễn Quy hoạch tuyến tính phận quan trọng quy hoạch toán học Đây mô hình toán nhiều toán thực tế thuộc lĩnh vực khác kinh tế, viễn thông, xây dựng… Các thuật toán giải quy hoạch tuyến tính công cụ để giải toán phức tạp quy hoạch phi tuyến, tối ưu đa mục tiêu… Ngày nay, với trợ giúp cách mạng công nghiệp thông tin, quy hoạch toán học ngày phát triển mạnh mẽ Các phương pháp tối ưu ứng dụng rộng rãi lĩnh vực khoa học, kỹ thuật, công nghệ, quản lý, kinh tế, khai thác liệu (data mining), viễn thông, v.v… Chương – KIẾN THỨC CHUẨN BỊ 1.1 Một số khái niệm kết từ giải tích lồi 1.1.1 Không gian Euclid  n 1.1.1.1 Điểm hay vectơ  n Trong luận văn làm việc không gian Euclid  n Mỗi điểm không gian  n số thực có thứ tự viết dạng cột số  x1   2 x x :       n  x  Mỗi số xi , i  1, , n , gọi tọa độ thứ i điểm x Để thuận tiện viết, ta qui ước T x   x1 , x , , x n   x1   2 x :      n  x  1.1.1.2 Tích vô hướng T  n Tích vô hướng hai vectơ x   x1 , x , , x n  y  y , y , , y T  số thực, ký hiệu x, y , xác định sau x, y : x1 y1  x y   x n y n 1.1.1.3 Chuẩn Euclid vectơ Chuẩn Euclid (hay độ dài) vectơ x   n , ký hiệu x , số thực xác định x : 2  n x, x    x i   i 1  1.1.2 Hàm nhiều biến 1.1.2.1 Định nghĩa Hàm f từ  n  quy tắc ứng điểm x thuộc  n với số thực kí hiệu số thực f  x Cách viết f : X   , với X   n , nói f  x xác định với điểm x  X Khi đó, ta gọi X miền xác định hàm f Vì biến số phần tử thuộc  n nên có n thành phần xem biến độc lập Do người ta thường gọi hàm xác định  n , với n  , hàm nhiều biến 1.1.2.2 Gradient Cho hàm f xác định tập mở X   n Giả sử x0 , đạo hàm riêng hàm f theo biến tồn Khi đó, vectơ T  f  x0  f  x0   , ,   x n   x gọi Gradient f x0 ký hiệu f  x0  1.1.2.3 Đạo hàm theo hướng Cho hàm f xác định  n vectơ d   n \ 0 Giới hạn lim t 0 f  x0  td   f  x0  t , tồn (hữu hạn vô cùng), gọi đạo hàm theo hướng d hàm f điểm x0   n ký hiệu f   x0 , d  Mệnh đề 1.1 ([1, tr 14]) Cho hàm f xác định  n điểm x0   n Nếu f khả vi x0 f   x0 , d   f  x0  , d , d   n \ 0 1.1.2.4 Hàm tuyến tính, hàm afin Một hàm số f  x  xác định  n gọi tuyến tính f   x1   x2    f  x1    f  x2  62  D3   x   | x1  0.997886 Đặt k : k   chuyển sang Bước lặp k   Bước lặp o Bước 41 Giải quy hoạch tuyến tính  f  x   x1 | x  D3  lấy T nghiệm tối ưu đại diện x4   0.997886,  o Bước 42 Vì g1  x4   0.008441142  g  x4   3.0042183  nên T Dừng thuật toán lấy nghiệm xấp xỉ x*   0.9978886,  T Bằng biểu diễn hình học, dễ thấy nghiệm tối ưu toán  1,0  Ví dụ 2: Xét toán  x  x1  với điều kiện x  D  x    x2  1  D :  x   |  x1   x  0;      Bài toán có n  2, m  2 g1  x    x   x  0, g x   x   x2    x1   x1    g1  x     , g  x      2x2   1    + Bước chuẩn bị (Xây dựng đa diện xuất phát ( D1 ): Chọn ba điểm T T T x1   0,1 , x   0, 1 x3  1,0 T  Với x1   0,1 63 0 0 g1 x1    , g x1    ; g1 x1  1; g x1   1   2                 h11  x   g1 x1 , x  x1  g1 x1   x  0; h21  x   g x1 , x  x1  g x1  x   T  Với x   1,0     g1 x  1   2  3    , g x    ; g1 x  1; g x     1                  h12  x    g1 x , x  x  g1 x  2 x1  x   0; h22  x   g x , x  x  g x  1 x  0 T  Với x3  1,0  1  2 3 g1 x    , g x    ; g1 x  1; g x     1  0                 h13  x   g1 x , x  x  g1 x  x1  x   0; h23  x   g x , x  x  g x3  1 x  0 Ta tập đa diện D1  D xác định D1 :  x   | hij  x   0, i  1, 2, j  1, 2,3 Dễ thấy ràng buộc h22  0, h23  thừa Vì ta có D1   x   | x  0; x  1; x1  x   0; x1  x   0 Đặt k : Bước lặp k   Bước lặp 64 o Bước 11 : Giải quy hoạch tuyến tính nghiệm tối ưu x1   ,0  2  o Bước 12 : Vì g1  x1    f  x  x  x1 | x  D1 T  nên chuyển sang Bước 13 o Bước 13 : Ta có g1  x1   15  g  x1   nên chọn i1  đặt 16   D2 : D1  x   | g1  x1  , x  x1  g1  x1   1   D1   x   | x1  x    4  Đặt k : k   chuyển sang Bước lặp k   Bước lặp  o Bước 21 : Giải quy hoạch tuyến tính f  x   x  x1 | x  D2  tập nghiệm tối ưu đoạn  ,    T 3 1 , ,  4 2 T    1 1 Ta lấy nghiệm tối ưu đại diện x2   ,  2 4  T o Bước 22 : Vì g1  x2   g  x2     nên Dừng thuật toán 1 1   T lấy nghiệm xấp xỉ x*   ,  65 4.2 Phương pháp hàm phạt Xét toán  f  x   (P)   gi  x   0, i  1, , m f , g i ,  i  1, , m  hàm liên tục  n Gọi D   x | gi  x   0, i  1, , m tập chấp nhận toán Trong suốt mục này, ta giả sử D tập compac Nhận xét: Nói chung toán tối ưu không ràng buộc thường dễ xử lý toán có ràng buộc Một ý tưởng nảy sinh chuyển toán có ràng buộc toán không ràng buộc Kỹ thuật để thực ý tưởng hàm phạt Hai vấn đề phải giải xây dựng hàm phạt toán phụ cho xấp xỉ lời giải toán ban đầu từ lời giải toán phụ Phương pháp hàm phạt có hai loại: Hàm phạt điểm hàm phạt điểm 4.2.1 Phương pháp hàm phạt điểm Hàm phạt p  x  thỏa mãn tính chất: i/ không âm liên tục tập int D   x   n : g i  x   0, i  1, , m ii/ p  x    gi  x   0 Hai hàm phạt điểm sử dụng nhiều, Fiacco McCormick đưa ra, là: m m i 1 i 1 p  x    ln   gi  x   p  x    gi  x  Phương pháp hàm phạt điểm xuất phát từ điểm x1 tập chấp nhận D , giải dãy toán tối ưu không ràng buộc min  x,  k  với điều kiện x   n   x,  k   f  x    k p  x  ,  k  dãy tham số dương, giảm đơn điệu 4.2.1.1 Thuật toán 66 Bước chuẩn bị: Cho số   đủ bé (để kiểm tra điều kiện dừng thuật toán) Chọn điểm x1  D thỏa mãn gi  x1   0, i  1, , m Chọn tham số phạt 1  số    0,1 Đặt k  Bước lặp k ,  k  1, 2,  Bước k1 Xuất phát từ xk , giải toán tối ưu không ràng buộc min  x,  k  với điều kiện x   n nhận nghiệm xk 1 Bước k2 : If  k p  xk 1    Then Dừng thuật toán (lấy xk 1 nghiệm tối ưu toán (P)) Else chuyển Bước k3 ; Bước k3 : Đặt tham số phạt  k 1 :  k Đặt k : k  Chuyển Bước lặp k 4.2.1.2 Sự hội tụ Định lý 4.2: Thuật toán 4.2.1.1 có tính chất sau: i/  x,   f  x  với   , x chấp nhận được; ii/ gi  xk   0, i  1, , m   iii/   x k , k  hội tụ đến giá trị tối ưu toán (P)  k   điểm tụ  xk  nghiệm tối ưu (P) Chứng minh: i/ Hiển nhiên ii/ Nếu gi  xk   , ta suy   xk , k  không hữu hạn iii/ Nếu x* nghiệm (P), ta cần chứng minh lim  xk , k   f  x*  k  67 Vì  k   0,   xk , k  dãy giảm  bị chặn (do f liên tục tập compac D nên bị chặn  k p  x   ), nên ta suy   xk , k  hội tụ tới   f  x*  Thật vậy:   xk , k   f  xk   f  x*    lim   xk , k   f  x*  k  Nếu     f  x*   , f liên tục nên ta chọn x cho  f x  f  x*     g i x  0, i  1, , m Lấy k đủ lớn cho m  k  i 1   gi x    xk , k     Ta có     xk , k  (vì   xk , k  dãy giảm hội tụ tới  )    x, k  (vì   x, k  đạt cực tiểu xk ) m   f x k  i 1   f x   gi x   f  x*     f  x*     f  x*    68 Điều mâu thuẫn Do đó,     f  x*   , tức   f  x*  Tiếp theo, x điểm tụ  xk  , tức x  lim xki , theo i/   xki ,  ki   f  xki  Từ suy    lim  xki , ki   f x Điều có nghĩa x cực tiểu □ Chú ý: Điểm cực tiểu   x,   nằm tập chấp nhận 4.2.1.3 Ví dụ Trước xét ví dụ, ta nhắc lại vài tính chất hàm lồi Mệnh đề 4.2: Cho hàm lồi f :  n   tập lồi khác rỗng D   n Xét toán  f  x  | x  D Khi đó: i/Nếu x* nghiệm tối ưu địa phương toán x* nghiệm tối ưu toàn cục; ii/Nếu x* nghiệm tối ưu địa phương chặt f hàm lồi chặt x* nghiệm tối ưu toàn cục toán Chứng minh: Xem [1, tr 55] Định lý 4.3: Giả sử f hàm lồi khả vi  n Khi x*   n nghiệm cực tiểu toàn cục toán  f  x  | x   n  , f :  n   hàm phi tuyến f  x*   Chứng minh: Xem [1, tr 222] Ví dụ: Xét toán 69  f  x   x1  x    2  g1  x   1  x   x     g  x    x  Ta chọn hàm phạt có dạng   p  x    ln   g1  x    ln   g  x     ln  x1   x   ln  x  Khi đó, hàm mục tiêu toán không ràng buộc tương ứng     x,    x1  x   ln  x1   x    ln  x  Trong  tham số dương giảm đơn điệu Với tham số  cụ thể, toán không ràng buộc   x,   | x   n  có hàm mục tiêu   x,   hàm lồi chặt theo biến x nên, theo mệnh đề định lý trên, có nghiệm cực tiểu điểm dừng Ta có  x1  x,      x1   x   x  x,    2  , 2 x  x1   x 2    x2 T Các điểm x   x1 , x  thuộc phần tập chấp nhận có  x1   x   x2  Giải hệ phương trình   0  1 2  x  x     x  x,     2 x   2   0  x  x1   x   T Ta nhận điểm dừng x*     x*1   , x*2    với 70  2  3    2 , x*    2 x*1    Vì tham số   giảm dần nên lim x*1     0  30   0 lim x*2       0  0  0, 1 T Ta có điểm  0,1 thuộc tập chấp nhận nghiệm tối ưu toán cần giải Bảng 4.1 trình bày giá trị x*1   x*2   với tham số 1  1,   101  k 1   k , k  1,2, ,6  Rõ ràng x*  k    x*1  k  , x*2  k   T  hội tụ đến nghiệm tối ưu toán xét Bảng 4.1 Kết tính k k x*1  k  x*2  k  100 1.8660254 1.3660254 10 1 0.1977226 1.0477226 10 2 0.0199752 1.0049752 10 3 0.0019998 1.0004998 10 4 0.0002000 1.0000500 10 5 0.0000200 1.0000050 4.2.2 Phương pháp hàm phạt điểm Hàm phạt p  x  định nghĩa m p  x     g i  x   i 1  hàm biến liên tục thỏa mãn 71   y   y    y   y  Thông thường, hàm phạt p  x  có dạng m m i 1 i 1 p  x  :  max 0, gi  x  p  x  :   max 0, gi  x  Hàm mục tiêu dãy toán tối ưu không ràng buộc tương ứng với toán (P)   x, k   f  x    k p  x  Trong dãy tham số  k  dãy số dương, đơn điệu tăng đến  Đại lượng  k p  x  lượng phạt Dễ thấy x phương án chấp nhận toán (P) không bị phạt (lượng phạt 0), x phương án chấp nhận phải chịu lượng phạt  k p  x  4.2.2.1 Thuật toán Bước chuẩn bị: Cho số   đủ bé (để kiểm tra điều kiện dừng thuật toán) Chọn điểm x1   n Chọn tham số phạt 1  số   Đặt k  ; Bước lặp k ,  k  1, 2,  Bước k1 : Giải toán tối ưu không ràng buộc   x, k  với điều kiện x   n Gọi nghiệm toán xk 1 Bước k2 If  k p  xk 1    Then Dừng thuật toán (lấy xk 1 nghiệm tối ưu toán (P)) Else Chuyển Bước k3 ; Bước k3 Đặt tham số phạt  k 1 :  k Đặt k  k  chuyển Bước k 4.2.2.2 Sự hội tụ 72 Định lý 4.4: Thuật toán 4.2.2.1 có tính chất sau: i/   x,    f  x  với   0, x   n ii/   xk , k  tăng,  k tăng; iii/   xk , k  hội tụ tới giá trị tối ưu (P),  k   , điểm tụ  xk  nghiệm tối ưu (P) Chứng minh Chứng minh tương tự định lý 4.2 4.2.2.3 Ví dụ Xét toán  f  x   x1  x    g1  x    x1   x     g  x    x1   x   8  g  x  x    Ta chọn hàm phạt có dạng p  x   max 0, g1  x   max 0, g  x   max 0, g  x  Khi hàm mục tiêu toán không ràng buộc tương ứng   x, k   f  x    k p  x  2   7    x1  x    k  max 0,  x1   x  max 0,  x1   x    max 0, x  3  8      Trong dãy tham số  k  dãy số dương, đơn điệu tăng đến  Bước chuẩn bị: Chọn giá trị   0.0005 x1   0,  1  10   10 73 Giải toán   x, 1  với điều kiện x  2 thuật toán 3.3.2.1 x11   -0.8, 1.2 x12   -0.96, 1.04 x13   -0.992, 1.008 x14   -0.9984, 1.0016 x15   -0.99968, 1.00032 x16   -0.999936, 1.000064 x1*   -0.9999872, 1.0000128  x2 Vì 1 p  x2     nên dừng thuật toán Lấy x*   -0.9999872, 1.0000128 nghiệm tối ưu toán f  x*   2.56.105  74 KẾT LUẬN VÀ KIẾN NGHỊ Luận văn trình bày vài phương pháp giải toán tối ưu Cụ thể sau: Chương – BÀI TOÁN QUY HOẠCH TUYẾN TÍNH (Tham khảo [2]) Mục Một vài kết 2.1 Thuật toán đơn Thuật toán: Trình bày lại cho gọn hơn, kết hợp thuật toán hình dạng bảng đơn hình hai pha để tìm phương án sở chấp nhận xuất phát Chương trình: Viết lại Matlab, bổ sung phần xử lý trường hợp xảy chuyển từ pha thứ sang pha thứ hai 2.2 Thuật toán đơn Thuật toán: Trình bày lại cho gọn hình dạng ma trận Chương trình: Dựa theo thuật toán mô tả để lựa chọn nghịch đảo cấu trúc liệu phù hợp viết chương trình mô 2.3 Thuật toán đơn Matlab hình đối ngẫu Chương – BÀI TOÁN QUY HOẠCH PHI TUYẾN KHÔNG RÀNG BUỘC (Tham khảo [4]) Mục Một vài kết 3.1 Phương pháp Chương trình: Gradient Viết lại Matlab nên mở rộng việc tính đạo hàm, 3.2 Phương pháp Matlab hỗ trợ đại số máy tính, phạm vi áp dụng Newton chương trình rộng [4] Matlab hỗ trợ phép tính ma trận nên chương trình gọn 3.3 Phương pháp Phần học viên xây dựng thuật toán cho trường hợp n = 1, tìm kiếm trực tiếp n=2 để giải ví dụ phương pháp hàm phạt điểm chương 75 Chương – BÀI TOÁN QUY HOẠCH PHI TUYẾN CÓ RÀNG BUỘC (Tham khảo [1], [3], [5]) Một vài kết + nhận xét Mục 4.1 Phương pháp Chứng minh: toán (P1), (P2) tương đương, dãy X k tuyến tính hóa giảm chứa X , X tập lồi đa diện X  X , hội tụ thuật toán Xây dựng thêm ví dụ Bài toán phụ mục quy hoạch tuyến tính nên ta sử dụng thuật toán chương để giải 4.2 Phương pháp Tìm hiểu, so sánh, tổng hợp từ [1], [3], [5] trình bày lại cách chi tiết Xây dựng thêm ví dụ hàm phạt 4.2.1 Phương pháp Bài toán phụ mục quy hoạch phi tuyến không ràng hàm phạt điểm buộc nên ta sử dụng thuật toán chương để giải 4.2.2 Phương pháp hàm phạt điểm Do thời gian có hạn trình độ thân người thực luận văn hạn chế nên luận văn trình bày vài phương pháp để giải toán quy hoạch tuyến tính phi tuyến Ngoài nhiều phương pháp khác mà có điều kiện, tiếp tục nghiên cứu trình bày dịp khác Chúng mong góp ý từ Quý Thầy, Cô hội đồng chấm luận văn để luận văn hoàn thiện 76 TÀI LIỆU THAM KHẢO Tiếng Việt Nguyễn Thị Bạch Kim (2008), Giáo trình Các Phương pháp Tối ưu Lý thuyết Thuật toán, Nxb Bách khoa, Hà Nội Nguyễn Đức Nghĩa (1999), Tối ưu hóa (Quy hoạch tuyến tính rời rạc), Nxb Giáo dục, Hà Nội Nguyễn Xuân Tấn, Nguyễn Bá Minh (2007), Lý thuyết tối ưu không trơn, Nxb Đại học Quốc gia Hà Nội, Hà Nội Mai Đức Thanh (2002), Một vài phương pháp số giải toán tối ưu phi tuyến, Luận văn tốt nghiệp đại học chuyên ngành Toán, Đại học Sư phạm Thành phố Hồ Chí Minh, Hồ Chí Minh Tiếng Anh B.N.Pshenichny and Yu.M.Danilin (1982), Numerical methods in extremal problems, Mir Publishers, Moscow ... j  0, j  1,8    Đầu vào: [nofeas,nobounded,x,fx]=QHTT_Bang(3,8,[ -2 -2 -1 0 ;-1 -3 -1 -1 0;3 -3 -4 -1 0 1],[1 3], [-2 -1 0 0]) 26 Đầu ra: nobounded = Kết luận: hàm mục tiêu không bị chặn Ví... j  1,    Đầu vào: [nofeas,nobounded,x,fx]=QHTT_Bang(3,4,[1 -1 2;2 -2 ;3 0],[1 3], [-1 -1 ]) Đầu ra: nofeas = nobounded = x= f= -1 0 0 Kết luận: Phương án tối ưu x*  1, 0,0,  Giá trị tối... 18   x j  0, j  1,    Đầu vào: [nofeas,nobounded,x,fx]=QHTT_Bang(3,4,[ -2 0;2 0;1 -1 -1 1],[27 50 18],[2 -2 0]) Đầu ra: nofeas = Kết luận: toán phương án Ví dụ 2: Giải toán quy hoạch tuyến