PHƯƠNG PHÁP CHIẾU GIẢI BÀI TOÁN BẤT ĐẲNG THỨC BIẾN PHÂN GIẢ ĐƠN ĐIỆU MẠNH

ĐẠI HỌC QUỐC GIA HÀ NỘI TRƯỜNG ĐẠI HỌC KHOA HỌC TỰ NHIÊN ————————— NGÔ THỊ THO PHƯƠNG PHÁP CHIẾU GIẢI BÀI TOÁN BẤT ĐẲNG THỨC BIẾN PHÂN GIẢ ĐƠN ĐIỆU MẠNH LUẬN VĂN THẠC SĨ TOÁN HỌC Hà Nội - 2015 ĐẠI HỌC QUỐC GIA HÀ NỘI TRƯỜNG ĐẠI HỌC KHOA HỌC TỰ NHIÊN ————————— NGÔ THỊ THO PHƯƠNG PHÁP CHIẾU GIẢI BÀI TOÁN BẤT ĐẲNG THỨC BIẾN PHÂN GIẢ ĐƠN ĐIỆU MẠNH Chuyên ngành: Toán ứng dụng Mã số: 60460112 LUẬN VĂN THẠC SĨ TOÁN HỌC NGƯỜI HƯỚNG DẪN KHOA HỌC: GS.TSKH LÊ DŨNG MƯU Hà Nội - 2015 Mục lục Lời cảm ơn Lời mở đầu Chương Bài toán bất đẳng thức biến phân 1.1 Kiến thức chuẩn bị 1.1.1 Hội tụ mạnh yếu không gian Hilbert 1.1.2 Toán tử chiếu 1.1.3 Tính liên tục hàm lồi 1.1.4 Đạo hàm vi phân hàm lồi 1.2 Bài toán bất đẳng thức biến phân 8 1.2.1 Các khái niệm 1.2.2 Sự tồn nghiệm 11 Chương Phương pháp chiếu giải toán bất đẳng thức biến phân giả đơn điệu mạnh 12 2.1 Phương pháp chiếu đạo hàm tăng cường 13 2.2 Phương pháp chiếu cải biên 16 Kết luận 22 Tài liệu tham khảo 23 LỜI CẢM ƠN Lời đầu tiên, em xin chân thành cảm ơn GS.TSKH Lê Dũng Mưu Thầy người hướng dẫn khóa luận tốt nghiệp hướng dẫn luận văn thạc sĩ cho em Hai chặng đường qua, thầy tận tình hướng dẫn bảo nghiêm khắc, thầy cung cấp nhiều tài liệu quan trọng giành nhiều thời gian giải đáp thắc mắc suốt trình làm việc thầy Em xin gửi tới thầy, cô Khoa Toán - Cơ - Tin học, Trường Đại học Khoa Học Tự Nhiên, Đại học Quốc Gia Hà Nội, thầy cô giảng dạy lớp Cao học Toán khóa 2013 - 2015, lời cảm ơn chân thành công lao dạy dỗ thầy, cô hai năm qua Đặc biệt, em muốn gửi lời cảm ơn tới thầy dạy chuyên ngành nhóm Toán Ứng Dụng Mặc dù nhóm có tám thành viên thầy lên lớp với nhiệt huyết chuyên đề hay, sâu sắc Cuối em xin gửi lời cảm ơn tới gia đình, bạn, anh, chị lớp cao học Toán khóa 2013 - 2015 giành riêng lời cảm ơn cho gia đình Toán Ứng Dụng Là em út nhóm, nên người quan tâm nhiều Thời gian học anh chị cho em kỷ niệm đẹp, học điều hay kiến thức thú vị Mặc dù có nhiều cố gắng, luận văn không tránh khỏi thiếu sót Em mong nhận ý kiến đóng góp thầy, cô bạn đọc để luận văn hoàn thiện Hà Nội, ngày tháng 10 năm 2015 Học viên Ngô Thị Tho LỜI MỞ ĐẦU Năm 1966, Hatman Stampacchia công bố nghiên cứu toán bất đẳng thức biên phân, liên quan tới việc giải toán biến phân, toán điều kiển tối ưu toán biên có dạng phương trình đạo hàm riêng Năm 1980, Kinderlehrer Stampacchia cho xuất sách "An Introduction to Variational Inequalities and Their Applications", giới thiệu toán biến phân không gian vô hạn chiều ứng dụng Năm 1984, sách "Variational and Quasivariational Inequalities: Applications to Free Boundary Problems" C Baiocci A Capelo áp dụng bất đẳng thức biến phân tựa biến phân để giải toán biên Hiện toán bất đẳng thức biến phân phát triển thành nhiều dạng khác nhau,như là: bất đẳng thức biến phân vectơ, tựa bất đẳng thức biến phân, giả bất đẳng thức biến phân, bất đẳng thức biến phân ẩn, bất đẳng thức biến phân suy rộng Bài toán thu hút quan tâm nhiều nhà toán học Vì mô hình chứa nhiều toán quan trọng số lĩnh vực toán học thực tế tối ưu hóa, toán bù, lý thuyết trò chơi, cân Nash, cân mạng giao thông, cân di trú Một hướng nghiên cứu quan trọng bất đẳng thức biến phân việc xây dựng phương pháp giải Dựa tính chất kiểu đơn điệu G Cohen nghiên cứu phương pháp nguyên lý toán phụ Ngoài có phương pháp hiệu chỉnh Tikhonov, phương pháp chiếu, phương pháp điểm Những phương pháp hiệu quả, dễ thực máy tính hội tụ chúng đảm bảo sở giả thiết khác tính chất đơn điệu Có nhiều phương pháp chiếu khác nhau, là: phương pháp chiếu bản, phương pháp chiếu đạo hàm, phương pháp chiếu siêu phẳng Mỗi phương pháp giải lớp toán bất đẳng thức biến phân định Do hội tụ thuật toán đảm bảo Luận văn trình bày phương pháp chiếu đạo hàm tăng cường chiếu cải biên để giải toán bất đẳng thức biến phân giả đơn điệu mạnh Các phương pháp tạo dãy hội tụ điểm lặp dễ dàng tính Chúng hội tụ tới nghiệm toán Luận văn gồm hai chương: Chương 1: Bài toán bất đẳng thức biến phân, chia làm hai phần: • Phần 1: Nhắc lại số kiến thức Giải tích hàm Giải tích lồi, là: hội tụ mạnh yếu không gian Hilbert, toán tử chiếu, tính liên tục hàm lồi, đạo hàm vi phân hàm lồi • Phần 2: Phát biểu toán, trình bày số khái niệm mô hình minh họa cho toán Sau đó, chứng minh tồn tính nghiệm toán Chương 2: Phương pháp chiếu giải toán bất đẳng thức biến phân giả đơn điệu mạnh Nội dung chương trình bày hai thuật toán chiếu đạo hàm tăng cường thuật toán chiếu cải biên để giải toán V I(K, F) Phát biểu chứng minh định lý hội tụ dãy lặp tạo thuật toán Đưa số ví dụ chứng minh điều kiện định lý tồn nghiệm cần thiết Nếu bỏ điều kiện đó, dãy lặp không hội tụ tới nghiệm toán Chương Bài toán bất đẳng thức biến phân Trong chương này, nhắc lại số kết Giải tích hàm có liên quan tới hội tụ mạnh hội tụ yếu dãy số Nhắc lại số khái niệm định lý Giải tích lồi, là: định nghĩa tính chất toán tử chiếu, tính liên tục, đạo hàm vi phân hàm lồi, Định lý tách, Định lý MoreauRockafellar Phần sau ta giới thiệu toán bất đẳng thức biến phân (VIP) nhấn mạnh toán bất đẳng thức biến phân giả đơn điệu mạnh Chỉ ví dụ toán bất đẳng thức biến phân thường gặp thực tế mô hình toán học Cuối chương phát biểu chứng minh định lý tồn tính nghiệm toán Nội dung chủ yếu trích dẫn từ tài liệu [1], [2], [3], [6], [10] Trong luận văn này, làm việc không gian Hilbert thực trang bị tô pô yếu, với tích vô hướng , chuẩn tương ứng ||.|| 1.1 Kiến thức chuẩn bị 1.1.1 Hội tụ mạnh yếu không gian Hilbert Định nghĩa 1.1.1 Giả sử H không gian tuyến tính thực, với x ∈ H xác định số gọi chuẩn x ( kí hiệu ||x||) thỏa mãn ba tiên đề sau: Xác định dương: ∀x ∈ H ||x|| ≥ 0; Thuần dương: ∀x ∈ H; ∀λ ∈ R Bất đẳng thức tam giác: ∀x, y ∈ H ||x|| = ⇔ x = ||λ x|| = |λ | ||x|| ||x + y|| ≤ ||x|| + ||y|| Định nghĩa 1.1.2 Giả sử H không gian tuyến tính thực, cặp (H, , ) với , : H ×H → R (x, y) → x, y thỏa mãn điều kiện: Xác định dương: x, x ≥ ∀x ∈ H; x, x = ⇔ x = Đối xứng: x, y = y, x ∀x, y ∈ H Song tuyến tính: αx + β y, z = α x, z + β y, z ∀α, β ∈ R, ∀x, y, z ∈ H gọi không gian tiền Hilbert Không gian tiền Hilbert, đầy đủ gọi không gian Hilbert, kí hiệu H Định nghĩa 1.1.3 1) Ta nói dãy {xk } hội tụ mạnh đến x ( kí hiệu xk → x) lim ||xk − x|| = k→∞ 2) Dãy {xk } hội tụ yếu đến x ( kí hiệu xk tức ∀ f ∈ H∗ x) {xk } hội tụ x theo tô pô yếu σ f (xk ) → f (x) Khi H không gian hữu hạn chiều tô pô yếu tô pô thông thường H trùng Đặc biệt, dãy hội tụ mạnh hội tụ yếu 1.1.2 Toán tử chiếu Định nghĩa 1.1.4 Giả sử C tập lồi, khác rỗng không gian Hilbert thực H x0 ∈ C Nón pháp tuyến (ngoài) C x0 kí hiệu NC (x0 ) định nghĩa bởi: NC (x0 ) := {ω ∈ H| ω T (x − x0 ) ≤ ∀x ∈ C} Tập −NC (x0 ) gọi nón pháp tuyến (trong) C x0 Nón pháp tuyến ε C x0 định nghĩa bởi: NCε (x0 ) := {ω ∈ H| ω T (x − x0 ) ≤ ε ∀x ∈ C} Định nghĩa 1.1.5 Giả sử C = 0/ (không thiết lồi) tập không gian Hilbert H y véc-tơ bất kỳ, khoảng cách từ y đến C định nghĩa dC (y) := inf ||x − y|| x∈C Nếu tồn π ∈ C cho dC (y) := ||π − y||, ta nói π hình chiếu (khoảng cách) y C, kí hiệu π = pC (y) Mệnh đề 1.1.1 Cho C tập lồi đóng khác rỗng Khi đó: Với y ∈ H, π ∈ C hai tính chất sau tương đương: a) π = pC (y), b) y − π ∈ NC (π) Với y ∈ H, hình chiếu pC (y) y C tồn Nếu y ∈ / C, pC (y) − y, x − pC (y) = siêu phẳng tựa C pC (y) tách hẳn y khỏi C, tức pC (y) − y, x − pC (y) ≥ ∀x ∈ C, pC (y) − y, y − pC (y) < Ánh xạ y → pC (y) có tính chất sau: a) ||pC (x) − pC (y)|| ≤ ||x − y|| ∀x, ∀y (tính không giãn), b) pC (x) − pC (y), x − y ≥ ||pC (x) − pC (y)||2 , (tính đồng bức) Toán tử chiếu công cụ hữu hiệu nhằm giải toán cân trường hợp đặc biệt như: Bài toán tối ưu, bất đẳng thức biến phân, điểm bất động, toán điểm yên ngựa Trong luận văn này, ta vận dụng giải toán bất đẳng thức biến phân giả đơn điệu mạnh 1.1.3 Tính liên tục hàm lồi Cho C ⊆ H tập lồi f : C → R ∪ {+∞}, ta kí hiệu: dom f := {x ∈ C : f (x) < +∞} Tập dom f gọi miền hữu dụng tập f Tập epi f := {(x, µ) ∈ C × R : f (x) ≤ µ}, gọi đồ thị hàm f Hàm f gọi thường dom f = 0/ f (x) > −∞ với x ∈ dom f Định nghĩa 1.1.6 Hàm f gọi lồi epi f tập lồi Hàm f hàm lõm − f hàm lồi Nếu f vừa lồi vừa lõm ta nói f hàm afin Định nghĩa 1.1.7 Cho f : H → R, hàm f gọi nửa liên tục x0 ∈ H ∀{xk } ⊂ H : xk → x0 ⇒ limk→∞ f (xk ) ≥ f (x0 ) Hàm f gọi nửa liên tục D ⊆ H liên tục x ∈ D Hàm f nửa liên tục − f nửa liên tục Nếu hàm f vừa liên tục vừa liên tục liên tục 1.1.4 Đạo hàm vi phân hàm lồi Tính khả vi hàm lồi đóng vai trò quan trọng phương pháp tối ưu hóa Lớp hàm lồi có tính chất đẹp mà lớp hàm khác Giả sử f : H → R hàm lồi Ta có khái niệm sau Định nghĩa 1.1.8 Vectơ ω ∈ H ∗ gọi đạo hàm f x0 ∈ H nếu: ω, x − x0 ≤ f (x) − f (x0 ), ∀x ∈ H Tập hợp tất đạo hàm hàm f x0 gọi vi phân hàm f x0 , kí hiệu ∂ f (x0 ) := {ω ∈ H ∗ : ω, x − x0 ≤ f (x) − f (x0 ), ∀x ∈ H} Hàm f gọi khả vi phân x0 ∂ f (x0 ) = / Định nghĩa 1.1.9 Cho ε > 0, vectơ ω ∈ H ∗ gọi ε-dưới đạo hàm f x0 ∈ H nếu: ω, x − x0 ≤ f (x) − f (x0 ) + ε, ∀x ∈ H Tập hợp tất ε-dưới đạo hàm hàm f x0 gọi ε-dưới vi phân hàm f x0 , kí hiệu ∂ε f (x0 ) := {ω ∈ H ∗ : ω, x − x0 ≤ f (x) − f (x0 ) + ε, ∀x ∈ H} Hàm f gọi ε-khả vi phân x0 ∂ε f (x0 ) = / Định nghĩa 1.1.10 Giả sử x ∈ H, d ∈ H\{0}, hàm f gọi là: a) Khả vi Frechet x0 tồn ω ∈ H ∗ cho f (x) − f (x0 ) − ω, x − x0 = 0, ∀x ∈ H lim ||x − x0 || x→x0 Một điểm ω tồn tại, gọi đạo hàm f x0 , kí hiệu f (x0 ) ∇ f (x0 ) b) Có đạo hàm theo hướng d x0 tồn giới hạn lim t→0+ f (x0 + td) − f (x0 ) t ta gọi giới hạn đạo hàm theo hướng d f x0 , kí hiệu f (x0 , d) Nói chung hàm lồi không thiết khả vi điểm Dưới vi phân khái niệm mở rộng đạo hàm trường hợp hàm không khả vi Trong trường hợp ∂ f (x∗ ) gồm điểm f khả vi x∗ 1.2 Bài toán bất đẳng thức biến phân 1.2.1 Các khái niệm Định nghĩa 1.2.1 Cho K ⊂ H tập đóng, khác rỗng, F : K → H ánh xạ đơn trị Bài toán bất đẳng thức biên phân (đơn trị) toán Tìm x∗ ∈ K cho F(x∗ ), y − x∗ ≥ 0, ∀y ∈ K Tập nghiệm toán kí hiệu S(K, F) (VIP) Định nghĩa 1.2.2 Giả sử K ⊂ H tập lồi đóng, khác rỗng toán tử F : K → H gọi a) đơn điệu mạnh K tồn γ > cho F(x) − F(y), x − y ≥ γ||x − y||2 ∀x, y ∈ K, b) đơn điệu K F(x) − F(y), x − y ≥ ∀x, y ∈ K, c) giả đơn điệu mạnh K tồn γ > cho F(x), y − x ≥ ⇒ F(y), y − x ≥ γ||y − x||2 ∀x, y ∈ K, d) giả đơn điệu K F(x), y − x ≥ ⇒ F(y), y − x ≥ ∀x, y ∈ K Theo định nghĩa kéo theo (a) ⇒ (b), (a) ⇒ (c), (c) ⇒ (d), (b) ⇒ (d), hiển nhiên Chú ý toán tử giả đơn điệu mạnh không đơn điệu Định nghĩa 1.2.3 Giả sử toán tử F toán tử giả đơn điệu mạnh toán bất đẳng thức biến phân (VIP) trở thành toán bất đẳng thức biến phân giả đơn điệu mạnh Kí hiệu VI(K, F) Mệnh đề 1.2.1 Cho K ⊂ H tập lồi đóng toán tử F : K → H liên tục a) Nếu F giả đơn điệu mạnh VI(K, F) có có nghiệm b) Nếu F giả đơn điệu S(K, F) tập lồi Định nghĩa 1.2.4 Ánh xạ F gọi liên tục Lipschitz K tồn số L > cho: ||F(u) − F(v)|| ≤ L||u − v||, ∀u, v ∈ K PK (.) ánh xạ liên tục Lipschitz với số Lipschitz L=1 Ví dụ 1.2.1 Bài toán điểm bất động Brouwer 10 Cho K tập lồi, đóng, khác rỗng, compăc yếu H ánh xạ T : K → K ánh xạ liên tục Bài toán điểm bất động phát biểu sau Tìm x∗ ∈ K : x∗ = T (x∗ ) Bài toán điểm bất động đưa toán bất đẳng thức biến phân (VIP) thông qua mệnh đề sau Mệnh đề 1.2.2 Giả sử ánh xạ F xác định F(x) = x − T (x), ∀x ∈ K Khi đó, nghiệm toán điểm bất động trùng với nghiệm toán bất đẳng thức biến phân (VIP) Tức toán bất đẳng thức biến phân tương đương với toán điểm bất động 1.2.2 Sự tồn nghiệm Trong phần này, trình bày chứng minh theo tài liệu [9] toán bất đẳng thức biến phân giả đơn điệu mạnh tồn nghiệm nghiệm nhấtcủa bái toán Bổ đề 1.2.1 Giả sử K tập lồi đóng, khác rỗng không gian Hilbert H Cho F : K → H toán tử cho F(∗), y − ∗ nửa liên tục với y ∈ K Giả sử ∃ tập compăc W : ∀x ∈ K\W, ∃y ∈ K : F(x), y − x < Thì toán bất đẳng thức biến phân (VI) có nghiệm Mệnh đề 1.2.3 Giả sử F β − giả đơn điệu mạnh K Nếu F(∗), y − ∗ nửa liên tục với y ∈ K toán bất đẳng thức biến phân (VI) có nghiệm 11 Chương Phương pháp chiếu giải toán bất đẳng thức biến phân giả đơn điệu mạnh Chương này, trình bày thuật toán chiếu giải toán bất đẳng thức biến phân giả đơn điệu mạnh V I(K, F) Phần đầu, trình bày thuật toán chiếu đạo hàm tăng cường (mỗi bước lặp có hai lần chiếu) với độ dài bước lấy từ khoảng đóng số thực dương, không yêu cầu phụ thuộc vào số Lipschitz Dãy lặp tổng quát thu từ thuật toán hội tụ tới nghiệm toán Phần sau, trình bày thuật toán chiếu cải biên (mỗi bước lặp có lần chiếu) với độ dài bước lấy tùy ý từ khoảng đóng cố định số thực dương Trình bày chứng minh dãy lặp tổng quát thu từ thuật toán hội tụ tuyến tính tới nghiệm toán Phương pháp đòi hỏi số Lipschitz môđun giả đơn điệu mạnh việc lựa chọn khoảng đóng cố định chứa độ dài bước Nội dung chủ yếu chương trích dẫn từ tài liệu [7], [8] 12 Đề tiện cho việc theo dõi, em xin phát biểu lại toán bất đẳng thức biến phân giả đơn điệu manh sau: Định nghĩa 2.0.5 Cho K ⊂ H tập đóng, khác rỗng, F : K → H toán tử giả đơn điệu mạnh K Bài toán bất đẳng thức biến phân giả đơn điệu mạnh toán Tìm x∗ ∈ K cho F(x∗ ), y − x∗ ≥ 0, ∀y ∈ K (VI) Tập nghiệm toán ký hiệu VI(K, F) 2.1 Phương pháp chiếu đạo hàm tăng cường Thuật toán 2.1.1 Chọn điểm đầu u0 ∈ K dãy độ dài bước {λk }∞ k=0 ⊂ R+ với ∞ ∑ λk = +∞, k=0 lim λk = k→∞ Bước 1: Đặt k=0 Bước 2: Tính uk = PK (uk − λk F(uk )), uk+1 = PK (uk − λk F(uk )) Bước 3: Nếu uk = uk dừng lại Ngược lại đặt k + = k quay lại bước Nếu thuật toán dừng bước thứ k, ta đặt uk = uk với k ≥ k + Vì thế, Thuật toán 2.1.1 tạo dãy lặp vô hạn Định lý 2.1.1 Cho K tập lồi, đóng, khác rỗng không gian Hilbert thực H Giả sử F : K → H liên tục Lipschitz giả đơn điệu mạnh K, dãy lặp {uk } tạo Thuật toán 2.1.1 hội tụ mạnh tới nghiệm u∗ toán Hơn nữa, tồn số k0 ∈ N cho γλk < với k ≥ k0 ||u k+1 k ∗ ∏ (1 − γλ j )||uk0 − u∗ ||, − u || ≤ j=k0 γ > số giả đơn điệu mạnh F Ngoài ra, k lim k→∞ ∏ (1 − γλ j ) = j=k0 13 Trong phần tiếp theo, ta chứng minh giả thiết hàm F giả đơn điệu mạnh ∞ hai điều kiện ∑ λk = +∞, lim λk = cần thiết cho khẳng định Định lý 2.1.1 k→∞ k=0 Ví dụ 2.1.1 Cho K = R F(u) = u Dễ thấy, F liên tục Lipschitz, đơn điệu mạnh K S(K, F) = {0} Chọn u0 = ∈ K định nghĩa dãy {λk } cách đặt λk = (k + 1)2 ∀k ∈ N ∞ Từ ∑ λk < +∞, dãy lặp {uk } tạo Thuật toán 2.1.1 cho k=0 uk+1 = PK (uk − λk F(PK (uk − λk F(uk )))) = uk − λk F(uk − λk F(uk )) = uk − λk (uk − λk uk ) = (1 − λk + λk2 )uk Vì u0 = 1, ta có k k+1 u = ∏ j=1 = = ∏ 1− ∏ k ≥ k (1 − λ j + λ j2 ) j=0 j=0 1− ( j + 1)2 k+2 2(k + 1) 1 + ( j + 1) ( j + 1)4 k j( j + 2) j=1 ( j + 1) =∏ ∀k ∈ N Do đó, {uk } dãy giảm bị chặn nên hội tụ k+2 = k→∞ 2(k + 1) lim uk ≥ lim k→∞ Ta lim uk = u∗ với u∗ ≥ Do vậy, dãy {uk } không hội tụ tới nghiệm k→∞ toán VI(K, F) Ví dụ 2.1.2 Cho K, F, u0 Ví dụ 2.1.1 λk = với k ∈ N Do k−1 ∞ ∑ λk = +∞ Tuy nhiên tính toán uk = ∏ (1 − λ j + λ j2 ) = với j=0 k=0 k ∈ N Do lim k→∞ uk = không hội tụ tới nghiệm toán Ví dụ 2.1.3 Cho K = R2 F(u) = (−u2 , u1 ) với u = (u1 , u2 ) ∈ K Dễ thấy hàm F liên tục Lipschitz đơn điệu K, nghiệm toán VI(K, F) (0, 0)T 14 Cho u0 = (u01 , u02 )T ∈ K\{(0, 0)T } {λk } cho λk = k+1 ∀k ∈ N ∞ Khi ∑ λk = +∞, lim λk = Để chứng minh F không giả đơn điệu mạnh K, k→∞ ta cần chọn u = (1, 0)T , v = (2, 0)T , k=0 Dãy lặp {uk } F(u), v − u = F(v), v − u = tạo Thuật toán 2.1.1 cho    u0 = (u01 , u02 )T     k+1 1 k u1 = 1− uk1 + u (k + 1) k+1    1 k   = 1− uk2 − u uk+1 2 (k + 1) k+1 Ta có ||uk+1 || k+1 2 = (uk+1 ) + (u2 ) = (uk1 )2 + (uk2 )2 =||uk || − =||u || ∏ 1 + (k + 1)2 (k + 1)4 1 + (k + 1)2 (k + 1)4 k 1− 1− j=0 1 + ( j + 1) ( j + 1)4 Khi k k lim ||u || = ||u || lim k→∞ k→∞ ≥ ||u0 || lim k→∞ ∏ 1− j=0 1 + ( j + 1)2 ( j + 1)4 k+2 2(k + 1) √ ≥ ||u || √ = µ||u0 || (µ ≥ ) Do đó, {uk } không hội tụ đến nghiệm toán VI(K, F) Qua ví dụ ta thấy rằng, bỏ ba kiện dãy lặp {uk } không hội tụ tới nghiệm toán VI(F, K) 15 2.2 Phương pháp chiếu cải biên Thuật toán 2.2.1 Cho trước u0 ∈ K {λk } ⊂ (0, +∞) Bước 1: Đặt k=0 Bước 2: Tính u k = PK (uk − λk F(uk )), u k = uk dừng lại Nếu không chuyển sang bước Bước 3: Đặt uk+1 = u k , sau quay lại bước Nếu thuật toán chấm dứt bước thứ k, ta đặt uk = uk với k ≥ k + Như vậy, dãy độ dài bước thay đổi {λk } ⊂ (0, +∞), Thuật toán 2.2.1 tạo cho điểm ban đầu u0 ∈ K dãy lặp vô hạn {uk } Mệnh đề 2.2.1 Cho K tập lồi, đóng, khác rỗng không gian Hilbert thực H ánh xạ F : K → H giả đơn điệu mạnh K với mô-đun γ liên tục Lipschitz K với số L Cho {uk } dãy lặp tạo Thuật toán 2.2.1 u∗ nghiệm toán VI(K, F) [1 + λk (2γ − λk L2 )] ||uk+1 − u∗ ||2 ≤ ||uk − u∗ ||2 , ∀k ∈ N Bài toán VI(K, F) có nghiệm, dãy lặp tạo Thuật toán 2.2.1, với độ dài bước chọn từ khoảng đóng số thực dương, hội tụ tuyến tính tới nghiệm Cụ thể, ta có định lý sau: Định lý 2.2.1 Cho K tập lồi, đóng, khác rỗng không gian Hilbert thực H, F : K → H ánh xạ giả đơn điệu mạnh K với môđun γ liên tục Lipschitz K với số L Giả sử < a ≤ λk ≤ b < 2γ , ∀k ∈ N, L2 (2.4) a, b số dương Cho {uk } dãy lặp tạo Thuật toán 2.2.1 Khi đó, dãy {uk } hội tụ tuyến tính tới nghiệm u∗ toán Hơn nữa, sai số tiên nghiệm hậu nghiệm ||uk+1 − u∗ || ≤ µ k+1 ||u − u0 ||, 1−µ ||uk+1 − u∗ || ≤ µ ||uk+1 − uk ||, 1−µ 16 với k ∈ N Ở µ= 1 + a(2γ − bL2 ) ∈ (0, 1) Chú ý 2.2.1 Khi a = b = λ , độ dài bước cố định Do phương pháp chiếu cải biên trở thành phương pháp chiếu µ trở thành µ= 1 + λ (2γ − λ L2 ) 2γ Các ước tính sai số chặt chẽ µ nhỏ Xét L2 2γ µ hàm λ ∈ 0, , ta tìm giá trị nhỏ µ L L γ µ ∗ := điểm λ∗ := L L2 + γ Ta có λ ∈ 0, Chú ý 2.2.2 Ngoài giá trị µ xem hàm µ = µ(a, b) biến (a, b) thuộc miền (a, b) ∈ R2 : < a ≤ b < 2γ L2 Đặt b = ta, với t ∈ [1, + ∞) cố định, tương tự ta tính hàm γ a = Từ µ(a, b) = µ(a,ta) đạt giá trị nhỏ tL γ2 1+ tL         L : ≤ t < +∞ = 2 + γ2   L γ    1+  tL Vậy suy µ(a, b) : < a ≤ b < Do đó, giá trị nhỏ µ µ ∗ = (a∗ , b∗ ) = ( γ γ , ) L2 L2 17 2γ L2 L L2 + γ = L L2 + γ đạt điểm Chú ý 2.2.3 Ước lượng sai số Định lý 2.2.1 hữu ích việc áp dụng Thuật toán 2.2.1 để giải toán bất đẳng thức biến phân giả đơn điệu mạnh.Ví dụ, µ k+1 ||u −u0 || cho phép ước tính số lượng công thức ||uk+1 −u∗ || ≤ 1−µ bước lặp cần để đạt độ xác định Cụ thể, với ε > bất µ k+1 kỳ, ||u − u0 || ≤ ε ta có ||uk+1 − u∗ || ≤ ε 1−µ Hệ 2.2.1 Trong kí hiệu Định lý 2.2.1, F đơn điệu mạnh liên tục Lipschitz K dãy {uk } tạo Thuật toán 2.2.1 hội tụ tuyến tính tới nghiệm toán VI(K,F) ước tính sai số nêu thỏa mãn Ví dụ 2.2.1 Cho H = l2 không gian Hilbert thực mà thành phần dãy bình phương khả tổng vô hướng thực Ví dụ H = {u = (u1 , u2 , · · · , ui , · · · ) : ∞ ∞ i=1 i=1 ∑ |ui |2 < +∞} Định nghĩa u, v = ∑ ui vi ||u|| = u, u tích vô hướng chuẩn hai vectơ u, v H, với u = (u1 , u2 , · · · , ui , · · · ),v = (v1 , v2 , · · · , vi , · · · ) ∈ H β Cho α, β ∈ R cho β > α > > Đặt Kα = {u ∈ H : ||u|| ≤ α}, Fβ (u) = (β − ||u||)u, α, β tham số Dễ thấy S(Kα , Fβ ) = {0} Hàm Fβ liên tục Lipschitz giả đơn điệu mạnh Kα Thật vậy, với u, v ∈ Kα bất kỳ, ||Fβ (u) − Fβ (v)|| = ||(β − ||u||)u − (β − ||v||)v|| = ||β (u − v) − ||u||(u − v) − (||u|| − ||v||)v|| ≤ β ||u − v|| + ||u|| ||u − v|| + | ||u|| − ||v|| | ||v|| ≤ β ||u − v|| + α ||u − v|| + α||u − v|| = (β + 2α)||u − v|| Do Fβ liên tục Lipschitz Kα với số Lipschitz L := β + 2α Cho u, v ∈ Kα cho Fβ (u), v − u ≥ Theo giả thiết ||u|| ≤ α < β Suy u, v − u ≥ Do Fβ (v), v − u = (β − ||v||) v, v − u ≥ (β − ||v||) v, v − u − u, v − u ≥ (β − α)||u − v||2 = γ||u − v||2 , 18 γ := β − α > Suy Fβ giả đơn điệu mạnh Kα Hơn Fβ β đơn điệu mạnh đơn điệu Kα Thật vậy, ta chọn u = , 0, · · · , 0, · · · , v = (α, 0, · · · , 0, · · · ) ∈ Kα Ta có Fβ (u) − Fβ (v), u − v = β −α < 2γ 2(β − α) = 0, tùy ý, đặt λk = λ với L2 (β + 2α)2 k ∈ N Theo Định lý 2.2.1, dãy {uk } tạo Thuật toán 2.2.1 hội tụ tuyến tính tới Hơn nữa, Lấy u0 ∈ Kα bất kỳ, λ ∈ 0, ||uk+1 − 0|| ≤ µ µ k+1 ||u − u0 || ||uk+1 − 0|| ≤ ||uk+1 − uk || 1−µ 1−µ với k ∈ N, µ= 1 + λ [2(β − α) − λ (β + 2α)2 ] Theo Chú ý 2.2.1 giá trị nhỏ µ µ∗ = λ = λ∗ = β + 2α (β − α)2 + (β + 2α)2 ] điểm β −α (β + 2α)2 Nếu độ dài bước tạo thành dãy không khả tổng số thực dương Thuật toán 2.2.1 tạo dãy lặp hội mạnh tới nghiệm toán Ta có: Định lý 2.2.2 Cho K tập lồi đóng khác rỗng không gian Hilbert thực H F : K → H ánh xạ giả đơn điệu mạnh K với mô-đun γ liên tục Lipschitz K với số L Giả sử {λk } dãy vô hướng dương với ∞ ∑ λk = +∞, k=0 lim λk = (2.5) k→∞ Dãy lặp {uk } tạo thành từ Thuật toán 2.2.1 hội tụ mạnh tới u∗ nghiệm toán VI(K, F) Hơn nữa, tồn số k0 ∈ N cho với k ≥ k0 , λk (2γ − λk L2 ) > 0, ||uk+1 − u∗ || ≤ ∏ki=k0 [1 + λi (2γ − λi L2 ] 19 ||uk0 − u∗ || Hệ 2.2.2 Cho {λk } Định lý 2.2.2 Cho F đơn điệu mạnh K với mô-đun γ liên tục Lipschitz K với số L Thì dãy {uk } tạo từ Thuật toán 2.2.1 hội tụ theo chuẩn tới nghiệm toán VI(K, F) tồn số k0 ∈ N cho ||uk+1 − u∗ || ≤ ||uk0 − u∗ ||, ∏ki=k0 [1 + λi (2γ − λi L2 ] với k ≥ k0 Ví dụ 2.2.2 Cho H = l2 , α, β ∈ R cho β > α > Kα = {u ∈ H : ||u|| ≤ α}, β > Đặt Fβ (u) = (β − ||u||)u, ∞ , ∀k ∈ N Khi ∑ λk = +∞, lim λk = k→∞ k+1 k=0 k Cho {u } dãy lặp tạo Thuật toán 2.2.1 Suy dãy {uk } hội tụ (β + 2α)2 mạnh tới nghiệm toán VI(Kα , Fβ ) Đặt k0 = , 2(β − α) λk (2γ − λk L2 ) > 0, với k ≥ k0 Ta có α, β tham số, λk = ||uk+1 − 0|| ≤ ∏ki=k0 + ||uk0 − 0||, + 2α)2 ∀k ≥ k0 2(β − α) (β − i+1 (i + 1)2 Chúng ta xét xem điều xảy bỏ điều kiện (2.4) (2.5) đưa Định lý 2.2.1 Định lý 2.2.2 thông qua ví dụ sau: Ví dụ 2.2.3 Đặt K = R F(u) = u Rõ ràng F liên tục Lipschitz, đơn điệu mạnh K S(K, F) = Chọn u0 = ∈ K λk = , (k + 2)2 ∀k ∈ N ∞ Từ lim λk = ∑ λk < +∞, hai điều kiện (2.4) (2.5) bị lược bỏ Dãy k→∞ lặp uk k=0 tạo Thuật toán 2.2.1 với u0 = cho uk+1 = PK (uk − λk F(uk )) = uk − λk uk = (1 − λk )uk 20 Do đó, k k uk+1 = ∏(1 − λi ) = ∏ 1− i=0 i=0 (i + 2)2 k k+3 (i + 1)(i + 3) = , (i + 2)2 2(k + 2) i=0 =∏ ∀k ∈ N Nghĩa {uk } không hội tụ đến nghiệm toán k→∞ VI(K, F) Vậy điều kiện (2.4) (2.5) bỏ đi, không dãy lặp không hội tụ tới nghiệm toán cần tìm Ví dụ 2.2.4 Cho K = R F(u) = u, u0 ∈ R\{0} Cho λk ⊂ (0, +∞) thỏa mãn Vậy lim uk = ∞ điều kiện ∑ λk = +∞, lim λk = λk = 1, ∀k ∈ N Khi dãy lặp {uk } trở thành k=0 k→∞ uk+1 = PK (uk − λk F(uk )) = uk − λk uk = (1 − λk )uk Để chứng minh {uk } hội tụ tuyến tính đến 0, ta cần chứng minh: ||uk+1 − 0|| lim = µ, với µ ∈ (0, 1) k→∞ ||uk − 0|| Vì lim λk = uk = với k ∈ N, ta có k→∞ ||uk+1 − 0|| = lim |1 − λk | = k→∞ ||uk − 0|| k→∞ lim Vậy {uk } không hội tụ tuyến tính tới nghiệm toán VI(K, F) Ví dụ cho thấy dãy {uk } xét Định lý 2.2.2 không hội tụ tuyến tính tới nghiệm toán VI(K, F) Mặt khác, so sánh với dãy lặp tạo Định lý 2.2.1, công thức lặp Định lý 2.2.2 có tốc độ hội tụ chậm Như vậy, bên cạnh ưu điểm nêu trên, phương pháp chiếu cải biên số tiên nghiệm có nhược điểm tốc độ hội tụ 21 KẾT LUẬN Sau Hatman Stampacchia giới thiệu lần đầu vào năm 1966, trải qua 50 năm phát triển không ngừng, toán bất đẳng thức biến phân trở thành công cụ hữu hiệu, để nghiên cứu giải toán cân kinh tế tài chính, vận tải, lý thuyết trò chơi nhiều toán khác Gần đây, toán bất đẳng thức biến phân nhiều nhà toán học quan tâm, có nhiều kết quan trọng Người ta tìm nhiều phương pháp để giải toán Bản luận văn nhằm mục đích giới thiệu toán bất đẳng thức biến phân giả đơn điệu mạnh Cụ thể là, sau tổng hợp lại số kiến thức Giải tích hàm Giải tích lồi, trình bày toán bất đẳng thức biến phân với ví dụ minh họa Sau đó, luận văn trình bày tồn nghiệm toán Tiếp đến, luận văn giới thiệu hai thuật toán chiếu để giải lớp toán này, đồng thời xét đến hội tụ thuật toán không gian Hilbert TÀI LIỆU THAM KHẢO Tiếng Việt Hoàng Tụy (2005), Hàm thực giải tích hàm, Nhà xuất Đại học Quốc gia Hà Nội Lê Dũng Mưu, Nguyễn Văn Hiền, Nguyễn Hữu Điển (2015), Giải tích lồi ứng dụng, Nhà xuất Đại học Quốc gia Hà Nội Phạm Kỳ Anh, Trần Đức Long (2001), Hàm thực giải tích hàm, Nhà xuất Đại học Quốc gia Hà Nội Tiếng Anh D Kinderlehrer and G Stampacchia (1980), An Introduction to Variational Inequalities and Their Applications, Academic Press, New York Fan Ky (1972), A minimax inequalities and applications In: Shisha O (Ed): Inequalities, Academic Press, New York Igor Konnov (2001), Combined Relaxation Methods for Variational Inequalities, Springer Pham Duy Khanh (2012), ”A new extragradient method for strongly pseudomonotone variational inequalities”, Submitted Pham Duy Khanh, Phan Tu Vuong (2014), ”Modified projection method for strongly pseudomonotone variational inequalities”, Journal of Global Optimization,58, no 2, 341 - 350 Phung M Duc, Le D Muu, and Nguyen V Quy (2014), ”Solution - existence and algorithms with their convergence rate for strongly pseudomonotone equilibrium problems”, Pracific Journal Mathematics, Pacific J Mathematics, To appear 23 [...]... Thì bài toán bất đẳng thức biến phân (VI) có một nghiệm Mệnh đề 1.2.3 Giả sử F là β − giả đơn điệu mạnh trên K Nếu F(∗), y − ∗ là nửa liên tục trên với mỗi y ∈ K thì bài toán bất đẳng thức biến phân (VI) có một nghiệm duy nhất 11 Chương 2 Phương pháp chiếu giải bài toán bất đẳng thức biến phân giả đơn điệu mạnh Chương này, trình bày các thuật toán chiếu giải bài toán bất đẳng thức biến phân giả đơn điệu. .. K, d) giả đơn điệu trên K nếu F(x), y − x ≥ 0 ⇒ F(y), y − x ≥ 0 ∀x, y ∈ K Theo định nghĩa trên các kéo theo (a) ⇒ (b), (a) ⇒ (c), (c) ⇒ (d), (b) ⇒ (d), là hiển nhiên Chú ý rằng một toán tử giả đơn điệu mạnh có thể không đơn điệu Định nghĩa 1.2.3 Giả sử toán tử F là toán tử giả đơn điệu mạnh thì bài toán bất đẳng thức biến phân (VIP) trở thành bài toán bất đẳng thức biến phân giả đơn điệu mạnh Kí hiệu... tục Bài toán điểm bất động được phát biểu như sau Tìm x∗ ∈ K : x∗ = T (x∗ ) Bài toán điểm bất động được đưa về bài toán bất đẳng thức biến phân (VIP) thông qua mệnh đề sau Mệnh đề 1.2.2 Giả sử ánh xạ F được xác định bởi F(x) = x − T (x), ∀x ∈ K Khi đó, nghiệm của bài toán điểm bất động trùng với nghiệm của bài toán bất đẳng thức biến phân (VIP) Tức là bài toán bất đẳng thức biến phân tương đương với bài. .. nhằm mục đích giới thiệu bài toán bất đẳng thức biến phân giả đơn điệu mạnh Cụ thể là, sau khi tổng hợp lại một số kiến thức cơ bản về Giải tích hàm và Giải tích lồi, trình bày về bài toán bất đẳng thức biến phân với các ví dụ minh họa Sau đó, luận văn trình bày về sự tồn tại và duy nhất nghiệm của bài toán Tiếp đến, luận văn giới thiệu hai thuật toán chiếu để giải lớp bài toán này, đồng thời xét đến... triển không ngừng, bài toán bất đẳng thức biến phân đã trở thành một công cụ hữu hiệu, để nghiên cứu và giải các bài toán cân bằng trong kinh tế tài chính, vận tải, lý thuyết trò chơi và trong nhiều bài toán khác Gần đây, các bài toán bất đẳng thức biến phân đã được nhiều nhà toán học quan tâm, và có được nhiều kết quả quan trọng Người ta đã tìm ra nhiều phương pháp để giải bài toán này Bản luận văn... như sau: Định nghĩa 2.0.5 Cho K ⊂ H là một tập đóng, khác rỗng, F : K → H là toán tử giả đơn điệu mạnh trên K Bài toán bất đẳng thức biến phân giả đơn điệu mạnh là bài toán Tìm x∗ ∈ K sao cho F(x∗ ), y − x∗ ≥ 0, ∀y ∈ K (VI) Tập nghiệm của bài toán được ký hiệu là VI(K, F) 2.1 Phương pháp chiếu dưới đạo hàm tăng cường Thuật toán 2.1.1 Chọn một điểm đầu u0 ∈ K và một dãy các độ dài bước {λk }∞ k=0 ⊂ R+... tổng quát thu được từ thuật toán này hội tụ tuyến tính tới nghiệm duy nhất của bài toán Phương pháp này đòi hỏi hằng số Lipschitz và môđun của giả đơn điệu mạnh trong việc lựa chọn khoảng đóng cố định chứa độ dài bước Nội dung chủ yếu của chương được trích dẫn từ tài liệu [7], [8] 12 Đề tiện cho việc theo dõi, em xin phát biểu lại bài toán bất đẳng thức biến phân giả đơn điệu manh như sau: Định nghĩa... việc áp dụng Thuật toán 2.2.1 để giải bài toán bất đẳng thức biến phân giả đơn điệu mạnh. Ví dụ, µ k+1 1 ||u −u0 || cho phép chúng ta ước tính số lượng công thức ||uk+1 −u∗ || ≤ 1−µ bước lặp cần để đạt được một độ chính xác nhất định Cụ thể, với ε > 0 bất µ k+1 1 kỳ, nếu ||u − u0 || ≤ ε thì ta có ||uk+1 − u∗ || ≤ ε 1−µ Hệ quả 2.2.1 Trong các kí hiệu của Định lý 2.2.1, nếu F là đơn điệu mạnh và liên tục... khả vi tại mọi điểm Dưới vi phân là một khái niệm mở rộng của đạo hàm trong trường hợp hàm không khả vi Trong trường hợp ∂ f (x∗ ) chỉ gồm duy nhất một điểm thì f khả vi tại x∗ 1.2 Bài toán bất đẳng thức biến phân 1.2.1 Các khái niệm Định nghĩa 1.2.1 Cho K ⊂ H là một tập đóng, khác rỗng, F : K → H là một ánh xạ đơn trị Bài toán bất đẳng thức biên phân (đơn trị) là bài toán Tìm x∗ ∈ K sao cho F(x∗... đương với bài toán điểm bất động 1.2.2 Sự tồn tại nghiệm Trong phần này, trình bày chứng minh theo tài liệu [9] rằng bài toán bất đẳng thức biến phân giả đơn điệu mạnh luôn tồn tại một nghiệm và đó là nghiệm duy nhấtcủa bái toán Bổ đề 1.2.1 Giả sử K là một tập lồi đóng, khác rỗng trong không gian Hilbert H Cho F : K → H là một toán tử sao cho F(∗), y − ∗ là nửa liên tục trên với mỗi y ∈ K Giả sử ∃ tập