ĐẠI HỌC THÁI NGUYÊN TRƯỜNG ĐẠI HỌC KHOA HỌC PHẠM THỊ LINH PHƯƠNG PHÁP DƯỚI ĐẠO HÀM CHO BÀI TOÁN CÂN BẰNG VÀ ÁNH XẠ GIẢ CO CHẶT Chuyên nghành: TOÁN ỨNG DỤNG Mã số: 60.46.01.12 LUẬN VĂN THẠC SĨ TOÁN HỌC Người hướng dẫn khoa học: PGS.TS. PHẠM NGỌC ANH Thái Nguyên - 2012 1Số hóa bởi Trung tâm Học liệu – Đại học Thái Nguyên http://www.lrc-tnu.edu.vn i Mục lục Mục lục . . . . . . . . . . . . . . . . . . . . . . . . . . . . . i Lời cảm ơn . . . . . . . . . . . . . . . . . . . . . . . . . . . . ii Những kí hiệu và chữ viết tắt . . . . . . . . . . . . . . . . iii Lời nói đầu . . . . . . . . . . . . . . . . . . . . . . . . . . . 1 Chương 1. Một số Khái niệm Cơ bản 2 1.1. Tập lồi và các phép toán cơ bản . . . . . . . . . . . . . . 2 1.2. Hàm lồi . . . . . . . . . . . . . . . . . . . . . . . . . . . 3 1.3. Bài toán cân bằng . . . . . . . . . . . . . . . . . . . . . 8 1.4. Ánh xạ giả co chặt và các tính chất . . . . . . . . . . . . 12 Chương 2. Phương pháp dưới đạo hàm cho bài toán cân bằng và ánh xạ giả co chặt 20 2.1. Cách tiếp cận . . . . . . . . . . . . . . . . . . . . . . . . 20 2.2. Thuật toán 2.1 . . . . . . . . . . . . . . . . . . . . . . . 21 2.3. Định lý hội tụ mạnh 2.1 . . . . . . . . . . . . . . . . . . 22 Chương 3. Ứng dụng của phương pháp dưới đạo hàm cho bài toán cân bằng 32 3.1. Thuật toán 3.1 . . . . . . . . . . . . . . . . . . . . . . . 32 3.2. Định lý hội tụ mạnh 3.1 . . . . . . . . . . . . . . . . . . 33 3.3. Ví dụ minh họa và các kết quả tính toán . . . . . . . . . 37 Kết luận . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 41 Tài liệu tham khảo . . . . . . . . . . . . . . . . . . . . . . . 42 2Số hóa bởi Trung tâm Học liệu – Đại học Thái Nguyên http://www.lrc-tnu.edu.vn ii Lời cảm ơn Luận văn này được hoàn thành tại trường Đại học Khoa học - Đại học Thái Nguyên. Tác giả xin bày tỏ lòng biết ơn sâu sắc đối với thầy PGS.TS Phạm Ngọc Anh (Học viện Công nghệ Bưu chính Viễn thông), thầy đã trực tiếp hướng dẫn tận tình và động viên tác giả trong suốt thời gian nghiên cứu vừa qua. Xin chân thành cảm ơn các quý thầy giáo, cô giáo đã trực tiếp giảng dạy lớp Cao học Toán K4C, các bạn học viên, và các bạn đồng nghiệp đã tạo điều kiện thuận lợi, động viên giúp đỡ tác giả trong quá trình học tập và nghiên cứu tại trường. Tác giả cũng xin bày tỏ lòng biết ơn sâu sắc tới gia đình và người thân luôn khuyến khích động viên tác giả trong suốt quá trình học cao học và viết luận văn này. Mặc dù có nhiều cố gắng nhưng luận văn khó tránh khỏi những thiếu sót và hạn chế. Tác giả mong nhận được những ý kiến đóng góp của quý thầy cô và bạn đọc để luận văn được hoàn thiện hơn. Xin chân thành cảm ơn! Thái Nguyên, tháng 10 năm 2012 Tác giả Phạm Thị Linh 3Số hóa bởi Trung tâm Học liệu – Đại học Thái Nguyên http://www.lrc-tnu.edu.vn iii Những ký hiệu và chữ viết tắt R : Tập hợp các số thực R + : Tập hợp các số thực không âm R n : Không gian số thực n - chiều R n + : Không gian số thực không âm n - chiều x ∈ C : x thuộc tập C x ∈ C : x không thuộc tập C ∀x : Với mọi x ∃x : Tồn tại x ∅ : Tập hợp rỗng ∩ : Phép giao các tập hợp ∪ : Phép hợp các tập hợp x := y : x được định nghĩa bằng y x, y : Tích vô hướng của x và y x k  x : Dãy  x k  hội tụ yếu tới x x k → x : Dãy  x k  hội tụ mạnh tới x I : Ánh xạ đồng nhất x : Chuẩn của véc tơ x [x, y] : Đoạn thẳng nối hai điểm x và y 4Số hóa bởi Trung tâm Học liệu – Đại học Thái Nguyên http://www.lrc-tnu.edu.vn 1 Lời nói đầu Bài toán cân bằng được mô tả dưới dạng một bất đẳng thức, gọi là bất đẳng thức Ky Fan, lần đầu được áp dụng để nghiên cứu các mô hình cân bằng kinh tế theo khái niệm cân bằng do J. Nash, nhà toán học Mỹ đoạt giải Nobel kinh tế trong những công trình nghiên cứu về cân bằng đưa ra vào năm 1994. Về mặt lý thuyết, như sự tồn tại nghiệm, nhiều kết quả quan trọng đã đạt được cho bài toán cân bằng tổng quát trên các không gian trừu tượng. Tuy nhiên, về mặt thuật toán và các ứng dụng, các kết quả còn hạn chế. Luận văn này trình bày một số thuật toán để giải bài toán tìm nghiệm chung của tập nghiệm của bài toán cân bằng và tập các điểm bất động của một họ hữu hạn các ánh xạ giả co chặt. Luận văn gồm mục lục, ba chương, phần kết luận và tài liệu tham khảo. Chương 1 sẽ nhắc lại các kiến thức cơ bản nhất của tập lồi và hàm lồi, mà các kết quả này sẽ được sử dụng ở các chương sau. Phần cuối của chương sẽ giới thiệu về bài toán cân bằng, một số ví dụ và cuối cùng sẽ trình bày về ánh xạ giả co chặt, phép chiếu trực giao với các tính chất. Chương 2 sẽ trình bày thuật toán để giải bài toán tìm nghiệm chung của tập nghiệm của bài toán cân bằng và tập các điểm bất động của một họ hữu hạn các ánh xạ giả co chặt với định lý hội tụ mạnh. Chương 3 là phần ứng dụng. Phần này trình bày về việc áp dụng thuật toán để giải một số bài toán cân bằng với các kết quả tính toán cụ thể. Đây cũng là những đóng góp mới được ứng dụng để giải bài toán cân bằng thông qua sự gắn kết giữa phương pháp dưới đạo hàm và các kỹ thuật điểm bất động. 5Số hóa bởi Trung tâm Học liệu – Đại học Thái Nguyên http://www.lrc-tnu.edu.vn 2 Chương 1 Một số Khái niệm Cơ bản Trong luận văn này, ta xét bài toán cân bằng và bài toán điểm bất động của một họ hữu hạn các ánh xạ giả co chặt trên không gian Hilbert thực H. Dưới đây, ta nhắc lại một số khái niệm và các tính chất cơ bản của giải tích lồi như: Tập lồi, hàm lồi, dưới vi phân, và một số kiến thức liên quan đến bài toán cân bằng, ánh xạ giả co chặt, phép chiếu trực giao cùng với các tính chất tương ứng. Các kiến thức trong chương này được lấy chủ yếu từ các tài liệu [1],[3],[5],[6]. 1.1. Tập lồi và các phép toán cơ bản Định nghĩa 1.1. Cho C là tập con, khác rỗng của không gian Hilbert thực H. Tập C được gọi là lồi (convex) nếu với mọi x, y ∈ C và λ ∈ [0, 1] , ta có λx + (1 − λ) y ∈ C. Đặc biệt, H và ∅ là các tập lồi. Ví dụ 1.1. Các nửa không gian là các tập lồi. Các tam giác và hình tròn trong mặt phẳng là các tập lồi. Hình cầu đơn vị trong không gian Banach là tập lồi Định nghĩa 1.2. 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 gọi là tập lồi đa diện (polyhedral convex set) hay khúc lồi. Định nghĩa 1.3. Tập con C khác rỗng trong không gian Hilbert thực H được gọi là nón (cone) nếu λx ∈ C, ∀ x ∈ C, ∀λ > 0. 6Số hóa bởi Trung tâm Học liệu – Đại học Thái Nguyên http://www.lrc-tnu.edu.vn 3 Tập con C khác rỗng trong không gian Hilbert thực H được gọi là nón lồi nếu nó vừa là nón vừa là lồi. Điều đó có nghĩa là λx + µy ∈ C, ∀ x, y ∈ C, ∀λ, µ > 0. Ví dụ 1.2. Tập R n + là nón lồi trong R n . Định nghĩa 1.4. Cho không gian Hilbert thực H, tập C ⊆ H lồi, khác rỗng và điểm x ∗ ∈ C. Khi đó, nón pháp tuyến ngoài của C tại x ∗ (hay còn gọi là nón lồi đóng), kí hiệu N C (x ∗ ) , được xác định bởi N C (x ∗ ) : = {p ∈ H : p, x − x ∗   0, ∀x ∈ C} . 1.2. Hàm lồi Định nghĩa 1.5. Cho không gian Hilbert thực H, tập C ⊆ H và hàm f : C → R ∪ {+∞, −∞}. Khi đó, các tập hợp dom f : = {x ∈ C : f (x) < +∞} , epif : = {(x, α) ∈ C × R : f (x)  α} , tương ứng, được gọi là miền hữu hiệu (effective domain) và trên đồ thị (epigraph) của f. Hàm f được gọi là chính thường (proper) trên C nếu dom f = ∅, f (x) > −∞, ∀ x ∈ C. Định nghĩa 1.6. Cho không gian Hilbert thực H và tập C ⊆ H. Hàm f : C → R ∪ {+∞, −∞} được gọi là lồi (convex) trên C nếu trên đồ thị của nó là tập con lồi của H × R. Hàm f được gọi là lõm (concave) nếu −f là lồi. Bổ đề 1.1. Cho không gian Hilbert thực H và tập C ⊆ H. Nếu hàm f : C → R ∪ {+∞, −∞} lồi trên C thì miền hữu hiệu của f là tập lồi. 7Số hóa bởi Trung tâm Học liệu – Đại học Thái Nguyên http://www.lrc-tnu.edu.vn 4 Mệnh đề 1.1. Cho không gian Hilbert thực H và tập C ⊆ H. Khi đó, hàm f : C → R ∪ {+∞} được gọi là lồi trên C nếu và chỉ nếu với mọi x, y ∈ C, λ ∈ [0, 1] , ta có f (λx + (1 − λ) y)  λf(x) + (1 − λ) f(y). Chứng minh. Giả sử f là hàm lồi, không mất tính tổng quát có thể coi λ ∈ (0, 1). Không thể xảy ra trường hợp f (x) < +∞, f (y) < +∞ mà f (λx + (1 − λ) y) = +∞, bởi vì dom f lồi. Hơn nữa, với mọi x, y ∈ dom f, thì [x, y] ⊂ dom f. Vì λ ∈ (0, 1) nên f (x) = +∞, suy ra λf (x) = +∞. Nếu x hoặc y /∈ dom f thì f (x) = +∞ hoặc f (y) = +∞. Mặt khác, vì epi f lồi nên với mọi (x, α) , (y, β) ∈ epif, λ ∈ (0, 1) , ta có λ (x, α) + (1 − λ) (y, β) = (λx + (1 − λ) y, λα + (1 − λ) β) ∈ epi f ⇒ f (λx + (1 − λ) y)  λα + (1 − λ) β ⇒ f (λx + (1 − λ) y)  λf (x) + (1 − λ) f (y) (lấy α = f (x) , β = f (y)). Ngược lại, giả sử với mọi x, y ∈ C, λ ∈ [0, 1] , ta có f (λx + (1 − λ) y)  λf(x) + (1 − λ) f(y). Lấy (x, α) , (y, β) ∈ epif, λ ∈ [0, 1] , ta có f (λx + (1 − λ) y)  λf (x) + (1 − λ) f (y)  λα + (1 − λ) β ⇒ (λx + (1 − λ) y, λα + (1 − λ) β) ∈ epi f ⇔ λ (x, α) + (1 − λ) (y, β) ∈ epi f. Ví dụ 1.3. Cho C là tập con lồi, khác rỗng của không gian Hilbert thực H. Khi đó, hàm chỉ (indicator funtion) trên C δ C (x) =  0 nếu x ∈ C; +∞ nếu x /∈ C; là một hàm lồi. 8Số hóa bởi Trung tâm Học liệu – Đại học Thái Nguyên http://www.lrc-tnu.edu.vn 5 Định nghĩa 1.7. Cho không gian Hilbert thực H và tập C ⊆ H. Hàm f : C → R ∪ {+∞} được gọi là lồi chặt (strict convex) trên C nếu với mọi x, y ∈ C, x = y, λ ∈ (0, 1) , ta có f (λx + (1 − λ) y) < λf(x) + (1 − λ) f(y). Ví dụ 1.4. Cho không gian Hilbert thực H. Khi đó, với mọi x ∈ H hàm x 2 = x, x là một hàm lồi chặt trên H. Định nghĩa 1.8. Cho không gian Hilbert thực H, tập C ⊆ H và hàm lồi f : C → R ∪ {+∞} . Khi đó, dưới vi phân (subdifferential) của f tại x ∗ , ký hiệu là ∂f (x ∗ ), được xác định bởi ∂f (x ∗ ) : = {p ∈ H : f(x) − f(x ∗ )  p, x − x ∗  , ∀x ∈ C} . Hàm f được gọi là khả dưới vi phân (subdifferentiable) trên C nếu ∂f (x) = ∅, ∀ x ∈ C. Ví dụ 1.5. (Dưới vi phân của hàm chỉ) Cho C là tập con lồi, khác rỗng của không gian Hilbert thực H. Xét hàm chỉ trên C δ C (x) =  0 nếu x ∈ C; +∞ nếu x /∈ C. Khi đó, ∂δ C (x ∗ ) = N C (x ∗ ) , ∀x ∗ ∈ C. Chứng minh. Nếu x ∗ ∈ C thì ∂δ C (x ∗ ) = 0 và ∂δ C (x ∗ ) = {p ∈ H : ∂δ C (x)  p, x − x ∗  , ∀x ∈ C} . = {p ∈ H : 0  p, x − x ∗  , ∀x ∈ C} = N C (x ∗ ) . 9Số hóa bởi Trung tâm Học liệu – Đại học Thái Nguyên http://www.lrc-tnu.edu.vn 6 Ví dụ 1.6. (Dưới vi phân của hàm lồi thuần nhất dương) Cho f : R n → R là hàm lồi thuần nhất dương, tức là hàm lồi f thỏa mãn f (λx) = λf (x) , ∀λ > 0, ∀x ∈ R n . Khi đó, ∂f (x ∗ ) = {p ∈ R n : p, x ∗  = f (x ∗ ) , p, x  f (x) , ∀x ∈ C} . Chứng minh. Nếu p ∈ ∂f (x ∗ ) thì p, x − x ∗   f (x) − f (x ∗ ) , ∀x ∈ C. (1.1) Thay x = 2x ∗ vào (1.1), ta có p, x ∗   f (x ∗ ) . (1.2) Thay x = 0 vào (1.1), ta nhận được − p, x ∗   −f (x ∗ ) . (1.3) Kết hợp (1.2) và (1.3), ta có p, x ∗  = f (x ∗ ) . Hơn nữa, p, x − x ∗  = p, x − p, x ∗  = p, x − f (x ∗ ) . Do đó p, x  f (x) , ∀x ∈ C. Ngược lại, nếu x ∗ ∈ R n thỏa mãn p, x ∗  = f (x ∗ ) và p, x  f (x) , ∀x ∈ C thì p, x − x ∗  = p, x − p, x ∗   f (x) − f (x ∗ ) , ∀x ∈ C. Vậy p ∈ ∂f (x ∗ ). 10Số hóa bởi Trung tâm Học liệu – Đại học Thái Nguyên http://www.lrc-tnu.edu.vn [...]... http://www.lrc-tnu.edu.vn 20 Chương 2 Phương pháp dưới đạo hàm cho bài toán cân bằng và ánh xạ giả co chặt 2.1 Cách tiếp cận Cho C là tập con lồi, đóng, khác rỗng của không gian Hilbert thực H Với mỗi i = 1, , p, các ánh xạ Si : C → C giả co chặt trên C với Li ∈ [0, 1) Khi đó, bài toán tìm điểm bất động chung của một họ hữu hạn các ánh xạ {Si }p (p i=1 1) giả co chặt trên C được mô tả như sau: Tìm x∗ ∈ C sao cho x∗ ∈ ∩p F... của phương pháp dưới đạo hàm cho bài toán cân bằng 3.1 Thuật toán 3.1 Trong trường hợp đặc biệt ánh xạ đồng nhất chính là ánh xạ giả co chặt trên một tập con lồi, đóng, khác rỗng của không gian Hilbert thực H cho trước với L bất kỳ trong nửa đoạn [0, 1) Vì vậy để minh họa cho thuật toán 2.1 trong chương này ta chọn họ hữu hạn các ánh xạ giả co chặt chính là họ các ánh xạ đồng nhất Khi đó, thuật toán. .. bước 3, ta thấy các dãy y k , z k và tk cũng hội tụ mạnh về điểm x∗ Kết luận chương Chương này đã trình bày các nội dung cơ bản về phương pháp dưới đạo hàm cho bài toán cân bằng và ánh xạ giả co chặt Trong đó có thuật toán và định lý hội tụ mạnh dùng đề tìm nghiệm chung của tập nghiệm của bài toán cân bằng và tập các điểm bất động của một họ hữu hạn các ánh xạ giả co chặt 35Số hóa bởi Trung tâm Học... động của ánh xạ Si (i = 1, 2, , p) Trong luận văn này, chúng ta quan tâm đến bài toán tìm nghiệm chung của tập nghiệm của bài toán cân bằng và tập các điểm bất động của một họ hữu hạn các ánh xạ giả co chặt trên C Điều này có nghĩa là Tìm x∗ ∈ ∩p F ix (Si ) ∩ Sol (f, C) i=1 (2.1) Nếu f (x, y) = F (x) , y − x thì bài toán (2.1) trở thành bài toán tìm một nghiệm chung của tập nghiệm của bài toán bất... Hilbert thực H và song hàm f : C × C → R sao cho f (x, x) = 0, ∀ x ∈ C Bài toán cân bằng, viết tắt là EP (f, C), được phát biểu như sau: Tìm x∗ ∈ C sao cho f (x∗ , y) 0, ∀ y ∈ C Ký hiệu tập nghiệm của bài toán EP (f, C) là Sol(f, C) Cho ánh xạ F : C → H, ta định nghĩa f (x, y) : = F (x) , y − x , ∀x, y ∈ C Khi đó, bài toán EP (f, C) trở thành bài toán bất đẳng thức biến phân: Tìm x∗ ∈ C sao cho F (x∗ )... lý 1.2 Cho C là tập con lồi, đóng, khác rỗng trong không gian Hilbert thực H và hàm f : C → R lồi, khả dưới vi phân trên C Khi đó, điều kiện cần và đủ để điểm x∗ ∈ C là cực tiểu của f trên C là 0 ∈ ∂f (x∗ ) + NC (x∗ ) , trong đó, ∂f (x∗ ) là ký hiệu dưới vi phân của hàm f tại x∗ còn NC (x∗ ) là ký hiệu của nón pháp tuyến ngoài của C tại x∗ 1.3 Bài toán cân bằng Phát biểu bài toán Cho C là tập con lồi,... của bài toán EP (f, C) Ví dụ 1.8 Bài toán bù phi tuyến Cho C ⊂ Rn là nón lồi đóng, C ∗ = {x ∈ Rn : x, y 0, ∀y ∈ C} là nón đối ngẫu của nón C và ánh xạ liên tục S : C → Rn Khi đó, bài toán bù phi tuyến được phát biểu như sau: Tìm x∗ ∈ C sao cho S (x∗ ) ∈ C ∗ và S (x∗ ) , x∗ = 0 Nếu đặt f (x, y) := S (x) , y − x , ∀x, y ∈ C thì bài toán bù phi tuyến tương đương với bài toán EP (f, C) Chứng minh Giả. .. λ ∈ (0, 1) và L ∈ [0, 1) nên V1 x − V2 x = 0 Suy ra, S1 x = S2 x = x Do đó (1 − λ) S1 x + λS2 x = x Vậy x ∈ F ix (S1 ) ∩ F ix (S2 ) Kết luận chương Trong chương này, ta đã nhắc lại các kết quả quan trọng của giải tích lồi như tập lồi, hàm lồi, dưới vi phân, bài toán cân bằng Đồng thời, ta cũng trình bày các khái niệm về ánh xạ đơn điệu, giả đơn điệu Đặt biệt là khái niệm ánh xạ giả co chặt, phép chiếu... Lipschitz trên C với hai hằng số c1 = c2 = 1.4 B−M 2 Ánh xạ giả co chặt và các tính chất Định nghĩa 1.11 Cho C là tập con lồi, đóng, khác rỗng của không gian Hilbert thực H Ánh xạ S : C → C được gọi là 16Số hóa bởi Trung tâm Học liệu – Đại học Thái Nguyên http://www.lrc-tnu.edu.vn 13 (i) giả co chặt (strict pseudocontraction) trên C nếu tồn tại L ∈ [0, 1) sao cho với mọi x, y ∈ C, ta có S (x) − S (y) 2 x−y... đẳng thức biến phân và tập nghiệm của bài toán điểm bất động Toàn bộ chương 2 của luận văn này sẽ trình bày về một sơ đồ lặp mới được các tác giả trong [2] đề cập để giải bài toán (2.1) với các giả thiết: Giả thiết 2.1 Song hàm f thoả mãn các điều kiện: (i) f giả đơn điệu và liên tục yếu trên C; (ii) f liên tục kiểu Lipschitz trên C; (iii) với mỗi x ∈ C, hàm f (x, ) lồi và khả dưới vi phân trên C 24Số . ĐẠI HỌC KHOA HỌC PHẠM THỊ LINH PHƯƠNG PHÁP DƯỚI ĐẠO HÀM CHO BÀI TOÁN CÂN BẰNG VÀ ÁNH XẠ GIẢ CO CHẶT Chuyên nghành: TOÁN ỨNG DỤNG Mã số: 60.46.01.12 LUẬN VĂN THẠC SĨ TOÁN HỌC Người hướng dẫn khoa. . . . 12 Chương 2. Phương pháp dưới đạo hàm cho bài toán cân bằng và ánh xạ giả co chặt 20 2.1. Cách tiếp cận . . . . . . . . . . . . . . . . . . . . . . . . 20 2.2. Thuật toán 2.1 . . . . toán để giải một số bài toán cân bằng với các kết quả tính toán cụ thể. Đây cũng là những đóng góp mới được ứng dụng để giải bài toán cân bằng thông qua sự gắn kết giữa phương pháp dưới đạo hàm