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 TỐN ĐẶC BIỆT TRONG VIỆC
GIẢI CÁC BÀI TỐN BIÊN

Tổ bộ 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 quả như hơm nay, trong q trình bắt đầu và hồn thiện
em đã nhận được rất nhiều sự giúp đỡ từ quý thầy cơ, bạn bè và gia đình. Em xin
gửi lời cảm ơn chân thành đến:

Đầu tiên là thầy Lương Lê Hải - giảng viên định hướng và trực tiếp hướng dẫn
em trong suốt q trình làm khóa luận. Thầy ln đồng hành giúp đỡ, động viên,
chỉ dẫn tận tâm khi em gặp vấn đề khó hiểu. Ngồi ra, em cịn nhận được từ thầy
sự tự tin, kinh nghiệm sống và niềm đam mê nghiên cứu khoa học.
Thứ hai, các thầy, cô trong Khoa Vật lý đã giảng dạy, truyền cho em những
kiến thức chuyên môn nền tảng, kĩ năng, phương pháp để em có thể vững bước
vào nghề trong tương lai.
Cùng với đó là gia đình và bạn bè thân thiết luôn bên cạnh và giúp đỡ em trong
thời gian qua.
Một lần nữa 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

2

Mở đầu

3

1 Hệ thống một số hàm toán đặc biệt
5
1.1 Hàm Bessel vàcác hàm trụ ......................................................................... 5
1.1.1 Định nghĩa và các tính chất
cơ bản của hàm Bessel.............
6

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 và các tính chất
cơ bản của hàm cầu...................... 19
1.3.2 Hàm riêng của quả cầu.................................................................... 23
2 ưng dụng của các hàm toán đặc biệt trong việc giải các bài toán
biên
25
2.1 Bài toán về sự làm nguội của hình trụ trịn dài vơ
hạn........................ 25
2.2 Bài toán khảo sát sự rung động của bề mặt trống..................................... 29
2.3 Bài tốn tán xạ vơ hướng trên phỏng cầu dài........................................... 34
Kết luận và 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

2


Mở đầu
1. Lý do chọn đề tài
Ngày nay, với sự phát triển mạnh mẽ của vật lý lý thuyết và vật lý toán việc
sử dụng các hàm toán đặc biệt đã trở nên rất cần thiết [1]. Tầm quan trọng của
các hàm đặc biệt này có liên quan đến hai yếu tố cơ bản. Thứ nhất, khi khảo sát
các mô hình tốn học của các hiện tượng vật lý xảy ra trong tự nhiên, ban đầu
chúng ta cần khảo sát những bài tốn đã được đơn giản hóa, tức là những bài tốn
mà nghiệm của chúng có thể tìm được ở dạng giải tích (nghiệm chính xác). Thứ
hai, những bài tốn đã được đơn giản hóa này có thể được sử dụng như là phép
thử (hàm cơ sở) hữu hiệu cho việc lựa chọn những thuật toán số học để giải quyết
những bài toán vật lý phức tạp hơn.
Trong quá trình khảo sát những bài tốn vật lý lý thuyết hay vật lý toán chúng
ta thường sử dụng những hàm đặc biệt khác nhau. Nghiệm của nhiều bài toán vật
lý quan trọng có liên quan đến các vấn đề như nghiên cứu các quá trình truyền
nhiệt và tương tác bức xạ với các chất [2], sự lan truyền của các sóng điện từ và
sóng âm [3], khảo sát lý thuyết phản ứng hạt nhân và cấu trúc bên trong của các
sao, dẫn đến việc tìm hàm riêng của bài tốn Sturm - Liouville chứa phương trình
Laplace hay Helmholts, mà có thể được tìm thấy ở dạng giải tích chỉ đối với một
số lượng nhỏ các miền khảo sát [4].
Trong các trường hợp các miền khảo sát có dạng đơn giản nhất, như là đoạn
thẳng, hình chữ nhật hay hình bình hành thì các nghiệm hàm này được biểu diễn
thơng qua các hàm sơ cấp cơ bản. Đối với những miền có dạng hình trịn, hình trụ,
hình cầu hay những miền phức tạp hơn thì các hàm riêng được biểu diễn thông
qua các hàm đặc biệt [5]. Trong thực tiễn những hàm đặc biệt thường đóng vai trị
như là nghiệm của những phương trình vi phân khác nhau của các bài tốn vật lý.
Từ đó, có thể thấy các hàm đặc biệt có ứng dụng vơ cùng to lớn trong các ngành

