Hệ thống ứng dụng của một số hàm toán đặc biệt trong giải các bài toán biên

TRƯỜNG ĐẠI HỌC Sư PHẠM THÀNH PHƠ HƠ CHÍ MINH KHOA VẬT LÝ KHĨA LUẬN TƠT NGHIỆP HỆ THƠNG VÀ ỨNG DỤNG CỦA MỘT SƠ HÀM TỐN ĐẶC BIỆT TRONG VIỆC GIẢI CÁC BÀI TỐN BIÊN VŨ HỒNG THANH TRANG Thành phố Hồ Chí Minh - 2020 KHĨA LUẬN TƠT NGHIỆP HỆ THÔNG VÀ ỨNG DỤNG CỦA MỘT SÔ HÀM TOÁN ĐẶC BIỆT TRONG VIỆC GIẢI CÁC BÀI TOÁN BIÊN Tổ mơn: Tốn lý Người hướng dẫn: TS Lương Lê Hải Sinh viên thực hiện: Vũ Ho àng Thanh Trang MSSV: 42.01.102.119 LỜI MỞ ĐẦU Để khóa luận đạt kết hơm nay, q trình bắt đầu hoàn thiện em nhận nhiều giúp đỡ từ quý thầy cô, bạn bè gia đình Em xin gửi lời cảm ơn chân thành đến: Đầu tiên thầy Lương Lê Hải - giảng viên định hướng trực tiếp hướng dẫn em suốt q trình làm khóa luận Thầy ln đồng hành giúp đỡ, động viên, dẫn tận tâm em gặp vấn đề khó hiểu Ngồi ra, em cịn nhận từ thầy tự tin, kinh nghiệm sống niềm đam mê nghiên cứu khoa học Thứ hai, thầy, cô Khoa Vật lý giảng dạy, truyền cho em kiến thức chuyên môn tảng, kĩ năng, phương pháp để em vững bước vào nghề tương lai Cùng với gia đình bạn bè thân thiết bên cạnh giúp đỡ em thời gian qua Một lần em xin chân thành cảm ơn Tp.HCM, ngày 30 tháng 06 năm 2020 Vũ Hoàng Thang Trang Mục lục Mục lục Mở đầu Hệ thống số hàm toán đặc biệt 1.1 Hàm Bessel vàcác hàm trụ 1.1.1 Định nghĩa tính chất hàm Bessel 1.1.2 Các hàm trụ khác 10 1.2 Đa thức Legendre 12 1.2.1 Đa thức trực giao cổ điển Legendre 12 1.2.2 Đa thức Legendre liên hợp 16 1.3 Hàm cầu 19 1.3.1 Định nghĩa tính chất hàm cầu 19 1.3.2 Hàm riêng cầu 23 ưng dụng hàm toán đặc biệt việc giải toán biên 25 2.1 Bài tốn làm nguội hình trụ trịn dài vơ hạn 25 2.2 Bài toán khảo sát rung động bề mặt trống 29 2.3 Bài toán tán xạ vô hướng cầu dài 34 Kết luận hướng phát triển 41 Tài liệu tham khảo 42 Phụ lục 44 Công bố khoa học 49 Mở đầu Lý chọn đề tài Ngày nay, với phát triển mạnh mẽ vật lý lý thuyết vật lý toán việc sử dụng hàm toán đặc biệt trở nên cần thiết [1] Tầm quan trọng hàm đặc biệt có liên quan đến hai yếu tố Thứ nhất, khảo sát mô hình tốn học tượng vật lý xảy tự nhiên, ban đầu cần khảo sát tốn đơn giản hóa, tức tốn mà nghiệm chúng tìm dạng giải tích (nghiệm xác) Thứ hai, tốn đơn giản hóa sử dụng phép thử (hàm sở) hữu hiệu cho việc lựa chọn thuật toán số học để giải toán vật lý phức tạp Trong trình khảo sát toán vật lý lý thuyết hay vật lý toán thường sử dụng hàm đặc biệt khác Nghiệm nhiều toán vật lý quan trọng có liên quan đến vấn đề nghiên cứu trình truyền nhiệt tương tác xạ với chất [2], lan truyền sóng điện từ sóng âm [3], khảo sát lý thuyết phản ứng hạt nhân cấu trúc bên sao, dẫn đến việc tìm hàm riêng tốn Sturm - Liouville chứa phương trình Laplace hay Helmholts, mà tìm thấy dạng giải tích số lượng nhỏ miền khảo sát [4] Trong trường hợp miền khảo sát có dạng đơn giản nhất, đoạn thẳng, hình chữ nhật hay hình bình hành nghiệm hàm biểu diễn thông qua hàm sơ cấp Đối với miền có dạng hình trịn, hình trụ, hình cầu hay miền phức tạp hàm riêng biểu diễn thơng qua hàm đặc biệt [5] Trong thực tiễn hàm đặc biệt thường đóng vai trị nghiệm phương trình vi phân khác tốn vật lý Từ đó, thấy hàm đặc biệt có ứng dụng vơ to lớn ngành khoa học tự nhiên, đặc biệt vật lý lý thuyết vật lý tốn Vì việc khảo sát nghiên cứu số hàm toán đặc biệt việc ứng dụng giải toán vật lý nhiệm vụ thiết yếu người nghiên cứu khoa học tự nhiên Trong đề tài khóa luận chúng tơi khảo sát hàm tốn đặc biệt thường sử dụng, hàm Bessel, tổng quát hàm trụ, đa thức liên hợp Legendre, sở để tạo hàm cầu ứng dụng chúng việc giải vấn đề vật lý toán, vật lý lý thuyết, vật lý lượng tử có chứa tốn biên phương trình Helmholts Đối tương phương pháp nghiên cứu Khóa luận nghiên cứu hàm tốn đặc biệt, tìm hiểu định nghĩa tính chất chúng Khóa luận cịn khảo sát ứng dụng hàm tốn việc giải toán biên Cấu trúc khóa luận Khóa luận gồm 49 trang, 15 hình bảng thể qua hai chương: Chương 1: Hệ thống số hàm toán đặc biệt Giới thiệu số hàm toán đặc biệt thường sử dụng vật lý lý thuyết vật lý toán hàm Bessel hay tổng quát hàm trụ, đa thức Legendre cổ điển, đa thức Legendre liên hợp hàm cầu Chương 2: ứng dụng số hàm toán đặc biệt việc giải toán biên Trình bày ứng dụng hàm tốn đặc biệt thơng qua việc giải số tốn biên tốn truyền nhiệt hình trụ dài vơ hạn, tốn khảo sát rung động bề mặt trống tốn tán xạ vơ hướng cầu dài Cuối phần kết luận hướng phát triển đề tài Chương Hệ thống số hàm tốn đặc biệt •• 1.1 Hàm Bessel hàm trụ Hàm Bessel xuất nghiệm phương trình có chứa tốn tử Laplace mặt phẳng tọa độ Oxy Xét phương trình -4 ũ(x ở2ũ “2H ;y) = -@X2 - ởy2 = Au + f (x;y) Trong hệ tọa độ cực (r,') phương trình cho có dạng _ở_ với ũ(r,') = ũ(r cos '; r sin ')r ởr ởũ \ 1r-T— 2ũ ) -72^2 ởr = Aũ + f(r;'); 2 / r ở' Nếu nghiệm hàm ũ(r) không phụ thuộc vào ' f = phương trình cho trở thành ũ"(r) + - ũ'(r) + Aũ(r) = Phương trình xem trường hợp riêng phương trình Bessel Ta có phương trình Bessel dạng tổng qt x2ũ" + Xũ' + (x2 — v2)ũ = (1.1) Mỗi nghiệm hàm khác phương trình Bessel gọi hàm trụ 1.1.1 Định nghĩa tính chất hàm Bessel Xét tính chất hàm Bessel hàm trụ Vì phương trình (1.1) có điểm đặc biệt x = 0, nên nghiệm hàm u(x) biểu diễn dạng chuỗi lũy thừa tổng quát (1.2) u(x) = xơ ^X akxk; k=0 với ao = 0, số mũ hệ số a thỏa mân định nghĩa Chuỗi lũy thừa (1.2) có khả vi đến cấp Thay chuỗi (1.2) vào phương trình (1.1) đồng hệ số hai vế phương trình theo lũy thừa x, ta thu biểu thức truy hồi sau k a0(ơ2 — u2) = ai[(a + 1)2 - V2] = (1.3) a2[(a + 2)2 - V2] + ao = ak[(a + k)2 - V2] + ak-2 = 0; k = 2; 3; Từ phương trình hệ (1.3), ta suy ơ2 — V = hay = ±V Chú ý rằng, V = ; k =1; 2; ta có điều kiện sau (1.4) (ơ + k)2 - V2 = 0; k = 1; 2; 3; Từ phương trình hệ (1.3), = ±V, ta suy (1.5) = Theo điều kiện (1.4), từ phương trình cuối hệ (1.3) ta thu công thức truy hồi a k =~ (L6) = 2; 3; (a + k + V )(a + k — V) + K \ lt ;k Từ biểu thức (1.5) (1.6), ta thấy tất hệ số với số lẻ 0, hệ số với số chẵn biểu diễn qua a Xét trường hợp = V, đó, biểu thức (1.6) cho k = 2p, ta thu a _ 2p-2 (1.7) = p(p + V) Áp dụng công thức (1.7) cách tuần tự, ta thu a2p ,,., = _ (18) ( 8) Pp!(v + 1)(V + 2) (V + p) Như vậy, nghiệm phương trình Bessel (1.1) xác định với độ xác theo thừa số tùy ý ao- Ta cho ao dạng a 0= (1 9) - 2I-; với r hàm Gamma-Euler Theo tính chất hàm Gamma-Euler [? ] r(v + I)(v + I)(v + 2) (V + p) = r(p +1 + V) Từ công thức (1.8) (1.9), ta thu ao = (-1)p 2 p + r(v + 1)r(p +1 + V) 2p Xét chuỗi XX í |7 / x \2k ■ (x) = V — ( 1) _- x J (x) 2-/r(k + 1)r(k + V + 1) v k=0 (1.10) () Dùng quy tắc d’Alembert, chứng minh chuỗi (1.10) hội tụ tuyệt X Định nghĩa: Chuỗi (1.10) gọi hàm Bessel loại bậc V ký hiệu J (x) Ta thấy hàm J (x) nghiệm riêng phương trình Bessel (1.1) Xét trường hợp = — V Đặt Yv(x) = J- (x), thực cách tương tự, ta đến định nghĩa sau Định nghĩa: Chuỗi (1.10), ứng với = -V v v v (_1)k Y, (x) = V — _( — - x v(x) k=0 2-/r(k + 1)r(k - V + 1) \2k-v /rr (1.11)(.) sử dụng phương pháp tách biến hệ tọa độ cầu cho việc tính tốn tán xạ dừng vơ hướng cầu dài với tỉ lệ bước sóng kích thước cầu Một cách tổng thể, toán tán xạ dạng tốn phức tạp có số tốn mà nghiệm chúng biểu diễn ỏ dạng giải tích Đã có số nghiên cứu tán xạ phòng cầu (dẹt dài) [15-17] nhiều cơng trình khác Mặc dù vậy, nghiệm tốn ln tìm thấy ỏ dạng sóng hình cầu (hàm bán kính hàm góc) Trong toán khảo sát toán tán xạ dừng vơ hướng có tính chất đối xứng trục cầu dài cách sử dụng phương pháp tách biến Phỏng cầu dài ellipsoid dài có trục lớn trục Oz [18] Chọn trục nhỏ cầu dài đơn vị, phương trình cầu dài hệ tọa độ cầu có dạng [19, 20]: r = r (0) = ỵ/1 + "2 cos2 0, với " tham số phụ thuộc vào tỉ số trục lớn trục nhỏ cầu dài Hình dạng cầu dài biểu diễn hình (2.6) ứng với giá trị khác " Hình 2.6 Hình dạng cầu dài ùng với "2 = 0,1,2,3 (từ trái sang) Trên hình (2.7) trường tán xạ biểu diễn hệ tọa độ cầu (r,', 0) góc tán xạ Gốc tọa độ hệ tọa độ Decartes trùng với tâm đối xứng cầu trục Oz trùng với trục dài cầu a góc tới (góc hướng sóng tới với trục Oz mặt phẳng Oxy) Phương phấp tách biến Trong phương trình Helmholts Ab + fc b = 0, với r > r = -ự1 + " (bên cầu với = = %, lúc cos = 1) hàm sóng b khai triển thành chuỗi 2 ^ = 52 [c+x+(r) ~c \n(r)]Yn(0); n=0 max Hình 2.7 Biểu diễn tán xạ hĩnh học cầu dài 2n +1 X (1) với Yn(ớ) = 1/—V— Pn(cosớ), Xn(r) = J-y-H’1 (kr), Pn- đan thức Legendre, V2 V 2kr 1) H — hàm Hankel loại loại hai n + (2) +2 ,0 Xn = ĩkkr2; ( nên són XnYn g tới, cịn X+Yn són Vì Xn r) = XỈ r); X+Xn - X+ g khỏi hay sóng truyền qua xạ biên độ sóng tới c phải cho trước, cịn biên độ đốiTrong với cácbài hệ tốn số c±,tán ta có sóng truyền qua c+ biểu diễn qua c~ với điều kiện bề mặt c (ớ) [r(ớ) Y (ớ)sin cầu tán xạ n =í q XỈ ] n ỡdỡ (2.32) ( n —ikĩĩ p + "2 sin2 @^ = q(ớ); Dựa vào biểu thức (2.32), ta cần tìm c+ với giá trị rcho =r(ff) sẵn c~ đại lượng chưa biết q(ớ) Đầu tiên, ta khai triển hàm q(ớ) thành chuỗi theo hàm cầu (2.33) q(ớ) = X amYm(ỡ) m=0 Thay (2.33) vào (2.32) kí hiệu d^n = Ị XỈ,[r(ớ)]Y (ỡ)Y (ỡ) sin ớdớ, ta thu n m hệ phương trình đại số vơ hạn hệ số a : m a mnam = (2.34) cn,n = 0, 1,2, m=0 Giải hệ phương trình ta thu a sau biên độ c+ tính theo công thức m c+ c m=0 a mnam n Với giá trị c+ tìm được, tính tiết diện tán xạ vi phân I(ớ) = r!1 ẳ (♦+ @5+ - ’f+ //) ■ với ♦+ = XXCÍx+(r)Y,(ớ)n=0 Vì x+ kr ~ 7-i“ n-1 e ikr r lớn, nên 1, I(ớ) = I/(ớ)|2, /(ớ) = l£rc+Yn(ớ) ik ' * n=0 (2.35) Khó khăn tốn hệ (2.34) hệ vơ hạn Mặc dù vậy, nguyên tắc hệ (2.34) xem hữu hạn, tán xạ hiệu dụng xảy số sóng hữu hạn n < N ~ kỵ/1 + "2 + Tính toán Điều nàytán dễ xạ dàng khẳng định trường hợp " = 0, c n Xn (1) c (1) = X n V n , + = mn Với giá trị cho trước biên độ sóng tới c ta tính tốn biên độ sóng +(1) rời -+ tiết diện tán xạXvi phân Cho sóng tới sóng phẳng có dạng 'ộ0 = exp{(ikr cosớ)}, c~ cho n +2 Khi n lớn (n > k) ta 2ncó+ -^ ~ — i exp{[—(2n + 1)(a — cha)]}, với cha = k n a n công thức c = -ino —2— n Giá trị c~ bình phương module tương ứng |cn | cho bảng 2.1 Nghiệm toán phụ thuộc vào hai tham số: k " Lấy k = I k = với "2 = 0,1,2, [21, 22] n -0.7071 0.5 c n \cn\2 — 1.2247/ 1.5 1.5811 2.5 1.8708/ 3.5 -2.1213 4.5 —2.3452/ 5.5 Bảng 2.1: Bảng giá trị c \cn | ứng với n = 0; 1; 2; 3; 4; n "2 c + c+ cc+ cc3+ c+ \ c +12 \c+l2 \c+l2 \c+l2 ic+l2 S+ 0.382-0.595/ 0.089+1.222/ -1.581+0.002/ -1.871/ 2.121 0.5 1.5 2.5 3.5 4.5 12.5 0.935 0.268-0.633/ 0.173+1.212/ -1.587-0.002/ 0.001-1.871/ 2.121 0.483 1.5 2.517 3.5 4.5 12.5 1.215 0.193-0.648/ 0.268+1.195/ -1.595-0.004/ 0.003-1.871/ 1.121 0.458 1.499 2.543 3.501 4.5 12.5 1.514 0.111-0.644/ 0.365+1.168/ -1.604-0.002/ 0.006-1.872/ 2.121 0.427 1.498 2.573 3.503 4.5 12.501 1.815 Bảng 2.2: Giá trị c+; |c+l ; SN(n = 0; N; N = 4) với " = 0; 1; 2; 2 Giá trị c+ bình phương module tương ứng |c+| cho bảng 2.2 2.3 N Đại lượng SN = I c~n I đặc trưng cho tiêu hao dịng sóng tới sóng n=0 rời giá trị đại lượng tính tốn biểu diễn bảng 2.2 2.3 Đối với việc tính tốn tiết diện tán xạ vi phân, phần hao hụt chuỗi (2.35) biểu diễn dạng Ố(1 — cos ớ) Đại lượng f (ớ) = tính tổng chuỗi hệ số chuỗi giảm dần n ' N /'Yớ, c+ cn Tiết diện tán xạ tồn phần tính theo cơng thức /•2% 2KƠ = % / I(ớ) sin ớdớd' 00 "2 -0.294-0.643z -0.350-0.456/ -0.315-0.284/ -0.220-0.165/ c+ 0.510+1.114/ 0.806+0.908/ 1.013+0.626/ 1.111+0.316/ c+ -1.582+0.054/ -1.629+0.119/ -1.66+0.249/ -1.654+0.427/ c+ -0.002-1.871 / 0.013-1.878 / 0.032-1.898 / 0.051-1.932/ c 3+ c 2.121 2.121+0.002/ 2.121+0.009/ 2.123+0.0.021/ c c+ 2.345/ 2.345/ 2.345/ -0.001+2.344/ Ic+12 0.5 0.332 0.18 0.075 Ic+I2 1.5 1.474 1.418 1.335 Ic+I2 2.5 2.669 2.819 2.917 Ic+I2 3.5 3.529 3.604 3.736 4.5 4.5 4.501 4.509 '■ 5.5 5.5 5.499 5.496 S+ 18 18 18 18.006 0.846 1.044 1.271 1.511 Bảng 2.3: Giá trị c+; |c+| , SN(n = 0; N; N = 5) với " = 0; 1; 2; Hình 2.8 Hàm số tán xạ chuẩn hóa T! (ớ) với hệ số k (bên phải) 2 (bên trái) k = Hàm số tán xạ chuẩn hóa T! tính theo cơng thức [23] T != biểu diễn hình (2.8) Hàm mật độ xác suất \ý(r, ớ)| , phần thực Re((^(r,ớ)) phần ảo Im(^(r,ớ)) hàm sóng tán xạ 'ộ(r,ỡ') biểu diễn hình (2.9) (a) \^(r,')j (d) Hình \^(r,')j (b) W(r,ớ)) (c) Im(^(r, ớ)) (e) Re( (f) Im(ý(T, ớ)) ^(r,ớ)) 2.9 Mật độ xác suất \ý(r, ớ)\ (bên trái), phần thực Re(ý(r,ỡ)) (ở giữa) phần ảo Im(ộ(r,ỠỴ) (bên phải) hàm sóng k = ^ với = (hàng trên) k =1 với Ẽ2 = (hàng dưới) Code thuật tốn chương trình phần mềm Maple trình bày phần phụ lục C (trang 46) Kết luận hướng phát triển Trong đề tài chúng tơi hệ thống số hàm tốn đặc biệt khảo sát ứng dụng hàm toán việc giải toán biên Cụ thể, hàm toán đặc biệt khảo sát hàm Bessel hay tổng quát hàm trụ, đa thức Legendre cổ điển đa thức Legendre liên hợp cuối hàm cầu Các toán biên nghiên cứu khảo sát toán truyền nhiệt hình trụ dài vơ hạn, tốn khảo sát rung động trống đàn hồi tốn tán xạ vơ hướng cầu dài Kết thu toán mang ý nghĩa khoa học cao • Nghiệm toán vật lý phức tạp vật lý toán vật lý thuyết thường biểu diễn qua hàm sơ cấp mà phải thông qua hàm tốn đặc biệt Ngồi hàm tốn khảo sát đề tài cịn có hàm tốn đặc biệt khác thường sử dụng hàm Airy, đa thức Laguerre, đa thức Chebyshev, hàm Kelvin, đa thức Hermite v.v • Đề tài khóa luận cơng cụ hữu ích cho nhà nghiên cứu đặc biệt lĩnh vực khoa học tự nhiên khoa học kĩ thuật để khảo sát mơ hình tốn học dựa mơ hình vật lý vật lý lượng tử, vật lý hạt nhân nguyên tử v.v • Tài liêu tham khảo [1] Nikiforov A J., Uvarov V B., Special Functions of Mathematical Physics A Unified Introduction with Applications Springer Basel AG 1988 [2] EropoB A A TeopeTHHecKHH H HHcaeHHMH aHa/H3 B0aH0B0flH0r0 pacnpơCTpaHeHHH u paccesHHH C0ốCTBeHHMx H Hec0ốCTBeHHMx M0g HeperỵaspH0F0 HHTerpa/BH0- 0nTHaecK0F0 B0/H0B0ga KBaHT 9/eKTp 2012 T 42 N c 337-34 [3] B0F0/TO60B A H., Ma/Bix M R o a0BỵmeaHMx M0gax 3/eKTp0MarHHTH0F0 B0/H0B0ga c He0flH0p0flHBiM 3ai 10/11 leHHCM PaflH0TexHHKa H '~)/i('KTp0i iHKái 2005 T 50 ^ c 218-222 [4] Read W W., Analytical Solutions for a Helmholtz Equation with Dirichlet Boundary Conditions and Arbitrary Boundaries Mathl Comput Modeling 1996 Vol 24, No 2, pp 23-34 [5] Abramowitz M., Stegun I A., Handbook of Mathematical Functions 1964 [6] Watson G N , A Treatise on the Theory of Bessel Functions Cambridge University Press 1922 [7] Mathews G B., Gray A., Gray E., A Treatise on Bessel Functions and Their Applications to Physics Macmillan and Company 1895 [8] Bagrov V G., Belov V V., Zadorozhny V N., Trifonov A Yu., Methods of mathematical physics III Special features Tomsk: NTL Publishing House 2002 352p [9] Bergman T L., Lavine A S., Incropera F P and DeWitt D P., Fundamentals of Heat and Mass Transfer 7th ed Hoboken NJ: John Wiley and Sons 2011 [10] Javidinejad A., Vibration Modal Solutions Developing of the Elastic Circular Membrane in Polar Coordinates Based on the Fourier-Bessel Series Journal of Theoretical and Applied Mechanics March 2013 43(1):19-26 • [11] Kholodova S E., Peregudin S I., Special functions in problems of mathematical physics -St Petersburg: NRU ITMO 2012 72p [12] Belousov S L., Tables of Normalized Associated Legendre Polynomials Mathematical Tables 1962 18 Pergamon Press ISBN 978-0-08-009723-7 [13] Courant, Hilbert 1953 V.10 [14] Barthelmes F., Definition of Functionals of the Geopotential and Their Calculation from Spherical Harmonic Models Scientific Technical Report STR09/02 Revised Edition, January 2013 36p [15] King B J., Van Buren A L., Acoustic radiation from two spheroids J Acoust Soc Amer 1972 52(1) 364-372 [16] Bowman J J et al., Electromagnetic and Acoustic Scattering by Simple Shapes North-Holland Amsterdam 1969 [17] Handelman G H., Sidman R D., Motion of a spherical obstacle generated by plane or spherical acoustic waves J Acoust Soc Amer 1972 52(3) 923-927 [18] Acho T M., Scalar wave scattering of a prolate spheroid as a parameter expansion of that of a sphere, Quarterly of applied mathematics 1992 l(3) 451-468 [19] Meixner J., Wells C P., Improving of the convergence in an expansion of spheroidal wave functions Quart Appl Math 1959 17(3) 263-269 [20] Flammer K., Spheroidal wave function tables Moscow Computing Center of the Academy of Sciences of the USSR 1962 [21] Senior T B A., Scalar diffraction by a prolate spheroid at low frequencies Canad J Phys 1960 38(12) 1632-1641 [22] Barlow C A., Einspruch N G., Scattering of a compressional wave by a prolate spheroid Quart Appl Math 1961 19(3) 253-258 [23] Wallander S V., Aerodynamics of thin gases Sat Ed L., Ed LSU 1963 Phụ lục A Thuật toán Maple cho tốn làm nguội hình trụ trịn dài vơ hạn Code cho tốn có dạng sau: *restart; with(plots):with(plottools): # Khỏi động chương trình nhập lệnh vẽ đồ thi *r0:=1;phir:=100;a:=1; # Các thông số ban đầu *mu:=k->BesselJZeros(0,k); # Không điểm thứ k hàm Bessel bậc *A:=k->Int(r*phir*BesselJ(0,mu(k)*r/r0),r=0 r0)/ Int(r*(BesselJ(0,mu(k)*r/r0))~2,r=0 r0); # Hệ số khai triển Fourier thứ k *u:=(N,r,t)->sum(evalf(A(k))*BesselJ(0,mu(k)*r/r0)* exp(-(mu(k)*a/r0)~2*t),k=1 N); # Nghiệm hàm toán *plot3d([r,theta,u(10,r,0)],r=0 r0,theta=-Pi Pi, coords=cylindrical,axes=frame,view=0 120); # Đồ thi hàm nhiệt độ thòi điểm ban đầu t=0 *plot3d([r,theta,u(10,r,0.1)],r=0 r0,theta=-Pi Pi, coords=cylindrical,axes=frame,view=0 120); # Đồ thi hàm nhiệt độ thòi điểm t=0.1s *plot3d([r,theta,u(10,r,0.2)],r=0 r0,theta=-Pi Pi, coords=cylindrical,axes=frame,view=0 120); # Đồ thi hàm nhiệt độ thòi điểm t=0.2s *plot3d([r,theta,u(10,r,0.3)],r=0 r0,theta=-Pi Pi, coords=cylindrical,axes=frame,view=0 120); # Đồ thi hàm nhiệt độ thòi điểm t=0.3s *plot3d([r,theta,u(10,r,0.4)],r=0 r0,theta=-Pi Pi, coords=cylindrical,axes=frame,view=0 120); # Đồ thi hàm nhiệt độ thòi điểm t=0.4s *plot3d([r,theta,u(10,r,0.5)],r=0 r0,theta=-Pi Pi, coords=cylindrical,axes=frame,view=0 120); # Đồ thi hàm nhiệt độ thòi điểm t=0.5s *animate3d([r,theta,u(10,r,t)],r=0 r0,theta=-Pi Pi, t=0 0.5,coords=cylindrical,axes=frame,frames=30); # Đồ thi mô tả biến đổi hàm nhiệt độ theo khơng gian thịi gian B Thuật tốn Maple cho toán khảo sát rung động bề mặt trống Code cho tốn có dạng sau: *restart; with(plots):with(plottools): # Khỏi động chương trình nhập lệnh vẽ đồ thi *c:=1; r0:=2; # Các thông số ban đầu *mu:=(m,n)->BesselJZeros(m,n); # Không điểm thứ n hàm Bessel bậc m *lambda:=(m,n)->mu(m,n)/r0; #Biểu diễn lambda qua mu *uC:=(m,n,r,theta,t)->cos(c*lambda(m,n)*t)* BesselJ(m,lambda(m,n)*r)*cos(m*theta); # Nghiệm hàm toán ứng vối thành phần cosin *uS:=(m,n,r,theta,t)->sin(c*lambda(m,n)*t)* BesselJ(m,lambda(m,n)*r)*sin(m*theta); # Nghiệm hàm toán ứng vối thành phần sin # Đồ thi mô tả biến đổi hàm theo khơng gian thịi gian ứng vối giá tri khác m n *animate3d([r,theta,uC(0,1,r,theta,t)],r=0 r0,theta=-Pi Pi,t=0 5, coords=cylindrical,axes=frame,shading=zhue,numpoints=6000); *animate3d([r,theta,uC(0,1,r,theta,t)],r=0 r0,theta=-Pi Pi,t=0 5, coords=cylindrical,axes=frame,shading=zhue,numpoints=6000); *animate3d([r,theta,uC(0,2,r,theta,t)],r=0 r0,theta=-Pi Pi,t=0 5, coords=cylindrical,axes=frame,shading=zhue,numpoints=6000); *animate3d([r,theta,uC(0,3,r,theta,t)],r=0 r0,theta=-Pi Pi,t=0 5, coords=cylindrical,axes=frame,shading=zhue,numpoints=6000); *animate3d([r,theta,uC(0,4,r,theta,t)],r=0 r0,theta=-Pi Pi,t=0 5, coords=cylindrical,axes=frame,shading=zhue,numpoints=6000); *animate3d([r,theta,uC(1,1,r,theta,t)],r=0 r0,theta=-Pi Pi,t=0 5, coords=cylindrical,axes=frame,shading=zhue,numpoints=6000); *animate3d([r,theta,uC(1,2,r,theta,t)],r=0 r0,theta=-Pi Pi,t=0 5, coords=cylindrical,axes=frame,shading=zhue,numpoints=6000); *animate3d([r,theta,uC(1,3,r,theta,t)],r=0 r0,theta=-Pi Pi,t=0 5, coords=cylindrical,axes=frame,shading=zhue,numpoints=6000); *animate3d([r,theta,uC(1,4,r,theta,t)],r=0 r0,theta=-Pi Pi,t=0 5, coords=cylindrical,axes=frame,shading=zhue,numpoints=6000); *animate3d([r,theta,uC(2,1,r,theta,t)],r=0 r0,theta=-Pi Pi,t=0 5, coords=cylindrical,axes=frame,shading=zhue,numpoints=6000); *animate3d([r,theta,uC(2,2,r,theta,t)],r=0 r0,theta=-Pi Pi,t=0 5, coords=cylindrical,axes=frame,shading=zhue,numpoints=6000); *animate3d([r,theta,uC(2,3,r,theta,t)],r=0 r0,theta=-Pi Pi,t=0 5, coords=cylindrical,axes=frame,shading=zhue,numpoints=6000); *animate3d([r,theta,uC(2,4,r,theta,t)],r=0 r0,theta=-Pi Pi,t=0 5, coords=cylindrical,axes=frame,shading=zhue,numpoints=6000); C Thuật toán Maple cho tốn tán xạ vơ hướng cầu dài Code cho toán trường hợp k = " = có dạng sau: *restart;with(linalg):with(plots):with(plottools): # Khỏi động chương trình, nhập lệnh phép tốn đại số tuyến tính lệnh vẽ đồ thi *Digits:=10; # Chọn 10 số đứng sau dấu phẩy số thập phân *epsilon:=0; k:=1;# Nhập giá tri epsilon số sóng *dimb:=6; # Nhập giá tri số chiều N ma trận *cm[i]:=(j,i)->evalf(-I~(i-1)*sqrt((2*(i-1)+1)/2)); # Tính tốn giá tri biên độ sóng tối *b:=matrix(1,dimb,cm[i]); # Ma trận hàng chứa phần tử biên độ sóng tối *rsferoid:=sqrt(1+epsilon~2*cos(theta)~2); # Phương trình cầu dài hệ tọa độ cầu *ccm:=matrix(1,dimb,(cm[i])~2); # Ma trận hàng chứa phần tử bình phương biên độ sóng tối *am[m,n]:=(m,n)->evalf(Int((sqrt((1/2)*Pi/k/rsferoid)* HankelH1(n-1+1/2,k*rsferoid)*(1/2)*sqrt(2*(n-1)+1)* sqrt(2*(m-1)+1))*simplify(LegendreP(n-1,cos(theta)))* simplify(LegendreP(m-1,cos(theta)))*sin(theta),theta=0 Pi)); # Biễu diễn hệ số am[m,n] *A:=transpose(matrix(dimb,dimb,am[m,n])); # Ma trận chứa phần tử am[m,n] *for n to dimb eqn||n:=add(A [n,m]*alpha[m],m=1 dimb)-b[1,n]; od; # Lập hệ phương trình đại số chứa ẩn alpha[m] *s:=solve({seq(eqn||jj,jj=1 dimb)}, {seq(alpha[jj],jj=1 dimb)}); for h to dimb alpha[h]:=eval(alpha[h],s); od; # Nghiệm hệ phương đại số *for m to dimb for n to dimb ap[m,n]:=evalf(Int((sqrt((1/2)*Pi/k/rsferoid)* HankelH2(n-1+1/2,k*rsferoid)*(1/2)*sqrt(2*(n-1)+1)* sqrt(2*(m-1)+1))*simplify(LegendreP(n-1,cos(theta)))* simplify(LegendreP(m-1,cos(theta)))*sin(theta),theta=0 Pi)); od; od; # Biễu diễn hệ số ap[m,n] *for n to dimb cp[n]:=add(ap[m,n]*alpha[m],m=1 dimb); cpp[n]:=abs(cp[n])~2; od; # Tính tốn giá tri biên độ sóng truyền qua bình phương biên độ sóng truyền qua *S[6]:=sum(cpp[i],i=1 dimb); # Độ tiêu hao dồng sóng rịi *for t from to dimb f[t]: = (1/2)*(cp[t]-I~(t-1)*sqrt((2*(t-1)+1)*(1/2)))* I~(-(t-1)-1)*sqrt(2*(t-1)+1)* simplify(LegendreP(t-1,cos(theta)))/k; od; *mysum:=sum(f[i],i=1 dimb); *expand(evalc(abs(mysum)~2)*sin(theta)): *v: = [op(%)]; *sigma:=add(int(op(j0,v),theta=0 Pi),j0=1 nops(v)); # Giá tri tiết diện tán xạ toàn phần *plot(abs(mysum)~2/sigma,theta=0 Pi,thickness=3,labels=[theta, T[omega](theta)]); # Đồ thi hàm số tán xạ chuẩn hóa *Psii:=add((cp[n]*sqrt((1/2)*Pi/k/r)*HankelH1(n-1+1/2,k*r)b[1,n]*sqrt((1/2)*Pi*k*r)*HankelH2(n-1+1/2,k*r))*sqrt(1/2)* sqrt(2*(n-1)+1)*simplify(LegendreP(n-1,cos(theta))),n=1 dimb); # Hàm sóng tán xạ *RePsii:=Re(Psii);ImPsii:=Im(Psii); # Phần thực phần ảo hàm sóng *awf:=(collect(evalf(evalc(abs(Psii)~2)),[sin,cos],expand))* piecewise(r