(Luận văn) phương pháp số giải bài toán quy hoạch lồi và ứng dụng

ĐẠI HỌC THÁI NGUYÊN TRƯỜNG ĐẠI HỌC KHOA HỌC  - TRƯƠNG TUẤN HƯNG lu an n va p ie gh tn to PHƯƠNG PHÁP SỐ GIẢI BÀI TOÁN QUY HOẠCH LỒI VÀ ỨNG DỤNG d oa nl w ll u nf va an lu oi m LUẬN VĂN THẠC SĨ TOÁN HỌC z at nh z m co l gm @ an Lu THÁI NGUYÊN - 2018 n va ac th si ĐẠI HỌC THÁI NGUYÊN TRƯỜNG ĐẠI HỌC KHOA HỌC  - TRƯƠNG TUẤN HƯNG lu an n va p ie gh tn to PHƯƠNG PHÁP SỐ GIẢI BÀI TOÁN QUY HOẠCH LỒI VÀ ỨNG DỤNG d oa nl w Chuyên ngành: Toán ứng dụng Mã số : 8460112 an lu ll u nf va LUẬN VĂN THẠC SĨ TOÁN HỌC oi m NGƯỜI HƯỚNG DẪN KHOA HỌC (Xác nhận) z at nh z m co l gm @ TS Vũ Vinh Quang an Lu THÁI NGUYÊN - 2018 n va ac th si iii Mục lục v Bảng ký hiệu Mở đầu lu Lời cảm ơn an n va gh tn to p ie Chương MỘT SỐ KIẾN THỨC CƠ BẢN 1.1 Mơ hình tổng qt tốn quy hoạch tuyến tính Mơ hình tổng quát Phân loại toán tối ưu Bài tốn quy hoạch tuyến tính Một số phương pháp giải d va an lu 1.2 1.3 oa nl w 1.1.1 1.1.2 Thuật tốn hình học Thuật tốn đơn hình 1.3.3 1.3.4 Thuật tốn đơn hình mở rộng 15 Phương pháp giải toán quy hoạch tuyến tính oi lm ul nf 1.3.1 1.3.2 z at nh tổng quát phần mềm MATLAB 16 z gm @ Chương BÀI TOÁN QUY HOẠCH LỒI, CÁC THUẬT 2.1.1 Khái niệm tập lồi, hàm lồi 18 2.1.2 2.1.3 Khái niệm Gradient đạo hàm theo hướng 20 Bài toán quy hoạch lồi tổng quát, điều kiện tối ưu 21 an Lu va 2.2 m co l TỐN 18 2.1 Mơ hình toán quy hoạch lồi tổng quát 18 n Cực tiểu hàm lồi biến 22 2.2.1 Thuật tốn chia đơi 22 ac th si iv 2.3 2.2.2 Thuật toán mặt cắt vàng 24 Mơ hình tốn quy hoạch lồi với ràng buộc tuyến tính 26 2.3.1 2.3.2 2.4 Mơ hình tổng qt 26 Thuật toán Frank-Wolfe 26 Mơ hình tốn quy hoạch lồi với ràng buộc phi tuyến 29 2.4.1 Mơ hình tổng quát 29 2.4.2 Thuật toán Gradient 29 Chương MỘT SỐ ỨNG DỤNG THIẾT KẾ TỐI ƯU lu an 3.1 3.2 32 va Mơ hình toán sản xuất sản phẩm 32 Mô hình tốn xác định thiết diện tối ưu giàn chịu n lực 36 39 p ie gh tn to Kết luận 40 d oa nl w Tài liệu tham khảo oi lm ul nf va an lu z at nh z m co l gm @ an Lu n va ac th si v Lời cảm ơn lu Trước hết, tơi xin bày tỏ lịng kính trọng lịng biết ơn sâu sắc tới thầy giáo TS Vũ Vinh Quang, người thầy tận tình hướng dẫn, bảo an n va cung cấp tài liệu hữu ích để tơi hồn thành luận văn gh tn to Xin cảm ơn lãnh đạo Trường Đại học Khoa học - Đại học Thái nguyên tạo điều kiện giúp đỡ mặt suốt trình học tập p ie thực luận văn oa nl w Tôi xin bày tỏ lịng biết ơn tới thầy, giáo giảng dạy lớp K10Y truyền đạt kiến thức, phương pháp nghiên cứu khoa học d suốt năm học vừa qua Xin chân thành cảm ơn anh chị em học viên cao học K10Y an lu ul nf va bạn đồng nghiệp động viên, khích lệ tơi q trình học tập, nghiên cứu oi lm Tơi xin bày tỏ lịng biết ơn sâu sắc đến gia đình, người thân, người ln động viên, khuyến khích giúp đỡ mặt để tơi z at nh hồn thành cơng việc nghiên cứu z m co l gm @ an Lu n va ac th si vi Lời cam đoan lu Tôi xin cam đoan: Những nội dung luận văn thực hướng dẫn trực tiếp thầy giáo hướng dẫn TS Vũ Vinh an n va Quang Mọi tham khảo dùng luận văn trích dẫn rõ ràng tác giả, tên cơng trình, thời gian, địa điểm cơng bố p ie gh tn to Tôi xin chịu trách nhiệm với lời cam đoan d oa nl w oi lm ul nf va an lu z at nh z m co l gm @ an Lu n va ac th si Bảng ký hiệu an n va X Vectơ không gian Rn f (X) → max f (X) → Bài tốn tìm cực đại Bài tốn tìm cực tiểu CT D Vectơ chuyển vị C Miền phương án ∇f |x| Vectơ đạo hàm hướng Giá trị tuyệt đối x gh tn to Tập số thực Không gian vectơ thực n chiều ie lu R Rn p oa nl w x không thuộc D Ma trận chuyển vị ma trận A ul nf va an lu x∈ /D AT Số phần tử tập X x thuộc D d |X| x∈D Là điểm để hàm f (x) đạt giá trị tối ưu f val Exitf lag Giá trị hàm mục tiêu Số ngun thơng báo kết thúc tính toán lb(lower bound) ub(upper bound) Giới hạn Giới hạn linprog bintprog Lệnh để lấy nghiệm không âm Lệnh để lấy nghiệm nguyên có giá trị QHT T Quy hoạch tuyến tính oi lm xopt z at nh z m co l gm @ an Lu n va ac th si Mở đầu lu Mơ hình tốn quy hoạch phi tuyến tính nói chung quy hoạch lồi nói riêng mơ hình quan trọng lớp tốn tối ưu an n va hóa, có nhiều ứng dụng toán học vật lý Về mặt lý thuyết, có nhiều tài liệu trình bày thuật tốn lý gh tn to thuyết giải mơ hình tốn mơ hình tổng quát Tuy nhiên việc nghiên cứu cài đặt chi tiết thuật toán ứng dụng vào p ie số mơ hình tốn cụ thể học vật lý chưa oa nl w nhiều người đề cập đến Nội dung luận văn nghiên cứu sở toán học d thuật toán giải tốn quy hoạch lồi có ràng buộc, tìm hiểu chi tiết bước mơ tả thuật tốn, xây dựng sơ đồ khối cài đặt an lu ul nf va thuật tốn ngơn ngữ lập trình cụ thể Trên sở thuật toán nghiên cứu cài đặt, luận văn xây dựng mơ hình ứng dụng số oi lm toán học vật lý xác định mơ hình tối ưu thiết kế Nội dung luận văn dự kiến gồm có chương, phần phụ lục z at nh cấu trúc sau: Chương 1: Trình bày số kiến thức bao gồm mơ hình tổng z @ qt tốn quy hoạch tuyến tính, thuật tốn: Hình học, đơn l gm hình, đơn hình mở rộng phương pháp giải tốn quy hoạch tuyến tính tổng quát phần mềm MATLAB m co Chương 2: Trình bày kiến thức thuật toán liên quan đến toán quy hoạch lồi bao gồm mơ hình tốn quy hoạch lồi tổng quát, an Lu n va thuật toán giải tốn cực tiểu hàm lồi biến, mơ hình tốn quy hoạch lồi với ràng buộc tuyến tính, thuật tốn Frank−Wolfe ac th Mơ hình tốn quy hoạch lồi với ràng buộc phi tuyến, thuật toán si Gradient Chương 3: Trình bày số mơ hình tốn ứng dụng thực tế: Mơ hình tốn sản xuất sản phẩm mơ hình tốn xác định thiết diện tối ưu giàn chịu lực Phần phụ lục đưa số chương trình nguồn mơi trường MATLAB giải tốn cực trị hàm lồi Thái Nguyên, tháng năm 2018 Tác giả luận văn lu an n va p ie gh tn to Trương Tuấn Hưng d oa nl w oi lm ul nf va an lu z at nh z m co l gm @ an Lu n va ac th si Chương MỘT SỐ KIẾN THỨC CƠ BẢN lu an n va Nội dung chương trình bày mơ hình tổng qt tốn quy hoạch tuyến tính, thuật tốn: Hình học, đơn hình, đơn gh tn to hình mở rộng, phương pháp giải tốn quy hoạch tuyến tính tổng p ie quát phần mềm MATLAB Các kết kiến thức quan trọng ứng dụng chương sau luận văn Các kiến thức Mơ hình tổng qt tốn quy hoạch tuyến nf Mơ hình tổng qt oi lm ul 1.1.1 va tính an lu 1.1 d oa nl w tham khảo tài liệu [1, 2] Tối ưu hóa lĩnh vực quan trọng tốn học có z at nh ảnh hưởng đến hầu hết lĩnh vực khoa học, công nghệ kinh tế xã hội Việc tìm giải pháp tối ưu cho tốn thực tế chiếm z gm @ vai trò quan trọng việc tiến hành lập kế hoạch sản xuất hay thiết kế hệ thống điều khiển trình Nếu sử dụng kiến n va Bài toán tối ưu tổng quát phát biểu sau: Cực đại hóa (cực tiểu hóa) hàm: an Lu vấn đề đặt thực tế m co l thức tảng toán học để giải toán cực trị, người ta đạt hiệu kinh tế cao Điều phù hợp với mục đích ac th f (X) → max(min) si 31 lu an n va p ie gh tn to d oa nl w oi lm ul nf va an lu z at nh z m co l gm @ an Lu n va ac th si 32 Chương MỘT SỐ ỨNG DỤNG THIẾT lu an KẾ TỐI ƯU n va p ie gh tn to Nội dung chương đưa số mơ hình tốn thực tế mơ hình hóa dạng tốn quy hoach lồi tổng qt Lời giải Mơ hình tốn sản xuất sản phẩm va an lu 3.1 d oa nl w toán tìm từ thuật tốn trình bày chương nf Bài toán: Một xưởng sản xuất hai loại sản phẩm A B với hai nguồn oi lm ul nguyên liệu Lượng tối đa loại nguyên liệu MA loại nguyên liệu MB đơn vị ngày Để sản xuất đơn vị A cần pA1 đơn z at nh vị nguồn pA2 đơn vị nguồn Một đơn vị B cần pB1 z đơn vị nguồn pB2 đơn vị nguồn Giá đơn vị nguồn nguồn cho quan hệ tương ứng (a-bu1) (c-du2), @ sB = fb − gb xA − hb xB an Lu sA = fa − ga xA − xB ; m co l gm ui số đơn vị nguồn i sử dụng Giá bán đơn vị sản phẩm A B là: n va ac th si 33 Trong xA , xB số đơn vị A B bán Hãy xây dựng phương án sản xuất cho lợi nhuận xưởng lớn nhất, giả thiết sản phẩm bán hết Xây dựng mơ hình tốn học Gọi xA , xB số đơn vị A B sản xuất ngày Yêu cầu ngày nguồn nguồn tương ứng : (pA1xA + pA2xB ); (pB1xA + pB2xB ) lu an Khi theo điều kiện tốn, ta có hệ ràng buộc: va n    pA1xA + pA2xB ≤ M A;      pB1x + pB2x ≤ M B; A B   xA ≥ 0;      x ≥ 0; B p ie gh tn to nl w d oa Chi phí cho nguồn nguồn ngày là: an lu (pA1xA + pA2xB ) [(a − b (pA1xA + pA2xB )] nf va + (pB1xA + pB2xB ) [c − d (pB1xA + pB2xB )] oi lm ul Tiền bán sản phẩm ngày: xA (fa − ga xA − xB ) + xB (fb − gb xA − hb xB ) z at nh Lợi nhuận tổng cộng tiền bán sản phẩm trừ chi phí Hàm mục tiêu xác định công thức z @ gm f (x) =(pA1xA + pA2xB ) [(a − b (pA1xA + pA2xB )] l + (pB1xA + pB2xB ) [c − d (pB1xA + pB2xB )] an Lu → m co − xA (fa − ga xA − xB ) − xB (fb − gb xA − hb xB ) va Nhận xét: Bài tốn dạng tốn cực trị hàm lồi với n ràng buộc tuyến tính, lời giải toán xác định thuật toán Frank-Wolf Sau kết toán với số liệu cụ thể: ac th si 34 Bài toán: Một xưởng sản xuất hai loại sản phẩm A B với hai nguồn nguyên liệu Lượng tối đa nguồn nguồn tương ứng 1000 250 đơn vị ngày Để sản xuất đơn vị A cần đơn vị nguồn 0,2 đơn vị nguồn Một đơn vị B cần 0,5 đơn vị nguồn 0,5 đơn vị nguồn Giá đơn vị nguồn nguồn cho quan hệ tương ứng (0,375-0,00005u1)và (0,75-0,0001u2), ui số đơn vị nguồn I sử dụng Giá bán đơn vị sản phẩm A B là: lu pA = 2, 00 − 0, 0005xA − 0, 00015xB ; an n va pB = 3, 50 − 0, 0002xA − 0, 0015xB p ie gh tn to Trong xA , xB số đơn vị A B bán Hãy xây dựng phương án sản xuất cho lợi nhuận xưởng lớn nhất, giả thiêt nl w sản phẩm bán hết Xây dựng mơ hình tốn học Gọi xA , xB số đơn vị A B sản xuất ngày Yêu cầu ngày nguồn d oa nguồn tương ứng : (xA + 0, 5xB ) (0, 2xA + 0, 5xB ) Các ràng buộc là:    xA + 0, 5xB ≤ 1000;      0, 2x + 0, 5x ≤ 250; B ul nf va an lu A oi lm   xA ≥ 0;      x ≥ B z at nh Chi phí cho nguồn nguồn ngày là: z gm @ (xA +0, 5xB [0, 375 − 0, 00005(xA + 0, 5xB )] + (0, 2xA + 0, 5xB [0, 75 − 0, 0001(0, 2xA + 0, 5xB )] m co l Tiền bán sản phẩm ngày: an Lu xA (2, 00 − 0, 0005xA − 0, 00015xB )+xB (3, 50 − 0, 0002xA − 0, 00015xB ) Lợi nhuận tổng cộng tiền bán sản phẩm trừ chi phí n va ac th si 35 Hàm mục tiêu xác định công thức f (x) = (xA + 0, 5xB ) [0, 375 − 0, 00005(xA + 0, 5xB )] + (0, 2xA + 0, 5xB ) [0, 75 − 0, 0001(0, 2xA + 0, 5xB )] − xA (2, 00 − 0, 0005xA − 0, 00015xB ) − xB (3, 50 − 0, 0002xA − 0, 00015xB ) → Như ta thu toán lu an f (xA , xB ) = c1 x2A + c2 xA xB + c3 x2B + c4 xA + c5 xB →    xA + 0, 5xB ≤ 1000;      0, 2x + 0, 5x ≤ 250; n va tn to A B gh p ie   xA ≥ 0;      x ≥ B w d oa nl Nghiệm tối ưu tốn tìm từ chương trình Frank_wolfee.m oi lm ul nf va an lu Xopt = (1000; 0); fmin = −69995 z at nh z m co l gm @ an Lu n va ac th si 36 3.2 Mơ hình toán xác định thiết diện tối ưu giàn chịu lực Xét giàn chịu tác dụng lực P (Hình vẽ) lu an n va p ie gh tn to d oa nl w an lu ul nf va Khoảng cách a=10 Gọi vectơ chuyển vị nút z = [z1 ; z2 ]T , oi lm Kí hiệu độ dài b1 , b2 , b3 ; b : [b1 ; b2 ; b3 ] Kí hiệu E mơđun đàn hồi, σia ứng suất cho phép i=1,2,3 z at nh Yêu cầu: xác định thiết diện ngang độ dài b1 , b2 , b3 (Các biến thiết kế) cho trọng lượng giàn nhỏ nhất, đồng thời z @ thỏa mãn điều kiện ứng suất m co l gm Mơ hình tốn học Theo lý thuyết sức bền vật liệu, trọng lượng giàn xác định theo công thức : an Lu  √ √  ψ0 = ρg 10 2b1 + 10b2 + 10 2b3 n va Theo lý thuyết đàn hồi, phương trình trạng thái giàn có dạng: (3.1) ac th K(b)z − S = 0; si 37 Trong K(b) ma trận độ cứng Đây ma trận xác định dương S vectơ tải trọng Ma trận độ cứng vectơ trọng tải xác định công thức: √ " 2E (b1 + b3 ) K(b) = 40 (b1 − b3 ) " # P cos θ S= P sin θ (b1 − b3 ) √
b1 + b3 + 2b2 # (3.2) (3.3) lu an Các ứng suất xác định bởi: E(z1 + z2 ) Ez2 E(z2 + z1 ) σ1 = ; σ2 = ; σ3 = 20 20 20 Như vậy, có ràng buộc ứng suất n va tn to E |z1 + z2 | − σ1a ≤ 0; 20 E |z2 | − σ2a ≤ 0; g2 = 10 E g1 = |z2 − z1 | − σ3a ≤ 20 p ie gh g1 = (3.4) (3.6) d oa nl w (3.5) an lu Mơ hình tốn giải cách áp dụng phương ul nf va pháp giải toán quy hoạch phi tuyến bị ràng buộc Để đơn giản ta xét trường hợp: oi lm b1 = b3 ; ≤ θ ≤ 900 ; P = 20000lb; σ1a = σ2a = 20000lb/in2 Khi " phương#trình " # (3.1) có"dạng: # √ b1 z1 10 P cos θ √ = E b1 + 2b2 z2 P sin θ Giải phương trình ta được: z at nh z gm @ an Lu (s1 = P cos θ; s2 = P sin θ) m co l √ √ 10 2s1 10 2s2 √
z1 = ; z2 = b1 E b1 + 2b2 E va n Khi bỏ qua đại lượng 10σg = const , ta nhận hàm mục tiêu có ac th si 38 dạng: √ ψ0 = 2b1 + b2 (3.7) Các hệ ràng buộc: √   s1 s2 √ − 2000 ≤ 0; g1 = + b1 b1 + 2b2 √ 2s2 √ g2 = − 2000 ≤ b1 + 2b2 (3.8) (3.9) lu Chúng ta chứng minh mơ hình tốn (3.7)-(3.9) mơ hình toán quy hoạch lồi phi tuyến an va n Thậy vậy: + Hàm mục tiêu ψ0 tuyến tính biến b1 , b2 hàm p ie gh tn to mục tiêu hàm lồi miền xác định + Ma trận đạo hàm cấp g1 , g2 là:  √
s1 √  √ + + 2b2 /b1 2s2   s √
3 √ b1 + 2b2 2 # " √ √ 2 2s2 √ √
3 2 b1 + 2b2 d oa nl w oi lm ul nf va an lu z at nh Ma trận đầu xác định dương, ma trận hai bán xác định dương, hàm g1 , g2 lồi Nghiệm tối ưu tốn tìm từ thuật tốn gradient z Xopt = (7.8868; 4.8882); fmin = 26.3896 m co l gm @ an Lu n va ac th si 39 Kết luận lu Nội dung luận văn nghiên cứu mơ hình tốn quy hoạch tuyến tính quy hoạch lồi, kết luận văn đạt an n va được: Trình bày thuật tốn giải mơ hình tốn quy hoạch gh tn to tuyến tính bao gồm: Thuật tốn hình học, thuật tốn đơn hình Nghiên cứu mơ hình tốn quy hoạch lồi tổng qt p ie Nghiên cứu thuật tốn tìm nghiệm tối ưu hàm lồi biến oa nl w số, cài đặt thuật tốn ngơn ngữ lập trình MATLAB Nghiên cứu mơ hình tốn quy hoạch lồi với ràng buộc tuyến d tính, thuật tốn Frank-Wolfe tìm nghiệm tối ưu Cài đặt thuật tốn ngơn ngữ lập trình MATLAB an lu ul nf va Nghiên cứu mơ hình tốn quy hoạch lồi tổng quát, thuật toán Gradient giải toán sở phối hợp thuật toán Frank-Wolfe với oi lm thuật tốn tìm cực trị hàm biến Tìm hiểu số mơ hình tốn thực tế, xây dựng mơ z at nh hình tốn học giải toán việc sử dụng thuật tốn trình bày Hướng phát triển luận văn tiếp tục nghiên z m co l gm tốn tìm nghiệm tối ưu @ cứu mơ hình tốn cực trị phi tuyến tổng quát thuật an Lu n va ac th si 40 Tài liệu tham khảo Tiếng Việt lu an [1] Bùi Minh Trí (1999), Quy hoạch tốn học, NXB Khoa học Kĩ n va thuật khiển tối ưu, Nhà xuất Khoa học Kĩ thuật p ie gh tn to [2] Nguyễn Nhật Lệ (2009), Các tốn tối ưu hóa điều nl w Tiếng Anh d oa [3] J Cea (1978), Lectures on Optimization – Theory and Algorithms an lu Tata Institute of Fundamental Ressarch, Springer nf va [4] R Fletcher (2000), Practical Methods of Optimization, Wiley oi lm ul [5] C T Kelley (1999), Iterative Methods for Optimization, SIAM z at nh [6] V J Nocedal, S J Wright (1999), Numerical Optimization, Springer z m co l gm @ an Lu n va ac th si 41 PHẦN PHỤ LỤC (Một số chương trình nguồn Matlab) lu an Thuật tốn chia đơi n va p ie gh tn to d oa nl w oi lm ul nf va an lu z at nh z m co l gm @ an Lu n va ac th si 42 Thuật toán mặt cắt vàng lu an n va p ie gh tn to d oa nl w oi lm ul nf va an lu z at nh z m co l gm @ an Lu n va ac th si 43 Thuật toán Frank −Wolfe lu an n va p ie gh tn to d oa nl w oi lm ul nf va an lu z at nh z m co l gm @ an Lu n va ac th si 44 Thuật toán Gradient lu an n va p ie gh tn to d oa nl w oi lm ul nf va an lu z at nh z m co l gm @ an Lu n va ac th si 45 lu an n va p ie gh tn to d oa nl w oi lm ul nf va an lu z at nh z m co l gm @ an Lu n va ac th si