Một số phương pháp giải xấp xỉ phương trình tích phân tuyến tính fredholm

B Ộ G IÁ O D Ụ C V À Đ À O T Ạ O TRƯỜNG ĐẠI HỌC s PH ẠM HÀ NỘI ĐỖ TH Ị T H U H Ư Ờ N G M ỘT SỐ PH Ư Ơ N G P H Á P GIẢI X A P x ỉ PH Ư Ơ N G T R ÌN H TÍCH P H Â N T U Y E N t í n h FREDHOLM L U Ậ N V Ă N TH Ạ C s ĩ T O Á N HỌC H N ội-2016 B Ộ G IÁ O D Ụ C V À Đ À O T Ạ O TRƯỜNG ĐẠI HỌC s PH ẠM HÀ NỘI ĐỖ TH Ị T H U H Ư Ờ N G M ỘT SỐ PH Ư Ơ N G P H Á P GIẢI X A P x ỉ PH Ư Ơ N G T R ÌN H TÍCH P H Â N T U Y E N t í n h FREDHOLM L U Ậ N V Ă N TH Ạ C s ĩ T O Á N HỌC Chuyên ngành: Toán giải tích Mã số: 60 46 01 02 Người hướng dẫn khoa học: P G S TS K huất Văn N inh H N ội-2016 Lời cảm ơn Trước trình bày nội dung luận văn, em xin bày tỏ lòng biết ơn sâu sắc tới PGS TS Khuất Văn Ninh, người tận tình hướng dẫn để em hoàn thành luận văn Em xin bày tỏ lòng biết ơn chân thành tới toàn thể thầy cô giáo Phòng Sau đại học, Trường đại học sư phạm Hà Nội dạy bảo em tận tình suốt trình học tập Nhân dịp em xin gửi lời cảm ơn chân thành tới gia đình, bạn bè bên em, cổ vũ, động viên, giúp đỡ em suốt trình học tập hoàn thành luận văn Hà Nội, tháng 06 năm 2016 Học viên Đ ỗ T hị Thu Hường Lời cam đoan Tôi xin cam đoan, hướng dẫn PGS TS Khuất Văn Ninh, luận văn Thạc sĩ chuyên ngành Toán giải tích với đề tài “Một số phương pháp giải xấp xỉ phương trình tích phẫn tuyến tính Fredholm” hoàn thành nhận thức thân tác giả Trong trình nghiên cứu thực luận văn, tác giả kế thừa thành tựu nhà khoa học với trân trọng biết ơn Hà Nội, tháng 06 năm 2016 Tác giả Đ ỗ T hị Thu Hường lĩ M ục lục M ỏ đầu 1 K iến th ứ c ch u ấ n bị 1.1 Không gian metric 1.2 Không gian đinh chuẩn 1.3 Không gian Hilbert 1.4 Tống quan phương trình tích phân tuyến tính Fredholm 10 1.4.1 Phương trình tích phân tuyến tính Fredholm loại 10 1.4.2 Phương trình tích phân tuyến tính Fredholm loại 11 1.5 Công thức cầu phương 11 1.6 Công thức hình thang 12 1.7 Công thức Simpson 13 1.8 Tích phân phụ thuộc tham số 14 M ộ t số p h n g p h p giải p h n g tr ìn h tíc h p h â n tu y ế n tín h F re d h o lm loại 16 2.1 Các phương pháp giải tích 16 2.1.1 16 Phương pháp nhân suy biến iii 2.2 2.1.2 Phương pháp phân tích Adomian 24 2.1.3 Phương pháp xấp xỉ liên tiếp 28 2.1.4 Phương pháp chuỗi lũy thừa 35 Phương DhấD số 39 2.2.1 39 Phương pháp cầu phương M ộ t số p h n g p h p giải p h n g tr ìn h tíc h p h â n tu y ế n tín h F re d h o lm loại 48 3.1 Phương pháp quy hóa| 48 3.2 Phương pháp nhiễu đồng luân 51 K ế t lu ậ n 60 T ài liệu th a m k h ảo 61 M đầu Lý chọn đề tài Trong toán học phương trình tích phân phương trình hàm số chưa biết xuất dấu tích phân Phương trình tích phân xem công cụ toán học hữu ích nhiều lĩnh vực nên quan tâm nghiên cứu theo nhiều khía cạnh khác tồn nghiệm, xấp xỉ, tính chỉnh hay không chỉnh, nghiệm chỉnh hóa Nó có ứng dụng rộng rãi không toán học mà nhiều ngành khoa học khác, ví dụ nghiên cứu phương trình tích phân với điều kiện xác định để giải số vấn đề vật lý mà phương trình vi phân mô tả tượng khuếch tán, tượng truyền Vì việc nghiên cứu giải phương trình tích phân đóng vai trò quan trọng lý thuyết toán học Một số toán (rất ít) có lời giải dạng giải tích, thường gọi “nghiệm đóng”, đa phần không Vì để tìm câu trả lời cho toán, ta cần biết số phương pháp cho nghiệm không dạng giải tích mà dạng bảng số Kiểu phương trình tích phân phương trình tích phân Fredholm Dưới hướng dẫn PGS TS Khuất Văn Ninh với mong muốn nghiên cứu tìm hiểu sâu sắc phương trình tích phân tuyến tính Fredholm, chọn đề tài “M ộ t số p h n g p h p g iả i xấ p x ỉ p h n g tr ìn h tíc h p h â n tu y ế n tín h F red h o lm ” làm luận văn thạc sĩ M ục đích nghiên cứu • Nghiên cứu phương trình tích phân tuyến tính Fredholm hệ thống lại số phương pháp giải phương trình tích phân tuyến tính Fredholm • Vận dụng phương pháp giải phương trình tích phân tuyến tính Fredholm để giải số phương trình tích phân tuyến tính Fredholm cụ thể N hiệm vụ nghiên cứu • Tìm hiểu khái quát khái niệm giải tích hàm, phương trình tích phân tuyến tính Fredholm • Làm rõ số phương pháp giải phương trình tích phân tuyến tính Fredholm • Minh họa qua ví dụ, tập cụ thể Đ ối tượng phạm vi nghiên cứu • Đối tượng nghiên cứu: Phương trình tích phân tuyến tính Fredholm • Phạm vi nghiên cứu: Nghiên cứu số phương pháp giải phương trình tích phân tuyến tính Fredholm Phương pháp nghiên cứu • Phương pháp lý luận: Trước tiên đọc nghiên cứu tài liệu, giáo trình có liên quan đến phương trình tích phân tuyến tính Fredholm số phương pháp giải phương trình tích phân tuyến tính Fredholm Sau phân hóa, hệ thống kiến thức • Một số phương pháp công cụ giải tích hàm, phương pháp giải phương trình tích phân, phương pháp số Các kết dự kiến Luận văn nhằm hệ thống số phương pháp giải phương trình tích phân tuyến tính Fredholm ứng dụng phương pháp để giải số phương trình tích phân tuyến tính Fredholm cụ thể Aj = ^ (íj) - tị x(tị) Xị 0 1 ,1 1 ,0 0 ,0 0 ,2 ,9 9 ,0 0 ,3 ,9 9 0 ,0 0 6 ,4 ,9 9 ,0 0 ,5 ,9 9 ,0 0 ,6 ,9 9 ,0 0 ,7 ,9 9 0 ,0 0 ,8 ,9 9 ,0 0 ,9 ,9 9 ,0 0 1 10 ,0 1 ,0 3 ,0 3 Xị\ Kết luận: Ta thấy với n lớn độ xác nghiệm cao, tức sai số nhỏ 47 Chương M ột số phương pháp giải phương trìn h tích phân tu y ến tín h Fredholm loại 3.1 Phương pháp quy hóa Đ ịn h n g h ĩa 3.1 Giả sử A toán tử H lim H (không gian Hilbert) Nếu —J7—^77— = +oo A gọi toán tử í i||í o + o o ||/ / || Phương pháp quy hóa bao gồm biến đổi toán đặt không chỉnh sang toán đặt chỉnh Phương pháp quy hóa biến phương trình tích phân tuyến tính Fredholm loại (3.1) a thành phương trình tích phân tuyến tính Fredholm loại (3.2) a Nghiệm Up hội tụ tới nghiệm u(x) ịi —»• Ta có bổ đề sau 48 B ổ đề 3.1 Giả sử toán tử tích phân (]3.lị) liên tục toán tử không gian Hilbert, f ( x ) , u ( x ) UfJi(x) xấc định i Ịwm| bị chặn độc lập với ịi ii ịuự —u(x) I —>0 ¡1 —>0 Phương trình tích phân loại (3.2) giải phương pháp mà ta trình bày trước đó, ta tìm ufl(x) Nghiệm u(x) (3.1) thu u(x) — lim ufl(x) ịx— ^ Sau ta trình bày số ví dụ minh họa, dùng phương pháp quy hóa để biến phương trình loại thành loại Phương trình loại giải phương pháp mà ta dùng trước chương V í d ụ 3.1 Dùng phương pháp quy hóa phương pháp nhăn suy biến để giải phương trình tích phần tuyến tính Fredholm loại r /4 e^u (t)d t (3.3) = l Dùng phương pháp quy hóa, phương trình (3.3) trở thành '1 /4 u, u - ỉ ? - ì ^x —t U n(t)dt (3.4) l Ta phương trình tích phân tuyến tính Fredholm loại (3.4) Phương trình giải phương pháp nhân suy biến Phương trình (3.4) viết lại sau u^ x) = (i 49 e"’ (3.5) ãó a = fo^4e tutl(t)dt Để xác định a ta (3.4) vào (3.5), lấy tích phân vế ta Oi = 1 + 4//' Dẫn đến u^x) = + 4// Nghiệm u ( x ) phương trình u ( x ) = lim ưM(a;) = ex /¿-»0 V í d ụ 3.2 Dùng phương pháp quy hóa phương phấp nhẫn suy biến để giải phương trình Fredholm loại + = ị (4tex + 3)u(t)dt (3.6) J0 Dùng phương pháp quy hóa, phương trình (3.6) trở thành u, {x) = -e* + - í (4tex + 3)ưM(í)cỉí F ự l1 dd 0 (3.7) Ta phương trình Fredholm loại (3.7) Phương trình giải phương pháp nhân suy biến Phương trình (3.7) viết lại sau ịà u,{x) = Q - e + p ịi 3(3' - (3.8) a = / tu ^d t, p = / u^{t)dt 50 (3.9) Để xác định a , /3 ta (3.8) vào (3.9), lấy tích phân vế ta _ 3(e —3 —ịì) 2(6e — 18 —7fi — ụ?) ’ —2(e + + ịiè) 6e — 18 —7ịi — ịi2 Thế kết vào (|3.8)cho ta nghiệm gần Ui Ảx ) (1 + ụ)ex + (7 - 3e + n) 6(3 - e) + (7fi + ịi2) Nghiệm u(x) phương trình (3.6) thu 1• , X , - 3e u{x) = lim u ^ x ) = — - -.e + — ■ mT õ 6(3 - e) 6(3 - e) Dễ dàng nghiệm khác phương trình u (a;) = X2 Như nói trước đó, phương trình tích phân tuyến tính Fredholm loại toán đặt không chỉnh nên không tồn nghiệm, tồn nghiệm không 3.2 Phương pháp nhiễu đồng luân Sau ta trình bày phương pháp nhiễu đồng luân để giải phương trình tích phân Fredholm loại dạng f ( x ) = í K (x,t)u (t)d t (3.10) J a Bây ta định nghĩa phép toán Đặt L(u)(x) := f ( x ) — í K (x,t)u (t)d t J a 51 (3.11) Khi phương trình (3.10) có dạng L(u)(x) = Ta xác định đồng luân lồi có dạng H(u,p) = (1 —p)u(x) + pL(u)(x) = (3.12) Tham số nhúng p tăng đơn điệu từ đến Phương pháp nhiễu đồng luân dùng khai triển 00 «=Ệ pnuri (3.13) n= 00 Đặt u(x) = lim pnun p ^ n=0 Chuỗi (3.13) hội tụ tới nghiệm nị nghiệm tồn Thế (3.13) vào (3.12) làm ta thu quan hệ truy hồi u0(x) = 0, Ui(x) = f( x ) , un+i(x) K ( x , t)un(t)dt, n > (3.14) Nếu hạt nhân tách tức K ( x , t ) = g(x)h(t) điều kiện sau |1 - í K(t,t)dt\ < (3.15) Ja điều kiện đủ để dãy {un} hội tụ Ta xét trường hợp K ( x , t ) = g(x).h(t) Phương pháp nhiễu đồng luân dùng để giải phương trình tích phân Predholm loại sau 52 V í dụ 3.3 Dùng phương pháp nhiễu đồng luân để giải phương trình tích phẫn tuyến tính Fredholm loại sau f 1/3 —ex = / ex~tu(t)dt Jữ biết [ /3 1- / Jữ K (t, t)dt -!■ Dùng quan hệ truy hồi (3.14) ta tìm u {x) = 0, M x ) = ị e x, r /3 un+1 (x) = un(x) - J ex~tun(t)dt, n > Dẫn tới thành phần w (z) 0, Ml (a:) - e x, ’ u2{x) Ui{x) í* /3 / ex tUi{ì)dt = -e /3 u3(x) u {x) L ex tu (t)dt = — 27“ ’ Uị(x) u (x) j eX ex~tu (t)dt _= — 81 ’ /* /3 53 làm Do nghiệm gần cho „ + „ + ^ r+ ^ T + 27 81 1 '3 ' 1- u (x ) eX (3.16) Sự hội tụ dẫn đến nghiệm u {x) ex V í d ụ 3.4 Dùng phương phấp nhiễu đồng luân để giải phương trình tích phân tuyến tính Fredholm loại sau /•1 /4 _ = / et~xu(t)dt 4e J0 biết /•1 /4 1— / Jữ K ( t , t)dt Dùng quan hệ truy hồi (3.14) ta tìm u0(:r) = 0, M x) = e /•1 /4 un+i(x) = un(x) - J et~xun(t)dt, n > 54 Dẫn tới thành phần u, U0{X) = u l(x ) = /•1 /4 u2(x) = u^x) - J et~xu1(t)dt = Y^e~x, u {x) = u {x) - j Z*1/4 g et~xu (t)dt = ~^e~x, {x) = u {x) - é ~ xu {t)dt = ^ e " * , Uị làm Do nghiệm gần cho , / w \4 _ 1 = e~x - - 16 64 27 256 \ ) - ĩ = e Sự hội tụ dẫn đến nghiệm u(x) = e~x V í d ụ 3.5 Dùng phương phấp nhiễu đồng luân để giải phương trình tích phân tuyến tính Fredholm loại sau X I xtu(t)dt Jo biết — í K ( t,t) d t = Dẫn tới thành phần uữ{x) = 0, Ui{x) — X, r1 u2(x) = Uị(x) — / xtUị(t)dt = -X, J0 r1 u3(x) = u2(x) — / xtu2(t)dt — — X, J0 r1 g « 4(2 ) = w3 (a;) — / x tu (t)dt = — a ao 27 làm Do nghiệm cho U( * ) = * í i + ! + ĩ + £ + = X I - -2 = 3x Vậy u(x) = 3x V í d ụ 3.6 Dùng phương pháp nhiễu đồng luân để giải phương trình tích phân tuyến tính Fredholm loại sau = / xtu(t)dt J0 —X 56 biết — ị K( t , t) d t J0 = un+i{x) = un(x) Dẫn tới thành phần u ữ{ x ) = 0, Ui{x) = 5X, r1 u 2{x) u 3(x) = = U ị ( x ) —I J0 u ( x ) —/ U ị(x) — u 3(x) xtUị(t)dt r 1x t u ( t ) d t — / J0 = -X, y = x t u 3(t) d t — 10 27 20 -X, 8Ĩ; làm Do nghiệm gần cho ,, 5_ / u{x) = ^-x + j-x Ị + ^ ^ 9 5 — —X -ị— X - õ 1_ 15 = - X -ị - X 45 + 90 - _5 135 = - X 54 5x x 57 + 27 Sự hội tụ dẫn đến nghiệm u(x) —X = y ’ V í d ụ 3.7 Dùng phương pháp nhiễu đồng luân để giải phương trình tích phân tuyến tính Fredholm loại sau 4X = Ị xtu(t)dt Jữ biết — í K (t,t)d t J0 - Dùng quan hệ truy hồi (3.14) ta tìm u 0(z ) = , Uị(x) = ỉx, un+i{x) = un(x) - xtun(t)dt ,n > í! Dẫn tới thành phần u0(x) = 0, Uị(x) = 4*' u2{x) = Ui{x) Uz{x) = I xtUị(t)dt J0 u2{x) Uị(x) - Uz(x) I xtu2(t)dt J0 — / Jl0 58 6X’ 9X’ xtu 3(t)dt — — X 27 làm Do nghiệm gần cho yJ / V Vậy nghiệm 59 27 \ ) K ết luận Luận văn trình bày số phương pháp giải phương trình tích phân tuyến tính Fredholm loại loại Các kết trình bày luận văn bao gồm: Các phương pháp giải phương trình tích phân tuyến tính Fredholm loại là: Phương pháp nhân suy biến, Phương pháp phân tích Adomian, Phương pháp xấp xỉ liên tiếp, Phương pháp chuỗi lũy thừa ví dụ minh họa Một số phương pháp giải phương trình tích phân tuyến tính Fred­ holm loại bao gồm: Phương pháp quy hóa, Phương pháp nhiễu đồng luân ví dụ minh họa Mặc dù tác giả cố gắng, song kiến thức hạn chế nên luận văn không tránh khỏi thiếu sót, tác giả mong nhận ý kiến đóng góp Thầy, Cô giáo bạn 60 Tài liệu th am khảo [A] Tài liệu tiến g V iệt [1] Phạm Kỳ Anh (2005), Giải tích số, Nhà xuất ĐHQG Hà Nội [2] Phan Đức Chính (1978), Giải tích hàm, Nhà xuất Đại học trung học chuyên nghiệp, Hà Nội [3] Nguyễn Minh Chương, Ya.D.Mamedov, Khuất Văn Ninh (1992), Giải xấp xỉ phương trình toán tử, Nhà xuất khoa học kĩ thuật, Hà Nội [4] Nguyễn Minh Chương, Nguyễn Văn Khải, Khuất Văn Ninh, Nguyễn Văn Tuấn, Nguyễn Tường (2001), Giải tích số, Nhà xuất Giáo dục [5] Nguyễn Phụ Hy (2005), Giải tích hàm, Nhà xuất khoa học kĩ thuật Hà Nội [6] 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 61 [...]... này, chúng tôi trình bày một số kiến thức cơ bản của giải tích hàm, giải tích số như: Không gian metric, không gian định chuẩn, tổng quan về phương trình tích phân tuyến tính Predholm, Những kiến thức này được sử dụng để trình bày các phương pháp giải xấp xỉ phương trình tích phân tuyến tính Predholm Các khái niệm này ta có thể tìm thấy trong [2] và [6] 1.1 Không gian m etric Cho X là một tập tùy ý... (2.10) là (f(s) = —es + es + 1 = 1 23 2.1.2 Phương pháp phân tích A dom ian Phương pháp phân tích Adomian (ADM) là phương pháp tìm ra nghiệm dạng chuỗi vô hạn Thành phần U j,j > 0 dễ dàng tính được nếu biết số hạng không thuần nhất f ( x ) trong phương trình tích phân tuyến tính Fredholm u(x) = f ( x ) + A í K( x, t )u( t ) dt Ja là một đa thức của một hay hai số hạng Tuy nhiên, nếu hàm f ( x ) là tổ... D Khi đó tích phân dy phụ thuộc tham số a là một hàm khả vi Đ ịn h lý 1.7 (Tính khả tích của tích phân phụ thuộc tham số) Nếu f i x , y ) là hàm liên tục trong hình chữ nhật D = [a,b] X [c,đ\ thì ta có công thức I{y)dy = í \ f ( x , y ) d x } dy = Hay là 15 ư f(x,y)dy\dx Chương 2 M ột số phương pháp giải phương trìn h tích phân tu y ến tín h Fredholm loại 2 2.1 2.1.1 Các phương pháp giải tích P h... tr ìn h tíc h p h â n tu y ế n tín h F re d h o lm loại 1 Phương trình tích phân tuyến tính Fredholm loại 1 có dạng (1.3) a trong đó D là tập đóng, bị chặn trên trường số thực, hàm ẩn ip(s) chỉ xuất hiện trong dấu tích phân, hạt nhân K(s, t ) và hàm f ( s ) là các hàm 10 giá trị thực, À là tham số cho trước Phương trình tích phân tuyến tính Fredholm loại 1 là "đặt không chỉnh" nên nó có thể vô nghiệm... toán giải phương trình (2.2) với nhân suy biến gồm các bước sau: Bước 1: Tính các tích phân f i = Ị Bj (t)f(t)dt]aij = Ị B i(t)Aj dt](i,j = l, ,n) Ja Ja Bước 2 Giải hệ phương trình tuyến tính n Cị A ^ ^C j ữ ị j fijT 1,2, , Tỉ 3= 1 Nếu định thức của hệ khác không, thì tồn tại duy nhất Cj, ỉ = 1, , n Bước 3 Nghiệm cần tìm có dạng ip(s) = /(s) + A ciA i(s )ỉ—1 20 V í d ụ 2.1 Giải phương trình tích phân tuyến. .. loại 2 Phương trình tích phân tuyến tính Fredholm loại 2 có dạng f ( s ) = f ( s) + X ị K( s , t ) f ( t ) dt , s e D (1.4) Ja trong đó D là tập đóng, bị chặn trên trường số thực Hàm ẩn f ( s ) xuất hiện bên trong và bên ngoài dấu tích phân, hạt nhân K ( s , t ) và hàm f ( s ) là các hàm giá trị thực, A là tham số Khi f ( s ) = 0 phương trình được gọi là thuần nhất 1.5 Công thức cầu phương Phương pháp. .. phương pháp giải tích P h ư ơ n g p h á p n h â n su y b iế n Xét phương trình tích phân tuyến tính Predholm loại II ip(s) = f ( s ) + X ị K(s,t)íf(t)dt Ja (2.1) Ta xét (2.1) trong trường hợp K ( s , t ) là nhân suy biến hoặc khi K ( s , t ) có thể xấp xỉ bởi một nhân suy biến Đ ịn h n g h ĩa 2.1 Nếu nhân K ( s , t ) của phương trình tích phân được n biểu diễn dưới dạng K ( s , t ) = K n(s, t ) = Ỵ2 Ai(s)Bi(t)(a... ta có hệ Ci — 2c2A + — 2 á c2 = |ciA Ci c2 = 3 _ 4AA2 9 - 12A2 Thay Ci,c2 vào (2.9) ta được , 2As 4A2 v(s)- r = 4 v + r ì ỡ là nghiệm của phương trình đã cho 21 2 / 2 3' +s (A ^4 +J í2< V í dụ 2.2 Giải phương trình tích phân tuyến tính Fredholm loại 2 bằng phương pháp nhân suy biến (f(s) = s + í (st2 + s2t)(p(t)dt (2.10) (|2.10D o (f(s) = s + s ị t2(f(t)dt + s2 í t(f(t)dt J0 Jữ (2.11) Đặt ( 2 12)... quy tắc tính mà có các công thức cầu phương với các đại lượng A k, x k, R n tương ứng 1.6 Công thức hình thang Giả sử ta phải tính tích phân xác định của hàm f ( x ) được cho bằng bảng hoặc biểu thức giải tích nhưng không biết nguyên hàm của nó Do đó không thể áp dụng khái niệm nguyên hàm để tính Còn nếu ta dùng định nghĩa tích phân thì phải thực hiện nhiều các tính toán Trong trường hợp này ta tính gần... > 1 Điều này dẫn đến nghiệm đúng là u(x) = X + sin 2.1.3 ,1 + 1 —— 2 P h ư ơ n g p h á p x ấ p xỉ liên tiế p Cho phương trình Fredholm tuyến tính loại hai u(x) = f ( x ) + \ ị K (x,t)u(t)dt, (2.28) J a trong đó u(x) là hàm ẩn cần được xác định, K ( x , t ) là hạt nhân và A là một tham số Phương pháp xấp xỉ liên tiếp dựa trên quan hệ truy hồi w0(x) = hàm giá trị thực được chọn bất kỳ, fb un+i(x) = f