ĐẠI HỌC THÁI NGUYÊN TRƯỜNG ĐẠI HỌC KHOA HỌC HOÀNG THỊ NGA MỘT SỐ ỨNG DỤNG CỦA CÔNG THỨC NỘI SUY LAGRANGE VÀ HERMITE LUẬN VĂN THẠC SĨ TOÁN HỌC THÁI NGUYÊN 2018 Tai ngay!!! Ban co the xoa dong chu nay!![.]
ĐẠI HỌC THÁI NGUYÊN TRƯỜNG ĐẠI HỌC KHOA HỌC HOÀNG THỊ NGA MỘT SỐ ỨNG DỤNG CỦA CÔNG THỨC NỘI SUY LAGRANGE VÀ HERMITE LUẬN VĂN THẠC SĨ TOÁN HỌC THÁI NGUYÊN - 2018 Tai ngay!!! Ban co the xoa dong chu nay!!! ĐẠI HỌC THÁI NGUYÊN TRƯỜNG ĐẠI HỌC KHOA HỌC HOÀNG THỊ NGA MỘT SỐ ỨNG DỤNG CỦA CÔNG THỨC NỘI SUY LAGRANGE VÀ HERMITE Chuyên ngành: PHƯƠNG PHÁP TOÁN SƠ CẤP Mã số: 46 01 13 LUẬN VĂN THẠC SĨ TOÁN HỌC Người hướng dẫn khoa học: GS.TSKH Nguyễn Văn Mậu THÁI NGUYÊN - 2018 i MỤC LỤC MỞ ĐẦU Chương Nội suy Lagrange nội suy 1.1 Bài toán nội suy Lagrange 1.2 Bài toán nội suy Hermite 1.3 Bài toán nội suy Lagrange - Newton ii Hermite 18 Chương Ứng dụng nội suy tính nguyên hàm tích phân hàm phân thức 21 2.1 Nguyên hàm hàm phân thức với cực điểm đơn 21 2.2 Nguyên hàm hàm phân thức với cực điểm bậc tùy ý 26 Chương Một số dạng toán liên quan 3.1 Một số toán đa thức nhận giá trị nguyên 3.2 Một số toán xác định đa thức 3.2.1 Tìm đa thức biết nghiệm 3.2.2 Sử dụng công thức nội suy Lagrange để xác định hệ số đa thức 3.2.3 Một số toán xác định đa thức khác không liên quan đến công thức nội suy 43 43 50 50 53 56 KẾT LUẬN 59 TÀI LIỆU THAM KHẢO 59 ii MỞ ĐẦU Trong kì thi học sinh giỏi toán cấp, Olympic Toán sinh viên, toán liên quan tới đa thức thường xuyên đề cập Những dạng toán thường xem thuộc loại khó, phần kiến thức nội suy đa thức lại khơng nằm chương trình thức giáo trình Đại số Giải tích bậc trung học phổ thông Như ta biết, công thức nội suy Lagrange đề cập bậc phổ thơng Tuy nhiên cơng thức nội suy Hermite có tài liệu chun khảo Vì vậy, tơi chọn đề tài luận văn ”Một số ứng dụng công thức nội suy Lagrange Hermite” Luận văn nhằm cung cấp số ứng dụng công thức nội suy Lagrange Hermite để tìm nguyên hàm hàm phân thức Luận văn gồm phần mở đầu, kết luận chương Chương Nội suy Lagrange nội suy Hermite Chương Ứng dụng nội suy tính nguyên hàm tích phân hàm phân thức Chương Một số dạng toán liên quan Tiếp theo, chương trình bày hệ thống tập áp dụng giải đề thi HSG quốc gia Olympic liên quan Luận văn thực hoàn thành Trường Đại học Khoa học, Đại học Thái Nguyên hướng dẫn GS.TSKH Nguyễn Văn Mậu Xin gửi lời cảm ơn chân thành sâu sắc đến Thầy, người tận tình hướng dẫn đạo tác giả tập dượt nghiên cứu khoa học suốt q trình tìm hiểu tài liệu, viết hồn thiện Luận văn Đồng thời em xin chân thành cảm ơn q thầy Bộ mơn tốn, Khoa Khoa học Tự nhiên, Thầy Cô trường Đại học Khoa học Tự nhiên Hà Nội, Thầy Cơ Viện Tốn học tận tình giảng dạy, quan tâm tạo điều kiện thuận lợi thủ tục hành để em hồn thành khố học bảo vệ luận văn Thạc sĩ Tôi chân thành cảm ơn gia đình, bạn bè quan, đồn thể nơi tơi công tác Trường Trung học Phổ thông Thuỷ Sơn, Sở Giáo dục Đào tạo Hải Phòng, tạo điều kiện vật chất lẫn tinh thần trình học tập, nghiên cứu viết luận văn Thái Nguyên, tháng năm 2018 Tác giả Hoàng Thị Nga Chương Nội suy Lagrange nội suy Hermite Chương dành để trình bày toán nội suy Lagrange, toán nội suy Hermite tốn nội suy Lagrange-Newton, từ định lí, hệ số ví dụ tính tốn cụ thể 1.1 Bài toán nội suy Lagrange Trong số trường hợp, để tính tổng hữu hạn phân thức, người ta thường sử dụng số tính chất đa thức, đặc biệt công thức nội suy Lagrange Dưới số đồng thức áp dụng chúng Định lý 1.1 (Đồng thức Lagrange) Nếu x1 , x2 , , xm m giá trị tuỳ ý, đôi khác f (x) đa thức bậc nhỏ thua m ta có đồng thức sau (x − x2 )(x − x3 ) (x − xm ) + f (x) = f (x1 ) (x1 − x2 )(x1 − x3 ) (x1 − xm ) +f (x2 ) (x − x1 )(x − x3 ) (x − xm ) (x2 − x1 )(x2 − x3 ) (x2 − xm ) + · · · + f (xm ) (x − x1 )(x − x2 ) (x − xm−1 ) (xm − x1 )(xm − x2 ) (xm − xm−1 ) (1.1) Chứng minh Ta cần chứng minh công thức f (x) − f (x1 ) −f (x2 ) − · · · − f (xm ) (x − x2 )(x − x3 ) (x − xm ) − (x1 − x2 )(x1 − x3 ) (x1 − xm ) (x − x1 )(x − x3 ) (x − xm ) (x2 − x1 )(x2 − x3 ) (x2 − xm ) (x − x1 )(x − x2 ) (x − xm−1 ) ≡ (xm − x1 )(xm − x2 ) (xm − xm−1 ) Nhận xét vế trái công thức đa thức bậc không m − có m nghiệm phân biệt x1 , x2 , , xm Vậy đa thức phải đồng Hệ 1.1 Từ Định lý 1.1, ta thu đồng thức sau √ √ √ √ √ √ (x − 3)(x − 5)(x − 7) (x − 2)(x − 5)(x − 7) √ √ √ √ √ √ + √ √ √ √ √ √ ( − 3)( − 5)( − 7) ( − 2)( − 5)( − 7) √ √ √ √ √ √ (x − 2)(x − 3)(x − 7) (x − 2)(x − 3)(x − 5) √ √ √ √ √ + √ √ √ √ √ √ ≡ 1, + √ ( − 2)( − 3)( − 7) ( − 2)( − 3)( − 5) a2 (x − c)(x − a) (x − a)(x − b) (x − b)(x − c) + b2 + c2 ≡ x2 (a < b < c) (a − b)(a − c) (b − c)(b − a) (c − a)(c − b) Định lý 1.2 Giả sử f (x) đa thức bậc nhỏ thua m − (m > 2) x1 , x2 , , xm m giá trị đôi khác cho trước tuỳ ý Khi đó, ta có đồng thức f (x1 ) f (x2 ) + (x1 − x2 )(x1 − x3 ) (x1 − xm ) (x2 − x1 )(x2 − x3 ) (x2 − xm ) f (xm ) ≡ + ··· + (xm − x1 )(xm − x2 ) (xm − xm−1 ) Chứng minh Nhận xét vế trái đẳng thức cho hệ số hạng tử ứng với bậc m − cách viết tắc đa thức f (x) Đồng hệ số đồng bậc ta có điều phải chứng minh Dưới đây, ta xét số ứng dụng trực tiếp đồng thức Lagrange Ví dụ 1.1 Tính tổng cos 1o cos 2o S= + (cos 1o − cos 2o )(cos 1o − cos 3o ) (cos 2o − cos 1o )(cos 2o − cos 3o ) cos 3o + (cos 3o − cos 1o )(cos 3o − cos 2o ) Lời giải Áp dụng Định lý 1.2, với f (x) = x, x1 = cos 1o , x2 = cos 2o , x3 = cos 3o , ta thu S = Ví dụ 1.2 Ta có đồng thức c+d+a b+c+d + + (b − a)(c − a)(d − a)(x − a) (c − b)(d − b)(a − b)(x − b) d+a+b a+b+c + (d − c)(a − c)(b − c)(x − c) (a − d)(b − d)(c − d)(x − d) x−a−b−c−d ≡ (x − a)(x − b)(x − c)(x − d) Lời giải Thật vậy, ta cần chứng minh (a + b + c + d) − a (a + b + c + d) − b + + (a − b)(a − c)(a − d)(a − x) (b − a)(b − c)(b − d)(b − x) + (a + b + c + d) − d (a + b + c + d) − c + + (c − a)(c − b)(c − d)(c − x) (d − a)(d − b)(d − c)(d − x) + (a + b + c + d) − x = (x − a)(x − b)(x − c)(x − d) Ta có, với đa thức bậc f (y) = a + b + c + d − y, y1 = a, y2 = b, y3 = c, y4 = d, y5 = x, theo Định lý 1.2 ta thu điều phải chứng minh Định lý 1.3 Cho x1 , x2 , , xm m giá trị tuỳ ý đôi khác Đặt Sn = xn1 xn2 + (x1 − x2 )(x1 − x3 ) (x1 − xm ) (x2 − x1 )(x2 − x3 ) (x2 − xm ) xnm +··· + (xm − x1 )(xm − x2 ) (xm − xm−1 ) Khi a) Sn = ≤ n < m − 1, b) Sm−1 = 1, c) Sm+k tổng tích, tích có k + thừa số (giống khác nhau) lấy số x1 , x2 , , xm Chứng minh a) Theo Định lý 1.2, với f (x) = 1, x, x2 , , xm−2 , ta S0 = S1 = = Sm−2 = b) Để chứng minh Sm−1 = 1, ta cần thay f (x) Định lý 1.2 xm−1 , so sánh hệ số hạng tử bậc m − hai vế đồng thức vừa thu c) Để tính Sn n > m − ta làm sau: Giả sử x1 , x2 , , xm thoả mãn phương trình bậc m αm + p1 αm−1 + p2 αm−2 + · · · + pm−1 α + pm = 0, −p1 p = x1 + x2 + · · · + xm = x1 x2 + x1 x3 + · · · + xm−1 xm (−1)k p k = x1 x2 x3 xk + · · · Nhân hai vế phương trình với αk , ta αm+k + p1 αm+k−1 + p2 αm+k−2 + · · · + pm−1 αk+1 + pm αk = Thay α đẳng thức x1 , x2 , , xm ; chia đẳng thức thứ cho (x1 − x2 )(x1 − x3 ) (x1 − xm ), đẳng thức thứ hai cho (x2 − x1 )(x2 − x3 ) (x2 − xm ) , cộng vế với vế đẳng thức vừa nhận được, ta thu Sm+k + p1 Sm+k−1 + · · · + pm−1 Sk+1 + pm Sk = (1.2) Đặt k = 0, ta thu Sm + p1 Sm−1 = Do Sm = −p1 = x1 + x2 + · · · + xm Nhờ đẳng thức (1.2) ta tính tiếp biểu thức Sm+1 , Sm+2 , Ta đặt = α1 ; (x1 − x2 )(x1 − x3 ) (x1 − xm ) = α2 ; (x2 − x1 )(x2 − x3 ) (x2 − xm ) = αm (xm − x1 )(xm − x2 ) (xm − xm−1 ) Khi ta có Sn = xn1 α1 + xn2 α2 + · · · + xnm αm Xét α1 α2 αm + + ··· + − x1 z − x2 z − xm z Dùng công thức cấp số nhân với giả thiết z chọn cho P = |x1 z| < 1, |x2 z| < 1, , |xm z| < 1, ta khai triển tổng P thành chuỗi vô hạn sau: P = α1 (1 + x1 z + x21 z + · · · ) + α2 (1 + x2 z + x22 z + · · · )+ + · · · + αm (1 + xm z + x2m z + · · · ) hay P = (α1 + α2 + · · · + αm ) + (x1 α1 + x2 α2 + · · · + xm αm )z + +(x21 α1 + x22 α2 + · · · + x2m αm )z + · · · tức P = S0 + S1 z + S2 z + S3 z + · · · Để cho gọn, ta đặt (1 − x1 z)(1 − x2 z) (1 − xm z) = Q Khai triển Q theo luỹ thừa z , ta viết Q = − δ1 z + δ2 z + · · · + (−1)m δm z m , δ1 = x1 + x2 + · · · + xm , δ2 = x1 x2 + x1 x3 + · · · + xm−1 xm Tiếp theo, nhân hai vế đẳng thức thứ hai với (1 − x1 z)(1 − x2 z) (1 − xm z), ta có P Q = α1 (1 − x2 z)(1 − x3 z) (1 − xm z)+ α2 (1 − x1 z)(1 − x3 z) (1 − xm z)+ α3 (1 − x1 z)(1 − x2 z)(1 − x4 z) (1 − xm z) + · · · + αm (1 − x1 z)(1 − x2 z) (1 − xm−1 z) Như P Q đa thức bậc m − z Ta chứng minh z m−1 , tức có đồng thức P Q = z m−1 Thật vậy, biểu thức P Q − z m−1 triệt tiêu z= 1 , , , x1 x2 xm vì, chẳng hạn, với z= x1 1 x2 x3 xm α1 − 1− ··· − − m−1 = m−1 − m−1 = x1 x1 x1 x1 x1 x1 Vậy nên P Q − z m−1 = Do z m−1 =P Q hay z m−1 = S0 + S1 z + · · · + Sm−1 z m−1 + · · · m m − δ1 z + δ2 z − + (−1) δm z Nếu khai triển vế trái thành chuỗi vô hạn theo luỹ thừa z chuỗi bắt đầu hạng tử chứa z m−1 Vì vậy, hệ số hạng tử bậc 0, 1, 2, , m − vế phải không, tức S0 = S1 = = Sm−2 = Ngoài ra, hệ số hạng tử ứng với bậc m = vế trái Vậy Sm−1 = Bây đẳng thức cần chứng minh có dạng sau: z m−1 = z m−1 + Sm z m + Sm+1 z m+1 + · · · m m − δ1 z + δ2 z − · · · (−1) δm z Uớc lượng hai vế cho z m−1 , ta thu = + Sm z + Sm+1 z + · · · m m − δ1 z + δ2 z − · · · + (−1) δm z = (1 − δ1 z + δ2 z − · · · + (−1)m δm z m )(1 + Sm z + Sm+1 z + · · · ) (1 − qx ) < an < q q q q Pn (x) = [(2x − 1)n + 1] Ta thấy hệ số Pn (x) số nguyên 18 Bài toán 1.7 Chứng minh không tồn đa thức f (x) ∈ Z[x] mà f (2005) = 2005 f (2007) = 2008 Lời giải Giả sử tồn đa thức f (x) = an xn + an−1 xn−1 + · · · + a0 , ∈ Z ∀i ∈ {0, 1, , n} thỏa mãn điều kiện tốn Khi đó, ta có f (2007) − f (2005) = an (2007n − 2005n ) + an−1 (2007n−1 − 2005n−1 ) + · · · + a1 (2007 − 2005) chia hết cho Mặt khác f (2007) − f (2005) = 2008 − 2005 = không chia hết cho Vậy không tồn đa thức f (x) thỏa mãn điều kiện đề 1.3 Bài toán nội suy Lagrange - Newton Trước hết ta xét toán nội suy hỗn hợp Lagrange - Newton Bài toán 1.8 (Nội suy Lagrange - Newton) Cho xki , aki ∈ R, với xki 6= xkj ∀i 6= j; k = 0, 1, , n − 1; i, j = 1, , rk+1 ; r0 = 0, r0 + r1 + · · · + rk = sk , r0 + r1 + · · · + rn = sn = N Hãy xác định đa thức f (x) có bậc deg f (x) N − thỏa mãn điều kiện f (sk ) (xki ) = aki , ∀k = 0, 1, , n − 1, ∀i = 1, , rk+1 (1.7) Lời giải Ký hiệu rk+1 Lki (x) = x − xkj , (k = 0, 1, , n − 1) phép lấy nguyên hàm x − x ki kj j=1, j6=i Y Khi đó, dễ thấy deg Lki (x) rk+1 − 1, 1 i = j Lki (xkj ) = 0 i 6= j hay Lki (xkj ) = δij 19 Tiếp theo, ta đặt yn (x) ≡ xây dựng dãy hàm rk+1 X yk (x) = Rsk (s ) k [aki − yk+1 (xki )]Lki (x), (k = n − 1, , 1, 0) i=1 Dễ thấy deg yk (x) N − deg yk (x) < sk+1 Bây giờ, ta chứng minh đa thức f (x) = y0 (x) + y1 (x) + · · · + yn−1 (x) nghiệm toán nội suy Lagrange - Newton (1.7) Thật vậy, dễ thấy deg f (x) N − Ngồi ra, ta có f (sn−1 ) (x(n−1)i ) = (y0 + · · · + yn−2 + yn−1 )(sn−1 ) (x(n−1)i ) (s ) n−1 (x(n−1)i ) = yn−1 rn X = a(n−1)j L(n−1)j (x(n−1)i ) = j=1 rn X a(n−1)j δij = a(n−1)i j=1 f (sn−2 ) (x(n−2)i ) = (y0 + · · · + yn−3 + yn−2 + yn−1 )(sn−2 ) (x(n−2)i ) (s ) (s ) n−2 n−2 = yn−2 (x(n−2)i ) + yn−1 (x(n−2)i ), (sn−2 ) yn−2 (x(n−2)i ) rn−1 X (sn−2 ) = [a(n−2)j − yn−1 (x(n−2)j )].L(n−2)j (x(n−2)i ) j=1 rn−1 X (sn−2 ) = [a(n−2)j − yn−1 (x(n−2)j )].δij j=1 (s ) n−2 = a(n−2)i − yn−1 (x(n−2)i ) Vậy nên f (sn−2 ) (x(n−2)i ) = a(n−2)i Bằng cách tương tự, ta chứng minh f (sk ) (xki ) = aki , ∀k = 0, 1, , n − 1, ∀i = 1, , rk+1 20 Cuối cùng, có đa thức f∗ (x) với bậc deg f∗ (x) N − thỏa mãn điều kiện tốn (1.7), đa thức P (x) = f (x) − f∗ (x) có bậc deg P (x) N − thỏa mãn điều kiện P (sk ) (xki ) = 0, ∀k = 0, 1, , n − 1, ∀i = 1, , rk+1 Khi đó, theo cách xây dựng đa thức f (x) phần trên, ứng với trường hợp aki = 0, ∀k = 0, 1, , n − 1, ∀i = 1, , rk+1 , ta có yk = 0, ∀k = 0, 1, , n − Suy P (x) ≡ 0, f (x) = f∗ (x) Từ tốn ta có nhận xét sau: Nhận xét 1.2 Nếu k = 0, tương ứng r2 = r3 = · · · = rn = 0, sk = Khi tốn nội suy Lagrange - Newton tốn nội suy Lagrange quen biết Nếu i = 1, tương ứng r1 = r2 = · · · = rn = 1, sk = k Khi tốn nội suy Lagrange - Newton tốn nội suy Newton biết 21 Chương Ứng dụng nội suy tính nguyên hàm tích phân hàm phân thức Chương giới thiệu cho số thuật tốn tìm ngun hàm hàm số hữu tỉ cho trước theo thứ tự từ đơn giản đến phức tạp 2.1 Nguyên hàm hàm phân thức với cực điểm đơn Lớp hàm quen thuộc có nguyên hàm sơ cấp lớp hàm số hữu tỉ Kết luận Laplace chứng minh lần vào năm 1812 cách phân tích hàm số hữu tỉ thành tổng phân thức đơn giản Định lý 2.1 (Laplace, 1812) Nguyên hàm hàm số hữu tỉ hàm sơ cấp (đó hàm số hữu tỉ, tổng hàm số hữu tỉ số hữu hạn logarit hàm số hữu tỉ) Chứng minh Nhận xét đa thức Q(x) = b0 xn + b1 xn−1 + · · · + bn biểu diễn sau Q(x) = b0 (x − a1 )m1 (x − a2 )m2 (x − ar )mr , m1 , , mr số nguyên dương có tổng n a1 , , ar số thực số phức Do hàm số hữu tỉ R(x) có mẫu đa thức Q(x) biễu diễn dạng ( ) r X β β β s,1 s,2 s,ms A0 xp + A1 xp−1 + · · · + Ap + + ms + ··· + x − a (x − a ) (x − a ) s s s s=1 Từ đó, suy R xp+1 xp R (x) dx = A0 + A1 + · · · + Ap x + C p(+ p r P βs,2 + βs,1 ln (x − as ) − − ··· x − as s=1 ) βs,ms − (ms − 1) (x − as )ms −1 22 hàm số sơ cấp Xét dạng toán minh họa sau 2x + hàm số hữu tỉ Nguyên hàm x3 − 3x + có dạng tổng hàm số hữu tỉ logarit hàm số hữu tỉ.Thật vậy, ta có Z Z 2x + 2x + dx = dx x3 − 3x + (x − 1)2 (x + 2) Z Z Z dx − dx = dx + (x − 1) (x + 2) 3(x − 1) x − =− + ln +C (x − 1) x + Bài toán 2.1 Ta có f (x) = C số Nhận xét 2.3 Nếu m n số ngun, ta có kết luận sau R i) Tích phân (1 − xn ) m dx hàm số sơ cấp m = ±1 n = ±1 m = −n R ii) Tích phân sinm xcosn xdx hàm số sơ cấp với giá trị m n Thật vậy, m, n số nguyên nên ta có n−1 R R sinm x.cosn xdx = sinm x cos2 x cos xdx n−1 R = sinm x − sin2 x cos xdx n−1 R m = u − u2 du, (u = sin x) Tích phân hàm số sơ cấp m lẻ n lẻ hai n, m chẵn Định lý Laplace 2.1 cho ta kết luận hàm hữu tỉ có nguyên hàm sơ cấp Tuy nhiên, việc tìm ngun hàm khơng phải công việc đơn giản Sau số thuật toán giải toán Trong mục này, ta xét thuật tốn tìm ngun hàm hàm số hữu tỉ có P (x) dạng f (x) = Q P (x) đa thức, xi 6= xj với i 6= j n (x − xi ) i=1 Từ cơng thức trên, ta có phương pháp tìm nguyên hàm hàm số hữu tỉ mà mẫu thức bậc n có n nghiệm đơn phân biệt