khoa học tự nhiên, đặc biệt là vật lý lý thuyết và vật lý tốn. Vì vậy việc khảo sát
và nghiên cứu một số hàm toán đặc biệt trong việc ứng dụng giải các bài toán vật
lý là một nhiệm vụ thiết yếu của người nghiên cứu khoa học tự nhiên.
Trong đề tài khóa luận này chúng tơi sẽ khảo sát những hàm toán đặc biệt

3


thường được sử dụng, như hàm Bessel, tổng quát hơn là hàm trụ, đa
thức
liên
hợp
Legendre, là cơ sở để tạo ra hàm cầu và những ứng dụng của chúng trong
việc
giải
quyết các vấn đề trong vật lý toán, vật lý lý thuyết, vật lý lượng tử có chứa
bài
tốn biên đối với phương trình Helmholts.

2. Đối tương và phương pháp nghiên cứu
Khóa luận nghiên cứu các hàm tốn đặc biệt, tìm hiểu định nghĩa và tính chất
của các chúng.
Khóa luận cịn khảo sát ứng dụng của các hàm toán này trong việc giải các bài
tốn biên.
3. Cấu trúc khóa luận
Khóa luận gồm 49 trang, 15 hình và 3 bảng được thể hiện qua hai chương:
Chương 1: Hệ thống một số hàm toán đặc biệt
Giới thiệu một số hàm toán đặc biệt thường được sử dụng trong vật lý lý thuyết
và vật lý toán như hàm Bessel hay tổng quát hơn là hàm trụ, đa thức Legendre cổ
điển, đa thức Legendre liên hợp và hàm cầu.

Chương 2: ứng dụng của một số hàm toán đặc biệt trong việc giải các bài tốn
biên

Trình bày ứng dụng của các hàm tốn đặc biệt thơng qua việc giải một số bài
toán biên như bài toán truyền nhiệt trong một hình trụ dài vơ hạn, bài tốn khảo
sát sự rung động của bề mặt trống và bài toán tán xạ vô hướng trên phỏng cầu
dài.
Cuối cùng là phần kết luận và hướng phát triển của đề tài.

4


Chương 1

Hệ thống một số hàm tốn
đặc biệt
••
1.1. Hàm Bessel và các hàm trụ
Hàm Bessel xuất hiện trong nghiệm của các phương trình có chứa tốn tử
Laplace 4 trong 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,') thì phương trình đã cho có dạng
1
_ở_
với ũ(r,') = ũ(r cos '; r sin ')r .ởr

ởũ
\
1r-T—
ở 2ũ
) -72^2 ởr
= Aũ + f(r;');
2
2
/ r ở'

Nếu nghiệm hàm ũ(r) không phụ thuộc vào ' và f = 0 thì phương trình đã cho
trở thành

ũ"(r) + - ũ'(r) + Aũ(r) = 0.

Phương trình này được xem như là trường hợp riêng của phương trình Bessel.
Ta có phương trình Bessel ở dạng tổng quát
x2ũ" + Xũ' + (x2 — v2)ũ = 0.

(1.1)



Mỗi nghiệm hàm khác 0 của phương trình Bessel được gọi là hàm trụ.

1.1.1.

Định nghĩa và các tính chất cơ bản của hàm Bessel

Xét các tính chất cơ bản của hàm Bessel và các 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) của nó có thể được 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ũ ơ và các hệ số a thỏa mân định nghĩa. Chuỗi lũy thừa (1.2) có
khả vi đến cấp bất kì. Thay chuỗi (1.2) vào phương trình (1.1) và đồng nhất hệ số
hai vế của phương trình theo lũy thừa của x, ta thu được các biểu thức truy hồi
sau
k

a0(ơ2 — u2) = 0
ai[(a + 1)2 - V2] = 0

(1.3)

a2[(a + 2)2 - V2] + ao = 0

ak[(a + k)2 - V2] + ak-2 = 0; k = 2; 3;...

Từ phương trình đầu tiên của hệ (1.3), ta suy ra ơ2 — V = 0 hay ơ = ±V. Chú ý
2

rằng, khi V = ; k =1; 2;... thì ta có điều kiện sau
(1.4)

