Luận văn sư phạm Phương pháp tính gần đúng tích phân

Lời cảm ơn Trong suốt thời gian học tập Khoa Toán – Trường Đại học Sư Phạm Hà Nội 2, dạy dỗ bảo tận tình thầy giáo, cô giáo, em tiếp thu nhiều tri thức khoa học, kinh nghiệm phương pháp học tập mới, bước đầu làm quen với việc nghiên cứu khoa học Qua đây, em xin chân thành cảm ơn tồn thể thầy giáo, giáo Khoa Toán – người giúp đỡ, chăm lo dìu dắt chúng em trưởng thành hơm Đặc biệt, em xin gửi lời cảm ơn chân thành sâu sắc tới thầy: Tiến sỹ Nguyễn Văn Hùng, người trực tiếp hướng dẫn, bảo đóng góp nhiều ý kiến quý báu suốt thời gian em thực khóa luận Sinh viên Nguyễn Thị Ngân Lời cam đoan Khóa luận em hoàn thành hướng dẫn thầy giáo Nguyễn Văn Hùng với cố gắng thân em trình nghiên cứu thực khóa luận, em có tham khảo tài liệu số tác giả (đã nêu mục tài liệu tham khảo) Em xin cam đoan kết khóa luận kết nghiên cứu thân, không trùng với kết tác giả khác Nếu sai em xin hoàn toàn chịu trách nhiệm Sinh viên Nguyễn Thị Ngân Mục lục Lời cảm ơn Lời cam đoan Mục lục Lời nói đầu Chương 1: Các kiến thức chuẩn bị A SAI SỐ 1.1 Số gần sai số 1.2 Các quy tắc tính sai số 1.3 Sai số tính tốn sai số phương pháp B ĐA THỨC NỘI SUY 1.4 Đa thức nội suy Lagrange C TÍCH PHÂN 14 1.5 Tích phân 14 Chương 2: Giải gần tích phân 16 2.1 Mở đầu 16 2.2 Cơng thức hình thang 16 2.3 Công thức Simpson 21 2.4 Công thức Newton – Cotes 27 2.5 Công thức Chebysev 31 2.6 Công thức Gauss 34 2.7 Giải gần tích phân bội phương pháp Monte – Carlo 38 Chương 3: Ứng dụng 41 Kết luận 49 Tài liệu tham khảo 50 Lời nói đầu Tốn học nhu cầu giải tốn có nguồn gốc từ thực tiễn Cùng với thời gian, toán học ngày phát triển chia thành hai lĩnh vực: Toán học lý thuyết toán học ứng dụng Khi nói đến tốn học ứng dụng khơng thể khơng nhắc đến Giải tích số Giải tích số mơn khoa học nghiên cứu cách giải gần phương trình, tốn xấp xỉ hàm số, toán tối ưu Sự đời phát triển Giải tích số góp phần quan trọng việc tạo thuật giải toán thực tế như: tốn ngược lĩnh vực thăm dị, chuẩn đoán, nhận dạng… Ngày nay, với phát triển Tin học kiến thức Giải tích số trở nên cần thiết Chúng ta chứng kiến xu song song hóa diễn tất lĩnh vực Giải Tích số Để tiết kiệm nhớ máy tính, người ta đề xuất phương pháp hữu hiệu xử lí hệ lớn, thưa kĩ thuật ném ma trận, kĩ thuật tiền xử lí ma trận… Trong thực tế, nhiều ta phải tính tích phân xác định hàm số khơng biết ngun hàm nó, dùng định nghĩa độ xác đạt khơng cao mà phải thực khối lượng tính tốn lớn Ngồi ra, nhiều trường hợp, hàm số cho dạng bảng nên khái niệm nguyên hàm trở nên vơ nghĩa Tuy nhiên, Giải tích số cung cấp cho phương pháp đơn giản để tính gần tích phân xác định mà độ xác khơng Vì vậy, với niềm u thích mơn Giải tích số em lựa chọn đề tài cho khóa luận tốt nghiệp em “ Phương pháp giải gần tích phân” Khóa luận bao gồm chương: Chương 1: Các kiến thức chuẩn bị Chương 2: Giải gần tích phân Chương 3: Bài tập CHƯƠNG 1: CÁC KIẾN THỨC CHUẨN BỊ A SAI SỐ 1.1 Số gần sai số 1.1.1 Số gần – Sai số tương đối sai số tuyệt đối ∗ Ta gọi x số gần số với x ∗ Hiệu số ∆= | ∗ − | gọi sai số thực x khơng biết ∗ giá trị x không sai khác nhiều so nên xác định ∆ Mặt khác ta tìm số ∆ ≥ cho | ∆ = ọ ố ệ đố ∆ ọ | | Suy ∆ = | | ố ươ ủ đố ∗ − | ≤ ∆ Khi đó: ủ cơng thức thể mối liên hệ sai số tương đối sai số tuyệt đối 1.1.2 Quy tròn số sai số quy tròn a, Hiện tượng quy tròn số Khi gặp số có nhiều số đằng sau dấu phẩy, người ta bỏ vài chữ số cuối, việc làm gọi quy trịn số Mỗi quy tròn số, ta tạo sai số gọi sai số quy tròn tuyệt đối b, Sai số quy tròn tuyệt đối ∗ Gọi x sai số gần số cho | − Vì | ∗ − |≤ |≤| ∗ − |+| − Gọi x số gần |≤ số quy trịn x Thế gọi sai số quy trịn trịn số sai số tuyệt đối tăng thêm cho | − ∗ |≤∆ + nên ta thấy làm số quy trịn x Thế số 1.1.3 Cách viết số gần a, Chữ số có nghĩa Chữ số có nghĩa tất chữ số khác khơng, kể số khơng kẹp hai chữ số khác khơng đại diện cho hàng giữ lại Chẳng hạn 0,000014060 có năm chữ số có nghĩa 1;4;0;6;0 b, Chữ số đáng tin Mọi chữ số thập phân x biểu diễn dạng =± 10 số nguyên từ đến ∗ Gọi x chữ số gần với sai số tuyệt đối ∆ Thế gọi chữ số hay chữ số đáng tin ∆ ≤ 0,5 10 ∆ ≥ 0,5 10 chữ số đáng nghi c, Cách viết số gần ∗ Gọi x chữ số gần hai cách viết số gần x Cách 1: ± ∆ (1 ± với sai số tuyệt đối ∆ Thế có ) Cách 2: Viết theo quy ước chữ số có nghĩa x chữ số đáng tin 1.2 Các quy tắc tính sai số 1.2.1 Mở đầu Xét hàm số u hai biến x y có dạng số x;y Hãy lập cơng thức tính sai số u Ta kí hiệu ∆ , ∆ , ∆ số gia x;y;u , , vi phân x;y;u ∆ , ∆ , ∆ sai số tuyệt đối x;y;u = ( , ) Cho biết sai Vì | ∗ − | ≤ ∆ nên ta ln có |∆ | ≤ ∆ ∆ ≤∆ Ta phải tìm ∆ để có |∆ | ≤ ∆ 1.2.2 Sai số tổng = + Ta có ∆ = ∆ + ∆ suy ∆ ≤ |∆ | + |∆ | nên ∆ ≤ ∆ + ∆ Ta chọn ∆ = ∆ + ∆ để có |∆ | ≤ ∆ Do ta có quy tắc: Sai số tuyệt đối tổng tổng sai số tuyệt đối số hạng Chú ý ế = − ù ấ ℎì = ∆ +∆ ∆ = | − | | | Cho nên | − | bé sai số tương đối lớn Vì vậy, tính tốn người ta tìm cách để tránh phải trừ số gần 1.2.3 Sai số tích Ta có ∆ ≈ = = + ≈ ∆ + ∆ Suy |∆ | ≤ | ||∆ | + | ||∆ | ≤ | |∆ + | |∆ Suy ∆ = | |∆ + | |∆ Tức = | |∆ + | |∆ ∆ ∆ ∆ = = + = | | | | | | | || | = = + + Vậy sai số tương đối tích tổng sai số tương đối thừa số tích Đặc biệt ta có ( 1.2.4 Sai số thương ) = = với n nguyên dương , ≠ Tương tự trường hợp tích ta có quy tắc: Sai số tương đối thương tổng sai số tương đối số hạng: ⁄ = + 1.2.5 Công thức tổng quát ℎ = ( , = 1.3 = ,…, ) ∆ ∆ = = | | | | | | ∆ ó ∆ = ∆ = ∆ = ln Sai số tính tốn sai số phương pháp | | ∆ ∆ Khi giải gần toán phức tạp ta phải thay toán cho toán đơn giản giải thơng qua việc thực phép tính thơng thường tay hay máy tính điện tử Phương pháp thay toán phức tạp toán đơn giản gọi phương pháp gần Sai số phương pháp gần tạo gọi sai số phương pháp Để giải toán đơn giản ta phải thực phép tính thơng thường, ta ln phải quy trịn kết không gian Sai số tạo tất lần quy tròn gọi sai số tính tốn Sai số cuối tổng hợp hai loại sai số phương pháp sai số tính tốn Chú ý Sai số tổng hợp cuối có phần sai số phương pháp sai số tính tốn.Vì vậy, phải ý điều chỉnh cho sai số cuối nhỏ sai số cho phép B ĐA THỨC NỘI SUY Trong thực tế, nhiều ta phải tìm hàm y = f(x), biết giá trị y điểm x ∈ [a, b] (i = 0,1, … , n).Cũng có trường hợp biểu thức giải tích f(x) cho q cồng kềnh Khi dùng phép nội suy ta dễ dàng tính f nhiêu ∈ [ , ] mà độ xác khơng bao Ngồi ý nghĩa lịch sử ra, đa thức đại số thường dùng phép nội suy lí đơn giản sau: phép tốn cộng, trừ, nhân, đạo hàm, tích phân dễ dàng thực đa thức Hơn P(x) đa thức, c số P(cx) P(x + c) đa thức Bài toán đặt sau: Cho mốc nội suy a ≤ x < x < ⋯ < x ≤ b Hãy tìm đa thức bậc m: P (x) = a x cho P (x ) = y ≔ f(x ) (i = 0, n) Ý nghĩa hình học tốn nội suy là: xây dựng đường cong đại số y = P (x) qua điểm cho trước (x , y ) (i = 0, n) Như ta cần xác định (m + 1) hệ số a (i = 0, n) từ hệ phương trình tuyến tính sau: a x = y (i = 0, n) (1.1) Dễ thấy m < n (m > n) hệ nói chung vơ nghiệm (vơ định) Khi m = n, hệ (1.1) có định thức Vandermonde x x … x … x = △= x x … … … … … x x … x x −x ≠0 Suy phương trình (1.1) có nghiệm Để giải tốn cho ta phải giải hệ (1.1), khó khăn phức tạp Sau ta trình bày cách xây dựng đa thức nội suy mà không cần giải hệ (1.1), gọi đa thức nội suy Lagrange 1.4 Đa thức nội suy Lagrange 1.4.1 Công thức nội suy Lagrange Tìm đa thức ( ) có bậc n, cho Ta có: Đặt ( )= ( , = 0, ) = ( − )…( − ( − )…( − )( − )( − ( )= ( ), = = )…( − ) )…( − ) ( = 0, ) Như vậy, P(x) đa thức nội suy (duy nhất) cần tìm Nếu mốc nội suy cách đều, tức ℎì đặ ( )= ( ( ≔ − − ℎ = ℎ ( = 0, − 1) ℎ + ℎ) = + ℎ) = 1.4.2 Sai số phép nội suy (−1) − = + ℎ, đượ ( − 1) … ( − ) ! ( − 1) … ( − ) ! (−1) − a) Sai số phương pháp Giả sử, P(x) đa thức nội suy bậc n hàm f(x), tức ( ) = ( ) ( = 0, ) Ta cố định giá trị cách ước lượng sai số ≠ ( = 0, ) ∈ [ , ] tùy ý tìm ( ) = ( ) − ( ) Dĩ nhiên cần xét ( ) = ( = 0, ) 10 Công thức cầu phương Gauss (n = 3) có dạng: f(t)dt ⋍ 5f − + 8f(0) + 5f b−a b+a t + 2 Khi cơng thức (2.10) có dạng: Nếu [a, b] , đặt x = f(x)dx ≃ x = b−a b+a t + , 2 b−a B f(x ), t nghiệm đa thức Legendre bậc n Sai số công thức Gauss là: (b − a) (n!) f ( ) (ξ) R = [(2n!)] (2n + 1) Nhận xét: Cơng thức Gauss có độ xác cao Nhược điểm phương pháp hệ số , nói chung vơ tỉ, tính tốn theo cơng thức (2.10) phức tạp Bây ta nói sai số tính tốn cơng thức cầu phương Giả sử ta có cơng thức cầu phương sau: f(x)dx = A f(x ) + R(f), (2.14) x ∈ [a, b] (i = 1, n) điểm chia, A (i = 1, n) số dương R(f) sai số phương pháp 39 Nếu tính cịn y = f(x ) tính với sai số ε, tức thay cơng thức (2.14) ta có: A y + R(f) (2.15) f(x)dx = |y − y | ≤ ε (i = 1, n) Dễ thấy R(f) = Và A (y − y ) + R(f) |R(f)| ≤ ε A + |R(f)| Vì cơng thức (2.14) cho đa thức bậc không, nên với f(x) ≡ 1, ta có: b−a= dx = A Như sai số tính tốn (2.15) là: |R(f)| ≤ (b − a)ε + |R(f)|, R(f) sai số phương pháp công thức cầu phương (2.14) Dưới bảng yếu tố công thức Gauss với n = 1,2,3,4 N 2 -0,577 -0,775 0,556 0,577 0,889 0,775 0,556 -0,861 0,348 -0,340 0,652 0,340 0,652 40 0,861 0,34888 2.5.2 Các tốn Bài tốn 10 Bằng cơng thức Gauss, với n = 4, giải gần tích phân sau: = ( )= i 1+ Giải 1+ ậ đượ , ( ) ả : ( ) -0,861 0,574276655 0,348 0,199848275 -0,340 0,89637863 0,652 0,584438866 0,340 0,89637863 0,652 0,584438866 0,861 0,574276655 0,348 0,199848275 = 1,568574282 Áp dụng cơng thức Gauss ta được: ≈ 1,568574282 2.7 Tính gần tích phân bội phương pháp Monte – Carlo Nói cách nơm na, phương pháp Monte – Carlo phương pháp giải toán cách sử dụng nhiều phép thử ngẫu nhiên Giả sử ta có hàm y = f(x , … , x ) liên tục giới nội miền W Ta cần tính tích phân bội sau: I= W f(x)dx = W … f(x , … , x )dx … dx Giả sử miền W nằm gọn khối hộp chữ nhật Π: W ⊂ Π ≔ {x ∈ : a ≤ x ≤ A (i = 1, m)} 41 Bằng phép đổi biến x = a + (A − a )ξ (i = 1, m), khối hộp chữ nhật  biến thành khối hộp đơn vị Ta có I = Π ={ ∈ W f(x)dx = : ≤ ξ ≤ (i = 1, n)} F(ξ)dξ , dξ = dξ … dξ , x = ξ F(ξ) = f(a + (A − a )ξ , … , a + (A − a )ξ ) det (A − a )f(a + (A − a )ξ , … , a + (A − a )ξ ) = a) Xét trường hợp ( ) ≥ Gọi V = {(ξ, y) ∈ Giả sử ≤ F(ξ) ≤ B Đặt v = Dễ thấy : ξ ∈ ; ≤ y ≤ F(ξ)}, ta có I = y v = ( , ) ∈ B :ξ ∈ ;0 ≤ V ⊂ {( , ): ≤ ξ ≤ 1; ≤ ≤ dξdy F(ξ) B ≤ 1, i = 1, m} Chọn (m +1) dãy số ngẫu nhiên phân bố [0,1]: ξ( ) ,…, ξ ( ) , {ν }(i = 1,2, … ) xét điểm M ξ ( ) … , ξ ( ) , ν Giả sử N (N V Khi đó: I= dξdy = B 1) điểm ngẫu nhiên M , có n điểm nằm dξd = BP(M ∈ V) ≃ B b) Xét trường hợp hàm f đổi dấu ặt F(ξ) = −b + (B + b)F(ξ) 42 n N ≤ ( )≤ (2.17) Ta có: = −bvol( ) + (B + b) F(ξ) F (ξ)dξ, F (ξ)dξ tính theo sơ đồ ≤ F ≤ tích phân nêu phần (a) Để đánh giá sai số công thức (2.17), coi B = n I = P(M ∈ V) ≃ N Áp dụng bất đẳng thức Chebysev, ta được: P n −I N < ≥1− I (1 − I ) ≥1− ε N 4ε N n = δ P −I < ≥ − δ 4ε N N 1 = hi δ cố định ta thấy ε = 2√δN √N Điều chứng tỏ phương pháp Monte – Carlo hội tụ chậm.Để Nếu chọn tăng độ xác lên 10 lần, khối lượng phép thử phải tăng lên 100 lần 4ε δ N 25 triệu!! Nếu cho trước ε, δ số phép thử N = Ví dụ ε = 10 ; = 10 Tuy nhiên, số phép thử không phụ thuộc vào số lớp lấy tích phân, nên phương pháp Monte – Carlo đặc biệt có lợi việc tính tích phân nhiều chiều, phương pháp tất định tỏ bất lực Ví dụ, để tính tích phân 10 lớp công thức cầu phương khối hộp đơn vị với bước h = 0.1,ta phải tính tổng có chứa tới 10 số hạng Trong phương pháp Monte – Carlo dễ dàng cho ta lời giải thô tốn 43 Chương 3: Ứng dụng A.TÍNH GẦN ĐÚNG MỘT SỐ TÍCH PHÂN Bài 1: Bằng phương pháp hình thang, với việc chia đoạn [3;5] thành 10 phần nhau, tính: = Bài 2: Bằng phương pháp Simpson, với việc chia đoạn [0;1] thành 10 phần nhau, tính: Bài 3: = ℎ = +2 Nếu tính I theo cơng thức hình thang cần chia đoạn [0;1] điểm chia (n =?) để sai số nhỏ 10 Bài 4: ℎ = +3 Tính I theo cơng thức biết với h = 0,1 44 Bài 5: Giải gần tích phân sau công thức Newton – Cotes với n = = Bài 6: 1+ Giải gần tích phân sau phương pháp Gauss với n = 1+ 1+ = Bài 7: Giải gần tích phân sau phương pháp Chebysev với n = = Bài 8: 1+ Giải gần tích phân sau phương pháp Monte-Carlo = = {( , )| Bài 9: ( + + ) , ≤1; ≥0; ≥ 0} Giải gần tích phân sau phương pháp Monte-Carlo = = ( , )| ( − 45 ) , + − ≤ B HƯỚNG DẪN GIẢI Bài 1: Tính tích phân cơng thức hình thang: = ℎ= 5−3 = 0,2 10 Chia đoạn [3;5] thành 10 phần Ta bảng sau: 3,0 0,708648150 3,2 0,954116007 3,4 0,974218425 3,6 0,716702360 3,8 0,306467225 4,0 0,019738581 4,2 0,103543674 4,4 0,533616051 4,6 0,944659613 4,8 0,912853014 5,0 0,425999184 Áp dụng công thức hình thang: = 0,2 [ + + 2( = 1,206647723 46 + + + ⋯+ )] Bài 2: Tính tích phân cơng thức Simpson: = Chia đoạn [0;1] thành 10 phần ℎ= Ta tính tốn bảng sau: 1−0 = 0,1 10 sin 0,0 0,000000000 0,1 0,003999989 0,2 0,015999317 0,3 0,035992224 0,4 0,063956318 0,5 0,099833416 0,6 0,143502851 0,7 0,194747485 0,8 0,253212945 0,9 0,318360975 1,0 0,389418342 Áp dụng công thức Simpson: = 0,1 [ + + 4( + 2( = 0,131816585 + 47 + + + + + + ) )] Bài 3: Phải chia [0;1] thành phần để = Dễ dàng tính được: Nên ≤ 10 ? +2 ( + 2) − 2( + 1)( + 2) ( + 2) ( )= ) ≤ (∀ ∈ [0,1]) ℎư ậ : ( − ) để = ℎ ≤ 10 ℎì: 12 ′′ ( ℎ≤ 4,8 10 ℎ ≤ 4,8 10 ≈ 0,069282032 ℎ= − ⇒ = − ℎ ≥ 14,43375679 Vậy với n=15,16,….thì thỏa mãn u cầu tốn Bài 4: Tính tích phân công thức biết: ℎ ℎ = 0,1 = = +3 − ℎ 48 = 1−0 = 10 0,1 Chia đoạn [0;1] thành 10 phần ta bảng sau: +3 0,0 0,333333333 0,1 0,3125 0,2 0,294117647 0,3 0,277777777 0,4 0,263157894 0,5 0,25 0,6 0,238095238 0,7 0,227272727 0,8 0,217391304 0,9 0,208333333 1,0 0,166666666 Áp dụng cơng thức hình thang: = 0,1 [ + + 2( + = 0,152588383 + )] + ⋯+ Áp dụng công thức Simpson: = Bài 5: 0,1 [ + + 4( + 2( = 0,254301985 + + + + + + + ) )] Tính tích phân cơng thức Newton – Cotes với n =6 : = 1+ 49 ℎ= 1 ⁄6 ⁄2 ⁄3 1−0 = 6 41 41 216⁄217 216 215,0046803 27 26, 03571429 8⁄9 272 241,7777778 27 20,82857143 216⁄341 216 136,8211144 41 20,5 27⁄28 ⁄3 27⁄35 ⁄6 1⁄2 = Bài 6: = 701,9677862 840 = 0,835675936 Giải gần tích phân sau phương pháp Gauss với n = = Ta tính tốn bảng sau: ( ) 1+ 1+ ( ) -0,861 1,10487105 0,348 0,384495125 -0,340 1,011800326 0,652 0,659693812 0,340 1,011800326 0,652 0,659693812 0,861 1,10487105 0,348 0,384495125 = 2,088377875 = 2,088377875 50 Bài 7: Giải gần tích phân sau phương pháp Chebysev với n = = 1+ ( )= sin 1+ -0,832 -4,400474083 -0,374 -0,58361316 0 0,374 0,265896534 0,832 0,272194812 = = −4,445995896 10 = −1,396750805 51 Kết luận Trên em trình bày xong tồn khóa luận “ phương pháp giải gần tích phân” Khóa luận cung cấp số phương pháp tính tích phân cách nhanh chóng dễ dàng Cùng với ví dụ minh họa cụ thể chọn lọc kĩ lưỡng, khóa luận giúp bạn đọc tiếp cận với mơn “ Giải tích số” coi tài liệu tham khảo Tuy nhiên, thời gian nghiên cứu hạn chế, phạm vi nghiên cứu tương đối rộng nên khóa luận khơng thể tránh khỏi thiếu xót Rất mong nhận ý kiến đóng góp q thầy bạn đọc Một lần cho phép em gửi lời cảm ơn tới tất thầy cô giảng viên trường, cán thư viện nhà trường đặc biệt thầy giáo Nguyễn Văn Hùng tận tình giúp đỡ em hoàn thành xong đề tài Em xin chân thành cảm ơn 52 Tài liệu tham khảo Phạm Kỳ Anh – Giải tích số – NXB Đại học Quốc Gia Hà Nội – 1996 Phạm Phú Chiêm, Nguyễn Bường – Giải tích số – NXB Đại học Quốc Gia Hà Nội – 2000 Nguyễn Minh Chương, Nguyễn Văn Khải, Khuất Văn Ninh, Nguyễn Văn Tuấn, Nguyễn Tường – Giáo trình Giải tích số – NXB Giáo dục – 2000 Tạ Văn Đĩnh – Phương pháp tính – NXB Giáo dục – 1998 Phan Văn Hạp, Hồng Đức Ngun, Lê Đình Thịnh – Bài tập phương pháp tính – NXB Khoa Học Kĩ Thuật Hà Nội – 1996 Hoàng Xuân Huấn – Giáo trình phương pháp số – NXB Đại học Quốc Gia Hà Nội Dương Thủy Vỹ – Giáo trình phương pháp tính – NXB Khoa Học Kĩ Thuật Hà Nội 53 ... số lớp lấy tích phân, nên phương pháp Monte – Carlo đặc biệt có lợi việc tính tích phân nhiều chiều, phương pháp tất định tỏ bất lực Ví dụ, để tính tích phân 10 lớp công thức cầu phương khối... Giải gần tích phân sau phương pháp Chebysev với n = = Bài 8: 1+ Giải gần tích phân sau phương pháp Monte-Carlo = = {( , )| Bài 9: ( + + ) , ≤1; ≥0; ≥ 0} Giải gần tích phân sau phương pháp Monte-Carlo... Giải tích số cung cấp cho phương pháp đơn giản để tính gần tích phân xác định mà độ xác khơng Vì vậy, với niềm u thích mơn Giải tích số em lựa chọn đề tài cho khóa luận tốt nghiệp em “ Phương pháp