(LUẬN văn THẠC sĩ) một số thuật toán chiếu điểm gần kề giải phương trình với toán tử đơn điệu

ĐẠI HỌC QUỐC GIA HÀ NỘI TRƯỜNG ĐẠI HỌC KHOA HỌC TỰ NHIÊN NGUYỄN THỊ HUỆ MỘT SỐ THUẬT TOÁN CHIẾU - ĐIỂM GẦN KỀ GIẢI PHƯƠNG TRÌNH TỐN TỬ VỚI TOÁN TỬ ĐƠN ĐIỆU LUẬN VĂN THẠC SĨ KHOA HỌC Hà Nội - 2011 TIEU LUAN MOI download : skknchat@gmail.com ĐẠI HỌC QUỐC GIA HÀ NỘI TRƯỜNG ĐẠI HỌC KHOA HỌC TỰ NHIÊN NGUYỄN THỊ HUỆ MỘT SỐ THUẬT TỐN CHIẾU - ĐIỂM GẦN KỀ GIẢI PHƯƠNG TRÌNH VỚI TỐN TỬ ĐƠN ĐIỆU Chun ngành: Tốn học tính tốn Mã số: 604630 LUẬN VĂN THẠC SĨ KHOA HỌC Giáo viên hướng dẫn: GS TSKH Phạm Kỳ Anh Hà Nội - 2011 TIEU LUAN MOI download : skknchat@gmail.com Mục lục Lời cảm ơn i Bảng ký hiệu iv Lời nói đầu v Kiến thức chuẩn bị 1.1 Một số khái niệm 1.2 Phép chiếu không gian Hilbert Phương pháp chiếu-điểm gần kề 2.1 Giới thiệu 2.2 Sự hội tụ mạnh phương pháp chiếu-điểm gần kề 10 2.2.1 Thuật toán chiếu-điểm gần kề 10 2.2.2 Định lý hội tụ 13 2.3 Phương pháp chiếu-điểm gần kề song song 18 Phương pháp CQ 23 3.1 Các bổ đề quan trọng 23 3.2 Một số thuật toán CQ không gian Banach 25 3.3 Một số thuật toán CQ không gian Hilbert 37 Áp dụng 41 4.1 Bài tốn khơi phục ảnh 41 4.1.1 Phương pháp chiếu-điểm gần kề song song 42 4.1.2 Thử nghiệm số 44 4.2 Giải hệ phương trình đại số tuyến tính 45 4.2.1 Phương pháp chiếu-điểm gần kề song song 46 ii TIEU LUAN MOI download : skknchat@gmail.com 4.2.2 Phương pháp CQ song song 48 Kết luận 49 Tài liệu tham khảo 51 TIEU LUAN MOI download : skknchat@gmail.com (LUAN.van.THAC.si).mot.so.thuat.toan.chieu.diem.gan.ke.giai.phuong.trinh.voi.toan.tu.don.dieu(LUAN.van.THAC.si).mot.so.thuat.toan.chieu.diem.gan.ke.giai.phuong.trinh.voi.toan.tu.don.dieu(LUAN.van.THAC.si).mot.so.thuat.toan.chieu.diem.gan.ke.giai.phuong.trinh.voi.toan.tu.don.dieu(LUAN.van.THAC.si).mot.so.thuat.toan.chieu.diem.gan.ke.giai.phuong.trinh.voi.toan.tu.don.dieu Bảng ký hiệu p φ(·, ·) ∂f E∗ F (T ) Fb(T ) Khoảng cách suy rộng E Gradient f Không gian đối ngẫu E Tập điểm bất động T Tập điểm bất động tiệm cận T I Ánh xạ đồng J PC (·) ri(D) h·, ·i ⇀ S T[k]  Ánh xạ đối ngẫu chuẩn tắc Hình chiếu lên tập C Tập điểm tương đối D Tích đối ngẫu tích vơ hướng Hội tụ yếu Tập không điểm hay tập nghiệm Tk(modN ) Kết thúc chứng minh TIEU LUAN MOI download : skknchat@gmail.com (LUAN.van.THAC.si).mot.so.thuat.toan.chieu.diem.gan.ke.giai.phuong.trinh.voi.toan.tu.don.dieu(LUAN.van.THAC.si).mot.so.thuat.toan.chieu.diem.gan.ke.giai.phuong.trinh.voi.toan.tu.don.dieu(LUAN.van.THAC.si).mot.so.thuat.toan.chieu.diem.gan.ke.giai.phuong.trinh.voi.toan.tu.don.dieu(LUAN.van.THAC.si).mot.so.thuat.toan.chieu.diem.gan.ke.giai.phuong.trinh.voi.toan.tu.don.dieu (LUAN.van.THAC.si).mot.so.thuat.toan.chieu.diem.gan.ke.giai.phuong.trinh.voi.toan.tu.don.dieu(LUAN.van.THAC.si).mot.so.thuat.toan.chieu.diem.gan.ke.giai.phuong.trinh.voi.toan.tu.don.dieu(LUAN.van.THAC.si).mot.so.thuat.toan.chieu.diem.gan.ke.giai.phuong.trinh.voi.toan.tu.don.dieu(LUAN.van.THAC.si).mot.so.thuat.toan.chieu.diem.gan.ke.giai.phuong.trinh.voi.toan.tu.don.dieu Lời nói đầu Nhiều vấn đề khoa học kỹ thuật tốn chấp nhận lồi có ứng dụng lý thuyết tối ưu, khôi phục ảnh, phương pháp xử lý xạ, , đưa tốn tìm nghiệm hệ phương trình với tốn tử đơn điệu tìm điểm bất động họ hữu hạn ánh xạ không giãn (tương đối) Thuật tốn điểm gần kề tìm khơng điểm tốn tử đơn điệu cực đại Rockafellar [8] đề xuất vào năm 1976 trải qua nhiều lần cải biên Tuy nhiên thuật tốn nói chung cho kết hội tụ yếu Năm 2000, M V Solodov B F Svaiter [11] kết hợp thuật toán điểm gần kề với phép chiếu đơn giản lên giao nửa không gian để thu kết hội tụ mạnh Gần đây, P K Anh C V Chung [3] thực song song hóa thuật tốn chiếu-điểm gần kề để tìm nghiệm hệ phương trình toán tử đơn điệu Cũng dựa ý tưởng lai ghép, Nakajo Takahashi thu định lý hội tụ mạnh cho ánh xạ không giãn tương đối không gian Hilbert S Matsushita [7] tổng quát kết cho không gian Banach Năm 2011, Liu [6] cải biên thuật toán CQ Qin Su [9] để thu định lý hội tụ mạnh cho thuật tốn lặp xoay vịng tìm điểm bất động chung họ hữu hạn ánh xạ không giãn tương đối không gian Banach Đến năm 2011, P K Anh C V Chung đề xuất phương pháp CQ song song thuật tốn tốt thuật tốn lặp xoay vịng Liu chạy chế độ Bản luận văn tập trung trình bày thuật toán lai ghép với phép chiếu để giải hệ phương trình với tốn tử đơn điệu khơng gian Hilbert tìm điểm bất động họ ánh xạ khơng giãn tương đối khơng gian Banach Ngồi phần Mở đầu, kết luận tài liệu tham khảo, luận văn gồm chương: Chương 1: "Kiến thức chuẩn bị" trình bày số khái niệm tốn tử đơn TIEU LUAN MOI download : skknchat@gmail.com (LUAN.van.THAC.si).mot.so.thuat.toan.chieu.diem.gan.ke.giai.phuong.trinh.voi.toan.tu.don.dieu(LUAN.van.THAC.si).mot.so.thuat.toan.chieu.diem.gan.ke.giai.phuong.trinh.voi.toan.tu.don.dieu(LUAN.van.THAC.si).mot.so.thuat.toan.chieu.diem.gan.ke.giai.phuong.trinh.voi.toan.tu.don.dieu(LUAN.van.THAC.si).mot.so.thuat.toan.chieu.diem.gan.ke.giai.phuong.trinh.voi.toan.tu.don.dieu (LUAN.van.THAC.si).mot.so.thuat.toan.chieu.diem.gan.ke.giai.phuong.trinh.voi.toan.tu.don.dieu(LUAN.van.THAC.si).mot.so.thuat.toan.chieu.diem.gan.ke.giai.phuong.trinh.voi.toan.tu.don.dieu(LUAN.van.THAC.si).mot.so.thuat.toan.chieu.diem.gan.ke.giai.phuong.trinh.voi.toan.tu.don.dieu(LUAN.van.THAC.si).mot.so.thuat.toan.chieu.diem.gan.ke.giai.phuong.trinh.voi.toan.tu.don.dieu vi Lời nói đầu điệu, hình học khơng gian Banach vài tính chất quan trọng phép chiếu dùng luận văn Chương 2: "Phương pháp chiếu-điểm gần kề" trình bày hai thuật tốn lai ghép phương pháp chiếu phương pháp điểm gần kề khơng gian Hilbert, phần Chương trình bày thuật tốn chiếu-điểm gần kề tìm khơng điểm toán tử đơn điệu cực đại (đơn trị đa trị) định lý hội tụ mạnh; phần Chương trình bày thuật tốn chiếu-điểm gần kề song song tìm nghiệm hệ phương trình tốn tử đơn điệu định lý hội tụ mạnh Chương 3: "Phương pháp CQ" trình bày định lý hội tụ mạnh phương pháp CQ cho tốn tử khơng giãn tương đối, phương pháp CQ xoay vòng phương pháp CQ song song cho họ tốn tử khơng giãn tương đối khơng gian Banach Phần Chương trình bày ứng dụng phương pháp CQ cải biên CQ song song không gian Hilbert Chương 4: "Áp dụng", áp dụng phương pháp song song luận văn để giải tốn khơi phục ảnh khơng gian Hilbert giải hệ phương trình đại số tuyến tính Trong chương chúng tơi đưa ví dụ số thực mơi trường Matlab minh họa hình học cho phương pháp chiếu-điểm gần kề song song Hà nội, ngày tháng 12 năm 2011 TIEU LUAN MOI download : skknchat@gmail.com (LUAN.van.THAC.si).mot.so.thuat.toan.chieu.diem.gan.ke.giai.phuong.trinh.voi.toan.tu.don.dieu(LUAN.van.THAC.si).mot.so.thuat.toan.chieu.diem.gan.ke.giai.phuong.trinh.voi.toan.tu.don.dieu(LUAN.van.THAC.si).mot.so.thuat.toan.chieu.diem.gan.ke.giai.phuong.trinh.voi.toan.tu.don.dieu(LUAN.van.THAC.si).mot.so.thuat.toan.chieu.diem.gan.ke.giai.phuong.trinh.voi.toan.tu.don.dieu (LUAN.van.THAC.si).mot.so.thuat.toan.chieu.diem.gan.ke.giai.phuong.trinh.voi.toan.tu.don.dieu(LUAN.van.THAC.si).mot.so.thuat.toan.chieu.diem.gan.ke.giai.phuong.trinh.voi.toan.tu.don.dieu(LUAN.van.THAC.si).mot.so.thuat.toan.chieu.diem.gan.ke.giai.phuong.trinh.voi.toan.tu.don.dieu(LUAN.van.THAC.si).mot.so.thuat.toan.chieu.diem.gan.ke.giai.phuong.trinh.voi.toan.tu.don.dieu Chương Kiến thức chuẩn bị 1.1 Một số khái niệm Cho H không gian Hilbert thực, E không gian Banach E ∗ đối ngẫu E Tập lồi: Tập C ⊂ H (hoặc E) gọi lồi C chứa đoạn thẳng nối hai điểm Tức C lồi ∀x, y ∈ C, ∀α ∈ [0, 1] ⇒ αx + (1 − α)y ∈ C Toán tử đơn điệu: Ánh xạ T : E → E ∗ gọi đơn điệu hT x − T y, x − yi ≥ ∀x, y ∈ E, ký hiệu hf, xi tích đối ngẫu Trường hợp E = E ∗ = H ta có tích vơ hướng H Bổ đề 1.1.1 Giả sử S tập không điểm T , tức S = {x ∈ H | T x = θ} Nếu T toán tử đơn điệu S 6= ∅ S tập lồi Chứng minh Với x1 , x2 ∈ S t ∈ (0, 1) đặt x = tx1 + (1 − t)x2 Vì T tốn tử đơn điệu nên ta có ≤ hT x − T x1 , x − x1 i = hT x, tx1 + (1 − t)x2 − x1 i = (1 − t)hT x, x2 − x1 i TIEU LUAN MOI download : skknchat@gmail.com (LUAN.van.THAC.si).mot.so.thuat.toan.chieu.diem.gan.ke.giai.phuong.trinh.voi.toan.tu.don.dieu(LUAN.van.THAC.si).mot.so.thuat.toan.chieu.diem.gan.ke.giai.phuong.trinh.voi.toan.tu.don.dieu(LUAN.van.THAC.si).mot.so.thuat.toan.chieu.diem.gan.ke.giai.phuong.trinh.voi.toan.tu.don.dieu(LUAN.van.THAC.si).mot.so.thuat.toan.chieu.diem.gan.ke.giai.phuong.trinh.voi.toan.tu.don.dieu (LUAN.van.THAC.si).mot.so.thuat.toan.chieu.diem.gan.ke.giai.phuong.trinh.voi.toan.tu.don.dieu(LUAN.van.THAC.si).mot.so.thuat.toan.chieu.diem.gan.ke.giai.phuong.trinh.voi.toan.tu.don.dieu(LUAN.van.THAC.si).mot.so.thuat.toan.chieu.diem.gan.ke.giai.phuong.trinh.voi.toan.tu.don.dieu(LUAN.van.THAC.si).mot.so.thuat.toan.chieu.diem.gan.ke.giai.phuong.trinh.voi.toan.tu.don.dieu Chương Kiến thức chuẩn bị ≤ hT x − T x2 , x − x2 i = hT x, tx1 + (1 − t)x2 − x2 i = thT x, x1 − x2 i Từ điều suy hT x, x1 − x2 i = Vậy T x = θ nghĩa x ∈ S  Toán tử ngược-đơn điệu mạnh: T gọi đồng với số c > hay c-ngược đơn điệu mạnh H hT x − T y, x − yi ≥ ckT x − T yk2 ∀x, y ∈ H Toán tử đơn điệu cực đại: Toán tử T gọi đơn điệu cực đại đơn điệu đồ thị khơng phải tập thực đồ thị toán tử đơn điệu khác Toán tử J : E → E ∗ xác định J(x) = {f ∈ E ∗ | hf, xi = kxk2E = kf k2E ∗ } gọi ánh xạ đối ngẫu chuẩn tắc Hàm số φ : E × E → R cho φ(y, x) = kyk2 − 2hy, Jxi + kxk2 với x, y ∈ E, J ánh xạ đối ngẫu chuẩn tắc từ E vào E ∗ Đại p lượng φ(x, y) gọi khoảng cách suy rộng E Ánh xạ không giãn ánh xạ không giãn tương đối: Cho C tập lồi đóng E T ánh xạ từ C vào Ký hiệu F (T ) tập điểm bất động T Một điểm p C gọi điểm bất động tiệm cận T C chứa dãy {xn } hội tụ yếu tới p cho lim (xn − T xn ) = Ký n→∞ hiệu tập tất điểm bất động tiệm cận T Fb(T ) Ánh xạ T gọi không giãn kT x − T yk ≤ kx − yk ∀x, y ∈ C T ánh xạ không giãn tương đối Fb(T ) = F (T ) φ(p, T x) ≤ φ(p, x) với x ∈ C p ∈ F (T ) Nếu T tốn tử khơng giãn khơng gian Hilbert A = I − T tốn tử đơn điệu, I tốn tử đồng Thật vậy, với x, y ∈ H ta có kT x − T yk ≤ kx − yk TIEU LUAN MOI download : skknchat@gmail.com (LUAN.van.THAC.si).mot.so.thuat.toan.chieu.diem.gan.ke.giai.phuong.trinh.voi.toan.tu.don.dieu(LUAN.van.THAC.si).mot.so.thuat.toan.chieu.diem.gan.ke.giai.phuong.trinh.voi.toan.tu.don.dieu(LUAN.van.THAC.si).mot.so.thuat.toan.chieu.diem.gan.ke.giai.phuong.trinh.voi.toan.tu.don.dieu(LUAN.van.THAC.si).mot.so.thuat.toan.chieu.diem.gan.ke.giai.phuong.trinh.voi.toan.tu.don.dieu (LUAN.van.THAC.si).mot.so.thuat.toan.chieu.diem.gan.ke.giai.phuong.trinh.voi.toan.tu.don.dieu(LUAN.van.THAC.si).mot.so.thuat.toan.chieu.diem.gan.ke.giai.phuong.trinh.voi.toan.tu.don.dieu(LUAN.van.THAC.si).mot.so.thuat.toan.chieu.diem.gan.ke.giai.phuong.trinh.voi.toan.tu.don.dieu(LUAN.van.THAC.si).mot.so.thuat.toan.chieu.diem.gan.ke.giai.phuong.trinh.voi.toan.tu.don.dieu Chương Kiến thức chuẩn bị Do ≤ kAx − Ayk2 = kx − yk2 − 2hx − y, T x − T yi + kT x − T yk2
≤ kx − yk2 − hx − y, T x − T yi = 2hAx − Ay, x − yi Không gian Banach lồi đều: Không gian Banach E gọi lồi với ε > tồn δ(ε) > cho với x, y ∈ E, kxk ≤ 1, kyk ≤ kx − yk ≥ ε ta có kx + yk < 2(1 − δ(ε)) Chú ý Có thể thay đổi sau: với x, y ∈ E, d > tồn δ( dε ), kxk ≤ d, kyk ≤ d, kx − yk ≥ ε ta có kx + yk ≤ 2d(1 − δ( dε )) Chú ý Mọi không gian Hilbert không gian lồi Thật vậy, giả sử ε > 0, kxk ≤ 1, kyk ≤ cho kx − yk ≥ ε Từ bất đẳng thức hình bình hành suy Do
kx + yk2 = kxk2 + kyk2 − kx − yk2 ≤ − ε2 √ kx + yk ≤ − ε2 = − 1− r − ε2 !! Chú ý Mọi không gian Banach lồi không gian phản xạ, tức E ∗∗ = E Không gian Banach lồi chặt: Không gian Banach E gọi lồi chặt với x, y ∈ E, x 6= y, kxk ≤ 1, kyk ≤ ta có kx + yk < Không gian Banach trơn trơn đều: Đặt U = {x ∈ E | kxk = 1} hình cầu đơn vị E Khi không gian Banach E gọi trơn giới hạn kx + tyk − kxk t→0 t lim tồn với x, y ∈ U E gọi trơn giới hạn với x, y ∈ U Nếu E trơn J liên tục chuẩn-chuẩn tập bị chặn E Tính chất Kadec-Klee: Khơng gian Banach E có tính chất Kadec-Klee xn ⇀ x kxn k → kxk xn → x n → ∞ Bổ đề 1.1.2 Mọi không gian Banach lồi đều có tính chất Kadec-Klee TIEU LUAN MOI download : skknchat@gmail.com (LUAN.van.THAC.si).mot.so.thuat.toan.chieu.diem.gan.ke.giai.phuong.trinh.voi.toan.tu.don.dieu(LUAN.van.THAC.si).mot.so.thuat.toan.chieu.diem.gan.ke.giai.phuong.trinh.voi.toan.tu.don.dieu(LUAN.van.THAC.si).mot.so.thuat.toan.chieu.diem.gan.ke.giai.phuong.trinh.voi.toan.tu.don.dieu(LUAN.van.THAC.si).mot.so.thuat.toan.chieu.diem.gan.ke.giai.phuong.trinh.voi.toan.tu.don.dieu (LUAN.van.THAC.si).mot.so.thuat.toan.chieu.diem.gan.ke.giai.phuong.trinh.voi.toan.tu.don.dieu(LUAN.van.THAC.si).mot.so.thuat.toan.chieu.diem.gan.ke.giai.phuong.trinh.voi.toan.tu.don.dieu(LUAN.van.THAC.si).mot.so.thuat.toan.chieu.diem.gan.ke.giai.phuong.trinh.voi.toan.tu.don.dieu(LUAN.van.THAC.si).mot.so.thuat.toan.chieu.diem.gan.ke.giai.phuong.trinh.voi.toan.tu.don.dieu 38 Chương Phương pháp CQ Vì T khơng giãn nên với x, y ∈ H ta có kT x − T yk ≤ kx − yk tương đương với φ(T x, T y) ≤ φ(x, y) Do T ánh xạ không giãn tương đôi theo Định lý 3.2.1 ta thu xn → PF (T )(x) n → ∞  Cho Ti (i = 1, 2, , N ) họ ánh xạ không giãn từ C và N T F (Ti ) khác rỗng Chúng ta xét cải tiến thuật toán CQ giả sử F = i=1 song song khơng gian Hilbert H mà tập Cn Qn nửa không gian H Khi xn+1 = PCn ∩Qn (x0 ) tính tường minh Thuật tốn 3.3.1 ([2]) Cho x0 ∈ C giá trị ban đầu dãy {αn } ⊂ (0, 1) Với n = 0, 1, 2, • Đặt zk = PC xk • Tính yki = αk x0 + (1 − αk )Ti zk , • Tìm i = 1, 2, , N ik := argmax {kyki − xk k} i=1,2, ,N • Xác định Ck = {z ∈ H | kz − ykik k2 ≤ αk kz − x0 k2 + (1 − αk )kz − xk k2 }, Qk = {z ∈ H | hx0 − xk , xk − zi ≥ 0} • Tính xk+1 = PCk ∩Qk (x0 ) Trong Thuật tốn cải tiến để Ck Qk nửa khơng gian trùng với H nên xk+1 = PCk ∩Qk (x0 ) không thuộc C TIEU LUAN MOI download : skknchat@gmail.com (LUAN.van.THAC.si).mot.so.thuat.toan.chieu.diem.gan.ke.giai.phuong.trinh.voi.toan.tu.don.dieu(LUAN.van.THAC.si).mot.so.thuat.toan.chieu.diem.gan.ke.giai.phuong.trinh.voi.toan.tu.don.dieu(LUAN.van.THAC.si).mot.so.thuat.toan.chieu.diem.gan.ke.giai.phuong.trinh.voi.toan.tu.don.dieu(LUAN.van.THAC.si).mot.so.thuat.toan.chieu.diem.gan.ke.giai.phuong.trinh.voi.toan.tu.don.dieu (LUAN.van.THAC.si).mot.so.thuat.toan.chieu.diem.gan.ke.giai.phuong.trinh.voi.toan.tu.don.dieu(LUAN.van.THAC.si).mot.so.thuat.toan.chieu.diem.gan.ke.giai.phuong.trinh.voi.toan.tu.don.dieu(LUAN.van.THAC.si).mot.so.thuat.toan.chieu.diem.gan.ke.giai.phuong.trinh.voi.toan.tu.don.dieu(LUAN.van.THAC.si).mot.so.thuat.toan.chieu.diem.gan.ke.giai.phuong.trinh.voi.toan.tu.don.dieu 39 Chương Phương pháp CQ Do để bước tính yki có nghĩa bước đặt zk = PC (xk ) cần thiết Trong trường hợp C = H zk = xk Vì Ck Qk nửa không gian nên Ck ∩ Qk khác rỗng có cơng thức hiển tính PCk ∩Qk (x0 ) tương tự Thuật toán 2.2.1 Định lý sau cho kết hội tụ dãy {xk } Thuật toán 3.3.1 Định lý 3.3.2 Giả sử {xk } dãy sinh Thuật toán 3.3.1, Ti (i = 1, 2, , N ) liên tục lim αk = Khi xk hội tụ mạnh đến PF (x0 ) k → ∞ k→∞ Chứng minh.1 Chứng minh tương tự Định lý 3.2.4 với ý zk ∈ C ta có F ⊂ Ck ∩ Qk với k ≥ Và F 6= ∅ nên dãy {xk } hoàn toàn xác định Sử dụng xk = PCk (x0 ), xk+1 = PCk ∩Qk (x0 ) ∈ Ck F ⊂ Ck ∩ Qk với k ≥ 0, tương tự Định lý 3.2.4 ta thu kết sau kx0 − xk k ≤ kx0 − PF (x0 )k (3.11) ∀k; kxk − xk+1 k → k → ∞; kxk − yki k → k → ∞ với i = 1, 2, , N (3.12) Từ (3.11) suy dãy {xk } {yki } bị chặn với i = 1, 2, , N Hơn ta có kzk − yki k = kPC (xk ) − PC (yki )k ≤ kxk − yki k Nên sử dụng (3.12) ta có lim kzk − yki k = k→∞ lim kzk − xk k = (vì kzk − xk k ≤ kzk − yki k + kxk − yki k) k→∞ (3.13) Từ suy {zk } bị chặn Từ Thuật toán suy kTi zk − yki k = kTi zk − (αk x0 + (1 − αk )Ti (zk ))k = αk kx0 − Ti zk k, i = 1, 2, , N Với p ∈ F cố định, ta có kx0 − Ti zk k ≤ kx0 − pk + kp − Ti zk k ≤ kx0 − pk + kp − zk k với i = 1, 2, , N Bất đẳng thức cuối suy {kx0 − Ti zk k} bị chặn Và lim αk = nên ta có k→∞ lim kTi zk − yki k = k→∞ TIEU LUAN MOI download : skknchat@gmail.com (LUAN.van.THAC.si).mot.so.thuat.toan.chieu.diem.gan.ke.giai.phuong.trinh.voi.toan.tu.don.dieu(LUAN.van.THAC.si).mot.so.thuat.toan.chieu.diem.gan.ke.giai.phuong.trinh.voi.toan.tu.don.dieu(LUAN.van.THAC.si).mot.so.thuat.toan.chieu.diem.gan.ke.giai.phuong.trinh.voi.toan.tu.don.dieu(LUAN.van.THAC.si).mot.so.thuat.toan.chieu.diem.gan.ke.giai.phuong.trinh.voi.toan.tu.don.dieu (LUAN.van.THAC.si).mot.so.thuat.toan.chieu.diem.gan.ke.giai.phuong.trinh.voi.toan.tu.don.dieu(LUAN.van.THAC.si).mot.so.thuat.toan.chieu.diem.gan.ke.giai.phuong.trinh.voi.toan.tu.don.dieu(LUAN.van.THAC.si).mot.so.thuat.toan.chieu.diem.gan.ke.giai.phuong.trinh.voi.toan.tu.don.dieu(LUAN.van.THAC.si).mot.so.thuat.toan.chieu.diem.gan.ke.giai.phuong.trinh.voi.toan.tu.don.dieu 40 Chương Phương pháp CQ Vậy ta có lim kzk − Ti zk k = (vì kzk − Ti zk k ≤ kzk − yki k + kTi zk − yki k) k→∞ Từ tính bị chặn {zk } nên tồn dãy {zkj } cho zkj ⇀ zb lim kzkj − j→∞ Ti zkj k = Từ tính khơng giãn tương đối Ti suy zb ∈ Fb(Ti ) = F (Ti ) với i = 1, 2, , N Vì zb ∈ F Sử dụng (3.11) ta có kzkj − PF (x0 )k2 = kPC (xkj ) − PC (PF (x0 ))k2 = kxkj − x0 k2 + kx0 − PF (x0 )k2 − 2hxkj − x0 , PF (x0 ) − x0 i ≤ 2kx0 − PF (x0 )k2 − 2hxkj − x0 , PF (x0 ) − x0 i = 2(kx0 − PF (x0 )k2 − hzkj − x0 , PF (x0 ) − x0 i + hzkj − xkj , PF (x0 ) − x0 i) Vì từ (3.13) zkj → zb ta suy lim sup kzkj − PF (x0 )k2 ≤ 2(kx0 − PF (x0 )k2 − hb z − x0 , PF (x0 ) − x0 i (3.14) j→∞ Vì zb ∈ F F lồi nên ta có kx0 − PF (x0 )k2 − hb z − x0 , PF (x0 ) − x0 i = hx0 − PF (x0 ), zb − PF (x0 )i ≤ Kết hợp với (3.14) suy lim kzkj − PF (x0 )k = j→∞ zkj → PF (x0 ) j → ∞ Vì nên PF (x0 ) điểm tụ yếu {zk } Rõ ràng dãy hội tụ yếu {zk } hội tụ mạnh PF (x0 ) zk → PF (x0 ) k → ∞ Từ lim kzk − xk k = ta có lim xk = PF (x0 ) k→∞ k→∞  Trong tài liệu [2], tác giả thử nghiệm Thuật tốn CQ xoay vịng Liu Thuật tốn CQ song song đề tìm điểm bất động hai tốn tử tích phân phi tuyến khơng gian Hilbert H = L2 [0, 1] Tính tốn cho thấy, chế độ tuần tự, tức sử dụng xử lý thuật toán CQ song song hiệu thuật tốn CQ xoay vịng Cịn chế độ chạy song song với hai xử lý, để đạt độ xác thuật tốn CQ song song tốn thời gian hẳn thuật tốn CQ xoay vòng TIEU LUAN MOI download : skknchat@gmail.com (LUAN.van.THAC.si).mot.so.thuat.toan.chieu.diem.gan.ke.giai.phuong.trinh.voi.toan.tu.don.dieu(LUAN.van.THAC.si).mot.so.thuat.toan.chieu.diem.gan.ke.giai.phuong.trinh.voi.toan.tu.don.dieu(LUAN.van.THAC.si).mot.so.thuat.toan.chieu.diem.gan.ke.giai.phuong.trinh.voi.toan.tu.don.dieu(LUAN.van.THAC.si).mot.so.thuat.toan.chieu.diem.gan.ke.giai.phuong.trinh.voi.toan.tu.don.dieu (LUAN.van.THAC.si).mot.so.thuat.toan.chieu.diem.gan.ke.giai.phuong.trinh.voi.toan.tu.don.dieu(LUAN.van.THAC.si).mot.so.thuat.toan.chieu.diem.gan.ke.giai.phuong.trinh.voi.toan.tu.don.dieu(LUAN.van.THAC.si).mot.so.thuat.toan.chieu.diem.gan.ke.giai.phuong.trinh.voi.toan.tu.don.dieu(LUAN.van.THAC.si).mot.so.thuat.toan.chieu.diem.gan.ke.giai.phuong.trinh.voi.toan.tu.don.dieu Chương Áp dụng 4.1 Bài tốn khơi phục ảnh Bài tốn Tìm vector x ∈ H, biết hướng quan sát vi ∈ H hình chiếu hx, vi i = µi , i = 1, 2, , n Gọi Pi (x) phép chiếu trực giao x lên siêu phẳng Hi = {ξ ∈ H | hξ, vi i = µi } Ta có Pi (x) = x − Với x, y ∈ H ta có i = 1, 2, , n hx, vi i − µi vi kvi k2 hx, vi i − hy, vi i vi kvi k2 hx − y, vi i ⇒ x − y = Pi (x) − Pi (y) + vi kvi k2 hx − y, vi i2 ⇒ kx − yk2 = kPi (x) − Pi (y)k2 + kvi k2 Pi (x) − Pi (y) = x − y − ⇒ kPi (x) − Pi (y)k ≤ kx − yk Do Pi ánh xạ khơng giãn Đặt Ai = I − Pi , I tốn tử đồng Khi ta có Ai (x) = hx, vi i − µi vi kvi k2 TIEU LUAN MOI download : skknchat@gmail.com (LUAN.van.THAC.si).mot.so.thuat.toan.chieu.diem.gan.ke.giai.phuong.trinh.voi.toan.tu.don.dieu(LUAN.van.THAC.si).mot.so.thuat.toan.chieu.diem.gan.ke.giai.phuong.trinh.voi.toan.tu.don.dieu(LUAN.van.THAC.si).mot.so.thuat.toan.chieu.diem.gan.ke.giai.phuong.trinh.voi.toan.tu.don.dieu(LUAN.van.THAC.si).mot.so.thuat.toan.chieu.diem.gan.ke.giai.phuong.trinh.voi.toan.tu.don.dieu (LUAN.van.THAC.si).mot.so.thuat.toan.chieu.diem.gan.ke.giai.phuong.trinh.voi.toan.tu.don.dieu(LUAN.van.THAC.si).mot.so.thuat.toan.chieu.diem.gan.ke.giai.phuong.trinh.voi.toan.tu.don.dieu(LUAN.van.THAC.si).mot.so.thuat.toan.chieu.diem.gan.ke.giai.phuong.trinh.voi.toan.tu.don.dieu(LUAN.van.THAC.si).mot.so.thuat.toan.chieu.diem.gan.ke.giai.phuong.trinh.voi.toan.tu.don.dieu 42 Chương Áp dụng hAi (x) − Ai (y), x − yi = kx − yk2 − hPi (x) − Pi (y), x − yi Hơn nữa, kAi (x) − Ai (y)k2 = kx − yk2 + kPi (x) − Pi (y)k2 − 2hx − y, Pi (x) − Pi (y)i Vì kPi (x) − Pi (y)k ≤ kx − yk2 nên hAi (x) − Ai (y), x − yi ≥ kAi (x) − Ai (y)k2 ≥ Vậy Ai toán tử ngược đơn điệu mạnh 4.1.1 Phương pháp chiếu-điểm gần kề song song Bài tốn khơi phục ảnh phát biểu lại thành tốn giải hệ phương trình với toán tử đơn điệu Ai (x) = i = 1, 2, , n (4.1) Áp dụng Thuật toán chiếu - điểm gần kề song song để giải hệ phương trình (4.1) Cho trước x0 ∈ H µik > Bước Giả sử bước lặp thứ k ≥ ta có xk Tìm yki ∈ H thỏa mãn phương trình Ai (yki ) + µik (yki − xk ) = hay i = 1, 2, , n, hyki , vi i − µi vi + µik (yki − xk ) = kvi k2 i = 1, 2, , n Nhân vô hướng hai vế phương trình với vi ta hyki , vi i − µi + µik hyki , vi i − µik hxk , vi i =
⇒ + µik hyki , vi i = µi + µik hxk , vi i
⇒ hyki , vi i = µi + µik hxk , vi i i + µk Đặt αi = hyki , vi i, ta thu yki = xk − αi − µi vi µik kvi k2 TIEU LUAN MOI download : skknchat@gmail.com (LUAN.van.THAC.si).mot.so.thuat.toan.chieu.diem.gan.ke.giai.phuong.trinh.voi.toan.tu.don.dieu(LUAN.van.THAC.si).mot.so.thuat.toan.chieu.diem.gan.ke.giai.phuong.trinh.voi.toan.tu.don.dieu(LUAN.van.THAC.si).mot.so.thuat.toan.chieu.diem.gan.ke.giai.phuong.trinh.voi.toan.tu.don.dieu(LUAN.van.THAC.si).mot.so.thuat.toan.chieu.diem.gan.ke.giai.phuong.trinh.voi.toan.tu.don.dieu (LUAN.van.THAC.si).mot.so.thuat.toan.chieu.diem.gan.ke.giai.phuong.trinh.voi.toan.tu.don.dieu(LUAN.van.THAC.si).mot.so.thuat.toan.chieu.diem.gan.ke.giai.phuong.trinh.voi.toan.tu.don.dieu(LUAN.van.THAC.si).mot.so.thuat.toan.chieu.diem.gan.ke.giai.phuong.trinh.voi.toan.tu.don.dieu(LUAN.van.THAC.si).mot.so.thuat.toan.chieu.diem.gan.ke.giai.phuong.trinh.voi.toan.tu.don.dieu 43 Chương Áp dụng µi + µik hxk , vi i − µi − µi µik = xk − vi (1 + µik )µik kvi k2 hxk , vi i − µi = xk − vi (1 + µik )kvi k2 Bước Tìm số tối ưu ik : Hki = {z ∈ Rn | hz − yki , Ai (yki )i ≤ 0} = {z ∈ Rn | hz − yki , xk − yki i ≤ 0} Với z ∈ Hki ta có hz − yki , xk − yki i ≤ ⇒ hz − xk + xk − yki , xk − yki i ≤ ⇒ kxk − yki k2 ≤ hxk − z, xk − yki i ≤ kxk − zkkxk − yki k ⇒ kxk − yki k ≤ kxk − zk, suy PHki (xk ) = yki xk − PHki (xk ) = xk − yki = hxk , vi i − µi vi (1 + µik )kvi k2 |hxk , vi i − µi | ta có (1 + µik )kvi k   |hxk , vi i − µi | ik = argmax (1 + µik )kvi k i=1,2, ,n suy kxk − PHki (xk )k = Bước Tính xk+1 Đặt Wk = {z ∈ H | hz − xk , x0 − xk i ≤ 0} Ta có hAik (ykik ), x0 − ykik i PH ik (x0 ) = x0 − Aik (ykik ) i k k kAik (yk )k ik µk hxk − ykik , x0 − ykik i ik = x0 − µk (xk − ykik ) ik ik (µk ) kxk − yk k hxk − ykik , x0 − xk i + kxk − ykik k2 (xk − ykik ) = x0 − ik kxk − yk k TIEU LUAN MOI download : skknchat@gmail.com (LUAN.van.THAC.si).mot.so.thuat.toan.chieu.diem.gan.ke.giai.phuong.trinh.voi.toan.tu.don.dieu(LUAN.van.THAC.si).mot.so.thuat.toan.chieu.diem.gan.ke.giai.phuong.trinh.voi.toan.tu.don.dieu(LUAN.van.THAC.si).mot.so.thuat.toan.chieu.diem.gan.ke.giai.phuong.trinh.voi.toan.tu.don.dieu(LUAN.van.THAC.si).mot.so.thuat.toan.chieu.diem.gan.ke.giai.phuong.trinh.voi.toan.tu.don.dieu (LUAN.van.THAC.si).mot.so.thuat.toan.chieu.diem.gan.ke.giai.phuong.trinh.voi.toan.tu.don.dieu(LUAN.van.THAC.si).mot.so.thuat.toan.chieu.diem.gan.ke.giai.phuong.trinh.voi.toan.tu.don.dieu(LUAN.van.THAC.si).mot.so.thuat.toan.chieu.diem.gan.ke.giai.phuong.trinh.voi.toan.tu.don.dieu(LUAN.van.THAC.si).mot.so.thuat.toan.chieu.diem.gan.ke.giai.phuong.trinh.voi.toan.tu.don.dieu 44 Chương Áp dụng hxk − ykik , x0 − xk i = x0 − (xk − ykik ) − (xk − ykik ) kxk − ykik k2 hxk , vik i − µik hvi , x0 − xk i − v vik = x0 − k i k kvik k2 (1 + µikk )kvik k2 (1 + µikk )hvik , x0 i − µikk hvik , xk i − µik = x0 − vik (1 + µikk )kvik k2 Nếu PH ik (x0 ) ∈ Wk , tức k hPH ik (x0 ) − xk , x0 − xk i ≤ k (1 + µikk )hvik , x0 i − µikk hvik , xk i − µik hvik , x0 − xk i ⇔ kx0 − xk k ≤ (1 + µikk )kvik k2 xk+1 = PH ik (x0 ) k Nếu khơng xk+1 xác định sau xk+1 = PH ik ∩W (x0 ) = x0 + λ1 Aik (ykik ) + λ2 (x0 − xk ), k k λ1 , λ2 thỏa mãn hệ phương trình  λ1 kAi (y ik )k2 + λ2 hAi (y ik ), x0 − xk i = −hx0 − y ik , Ai (y ik )i, k k k k k k k λ hA (y ik ), x − x i + λ kx − x k2 = −kx − x k2 4.1.2 ik k k 0 k k Thử nghiệm số Xét hệ phương trình hx, vi i = µi , H = L2 [0, 1], vi = ti , µi = i+2 i = 1, 2, 3, 4, nghiệm xác tốn x∗ (t) = t Ta có αi = hyki , vi i = kvi k2 = Z1   + µik + µik + i t2i dt = , + 2i Z1  xk ti dt , TIEU LUAN MOI download : skknchat@gmail.com (LUAN.van.THAC.si).mot.so.thuat.toan.chieu.diem.gan.ke.giai.phuong.trinh.voi.toan.tu.don.dieu(LUAN.van.THAC.si).mot.so.thuat.toan.chieu.diem.gan.ke.giai.phuong.trinh.voi.toan.tu.don.dieu(LUAN.van.THAC.si).mot.so.thuat.toan.chieu.diem.gan.ke.giai.phuong.trinh.voi.toan.tu.don.dieu(LUAN.van.THAC.si).mot.so.thuat.toan.chieu.diem.gan.ke.giai.phuong.trinh.voi.toan.tu.don.dieu (LUAN.van.THAC.si).mot.so.thuat.toan.chieu.diem.gan.ke.giai.phuong.trinh.voi.toan.tu.don.dieu(LUAN.van.THAC.si).mot.so.thuat.toan.chieu.diem.gan.ke.giai.phuong.trinh.voi.toan.tu.don.dieu(LUAN.van.THAC.si).mot.so.thuat.toan.chieu.diem.gan.ke.giai.phuong.trinh.voi.toan.tu.don.dieu(LUAN.van.THAC.si).mot.so.thuat.toan.chieu.diem.gan.ke.giai.phuong.trinh.voi.toan.tu.don.dieu 45 Chương Áp dụng (2i + 1)(αi − µi ) i t, µik (2i + 1)(αi − µi ) i t, xk − yki = µik yki = xk − Tính tốn thực chế độ với x0 = 104 × t µik = với i, k, sai số εk tính theo cơng thức ε2k ∗ = kxk − x k = Z1 (xk − t)2 dt Bảng 4.1: Thuật toán chiếu-điểm gần kề song song 4.2 Số lần lặp (k) Sai số (εk ) 2.89 × 103 10 30 5.6376 5.36 × 10−6 40 5.25 × 10−9 Giải hệ phương trình đại số tuyến tính Xét hệ phương trình Ax = b, A ma trận thực kích thước m × n b vectơ không gian Rm Đặt aTi = (ai1 , , ain ) vector hàng thứ i ma trận A Viết lại hệ phương trình dạng  aT x = bi i = 1, , m, i a = (a , , a )T , x ∈ Rn , b ∈ R i i1 in i (4.2) Gọi Pi (x) phép chiếu trực giao x lên siêu phẳng Hi = {z ∈ Rn | aTi z = bi } Ta có Pi (x) = x − i = 1, , m aTi x − bi kai k2 TIEU LUAN MOI download : skknchat@gmail.com (LUAN.van.THAC.si).mot.so.thuat.toan.chieu.diem.gan.ke.giai.phuong.trinh.voi.toan.tu.don.dieu(LUAN.van.THAC.si).mot.so.thuat.toan.chieu.diem.gan.ke.giai.phuong.trinh.voi.toan.tu.don.dieu(LUAN.van.THAC.si).mot.so.thuat.toan.chieu.diem.gan.ke.giai.phuong.trinh.voi.toan.tu.don.dieu(LUAN.van.THAC.si).mot.so.thuat.toan.chieu.diem.gan.ke.giai.phuong.trinh.voi.toan.tu.don.dieu (LUAN.van.THAC.si).mot.so.thuat.toan.chieu.diem.gan.ke.giai.phuong.trinh.voi.toan.tu.don.dieu(LUAN.van.THAC.si).mot.so.thuat.toan.chieu.diem.gan.ke.giai.phuong.trinh.voi.toan.tu.don.dieu(LUAN.van.THAC.si).mot.so.thuat.toan.chieu.diem.gan.ke.giai.phuong.trinh.voi.toan.tu.don.dieu(LUAN.van.THAC.si).mot.so.thuat.toan.chieu.diem.gan.ke.giai.phuong.trinh.voi.toan.tu.don.dieu 46 Chương Áp dụng Đặt Ai = I − Pi Tương tự Bài tốn khơi phục ảnh, ta có Pi ánh xạ khơng giãn Ai toán tử ngược đơn điệu mạnh Bài toán giải hệ phương trình (4.2) tương đương với tốn giải hệ phương trình với tốn tử đơn điệu Ai (x) = 4.2.1 i = 1, , m Phương pháp chiếu-điểm gần kề song song Đây Bài tốn khơi phục ảnh với H = Rn , vi = µi = bi Áp dụng Thuật toán chiếu - điểm gần kề song song ta thu thuật tốn giải hệ phương trình (4.2) sau Cho trước x0 ∈ Rn µik > lần lặp thứ k ≥ ta có xk • Tính (đồng thời) yki = xk − • Đặt hxk , i − bi (1 + µik )kai k2 σki := Tìm i = 1, , m (4.3) |hxk , i − bi | (1 + µik )kai k  ik = argmax σki i=1,2, ,n • Tính PH ik (x0 ) = x0 − k (1 + µikk )haik , x0 i − µikk haik , xk i − bik aik (1 + µikk )kaik k2 Nếu (1 + µikk )haik , x0 i − µikk haik , xk i − bik kx0 − xk k ≤ haik , x0 − xk i (1 + µikk )kaik k2 xk+1 = PH ik (x0 ) k Nếu khơng tìm λ1 λ2 thỏa mãn hệ phương trình  λ1 kAi (y ik )k2 + λ2 hAi (y ik ), x0 − xk i = −hx0 − y ik , Ai (y ik )i, k k k k k k k λ hA (y ik ), x − x i + λ kx − x k2 = −kx − x k2 ik k k k k TIEU LUAN MOI download : skknchat@gmail.com (LUAN.van.THAC.si).mot.so.thuat.toan.chieu.diem.gan.ke.giai.phuong.trinh.voi.toan.tu.don.dieu(LUAN.van.THAC.si).mot.so.thuat.toan.chieu.diem.gan.ke.giai.phuong.trinh.voi.toan.tu.don.dieu(LUAN.van.THAC.si).mot.so.thuat.toan.chieu.diem.gan.ke.giai.phuong.trinh.voi.toan.tu.don.dieu(LUAN.van.THAC.si).mot.so.thuat.toan.chieu.diem.gan.ke.giai.phuong.trinh.voi.toan.tu.don.dieu (LUAN.van.THAC.si).mot.so.thuat.toan.chieu.diem.gan.ke.giai.phuong.trinh.voi.toan.tu.don.dieu(LUAN.van.THAC.si).mot.so.thuat.toan.chieu.diem.gan.ke.giai.phuong.trinh.voi.toan.tu.don.dieu(LUAN.van.THAC.si).mot.so.thuat.toan.chieu.diem.gan.ke.giai.phuong.trinh.voi.toan.tu.don.dieu(LUAN.van.THAC.si).mot.so.thuat.toan.chieu.diem.gan.ke.giai.phuong.trinh.voi.toan.tu.don.dieu 47 Chương Áp dụng tính xk+1 = x0 + λ1 Aik (ykik ) + λ2 (x0 − xk ) Ý nghĩa hình học • Ta có Ai (yki ) + µik (yki − xk ) = ⇔ yki − Pi (yki ) + µik (yki − xk ) = ⇔ (1 + µik )yki = Pi (yki ) + µik xk ⇔ yki = µik i P (y ) + xk i k + µik + µik (4.4) Do µik > nên phương trình (4.4) cho thấy yki nằm đoạn thẳng [xk , Pi (yki )] Hơn nữa, từ (4.4) (4.3) ta thu Pi (yki ) = (1 + µik )yki − µik xk   hxk , i − bi i − µik xk = (1 + µk ) xk − (1 + µik )kai k2 hxk , i − bi = xk − kai k2 = Pi (xk ) Như vậy, điểm yki nằm đoạn thẳng [xk , Pi (xk )] • Xét siêu phẳng ′ Hki = {z ∈ Rn | hz − yki , xk − yki i = 0} = {z ∈ Rn | hz, xk − yki i = hyki , xk − yki i}   (1 + µik )kai k2 i i n hy , xk − yk i = z ∈ R | hz, i = hxk , i − µik k Điều cho thấy Hki song song với Hi ′ • Xét siêu phẳng Wk′ = {z ∈ Rn | hz − xk , x0 − xk i = 0} = {z ∈ Rn | hz, x0 − xk i = hxk , x0 − xk i} Vậy Wk′ siêu phẳng qua điểm xk vng góc với vector x0 − xk TIEU LUAN MOI download : skknchat@gmail.com (LUAN.van.THAC.si).mot.so.thuat.toan.chieu.diem.gan.ke.giai.phuong.trinh.voi.toan.tu.don.dieu(LUAN.van.THAC.si).mot.so.thuat.toan.chieu.diem.gan.ke.giai.phuong.trinh.voi.toan.tu.don.dieu(LUAN.van.THAC.si).mot.so.thuat.toan.chieu.diem.gan.ke.giai.phuong.trinh.voi.toan.tu.don.dieu(LUAN.van.THAC.si).mot.so.thuat.toan.chieu.diem.gan.ke.giai.phuong.trinh.voi.toan.tu.don.dieu (LUAN.van.THAC.si).mot.so.thuat.toan.chieu.diem.gan.ke.giai.phuong.trinh.voi.toan.tu.don.dieu(LUAN.van.THAC.si).mot.so.thuat.toan.chieu.diem.gan.ke.giai.phuong.trinh.voi.toan.tu.don.dieu(LUAN.van.THAC.si).mot.so.thuat.toan.chieu.diem.gan.ke.giai.phuong.trinh.voi.toan.tu.don.dieu(LUAN.van.THAC.si).mot.so.thuat.toan.chieu.diem.gan.ke.giai.phuong.trinh.voi.toan.tu.don.dieu 48 Chương Áp dụng • Minh họa hình học trường hợp n = 2, m = µik = với k, i giả sử hệ phương trình có nghiệm Trong trường hợp này, Hi đường thẳng với µik = yki trung điểm đoạn thẳng [xk , Pi (xk )] (với k, i) Hình 4.1: Minh họa hình học thuật tốn chiếu-điểm gần kề song song 4.2.2 Phương pháp CQ song song Giải hệ phương trình (4.2) tương đương với tốn tìm điểm bất động họ tốn tử khơng giãn Pi (x) từ Rn vào Với k ≥ 0, hình chiếu xk lên C = Rn nó, tức zk = PRn (xk ) = xk Hơn nữa, ta có yki  aTi xk − bi = αk x0 + (1 − αk )Pi (xk ) = αk x0 + (1 − αk ) xk − kai k2  Viết lại Ck dạng Ck = {v ∈ Rn | hv, xk − Pi (xk )i ≤ µk }, xk − Pi (xk ) = µk = aTi xk − bi kai k2
αk kx0 k2 + (1 − αk )kxk k2 − kyki k2 2(1 − αk ) TIEU LUAN MOI download : skknchat@gmail.com (LUAN.van.THAC.si).mot.so.thuat.toan.chieu.diem.gan.ke.giai.phuong.trinh.voi.toan.tu.don.dieu(LUAN.van.THAC.si).mot.so.thuat.toan.chieu.diem.gan.ke.giai.phuong.trinh.voi.toan.tu.don.dieu(LUAN.van.THAC.si).mot.so.thuat.toan.chieu.diem.gan.ke.giai.phuong.trinh.voi.toan.tu.don.dieu(LUAN.van.THAC.si).mot.so.thuat.toan.chieu.diem.gan.ke.giai.phuong.trinh.voi.toan.tu.don.dieu (LUAN.van.THAC.si).mot.so.thuat.toan.chieu.diem.gan.ke.giai.phuong.trinh.voi.toan.tu.don.dieu(LUAN.van.THAC.si).mot.so.thuat.toan.chieu.diem.gan.ke.giai.phuong.trinh.voi.toan.tu.don.dieu(LUAN.van.THAC.si).mot.so.thuat.toan.chieu.diem.gan.ke.giai.phuong.trinh.voi.toan.tu.don.dieu(LUAN.van.THAC.si).mot.so.thuat.toan.chieu.diem.gan.ke.giai.phuong.trinh.voi.toan.tu.don.dieu 49 Chương Áp dụng Nếu PCk ∩Qk (x0 ) 6∈ Qk hình chiếu x0 lên Ck ∩ Qk viết dạng PCk ∩Qk (x0 ) = x0 + λ1 (xk − Pik (xk )) + λ2 (x0 − xk ), λ1 , λ2 thỏa mãn hệ phương trình  λ1 kxk − Pi (xk )k2 + λ2 hxk − Pi (xk ), x0 − xk i = µk − hx0 , xk − Pi (xk )i, k k k λ hx − P (x ), x − x i + λ kx − x k2 = −kx − x k2 k ik k k k k (4.3) Tóm lại, áp dụng Thuật toán CQ song song giải hệ phương trình (4.2) sau Cho x0 ∈ Rn dãy {αk } ⊂ (0, 1) Giả sử lần lặp thứ k ≥ ta có xk • Tính (đồng thời) yki • Tìm   aTi xk − bi = αk x0 + (1 − αk ) xk − kai k2 ik = argmax{kyki − xk k} • Tìm xk+1 :    xk xk ∈ Ck ,    xk+1 = PCk (x0 ) PCk (x0 ) ∈ Qk ,     x0 + λ1 (xk − Pi (xk )) + λ2 (x0 − xk ) trường hợp cịn lại, k λ1 , λ2 nghiệm hệ phương trình (4.3) PCk (x0 ) = x0 − hx0 , xk − Pik (xk )i − µk (xk − Pik (xk )) kxk − Pik (xk )k2 TIEU LUAN MOI download : skknchat@gmail.com (LUAN.van.THAC.si).mot.so.thuat.toan.chieu.diem.gan.ke.giai.phuong.trinh.voi.toan.tu.don.dieu(LUAN.van.THAC.si).mot.so.thuat.toan.chieu.diem.gan.ke.giai.phuong.trinh.voi.toan.tu.don.dieu(LUAN.van.THAC.si).mot.so.thuat.toan.chieu.diem.gan.ke.giai.phuong.trinh.voi.toan.tu.don.dieu(LUAN.van.THAC.si).mot.so.thuat.toan.chieu.diem.gan.ke.giai.phuong.trinh.voi.toan.tu.don.dieu (LUAN.van.THAC.si).mot.so.thuat.toan.chieu.diem.gan.ke.giai.phuong.trinh.voi.toan.tu.don.dieu(LUAN.van.THAC.si).mot.so.thuat.toan.chieu.diem.gan.ke.giai.phuong.trinh.voi.toan.tu.don.dieu(LUAN.van.THAC.si).mot.so.thuat.toan.chieu.diem.gan.ke.giai.phuong.trinh.voi.toan.tu.don.dieu(LUAN.van.THAC.si).mot.so.thuat.toan.chieu.diem.gan.ke.giai.phuong.trinh.voi.toan.tu.don.dieu Kết luận Luận văn trình bày lại cách hệ thống thuật toán hội tụ mạnh giải phương trình, hệ phương trình tốn tử đơn điệu tìm điểm bất động họ hữu hạn ánh xạ không giãn tương đối không gian Hilbert dựa phương pháp lai ghép thuật toán điểm gần kề với phép chiếu lên giao tập lồi nửa không gian Luận văn đề cập đến vấn đề sau: Hệ thống lại số phương pháp điểm gần kề cải biên khơng gian Hilbert Khái niệm nghiệm gần với sai số cho phép số tính chất quan trọng; Trình bày thuật tốn chiếu-điểm gần kề để tìm khơng điểm tốn tử đơn điệu cực đại định lý hội tụ mạnh đến điểm tập khơng điểm; Trình bày thuật toán song song chiếu-điểm gần kề CQ tương ứng để giải hệ phương trình với tốn tử đơn điệu tìm điểm bất động họ ánh xạ không giãn tương đối không gian Hilbert Banach; Trình bày định lý hội tụ mạnh cho tốn tử khơng gian tương đối khơng gian Banach định lý hội tụ mạnh phương pháp CQ xoay vòng CQ song song cho họ hữu hạn ánh xạ không giãn tương đối khơng gian Banach; Áp dụng thuật tốn chiếu-điểm gần kề song song cho tốn khơi phục ảnh khơng gian Hilbert ví dụ số minh họa cho hội tụ thuật toán này; Áp dụng thuật toán chiếu-điểm gần kề song song phương pháp CQ song song để giải hệ phương trình với số ẩn số phương trình khơng thiết Minh họa hình học cho thuật tốn chiếu-điểm gần kề song song R2 TIEU LUAN MOI download : skknchat@gmail.com (LUAN.van.THAC.si).mot.so.thuat.toan.chieu.diem.gan.ke.giai.phuong.trinh.voi.toan.tu.don.dieu(LUAN.van.THAC.si).mot.so.thuat.toan.chieu.diem.gan.ke.giai.phuong.trinh.voi.toan.tu.don.dieu(LUAN.van.THAC.si).mot.so.thuat.toan.chieu.diem.gan.ke.giai.phuong.trinh.voi.toan.tu.don.dieu(LUAN.van.THAC.si).mot.so.thuat.toan.chieu.diem.gan.ke.giai.phuong.trinh.voi.toan.tu.don.dieu (LUAN.van.THAC.si).mot.so.thuat.toan.chieu.diem.gan.ke.giai.phuong.trinh.voi.toan.tu.don.dieu(LUAN.van.THAC.si).mot.so.thuat.toan.chieu.diem.gan.ke.giai.phuong.trinh.voi.toan.tu.don.dieu(LUAN.van.THAC.si).mot.so.thuat.toan.chieu.diem.gan.ke.giai.phuong.trinh.voi.toan.tu.don.dieu(LUAN.van.THAC.si).mot.so.thuat.toan.chieu.diem.gan.ke.giai.phuong.trinh.voi.toan.tu.don.dieu Tài liệu tham khảo [1] Lê Dũng Mưu Nguyễn Văn Hiền (2009), Nhập môn giải tích lồi ứng dụng, Nhà xuất Khoa học tự nhiên Công nghệ, Hà Nội [2] P K Anh and C V Chung (2011), A parallel CQ method for a finite family of relatively nonexpansive mappings (Submitted for publication) [3] P K Anh and C V Chung (2011), On strongly convergent parallel proximal point algorithms, Journal of Science, VNU, 27(2) [4] O Guler (1991), On the convergence of the proximal point algorithm for convex minimization, SIAM J Optim., 2, 649-664 [5] S Kamimura and W Takahashi (2002), Strong convergence of a proximaltype algorithm in a Banach space, SIAM J Optim., 13, 938-945 [6] X F Liu (2011), Strong convergence theorems for a finite family of relatively nonexpansive mappings, Vietnam J Math., 39, 63-69 [7] S Matsushita and W Takahashi (2005), A strong convergence theorem for relatively nonexpansive mappings in a Banach space, J Approx Theory, 134, 257-266 [8] R T Rockafellar (1976), Monotone operators and proximal point algorithm, SIAM J Contr Optim., 14, 877-897 [9] X L Qin and Y F Su (2007), Strong convergence theorems for relatively nonexpansive mappings in a Banach space, Nonlinear Anal., 67, 1958-1965 [10] M V Solodov and B F Svaiter (1999), A hybrid projection-proximal point algorithm, J Conv Anal., 6, 59-70 TIEU LUAN MOI download : skknchat@gmail.com (LUAN.van.THAC.si).mot.so.thuat.toan.chieu.diem.gan.ke.giai.phuong.trinh.voi.toan.tu.don.dieu(LUAN.van.THAC.si).mot.so.thuat.toan.chieu.diem.gan.ke.giai.phuong.trinh.voi.toan.tu.don.dieu(LUAN.van.THAC.si).mot.so.thuat.toan.chieu.diem.gan.ke.giai.phuong.trinh.voi.toan.tu.don.dieu(LUAN.van.THAC.si).mot.so.thuat.toan.chieu.diem.gan.ke.giai.phuong.trinh.voi.toan.tu.don.dieu (LUAN.van.THAC.si).mot.so.thuat.toan.chieu.diem.gan.ke.giai.phuong.trinh.voi.toan.tu.don.dieu(LUAN.van.THAC.si).mot.so.thuat.toan.chieu.diem.gan.ke.giai.phuong.trinh.voi.toan.tu.don.dieu(LUAN.van.THAC.si).mot.so.thuat.toan.chieu.diem.gan.ke.giai.phuong.trinh.voi.toan.tu.don.dieu(LUAN.van.THAC.si).mot.so.thuat.toan.chieu.diem.gan.ke.giai.phuong.trinh.voi.toan.tu.don.dieu (LUAN.van.THAC.si).mot.so.thuat.toan.chieu.diem.gan.ke.giai.phuong.trinh.voi.toan.tu.don.dieu(LUAN.van.THAC.si).mot.so.thuat.toan.chieu.diem.gan.ke.giai.phuong.trinh.voi.toan.tu.don.dieu(LUAN.van.THAC.si).mot.so.thuat.toan.chieu.diem.gan.ke.giai.phuong.trinh.voi.toan.tu.don.dieu(LUAN.van.THAC.si).mot.so.thuat.toan.chieu.diem.gan.ke.giai.phuong.trinh.voi.toan.tu.don.dieu