(ơ + k)2 - V2 = 0; k = 1; 2; 3; ...

Từ phương trình đầu tiên của hệ (1.3), khi ơ = ±V, ta suy ra
(1.5)

ai = 0

Theo điều kiện (1.4), từ phương trình cuối cùng của hệ (1.3) ta thu được 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) và (1.6), ta thấy rằng tất cả các hệ số với chỉ số dưới lẻ đều
bằng 0, còn các hệ số với chỉ số dưới chẵn có thể được biểu diễn qua a . Xét trường
0


hợp ơ = V, khi đó, trong biểu thức (1.6) cho k = 2p, ta thu được
a

_
2

2p-2

(1.7)

= 2 p(p + V) .
2

Áp dụng công thức (1.7) một cách tuần tự, ta thu được
a2p

,,.,.. = ________ 1_____________ (18)

( 8)

2 Pp!(v + 1)(V + 2)...(V + p)
2

Như vậy, nghiệm của phương trình Bessel (1.1) được xác định với độ chính xác
theo thừa số tùy ý ao- Ta có thể cho ao ở dạng

a

0=

(1 9)

-

2I-;

với r là hàm Gamma-Euler. Theo tính chất của hàm Gamma-Euler [? ]
r(v + I)(v + I)(v + 2)...(V + p) = r(p +1 + V)

Từ công thức (1.8) và (1.9), ta thu được
ao. =
(-1)p
2
2 p + r(v + 1)r(p +1 + V).

2p

Xét chuỗi
XX
í |7
/ x \2k ■ 2
7 (x) = V —________( 1)___________- x .
J (x)
2-/r(k + 1)r(k + V + 1) 2
v

k=0

________

.

(1.10)
()

Dùng quy tắc d’Alembert, có thể chứng minh được chuỗi (1.10) hội tụ tuyệt đối
với mọi X.
Định nghĩa: Chuỗi (1.10) được gọi là hàm Bessel loại một bậc V và được ký
hiệu là J (x).
Ta thấy rằng hàm J (x) là một nghiệm riêng của phương trình Bessel (1.1).
Xét trường hợp khi ơ = — V. Đặt Yv(x) = J- (x), thực hiện một cách tương tự,
ta cũng sẽ đi đến định nghĩa sau
Định nghĩa: Chuỗi (1.10), ứng với ơ = -V
v

v

v

1

(_1)k

Y, (x) = V —_______( —________- x
v(x)

k=0

2-/r(k + 1)r(k - V + 1) 2

.\2k-v

/rr

(1.11)(.)


Hình 1.1. Hàm Bessel loại một J (x) ùng với n = 0; 1; 2.
n

là hàm Bessel loại hai bậc V.

Hình 1.2. Hàm Bessel loại hai Y (x) ùng với n = 0; 1; 2.
n

Khi V nhận giá trị không nguyên (V = ±1; ±2..) thì hàm Y/(x) là nghiệm thứ hai
của phương trình Bessel. Nghiệm này sẽ độc lập tuyến tính với hàm Bessel loại
một J (x) nên các hàm J (x) và Y (x) sẽ hình thành một hệ nghiệm cơ bản của
phương trình Bessel bậc V.
v

v

v

Nếu V = n - số nguyên, thì
Yn(x) = (-1)n Jn(x).

Khi đó các hàm Yn(x) và Jn(x) sẽ phụ thuộc tuyến tính và khơng hình thành hệ
nghiệm hàm cơ bản.
Ta chứng minh biểu thức (1.12). Vì r(—k) = 1, k = 0; 1,..., nên tổng chuỗi trong
công thức (1.11) bắt đầu từ giá trị k = n, ta có
/ )
-1
( —1)k (
x 2k-n
x
(
)
=
V
—__________
—)
:______Y x
(x)
z^r(k - n + 1)r(k + 1) 2
n

k=n

V
(-1)n+s
= 2^r(s + 1)r(s + n + 1)

/

s=0

x )2s+n
2

= (-1)n Jn(x).

Các hàm trụ Bessel với các chỉ số dưới liên tiếp nhau
đạo hàm liên hệ với nhau bằng hệ thức truy hồi

(V — 1, V

và V +1) cùng các
(1.13)

