(Đề tài NCKH) chỉnh hóa phương trình helmholtz có hiệu chỉnh

BỘ GIÁO DỤC VÀ ĐÀO TẠO TRƯỜNG ĐẠI HỌC SƯ PHẠM KỸ THUẬT THÀNH PHỐ HỒ CHÍ MINH BÁO CÁO TỔNG KẾT ÐỀ TÀI KH&CN CẤP TRƯỜNG TRỌNG ÐIỂM CHỈNH HĨA PHƯƠNG TRÌNH HELMHOL TZ CĨ HIỆU CHỈNH MÃ SỐ: T2013 - 159 SKC005723 Tp Hồ Chí Minh, tháng 02/2014 BỘ GIÁO DỤC VÀ ĐÀO TẠO TRƯỜNG ĐẠI HỌC SƯ PHẠM KỸ THUẬT THÀNH PHỐ HỒ CHÍ MINH BÁO CÁO TỔNG KẾT ĐỀ TÀI KH&CN CẤP TRƯỜNG CHỈNH HÓA PHƯƠNG TRÌNH HELMHOLTZ CĨ HIỆU CHỈNH Mã số: T2013 - 159 Chủ nhiệm đề tài: ThS NGUYỄN QUANG HUY TP HCM, 2/2014 Danh sách thành viên tham gia nghiên cứu đơn vị phối hợp Thành viên : 1) ThS Nguyễn Quang Huy – ĐH Sư phạm Kỹ thuật TpHCM MỤC LỤC Trang Danh mục bảng biểu Danh mục chữ viết tắt Thông tin kết nghiên cứu Tiếng Việt Tiếng Anh Mở đầu ………………………………………………………… ……… Chương Kiến thức chuẩn bị …………………………………………………… 10 Chương Các kết chỉnh hóa………………………………………………… 18 1.1 Giới thiệu toán……………………………………………………… 18 1.2 Biến đổi toán………………………………………………………… 18 1.3 Xây dựng nghiệm chỉnh hóa tốn …………………… 19 1.4.Tính ổn định nghiệm chỉnh hóa 20 1.5 Đánh giá sai số nghiệm chỉnh hóa nghiệm xác 23 Chương Ví dụ số ………………………………………………………….29 Kết luận …………………………………………………………………………… 35 Tài liệu tham khảo………………………………………………………………… 36 Danh mục bảng biểu Bảng 2.1 Bảng đánh giá sai số sai số nghiệm chỉnh hóa u ( f )(0,.) nghiệm xác uex ( fex )(0,.) Bảng 2.2 Bảng đánh giá sai số sai số nghiệm chỉnh hóa v ( f )(0,.) nghiệm xác uex ( fex )(0,.) Danh mục chữ viết tắt THƠNG TIN KẾT QUẢ NGHIÊN CỨU Thơng tin chung: - Tên đề tài: Chỉnh hóa phương trình Helmholtz có hiệu chỉnh - Mã số: T2013 - 159 - Chủ nhiệm: ThS Nguyễn Quang Huy - Cơ quan chủ trì: Đại học Sư phạm Kỹ thuật Thành phố Hồ Chí Minh - Thời gian thực hiện: Từ tháng năm 2013 đến tháng năm 2014 Mục tiêu: - Khảo sát chỉnh hóa nghiệm phương trình Helmholtz có hiệu chỉnh Đưa đánh giá sai số ví dụ minh họa Tính sáng tạo: - Dạng sai số đánh giá nghiệm chỉnh hóa nghiệm xác Kết nghiên cứu: - Xây dựng nghiệm chỉnh hóa tốn - Đánh giá tính ổn định nghiệm chỉnh hóa - Đánh giá sai số nghiệm chỉnh hóa nghiệm xác - Đưa ví dụ số minh họa cho phương pháp Sản phẩm: - Báo cáo hội thảo khoa học khoa, tóm tắt đăng trang web khoa Hiệu quả, phương thức chuyển giao kết nghiên cứu khả áp dụng: - Làm tài liệu tham khảo cho sinh viên học viên cao học ngành Toán Trưởng đơn vị Chủ nhiệm đề tài (ký, họ tên) (ký, họ tên) University of Technical Education HCM City Faculty of Foundation Sciences HCMC, Feb 17th , 2014 INFORMATION ON RESEARCH RESULT General information - Title of project: Regularization of a modified Helmholtz equation - Code number: T2013 – 159 - Chairman: Nguyen Quang Huy - Executive organization: University of Technical Education HCMC - Duration: From January 2013 to February 2013 Objective We studied and regularized the solution of a modified Helmholtz equation Besides, we gave the error estimates and numerical example Creativeness and innovativeness The error estimates is new Research results - Build the regularized solutions - Consider the stability of the regularized solutions - Estimate the errors of the regularized solutions and the exact solution - Give the numerical example for the theory Products - Report at the science seminar of Faculty, post the abstract on the website of Faculty Effects, transfer alternatives of research results and applicability The results established in this project are references for PhD and Master students who major in mathematics MỞ ĐẦU I) Giới thiệu tổng quan Bài tốn Cauchy cho phương trình Helmholtz có ý nghĩa nhiều ngành khoa học kỹ thuật Dạng tốn thường khơng chỉnh nhiều nhà Toán học nước nước quan tâm nghiên cứu , đưa nhiều phương pháp chỉnh hóa cho tốn Trong [4], tác giả N H Tuấn, P H Quân Đ Đ Trọng chỉnh hóa tốn u k 2u 0, u (0, y ) u ( , y ) 0, u y ( x ,0) f ( x ), u ( x ,0) g ( x ), Trong [5], tác giả P H Quân P T Hiếu khảo sát toán u k 2u 0, u y ( x ,0) ( x ), u ( x ,0) ( x ), Trong [6], H H Qin T Wei sử dụng phương pháp Tikhonov cải biên phương pháp “chặt cụt” để chỉnh hóa tốn u k 2u 0, ( x , y) (0, ) (0, 1) u (0, y ) u ( , y ) 0, y (0, 1) x (0, ) u y ( x ,0) 0, u ( x ,0) x (0, ( x ), ) Trong [7], tác giả R Shi, T Wei, H H Qin cho phương trình Helmholtz có hiệu chỉnh cách chỉnh hóa tốn Cauchy thêm “đạo hàm bậc 4” vào phương trình u k 2u u x (0, y ) u (0, y ) Trong đề tài này, chỉnh hóa tốn (III) theo phương pháp dạng sai số đánh giá khác so với [7] Kết chỉnh hóa trình bày chương Hơn nữa, chương 2, tơi đưa ví dụ số minh họa cho phương pháp II) Tính cấp thiết đề tài Hiện nay, toán ứng dụng xuất nhiều lĩnh vực vật lý, y học, khoa học công nghệ Đặc biệt tốn ngược có nhiều ý nghĩa Vật lý nhiều nhà Toán học quan tâm, nghiên cứu thời TÀI LIỆU THAM KHẢO [1] Andreas Kirsch, An introduction to the mathematical theory of inverse problem, Springer, 1966 [2] Đặng Đức Trọng, Giáo trình Giải tích thực, Đề tài nghiên cứu khoa học cấp trường , Đại học Khoa học tự nhiên Tp HCM, 2005 [3] Đặng Đình Áng, Trần Lưu Cường, Huỳnh Bá Lân, Nguyễn Văn Nhân, Phạm Hồng Qn, Biến đổi tích phân, Nhà xuất Giáo dục, 2007 [4] Nguyen Huy Tuan, Pham Hoang Quan and Nguyen Duc Trong, Regularization and new error estimates for a modified helmholtz equation, An St Univ Ovidius Constanta, Vol 18(2), 267-280, 2010 [5] Pham Hoang Quan and Phan Trung Hieu, A new regularization method for the Cauchy problem of the Helmholtz equation with nonhomogeneous Cauchy data, Acta Mathematica Vietnamica, Volume 36, Number 2, pp 419-430, 2011 [6] H H Qin, T.Wei, Two regularization methods for the Cauchy problems of the Helmhotz equation, Applied Mathematical Modelling 34 , 947-967, 2010 [7] R Shi, T Wei, H H Qin, Fourth – Order Modified for the Cauchy problem of the modified Helmholtz equation, Numer Math Theor Meth Appl, 2, 326-340, 2009 36 CHỈNH HĨA PHƯƠNG TRÌNH HELMHOLTZ CĨ HIỆU CHỈNH ThS Nguyễn Quang Huy Khoa Khoa học _ ĐH SPKT Tóm tắt Trong viết này, khảo sát tốn Cauchy cho phương trình Helmholtz có hiệu chỉnh Bài tốn tốn khơng chỉnh Để chỉnh hóa tốn này, chúng tơi xây dựng nghiệm chỉnh hóa , đánh giá tính ổn định nghiệm chỉnh hóa đánh giá sai số nghiệm chỉnh hóa nghiệm xác Cuối cùng, chúng tơi đưa ví dụ số minh họa cho phương pháp Abstract In this article, we consider a Cauchy problem of the Helmholtz equation This problem is illposed problem For regularizing this problem, we construct approximate solutions , consider stability characteristic and evaluate the errors between approximate solutions and exact solution I) Giới thiệu tổng quan Bài toán Cauchy cho phương trình Helmholtz có ý nghĩa nhiều ngành khoa học kỹ thuật Dạng toán thường khơng chỉnh nhiều nhà Tốn học nước nước quan tâm nghiên cứu , đưa nhiều phương pháp chỉnh hóa cho tốn Trong [4], tác giả N H Tuấn, P H Qn Đ Đ Trọng chỉnh hóa tốn u k 2u 0, ( x , y) (0, ) (0, 1) u (0, y ) u ( , y ) 0, y (0,1) ( x, y) (0, ) (0, 1) u y ( x ,0) f ( x ), u ( x ,0) x (0, ) g ( x ), Trong [5], tác giả P H Quân P T Hiếu khảo sát toán 37 u k 2u 0, u y ( x ,0) ( x ), u ( x ,0) ( x ), Trong [6], H H Qin T Wei sử dụng phương pháp Tikhonov cải biên phương pháp “chặt cụt” để chỉnh hóa tốn k 2u 0, u ( x , y) (0, ) (0, 1) u (0, y ) u ( , y ) 0, y (0, 1) x (0, ) u y ( x ,0) 0, u ( x ,0) ( x ), x (0, ) Trong [7], tác giả R Shi, T Wei, H H Qin chỉnh hóa tốn Cauchy cho phương trình Helmholtz có hiệu chỉnh cách thêm “đạo hàm bậc 4” vào phương trình u k 2u 0, u x (0, y ) 0, u (0, y ) f ( y ), Trong đề tài này, tơi chỉnh hóa toán (I) theo phương pháp dạng sai số đánh giá khác so với [7] Kết chỉnh hóa trình bày chương Hơn nữa, chương 2, tơi đưa ví dụ số minh họa cho phương pháp II) Các kết chỉnh hóa Xây dựng nghiệm chỉnh hóa tốn Ta tìm nghiệm xác tốn (I) (1) Theo đó, ta xây dựng nghiệm chỉnh hóa tốn sau u (x,y) (2) 38 v (x, y) với( ) (0, 1) tham số chỉnh hóa Bổ đề (Tính ổn định nghiệm chỉnh hóa cho (2)) Cho f1 , f L1 ( ) L2 ( ) tương ứng với liệu f1 , f2 x 0, , (0, 1) Bổ đề (Tính ổn định nghiệm chỉnh hóa cho (3)) Cho f1 , f L1 ( ) L2 ( ) v ( f1 ), v ( f2 ) tương ứng với liệu f1 , f2 ,ta x 0, , (0, 1), k Đánh giá sai số nghiệm chỉnh hóa nghiệm xác Định lý Cho f , f ex L1 ( ) L2 ( ) cho hai nghiệm cho (2) (1) tương ứng với liệu f , Giả sử Q e 39 Nếu chọn Q 1 Định lý Cho 0Min 1, v ( f ), u ex ( fex ) hai nghiệm cho (3) (1) tương ứng với liệu Giả sử Nếu chọn T III) k Ví dụ số Ta xét ví dụ cụ thể minh họa cho tính tốn lý thuyết chương Xét tốn tìm uex ( x , y) thỏa mãn phương trình u xx ( x , y ) u yy ( x, y ) đồng thời thỏa hai điều kiện sau y u (0, y ) f ex ( y ) e , u x (0, y) 4u ( x , y ) , x [0, 1], y 40 Ta có Xét liệu nhiễu Ta có f f Chọn ˆ u (x, v (x,) Từ định lý 1.1, ta có Từ định lý 1.2, ta có ex T k Tại x ,ta có bảng đánh giá sai số hình vẽ biến đổi Fourier nghiệm xác nghiệm xấp xỉ sau Đánh giá sai số 10 41 10 10 10 20 10 50 10 100 Bảng Bảng đánh giá sai số sai số nghiệm chỉnh hóa u ( f )(0,.) nghiệm xác uex ( fex )(0,.) Đánh giá sai số 10 10 10 10 20 10 50 10 100 Bảng Bảng đánh giá sai số sai số nghiệm chỉnh hóa v ( f )(0,.) nghiệm xác uex ( fex )(0,.) 42 Hình Hình vẽ biến đổi Fourier nghiệm xác uex (0,.) Hình Các hình vẽ biến đổi Fourier nghiệm xấp xỉ u (0,.) (k k 43 1,2,3,4,5) Hình Các hình vẽ biến đổi Fourier nghiệm xấp xỉ v (0,.) (k 1,2,3,4,5) k IV) Kết luận Trong đề tài này, chỉnh hóa nghiệm tốn Cauchy cho phương trình Helmholtz có hiệu chỉnh Trong chương 1, tơi tính khơng chỉnh tốn, xây dựng hai nghiệm chỉnh hóa tốn, đánh giá tính ổn định nghiệm chỉnh hóa đánh giá sai số nghiệm chỉnh hóa nghiệm xác Trong chương 2, chúng tơi đưa ví dụ số minh họa cho phương pháp lập trình phần mềm Matlab 44 Trong thời gian tới, cải thiện đánh giá sai số thực phương pháp chỉnh hóa khác cho tốn Hơn nữa, tơi khảo sát trường hợp phương trình khơng Tài liệu tham khảo [1] Andreas Kirsch, An introduction to the mathematical theory of inverse problem, Springer, 1966 [2] Đặng Đức Trọng, Giáo trình Giải tích thực, Đề tài nghiên cứu khoa học cấp trường , Đại học Khoa học tự nhiên Tp HCM, 2005 [3] Đặng Đình Áng, Trần Lưu Cường, Huỳnh Bá Lân, Nguyễn Văn Nhân, Phạm Hồng Qn, Biến đổi tích phân, Nhà xuất Giáo dục, 2007 [4] Nguyen Huy Tuan, Pham Hoang Quan and Nguyen Duc Trong, Regularization and new error estimates for a modified helmholtz equation, An St Univ Ovidius Constanta, Vol 18(2), 267-280, 2010 [5] Pham Hoang Quan and Phan Trung Hieu, A new regularization method for the Cauchy problem of the Helmholtz equation with nonhomogeneous Cauchy data, Acta Mathematica Vietnamica, Volume 36, Number 2, pp 419-430, 2011 [6] H H Qin, T.Wei, Two regularization methods for the Cauchy problems of the Helmhotz equation, Applied Mathematical Modelling 34 , 947-967, 2010 [7] R Shi, T Wei, H H Qin, Fourth – Order Modified for the Cauchy problem of the modified Helmholtz equation, Numer Math Theor Meth Appl, 2, 326-340, 2009 45 ... H Qin cho phương trình Helmholtz có hiệu chỉnh cách chỉnh hóa tốn Cauchy thêm “đạo hàm bậc 4” vào phương trình u k 2u u x (0, y ) u (0, y ) Trong đề tài này, chỉnh hóa tốn (III) theo phương pháp... sát chỉnh hóa nghiệm phương trình Helmholtz có hiệu chỉnh Đưa đánh giá sai số ví dụ minh họa Tính sáng tạo: - Dạng sai số đánh giá nghiệm chỉnh hóa nghiệm xác Kết nghiên cứu: - Xây dựng nghiệm chỉnh. .. tượng nghiên cứu Một phương trình Helmholtz có hiệu chỉnh b) Phạm vi nghiên cứu Chỉnh hóa nghiệm có đánh giá sai số VI) Nội dung nghiên cứu: Chương 1: Xây dựng nghiệm chỉnh hóa tốn , khảo sát tính