Toán tử không giãn trung bình và ứng dụng

ĐẠI HỌC QUỐC GIA HÀ NỘI TRƯỜNG ĐẠI HỌC KHOA HỌC TỰ NHIÊN KHOA TOÁN - CƠ - TIN HỌC Kiều Thị Thùy Linh TOÁN TỬ KHÔNG GIÃN TRUNG BÌNH VÀ ỨNG DỤNG LUẬN VĂN THẠC SĨ Chuyên ngành: Toán Giải tích Hà Nội - 2017 ĐẠI HỌC QUỐC GIA HÀ NỘI TRƯỜNG ĐẠI HỌC KHOA HỌC TỰ NHIÊN KHOA TOÁN - CƠ - TIN HỌC Kiều Thị Thùy Linh TOÁN TỬ KHÔNG GIÃN TRUNG BÌNH VÀ ỨNG DỤNG LUẬN VĂN THẠC SĨ Chuyên ngành: Toán Giải tích Mã số: 60460102 Cán hướng dẫn: GS.TSKH Phạm Kỳ Anh Hà Nội - 2017 LỜI CẢM ƠN Với lòng kính trọng biết ơn sâu sắc, xin bày tỏ lời cảm ơn chân thành tới GS.TSKH Phạm Kỳ Anh, giảng viên 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 trực tiếp giao đề tài cho bỏ công sức để hướng dẫn tận tình Thày cho kiến thức kinh nghiệm quý báu, tạo điều kiện thuận lợi cho trình thực hoàn thành luận văn hoàn thành chương trình học Thạc sĩ Tôi xin gửi lời cảm ơn tới thầy giáo, cô giáo 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 Những người Thày tràn đầy nhiệt huyết với nghề để truyền thụ cho sinh viên, học viên kiến thức suốt trình học tập khoa trường Và xin bày tỏ lòng biết ơn vô hạn tới gia đình bạn bè chỗ dựa tinh thần vững nguồn động viên sống trình học tập Hà nội, ngày 25 tháng năm 2017 Học viên Kiều Thị Thùy Linh Danh mục kí hiệu L( X, Y ) B( X ) En FixT K∗ −1 → N R R+ R++ Rn inf f sup f Ui i∈ I f u, v argmin f ∂f NC ( x ) bdC tập hợp toán tử tuyến tính liên tục từ X vào Y, tập hợp toán tử tuyến tính liên tục từ X vào X, ma trận đơn vị kích thước n × n, tập hợp điểm bất động toán tử T, toán tử liên hợp toán tử tuyến tính K, nghịch đảo số ma trận, hội tụ yếu, sụ hội tụ mạnh, tập hợp số tự nhiên, tập hợp số thực, tập hợp số thực không âm, tập hợp số thực dương, không gian thực n - chiều, infimum hàm f , supermum hàm f , giao tập hợp Ui , i ∈ I, gradient hàm f, tích vô hướng u v, chuẩn không gian Hilbert véctơ ma trận, điểm cực tiểu hàm f , vi phân hàm f , nón chuẩn tập lồi C x, biên tập C, kết thúc chứng minh Mục lục Lời mở đầu Toán tử không giãn 1.1 Toán tử không giãn 1.2 Toán tử không giãn vững 13 Toán tử không giãn trung bình 19 2.1 Toán tử không giãn trung bình 19 2.2 Toán tử chiếu mêtric 28 2.3 Phương pháp lai ghép 32 2.4 Phương pháp xấp xỉ gắn kết 35 Ứng dụng toán tử không giãn trung bình 43 3.1 Bài toán tối ưu có ràng buộc 43 3.2 Bài toán chấp nhận lồi 44 3.3 Kĩ thuật khôi phục đại số xử lý ảnh 45 3.4 Phương pháp ngoại suy tín hiệu với dải tần hữu hạn 46 3.5 Bài toán chấp nhận tách 48 Kết luận 51 Tài liệu tham khảo 52 LỜI MỞ ĐẦU Năm 1920, Stefan Banach chứng minh ánh xạ co T không gian mêtric đủ X có điểm bất động nhất, tức tồn x ∗ ∈ X cho Tx ∗ = x ∗ Hơn với x ∈ X, dãy quỹ đạo { T k x } hội tụ đến x ∗ Trên thực tế, nhiều vấn đề lại đưa toán tìm điểm bất động toán tử không giãn, hay điểm bất động chung họ toán tử không giãn Một ví dụ đơn giản toán tìm nghiệm hệ phương trình tuyến tính Ax = b, với ma trận chữ nhật A Mỗi nghiệm hệ coi điểm bất động chung toán tử chiếu trực giao siêu phẳng tương ứng với phương trình hệ Tuy nhiên trường hợp T toán tử không giãn ta phải bổ sung thêm số điều kiện cho không gian X Điều Browder, Gohde Kirk chứng minh vào năm 1965 Và theo định lý Banach trường hợp này, dãy lặp { T k x } nói chung không hội tụ nên người ta quan tâm đến việc xây dựng phương pháp lặp tìm điểm bất động toán tử không giãn Đầu tiên phương pháp lặp tìm nghiệm hệ phương trình đại số tuyến tính Kaczmarz đề xuất năm 1937 với phép chiếu xoay vòng Cimmino đưa năm 1938 với phép chiếu đồng thời lên siêu phẳng ứng với phương trình hệ Cả hai phương pháp hữu ích việc giải hệ phương trình cỡ lớn Vì chúng có vai trò ứng dụng quan trọng kĩ thuật chụp X quang cắt lớp máy tính Cả hai kết trở thành phép lặp tìm điểm bất động toán tử không giãn Một kết đưa năm 1950 John von Neumann phương pháp chiếu hai không gian không gian Hilbert Phương pháp cho hội tụ tới giao hai không gian Các kết Kaczmarz, Cimmino, Neumann tổng quát hóa nhiều thập kỉ Ngày nay, hội tụ phương pháp thiết lập mở rộng không cho toán tử chiếu trực giao siêu phẳng mà cho toán tử không giãn, toán tử tựa không giãn, toán tử không giãn vững, toán tử không giãn trung bình,vv Trong luận văn này, tác giả trình bày cách tổng quan toán tử không giãn, toán tử không giãn trung bình, toán tử chiếu mêtric số ứng dụng qua toán thực tế Các kiến thức tìm hiểu xét không gian Hibert thực H Bố cục luận văn bao gồm chương: • Chương luận văn trình bày khái niệm, tính chất toán tử không giãn, toán tử tựa không giãn, toán tử không giãn vững, toán tử đơn điệu mạnh ngược mối quan hệ loại toán tử Một số tính chất tập điểm bất động toán tử không giãn họ toán tử không giãn đưa • Chương luận văn trình bày toán tử không giãn trung bình, toán tử chiếu mêtric mối liên hệ toán tử với toán tử trình bày Chương Ngoài tính ổn định phép toán hợp thành, tổ hợp lồi toán tử không giãn trung bình thể rõ Cuối chương, tác giả giới thiệu hai phương pháp kết hợp để phép lặp hội tụ mạnh phương pháp lai ghép phương pháp xấp xỉ gắn kết Kĩ thuật lai ghép kết hợp với phép lặp KrasnoselskiMann cho phép bước xây dựng hai nửa không gian tách tập điểm bất động xấp xỉ ban đầu cho hình chiếu điểm ban đầu lên giao hai nửa không gian hội tụ mạnh tập điểm bất động toán tử Phương pháp xấp xỉ gắn kết mở rộng phương pháp lặp Halpern để tìm điểm bất động toán tử không giãn Phương pháp sử dụng tổ hợp lồi toán tử không giãn ánh xạ co với cách chọn trọng số thích hợp thu dãy lặp hội tụ mạnh đến điểm bất động toán tử không giãn • Chương luận văn trình bày số ứng dụng thực tế vấn đề trình bày Chương Chương Thứ toán tối ưu có ràng buộc sử dụng nhiều việc nhận tín hiệu xử lý ảnh; toán dùng hàm lồi, khả vi tập lồi đóng khác rỗng không gian Hilbert thực Tiếp theo toán chấp nhận lồi, toán chấp nhận tách có nhiều ứng dụng khoa học công nghệ mà điển hình phương pháp xạ trị với cường độ thay đổi IRMT (Intensively Modulated Radiation Therapy) Đây phương pháp xạ trị tiên tiến Ngoài tác giả giới thiệu toán ngoại suy giải tần hữu hạn xử dụng kĩ thuật xử lý tín hiệu, toán có sử dụng đến phép biến đổi Fourier Cuối kĩ thuật khôi phục đại số xử lý ảnh giúp khôi phục lại ảnh gốc từ hình chiếu theo nhiều phương pháp khác Các kiến thức tìm hiểu, tham khảo trình bày luận văn chủ yếu qua tài liệu số [1-4] [10-11] Do thời gian kiến thức có hạn, luận văn không tránh khỏi hạn chế sai sót Tác giả mong nhận góp ý quý thầy cô bạn đọc Tác giả xin chân thành cảm ơn! Hà Nội, ngày 25 tháng năm 2017 Học viên Kiều Thị Thùy Linh Chương Toán tử không giãn Cho H không gian Hilbert thực với tích , chuẩn tương ứng X tập khác rỗng H 1.1 Toán tử không giãn Định nghĩa 1.1 Toán tử T : X → H gọi (i) không giãn Tx − Ty ≤ x − y với x, y ∈ X (ii) không giãn chặt Tx − Ty < x − y x − y = Tx − Ty với x, y ∈ X (iii) co với hệ số α ∈ (0, 1) Tx − Ty ≤ α x − y với x, y ∈ X Mệnh đề 1.1 Cho X tập lồi, đóng khác rỗng H Khi tập điểm bất động toán tử không giãn T : X → H tập đóng lồi Chứng minh (i) Tính đóng Lấy dãy { xk } ⊂ Fix T ⊂ X xk → x Do X tập đóng nên x ∈ X Hơn nữa, T toán tử không giãn nên T liên tục X, ta có x = lim xk = lim Txk = Tx, tức Fix T tập đóng (ii) Tính lồi Lấy x, y ∈ Fix T, x = y z = (1 − λ) x + λy, λ ∈ (0, 1) Do tính không giãn T tính dương chuẩn, ta thu hai đẳng thức sau x − Tz = Tx − Tz ≤ x − z = λ x − y , Tz − y = Tz − Ty ≤ z − y = (1 − λ) x − y Theo bất đẳng thức tam giác, ta có x−y ≤ x − Tz + Tz − y ≤ λ x − y + (1 − λ ) x − y = x−y Do ( x − Tz) + ( Tz − y) = x − y = x − Tz + Tz − y Mặt khác, tính lồi chặt chuẩn, véc tơ x − Tz, Tz − y phải đồng tuyến phương, tức tồn số α > cho Tz − y = α( x − Tz) Suy Tz = α x+ y 1+α 1+α Hơn nữa, tính không giãn T, ta có x − y = x − Tz = Tx − Tz ≤ x − z = λ x − y , 1+α α x − y = Tz − y = Tz − Ty ≤ z − y = (1 − λ) x − y 1+α α = λ = − λ Từ hai bất đẳng thức suy 1+α 1+α Vì Tz = (1 − λ) x + λy = z hay z ∈ Fix T Do { xt } dãy bị chặn Từ suy { Txt }, { f ( xt )} dãy bị chặn Hơn ta có xt − Txt = t Txt − f ( xt ) → 0, t → Tiếp theo ta chứng minh { xt } dãy compact tương đối t → Thật vậy, giả sử {tn } ⊂ (0, 1) cho tn → n → ∞ Đặt xn := xtn Ta chứng minh dãy { xn } chứa dãy hội tụ mạnh tới x˜ ∈ Fix T nghiệm toán bất đẳng thức biến phân ˜ x − x˜ ≥ 0, x ∈ Fix T ( Id − f ) x, Thật vậy, dãy { xn } bị chặn nên ta giả sử dãy { xn } hội tụ yếu tới điểm x˜ ∈ Fix T Kết hợp với xn − Txn → 0, theo Bổ đề 2.5 ta suy x˜ ∈ Fix T Mặt khác xt − x˜ = (1 − t)( Txt − x˜ ) + t( f ( xt ) − x˜ ) nên ta có xt − x˜ ˜ xt − x˜ + t f ( xt ) − x, ˜ xt − x˜ =(1 − t) Txt − x, ≤(1 − t) xt − x˜ ˜ xt − x˜ + t f ( xt ) − x, Vì xt − x˜ ˜ xt − x˜ = f ( xt ) − f ( x˜ ), xt − x˜ + f ( x˜ ) − x, ˜ xt − x˜ ≤ f ( xt ) − x, ≤α xt − x˜ ˜ xt − x˜ + f ( x˜ ) − x, Suy xt − x˜ ≤ ˜ f ( x˜ ) − x˜ xt − x, 1−α | xn − x˜ ≤ ˜ f ( x˜ ) − x˜ xn − x, 1−α Nói riêng Do xn (2.6) ˜ x˜ từ (2.6) ta suy xn → x Tiếp theo ta chứng minh x˜ ∈ Fix T nghiệm toán bất đẳng thức biến phân xt nghiệm toán điểm bất động x = (1 − t) Tx + t f ( x ) Ta có ( Id − f ) xt = − 1−t ( Id − T ) xt t 39 Do với z ∈ Fix T 1−t ( Id − T ) xt , xt − z t 1−t ( Id − T ) xt − ( Id − T )z, xt − z , =− t ( Id − f ) xt , xt − z = − hay ( Id − f ) xt , xt − z ≤ Id − T toán tử đơn điệu Thay t tn cho n → ∞ ta x˜ nghiệm toán bất đẳng thức biến phân ˜ ta giả sử xsn → x, ˆ sn → Do xˆ ∈ Fix T, ta có Để chứng minh { xt } hội tụ tới điểm x, ˜ x˜ − xˆ ≤ ( Id − f ) x, Thay đổi vai trò x˜ xˆ ta ˆ xˆ − x˜ ≤ ( Id − f ) x, Cộng vế với vế hai bất đẳng thức ta ˆ ( Id − f ) x˜ − ( Id − f ) xˆ ≤ x˜ − x, ˜ Theo Bổ đề 2.6 ta suy xˆ = x Tiếp theo ta xét dãy lặp x0 ∈ C xn+1 = (1 − αn ) Txn + αn f ( xn ), với n ≥ {αn } ⊂ (0, 1) thỏa mãn điều kiện: (i) αn → 0, ∞ (ii) ∑ αn = ∞, n =0 ∞ α n +1 = n→∞ αn (iii) ∑ |αn+1 − αn | < ∞ lim n =0 40 (2.7) Định lý 2.8 Cho H không gian Hilbert thực, C tập lồi đóng H Giả sử T : C → C toán tử không giãn với Fix T = ∅ f : C → C toán tử co Xét dãy { xn } ˜ với x˜ nghiệm toán bất đẳng thức biến sinh (2.7) Khi xn → x, phân ˜ x − x˜ ≥ 0, x ∈ Fix T ( Id − f ) x, Chứng minh Trước hết ta chứng minh dãy { xn } dãy bị chặn Lấy p ∈ Fix T, ta có xn+1 − p ≤(1 − αn ) Txn − p + αn f ( xn ) − p ≤(1 − αn ) xn − p + αn ( f ( xn ) − f ( p) + f ( p) − p ) ≤(1 − αn ) xn − p + αn (α xn − p + f ( p) − p ) ≤ (1 − (1 − α ) α n ) x n − p + α n f ( p ) − p f ( p) − p ≤ max xn − p , 1−α Từ xn − p ≤ max x0 − p , f ( p) − p 1−α , ∀n ≥ Suy dãy { xn } bị chặn hai dãy { Txn }, { f ( xn )} bị chặn Tiếp theo ta chứng minh xn+1 − xn → Thật với số M > thích hợp ta có xn+1 − xn = (1 − αn )( Txn − Txn−1 ) + (αn − αn−1 ) ( f ( xn−1 ) − Txn−1 ) + α n ( f ( x n ) − f ( x n −1 ) ≤(1 − αn ) xn − xn−1 + M αn − αn−1 + ααn xn − xn−1 =(1 − (1 − α)αn ) xn − xn−1 + M|αn − αn−1 | Theo Bổ đề 2.4, ta có xn+1 − xn → Tiếp theo ta chứng minh xn − Txn → Ta có xn − Txn ≤ xn − xn+1 + xn+1 − Txn = xn − xn+1 + αn Txn − f ( xn ) → 41 Tiếp theo ta chứng minh lim x˜ − xn , x˜ − f ( x˜ ) ≤ n→∞ Thật vậy, lấy dãy { xnk } ⊂ { xn } cho lim sup x˜ − xn , x˜ − f ( x˜ ) = lim sup x˜ − xnk , x˜ − f ( x˜ ) n→∞ Ta giả sử xnk k→∞ ¯ Kết hợp Bổ đề 2.5 ta suy x¯ ∈ Fix T Vì x ¯ x˜ − f ( x˜ ) ≤ lim sup x˜ − xn , x˜ − f ( x˜ ) = x˜ − x, n→∞ ˜ Thật ta có Cuối chứng minh xn → x xn+1 − x˜ = (1 − αn )( Txn − x˜ ) + αn ( f ( xn ) − x˜ ) =(1 − αn )2 Txn − x˜ 2 + α2n f ( xn ) − x˜ ˜ f ( xn ) − x˜ + 2αn (1 − αn ) Txn − x, ≤(1 − 2αn + α2n ) xn − x˜ + α2n f ( xn ) − x˜ ˜ f ( xn ) − f ( x˜ ) + 2αn (1 − αn ) Txn − x, ˜ f ( x˜ ) − x˜ + 2αn (1 − αn ) Txn − x, ≤(1 − 2αn + α2n + 2ααn (1 − αn )) xn − x˜ ˜ f ( x˜ ) − x˜ + αn f ( xn ) − x˜ + αn 2(1 − αn ) Txn − x, =(1 − α˜n ) xn − x˜ 2 + α˜n β˜n , với α˜n = αn (2 − αn − 2α(1 − αn )) , ˜ f ( x˜ ) − x˜ + αn f ( xn ) − x˜ 2(1 − αn ) Txn − x, β˜n = − αn − 2α(1 − αn ) ∞ Ta chứng minh α˜n → 0, ∑ α˜n = ∞ lim sup β˜n ≤ n =1 ˜ Theo Bổ đề 2.4 suy xn → x 42 n→∞ Chương Ứng dụng toán tử không giãn trung bình 3.1 Bài toán tối ưu có ràng buộc Trong thực tế, toán tối ưu có ràng buộc sử dụng nhiều việc nhận tín hiệu xử lý ảnh Xét toán: Cho C tập khác rỗng không gian Hilbert thực H Tìm điểm cực tiểu địa phương hàm lồi khả vi f : H → R tập C Ở C := { x ∈ H : ci ( x ) ≤ 0, i ∈ I } tập lồi đóng, hàm ci : H → R, i ∈ I := {1, 2, , m} hàm lồi khả vi Năm 1977, Baillon Haddad đưa kết sau Định lý 3.1 Cho f : H → R hàm lồi, khả vi Fréchet λ > Khi đó, ∇ f λ−liên tục Lipschitz, tức với x, y ∈ H ∇ f ( x ) − ∇ f (y) ≤ λ x − y , toán tử ∇ f λ1 -đơn điệu mạnh ngược Nói riêng, ∇ f toán tử không giãn ∇ f toán tử không giãn vững Nếu γ ∈ 0, λ2 toán tử G = γ∇ f 2λ -đơn điệu mạnh ngược Khi đó, toán tử A = Id − γ∇ f toán tử PC A toán tử không giãn trung bình Từ định lý, ta có hệ sau 43 Hệ 3.1 Cho f hàm lồi khả vi tập mở D ⊆ H chứa tập lồi đóng C Giả sử f đạt cực tiểu C Nếu ∇ f toán tử λ−liên tục Lipschitz D γ ∈ 0, λ2 dãy lặp sinh x k+1 = PC x k − γ∇ f ( x k ) hội tụ yếu đến điểm cực tiểu hàm f tập C, với véctơ ban đầu x0 Chứng minh Trước hết, ta chứng minh x ∗ = argmin f ( x ) ⇔ x ∗ ∈ Fix ( PC ( Id − λ∇ f )) Thật ta có x ∗ = argmin f ( x ) = argmin f ( x ) + δC ( x ) ⇔ θ ∈ ∂( f + δC ( x ∗ )) Mặt khác θ ∈ ∇ f ( x ∗ ) + NC ( x ∗ ) ⇔ −∇ f ( x ∗ ) ∈ NC ( x ∗ ), hay −∇ f ( x ∗ ), x − x ∗ ≤ 0, ∀ x ∈ C Suy x ∗ ∈ Fix ( PC ( Id − λ∇ f )) Với λ ∈ 0, λ2 Id − λ∇ f toán tử không giãn trung bình Do PC ( Id − λ∇ f ) toán tử không giãn trung bình Hơn ∇ f toán tử λ−Lipschitz nên λ1 ∇ f toán tử không giãn Theo định lý Baillon-Haddat λ1 ∇ f toán tử 1−đơn điệu mạnh ngược Ta có γ∇ f = γλ γλ −đơn điệu mạnh ngược x ∗ ∈ Fix PC ( Id − γ∇ f ) toán tử yếu đến với γλ λ∇f > 21 Khi ta có dãy quỹ đạo { x k } hội tụ 3.2 Bài toán chấp nhận lồi Bài toán thực tế lấy lại cấu trúc liệu ảnh đưa giải toán chấp nhận lồi không gian Hilbert thực H với véctơ ảnh tương ứng phần tử không gian tập lồi biểu thị giá trị bắt buộc đặt lên cấu trúc ảnh Xét toán: Cho C1 , , C M tập lồi đóng khác rỗng không gian Hilbert thực H Tìm điểm x ∈ C := M i =1 Ci Một phương pháp giải toán ta sử dụng đến toán tử chiếu mêtric Ta xét hàm gần kề f : H → R liên kết với tập C1 , , C M f (x) = M PCm x − x 2M m∑ =1 Gradient f ∇ f (x) = x − M PCm x; M m∑ =1 44 Một điểm cực tiểu f không điểm ∇ f điểm bất động toán tử không giãn trung bình A cho A= M PCm M m∑ =1 Nếu tập C = ∅ Fix A = C dãy quỹ đạo { Ak x } hội tụ yếu tới phần tử C Nếu tập C = ∅ dãy hội tụ yếu tới điểm cực tiểu f Trong trường hợp giới hạn dãy không thiết phần tử tập Cm Nếu muốn điểm cực tiểu f gần với véctơ x tập K lồi đóng khác rỗng ta sử dụng dãy lặp định nghĩa x k+1 = PK ( Ax k ) Toán tử PK A toán tử không giãn trung bình nên với véctơ ban đầu x0 dãy lặp hội tụ tới điểm cực tiểu gần K Phương pháp thường gọi phương pháp đồng ta tính toán đồng thời toán tử mêtric tập Cm bước dãy lặp Mặt khác toán tử PCm toán tử không giãn trung bình nên tích chúng Do ta có dãy lặp sau x k+1 = PCm (k) x k , k = 0, 1, m(k ) = k (modM) + Nếu C = ∅ dãy { x k } hội tụ tới điểm thuộc C Còn C = ∅ dãy { x k } không hội tụ Vì tích toán tử PCm toán tử không giãn trung bình nên dãy { x jM+m j = 1, 2, } hội tụ yếu với m cố định tới điểm giới hạn với điều kiện toán tử tích có điểm bất động 3.3 Kĩ thuật khôi phục đại số xử lý ảnh Trong kỹ thuật chụp cắt lớp máy tính (computerized tomography), người ta cần khôi phục lại ảnh gốc từ hình chiếu theo số phương khác Ở dạng rời rạc, toán đặt sau: 45 Giả sử ảnh gốc biểu diễn dạng véctơ x R N ta biết hình chiếu theo hướng am ∈ R N , m = 1, , M, cho trước, tức biết am , x = bm Cần khôi phục lại ảnh gốc x Khi đó, ta có toán tìm nghiệm phương trình tuyến tính Ax = b, A ma trận thực cấp M × N b T = (b1 , , b M ) Với m = 1, , M ta đặt am cột thứ m ma trận A T Khi bm = am , x Đặt Cm = {w am , w = bm } Ta có PCm x = x + (bm − am , x ) am Thuật toán song song hay gọi phương pháp Cimmino có dạng x k +1 = x k + M bm − a m , x k ∑ M m =1 am x k +1 = x k + T A (b − Ax k ) M Thuật toán có dạng sau x k + = x k + bm ( k ) − a m ( k ) , x k am(k) , với m(k) = k mod ( M) + Phương pháp Kaczmarz đề xuất sau phát minh lại tên gọi kĩ thuật khôi phục đại số (ART) Khi phương trình Ax = b có nhiều nghiệm hai phương pháp hội tụ đến nghiệm gần với vectơ ban đầu x0 3.4 Phương pháp ngoại suy tín hiệu với dải tần hữu hạn Trong kỹ thuật xử lý tín hiệu, người ta nói tín hiệu f ∈ L2 (R) có giải tần bị chặn biến đổi Fourier khác không đoạn hữu hạn miền tần số Nếu biết f (t) đoạn hữu hạn đó, ta xác định biến đổi Fourier F (ω ) xác định f (t) đoạn đó, nói cách khác ta ngoại suy hàm f (t) "cửa sổ" Tuy nhiên thực tế, người ta biết giá trị f (t) số điểm "lấy mẫu" rời rạc tn = a + n∆ 46 Cho f (t) F (ω ) cặp biến đổi Fourier, với t ω biến thực, +∞ f (t)eitw dt, F (ω ) = −∞ f (t) = 2π +∞ F (ω )e−itω dω −∞ Giả sử F (ω ) = 0, với |ω | > Ω, Ω số dương Hàm f (t) gọi Ω−giải tần hữu hạn Tiếp theo, giả sử F (ω ) hàm có giá đoạn [−Ω, Ω], với Ω < π Dãy hệ số Fourier F kí hiệu f (n), n ∈ { M, M + 1, , N } Xét hàm G (Ω) hàm ΩG (ω ) = |ω | ≤ Ω, ngược lại G (ω ) Xét dãy hệ số Fourier g = g(n) đặt Dg có thứ tự hạng tử g(n) n ∈ { M, M + 1, , N } nhận giá trị trường hợp lại Đặt F g = G toán tử cho +∞ G (ω ) = ∑ g(n)einω , n=−∞ với ω ∈ (−π, π ) Xét không gian Hilbert H = L2 (−π, π ), C1 = L2 (Ω, Ω) C2 tập tất phần tử G (ω ) H, với hệ số Fourier thỏa mãn g(n) = f (n), n ∈ { M, M + 1, , N } Khi đó, toán tử chiếu mêtric hàm G (ω ) ∈ H lên C1 ΩG (ω ) toán tử chiếu G (ω ) lên C2 F ( D f + ( I − D )F − 1G ) Xét dãy lặp F0 (ω ) = 0, F k+1 = ΩF ( D f + ( I − D )F −1 F k ), với ω ∈ (−π, π ) vớik = 0, 1, Để thực thuật toán, ta cần tính toán phần tử dãy ( I − D )F −1 F k với số nguyên dương n ∈ / { M, M + 1, , N } Hơn ta có F k +1 − F k = Ω F D ( f − f k ) = Ω F a k , với D ( f − f k ) = ak = 0, n ∈ {0, 1, , M − 1} Do F0 = nên ta có F k = ΩF bk , với bk thỏa mãn bk (n) = 0, n ∈ / { M, M + 1, , N } Từ đó, ta suy giới hạn F ∞ có dạng N F ∞ (ω ) = Ω ∑ n= M 47 cn einω , với hệ số cn thích hợp Hệ số cn xác định cách giải hệ phương trình tuyến tính N f (m) = ∑ n= M cn sin Ω(m − n) , π (m − n) với m = M, , N 3.5 Bài toán chấp nhận tách Bài toán chấp nhận tách có nhiều ứng dụng khoa học công nghệ, mà điển hình phương pháp xạ trị với cường độ thay đổi IRMT (Intensively Modulated Radiation Therapy) Đây phương pháp xạ trị tiên tiến với độ xác cao Theo phương pháp IRMT, tia xạ với cường độ khác chiếu vào thể người bệnh, nhằm cung cấp đủ liều xạ tới vùng cần chữa trị (planned target volumes = PTV) hạn chế lượng xạ tới vùng khác (organs at risk = OAR) Giả sử toàn vùng xạ trị bệnh nhân với cấu trúc giải phẫu T gồm PTV OAR, chia thành I đơn vị thể tích, viết tắt voxels (volume elements) Kí hiệu {St } tập hợp voxels cấu trúc t ∈ T Bức xạ phát cách độc lập từ J-tia Gọi R J không gian cường độ chiếu xạ (intensity space) R I không gian liều lượng phóng xạ hấp thụ (dose space) Kí hiệu x = ( x j ) jJ=1 véctơ cường độ chiếu xạ, D = (dij ), dij lượng chiếu xạ hấp thụ voxel i đơn vị cường độ từ tia thứ j Gọi h = (hi )iI=1 véctơ liều lượng, thành phần hi tổng lượng chiếu xạ hấp thụ voxel i, tức hi = ∑ jJ=1 dij x j , hay h = Dx Đặt Hm , m = 1, , M, tập ràng buộc thứ m liều lượng, ví dụ lượng chiếu xạ không vượt mức tối đa không thấp mức tối thiểu, Xn , n = 1, , N, tập ràng buộc thứ n cường độ Ngoài ra, đặt X+ := { x ∈ R J : x j ≥ ∀ j ∈ J } Bài toán đặt tìm véctơ cường độ chiếu xạ x ∗ thỏa mãn điều kiện sau M x ∗ ∈ X+ ∩ (∩nN=1 xn ); h∗ := Dx ∗ ∈ ∩m =1 Hm Đây ví dụ thực tế toán chấp nhận tách đa tập Xét toán chấp nhận tách tổng quát: Tìm điểm c ∈ C với Ac ∈ Q tồn tại, A ma trận thực cấp M × N, C Q tập lồi đóng, khác rỗng RN 48 RM tương ứng Để giải toán, ta xét thuật toán CQ với bước lặp x k+1 = PC ( x k − γA T (( I − PQ ) Ax k ), với γ ∈ 0, ρ( A2T A) Thuật toán CQ hội tụ tới nghiệm toán chấp nhận tách với phần tử đầu vào x0 , toán có nghiệm Còn trường hợp toán nghiệm thuật toán hội tụ đến điểm cực tiểu hàm f (x) = PQ Ax − Ax Do thuật toán CQ giống với phương pháp chiếu gradient giải toán tối ưu có ràng buộc Hàm f ( x ) hàm lồi, khả vi R N đạo hàm ∇ f ( x ) = A T ( I − PQ ) Ax Mệnh đề 3.1 Toán tử đạo hàm ∇ f λ−liên tục Lipschitz với λ = ρ( A T A) Do ∇ f ν−đơn điệu mạnh ngược với ν = λ1 Chứng minh Ta có ∇ f ( x ) − ∇ f (y) = A T ( I − PQ ) Ax − A T ( I − PQ ) Ay ≤λ ( I − PQ ) Ax − ( I − PQ ) Ay =λ Ax − Ay 2 + λ PQ Ax − PQ Ay − 2λ PQ Ax − PQ Ay, Ax − Ay Do PQ toán tử không giãn vững nên PQ Ax − PQ Ay, Ax − Ay ≥ PQ Ax − PQ Ay Từ suy ∇ f ( x ) − ∇ f (y) ≤λ Ax − Ay ≤λ Ax − Ay ≤ λ2 x − y 49 2 − PQ Ax − PQ Ay Nếu γ ∈ 0, λ2 toán tử B := PC ( I − γA T ( I − PQ )) toán tử không giãn trung bình Khi B có điểm bất động quỹ đạo dãy { Bk x } hội tụ yếu đến điểm bất động B Giả sử z = PC (z − γA T ( I − PQ ) Az) Do đó, với điểm c ∈ C c − z, z − (z − γA T ( I − PQ ) Az ≥ hay c − z, A T ( I − PQ ) Az ≥ Suy z điểm cực tiểu f ( x ) C Thuật toán CQ sử dụng tham số γ ∈ (0, L2 ), với L giá trị riêng lớn ma trận A T A Việc lựa chọn giá trị tham số tốt nói chung không dễ dàng Vì người ta đề xuất phương pháp lặp thích nghi không sử dụng giá trị L 50 Kết luận Trong luận văn tác giả tìm hiểu trình bày số lớp toán tử không giãn với thủ tục, dãy lặp không gian Hilbert thực giải nhiều toán ứng dụng thực tế chụp cắt lớp ảnh thuộc ngành y học hay kỹ thuật khôi phục tín hiệu có giải tần bị chặn Như vậy, đóng góp luận văn gồm: Trình bày tổng quan toán tử không giãn, toán tử không giãn vững, toán tử đơn điệu mạnh ngược Tìm hiểu toán tử không giãn trung bình, toán tử chiếu mêtric hai phương pháp kết hợp cho hội tụ mạnh phép lặp phương pháp lai ghép phương pháp xấp xỉ gắn kết Một số ứng dụng thực tế lớp toán tử: toán tối ưu có ràng buộc, toán chấp nhận lồi, toán chấp nhận tách phương pháp xạ trị tiên tiến với độ xác cao, kĩ thuật khôi phục đại số xử lý ảnh, toán khôi phục tín hiệu có giải tần bị chặn kĩ thuật xử lý tín hiệu 51 Tài liệu tham khảo [1] Charles Byrne, A unified treatment of some iterative algorithms in signal processing and image reconstruction, Institute of Physics Publishing, 20 (2003), pp 103-120 [2] Kazuhide Nakajo and Wataru Takahashi, Strong convergence theorems for nonexpansive mappings and nonexpansive semigroups, J Math Anal Appl., 279 (2001), pp 372 - 379 [3] Hong-Kun Xu, Viscosity approximation methods for nonexpansive mappings, J Math Anal Appl., 298 (2004), pp 279-291 [4] Patrick L Combettes, Isao Yamadab, Compositions and convex combinations of averaged nonexpansive operators, J Math Anal Appl., 425 (2015), pp 55-70 [5] J B Baillon and G Haddad, Quelques propriétés des opératerus angle-bornés et ncycliquement monotones, Isarel J Math., 26 (1977), pp 137-150 [6] B Halpern, Fixed points of nonexpanding maps, Bull Amer Math Soc., 73 (1967), pp 957-961 [7] P L Lions, Approximation de points fixes de contractions, C R Acad Sci Sèr A-B Paris, 284 (1977), pp 1357-1359 [8] A Moudafi, Viscosity approximation methods for fixed-points problem, J Math Anal Appl., 241 (2000), pp 46-55 [9] Z Opial, Weak convergence of the sequence of successive approximations for nonexpansive mappings, Bull Amer Math Soc., 73 (1967), pp 591-597 52 [10] Andrzej Cegeielski, Iterative Methods for Fixed Pointed Problems in Hilbert Spaces, Springer (2014) [11] Heinz H Bauschke, Patrick L Combettes, Convex Analysis and Monotone Operator Theory in Hilbert Spaces, Springer (2010) 53 ... chất toán tử không giãn, toán tử tựa không giãn, toán tử không giãn vững, toán tử đơn điệu mạnh ngược mối quan hệ loại toán tử Một số tính chất tập điểm bất động toán tử không giãn họ toán tử không. .. pháp thiết lập mở rộng không cho toán tử chiếu trực giao siêu phẳng mà cho toán tử không giãn, toán tử tựa không giãn, toán tử không giãn vững, toán tử không giãn trung bình, vv Trong luận văn... tương ứng X tập khác rỗng H Định nghĩa 2.1 Cho T : X → H toán tử không giãn, α ∈ (0, 1) Toán tử T gọi toán tử không giãn trung bình với hệ số α hay α không giãn trung bình tồn toán tử không giãn