J^+1(x) = -Jv-1 (x) + ^jJ (x),
J'v (x) = -J^+1(x) + J (x),
Jv(x) = Jv—1 (x) - V.ỵ(x).

(1.14)
(1.15)

Chú ý rằng, các hàm trụ Bessel với chỉ số dưới là số bán nguyên được biểu diễn
thông qua các hàm sơ cấp. Thật vậy, ta có

và

1 • 3 • • • (2k + 1) p2k+1
vY

r /2 + k)

Từ phương trình (1.10), ta có
J (x)
2

(—1)k
= 12 k!r(3 + k)
k=0

,/2ĩ V dCnb x
V ~x (2k + 1)!
k=0

/

x )2k+2
2

í2 ■ x.
V ~x


Chú
Tương
ý rằng,
tự, tatừ
cũng
cơng

cóthức (1.13) và các biểu thức của các hàm số J1 (x) và Y1 (x) ta
có cơng thức tổng quát sau
2
Y1 (x) =

- - -cos x.
xx

2

J

n+1

(x)

=

với P (v) và Q (v) là những đa thức có bậc khơng vượt quá n phụ thuộc vào v,
n

n

ngoài ra Pn(0) = 1; Qn(0) = 0.
Khi x ! 1 thì ta có công thức tiệm cận của hàm Bessel
4 ) + O (x

Jv (x) =

(1.16)

Điểm không của hàm Bessel: Là những điểm mà tại đó hàm Bessel nhận giá
trị bằng 0 [7]. ở đây ta chỉ xét điểm không của hàm Bessel loại một có chỉ số dưới
là số nguyên. Cụ thể, xét phương trình J (x)=0. Phương trình này ln có một bộ
nghiệm dương và các nghiệm này được phân bố theo thứ tự tăng dần, tức là
n

T’
Thực hiện tính tốn các giá trị của 6 điểm khơng đầu tiên của hàm J (x) với độ
chính xác đến 4 chữ số thập phân, ta được:
0

//1' = 2.4048; 4° = 5.5201; 4° = 8.6537;
gỊ0 = 11.7015; $0 = 14.9309; 4° = 18.0711.

Chú ý rằng, ta có thể tìm khơng điểm của hàm Bessel bằng cách sử dụng công
thức tiệm cận (1.16), cụ thể cho J (x) = 0, ta suy ra
v

(ỳ)

3X

Xì'

Z

hk = 4 + 2 + kX; k 2 .

1.1.2.

Các h àm trụ khác

Cùng với hàm Bessel J (x) thì trong các bài tốn vật lý thường sử dụng các
hàm trụ khác [8]. Ta xét một số hàm trụ sau
v

Hàm Hankel loại một
//,,(1' (x) = ~~~—[Jv(x)e~ '
sin Xì

i KU

— Yv(x)];

ì = U; n 2 Z;


_ (__1)” @YÍ (x)
dv ____) @V

Hàm Hankel loại hai

”

(

V

= n;

—.------[Jv(x)e ' — Yv(x)]; V = n, n 2 Z;
H
(Xi) sin
= J”nv
(x) + ỴSJ'(x)
Một bộ
đôi hàm số bất kỳ dJvtừ(x) bộ các hàm số J (x); N (x); HÍ1)(x); HÍ2)(x) tạo
Hàm
Neumann
N”Jn(x)
(x) = 1- ỉ \ dJv (x)
= n;của V. Từ đó ta
thành hệ nghiệm cơ bản của phương
trình Bessel với mọi giáV trị
n
dv
Nv (x)2) = —-----[Jv(x)cosnV — Yv(x)]; V = n, n 2 Z;
H
(X) sin nV

”

i KÍ

1)

Í

Í

Hàm Infeld và hàm Macdonald
Tí\

(

n

AT/- \

ỉv(x) = exp í — ỶVÍ} Jv(ix); Kv(x)

ni
= 2 exp

T^/\

2 vộ HÍ1)(ix).

Hàm Hankel loại một và loại hai được biểu diễn qua hàm Bessel và Neumann
như sau
HÍ\x) = J (x) + ÍNV(x); HÍ2)(x) = J (x) - ÍNV(x).
Í

Í

(1.17)

Sử dụng cơng thức tiệm cận của các hàm Bessel J (x), ta có các biểu thức tiệm

cận của các hàm trụ như sau
v

n

HÍ1) (x) =

