(Luận văn) các phương pháp giải tích giải bài toán phương trình vật lý toán

BỘ GIÁO DỤC VÀ ĐÀO TẠO TRƯỜNG ĐẠI HỌC SƯ PHẠM THÀNH PHỐ HỒ CHÍ MINH KHOA VẬT LÝ  lu an n va p ie gh tn to HUỲNH TRÚC PHƯƠNG d oa nl w KHÓA LUẬN TỐT NGHIỆP va an lu CÁC PHƯƠNG PHÁP GIẢI TÍCH oi lm ul nf GIẢI BÀI TỐN PHƯƠNG TRÌNH z at nh VẬT LÝ - TOÁN z @ m co l gm Chuyên ngành: Sư phạm Vật lý an Lu TP Hồ Chí Minh, tháng 04 năm 2019 n va ac th si BỘ GIÁO DỤC VÀ ĐÀO TẠO TRƯỜNG ĐẠI HỌC SƯ PHẠM TP HỒ CHÍ MINH KHOA VẬT LÝ lu an n va CÁC PHƯƠNG PHÁP GIẢI TÍCH ie gh tn to GIẢI BÀI TỐN PHƯƠNG TRÌNH p VẬT LÝ - TỐN d oa nl w oi lm ul nf va an lu Sinh viên thực hiện: Huỳnh Trúc Phương z at nh Người hướng dẫn khoa học: ThS Nguyễn Vũ Thụ Nhân z m co l gm @ an Lu TP Hồ Chí Minh, tháng năm 2019 n va ac th si i LỜI CẢM ƠN Tôi xin chân thành cảm ơn thầy Nguyễn Vũ Thụ Nhân – người tận tình giúp đỡ hướng dẫn tơi q trình học tập, nghiên cứu hồn thiện khóa luận Tơi xin chân thành cảm ơn Trường, Phòng đào tạo, thầy cô khoa Vật lý, trường Đại học Sư phạm TP Hồ Chí Minh tạo điều kiện thuận lợi cho tơi thực khóa luận lu Qua đây, tơi xin bày tỏ lịng biết ơn gia đình, bạn bè người thân an n va giúp đỡ, động viên, hỗ trợ thời gian thực khóa luận gh tn to TP Hồ Chí Minh, ngày 29 tháng 04 năm 2019 p ie SINH VIÊN d oa nl w lu oi lm ul nf va an Huỳnh Trúc Phương z at nh z m co l gm @ an Lu n va ac th si ii MỤC LỤC LỜI CẢM ƠN i MỤC LỤC ii DANH MỤC BẢNG BIỂU v DANH MỤC HÌNH VẼ vi MỞ ĐẦU Lí chọn đề tài I lu II Mục đích nghiên cứu an III Đối tượng nghiên cứu va n IV Nhiệm vụ nghiên cứu VI Cấu trúc đề tài gh tn to V Phạm vi nghiên cứu p ie Chương CƠ SỞ LÝ THUYẾT CỦA ĐỀ TÀI NGHIÊN CỨU 1.1.1 oa 1.1.2 Hàm delta Dirac nl w 1.1 Một số hàm đặc biệt Hàm Heaviside d Hàm Bessel 1.1.4 Đa thức Legendre va an lu 1.1.3 ul nf 1.2 Các phép biến đổi tích phân Phép biến đổi Fourier 1.2.2 Phép biến đổi Fourier Sin Cos 1.2.3 Phép biến đổi Fourier phức 1.2.4 Phép biến đổi Laplace 10 oi lm 1.2.1 z at nh z Chương CÁC PHƯƠNG PHÁP GIẢI TÍCH GIẢI CÁC PHƯƠNG TRÌNH VẬT LÝ – TỐN 15 gm @ l 2.1 PHƯƠNG PHÁP TÁCH BIẾN 15 Giới thiệu phương pháp 15 2.1.2 Phương pháp tách biến việc giải phương trình truyền sóng 15 m co 2.1.1 an Lu 2.1.2.1 Truyền sóng dây hữu hạn dao động tự 15 n va 2.1.2.2 Truyền sóng dây hữu hạn dao động cưỡng 22 ac th si iii Phương pháp tách biến việc giải phương trình truyền nhiệt 25 2.1.3 2.1.3.1 Truyền nhiệt hữu hạn không chứa nguồn 25 2.1.3.2 Truyền nhiệt hữu hạn có chứa nguồn 31 2.1.4 Phương pháp tách biến việc giải phương trình Laplace 34 2.1.5 Phương pháp tách biến hệ tọa độ khác 38 2.2 PHƯƠNG PHÁP ĐA THỨC D’ALEMBERT 44 2.2.1 Giới thiệu phương pháp 44 2.2.2 sóng Phương pháp đa thức d’Alembert việc giải phương trình truyền 44 lu an 2.2.2.1 Truyền sóng dây dài vơ hạn 44 n va 2.2.2.2 Truyền sóng dây dài nửa vô hạn 46 2.3.1 Giới thiệu phương pháp 48 2.3.2 toán Phương pháp biến đổi tích phân việc giải phương trình vật lý – 48 p ie gh tn to 2.3 PHƯƠNG PHÁP BIẾN ĐỔI TÍCH PHÂN 48 lu Nghiệm hàm Green cho phương trình sóng độc lập với thời gian 56 an 2.4.3 Hàm Green 54 d 2.4.2 Giới thiệu phương pháp 54 oa 2.4.1 nl w 2.4 PHƯƠNG PHÁP HÀM GREEN 54 ul nf va 2.4.4 Nghiệm hàm Green cho phương trình sóng khơng không gian ba chiều 60 oi lm 2.4.5 Nghiệm hàm Green cho phương trình Maxwell tốn phụ thuộc vào thời gian 62 z at nh Chương ÁP DỤNG CÁC PHƯƠNG PHÁP GIẢI TÍCH TRONG VIỆC GIẢI CÁC PHƯƠNG TRÌNH VẬT LÝ - TỐN 68 z 3.1 ÁP DỤNG PHƯƠNG PHÁP TÁCH BIẾN 68 @ Giải tốn truyền sóng 68 3.1.2 Giải toán truyền nhiệt 75 3.1.3 Giải toán Laplace 81 3.1.4 Giải toán hệ tọa độ khác 88 m co l gm 3.1.1 an Lu 3.2 ÁP DỤNG PHƯƠNG PHÁP ĐA THỨC D’ALEMBERT 98 n va 3.3 ÁP DỤNG PHƯƠNG PHÁP BIẾN ĐỔI TÍCH PHÂN 101 ac th si iv 3.4 ÁP DỤNG PHƯƠNG PHÁP HÀM GREEN 115 KẾT LUẬN VÀ HƯỚNG PHÁT TRIỂN 126 TÀI LIỆU THAM KHẢO 127 lu an n va p ie gh tn to d oa nl w oi lm ul nf va an lu z at nh z m co l gm @ an Lu n va ac th si v DANH MỤC BẢNG BIỂU Bảng 1.1 Bảng biến đổi Laplace 13 Bảng 1.2 Bảng biến đổi Laplace mở rộng 14 lu an n va p ie gh tn to d oa nl w oi lm ul nf va an lu z at nh z m co l gm @ an Lu n va ac th si vi DANH MỤC HÌNH VẼ Hình 1.1 Chu tuyến l  L mặt phẳng phức 12 Hình 2.1 Đồ thị hàm số y  u ( x, t ) 100 lu an n va p ie gh tn to d oa nl w oi lm ul nf va an lu z at nh z m co l gm @ an Lu n va ac th si MỞ ĐẦU Lí chọn đề tài I Trong vật lý, việc giải phương trình đạo hàm riêng như: phương trình truyền sóng, phương trình truyền nhiệt,… mang lại ý nghĩa quan trọng Các nhà vật lý biết dao động dây, dao động sóng nước, nhờ việc giải phương trình truyền sóng, biết biến thiên nhiệt độ theo thời gian miền cho trước nhờ việc giải phương trình truyền nhiệt, [7],[5] Để giải phương trình này, nhà vật lý thường sử dụng số phương pháp tốn học: phương pháp số, phương pháp giải tích Phương lu pháp số giải nhiều tốn phức tạp, giải nghiệm gần [4] an va Cịn phương pháp giải tích giải nghiệm cách xác trở nên khó khăn n tốn phức tạp [3] Do đó, phương pháp giải tích thường sử dụng để gh tn to giảng dạy cho sinh viên tốn vật lý chương trình học sinh viên không ie phức tạp p Hiện nay, nhiều trường đại học, sinh viên chuyên ngành vật lý học nl w phương pháp giải tích để giải phương trình vật lý tốn: phương trình truyền sóng, d oa phương trình truyền nhiệt, phương trình Laplace Mỗi loại phương trình có nhiều dạng an lu khác nhau: phương trình truyền sóng dây dài hữu hạn vơ hạn, truyền sóng dây dao động cưỡng bức; phương trình truyền nhiệt dài hữu hạn chứa nguồn va ul nf không chứa nguồn, phương trình Laplace, Các phương pháp giải tích thường oi lm sử dụng để giải phương trình phương pháp tách biến phương pháp đa thức D’Alembert Hai phương pháp dùng phổ biến khơng địi hỏi sinh viên biết z at nh nhiều kiến thức tốn phức tạp Ngồi cịn có phương pháp tìm nghiệm dễ dàng nặng kiến thức toán phương pháp biến đổi tích phân, phương z gm @ pháp hàm Green Do có nhiều dạng phương trình, nhiều phương pháp giải tích để giải chúng nên việc hệ thống lại phương pháp giải tích giải phương trình vật lý tốn l m co cần thiết Nhờ đó, sinh viên xâu chuỗi lại kiến thức học, biết thêm phương pháp mới, giúp cho việc học trở nên dễ dàng an Lu Vì vậy, nhằm đáp ứng nhu cầu trên, hệ thống lại phương pháp giải tích n va để giải tốn phương trình vật lý - tốn đề tài ac th si Mục đích nghiên cứu II Đề tài hướng đến hai mục đích sau:  Đưa hệ thống phương pháp giải tích để giải phương trình đạo hàm riêng ứng dụng rộng rãi vật lý: phương trình truyền sóng, phương trình truyền nhiệt, phương trình Laplace  Đưa hệ thống giải tập phương trình đạo hàm riêng nói phương pháp giải tích: phương pháp tách biến, phương pháp đa thức D’Alembert, phương pháp biến đổi tích phân phương pháp hàm Green lu an Đối tượng nghiên cứu III va n  Các toán đạo hàm riêng ứng dụng vật lý tn to  Các phương pháp giải tích áp dụng giải tốn đạo hàm riêng vật lý Nhiệm vụ nghiên cứu ie gh IV p  Tìm hiểu tốn đạo hàm riêng thường gặp vật lý thông qua giáo w trình, sách, tài liệu liên quan oa nl  Phân tích ưu điểm, nhược điểm phương pháp giải tích áp dụng giải d toán vật lý – toán lu Phạm vi nghiên cứu va an V nf Đề tài tập trung nghiên cứu phương pháp giải tích để giải phương trình vật oi lm ul lý - toán thường gặp: phương trình truyền sóng dây, phương trình truyền nhiệt thanh, phương trình Laplace,… Cấu trúc đề tài z at nh VI Mở đầu: Phần tơi trình bày tổng quan đề tài nghiên cứu, bao gồm: lí chọn đề z m co l Chương Cơ sở lý thuyết đề tài nghiên cứu gm cứu cấu trúc đề tài @ tài, mục đích nghiên cứu, đối tượng nghiên cứu, nhiệm vụ nghiên cứu, phạm vi nghiên Trong chương này, tơi trình bày số hàm đặc biệt đề cập tới đề tài chương an Lu phép biến đổi tích phân để làm sở cho phương pháp biến đổi tích phân n va ac th si 113 {𝑢 (𝑥, 0) = {𝜕𝑢 , ∀𝑥 ∈ (0,2) (𝑥, 0) = sin 𝜋𝑥 𝜕𝑡 (3.3.34) 𝑢 (0, 𝑡) = 𝑢 (2, 𝑡) = 0, ∀𝑡 > (3.3.35) Và điều kiện biên: Giải lu Biến đổi Laplace hai vế phương trình (3.3.33) theo biến t, ta được: an n va 𝜕2𝑢 𝜕2 𝑢 ℒ ( 2) = ( ) 𝜕𝑥 𝜕𝑡 tn to 𝑑2𝒰 𝜕𝑢 (𝑥, 𝑝) = [𝑝 𝒰(𝑥, 𝑝) − 𝑝 𝑢(𝑥, 0) − (𝑥, 0)] 𝑑𝑥 𝜕𝑡 p ie gh ⇒ w Từ điều kiện đầu (3.3.34): d oa nl 𝑑2𝒰 ⇒ = 6𝑝 𝒰(𝑥, 𝑝) − sin 𝜋𝑥 𝑑𝑥 (3.3.36) oi lm Xét phương trình nhất: ul nf va an lu 𝑑2𝒰 ⇔ − 6𝑝 𝒰 = −6 sin 𝜋𝑥 𝑑𝑥 z at nh 𝑑2𝒰 − 6𝑝 𝒰 = 𝑑𝑥 z Nghiệm phương trình nhất: 𝒰∗ = 𝐴𝑒 𝑝√6𝑥 + 𝐵𝑒 −𝑝√6𝑥 𝒰𝑅 = 𝐴 sin 𝜋𝑥 + 𝐵 cos 𝜋𝑥 ⇒ 𝒰𝑅′ = 𝜋𝐴 cos 𝜋𝑥 − 𝜋𝐵 sin 𝜋𝑥 m co l gm @ Nghiệm riêng phương trình (3.3.36) có dạng: an Lu ⇒ 𝒰𝑅′′ = −𝜋 𝐴 sin 𝜋𝑥 − 𝜋 𝐵 cos 𝜋𝑥 n va ac th si 114 Thay 𝒰𝑅 , 𝒰𝑅 ′′ vào (3.3.36), ta có: −𝜋 𝐴 sin 𝜋𝑥 − 𝜋 𝐵 cos 𝜋𝑥 − 6𝑝 (𝐴 sin 𝜋𝑥 + 𝐵 cos 𝜋𝑥 ) = −6 sin 𝜋𝑥 ⇔ (𝜋 + 6𝑝 )𝐴 sin 𝜋𝑥 + (𝜋 + 6𝑝 )𝐵 cos 𝜋𝑥 = sin 𝜋𝑥 lu 𝐴= ⇒{ 𝜋 + 6𝑝 𝐵=0 sin 𝜋𝑥 ⇒ 𝒰𝑅 = 𝜋 + 𝑝2 an Nghiệm phương trình (3.3.36) là: n va sin 𝜋𝑥 𝜋2 + 𝑝2 (3.3.37) ie gh tn to 𝒰 = 𝒰∗ + 𝒰𝑅 = 𝐴𝑒 𝑝√6𝑥 + 𝐵𝑒 −𝑝√6𝑥 + p Từ điều kiện biên (3.3.35): w nl +∞ oa 𝒰(0, 𝑝) = ∫ 𝑒 −𝑝𝑥 𝑢(0, 𝑡)𝑑𝑡 = 0 +∞ d ⇒ 𝑒 𝑢(2, 𝑡)𝑑𝑡 = 0 nf va an lu { 𝒰(2, 𝑝) = ∫ (3.3.38) −𝑝𝑥 𝒰(0, 𝑝) = 𝐴 + 𝐵 (3.3.39) 𝒰(2,0) = 𝐴𝑒 𝑝2√6 + 𝐵𝑒 −𝑝2√6 z Từ (3.3.38) (3.3.39) ⇒ 𝐴 = 𝐵 = z at nh { oi lm ul Từ (3.3.37) (3.3.38), suy ra: @ sin 𝜋𝑥 𝜋2 + 𝑝2 m co l gm ⇒𝒰= Dùng phép biến đổi Laplace ngược, ta thu nghiệm phương trình (3.3.33): an Lu n va ac th si 115 𝑢(𝑥, 𝑡) = ℒ −1 (𝒰) = √6 √6 sin 𝜋𝑥 sin 𝑡 𝜋 𝜋 3.4 ÁP DỤNG PHƯƠNG PHÁP HÀM GREEN Bài 3.20 Bằng phương pháp biến đổi Laplace, tìm nghiệm tổng quát hàm Green cho phương trình: lu 𝜕2 ( + 𝑘 ) 𝑢(𝑥, 𝑘 ) = −𝑓 (𝑥 ), ∀𝑥 ∈ [0, +∞) 𝜕𝑥 (3.4.1) an n va 𝑘 số to gh tn Giải p ie Hàm Green nghiệm phương trình: (3.4.2) d oa nl w 𝜕2 ( + 𝑘 ) 𝑔(𝑥|𝑥0 , 𝑘) = −𝛿 (𝑥 − 𝑥0 ) 𝜕𝑥 an lu Biến đổi Laplace hai vế phương trình (3.4.2): oi lm ul nf va 𝜕2 ℒ [( + 𝑘 ) 𝑔 (𝑥|𝑥0 , 𝑘)] = −ℒ(𝛿) 𝜕𝑥 +∞ m co l Sau tính tốn, ta kết quả: gm @ +∞ +∞ 𝜕 𝑔 −𝑝𝑥 −𝑝𝑥 ∫ ∫ 𝑒 𝑑𝑥 + 𝑘 𝑔𝑒 𝑑𝑥 = − 𝛿 (𝑥 − 𝑥0 )𝑒 −𝑝𝑥 𝑑𝑥 (3.4.3) 𝜕𝑥 0 z ⇒∫ z at nh 𝜕2 ⇒ ℒ ( 𝑔) + ℒ (𝑘 𝑔) = −ℒ(𝛿) 𝜕𝑥 +∞ 𝑔 (𝑥|𝑥0 , 𝑘)𝑒 −𝑝𝑥 𝑑𝑥 n va với 𝑔̅ (𝑝|𝑥0 , 𝑘 ) = ∫0 an Lu 𝑝 𝑒 −𝑝𝑥0 ′ 𝑔̅ (𝑝|𝑥0 , 𝑘) = 𝑔(0|𝑥0 , 𝑘) + 𝑔 (0|𝑥0 , 𝑘 ) + 𝑝 + 𝑘2 𝑝 + 𝑘2 𝑝 + 𝑘2 ac th si 116 Biến đổi Laplace ngược hàm 𝑔̅ (𝑝|𝑥0 , 𝑘) để thu hàm 𝑔(𝑥|𝑥0 , 𝑘): 𝑝 𝑒 −𝑝𝑥0 ′ ⇒ ℒ [𝑔̅ (𝑝|𝑥0 , 𝑘)] = ℒ [ 𝑔(0|𝑥0 , 𝑘) + 𝑔 (0|𝑥0 , 𝑘 ) + ] 𝑝 + 𝑘2 𝑝 + 𝑘2 𝑝 + 𝑘2 𝑝 𝑒 −𝑝𝑥0 ′ ⇒ ℒ [𝑔̅ (𝑝|𝑥0 , 𝑘)] = ℒ [ 𝑔(0|𝑥0 , 𝑘)] + ℒ [ 𝑔 (0|𝑥0 , 𝑘 )] + ℒ [ ] 𝑝 + 𝑘2 𝑝 + 𝑘2 𝑝 + 𝑘2 1 ⇒ 𝑔(𝑥|𝑥0 , 𝑘 ) = 𝑔(0|𝑥0 𝑘) cos 𝑘𝑥 + 𝑔′ (0, 𝑥0 , 𝑘 ) sin 𝑘𝑥 + sin[𝑘 (𝑥 − 𝑥0 )] 𝐻(𝑥 𝑘 𝑘 − 𝑥0 ) lu an Nghiệm 𝑢(𝑥0 , 𝑘) phương trình (3.4.1) có dạng: va n +∞ tn to 𝑢(𝑥0 , 𝑘 ) = ∫ 𝑔 (𝑥|𝑥0 , 𝑘)𝑓(𝑥 )𝑑𝑥 gh +∞ 𝑔(0|𝑥0 𝑘 ) cos 𝑘𝑥0 𝑓(𝑥 )𝑑𝑥 + ∫ 0 p ie +∞ +∫ ′ 𝑔 (0, 𝑥0 , 𝑘 ) sin 𝑘𝑥0 𝑓 (𝑥 )𝑑𝑥 𝑘 sin[𝑘 (𝑥 − 𝑥0 )] 𝐻 (𝑥 − 𝑥0 )𝑓(𝑥 )𝑑𝑥 𝑘 d oa nl w +∞ 𝐴(𝑘 ) = ∫ ul nf va an lu Đặt +∞ =∫ 𝑔(0|𝑥, 𝑘 )𝑓(𝑥)𝑑𝑥 oi lm +∞ ′ 𝐵(𝑘 ) = ∫ 𝑔 (0|𝑥, 𝑘 )𝑓 (𝑥 )𝑑𝑥 𝑘 { z at nh z ⇒ 𝑢 (𝑥0 , 𝑘 ) = 𝐴(𝑘 ) cos 𝑘𝑥𝑜 + 𝐵(𝑘 ) sin 𝑘𝑥0 +∞ + ∫ sin[𝑘 (𝑥 − 𝑥0 )] 𝐻 (𝑥 − 𝑥0 )𝑓 (𝑥 )𝑑𝑥 𝑘 m co l Tìm nghiệm phương trình: gm @ Bài 3.21 an Lu n va ac th si 117 𝜕2 ( + 𝑘 ) 𝑢(𝑥, 𝑘 ) = 0, ∀ 𝑥 ∈ [0, 𝐿] 𝜕𝑥 (3.4.4) 𝑢(0, 𝑘) = 𝑢(𝐿, 𝑘 ) = (3.4.5) với điều kiện: { Giải lu Hàm Green nghiệm phương trình: an va n ( 𝜕2 + 𝑘 ) 𝑔(𝑥|𝑥0 , 𝑘) = −𝛿 (𝑥 − 𝑥0 ) 𝜕𝑥 (3.4.6) tn to ie gh Nghiệm hàm Green không gian chiều có dạng: p 𝑔(𝑥|𝑥0 , 𝑘) = − sin 𝑘|𝑥 − 𝑥0 | 𝑘 nl w d oa Với điều kiện (3.4.5), nghiệm riêng phương trình (3.4.4) có dạng: an lu 𝑢 (𝑟, 𝑘 ) = 𝐴 sin 𝑘𝑥 + 𝐵𝑠𝑖𝑛 𝑘(𝐿 − 𝑥) ul nf va Nghiệm tổng quát phương trình (2.4.4) là: oi lm 𝑢(𝑥, 𝑘 ) = − sin 𝑘|𝑥 − 𝑥0 | + 𝐴 sin 𝑘𝑥 + 𝐵 sin 𝑘(𝐿 − 𝑥) 𝑘 z at nh Từ điều kiện (3.4.5): 𝑢(0, 𝑘) = 𝑢(𝐿, 𝑘 ) = z gm @ { m co l − sin 𝑘𝑥0 + 𝐵 sin 𝑘𝐿 = 𝑘 ⇒{ − sin 𝑘(𝐿 − 𝑥0 ) + 𝐴 sin 𝑘𝐿 = 𝑘 an Lu n va ac th si 118 sin 𝑘𝑥0 𝑘 sin 𝑘𝐿 ⇒{ sin 𝑘(𝐿 − 𝑥0 ) 𝐴= 𝑘 sin 𝑘𝐿 𝐵= Vậy nghiệm phương trình (3.4.4) là: 1 sin 𝑘 (𝐿 − 𝑥0 ) sin 𝑘𝑥0 𝑢(𝑥, 𝑘 ) = − sin 𝑘|𝑥 − 𝑥0 | + sin 𝑘𝑥 + sin 𝑘(𝐿 − 𝑥) 𝑘 𝑘 sin 𝑘𝐿 𝑘 sin 𝑘𝐿 lu sin 𝑘 (𝐿 − 𝑥0 ) sin 𝑘𝑥0 ⇒ 𝑢(𝑥, 𝑘 ) = − [sin 𝑘|𝑥 − 𝑥0 | − sin 𝑘𝑥 − sin 𝑘 (𝐿 − 𝑥 )] 𝑘 sin 𝑘𝐿 sin 𝑘𝐿 an n va Bài 3.22 p ie gh tn to Chứng minh rằng, phương pháp biến đổi Laplace theo biến x phương trình (3.4.7) 𝑔(0|𝑥0 , 𝑘) = { 𝜕𝑔(𝑥|𝑥0 , 𝑘 ) [ ]| = 𝑔(1|𝑥0 , 𝑘) 𝜕𝑥 𝑥=0 (3.4.8) w 𝜕2 ( + 𝑘 ) 𝑔(𝑥|𝑥0 , 𝑘) = −𝛿 (𝑥 − 𝑥0 ) 𝜕𝑥 d oa nl với điều kiện biên: oi lm sin(𝑘𝑥 ) sin[𝑘 (1 − 𝑥0 )] sin[𝑘(𝑥 − 𝑥0 )] − 𝐻 (𝑥 − 𝑥0 ) 𝑘 (sin 𝑘 − 𝑘 ) 𝑘 z at nh 𝑔(𝑥|𝑥0 , 𝑘 ) = ul nf va an lu nghiệm hàm Green có dạng: z Giải gm @ Biến đổi Laplace hai vế phương trình (3.4.7), ta được: l m co 𝜕2𝑔 ℒ ( ) + 𝑘 ℒ(𝑔) = −ℒ(𝛿) 𝜕 𝑥 an Lu Sử dụng kết 2.20, ta được: n va ac th si 119 𝑔̅ (𝑝|𝑥0 , 𝑘 ) = 1 𝑒 −𝑝𝑥0 ′( ( ) ) 𝑔 0|𝑥 , 𝑘 + 𝑔 0|𝑥 , 𝑘 + 0 𝑝2 + 𝑘 𝑝2 + 𝑘 𝑝2 + 𝑘 Mà: 𝑔(0|𝑥0 , 𝑘 ) = 𝑒 −𝑝𝑥0 ′ ⇒ 𝑔̅ (𝑝|𝑥0 , 𝑘 ) = 𝑔 (0|𝑥0 , 𝑘) + 𝑝 + 𝑘2 𝑝 + 𝑘2 (3.4.9) Biến đổi Laplace ngược hàm 𝑔̅ (𝑝|𝑥0 , 𝑘), ta thu hàm 𝑔(𝑥|𝑥0 , 𝑘 ): lu 𝑔(𝑥|𝑥0 , 𝑘) = an 1 sin 𝑘𝑥 𝑔′ (0|𝑥0 , 𝑘 ) + sin[𝑘(𝑥 − 𝑥0 )] 𝐻(𝑥 − 𝑥0 ) 𝑘 𝑘 n va Mặt khác: to p ie gh tn 𝑔′ (0|𝑥0 , 𝑘) = 𝑔(1|𝑥0 , 𝑘) ⇒ 𝑔(𝑥|𝑥0 , 𝑘 ) = w 1 sin 𝑘𝑥 𝑔(1|𝑥0 , 𝑘 ) + sin[𝑘(𝑥 − 𝑥0 )] 𝐻 (𝑥 − 𝑥0 ) (3.4.10) 𝑘 𝑘 d oa nl Thay 𝑥 = vào phương trình (3.4.10): 1 sin 𝑘 𝑔(1|𝑥0 , 𝑘 ) + sin[𝑘(𝑥 − 𝑥0 )] 𝐻(1 − 𝑥0 ) 𝑘 𝑘 va an lu 𝑔(1|𝑥0 , 𝑘) = Thay (3.4.11) vào (3.4.10) (3.4.11) z at nh ⇒ 𝑔(𝑥|𝑥0 , 𝑘 ) = sin[𝑘(1 − 𝑥0 )] 𝐻(1 − 𝑥0 ) 𝑘 − sin 𝑘 oi lm ul nf ⇒ 𝑔(1|𝑥0 , 𝑘) = z sin 𝑘𝑥 sin[𝑘(1 − 𝑥0 )] + sin[𝑘(𝑥 − 𝑥0 )] 𝐻(𝑥 − 𝑥0 ) 𝑘(𝑘 − sin 𝑘 ) 𝑘 @ Trường điện U thỏa mãn phương trình: m co l gm Bài 3.23 𝜕2 𝑈 (𝑟, 𝑡) = −4𝜋𝜌 (𝑟)𝑒 𝑖𝜔𝑡 2 𝑐 𝜕𝑡 (3.4.12) an Lu ∇2 𝑈 (𝑟, 𝑡) − n va ac th si 120 với 𝜌 mật độ điện tích, 𝜔 tần số góc 𝑐 tốc độ lan truyền sóng điện từ chân khơng Sử dụng phương pháp hàm Green, tính biên độ trường điện cung cấp antenna phát bước sóng 𝜆 = 10𝑚, khoảng cách 1000m từ antenna, với 𝜌(𝑟) = 𝑟2 Giải Đặt lu 𝑈 (𝑟, 𝑡) = 𝑢(𝑟, 𝜔)𝑒 𝑖𝜔𝑡 an n va Phương trình (3.4.12) thành: to 𝜔2 + 𝑢(𝑟, 𝜔)𝑒 𝑖𝜔𝑡 = −4𝜋𝜌(𝑟)𝑒 𝑖𝜔𝑡 𝑐 ie gh tn ∇ 𝑢 (𝑟, 𝜔)𝑒 𝑖𝜔𝑡 p ⇔ ∇2 𝑢(𝑟, 𝜔) + w 𝜔2 𝑢(𝑟, 𝜔) = −4𝜋𝜌 𝑐2 (3.4.13) d oa nl ⇔ ∇2 𝑢 (𝑟, 𝜔) + 𝑘 𝑢 = −4𝜋𝜌 va an lu Nghiệm chứa hàm Green cho phương trình (3.4.13) điểm 𝑟0 : 𝑢(𝑟0 , 𝑘 ) = 4𝜋 ∫ 𝜌 (𝑟)𝑔 (𝑟|𝑟0 , 𝑘)𝑑 𝑟 oi lm ul nf (3.4.14) Hàm Green ứng với sóng khơng gian chiều: z at nh 𝑒 𝑖𝑘|𝑟−𝑟0| 𝑔(𝑟|𝑟0 , 𝑘) = 4𝜋|𝑟 − 𝑟0 | z gm 𝑖𝑘𝑟 −𝑖𝑘𝑛̂ 𝑟 𝑒 𝑒 4𝜋𝑟0 m co l 𝑔(𝑟|𝑟0 , 𝑘) = @ Sử dụng phương pháp tiệm cận, ta thu được: an Lu ⇒ 𝑢(𝑟0 , 𝑘 ) = 𝑒 𝑖𝑘𝑟0 ∫ 𝜌(𝑟) 𝑒 −𝑖𝑘𝑛̂0𝑟 𝑑 𝑟 n va ac th si 121 Đổi sang hệ tọa độ cực thay 𝜌(𝑟) = 𝑟2 𝑒 𝑖𝑘𝑟0 ∫ ∫ ∫ 𝑒 −𝑖𝑘𝑛̂0𝑟 𝑟 𝑑(cos 𝜃)𝑑𝑟𝑑𝜑 𝑢(𝑟𝑜 , 𝑘 ) = 𝑟0 𝑟 𝜋 2𝜋 𝑒 𝑖𝑘𝑟0 +∞ −𝑖𝑘𝑟 cos 𝜃 ( ) ∫ 𝑑𝑟 ∫ 𝑒 = 𝑑 cos 𝜃 ∫ 𝑑𝜑 𝑟0 0 = 2𝜋 𝜋 𝑒 𝑖𝑘𝑟0 +∞ ∫ 𝑑𝑟 ∫ 𝑒 −𝑖𝑘𝑟 cos 𝜃 𝑑 (cos 𝜃) 𝑟0 0 𝜋 lu Đặt 𝐼 = ∫ 𝑒 −𝑖𝑘𝑟 cos 𝜃 𝑑 (cos 𝜃) an n va gh tn to −𝑖𝑘𝑟 cos 𝜃 𝜋 (𝑒 𝑖𝑘𝑟 − 𝑒 −𝑖𝑘𝑟 ) = sin 𝑘𝑟 ⇒𝐼=− 𝑒 | = 𝑖𝑘𝑟 𝑖𝑘𝑟 𝑘𝑟 p ie ⇒ 𝑢(𝑟𝑜 , 𝑘 ) = 4𝜋 𝑖𝑘𝑟 +∞ sin(𝑘𝑟) 2𝜋 𝑖𝑘𝑟 𝑒 0∫ 𝑑𝑟 = 𝑒 𝑟0 𝑘𝑟 𝑟𝑜 𝑘 d Ta có: oa nl w  Tính biên độ trường điện thế: an lu 2𝜋 = 0,2𝜋 𝜆 nf va 𝑘= cấp antenna: oi lm ul Thay 𝑘 = 0,2 𝜋 𝑟0 = 1000𝑚 vào 𝑢, ta tìm biên độ trường điện cung z at nh 2𝜋 𝑖𝑘𝑟 2𝜋 𝑢= 𝑒 = 𝑒 0,2𝜋.1000.𝑖 = 0,01𝜋 𝑟𝑜 𝑘 1000.0,2𝜋 z @ gm Bài 3.24 m co l Tính tốn hàm Green khơng gian ba chiều cho tốn tử Klein – Gordon mơ tả trường sóng phụ thuộc vào thời gian với toán tử Klein – Gordon: an Lu n va 𝜕2 ∇ − 2 − 𝜎2 𝑐 𝜕𝑡 ac th si 122 Giải Theo đề bài, ta có phương trình: 𝜕2 (∇ − 2 − 𝜎 ) 𝐺 (𝑟|𝑟0 , 𝑡|𝑡0 ) = −𝛿 (𝑟 − 𝑟0 )𝛿 (𝑡 − 𝑡0 ) 𝑐 𝜕𝑡 Đặt { (3.4.15) 𝑅 = |𝑟 − 𝑟0 | 𝜏 = 𝑡 − 𝑡0 Đặt: lu +∞ ∫ 𝑔(𝑅, 𝜔)𝑒 𝑖𝜔𝜏 𝑑𝜔 𝐺 (𝑅, 𝜏) = 2𝜋 −∞ +∞ 𝑖𝜔𝜏 ∫ 𝑒 𝑑𝜔 𝛿 (𝜏) = 2𝜋 −∞ { an n va (3.4.16) tn to ie gh Thay (3.4.16) vào (3.4.15), ta được: p +∞ 𝜕 +∞ 𝑖𝜔𝜏 𝑖𝜔𝜏 3( ) ∫ (∇ − 2 − 𝜎 ) 𝑔 𝑒 𝑑𝜔 = −𝛿 𝑅 ∫ 𝑒 𝑑𝜔 2𝜋 −∞ 𝑐 𝜕𝑡 2𝜋 −∞ ∇ 𝑔𝑒 𝑖𝜔𝜏 +∞ 𝑑𝜔 − ∫ d oa −∞ nl ⇒∫ w +∞ +∞ 2) 𝑖𝜔𝜏 ( 𝑖𝜔 𝑔𝑒 𝑑𝜔 − ∫ 𝜎 𝑔𝑒 𝑖𝜔𝜏 𝑑𝜔 𝑐 −∞ lu −∞ +∞ (3.4.17) oi lm ∞ +∞ 𝜔2 𝑖𝜔𝜏 3( ) ∫ (∇ + − 𝜎 ) 𝑔𝑒 𝑑𝜔 = −𝛿 𝑅 𝑒 𝑖𝜔𝜏 𝑑𝜔 𝑐 −∞ ul ⇒∫ −∞ nf +∞ 𝑒 𝑖𝜔𝜏 𝑑𝜔 va an = −𝛿 (𝑅) ∫ @ m co l an Lu { gm +∞ ∫ 𝑔̅ (𝑢, 𝜔)𝑒 𝑖𝑢𝑅 𝑑 𝑢 (2𝜋) −∞ +∞ 3( ) ∫ 𝛿 𝑅 = 𝑒 𝑖𝑢𝑅 𝑑 𝑢 (2𝜋) −∞ 𝑔(𝑅, 𝜔) = (3.4.18) ⇒ (∇2 + (3.4.18) z Đặt z at nh 𝜔2 ⇒ (∇ + − 𝜎 ) 𝑔(𝑅, 𝑘 ) = −𝛿 (𝑅) 𝑐 n va +∞ +∞ 𝜔2 1 𝑖𝑢𝑅 ( ) ∫ ∫ − 𝜎 ) 𝑔̅ 𝑢, 𝜔 𝑒 𝑑 𝑢 = − 𝑒 𝑖𝑢𝑅 𝑑 𝑢 3 (2𝜋) −∞ (2𝜋) −∞ 𝑐 ac th si 123 +∞ +∞ 𝜔2 𝑖𝑢𝑅 (∇ + − 𝜎 ) 𝑔̅ (𝑢, 𝜔)𝑒 𝑑 𝑢 = − ∫ 𝑒 𝑖𝑢𝑅 𝑑 𝑢 𝑐 −∞ ⇒∫ −∞ 𝜔2 ⇒ (−𝑢 + − 𝜎 ) 𝑔̅ = −1 𝑐 ⇒ 𝑔̅ (𝑢, 𝜔) = 𝑢2 − 𝑘 + 𝜎 +∞ 𝑒 𝑖𝑢𝑅 ∫ ⇒ 𝑔(𝑅, 𝜔) = 𝑑3𝑢 (2𝜋)3 −∞ 𝑢2 − 𝑘 + 𝜎 lu an (2.4.17) ⇒ 𝐺 (𝑅, 𝜏) = +∞ ∫ (2𝜋)4 −∞ n va 𝑒 𝑖𝑢𝑅 𝑒 𝑖𝜔𝜏 𝑑 𝑢𝑑𝜔 𝜔 𝑢2 − + 𝜎 𝑐 (3.4.19) tn to Đặt +∞ gh 𝑒 𝑖𝑢𝑅 𝑑 𝑢 𝐼=∫ ie −∞ p w Đổi sang hệ tọa độ cầu: nl +∞ 𝑒 𝑑 𝑢=∫ −∞ d oa 𝐼=∫ +∞ 𝑖𝑢𝑅 𝑒 𝑖𝑢𝑅 𝑢2 sin 𝜃 𝑑𝑢𝑑𝜃𝑑𝜑 −∞ 𝜋 =∫ 𝑢 𝑑𝑢 ∫ 𝑒 −∞ 2𝜋 𝑖𝑢𝑅𝑐𝑜𝑠 𝜃 sin 𝜃𝑑𝜃 ∫ 𝑑𝜑 0 nf va an lu +∞ oi lm ul Đặt 𝑡 = cos 𝜃 ⇒ 𝑑𝑡 = − sin 𝜃 𝑑𝜃 +∞ 4𝜋 +∞ ∫ 𝑢 sin 𝑢𝑅 𝑑𝑢 ⇒ 𝐼 = 2𝜋 ∫ ∫ 𝑢 𝑒 𝑑𝑡𝑑𝑢 = 𝑖𝑅 −∞ −1 −∞ +∞ 𝑢 sin 𝑢𝑅 𝑒 𝑖𝜔𝜏 ⇒ 𝐺 (𝑅, 𝜏) = ∫ 𝑑𝑢𝑑𝜔 4𝜋 𝑖𝑅 −∞ 𝜔 2 𝑢 − +𝜎 𝑐 𝑖𝑢𝑅𝑡 z at nh z @ 𝜔2 𝑐2 − 𝜎2 an Lu Xét cực điểm dương 𝑧 = √ m co 𝑧𝑒 𝑖𝑅𝑧 ∮ 𝑑𝑧 𝐶 𝑧2 − 𝜔 + 𝜎2 𝑐2 l gm Để tìm hàm 𝐺(𝑅, 𝜏), ta tính tính tích phân chu tuyến C: n va ac th si 124 Theo định lý thặng dư: 𝑛 ∮𝑓(𝑧)𝑑𝑧 = 2𝜋𝑖 ∑ Res 𝑓(𝑧) 𝐶 𝑗=1 (𝑧 − √ Res ω2 z=√ −𝜎 𝑐 𝑓(𝑧) = lim 𝜔2 𝑧→√ −𝜎2 (𝑧 𝑐 z=a 𝜔2 − 𝜎 ) 𝑧𝑒 𝑖𝑅𝑧 𝑐2 𝜔2 𝜔2 − √ − 𝜎 ) (𝑧 + √ − 𝜎 ) 𝑐 𝑐 = 𝜔2 𝑖𝑅√ −𝜎 𝑒 𝑐 𝜔2 𝑧𝑒 𝑖𝑅𝑅𝑧 𝑖𝑅√ −𝜎2 𝑐 𝑑𝑧 = 𝜋𝑖𝑒 𝐶 𝑧2 − 𝜔 + 𝜎2 𝑐2 lu ⇒∮ an n va tn to ⇒ 𝐺 (𝑅, 𝜏) = +∞ 𝜔2 𝑖𝑅√ −𝜎2 𝑖𝜔𝜏 𝑐 ∫ 𝜋𝑖 𝑒 𝑒 𝑑𝜔 4𝜋 𝑖𝑅 −∞ ie gh Đặt 𝑖𝜔 = 𝑝 p Ta có: d oa nl w 𝜔2 𝜔2 √𝜎 𝑐 + 𝑝 √ − 𝜎 = 𝑖√ − 𝜎 = 𝑐2 𝑐2 𝑐 lu +∞ √𝜎 𝑐 +𝑝2 −𝑅 𝑐 ∫ 𝑒 𝑒 𝑝𝜏 𝑑𝑝 4𝜋 𝑅 −∞ +∞ √𝜎 𝑐 +𝑝2 𝜕 −𝑐 −𝑅 𝑐 ⇒ 𝐺 (𝑅, 𝜏) = ∫ [ 𝑒 ] 𝑒 𝑝𝜏 𝑑𝑝 4𝜋 𝑅 −∞ 𝜕𝑅 √𝜎 𝑐 + 𝑝 oi lm ul nf va an ⇒ 𝐺 (𝑅, 𝜏) = z at nh √𝜎 𝑐 +𝑝2 𝑐 𝜕 +∞ −𝑅 𝑐 ∫ − ⇒ 𝐺 (𝑅, 𝜏) = − 𝑒 𝑒 𝑝𝜏 𝑑𝑝 4𝜋 𝑅 𝜕𝑅 −∞ √𝜎 𝑐 + 𝑝 Để tìm hàm 𝐺 (𝑅, 𝜏), ta tìm biến đổi Laplace ngược hàm: z 𝑒 𝑅√𝜎 𝑐 +𝑝2 𝑐 (3.4.20) m co l √𝜎 𝑐 + 𝑝 − gm Mặt khác, ta có: @ 𝐹 (𝑝) = n va Có biến đổi ngược: an Lu 2 𝑓(𝑥 ) = {𝐽0 𝑎√𝑥 − 𝑏 , 𝑥 > 𝑏 ( 𝐽0 hàm Bessel bậc 0, a b số dương) 0, 𝑥 < 𝑏 ac th si 125 𝐹 (𝑝) = 𝑒 −𝑏√𝑝 +𝑎 √𝑝 + 𝑎2 𝑅 Áp dụng cho hàm 𝐹(𝑝) (3.1.20) với 𝑎 = 𝜎𝑐 𝑏 = , hàm Green cần tìm là: 𝑐 𝐺 (𝑅, 𝜏) = − 𝑐 𝜕 𝑅 𝑅 − ( ) (với 𝜏 > ) √ 𝐽 [𝜎𝑐 𝜏 4𝜋 𝑅 𝜕𝑅 𝑐 𝑐 lu an n va p ie gh tn to d oa nl w oi lm ul nf va an lu z at nh z m co l gm @ an Lu n va ac th si 126 KẾT LUẬN VÀ HƯỚNG PHÁT TRIỂN Về bản, hoàn thành mục tiêu ban đầu đề sau:  Hệ thống lại phương pháp giải tích: phương pháp tách biến, phương pháp D’Alembert, phương pháp biến đổi tích phân, phương pháp hàm Green để giải phương trình vật lý tốn: phương trình truyền sóng, truyền nhiệt, Laplace,  Hệ thống giải toán phương trình đạo hàm riêng vật lý theo phương pháp giải tích Tuy nhiên, thời gian có hạn nên đề tài chưa đề cập đến việc áp dụng phương pháp lu hàm Green cho số toán: giải phương trình Schrodinger phương trình an n va Helmholtz không nhất; tán xạ Rutherford, Rayleigh; nhiễu xạ Kirchhoff, Xuất phát từ hạn chế đề tài, đề tài phát triển lên thêm cách mở gh tn to Fraunhofer, Fresnel,…Các tập sử dụng hàm Green chưa nhiều đa dạng p ie rộng cho toán áp dụng phương pháp hàm Green quang học: nhiễu xạ, tán d oa nl w xạ,…; điện từ trường; vật lý thống kê,… oi lm ul nf va an lu z at nh z m co l gm @ an Lu n va ac th si 127 TÀI LIỆU THAM KHẢO Tiếng Việt [1] Vũ Văn Thanh, Nguyễn Nhật Khanh (2000) Phương trình đạo hàm riêng vật lý NXB Đại học Quốc gia TP Hồ Chí Minh Tiếng Anh [2] Duffy, D G (2004) Transform methods for solving partial differential equations Chapman and Hall/CRC lu an [3] Evans, G., Blackledge, J., & Yardley, P (2012) Analytic methods for partial va differential equations Springer Science & Business Media n gh tn to [4] Ruas, V (2016) Numerical Methods for Partial Differential Equations: An Introduction John Wiley & Sons p ie w [5] Heat equation Truy cập ngày 10 tháng 1, 2019 từ oa nl https://en.wikipedia.org/wiki/Heat_equationruy d [6] Laplace’s equation Truy cập ngày 11 tháng 1, 2019 từ an lu nf va https://en.wikipedia.org/wiki/Laplace%27s_equation oi lm ul [7] Wave equation Truy cập ngày tháng 1, 2019 từ https://en.wikipedia.org/wiki/Wave_equation z at nh z m co l gm @ an Lu n va ac th si