Luận án tiến sĩ toán học giải gần đúng một số bài toán biên phi tuyến cho phương trình vi phân cấp bốn

BỘ GIÁO DỤC VÀ ĐÀO TẠO VIỆN HÀN LÂM KHOA HỌC VÀ CÔNG NGHỆ VIỆT NAM HỌC VIỆN KHOA HỌC VÀ CÔNG NGHỆ NGUYỄN THANH HƯỜNG GIẢI GẦN ĐÚNG MỘT SỐ BÀI TỐN BIÊN PHI TUYẾN CHO PHƯƠNG TRÌNH VI PHÂN CẤP BỐN LUẬN ÁN TIẾN SĨ TOÁN HỌC HàNội - 2019 BỘ GIÁO DỤC VÀ ĐÀO TẠO VIỆN HÀN LÂM KHOA HỌC VÀ CÔNG NGHỆ VIỆT NAM HỌC VIỆN KHOA HỌC VÀ CÔNG NGHỆ NGUYỄN THANH HƯỜNG GIẢI GẦN ĐÚNG MỘT SỐ BÀI TOÁN BIÊN PHI TUYẾN CHO PHƯƠNG TRÌNH VI PHÂN CẤP BỐN Chun ngành: Tốn ứng dụng Mãsố: 46 01 12 LUẬN ÁN TIẾN SĨ TOÁN HỌC NGƯỜI HƯỚNG DẪN KHOA HỌC GS TS Đặng Quang Á TS Vũ Vinh Quang HàNội – 2019 LỜI CAM ĐOAN Tôi xin cam đoan công trình nghiên cứu tơi với hướng dẫn khoa học GS TS Đặng Quang Á TS Vũ Vinh Quang Những kết trình bày Luận án mới, trung thực chưa công bố cơng trình khác Các kết thực nghiệm kiểm tra chương trình tơi thiết kế thử nghiệm mơi trường MATLAB, số liệu hồn tồn trung thực Các kết công bố chung cán hướng dẫn đồng tác giả cho phép sử dụng Luận án Nghiên cứu sinh Nguyễn Thanh Hường i LỜI CẢM ƠN Trước hết, em xin bày tỏ lòng biết ơn chân thành sâu sắc tới Thầy hướng dẫn, GS TS Đặng Quang Á TS Vũ Vinh Quang Trong suốt trình học tập, nghiên cứu thực Luận án, Thầy kiên nhẫn, tận tình bảo, dìu dắt giúp đỡ em Chính niềm say mê khoa học, nghiêm khắc khoa học với quan tâm, động viên khích lệ Thầy động lực khiến em không ngừng nỗ lực, cố gắng vượt qua khó khăn, vất vả để hồn thành Luận án Em xin chân thành cảm ơn Thầy Cơ thành viên nhóm Seminar khoa học Phịng Phương pháp Tốn học Cơng nghệ thông tin, Viện Công nghệ Thông tin cán nghiên cứu Những ý kiến nhận xét đóng góp vơ q báu buổi báo cáo thảo luận giúp em hoàn thành tốt Luận án Em xin chân thành cảm ơn sở đào tạo - Viện Công nghệ Thông tin Học viện Khoa học Công nghệ Quý Viện Học viện tạo điều kiện thuận lợi để em hồn thành tốt q trình học tập nghiên cứu Em xin bày tỏ lòng biết ơn đến lãnh đạo Trường Đại học Khoa học - Đại học Thái nguyên, Ban chủ nhiệm Khoa Toán - Tin, bạn bè đồng nghiệp, gia đình người thân ln đồng hành, hỗ trợ, động viên, giúp đỡ em suốt trình học tập, nghiên cứu hoàn thành Luận án Xin chân thành cảm ơn! ii Danh mục chữ viết tắt ký hiệu Tập số thực Tập số thực không âm Tập số phức Không gian Euclide K chiều Khơng gian hàm có đạo hàm cấp k liên tục [a, b] C([0, ∞)) Không gian hàm liên tục [0, ∞) C(R) Không gian hàm liên tục R C([0, 1] × R) Khơng gian hàm liên tục [0, 1] × R C([a, b], K) Khơng gian hàm liên tục f : [a, b] → K C([a, b] × R4 , R) Khơng gian hàm liên tục f : [a, b] × R4 → R Ω Miền giới nội Γ Biên miền Ω Ω Bao đóng miền Ω C(Ω) Khơng gian hàm liên tục Ω C(Ω × R) Khơng gian hàm liên tục Ω × R C (Ω × R, R) Không gian hàm f : Ω × R → R có đạo hàm riêng cấp liên tục Ω × R C (Ω) Khơng gian hàm có đạo hàm riêng đến cấp hai liên tục Ω ∞ C (Γ) Không gian hàm khả vi vô hạn Γ ∞ C (Ω × R × R) Khơng gian hàm khả vi vơ hạn Ω × R × R ∆, ∆2 , ∇ Toán tử Laplace, toán tử song điều hịa, tốn tử Gradient q L (Ω) Khơng gian hàm khả tích bậc q Ω ∞ L (Ω) Không gian hàm bị chặn hầu khắp nơi Ω kxk Chuẩn phần tử x kxk2 Chuẩn không gian L2 phần tử x H (Ω) Khơng gian Sobolev hàm có đạo hàm suy rộng cấp hai thuộc L2 (Ω) H01 (Ω) Không gian Sobolev hàm triệt tiêu biên Ω, có đạo hàm suy rộng cấp thuộc L2 (Ω) O(h) Vô bé bậc cao h R R+ C RK C k [a, b] iii Danh sách hình vẽ Đồ Đồ Đồ Đồ Đồ Đồ Đồ Đồ Đồ Đồ Đồ Đồ Đồ Đồ Đồ Đồ Đồ Đồ Đồ Đồ Đồ Đồ 3.1 Đồ thị sai số e(m) Phương pháp DIM2 (trái) phương pháp nhanh Wang (phải) Ví dụ 3.1 với k = 1, σ = 2π Đồ thị sai số e(m) Phương pháp DIM2 (trái) phương pháp nhanh Wang (phải) Ví dụ 3.1 với k = 1, σ = 12π Đồ thị sai số e(m) Phương pháp DIM2 (trái) phương pháp nhanh Wang (phải) Ví dụ 3.2 với k = 0.1 Đồ thị sai số e(m) Phương pháp DIM2 (trái) phương pháp nhanh Wang (phải) Ví dụ 3.2 với k = 0.4 Đồ thị sai số e(m) tỉ số r(m) Ví dụ 3.3 3.2 3.3 3.4 3.5 thị thị thị thị thị thị thị thị thị thị thị thị thị thị thị thị thị thị thị thị thị thị e(K) Ví dụ 2.1 r(K) Ví dụ 2.1 e(K) Ví dụ 2.2 r(K) Ví dụ 2.2 nghiệm xấp xỉ Ví dụ 2.1.2 e(K) Ví dụ 2.3 r(K) Ví dụ 2.3 nghiệm xấp xỉ Ví dụ 2.3 e(K) Ví dụ 2.4 e(K) Ví dụ 2.5 nghiệm xấp xỉ Ví dụ 2.5 nghiệm xấp xỉ Ví dụ 2.6 nghiệm xấp xỉ Ví dụ 2.7 nghiệm xấp xỉ Ví dụ 2.10 nghiệm xấp xỉ Ví dụ 2.11 nghiệm xấp xỉ Ví dụ 2.15 nghiệm xấp xỉ Ví dụ 2.16 e(K) Ví dụ 2.17 e(K) Ví dụ 2.18 e(K) Ví dụ 2.19 e(K) Ví dụ 2.20 nghiệm xấp xỉ Ví dụ 2.20 2.1 2.2 2.3 2.4 2.5 2.6 2.7 2.8 2.9 2.10 2.11 2.12 2.13 2.14 2.15 2.16 2.17 2.18 2.19 2.20 2.21 2.22 iv 46 47 48 48 49 49 50 50 59 61 61 63 64 66 69 75 76 84 86 86 87 87 111 111 112 113 113 3.6 3.7 3.8 3.9 3.10 3.11 3.12 3.13 3.14 3.15 Đồ Đồ Đồ Đồ Đồ Đồ Đồ Đồ Đồ Đồ thị thị thị thị thị thị thị thị thị thị nghiệm xấp xỉ Ví dụ 3.3 sai số e(m) tỉ số r(m) Ví dụ nghiệm xấp xỉ Ví dụ 3.4 e(m) Ví dụ 3.5 e(m) Ví dụ 3.6 nghiệm xấp xỉ Ví dụ 3.6 e(m) Ví dụ 3.7 nghiệm xấp xỉ Ví dụ 3.7 e(m) Ví dụ 3.8 nghiệm xấp xỉ Ví dụ 3.8 v 3.4 114 114 114 123 124 124 124 125 126 126 Danh sách bảng 2.1 2.2 2.3 2.4 2.5 2.6 2.7 2.8 2.9 2.10 2.11 2.12 2.13 Sự Sự Sự Sự Sự Sự Sự Sự Sự Sự Sự Sự Sự hội hội hội hội hội hội hội hội hội hội hội hội hội tụ tụ tụ tụ tụ tụ tụ tụ tụ tụ tụ tụ tụ trong trong trong trong trong trong Ví Ví Ví Ví Ví Ví Ví Ví Ví Ví Ví Ví Ví dụ dụ dụ dụ dụ dụ dụ dụ dụ dụ dụ dụ dụ 2.1 2.4 2.13 với xấp xỉ đầu v0 = 2.13 với xấp xỉ đầu v0 = 2.14 với xấp xỉ đầu v0 = 2.15 với xấp xỉ đầu v0 = 2.16 với xấp xỉ đầu v0 = 2.17 2.21 với xấp xỉ đầu u0 = 0, v0 = 2.22 với xấp xỉ đầu u0 = 0, v0 = 2.23 với xấp xỉ đầu u0 = 0, v0 = 2.24 với xấp xỉ đầu u0 = 0, v0 = 2.25 với xấp xỉ đầu u0 = 0, v0 = 3.1 3.2 3.3 3.4 3.5 3.6 Sự Sự Sự Sự Sự Sự hội hội hội hội hội hội tụ tụ tụ tụ tụ tụ phương pháp lặp Ví dụ 3.1 phương pháp lặp Ví dụ 3.2 Ví dụ 3.5 Ví dụ 3.6 Ví dụ 3.7 Ví dụ 3.8 vi 46 59 73 74 74 75 75 84 95 96 96 97 98 lưới 65 × 65 nút110 112 122 123 125 125 Mục lục Lời cam đoan i Lời cảm ơn ii Danh mục chữ viết tắt ký hiệu iii Mở đầu Chương Kiến thức bổ trợ 11 1.1 Một số định lý điểm bất động 1.1.1 Giới thiệu chung 1.1.2 Định lý điểm bất động Schauder 1.1.3 Định lý điểm bất động Banach 11 11 12 14 1.2 Hàm Green số toán 16 1.3 Đạo hàm số, tích phân số với sai số cấp cao 1.3.1 Đạo hàm số 1.3.2 Tích phân số 21 21 22 1.4 Lược đồ sai phân với độ xác cấp bốn cho phương trình Poisson 24 1.5 Phương pháp giải hệ phương trình lưới 1.5.1 Phương pháp truy đuổi giải hệ phương trình ba điểm 1.5.2 Phương pháp rút gọn hoàn toàn giải hệ phương trình véc ba điểm 26 26 tơ 29 Chương Sự tồn nghiệm phương pháp lặp giải toán biên cho phương trình vi phân thường phi tuyến cấp bốn 34 2.1 Bài tốn biên cho phương trình vi phân thường phi tuyến cấp bốn địa phương 2.1.1 Trường hợp điều kiện biên tổ hợp 2.1.2 Trường hợp điều kiện biên Dirichlet 2.1.3 Trường hợp điều kiện biên phi tuyến 35 35 50 68 2.2 Bài tốn biên cho phương trình vi phân thường phi tuyến cấp bốn không địa phương 76 2.2.1 Trường hợp điều kiện biên dạng gối - tựa đơn giản 76 2.2.2 Trường hợp điều kiện biên phi tuyến 88 vii Chương Sự tồn nghiệm phương pháp lặp giải toán biên phi tuyến cho phương trình đạo hàm riêng cấp bốn 100 3.1 Bài toán biên phi tuyến cho phương trình song điều hịa 3.1.1 Sự tồn nghiệm 3.1.2 Phương pháp giải ví dụ số 100 102 106 3.2 Bài toán biên phi tuyến cho phương trình song điều hịa loại Kirchhoff 115 3.2.1 Sự tồn nghiệm 117 3.2.2 Phương pháp giải ví dụ số 120 Kết luận chung 128 Danh mục cơng trình công bố Luận án 129 Tài liệu tham khảo 130 viii Từ (2.2.7) ta có Z L u (x) = G0x (x, t)ϕ(t)dt, (2.2.8) G0x (x, t) = 6L ( (t − L)(t2 + 3x2 − 2tL), ≤ x ≤ t ≤ L, t(t2 + 3x2 − 6xL + 2L2 ), ≤ t ≤ x ≤ L Công thức (2.2.8) đảm bảo hàm G(x, t) G0x (x, t) liên tục hình vng [0, L] × [0, L] Chú ý Z L Z L L3 5L4 (2.2.9) , max |G0x (x, t)|dt = max |G(x, t)|dt = 0≤x≤L 0≤x≤L 384 24 Bây giờ, ta đặt v(x) = u00 (x) tốn (2.2.5) đưa hai toán cấp hai ( v 00 (x) = ϕ(x), < x < L, (2.2.10) v(0) = v(L) = 0, ( u00 (x) = v(x), < x < L, (2.2.11) u(0) = u(L) = Lúc toán tử A xác định (2.2.4) biểu diễn dạng c(kyk2 )v(x) + f (x, u(x), y(x), v(x), z(x)) (Aϕ)(x) := M (2.2.12) với y(x) = u0 (x), z(x) = v (x) (2.2.13) Để A tốn tử co miền thích hợp, với số dương R, ta xét miền n 5L4 R L3 R L2 R LR o DR := (x, u, y, v, z) | ≤ x ≤ L, |u| ≤ , |y| ≤ , |v| ≤ , |z| ≤ 384 24 Cũng mục trước, ta kí hiệu B[O, R] hình cầu đóng tâm O bán kính R khơng gian C[0, L] với chuẩn kϕk = max0≤x≤L |ϕ(x)| Bổ đề 2.4 Giả sử tồn số R > 0, ≤ m ≤ , λM c, K1 , K2 , K3 , K4 ≥ L cho   mL2 c |M (s)| ≤ m, |f (x, u, y, v, z)| ≤ R − , (2.2.14) R2 L7 với (x, u, y, v, z) ∈ DR ≤ s ≤ Khi đó, tốn tử A ánh xạ B[O, R] 576 vào Ngoài ra, c(s2 ) − M c(s1 )| ≤ λ c|s2 − s1 |, |M M (2.2.15) |f (x, u2 , y2 , v2 , z2 ) − f (x, u1 , y1 , v1 , z1 )| ≤ K1 |u2 − u1 | + K2 |y2 − y1 | + K3 |v2 − v1 | + K4 |z2 − z1 |, 79 (2.2.16) với (x, ui , yi , vi , zi ) ∈ DR , ≤ si ≤ R2 L7 (i = 1, 2) 576 L3 L2 L mL2 λM 5L4 cR L + K2 + K3 + K4 + +