nV-----7
2
n

HÍ2)(x) =

N

Í

(x)

— V----7
2
n
n\ .
_3\
+Ox2;
\
)
— V 47
.x
--------O(x 1

7 2nx .= (1 + ~ ));
2
K (x) = 4 ỉỉe~x (1 + O(x-1)) .
2x
Í

(1.18)


có cơng thức nghiệm tong qt của phương trình Bessel là
u(x) = C1 Jv (x) + C2N., (x)

hoặc
u(x) = CiHÍ1)(x) + C2H2 (x).

Nếu trong phương trình Bessel, ta thay x bởi ix, thì các hàm Infeld và Macdonald
sẽ tạo thành hệ nghiệm cơ bản của phương trình, từ đó ta có nghiệm tổng quát

u(x) = C1IV (x) + C2K., (x).

Khi x ! 0, ta có biểu thức mơ tả trạng thái của các hàm trụ
•Ạ'(x) - wn d
1(V + 1) \x

I

< 2 inx V = 0;
N(x)
- < 7rr(V) /x \-^
I.I;V

l % \2>
’
±i2 In

HÍ1;2)(x)
±l

2;V=0;
2
^
r(v) /x Y~v

~

82
K (x) « < -. x’

> 0;

’ V > 0;

,x’
> 0.

1.2. Đa thức Legendre
1.2.1. Đa thức trực giao cổ điển Legendre
Trong các bài toán vật lý toán và vật lý lý thuyết, đa thức Legendre là một hệ
đa thức hoàn chỉnh và trực giao. Đa thức Legendre có thể được định nghĩa theo
nhiều cách khác nhau và mỗi cách định nghĩa làm nổi bật các tính chất cũng như
những ứng dụng khác nhau trong việc giải các bài tốn biên.

Ta sẽ tìm hiểu định nghĩa và các tính chất của đa thức Legendre thơng qua



phương trình vi phân của bài tốn Sturm-Liouville
Phương trình (1.19)d được gọi
là phương trình vi phân Legendre. Ta sẽ tìm nghiệm
(1 - x ) ì
Ay điểm
= 0; —1
1;
(1.19)
của phương trình (1.19dx
) trong lân cậndxcủa
x =<0xở< dạng
chuỗi lũy thừa
2

(1.20)
(1.21)

|y(-l)| < 1; |y(1)I < 1.

y(x) = X ak xk.
k=0

Thay biểu thức (1.21) vào phương trình (1.19) và thực hiện một số phép biến
đổi cơ bản, ta thu được
X[(k + 2)(k + 1)ak+2 + (A - k2 - k)ak]xk = 0;
k=0

suy ra
[(k + 2)(k + 1)ak+2 + (A — k2 — k)ak]xk = 0; k = 0; 1.

Từ đó, ta có cơng thức truy hồi
A - k2 - k
ak+2
= _ (k + 2)(k + 1)ak.

(1.22)

Công thức này cho phép ta biểu diễn các hệ số chẵn qua hệ số a và các hệ số
0

lẻ qua a .
1

Khi a = 0; a = 0, ta có nghiệm riêng chỉ chứa lũy thừa bậc chẵn của x
0

1

yi(x) = X a2px2 .
p

(1.23)

p=0

Khi a = 0; a = 0 - nghiệm riêng chỉ chứa lũy thừa bậc lẻ của x

0

1

y2(x) = X a2p+ix2 +1.
p

(1.24)

p=0

Ta thấy rằng, các chuỗi y (x) và y (x) hội tụ trên đoạn [—1; 1]. Nếu A = n(n + 1),
1

thì j

y(-1)j < 1; jy(1)| < 1.

2


Thật vậy, từ công thức truy hồi (1.22), khi A = n(n + 1), ta có
a

n+2 — an+4 — ■■■

— an+2p — ■■■ — 0;

nghĩa là một trong các chuỗi (1.23) hoặc (1.24) sẽ triệt tiêu và sẽ tạo thành đa
thức bậc n.

y(x) — Pn(x), n — 0,1.

Đa thức trên là hàm riêng của bài toán (1.19), (1.20). Các đa thức này được gọi
là đa thức Legendre.
Xét một số tính chất của đa thức Legendre [11]
1. Các đa thức Legendre trực giao trên đoạn [—1; 1] với trọng số p(x) — 1, tức là
Ị Pn(x)Pm(x)dx

