Hàm riêng của toán tử sturm liouville trên khoảng hữu hạn và trên khoảng vô hạn

ĐẠI HỌC QUỐC GIA HÀ NỘI TRƯỜNG ĐẠI HỌC KHOA HỌC TỰ NHIÊN ——————- Nguyễn Viết Đại HÀM RIÊNG CỦA TOÁN TỬ STURM-LIOUVILLE TRÊN KHOẢNG HỮU HẠN VÀ TRÊN KHOẢNG VÔ HẠN LUẬN VĂN THẠC SĨ Hà Nội - 2019 ĐẠI HỌC QUỐC GIA HÀ NỘI TRƯỜNG ĐẠI HỌC KHOA HỌC TỰ NHIÊN ——————- Nguyễn Viết Đại HÀM RIÊNG CỦA TOÁN TỬ STURM-LIOUVILLE TRÊN KHOẢNG HỮU HẠN VÀ TRÊN KHOẢNG VƠ HẠN Chun ngành: Tốn giải tích Mã số: 8460101.02 LUẬN VĂN THẠC SĨ KHOA HỌC Cán hướng dẫn: TS Đặng Anh Tuấn Hà Nội - 2019 LỜI CẢM ƠN Trước trình bày nội dung khóa luận, em xin bày tỏ lịng biết ơn tới thầy Đặng Anh Tuấn Thầy tận tình hướng dẫn để em hồn thành luận văn Thầy không hướng dẫn em mặt chuyên môn tốn, thầy cịn dạy em nhiều điều sống Những lời dạy bảo thầy giúp em nhìn chuyện, giúp em vượt khúc mắc, yếu đuối mặt tâm lý mà tưởng chừng vượt qua Em xin lỗi thầy nhiều yếu đuối muốn bỏ cuộc, em ngắt liên lạc với thầy, thầy không bao dung ln quan tâm đến em em tiếp tục Em xin bày tỏ lịng biết ơn chân thành tới tồn thể thầy giáo khoa Tốn - Cơ - Tin học, Đại học Khoa Học Tự Nhiên, Đại Học Quốc Gia Hà Nội dạy bảo em tận tình suốt trình học tập khoa Nhân dịp em xin cám ơn tới ông nội , ông bà ngoại, bố mẹ cậu mợ em người thương yêu , quan tâm che chở cho em Ngoài em xin cảm ơn trung tâm anh ngữ ViViAn tận tình dạy cho em để thi tiếng anh B1 Em xin cảm ơn anh Đỗ Duy Hiếu nhận em vào làm trung tâm anh để có tiền trang trải sống suốt thời gian Hà Nội, em xin lỗi bỏ mà khơng nói lời Em xin cảm ơn Viện Tốn kí hợp đồng với em tháng, khơng có hợp đồng làm động lực để quay lại em khơng thể vượt qua tiếng anh B1 Cuối em xin cảm ơn bạn Tô Thị Vân Anh bạn Nguyễn Đức Ngà, bạn Vân Anh liên lạc gọi em lại học tiếng anh chuyên ngành, bạn Ngà hướng dẫn em bước làm thủ tục bảo vệ Hà Nội, ngày tháng 12 năm 2018 Học viên Nguyễn Viết Đại LỜI MỞ ĐẦU Từ đại số tuyến tính hữu hạn chiều, cho ma trận đối xứng ta tìm thấy sở trực chuẩn khơng gian gồm tồn vectơ riêng ma trận Khi ta có khai triển n v = å(v, vk)vk k=1 với vk vectơ riêng chuẩn hóa ma trận A Ngồi ta có đẳng n thức Pythagoras jvj = å k=1 j(v, vk)j Từ lý thuyết chuỗi Fourier hàm f tuần hồn chu kì 2p khả vi liên tục R có khai triển +¥ f (x) = vk(x) = e ikx p å ( f , vk)vk(x) k= ¥ / 2p hàm riêng chuẩn hóa ứng với giá trị riêng k toán tử vi phân thường +¥ d2 Ngồi ta có đẳng thức Parseval dx2 jj f jj = å ¥ j( f , vk)j Chuỗi Fourier xuất ta giải phương trình truyền nhiệt, dao động sợi dây, dao động màng mỏng, phương pháp tách biến Sự tương tự vấn đề đại số tuyến tính lý thuyết phương trình nhà toán học thấy từ lâu trước Tuy nhiên D.Hilbert người hệ thống lại tương tự việc làm lý thuyết phương trình tích phân, xem [5] Một kết việc làm làm nảy sinh khơng gian Hilbert l sau khơng gian Hilbert tổng qt Xây dựng tốn học cho khơng gian l không gian Hilbert trừu tượng dẫn đường cho phát triển mạnh mẽ lý thuyết phổ tốn tử tự liên hợp khơng gian Hilbert Lý thuyết phổ trừu tượng hoàn thiện, định lý sở toàn lý thuyết định lý khai triển phổ Một tốn tử tự liên hợp khơng gian Hilbert khai triển thông qua phép chiếu phổ El (còn gọi họ phổ giải thức đơn vị) Tuy nhiên trường hợp tốn tử cụ thể thơng tin tiệm cận giá trị riêng, hàm riêng họ phổ Trong luận văn em đọc hiểu trình bày chi tiết lại kết khai triển hàm riêng toán tử Sturm-Liouville cho hai trường hợp khoảng hữu hạn nửa đường thẳng Nội dung luận văn gồm chương Chương 1: kiến thức chuẩn bị Chương 2: khai triển khoảng hữu hạn Chương 3: khai triển nửa đường thẳng Nội dung chương trình bày cơng thức tiệm cận giá trị riêng hàm riêng toán tử Sturm-Liouville, chứng minh tồn dãy đếm giá trị riêng cách khác nhau: sử dụng định lý Rouche, lý thuyết dao động Sturm, phương pháp phương trình tích phân Ngồi chương có cách chứng minh khác cho định lý khai triển hàm riêng : phương pháp phương trình tích phân, phương pháp thặng dư Cauchy Ở cuối chương định lý , hội tụ điểm khai triển hàm riêng Sturm-Liouville giống hội tụ điểm chuỗi Fourier thông thường Nội dung chương 3, xây dựng hàm phổ r(l) (cịn gọi độ đo phổ) từ định nghĩa biến đổi Fourier tổng quát thu đẳng thức Parseval định lý khai triển dạng tương tự chương Đồng thời chương trình bày phân loại giới hạn điểm, giới hạn trịn tốn tử Sturm-Liouville nhiên em chưa tìm hiểu xuất phát điểm vật lý khái niệm Ngoài chương trình bày biểu diễn tích phân giải thức, rõ họ phổ E l toán tử Sturm-Liouville Ở cuối chương ánh xạ f (x) 7!F(l) đặt tương ứng hàm f (x) L (0, ¥) với biến đổi Fourier tổng quát F(l) 2 L ( ¥, +¥) ánh xạ Unitary ( song ánh bảo toàn chuẩn) r(l) Các kết mục 2.2 tham khảo [7] [9], mục 2.3 tham khảo [4] [11], mục 2.4 2.5 tham khảo [4] [9], mục 2.6 tham khảo [8] [9], chương tham khảo [9], họ phổ E l trình bày trừu tượng tìm đọc [6] phụ lục [9] Hà Nội, ngày tháng 12 năm 2018 Học viên Nguyễn Viết Đại Mục lục Lời mở đầu Kiến thức chuẩn bị 1.1 Tính trù mật 1.2 Một số định lý phương trình vi ph 1.3 Một số định lý giải tích phức 1.4 Một số kết tích phân Khai triển khoảng hữu hạn 2.1 Giới thiệu số tính chất 2.2 Công thức tiệm cận cho giá trị riê 2.3 Phân bố không điểm hàm riê 2.4 Hàm Green, toán tử compact đối xứn 2.5 Định lý khai triển đẳng thức Parse 2.6 Chứng minh định lý khai triển tí 2.7 Hội tụ điểm khai triển hàm riêng Khai triển nửa đường thẳng 3.1 Đẳng thức Parseval với nửa đường t 3.2 Giới hạn điểm, giới hạn tròn 3.3 Biểu diễn tích phân giải thức 3.4 Tính trực giao khai triển Tài liệu tham khảo 93 Chương Kiến thức chuẩn bị 1.1 Tính trù mật Ký hiệu C[a, b] khơng gian hàm giá trị phức, liên tục khoảng mở hữu hạn (a, b), liên tục phải a liên tục trái b Khơng gian C[a, b] có tích vơ hướng cho bởi: Zb ( f , g) = a f (x)g (x)dx, f , g C[a, b], g (x) liên hợp phức g(x) Cho D(L) tập C[a, b] xác định D(L) = fy(x) C [a, b] : BCa(y) = BCb(y) = 0g, BCa(y) = y(a)cos(a) + y0 (a)sin(a) = 0, BCb(y) = y(b)cos(b) + y0 (b)sin(b) = (a, b R), C [a, b] không gian hàm giá trị phức, khả vi liên tục cấp hai (a, b), khả vi liên tục cấp hai bên phải a bên trái b Khi ta có khẳng định sau: Bổ đề 1.1.1 ([2]) D(L) trù mật khơng gian C[a, b] với chuẩn cảm sinh từ tích vô hướng 1.2 Một số định lý phương trình vi phân thường 1.2 Một số định lý phương trình vi phân thường Bổ đề 1.2.1 (Cơng thức Liouville) ([4]) Xét phương trình y00 (x) + p(x)y0 (x) + q(x)y(x) = 0, với p(x), q(x) C[a, b] Giả sử y1(x) y2(x) hai nghiệm phương trình Khi định thức Wfy1, y2 g(x) = y1(x)y20(x) y10(x)y2(x) Wronskian y1(x) y2(x) cho công thức Liouville: Zx Wfy1, y2g(x) = c.exp p(t)dt 8x [a, b], (1.2.1) a với c số Bổ đề 1.2.2 (Bất đẳng thức Gronwall-dạng vi phân) ([4]) Cho h(.) hàm không âm, liên tục [0, T] thỏa mãn bất đẳng thức vi phân h0 (t) f(t)h(t) + y(t), 8t [0, T], f(t), y(t) hàm không âm liên tục [0, T] Khi đó: Zt R t f(r)dr [h(0) + y(s)ds], 8t [0, T] h(t) e 0 Bổ đề 1.2.3 (bất đẳng thức Gronwall-dạng tích phân) ([4]) Cho x(t) hàm không âm, liên tục [0, T] thỏa mãn theo t bất đẳng thức tích phân: x(t) với số C1, C2 x(t) Z C1 t x(s)ds + C2, Khi đó: C2(1 + C1te Ct ), t T Định lý 1.2.1 (định lý tồn nghiệm ([9]) ) Nếu q(x) hàm liên tục [a, b], với a R, l C toán Cauchy: y00 (x) + (l q(x))y(x) = 0, (1.2.2) cos(a), (x0 [a, b] cố định ) có nghiệm j(x, l), x [a, b] Với x cố định thuộc [a, b] hàm j(x, l) hàm nguyên l, tức hàm chỉnh hình toàn mặt phẳng phức C j(x0, l) = sin(a), jx(x0, l) = 3.4 Tính trực giao khai triển Chứng minh Đầu tiên ta xét f (x) = fn(x) triệt tiêu bên ngồi đoạn [0, n] Ta có Zb Zn Gb(x, t, z) = Z fn(u) j(t, l)j(u, l)drb(l)du Z Zb Zn Gb(x, t, z)j(t, l)dt D dt D fn(u)j(u, l)du drb(l) l l lD k,b k,b + a = l lk,b l+D a k,b = j x, Z = D z l j( x Z z D Z b Gb(x, t, z) Z n fn(u) Z D j(t, l)j(u, l)drb(l)du dt = + Zb a Zn Gb(x, t, z) Z fn(u) j(t, l)j(u, l)drb(l)du dt D Sử dụng bổ đề 3.2.2, với a > cố định, cho b ! ¥ số hạng vế phải tiến tới Za Zn G(x, t, z) Z j(t, l)j(u, l)dr(l)du fn(u) dt D Ta chứng minh số hạng thứ hai tiến tới khơng a ! ¥, b > a Để làm Z a Gb(x, 3.4 Tính trực giao khai triển điều ta sử dụng bất đẳng thức Cauchy-Schwarz: Gb(x, t, z) b Z a b Z j a b Z j a Gb(x, t, z) Gb(x, t, z) b Z j a Gb(x, t, z) b Z jGb(x, t, z)j dt a Ta chứng minh Zb a tiến tới cho a, b đủ lớn cho ! 2 jGb(x, t, z)j dt ¥, b > a Thật vậy, cho trước e > bất kỳ, chọn a Z ¥ jy(t, z)j dt e a Từ bổ đề 3.2.2 ta chọn b đủ lớn cho Z Z )j dt )j dt e, Khi ta thu Z a b + jy(t, z)j dt 3e e 3.4 Tính trực giao khai triển Như cho b ! ¥ (3.4.3) ta chứng minh (3.4.2) Với trường hợp f (x) L (0, ¥) bất kỳ, ta đặt fn(x) = f (x) x [0, n] fn(x) = x > n Khi ta có với Fn(l) = Cho n ! ¥ trong(3.4.4) ta thu Từ hồn thành chứng minh Nhận xét: Cho g(x) L (0, ¥) Nhân hai vế (3.4.2) với g(x) sau lấy tích phân từ đến n ta Z0 Cho n ! ¥ ta Z Zn G(l) = l.i.mn!¥ g(x)j(x, l)dx Họ tốn tử ED có tính chất sau Mệnh đề 3.4.1 ([9]) Cho f , g L (0, ¥) Khi (i) tự liên hợp: (ED f , g) = ( f , EDg), (ii) đơn điệu: D0 D (ED0 f , f ) (ED f , f ), (iii) đầy đủ: (E( ¥,+¥) f , g) = ( f , g) = (iv) trực giao: EDED0 = ED D0 R ¥ f (x)g(x)dx, 3.4 Tính trực giao khai triển D D0 giao D0 với D Chứng minh Nhân hai vế (3.4.1) với g(x) sau lấy tích phân theo x từ đến n ta Zn Z¥ g(x) Z ED(x, t) f (t)dt dx = D Gn(l)F(l)dr(l) Cho n ! ¥ ta ( ED f , g ) = Z Từ (3.4.6) tính đối xứng hàm E D(x, t) ta thu (i) Lấy g = f (3.4.6) trở thành Z (ED f , f ) = F (l)dr(l) D Từ đẳng thức ta thu (ii) Từ (3.4.6) đẳng thức Parseval ta thu (iii) Đặt El = E( Thật vậy, D0 E E = E nE D D D =E 0 ( ¥,l +D ) 0 D ( ¥,l +D ) Từ (3.4.6) lấy D = ( ¥, l) ta Z (El f , g) = ¥ l F(l)G(l)dr(l) Từ tính chất (ii) ta thu hàm số l 7!(El f , g) hàm đơn điệu tăng theo l, dl(El f , g) = F(l)G(l)dr(l) Khi (3.4.5) viết lại thành d (E l (RzED f , g) = Z ( ¥,l) D f , g) z l (3.4.7) 3.4 Tính trực giao khai triển Mặt khác (3.3.12) thay f E D f ta Z0 Từ (3.4.7) (3.4.8) tính độ đo Stieltjes ta (E ( ¥,l) D f , g) (E ( ¥,l) E f , g) = c D Lấy l = +¥, từ tính đầy đủ ta c = Do f , g L (0, ¥) nên E ( ¥,l) D =E E ( ¥,l) D Bổ Z ED(x, u)ED0 (u, t)du = D D0 j(x, l)j(t, l)dr(l) = ED D0 (x, t) (3.4.9) Chứng minh Cho x cố định thuộc (0, ¥), với f (t) L (0, ¥) ta có: Z¥ Z¥ Z¥ f (t)dt = ED(x, u)ED0 (u, t)du Z¥ ED(x, u) Do E ED0 (u, t) f (t)dt du (x, t) DD thức xảy Z ¥ ED(x, u)ED0 (u, t)du = ED D0 (x, t) Do x lấy nên ta điều phải chứng minh Bổ đề 3.4.6 ([9]) Cho ED(x) = E ( ) x Khi với hầu khắp nơi D FD(l) = l 2/ D Chứng minh Xét trường hợp sin(a) 6= Lấy x = t = (3.4.9) ta Z¥ Z E (t)E D D0 (t)dt = ( ) D D0 r l (3.4.10) 3.4 Tính trực giao khai triển Với trường hợp sin(a) = ta thu đẳng thức cách đạo hàm (3.4.9) theo x, t lấy x = t = Chia hai vế (3.4.10) cho r(l Cho D0 ! từ bổ đề 3.4.4 vế trái đẳng thức tiến tới F D(l) hầu khắp nơi, vế phải tiến tới l D, tiến tới l 6= D Bổ đề 3.4.7 ([9]) Với x cố định, z khơng thực biến đổi Fourier hàm Green G(x,t,z) z Z ¥ j(x, l) l ¥ G(x, t, z)ED(t)dt = z l0 Z ¥ j(x, l) FD(l)dr(l) Chia hai vế cho r(D) ta Cho D tiến 0, sử dụng bổ đề 3.4.4 vế trái đẳng thức tiến tới biến đổi Fourier hàm Green G(x,t,z), vế phải tiến tới j(x, l) , từ ta z l điều phải chứng minh Bây ta quay lại chứng minh định lý mục Định lý 3.4.2 ([9]) Cho r(l) xây dựng định lý 3.1.1 Khi với hàm F(l) L Z f (x) = l.i.mN!¥ N N F(l)j(x, l)dr(l) thỏa mãn đẳng thức Parseval: Z¥ Z +¥ ¥ F2(l)dr(l) = f 2(x)dx 3.4 Tính trực giao khai triển Chứng minh Đầu tiên ta xem xét FN (l) hàm liên tục, có giá nằm [ N, N] Ta đặt ZN gN (x) = N FN (l)j(x, l)dr(l) Ta chứng minh Z ¥ g N (x)dx ¥ (3.4.11) Cho l0 = N, l1, l2, , ln = N phân hoạch khoảng ( N, N), cho xi [li 1, li] Ta đặt n g N (x) = Z n å FN (xi) l i=1 i li j(x, l)dr(l) Nếu max(li đoạn hữu hạn Từ bổ đề 3.4.5 ta có ¥ Z f n g N (x) Cho a > số cố định Từ (3.4.12) tồn số K cho với n ta có Cho n ! ta thu minh (3.4.11) Tiếp theo ta chứng minh biến đổi Fourier hàm g N (x) trùng với FN (l) hầu khắp nơi (theo độ đo r(l)) Từ (3.4.10) ta Z (3.4.13) ta tiến qua giới hạn dấu tích phân (3.4.14) để thu Z¥ Z gN (x)ED(x)dx = D FN (l)dr(l) (3.4.15) Chia hai vế cho r(D), sau cho D ! từ bổ đề 3.4.4 vế trái tiến tới biến đổi Fourier gN (x) hầu khắp nơi (theo độ đo r(l) ) vế phải tiến tới FN (l) Sử dụng đẳng thức Parseval ta Z ¥ g N (x)dx = Z 3.4 Tính trực giao khai triển Bây giờ, cho F(l) L liên tục FN (l) có giá nằm F(l) Cho ¥, +¥) bất kỳ, ta tìm dãy hàm [ N, N] hội tụ bình phương trung bình tới ZN gN (x) = Từ (3.4.16) với N0 > N ta có Z ¥ gN0 (x) gN (x) dx = Z Do gN (x) phải hội tụ bình phương trung bình tới hàm f (x) ta có Parseval Z¥ f (x)dx = Z +¥ F (l)dr(l) ¥ Cuối ta cần chứng minh hội tụ bình phương trung bình tới f (x) Cho H(l) Lr(l)( f (x) xây dựng từ F(l) Khi ta có Z ¥ f f (x) h(x)g2 dx = Cho H(l) = F(l) với Z ¥ ff (x) fN (x)g dx = Z N Z¥ F (l)dr(l) + ¥ F (l)dr(l) N Cho N ! ¥ ta fN (x) hội tụ bình phương trung bình tới f (x) 3.4 Tính trực giao khai triển KẾT LUẬN Khóa luận khơng đưa kết mới, nội dung khóa luận gồm đọc hiểu trình bày chi tiết lại kết khai triển hàm riêng toán tử Sturm-Liouville đoạn hữu hạn nửa đường thẳng Khóa luận trình bày với chương Ở chương kiến thức chuẩn bị phục vụ cho chứng minh chương chương Ở chương trình bày cơng thức tiệm cận giá trị riêng hàm riêng toán tử Sturm-Liouville, chứng minh tồn dãy đếm giá trị riêng cách khác nhau: sử dụng định lý Rouche (xem định lý 2.2.1), lý thuyết dao động Sturm (xem định lý 2.3.3), phương pháp phương trình tích phân ( xem bổ đề 2.4.1 chứng minh tồn dãy đếm giá trị riêng giải thức ) Ngồi chương có cách chứng minh khác cho định lý khai triển hàm riêng : phương pháp phương trình tích phân ( xem định lý 2.5.2), phương pháp thặng dư Cauchy (xem định lý 2.6.2 định lý 2.6.3) Ở cuối chương định lý , hội tụ điểm khai triển hàm riêng Sturm-Liouville giống hội tụ điểm chuỗi Fourier thông thường ( xem định lý 2.7.1 nhận xét cuối mục 2.7 ) Ở chương 3, xây dựng hàm phổ r(l) (cịn gọi độ đo phổ) từ định nghĩa biến đổi Fourier tổng quát thu đẳng thức Parseval định lý khai triển dạng tương tự chương (xem định lý 3.1.1, định lý 3.1.3 định lý 3.3.4 ) Đồng thời chương trình bày phân loại giới hạn điểm, giới hạn trịn tốn tử Sturm-Liouville nhiên em chưa tìm hiểu xuất phát điểm vật lý khái niệm Ngoài chương trình bày biểu diễn tích phân giải thức (xem định lý 3.3.3), rõ họ phổ E l toán tử Sturm-Liouville ( xem bổ đề 3.4.3 mệnh đề 3.4.1) Ở cuối chương ánh xạ f (x) 7!F( l) đặt tương ứng hàm f (x) L (0, ¥) biến đổi Fourier tổng quát F(l) L ( ¥, +¥) ánh xạ Unitary ( song ánh bảo toàn r(l) chuẩn) Mặc dù cố gắng, nhiên luận văn khơng tránh khỏi sai sót, mong nhận góp ý q thầy bạn đọc Tài liệu tham khảo [1] Apostol T.M (1974), Mathematical analysis, 2nd, Pearson [2] Al-Gwaiz M.A (2008), Sturm-Liouville and its applications, Springer [3] Brown J.W and Churchill R.V (2013), Complex variables and applica-tions, 9th, McGraw-Hill [4] Coddington E.A and Levinson N (1955), Theory of ordinary differen-tial equations, McGraw-Hill [5] Courant R and Hilbert D (1989), Methods of mathematical physics, vol I, Wiley-VCH [6] EidelmanY., Milman V and Tsolomitis A (2004), Functional analysis an introduction, AMS [7] Freiling G and Yurko V (2001), Inverse Sturm-Liouville problems and their applications, Nova Science [8] Folland G.B (1992), Fourier analysis and its applications, Wadsworth Brooks/ cole [9] Levitan B.M and Sargsjan I.S (1975), Introduction to spectral theory: selfadjoint ordinary differential operators, AMS [10] Levitan B.M and Sargsjan I.S (1991), Sturm-Liouville and Dirac oper-ators, Springer [11] Lebovitz N (2019), Ordinary differential equations, Cengage Learning [12] Stein E and Shakarchi R (2005), Real analysis, Princeton university press TÀI LIỆU THAM KHẢO [13] Teschl G (2012), Ordinary differential equations and Dynamical sys-tems, AMS [14] Teschl G (2014), Mathematical methods in quantum mechanics with applications to Schrodinger operators, 2nd, AMS Titchmarsh E.C (1950), Eigenfunction expansions associated with second-order differential equations part I, Clarendon Oxford [15] ... TRƯỜNG ĐẠI HỌC KHOA HỌC TỰ NHIÊN ——————- Nguyễn Viết Đại HÀM RIÊNG CỦA TOÁN TỬ STURM- LIOUVILLE TRÊN KHOẢNG HỮU HẠN VÀ TRÊN KHOẢNG VƠ HẠN Chun ngành: Tốn giải tích Mã số: 8460101.02 LUẬN VĂN THẠC... tử cụ thể thơng tin tiệm cận giá trị riêng, hàm riêng họ phổ Trong luận văn em đọc hiểu trình bày chi tiết lại kết khai triển hàm riêng toán tử Sturm- Liouville cho hai trường hợp khoảng hữu hạn. .. giá trị riêng tăng vơ tốn tử Sturm- Liouville 2.4 Hàm Green, toán tử compact đối xứng Ở mục trước, ta chứng minh tồn giá trị riêng toán tử Sturm- Liouville định lý Rouche định lý so sánh Sturm Mục