Phương pháp hiệu chỉnh Browder-Tikhonov cho phương trình phi tuyến không chỉnh loại J-đơn điệu

ĐẠI HỌC QUỐC GIA HÀ NỘI TRƯỜNG ĐẠI HỌC KHOA HỌC TỰ NHIÊN VŨ XUÂN QUỲNH PHƯƠNG PHÁP HIỆU CHỈNH BROWDER TIKHONOV CHO PHƯƠNG TRÌNH PHI TUYẾN KHƠNG CHỈNH LOẠI J - ĐƠN ĐIỆU LUẬN VĂN THẠC SĨ KHOA 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 VŨ XUÂN QUỲNH PHƯƠNG PHÁP HIỆU CHỈNH BROWDER TIKHONOV CHO PHƯƠNG TRÌNH PHI TUYẾN KHÔNG CHỈNH LOẠI J - ĐƠN ĐIỆU Chuyên ngành: Toán ứng dụng Mã số: 60.46.01.12 LUẬN VĂN THẠC SĨ KHOA HỌC NGƯỜI HƯỚNG DẪN KHOA HỌC GS TS Nguyễn Bường Hà Nội - 2015 LỜI CẢM ƠN Trước trình bày nội dung luận văn, tơi xin bày tỏ lòng biết ơn chân thành sâu sắc tới GS TS Nguyễn Bường-Viện Công nghệ thông tin-Viện Hàn lâm Khoa học Công nghệ Việt nam, người tận tình hướng dẫn, bảo tơi hồn thành luận văn Tôi xin gửi lời cảm ơn tới thầy, cô giáo công tác trường Đại học Khoa học Tự nhiên-Đại học Quốc gia Hà nội truyền đạt kiến thức cho tơi suốt q trình học tập trường Cuối cùng, xin gửi lời cảm ơn tới lãnh đạo Viện Công nghệ thông tin, bạn đồng nghiệp gia đình tạo điều kiện thuận lợi để tơi hồn thành luận văn Hà Nội, tháng 10 năm 2015 Học viên Vũ Xuân Quỳnh Mục lục Mở đầu Khái niệm 1.1 Không gian Banach 1.2 Bài tốn đặt khơng chỉnh 1.2.1 Khái niệm tốn đặt khơng chỉnh 1.2.2 Khái niệm thuật toán hiệu chỉnh 10 1.3 Phương trình với tốn tử loại J-đơn điệu 12 1.3.1 Một số khái niệm 12 1.3.2 Phương trình với tốn tử loại J-đơn điệu 16 Phương pháp hiệu chỉnh Browder-Tikhonov 19 2.1 Phương pháp Browder-Tikhonov với toán tử loại J-đơn điệu 19 2.2 Phương pháp hiệu chỉnh Newton-Kantorovich 37 Kết luận 44 Tài liệu tham khảo 45 Mở đầu Trong lớp toán nảy sinh từ khoa học, kỹ thuật nghành kinh tế quốc dân tồn lớp tốn mà nghiệm khơng ổn định với kiện ban đầu Khi kiện ban đầu thay đổi chút phương trình khơng có nghiệm có nghiệm tương ứng lại cách xa nghiệm xác nhiều Người ta nói tốn đặt khơng chỉnh đặt yêu cầu tìm phương pháp giải ổn định tốn Ta xét tốn đặt khơng chỉnh dạng phương trình tốn tử A(x) = f, f ∈ X, (1) A tốn tử từ không gian Banach X vào không gian Banach Y Khi tốn hiệu chỉnh phương pháp cực tiểu phiếm hàm làm trơn Tikhonov h Fδ,α (x) = ||Ah (x) − fδ ||2 + αΩ(x), x ∈ D(Ah ) = D(A), với việc chọn tham số α = α(h, δ) thích hợp, (Ah , fδ ) xấp xỉ (A, f ), α > tham số hiệu chỉnh, Ω(x) phiếm hàm ổn định Tuy nhiên toán phi tuyến việc tìm phần tử cực tiểu phiếm hàm Tikhonov trở nên khó khăn Do để giải toán trường hợp phi tuyến, A : X → X ∗ toán tử đơn điệu, [7] Browder đề xuất dạng phương pháp hiệu chỉnh Tikhonov cách sử dụng tốn tử có tính chất h-liên tục đơn điệu mạnh Tiếp tục tư tưởng này, Alber [3] sử dụng ánh xạ đối ngẫu tổng quát để hiệu chỉnh tốn Để tìm nghiệm cho tốn (1), chúng tơi xem xét phương pháp hiệu chỉnh Browder-Tikhonov có dạng A(x) + α(x − x+ ) = fδ , (2) A : X → X tốn tử loại J-đơn điệu khơng gian Banach X có tính chất xấp xỉ Khi ánh xạ đối ngẫu chuẩn tắc J liên tục yếu theo dãy liên tục mạnh (2) có nghiệm xδα hội tụ tới x0 nghiệm (1) Ta hội tụ J tính liên tục yếu theo dãy bổ sung thêm hai điều kiện ||A(x) − A(x0 ) − J ∗ A (x0 )∗ J(x − x0 )|| ≤ τ ||A(x) − A(x0 )||, (3) x ∈ X , τ > 0, J ∗ ánh xạ đối ngẫu chuẩn tắc X ∗ , x0 nghiệm (1) tồn z ∈ X cho A (x0 )z = x+ − x0 (4) Cuối cùng, J không liên tục yếu theo dãy khơng thỏa mãn hai điều kiện (3), (4) nghiệm hiệu chỉnh phương pháp hội tụ tới nghiệm toán Nội dung luận văn trình bày hai chương Chương giới thiệu tốn đặt khơng chỉnh, phương trình với tốn tử loại J-đơn điệu số khái niệm dùng tồn luận văn Chương trình bày phương pháp hiệu chỉnh Browder-Tikhonov cho phương trình phi tuyến khơng chỉnh với tốn tử loại J-đơn điệu phương pháp lặp Newton-Kantorovich kết hợp với phương pháp hiệu chỉnh Luận văn chắn tránh khỏi sai sót Em mong nhận góp ý bảo thầy Em xin chân thành cảm ơn! Chương Khái niệm Chương gồm ba mục Mục 1.1 trình bày khái niệm số ví dụ khơng gian Banach Mục 1.2 giới thiệu tốn đặt khơng chỉnh thuật tốn hiệu chỉnh Trong mục 1.3, chúng tơi trình bày số khái niệm giải tích hàm có liên quan tới luận văn phương trình với tốn tử loại J-đơn điệu Các kiến thức tham khảo từ tài liệu [1], [4] [7] 1.1 Không gian Banach Định nghĩa 1.1 Cho (X, d) không gian metric Dãy {xn } ⊂ (X, d) gọi dãy ∀ > ∃ N = N ( ), ∀ m, n ≥ N ⇒ d(xm , xn ) < (X, d) gọi không gian metric đủ, dãy có giới hạn X Định nghĩa 1.2 Cho X khơng gian tuyến tính Ta nói X khơng gian tuyến tính định chuẩn, với x ∈ X xác định số, gọi chuẩn x (kí hiệu ||x||) thỏa mãn ba tiên đề sau: a) Xác định dương: ∀x ∈ X, ||x|| ≥ Đẳng thức xảy x = 0; b) Thuần dương: ∀x ∈ X, ∀λ ∈ R ||λx|| = |λ|.||x||; c) Bất đẳng thức tam giác: ∀x, y ∈ X ||x + y|| ≤ ||x|| + ||y|| Định nghĩa 1.3 Khơng gian Banach khơng gian tuyến tính định chuẩn đầy đủ Ví dụ 1.1 Khơng gian Rn với chuẩn Euclid khoảng cách xác định sau: n |ξi |2 ) , ||x|| = ( i=1 d(x, y) = ||x − y||, với x = (ξ1 , ξ2 , , ξn ) ∈ Rn , y ∈ Rn khơng gian Banach Ví dụ 1.2 Không gian hàm thực liên tục C[a,b] với chuẩn khoảng cách xác định sau: ||x|| = max |x(t)|, a≤t≤b d(x, y) = max |x(t) − y(t)|, a≤t≤b với x(t), y(t) ∈ C[a, b] không gian Banach Ví dụ 1.3 Khơng gian lp ( p ≥ 1), tập dãy số η1 , η2 , , ηn , thỏa mãn ∞ |ηn |p < ∞ n=1 không gian Banach Thật vậy, với x = (η1 , η2 , ), y = (ν1 , ν2 , ) ∈ lp (p ≥ 1) chuẩn khoảng cách xác định sau: ∞ |ηn |p ) p , ||x|| = ( n=1 ∞ |ηn − νn |p ) p ||x − y|| = ( n=1 (m) (m) Giả sử {xm }∞ m=1 dãy Cauchy lp , xm = (η1 , η2 , ) Do với > ∃m0 , ∀m ≥ m0 , ∀r nguyên dương, ta có ||xm − xm+r || ≤ hay ∞ |ηn(m) − ηn(m+r) |p ) p ≤ , ( n=1 (m) (m+r) (m) − ηn | ≤ ∀n, ∀m ≥ m0 , dó ∀n ∃ limm→∞ ηn = ηn0 (m+r) p p1 (m) Như vậy, ( N − η | ) ≤ ∀N Cho r → ∞ ta |η n n n=1 suy |ηn N |ηn(m) − ηn0 |p ) p ≤ ∀N ( n=1 Ta cho N → ∞ ∞ |ηn(m) − ηn0 |p ) p ≤ , ( n=1 (m) (m) − η10 , η2 − η20 , ) ∈ lp Suy x0 = (η10 , η20 , ) = (m) (m) (m) (m) (η1 , η2 , ) − (η1 − η10 , η2 − η20 , ) = xm − z ∈ lp Nên ||z|| = ||xm − x0 || ≤ ∀m ≥ m0 , hay x0 = limm→∞ xm Vậy lp không gian Banach từ z := (η1 Ví dụ 1.4 Khơng gian co tập tất dãy số µ1 , µ2 , , µn , hội tụ tới không gian Banach Trong không gian này, chuẩn khoảng cách xác định bởi: ||x|| = sup |µn |, n d(x, y) = ||x − y||, (k) (k) x, y ∈ co Cho {xk } dãy Cauchy co , xk = (µ1 , µ2 , ) Khi ∀ > ∃k0 , ∀k ≥ k0 , ∀p nguyên dương ||xk − xk+p || ≤ Vì (k) (k) (k+p) ||xk || = supn |µn | nên |µn − µn | ≤ ∀n Do n cố định (k) (k) {µn } dãy Cauchy, tồn µ0n để limk→∞ µn = µ0n Ta suy (k) (k) |µn − µ0n | ≤ ∀n Với k cố định limn→∞ µn = nên ∃n0 ∀n ≥ n0 (k) (k) (k) để |µn | ≤ Từ |µ0n | ≤ |µn | + |µn − µ0n | ≤ ∀n ≥ n0 Do x0 = (µ01 , µ02 , ) ∈ co limk→∞ xk = x0 1.2 1.2.1 Bài toán đặt khơng chỉnh Khái niệm tốn đặt khơng chỉnh Cho phương trình tốn tử A(x) = f, (1.1) A : (X, d) → (Y, ρ)., X, Y khơng gian mêtric Phương trình (1.1) đặt chỉnh nếu: • Với f ∈ Y tồn nghiệm x(f ) ∈ X (1.1); • Nghiệm nhất; • Nghiệm phụ thuộc liên tục vào kiện (f, A) Phương trình (1.1) gọi đặt không chỉnh ba điều kiện khơng thỏa mãn, tức là, • Phương trình (1.1) khơng có nghiệm; • Phương trình (1.1) có nhiều nghiệm; • Phương trình (1.1) có nghiệm x = x(f ) không phụ thuộc liên tục vào kiện tốn Ví dụ 1.5 Cho phương trình tích phân Fredholm loại b K(t, s)x(s)ds = f0 (t), t ∈ [c, d], (1.2) a −∞ < a < b < +∞, −∞ < c < d < +∞, x(s) nghiệm, vế phải f0 (t) hàm số cho trước, nhân dK K(t, s) hàm liên tục Ta giả thiết nghiệm x(s) ∈ C[a, b] với dt khoảng cách hai hàm x1 x2 lớp ρC[a,b] (x1 , x2 ) = max |x1 (s) − x2 (s)| s∈[a,b] A m-J-đơn điệu nên α1 xδα11 − xδα22 , j(xδα11 − xδα22 ) + (α1 − α2 ) xδα22 − x+ , j(xδα11 − xδα22 ) ≤ fδ1 − fδ2 , j(xδα11 − xδα22 ) , từ α1 ||xδα11 −xδα22 ||2 ≤ (M1 +||x+ ||)|α1 −α2 | ||xδα11 −xδα22 ||+(δ1 +δ2 )||xδα11 −xδα22 ||, ta chọn M1 cho ||xδα22 || ≤ M1 Do vậy, ||xδα11 − xδα22 || ≤ (M1 + ||x+ ||) |α1 − α2 | δ1 + δ2 + α1 α1 ✷ Định lý 2.6 Cho X, A f định lý 2.4 cho S = ∅ Giả sử tồn phần tử v ∈ X để x+ − y∗ = A (y∗ )v đạo hàm Frechet A (.) liên tục Lipchitz địa phương hình cầu Br (y∗ ) = {x ∈ X : ||x − y∗ || ≤ ||x+ − y∗ ||} Khi đó, với α > 0, ta có ||xα − y∗ || ≤ 2(2L||v||2 + ||v||)α Chứng minh Đặt F = A (y∗ ), Rα = α(F + αI)−1 B = F Rα , αF − αB = αF − αF Rα = αF (I − Rα ) = αF [(F + αI) − αI](F + αI)−1 = F B Đặt zα = xα − y∗ , ta A(xα ) − A(y∗ + Bv) + α(zα − Bv) = A(y∗ ) − A(y∗ + Bv) + α(zα − Bv) − α(xα − x+ ) = A(y∗ ) − A(y∗ + Bv) + α(x+ − y∗ ) − αBv = A(y∗ ) − A(y∗ + Bv) + αF (v) − αBv = A(y∗ ) − A(y∗ + Bv) + F Bv, 32 Sử dụng tính liên tục Lipchitz A (.) bổ đề 2.1, ta có α||zα − Bv||2 ≤ Axα − A(y∗ + Bv) + α(zα − Bv), j(zα − Bv) ⇔ α||zα − Bv||2 ≤ A(y∗ ) − A(y∗ + Bv) + F Bv, j(zα − Bv) ⇔ ||zα − Bv|| ≤ ||A(y∗ ) − A(y∗ + Bv) + A (y∗ )Bv|| α ⇔ ||zα − Bv|| ≤ ||A (y∗ ) − A (¯ y∗ )||.||Bv||, y¯∗ ∈ Br (y∗ ) α L ⇔ ||zα − Bv|| ≤ ||Bv||2 α L ⇔ ||zα − Bv|| ≤ ||αF (F + αI)−1 ||2 ||v||2 α ⇔ ||zα − Bv|| ≤ 4Lα.||v||2 Mặt khác ||Bv|| ≤ 2α||v||, ||zα || = ||xα − y∗ || ≤ ||zα − Bv|| + ||Bv|| ≤ 2(2L||v||2 + ||v||)α ✷ Định lý 2.7 Cho X, A f định lý 2.6 Giả sử tồn phần tử v ∈ X để x+ − y∗ = A (y∗ )v hai điều kiện sau thỏa mãn: A (.) liên tục Lipchitz định lý 2.6 Tồn số k0 > cho k0 ||x+ − y∗ || < (A (x)−A (y∗ ))ω = A (y∗ )k(x, y∗ , ω), ||k(x, y∗ , ω)|| ≤ k0 ||ω|| ||x−y∗ ||, ∀x, ω ∈ Br˜(y∗ ), r˜ > r + δ/α √ Khi đó, α chọn cho α = O( δ), √ ||xδα − y∗ || ≤ O( δ) Chứng minh Nếu điều kiện định lý 2.6 thỏa mãn với kết định lý 2.4 ta thu ||xδα − y∗ || ≤ ||xδα − xα || + ||xα − y∗ || ≤ 33 δ + 2(2L||v||2 + ||v||)α α √ √ Do vậy, chọn α = O( δ) ||xδα − y∗ || ≤ O( δ) Bây ta xét trường hợp điều kiện sau thỏa mãn Đặt xt = y∗ + t(xα − y∗ ) với t ∈ (0, 1) Vì xα ∈ Br (y∗ ) nên ||xt − y∗ || = t||xα − y∗ || ≤ ||x+ − y∗ ||, suy xt ∈ Br˜(y∗ ) Mặt khác, từ (2.23) (A (xt ) − A (y∗ ))(xα − y∗ )dt A(xα ) − A(y∗ ) = A (y∗ )(xα − y∗ ) + = A (y∗ )(xα − y∗ ) + A (y∗ )k(xt , y∗ , xα − y∗ )dt suy + −α(xα − x ) = A (y∗ )(xα − y∗ ) + A (y∗ )k(xt , y∗ , xα − y∗ )dt Đẳng thức tương đương với + −α(xα −y∗ )−A (y∗ )(xα −y∗ ) = −α(x −y∗ )+ A (y∗ )k(xt , y∗ , xα −y∗ )dt Vì vậy, ta −1 xα − y∗ = α(F + αI) F v − (F + αI)−1 F k(xt , y∗ , xα − y∗ )dt Ngoài ra, với ||(F + αI)−1 F v|| ≤ 2||v||, ta có 1 ||(F + αI)−1 F k(xt , y∗ , xα − y∗ )||dt ≤ k0 ||xt − y∗ || ||xα − y∗ ||dt t||xα − y∗ ||2 dt ≤ 2k0 ≤ k0 ||xα − y∗ ||2 ≤ k0 ||x+ − y∗ || ||xα − y∗ || Từ đó, suy (1 − k0 ||x+ − y∗ ||) ||xα − y∗ || ≤ 2α||v|| Vì 2α||v|| δ + α − k0 ||x+ − y∗ || √ √ Do vậy, α = O( δ) ta thu ||xδα − y∗ || ≤ O( δ) Định lý chứng minh ✷ ||xδα − y∗ || ≤ ||xδα − xα || + ||xα − y∗ || ≤ 34 Trong [13] Ryazantseva xem xét thuật toán điểm gần kề kết hợp với phương pháp hiệu chỉnh Browder-Tikhonov có dạng: ck (A(xk+1 ) + αk xk+1 − fk ) + xk+1 = xk , x0 ∈ X Sau đó, ơng dãy {xk } sinh phương trình hội tụ mạnh tới nghiệm phương trình (2.1) J có tính chất liên tục liên tục yếu theo dãy, dãy {xk } bị chặn có thêm số điều kiện cho ck αk Ta tính bị chặn dãy {xk } sinh phương trình tổng quát sau ck (A(xk+1 ) + αk xk+1 − fk ) + xk+1 − xk = γk (xk − xk−1 ) (2.29) mà khơng cần tính liên tục yếu theo dãy J Bổ đề 2.3 Cho {ak }, {bk }, {ck } ba dãy số thực dương thỏa mãn điều kiện sau: ak+1 ≤ (1 − bk )ak + bk ck , bk < 1; ∞ n=0 bk = +∞, limk→+∞ ck = Khi limk→+∞ ak = Định lý sau nghiệm phương trình (2.29) hội tụ tới y∗ nghiệm (2.24) với điều kiện ràng buộc cho ck , αk γk Định lý 2.8 Cho X, A, f định lý 2.5 cho S = ∅ fk ∈ X thỏa mãn ||fk − f || ≤ δk , k → ∞ Giả sử tham số ck , αk γk thỏa mãn điều kiện sau: C0 > ck > c0 > 0, γk ≥ 0, C0 c0 số dương, ∞ k=1 τk = +∞, ∞ −1 k=1 γk τk ||xk τk = αk ck /(1 + αk ck ), − xk−1 || < +∞, αk > 0, limk→+∞ αk = 0, limk→+∞ Dk /τk = 0, Dk = M1 |αk − αk+1 | δk + δk+1 + αk αk 35 Khi đó, dãy {xk } xác định phương trình (2.29) hội tụ mạnh tới y∗ nghiệm (2.24), k → +∞ Chứng minh Ký hiệu xk nghiệm (2.22), α fδ thay αk fk x+ = Ta viết phương trình (2.29) (2.22) dạng tương đương sau λk (A(xk+1 ) − fk ) + xk+1 = βk xk + βk γk (xk − xk−1 ), λk (A(xk ) − fk ) + xk = βk xk , λk = βk ck βk = 1/(1 + ck αk ) Từ đó, ta suy λk A(xk+1 ) − A(xk ), j(xk+1 − xk ) + xk+1 − xk , j(xk+1 − xk ) = βk xk − xk , j(xk+1 − xk ) + βk γk xk − xk−1 , j(xk+1 − xk ) Do λk > A m-J-đơn điệu, nên ta thu ||xk+1 − xk ||2 ≤ (βk ||xk − xk || + βk γk ||xk − xk−1 ||)||xk+1 − xk || (2.30) Nếu xk+1 = xk theo định lý 2.5, lim xk+1 = lim xk = y∗ k→+∞ k→+∞ Không tổng quát, giả sử xk+1 = xk Từ (2.30) ta suy ||xk+1 − xk || ≤ βk ||xk − xk || + βk γk ||xk − xk−1 || Kết hợp bất đẳng thức với định lý 2.5 dẫn đến ||xk+1 − xk+1 || ≤ ||xk+1 − xk || + ||xk+1 − xk || ≤ βk ||xk − xk || + βk γk ||xk − xk−1 || + Dk = (1 − τk )||xk − xk || + dk , dk = βk γk ||xk − xk−1 || + Dk ≤ γk ||xk − xk−1 || + Dk Theo điều kiện thứ hai định lý ta suy γk τk −1 ||xk − xk−1 || → k → +∞ Kết hợp điều với điều kiện thứ ba định lý, ta có Dk dk ≤ lim (γk τk −1 ||xk − xk−1 || + ) = 0 ≤ lim k→+∞ k→+∞ τk τk 36 Do vậy, dk = k→+∞ τk lim Áp dụng bổ đề 2.3 với ak = ||xk − xk ||, bk = τk , ck = dk , τk lim ||xk − xk || → k→+∞ Vì ||xk − y∗ || → k → +∞, nên từ ||xk − y∗ || ≤ ||xk − xk || + ||xk − y∗ || ta có xk → y∗ k → +∞ Định lý chứng minh ✷ 2.2 Phương pháp hiệu chỉnh Newton-Kantorovich Trong phần tìm nghiệm cho phương trình phi tuyến khơng chỉnh (2.1) với A toán tử m-J-đơn điệu, dựa vào phương pháp hiệu chỉnh Newton-Kantorovich Ta giả sử tập nghiệm (2.1), ký hiệu S khác rỗng, thay cho f, ta biết xấp xỉ fδ thỏa mãn ||fδ − f || ≤ δ → (2.31) Nếu A khơng có tính chất J-đơn điệu mạnh phương trình (2.1) nói chung khơng chỉnh Để giải toán (2.1), ta dùng phương pháp ổn định phương pháp hiệu chỉnh Browder-Tikhonov (2.22) trình bày phần trước Tuy nhiên A tốn tử phi tuyến tốn phi tuyến, khó giải thực tế Để khắc phục điều này, phương pháp ổn định khác phương pháp hiệu chỉnh Newton-Kantorovich không gian Hilbert H Bukushinskii [5] đưa sau z0 ∈ H, A(zn ) + A (zn )(zn+1 − zn ) + αn zn+1 = fδ (2.32) Sau Bakushinskii Smirnova [6] tìm hội tụ (2.32) với quy tắc dừng hậu nghiệm ||A(zN ) − fδ ||2 ≤ τ δ < ||A(zn ) − fδ ||2 , ≤ n < N = N (δ), 37 (2.33) với điều kiện ||A (x)|| ≤ ∀x ∈ H A liên tục L-Lipchitz, τ > số cố định Phương pháp (2.32) Abel Ryazantseva [4] mở rộng từ không gian Hilbert vào không gian Banach với dạng sau z0 ∈ X, A(zn ) + A (zn )(zn+1 − zn ) + αn J s (zn+1 ) = fδ , với điều kiện ||A (x)|| ≤ ϕ(||x||), (2.34) ϕ(t) hàm khơng âm, khơng giảm với t ≥ 0, A ánh xạ đơn điệu từ X vào X ∗ Ta biết A khả vi Gateaux J-đơn điệu A (z)(.) ánh xạ J-đơn điệu cho z cố định thuộc X Hơn nữa, A (z)(.)là tuyến tính liên tục, liên tục yếu Từ đó, phương trình z0 ∈ X, A(zn ) + A (zn )(zn+1 − zn ) + αn (zn+1 − x+ ) = fδ (2.35) có nghiệm zn+1 , cho số nguyên n ≥ 0, với αn , δ > Để giải toán (2.1) với điều kiện (2.31), chứng minh định lý hội tụ mạnh cho phương pháp hiệu chỉnh Newton-Kantorovich (2.35) với quy tắc dừng hậu nghiệm (2.33) điều kiện (2.34) mà khơng có điều kiện nhiễu vế phải, tức ta xem xét phương pháp lặp sau x0 ∈ X, A(xn ) + A (xn )(xn+1 − xn ) + αn (xn+1 − x+ ) = f (2.36) Định lý hội tụ mạnh dãy {xn } Định lý 2.9 Cho X không gian Banach thực, phản xạ lồi chặt với chuẩn khả vi Gateaux cho A ánh xạ m-J-đơn điệu, khả vi Frechet hai lần X với điều kiện (2.34) Giả sử dãy {αn }, số thực d điểm ban đầu x0 thỏa mãn điều kiện sau: {αn } dãy đơn điệu giảm với < αn < tồn σ > thỏa mãn αn+1 ≥ σαn với n = 0, 1, ; ϕ0 ||x0 − xα0 || ≤ q < 1, ϕ0 = ϕ(d + γ), 2σα0 38 số dương γ tìm từ giới hạn 2σα0 ≤ γ, d ≥ 2||x∗ − x+ || + ||x+ ||, ϕ0 (2.37) x∗ nghiệm bất đẳng thức biến phân sau: x∗ ∈ S : x∗ − x+ , J(x∗ − y) ≤ ∀y ∈ S, (2.38) xα0 nghiệm phương trình A(x) + αn (x − x+ ) = f, (2.39) với n = 0; q − q2 ϕ0 d αn − αn+1 ≤ , c = αn2 c1 2σ Khi limn→+∞ ||xn − x∗ || = Chứng minh Gọi xαn nghiệm (2.39) Theo định lý 2.4, ta có dãy {xαn } hội tụ mạnh tới x∗ nghiệm (2.38) n → +∞ Ngoài ra, xαn ∈ Br (x∗ ) nên ||xαn − x+ || ≤ ||xαn − x∗ || + ||x∗ − x+ || ≤ 2||x∗ − x+ || điều dẫn đến ||xαn || ≤ 2||x∗ − x+ || + ||x+ || ≤ d ta suy ||xαn − xαn+1 || ≤ αn − αn+1 d, ∀n = 0, 1, αn (2.40) Từ giả thiết, xαn nghiệm phương trình (2.39) nên A(xαn )+αn (xαn − x+ ) = f Theo cơng thức Taylor A(xαn ) = A(xn ) + A (xn )(xαn − xn ) + A (θn )(xαn − xn )2 , 39 với θn = xαn + θ(xn − xαn ) < θ < Do vậy, A(xn )+A (xn )(xαn −xn )+ A (θn )(xαn −xn )2 +αn (xαn −x+ ) = f (2.41) Hơn nữa, từ (2.36) (2.41) ta suy A (xn )(xn+1 −xαn )− A (θn )(xαn −xn )2 +αn (xn+1 −xαn ), J(xn+1 −xαn ) = Do tính chất J-đơn điệu toán tử A , từ đẳng thức dẫn đến αn ||xn+1 − xαn || ≤ ||A (θn )|| ||xαn − xn ||2 , kết hợp bất đẳng thức với (2.34) ta có ||xn+1 − xαn || ≤ ϕ(rn ) ∆ , 2αn n (2.42) ∆n = ||xn − xαn || rn ≥ ||xαn || + ∆n ≥ ||θn || Từ điều kiện đầu (2.37), ta có σ < 2σαn ≤ γ ∀n ≥ ϕ0 (2.43) Ta thu ∆0 ≤ γ từ điều kiện thứ hai (2.37) với r0 = d + γ Bây ta chứng minh quy nạp bất đẳng thức sau ϕ0 ∆n ≤ q < 1, ∀n ≥ 2σαn (2.44) Với n = điều kiện thứ hai định lý Giả sử (2.44) với n = k rk = d + γ, ta chứng minh với n = k + Thật vậy, rõ ràng từ (2.40) (2.42) ta có ∆k+1 = ||xk+1 − xαk+1 || ≤ ||xk+1 − xαk || + ||xαk − xαk+1 || ϕ(rk ) αk − αk+1 ≤ ∆ + d 2αk k αk Kết sau suy từ bất đẳng thức điều kiện thứ ba định lý ϕ0 ∆k+1 ϕ0 ∆k+1 ϕ0 ≤ ≤ ∆k 2σαk+1 2σ αk 2σαk ≤ q + q − q = q < 40 + (αk − αk+1 )ϕ0 d 2σ αk Như phần chứng minh quy nạp hoàn thành Tiếp theo, từ bất đẳng thức (2.43), ta có ∆k+1 ≤ 2σαk+1 ≤ γ, ϕ0 theo định nghĩa rn , ta lấy rk+1 = d + γ Vì αn → n → +∞, từ (2.44) ta có ∆n → Mặt khác, từ ||xn − x∗ || ≤ ∆n + ||xαn − x∗ || ||xαn − x∗ || → 0, n → ∞, ta suy hội tụ mạnh dãy {xn } tới x∗ ✷ Định lý 2.10 Cho X, A, αn định lý 2.9 với toán tử A liên tục L-Lipchitz điều kiện thứ thỏa mãn Cho τ > (2.33) chọn cho √ 3dL ϕ||z ˜ − xα0 || ≤ q, < q < − , ϕ˜ = ϕ0 + 2L2 /˜ τ , τ˜ = ( τ − 1)2 , 2σα0 τ˜σ (2.45) d xác định định lý 2.9, số dương γ lấy từ giới hạn (2.37) với ϕ0 thay ϕ˜, αn − αn+1 d 2Lσ + ≤ q αn2 τ˜ ϕ˜ ˜τ (2.46) Khi Cho n = 0, 1, ., N (δ), ϕ||z ˜ n − xαn || ≤ q, 2σαn zn nghiệm (2.35) N (δ) chọn (2.33) Dãy {N (δ)} chấp nhận được, tức là, lim ||zN (δ) − y|| = 0, δ→0 y ∈ S Nếu N (δ) → ∞ δ → 0, y = x∗ 41 (2.47) Chứng minh Ta chứng minh (2.47) phương pháp quy nạp Thật vậy, với trường hợp n = bất đẳng thức (2.47) điều kiện (2.45) Lấy n < N (δ) giả sử cho k cho ≤ k ≤ n < N (δ) (2.47) Theo cơng thức Taylor, ta có A(xαn ) = A(zn ) + A (zn )(xαn − zn ) + A (θ˜n )(xαn − zn )2 , với θ˜n = xαn + θ(zn − xαn ) < θ < Do vậy, A(zn ) + A (zn )(xαn − zn ) + A (θ˜n )(xαn − zn )2 + αn (xαn − x+ ) = f Kết hợp điều với (2.35) ta A (zn )(zn+1 − xαn ) − A (θ˜n )(xαn − zn )2 + αn (zn+1 − xαn ) = fδ − f Do A có tính J-đơn điệu nên từ đẳng thức suy ||zn+1 − xαn || ≤ δ ||A (θ˜n )|| ||xαn − zn ||2 + 2αn αn Mặt khác, ||zn+1 − xαn+1 || ≤ ||zn+1 − xαn || + ||xαn − xαn+1 || δ rn ) ˜ αn − αn+1 ˜ n+1 ≤ ϕ(˜ d+ , ∆ ∆n + 2αn αn αn (2.48) ˜ n = ||zn − xα || r˜n ≥ ||xα || + ∆ ˜ n ≥ θ˜n Vì n < N (δ) từ ∆ n n (2.33) nên τ δ ≤ ||A(zn ) − fδ ||2 Ngoài ra, A có tính liên tục L-Lipchitz nên √ τ δ ≤ ||A(zn ) − fδ || ≤ ||A(zn ) − f || + ||f − fδ || ≤ ||A(zn ) − A(xαn ) − αn (xαn − x+ )|| + ||f − fδ || ≤ L||zn − xαn || + αn d + δ, đó, √ τ δ − δ ≤ L||zn − xαn || + αn d Không tổng quát, ta giả sử δ < Như vậy, (L||zn − xαn || + αn d)2 δ≤ , τ˜ 42 điều với (2.48) dẫn đến ˜ n d αn − αn+1 ϕ(˜ rn ) + 2L2 /˜ τ ˜ 2L∆ αn d2 ˜ ∆n+1 ≤ ∆n + + d+ 2αn τ˜ αn τ˜ Nhân hai vế bất phương trình với ϕ˜ > với lưu ý < σ < 1, 2σ αn sau kết hợp với (2.46) ta αn − αn+1 ϕ˜ d2 ϕ˜ 2dL γ˜n + d + τ˜σ αn2 2σ 2σ τ˜ 2dL dϕ˜ αn − αn+1 d ≤ γ˜n2 + γ˜n + ( + ) τ˜σ 2σ αn2 τ˜ 3dL ) ≤ q < 1, ≤ q(q + τ˜σ γ˜n+1 ≤ γ˜n2 + ˜ n /(2σαn ) Tương tự định lý 2.9, ta lấy γ˜n = ϕ˜∆ rn+1 = d + γ Hai trường hợp xảy sau: N (δ) = N0 ∀δ ≤ δ0 Khi đó, từ (2.33) ta có limδ→0 ||A(zN (δ) ) − fδ || = 0, ||A(zN0 ) − f || = hay zN0 nghiệm (2.1) N (δ) → ∞ ∀δ ≤ δ0 Rõ ràng, xαN (δ) → x∗ δ → Mặt khác, từ ˜ n ≤ (2σαn )/ϕ˜ → n → ∞ ∆ ˜ N (δ) + ||xα − x∗ ||, ||zN (δ) − x∗ || ≤ ∆ N (δ) ta suy dãy {zN (δ) } hội tụ mạnh tới x∗ ✷ 43 Kết luận Sau thời gian làm việc hướng dẫn GS TS Nguyễn Bường luận văn hoàn thành Luận văn giới thiệu tốn đặt khơng chỉnh phương trình với tốn tử loại J-đơn điệu Trình bày phương pháp hiệu chỉnh Browder-Tikhonov cho phương trình phi tuyến khơng chỉnh loại J-đơn điệu J có tính liên tục yếu theo dãy khơng có tính chất Cuối cùng, giới thiệu phương pháp hiệu chỉnh lặp Newton-Kantorovich Trong luận văn định lý chứng minh chi tiết Các vấn đề đưa luận văn dựa kết nghiên cứu gần GS.TS Nguyễn Bường đồng nghiệp 44 Tài liệu tham khảo [1] Phạm Kỳ Anh, Nguyễn Bường (2005), Bài tốn đặt khơng chỉnh, NXB ĐHQG Hà Nội [2] Alber Y (1975), On solution by the method of regularization for operator equation of the first kind involving accretive mappings in Banach spaces , Differential Equations SSSR, XI, 2242-2248 [3] Alber Y (1975), On solving nonlinear equations involving monotone operator in Banach space, Sibirian Mathematics Journal, 26, 3-11 [4] Alber Y., Ryazantseva I (2006), Nonlinear ill-posed problems of monotone type, Springer [5] Bakushinskii A (1976), Regularization Algorithm based on the NewtonKantorovich method for solving variational inequalities, Zh Vychisl Mat Fiz SSSR, 16(6) 1397-1404 [6] Bakushinskii, A B., Smirnova, A (2007), Iterative regularization and genaralized discrepancy principle for monotone operator equations, Numer Funct Anal and Optim 28(1-2) 13-25 [7] Browder F E (1966), Existence and approximation of solutions of nonlinear variational inequalities, Proc Nat Acad Sci USA 56(4) 1080-1086 [8] Nguyen Buong (2004), Convergence rates in regularization for nonlinear ill-posed equations under accretive perturbations, Zh Vychist Mat Mat Fiz, 44(3) 397-402 [9] Buong Ng., Phuong Ng.Th.H (2012) , Convergence rates in regularization for nonlinear ill-posed equations involving m-accretive mappings in Banach spaces, Applied Math Sciences, 6(63) 3109-3117 [10] Buong Ng., Phuong Ng.Th.H (2013), Regularization methods for nonlinear ill-posed equations involving m-accretive mappings in Banach spaces, Russian Math., 57(2) 58-64 [11] Nguyen Buong, Nguyen Duong Nguyen, and Nguyen Thi Thu Thuy (2015), Newton-Kantorovich iterative regularization and genaralized discrepancy principle for nonlinear ill-posed equations involving accretive 45 mappings, Russian Math., 59(5) 32-37 [12] Konyagin C.V (1980), On approximative properties of closed sets in Banach spaces and the characteristics of strongly convex spaces, Dokl Akad Nauk SSSR, 251(2) 276-280 [13] Ryazantseva I P (2002), Regularization proximal point algorithm for nonlinear equations of monotone type in Banach space, Zhurn Vychisl Matem i Matem Fiz., 42(9) 1295-1303 [14] Wang J., Li J., Liu Z (2008), Regularization methods for nonlinear ill-posed problems with accretive operators, Acta Math Scientia., 28b(1) 141-150 46 ... Luận văn giới thiệu tốn đặt khơng chỉnh phương trình với tốn tử loại J-đơn điệu Trình bày phương pháp hiệu chỉnh Browder-Tikhonov cho phương trình phi tuyến khơng chỉnh loại J-đơn điệu J có tính... thiệu toán đặt khơng chỉnh, phương trình với tốn tử loại J-đơn điệu số khái niệm dùng toàn luận văn Chương trình bày phương pháp hiệu chỉnh Browder-Tikhonov cho phương trình phi tuyến khơng chỉnh. .. lý 1.8 J-đơn điệu chặt, phương trình tốn tử tương ứng có nghiệm Chương sau chúng tơi trình bày phương pháp hiệu chỉnh BrowderTikhonov cho phương trình phi tuyến với toán tử loại J-đơn điệu 18