— 0, n — 'm

Thật vậy, theo phương trình Legendre, ta có các đồng nhất thức sau:
Nhân phương trình thứ nhất với
P d(x), và phương trình thứ hai với P (x), trừ vế
1
m

d
í - x2) ~Ề
+ n(n1]+, 1)P
n(x) = 0,
dxphân trên đoạn [—1;
theo vế và sau đó lấy tích
ta được

Ạ í1 - x2) d~Ề
- P„d
dPm
(1
- x2)
dx

dPm
m (1, -X2)
n,
+ m(m + 1)Pm(x) = 0^ dx
dx
dx
dx
Pm

— [m(m +1) — n(n + 1)] Ị Pn(x)Pm(x)dx

hay

í1 d ((1 - x2) [dPnPm - dPmPnlỊ dx
J_1 dx
[ dx dx J J
— (m — n)(m + n + 1) Ị Pn(x)Pm(x)dx■
Ị Pn(x)Pm(x)dx

Từ đó

n

dx


a
(m — n)(m + m + 1)

dPm

\
, Pny
dx

1

= 0;
-1

khi n = m.
2. Đa thức Legendre hoặc là đa thức chẵn hoặc là đa thức lẻ.
3. Đa thức Legendre thỏa mân công thức Rodrigues [12]
P

n(x)

dn
= CdXn(x2 - 1)n:

(1.25)

Thật vậy, hàm u(x) = C(x — 1) là nghiệm của phương trình
2

n

(x2 — 1)u0 — 2nxu = 0;

Lấy vi phân phương trình trên n +1 lần, áp dụng công thức Leibnits đối với đạo
hàm của tích hai hàm số, ta được

. d.u(n') ’
(x2 - 1) 'Ị
— n(n + 1)u(n) = 0;
dx

_d
_

dx

từ đây suy ra hàm u (x) là nghiệm của phương trình Legendre khi q = n(n + 1),
và u (x) là đa thức bậc n, trùng với đa thức Legendre.
4. Đa thức Ledendre P (x) có n khơng điểm khác nhau trong đoạn [—1; 1].
5. Trong biểu thức (1.25) đặt C = (n) ; ta được đa thức Legendre chuẩn hóa
(n)

(n)

n

Pn(x) = -7-^(x2 - 1)n
n
'
2(n)n! dxn-

(1.26)
v7

7

Sử dụng cơng thức (1.26), ta có thể tính được các biểu thức một số đa thức Legendre
đầu tiên

Po(x) = 1; P1(x) = x P2(x) = - 1 + 3 x2; P 3 (x) = 5 x3 - 2 x

P ( x)

4

(x)

= 3 + 35x4_ 15x2;
8+8 x 4 x ; P5(x)
P6(x) =
P7(x) =

(x)

16

16

= 63x5
8x

-Ẵ + 231 x6 - 3^x105
4 + -x
,22;
16
16

’

429 7
693
16x _ 16x +16x _ 16x;

315

5

3003
P s(x) = —— ++ 6435
128 126x _ ÍTx + 6Tx _ 62x ;
8

8V 7

35

3

3465

6

P _ 12155
6435 7
9009
9(x)
= 128 x _ 32 x + 64 x _ 32 x + 128x;

9

35 3 . 15
- 4 x + 8 x;

5

315

4

1155

3

2

315



-

63
1o(x)

46189
109395
45045
15015

= _ 256 + 256 x 256 x + 128 x _ 128 x + “256’x ■
10

8

6

4

3465

2

Sau này khi nói về đa thức Legendre chuẩn hóa thì ta sẽ dùng kí hiệu P (x),
thay cho P (x). Hình (1.3) biểu diễn 6 hàm đa thức Legendre chuẩn hóa đầu tiên.
n

n

Hình 1.3. Đồ thị đa thức Legendre P (x) ứng với n = 0; 1; 2; 3; 4; 5.
n

1.2.2.

Đa thức Legendre liên hợp

Trong toán học, đa thức Legendre liên hợp là nghiệm chính tắc của phương
trình Legendre tổng qt
(1
(1

_

x

2
p4Ĩ. p) dx2 Pn (x)

_ 2x—p
x

2x

dx

P (x)

n

+
+

ín,(n, +1)-

''"2 ,

n(n + 1)

x2

P,ỉm) (x)=0.

(1.27)

hoặc ta có thể viết lại ở dạng tương đương
Trong đó các hệ số n và m là các số nguyên. Phương
trình (1.27) (hay (1.28))
m2
d (1 — x ) — P( (x) + n( n + 1)
P,
í
= 0.
có các nghiệm khác khơng chỉ trên đoạn [—1; 1] khi n và m(x)nguyên
với 0 < m
< n.
(1.28)
1
—
x
x
x + n n +1
dx
(±
)
d'n
(
)
(

)
Khi m là số chẵn thì hàm Pn (x) là một đa thức, khi m = 0 và n ngun thì hàm
số chính là đa thức trực giao cổ điển Legendre (đã được khảo sát ở mục 1.2.1).
Phương trình vi phân Legendre thường được gặp trong các bài toán vật lý toán
và các bài tốn vật lý lý thuyết. Cụ thể, nó xuất hiện trong việc giải phương trình
Laplace (và phương trình vi phân đạo hàm riêng) trong hệ tọa độ cầu. Các đa thức
2

m)

m)

2

m)


Legendre liên hợp đóng vai trị quan trọng trong việc định nghĩa hàm cầu.

Định nghĩa các tham số nguyên không âm n và m
Đa thức Legendre liên hợp được biểu diễn qua đa thức Legendre cổ điển như
sau
(1.29)

FÍm)(x) = (—1)m(1 - x2)m-dm(Pn(x)).
n
\/
\
/\
/ J„m " n' 77

'

7

Các hàm số này trong biểu thức (1.29) thỏa mân phương trình vi phân Legendre
với các giá trị xác định của các tham số n và m bằng cách lấy vi phân m lần phương
trình Legendre đối với P (x) [13]
n

(1 _ x2)^uPn2 (x) - 2x^-Pn(x) + n(n + 1)Pn(x) = 0.
dx
dx

Mặt khác, theo công thức Rodrigues
1 dn
õpn. dxn

Hàm Pn (x) có thể được viếtPnlại
(x) =
m)

m

Pnm} (x) = ...
nv
'
2"n'

[(x2 - 1)n] .

n+m

(1 - x2)* Ặ

dxn+m v

(x2 - 1)n

(1.30)
(1.31)
'

Phương trình này cho phép mở rộng giới hạn của m thành: — n < m < n.
Ta có biểu thức liên hệ giữa Pn~ ~) (x) và Pn™) (x)
m

/> m (x) = (_i)m(^_________HIP(m)(x).
Pn
(x)
( 1}
(n + m)!Pn (x).

(1.32)

Tính trực giao
Giả sử 0 < m
giá trị m cố định

< n,

đa thức Legendre liên hợp thỏa mân điều kiện trực giao với

P^Pnr^dx

2(n + m)! x
(2n + 1)(n + m)! ỏfc’n’

trong đó ỏk, là ký hiệu Delta Kronecker.
Ngồi ra, đa thức Legendre liên hợp còn thỏa mân điều kiện trực giao khi cố
n

(1.33)


định n, tức là
0 với m = P;
(n + m)!
—-———-7 với m = p = 0;
m(m — m)!

P ( )P (p)
m

1 — X2

dx

1 với m = p = 0.

Một số hàm đa thức Legendre liên hợp

Hình 1.4. Da thức Legendre liên hợp với m =1 (bên trái) và m = 2 (bên phải).
Một số hàm đa thức Legendre liên hợp đầu tiên với các giá trị nguyên âm và
nguyên dương của m
P(0)
P

0

T3>(—2) Ị
P
(X)
2
p_(1)
P (X)
2

1 P(-1)
= 1; P1

(x)

= -3

\

1

D

(2) ( \

= 24P2

X(1

-

X

)

2

_____1

(x)

=2

(X)

; P (X)

2

;

22

=

3(1

P(1)

P

D(—1)/™\
P

2

-

(X)

1

(x);

p(o^ M = V P(1)D = —lơ — x2 2
P
(x)
P
X
1
= X;
1 ( )
= (1 x ) ;

1

=6

p(2)
3

D (1)D(0)
P
X
2 ( );
2\

X ); P

P

2

p(-3)

1
ÍQ ...2
= 2(,yx

___

(X)

(X)

=

;

1)

____1 p(3)
720P3 (X);

P„(~2)(X) = 1 P„(2) (X); P„(-1)(X) = zỉ P„(1)(X). P„(0)(X) = Ã^3 — 3X);
P
(X)
(X)
3
= 120P3 (X); P3
= 12 P3 (X); P3 (X) = 2(OX 3X);
p(1)
P3

______3

2

2

(X) — 2 (5X 1)(1 X ) 2 ;

2 p(2)

2
\ p(3)
P3 (X) — 15X (1 X );

P3 (X) —

2
2
15(1 X ) 2 .

Đồ thị của các hàm đa thức Legendre liên hợp với các giá trị khác nhau của m
và n được biểu diễn trên hình (1.4).


1.3. Hàm cầu
1.3.1. Định nghĩa và các tính chất cơ bản của hàm cầu
Hàm cầu là phần góc của họ các hàm trực giao trong phương trình Laplace, được
biểu diễn trong hệ tọa độ cầu. Hàm cầu được sử dụng rộng rãi trong việc nghiên
cứu các hiện tượng vật lý trong các miền không gian được giới hạn bởi các bề mặt
cầu và trong lời giải của các bài toán vật lý có tính đối xứng cầu. Hàm cầu đóng
vai trị quan trọng trong lý thuyết phương trình vi phân đạo hàm riêng và trong
vật lý lý thuyết, đặc biệt trong các bài tốn tính tốn quỹ đạo của electron trong
nguyên tử, trọng trường của Geoid [14], từ trường của các hành tinh và cường độ
bức xạ di tích.
Phương trình Laplace có dạng
Au = 0
d2u

d2u
dy-

, dx2

d 2u
dz2

’

Trong hệ tọa độ cầu phương trình trên được viết lại như sau
1 d í 2 du \
1dí.
r2 dr y dr ) + r2 sin ỡ dỡ ự

du \
1 d2u
2
2
dỡ J + r sin ỡ d'2

Nghiệm của phương trình này được tìm bằng phương pháp tách biến
u(r,ỡ,') = R(r)Y (ỡ,')_

Từ đó ta thu được hai phương trình vi phân, phương trình của hàm bán kính
1 d 2 dR
Rdr dr

và phương trình hàm góc
A 'Y + XY = 0;
Ớ

với Aỹ'— phần góc của tốn tử Laplace trong hệ tọa độ cầu, có dạng
1d
3ífi

d\
1d
2
sin ỡ dỡ ỵ dỡ J + sin ỡ d'2 ’
sm

(1.34)


Phương trình (1.34) là bài tốn Sturm-Liouville trên mặt cầu đơn vị: 0 < ỡ <
%. 0 < ' < 2%. với điều kiện tuần hoàn
Y (ớ.') = Y (ớ.' + 2%).

(1.35)

và điều kiện hữu hạn
|Y( ,

0 ')\ < 1. |-.') < 1.

(1.36)

Định nghĩa: Nghiệm hữu hạn của phương trình (1.34) trên mặt cầu đơn vị
thỏa mân điều kiện tuần hồn theo ' và có đạo hàm liên tục đến bậc hai được gọi
là hàm cầu. Nghiệm của bài toán (1.34) - (1.36) được tìm bằng phương pháp tách
biến ở dạng

Y (ớ.') = e(ớ)T(').

(1.37)

Thay (1.37) vào phương trình (1.34), ta được
T" + VT = 0. $(') = T(' + 2%).

(1.38)

Khi V = m nghiệm khơng tầm thường của phương trình này có dạng
2

M') = eim'. m = 0. ±1....

Đối với hàm 0(ớ), ta có phương trình
-7^7Ẫ

(sinớ@|ì +
sin ớ de de

(

À

(1.39)

--iì ỡ = 0.
sin ớ
2

0 < ớ < %. |0(ớ)| < 1. |0(%)| < 1.

Nếu thay đổi biến số x = cosớ. y(x) = y(cosớ). thì phương trình (1.39) có dạng
-d [(1 - x2) dy 1 +
dx p 7 dxj ỵ

(

À

ì y =0.
1 — x2 y

(1.40)

— 1 < x < 1. |y(±1)| < 1.

Bài toán (1.40) là bài toán Sturm-Liouville đối với hàm Legendre liên hợp, vì vậy
trị riêng có dạng
Àn = n(n + 1).