giáo trình phương pháp tính đại học dong a sinhvienzone com

BỘ GIÁO DỤC VÀ ĐÀO TẠO TRƯỜNG ĐẠI HỌC ĐÔNG Á ne C om ThS.PHẠM THỊ NGỌC MINH Si nh Vi en Zo GIÁO TRÌNH PHƯƠNG PHÁP TÍNH LƯU HÀNH NỘI BỘ Đà Nẵng, 2013 SinhVienZone.com https://fb.com/sinhvienzonevn Môn: Phương pháp tính CHƯƠNG.1 SAI SỐ 1.1 NHẬP MƠN PHƯƠNG PHÁP TÍNH 1.1.1 Giới thiệu mơn phương pháp tính Phương pháp tính mơn tốn học có nhiệm vụ giải đến kết số cho tốn, cung cấp phương pháp giải cho toán thực tế mà khơng có lời giải xác Mơn học cầu nối toán học lý thuyết ứng dụng thực tế .C om Trong thời đại tin học việc áp dụng phương pháp tính trở nên phổ biến nhằm tăng tốc độ tính tốn 1.1.2 Nhiệm vụ mơn học - Tìm phương pháp giải cho toán gồm: phương pháp (PP) phương pháp gần + Phương pháp: kết dạng biểu thức giải tích cụ thể Zo ne + Phương pháp gần đúng: thường cho kết sau q trình tính lặp theo quy luật đó, áp dụng trường hợp tốn khơng có lời giải có phức tạp - Xác định tính chất nghiệm en - Giải toán cực trị Vi - Xấp xỉ hàm: khảo sát, tính tốn hàm f(x) phức tạp, ta thay hàm f(x) hàm g(x) đơn giản cho g(x) ≈ f(x) Việc lựa chọn g(x) gọi phép xấp xỉ hàm Si nh - Đánh giá sai số: giải tốn phương pháp gần sai số xuất sai lệch giá trị nhận với nghiệm thực tốn Vì ta phải đánh giá sai số để từ chọn phương pháp tối ưu 1.1.3 Trình tự giải tốn phương pháp tính - Khảo sát, phân tích tốn - Lựa chọn phương pháp dựa vào tiêu chí sau: + Khối lượng tính tốn + Đơn giản xây dựng thuật tốn + Sai số bé + Khả thi - Xây dựng thuật tốn: sử dụng ngơn ngữ giả sơ đồ khối (càng mịn tốt) - Viết chương trình: sử dụng ngơn ngữ lập trình (C, C++, Pascal, Matlab,…) SinhVienZone.com https://fb.com/sinhvienzonevn Mơn: Phương pháp tính - Thực chương trình, thử nghiệm, sửa đổi hồn chỉnh 1.2 SAI SỐ 1.2.1 Khái niệm Giả sử x số gần x* (x* : số đúng), ∆ = x − x * gọi sai số thực x Vì khơng xác định ∆ nên ta xét đến loại sai số sau: - Sai số tương đối : δ x = ∆x x 1.2.2 Các loại sai số C om - Sai số tuyệt đối : Giả sử ∃∆x > đủ bé cho x − x * ≤ ∆x Khi ∆x gọi sai số tuyệt đối Dựa vào nguyên nhân gây sai số, ta có loại sau: ne - Sai số giả thiết: xuất việc giả thiết toán đạt số điều kiện lý tưởng nhằm làm giảm độ phức tạp toán Zo - Sai số số liệu ban đầu: xuất việc đo đạc cung cấp giá trị đầu vào khơng xác - Sai số phương pháp : xuất việc giải toán phương pháp gần en Vi - Sai số tính tốn : xuất làm tròn số q trình tính tốn, q trình tính nhiều sai số tích luỹ lớn nh 1.2.3 Sai số tính toán Si Giả sử dùng n số gần x i = (i = 1, n) để tính đại lượng y, với y = f ( xi ) = f ( x1 , x2 , , xn ) Trong : - f hàm khả vi liên tục theo đối số xi Khi sai số y xác định theo công thức sau : n - Sai số tuyệt đối : ∆y = ∑ i =1 n - Sai số tương đối : δy=∑ i =1 ∂f ∆xi ∂xi ∂ ln f ∆xi ∂xi - Trường hợp f có dạng tổng : y = f ( xi ) = ± x1 ± x2 ± ± xn SinhVienZone.com https://fb.com/sinhvienzonevn Mơn: Phương pháp tính n ∂f = ∀i suy ∆y = ∑ ∆xi ∂xi i =1 - Trường hợp f có dạng tích : y = f ( xi ) = x1 * x2 * * xn ln f = ln( x1x2 xn ) = (ln x1 + ln x2 + + ln xn ) n n ∆x ∂ ln f = ∀i suy δ y = ∑ i = ∑ δ xi ∂xi xi i =1 xi i =1 n Vậy δ y = ∑ δ xi - Trường hợp dạng thương: y = f ( x) = x1 x2 ne ∂y ∂y − x1 = ; = ∂x1 x2 ∂x2 x2 C om i =1 −x x ∆x + x ∆x ∆x1 + ∆x2 = 2 x2 x2 x22 ⇒δ y = ∆y ∆x1 ∆x2 = + = δ x1 + δ x2 y x1 x2 en Zo ⇒ ∆y = Vi - Trường hợp dạng lũy thừa : y = f ( x) = xα (α > 0) lny = lnf = αlnx Si nh ∂ ln f α = ∂x x Suy δ y = α ∆x = αδ x x Ví dụ 1.1: Cho a ≈ 10,25 ; b ≈ 0,324 ; c ≈ 12,13 Tính sai số : a3 y1 = ; b c y2 = a − b c Giải : SinhVienZone.com https://fb.com/sinhvienzonevn Si nh Vi en Zo ne C om Mơn: Phương pháp tính SinhVienZone.com https://fb.com/sinhvienzonevn Mơn: Phương pháp tính BÀI TẬP CUỐI CHƯƠNG 1 Nêu khái niệm sai số tuyệt đối Nêu khái niệm sai số tương đối Dựa vào nguyên nhân gây sai số, trình bày loại sai số Trình bày sai số tuyệt đối f hàm có dạng tổng Trình bày sai số tương đối f hàm có dạng tích Cho a ≈ 10.25, b ≈ 0.324, c ≈ 12.13 Tính sai số y = a2 b c C om Cho a ≈ 10.25, b ≈ 0.324, c ≈ 12.13 Tính sai số y = a − b c en Zo ne Tính thể tích khối cầu có đường kính d = 3.7cm π = 3.14 ± 0.0016 Một hình trụ có bán kính R = 2m, chiều cao h = 3m Hỏi ∆R ∆h để thể tích V có độ xác ∆V = 0.1m3 ? 10 Một hình cầu có bán kính đáy R = 5.87cm Tính thể tích hình cầu với độ xác 0.01cm3? 11 Xác định sai số tuyệt đối số gần sau biết sai số tương đối chúng: a = 35,72; δ a = 1% b = 0,896; δ b = 10% c = 231,44; δ c = 1% Si nh Vi 12 Khi đo số góc, ta nhận kết sau: a = 450; b = 75020’44” Hãy xác định sai số tương đối số gần đó, sai số tuyệt đối phép đo 1” 13 Xác định số chữ số đáng tin số gần sau biết sai số tuyệt đối chúng: a = 0,1132; ∆ a = 0,1.10−3 b = 2,325; ∆b = 0,1.10−1 c = 293,481; ∆ c = 0,1 14 Hãy xác định số chữ số đáng tin số gần sau biết sai số tương đối chúng là: a = 0,2218; δ a = 0,2.10−1 b = 0,02425; δ b = 0,5.10−2 c = 0,000135; δ c = 0,15 15 Quy tròn số gần với chữ số có nghĩa đáng tin xác định sai số tuyệt đối, sai số tương đối chúng: a) 1,255 b) -392,85 SinhVienZone.com https://fb.com/sinhvienzonevn Mơn: Phương pháp tính c) 0,1545 d) 625,55 16 Đường kính đường tròn đo xác tới 1mm d = 0,842m Tìm diện tích hình tròn 17 Tìm giá trị hàm u = xy2z3 nếu: x = 37,1 ∆ x = 0,3 y = 9,87 ∆ y = 0,11 18 Hãy xác định sai số tuyệt đối số xấp xỉ sau đây, cho biết sai số tương đối nó: b = 12627; δ b = 0, 2% Si nh Vi en Zo ne C om 19 Tính sai số tuyệt đối giới hạn sai số tương đối giới hạn thể tích hình cầu: V = π d3 cho đường kính d = 3,5 ± 0,03cm π = 3,14 ± 0,0016 SinhVienZone.com https://fb.com/sinhvienzonevn Môn: Phương pháp tính CHƯƠNG.2 GIẢI GẦN ĐÚNG PHƯƠNG TRÌNH 2.1 GIỚI THIỆU Để tìm nghiệm gần phương trình f(x) = cần tiến hành qua bước: - Tách nghiệm: xét tính chất nghiệm phương trình, phương trình có nghiệm hay khơng, có nghiệm, khoảng chứa nghiệm có Đối với bước này, ta dùng phương pháp đồ thị, kết hợp với định lý mà tốn học hỗ trợ - Chính xác hoá nghiệm: thu hẹp dần khoảng chứa nghiệm để hội tụ đến giá trị nghiệm gần với độ xác cho phép Trong bước ta áp dụng phương pháp: C om + Phương pháp chia đôi + Phương pháp lặp + Phương pháp tiếp tuyến + Phương pháp dây cung ne 2.2 TÁCH NGHIỆM Zo * Phương pháp đồ thị: Trường hợp hàm f(x) đơn giản en - Vẽ đồ thị f(x) Vi - Nghiệm phương trình hồnh độ giao điểm f(x) với trục x, từ suy số nghiệm, khoảng nghiệm Trường hợp f(x) phức tạp nh - Biến đổi tương đương f(x)=0 g(x) = h(x) - Vẽ đồ thị g(x), h(x) Si - Hoành độ giao điểm g(x) h(x) nghiệm phương trình, từ suy số nghiệm, khoảng nghiệm * Định lý 1: Giả sử f(x) liên tục (a,b) có f(a)*f(b) < Khi (a,b) tồn số lẻ nghiệm thực x∈(a,b) phương trình f(x) = Nghiệm f’(x) tồn khơng đổi dấu (a,b) Ví dụ 2.1: Tách nghiệm cho phương trình : x3 - x + = Giải : f(x) = x3 - x + = f’(x) = 3x2 - 1, f '( x) = x = ±1/ Bảng biến thiên : SinhVienZone.com https://fb.com/sinhvienzonevn Mơn: Phương pháp tính x −1 / -∞ f’(x) + +∞ 1/ - + yCĐ > f(x) +∞ -∞ CT Từ bảng biến thiên, phương trình có nghiệm x < −1 / f(-1)*f(-2) < 0, phương trình có nghiệm x ∈ (−2, −1) Ví dụ 2.2: Giải : 2x + x - = ⇔ 2x = -x + nh Vi en Zo ne Áp dụng phương pháp đồ thị : C om Tách nghiệm cho phương trình : 2x + x - = Từ đồ thị suy phương trình có nghiệm x ∈ (1, 2) Si * Định lý 2: Giả sử α nghiệm x nghiệm gần phương trình f(x) = 0, nằm khoảng nghiệm [a, b] f’(x) ≥ m ≥ a ≤ x ≤ b Khi f ( x) x −α ≤ m Ví dụ 2.3: Cho nghiệm gần phương trình x4 - x - = 1,22 Hãy ước lượng sai số tuyệt đối ? Giải : f(x) = f(1,22) = 1,224 - 1,22 -1 = -0,0047 < f(1,23) = 0,588 > ⇒ nghiệm phương trình x ∈ (1, 22;1, 23) SinhVienZone.com https://fb.com/sinhvienzonevn Mơn: Phương pháp tính f’(x) = 4x3 -1 ≥ 4*1,223 - = 6,264 = m ∀x ∈ (1, 22;1, 23) Theo định lý : ∆x = 0,0047/6,264 = 0,0008 (vì x − α ≤ 0,008 ) 2.3 TÁCH NGHIỆM CHO PHƯƠNG TRÌNH ĐẠI SỐ Xét phương trình đại số : f ( x) = a0 x n + a1x n −1 + + an −1 x + an = (1) * Định lý 3: Cho phương trình (1) có m1 = max { } i = 1, n m2 = max { } i = 0, n − x1 = an m ≤ x ≤ + = x2 m2 + an a0 * Định lý 4: C om Khi nghiệm x phương trình thỏa mãn : ne Cho phương trình (1) có a0 > 0, am hệ số âm Khi nghiệm dương phương trình ≤ N = + m a / a0 , với a = max { } i = 0, n cho < Zo Ví dụ 2.4: en Cho phương trình : 5x5 - 8x3 + 2x2 - x + = Tìm cận nghiệm dương phương trình Vi Giải : nh Ta có a2 = -8 hệ số âm đầu tiên, nên m = 2, a = max(8,1) = Vậy cận nghiệm dương : N = + / Si * Định lý 5: Cho phương trình (1), xét đa thức : ϕ1 ( x) = x n f (1 / x) = a0 + a1 x + + an x n ϕ2 ( x) = f (− x) = (−1) n ( a0 x n − a1x n−1 + a2 x n−2 − + (−1)n an ) ϕ3 ( x) = x n f (−1 / x) = (−1)n (a0 x n − an −1 x n−1 + an−2 x n−2 − + (−1) n a0 ) Giả sử N0, N1, N2, N3 cận nghiệm dương đa thức f(x), φ1(x), φ2(x), φ3(x) Khi nghiệm dương phương trình (1) nằm khoảng [1/N1, N0] nghiệm âm nằm khoảng [-N2,-1/ N3] Ví dụ 2.5: Xét phương trình : 3x2 + 2x - = → N = + / (định lý 4) SinhVienZone.com https://fb.com/sinhvienzonevn Mơn: Phương pháp tính Ta đưa dạng: Y = A + BX Giải hệ phương trình ta A, B => a = 10A, b = B Ví dụ 4.5: Cho biết cặp giá trị x y theo bảng sau: xi 0,65 0,75 0,85 0,95 1,15 yi 0,96 1,06 1,17 1,29 1,58 Lập công thức thực nghiệm y dạng aebx Giải Ta có: y = aebx Đặt Y = lny; A = lna; B = b; X = x 0,65 0,75 0,85 0,95 1,15 Yi = lnyi -0,04 0,06 0,18 0,25 0,46 ƩXi ƩXi2 ƩXiYi ƩYi 4,35 0,92 0,89 ne X i = xi Zo Ta đưa dạng: Y = A + BX C om Lấy Logarit số e hai vế: lny = lna + bx en 3,93 Phương pháp bình phương bé nhất: A, B nghiệm hệ phương trình: Si nh Vi n n  nA + B ∑ X i = ∑ Yi  i =1 i =1  n n n  A X + B X = ∑ i ∑ X iYi  ∑ i i =1 i =1  i =1 5 A + 4.35 B = 0.89   4.35 A + 3.93B = 0.92 Giải hệ phương trình ta được: A = -0.69, B = Suy ra: a = eA = 1/2, b = B =1 Vậy f(x) = 1/2 ex 43 SinhVienZone.com https://fb.com/sinhvienzonevn Mơn: Phương pháp tính BÀI TẬP CUỐI CHƯƠNG Nêu công thức nội suy Lagrange x 1,50 1,54 1,56 y 3,873 3,924 3,950 C om Nêu công thức nội suy Lagrange với mối cách Nêu công thức nội suy Newton Nêu công thức nội suy tổng quát Trình bày phương pháp bình phương bé với trường hợp f hàm tuyến tính bậc f(x) = ax + b Trình bày phương pháp bình phương bé với trường hợp f hàm tuyến tính bậc hai f(x) = ax2 + bx + c Trình bày phương pháp bình phương bé với trường hợp f hàm phi tuyến tính dạng f(x) = aebx Trình bày phương pháp bình phương bé với trường hợp f hàm phi tuyến tính dạng f(x) = axb Cho bảng giá trị hàm 1,60 1,63 1,70 4,000 4,037 4,123 ne 1,1 1,2 nh y 1,0 0,5652 0,6375 1,3 1,4 1,5 1,6 1,7 1,8 1,9 2,0 0,7147 0,7973 0,8861 0,9817 1,0848 1,1964 1,3172 1,4482 Si x Vi en Zo Sử dụng công thức nội suy Lagrăng tìm giá trị hàm điểm: a) 1,52 b) 1,55 c) 1,58 d) 1,61 10 Cho bảng giá trị hàm Sử dụng công thức nội suy Newton xác định giá trị hàm điểm: a) 1,0113 b) 1,0428 c) 1,9592 d) 1,9728 11 Cho giá trị hai đại lượng x, y bảng sau: i xi 0,56 0,84 1,14 2,44 3,16 yi -0,80 -0.97 -0,98 1,07 3,66 Tìm xấp xỉ hàm dạng bậc 2: y = ax2 + bx + c 44 SinhVienZone.com https://fb.com/sinhvienzonevn 1,5906 Mơn: Phương pháp tính 12 Cho bảng giá trị hàm x 19 22 25 28 32 35 y 0,66 0,367 0,223 0,14 0,084 0,06 C om Tìm hàm xấp xỉ phương pháp bình phương bé sau đánh giá sai số hàm xấp xỉ nếu: 13 Quan hệ y x tuyến tính: y = ax + b; 14 Quan hệ y x tam thức bậc hai: y = ax + bx + c; bx 15 Quan hệ y x hàm mũ: y = ae 16 Cho bảng số liệu: x y 16 26,5 211,5 x -1 0,5 1,5 2,5 en y Zo ne Hãy lập đa thức nội suy Lagrăng tương ứng 17 Hai đại lượng x y phụ thuộc theo quy luật y = ax + b Hãy xác định a, b phương pháp bình phương bé nhất, biết: x y -2 Si nh Vi 18 Cho hàm số f(x) thỏa mãn: Xây dựng hàm nội suy Newton 19 Cho hàm số f(x) thỏa mãn: x -1 y -2 Xây dựng hàm nội suy Newton 20 Cho hàm số f(x) thỏa mãn: x y 17 27,5 76 210,5 45 SinhVienZone.com https://fb.com/sinhvienzonevn Mơn: Phương pháp tính Xây dựng hàm nội suy Newton 21 Cho hàm số f(x) thỏa mãn: x y -1 Xây dựng hàm nội suy Lagrăng f(x), tính f(3,5) 22 Cho hàm số f(x) thỏa mãn: y 17,0 27,5 76 C om x 210,5 1970 Si nh Vi en Zo ne Xây dựng hàm nội suy Lagrăng f(x), tính f(5) 46 SinhVienZone.com https://fb.com/sinhvienzonevn Mơn: Phương pháp tính CHƯƠNG.5 TÍNH GẦN ĐÚNG ĐẠO HÀM VÀ TÍCH PHÂN XÁC ĐỊNH 5.1 TÍNH GẦN ĐÚNG ĐẠO HÀM Người ta thường dùng số phương pháp để tính gần đạo hàm hàm f(x) x phương pháp áp dụng đa thức nội suy thường dùng Giả sử người ta phải tính xấp xỉ đạo hàm hàm số f(x) đoạn (a,b) Trước hết người ta thay hàm f(x) đa thức nội suy P(x), sau lấy đạo hàm P'(x) coi xấp xỉ đạo hàm f'(x) Ví dụ 5.1: Giả sử ta xác định đa thức nội suy là: C om P3(x) = 8x3 - 29x + Khi đó, đạo hàm P3’(x) = 24x2 - 29 xem xấp xỉ f’(x) 5.2 TÍNH GẦN ĐÚNG TÍCH PHÂN XÁC ĐỊNH ne Xét hàm số f(x) liên tục [a,b], xác định ngun hàm F(x) ta có cơng thức tính tích phân: b Zo ∫a f ( x)dx = F (b) − F (a) Vi en Nhưng đa số trường hợp ta không xác định nguyên hàm không xác định biểu thức f(x) mà nhận giá trị điểm rời rạc Trong trường hợp ta sử dụng cơng thức gần sau để tính tích phân: - Cơng thức hình thang nh - Công thức Parabol - Công thức Newton - Cotet Si 5.2.1 Cơng thức hình thang Chia [a, b] thành n đoạn với khoảng cách h = (b - a)/n theo điểm chia: x0 = a, x1 = a + h, , xn = b b x1 x2 ∫a f ( x)dx = ∫x =a f ( x)dx + ∫x f ( x)dx + + ∫ xn xn −1 f ( x)dx = S S diện tích giới hạn đường cong f(x), x = a, x = b, trục x 47 SinhVienZone.com https://fb.com/sinhvienzonevn Mơn: Phương pháp tính Xét [x0, x1], ta xem đường cong f(x) đường thẳng: Tương tự: ne S2 ≈ h ( y1 + y2 ) C om S1 ≈ Sh.thang = h ( y0 + y1 ) … Zo Sn ≈ h ( yn−1 + yn ) b f ( x)dx ≈ h ( y0 + y1 + y2 + + yn−1 + yn ) Vi ∫a en Vậy: 5.2.2 Công thức Parabol Si nh Chia [a, b] thành 2n đoạn với khoảng cách h = (b - a)/2n theo điểm chia: x0 = a, x1 = a + h, , x2n = b b x2 ∫a f ( x)dx = ∫x f ( x)dx + ∫ x4 x2 f ( x)dx + + ∫ x2 n x2 n − f ( x)dx Xét [x0, x2] xem đường cong f(x) Parabol (nội suy bậc điểm x0, x1, x2): f ( x) ≈ L2 ( x) = y0 x2 ∫x ( x − x1 )( x − x2 ) + y ( x − x0 ) ( x − x2 ) + y ( x − x0 ) ( x − x1 ) ( x0 − x1 )( x0 − x2 ) ( x1 − x0 ) ( x1 − x2 ) ( x2 − x0 ) ( x2 − x1 ) x2 f ( x)dx ≈ ∫ L2 ( x)dx x0 Thay x0 = a, x1 = a + h, x2 = a + 2h vào ta có: x2 ∫x f ( x)dx ≈ h ( y0 + y1 + y2 ) 48 SinhVienZone.com https://fb.com/sinhvienzonevn Môn: Phương pháp tính Tương tự: x4 ∫x h ( y2 + y3 + y4 ) f ( x)dx ≈ x2 n ∫x n− f ( x)dx ≈ h ( y2n−2 + y2n−1 + y2n ) Vậy: b ∫a f ( x)dx ≈ h ( y0 + y1 + y2 + + y2n−2 + y2n−1 + y2n ) Ví dụ 5.2: dx 1 + x2 I =∫ Giải: dx 1+ x = arctan x = arctan − π ne Cách 1: I = ∫ C om Tính tích phân sau theo cách: ≈ 0.588 yi 1/2 1/5 1/10 1/17 1/26 en xi Zo Cách 2: Chia [1,5] thành đoạn thẳng (h = 1) với điểm chia: Theo cơng thức hình thang: Vi I ≈ (1/2 + 2/5 +2/10 + 2/17 +1/26)/2 ≈ 0.628 nh Cách 3: Công thức Parabol: I ≈ (1/2 + 4/5 +2/10 + 4/17 +1/26)/3 ≈ 0.591 Si 5.2.3 Công thức Newton-Cotet Chia [a, b] thành n đoạn với khoảng cách h = (b - a)/n theo điểm chia: x0 = a, x1 = a + h, , xn = b Đặt x = a + (b - a)t => dx = (b - a) dt xi a a+h a + 2h … b yi 1/n 2/n … Khi đó: b ∫a f ( x)dx = ( b − a ) ∫ f (a + (b − a )t )dt = (b − a) ∫ φ (t )dt 1 0 với Ф(t) = f(a + (b - a)t) Xem Ф(t) hàm nội suy Lagrange n + điểm: t0, t1, , tn 49 SinhVienZone.com https://fb.com/sinhvienzonevn Môn: Phương pháp tính    ( t − )  t −  ( t − 1)  t −  t −  ( t − 1) n  n   n φ (t ) ≈ Ln (t ) = y0  + y1 +    1      −  −  (−1)  −  −   − 1  n  n  n  n n   n  n −1 ( t − )  t −   t −  n   n  + yn n −1 (1 − ) 1 −  1 −  n   n  ∫0 Khi đó: φ (t )dt ≈ ∫ Ln (t )dt Đặt C om   i −  i +    t −  t −  ( t − 1) n  n  n  i  Pn = ∫ dt 0 i  i   i i −  i i +   i   −  −   −  −   − 1 n  n n  n  n  n n   n ∫a i =0 Zo Vậy; n f ( x)dx ≈ ( b − a ) ∑ yi pni b Xét n = (h = b - a) −1 dt = − 0 −1 1t en P10 = ∫ ne ( t − )  t − −0 dt = 01− y  h y f ( x )dx = ( b − a )  +  = ( y0 − y1 ) (Cơng thức hình thang) 2  nh b ∫a 1t Vi P11 = ∫ Si Giá trị Pni tra bảng sau: n Pni 1/2 1/2 1/6 4/6 1/6 1/8 3/8 3/8 1/8 9/71 16/45 2/15 16/45 … 9/70 19/288 25/95 25/144 25/144 25/95 … … … … 19/288 … … 50 SinhVienZone.com https://fb.com/sinhvienzonevn Mơn: Phương pháp tính xi 0,4 0,6 0,8 1,0 1,2 1,4 yi = f(xi) 0,4000 1,4848 2,6813 3,9975 5,3456 6,2465 C om BÀI TẬP CUỐI CHƯƠNG Nêu cơng thức tính gần đạo hàm trường hợp toán mốc cách Trình bày cơng thức hình thang tính gần tích phân xác định Trình bày cơng thức Parabol tính gần tích phân xác định Trình bày cơng thức Newton-Cotet tính gần tích phân xác định Tính giá trị đạo hàm cấp 1, cấp 2, giá trị hàm cho bảng sau: Tính gần y’(50) hàm số y = logx dựa vào bảng giá trị cho sau: 50 yi = log(xi) 1,6990 60 1,7404 1,7782 Zo en Cho hàm f(x) bảng sau: Vi xi 50 55 60 65 1,6990 1,7404 1,7782 1,8129 nh yi = log(xi) 55 ne xi Áp dụng đa thức nội suy tính gần đạo hàm hàm f(x) x = 50 sso Si sánh với kết tính trực tiếp Cho hàm f(x) bảng sau: xi 0,98 1,00 1,02 yi = f(xi) 0,7739332 0,7651977 0,7563321 Tính gần đạo hàm hàm f(x) x = Tính giá trị đạo hàm cấp cấp 2, giá trị hàm cho bảng sau: 51 SinhVienZone.com https://fb.com/sinhvienzonevn Mơn: Phương pháp tính xi 1,0 1,1 1,2 1,3 1,4 1,5 yi = f(xi) 1,266 1,326 1,393 1,469 1,553 1,647 10 Cho hàm y = f(x) dạng bảng sau: x y 12,3 11,1 7,2 4,1 6,3 8,8 9,2 10,8 13,1 Tính tích phân: I = ∫ f ( x )dx theo cơng thức hình thang .C om x t2 − 11 Trong kĩ thuật ta thường gặp tích phân xác suất: ϕ ( x) = e dt ∫ 2π ne Hãy tính φ(1) theo cơng thức hình thang chia khoảng tích phân thành 10 phần Cho bảng giá trị hàm dấu tích phân y = e y 0,1 0,2 0,3 0,4 0,9950 0,9802 0,9560 0,9231 0,5 0,6 0,7 0,8 0,9 1,0 0,8822 0,8353 0,7827 0,7261 0,6670 0,6065 Zo t2 en x − Vi 12 nh Cho tích phân I =∫ 1+ x dx Bằng cách phân hoạch đoạn [0, 1] thành đoạn nhau, tính gần tích Si phân theo cơng thức hình thang 13 Cho tích phân sin x dx x I =∫ a) Bằng cách phân hoạch [0, 1] thành đoạn Tính gần tích phân cho cơng thức hình thang công thức Simpson Đánh giá sai số? b) Tính gần tích phân cơng thức hình thang với sai số không 3.10-4 14 Cho hàm f(x) dạng bảng sau: 52 SinhVienZone.com https://fb.com/sinhvienzonevn Môn: Phương pháp tính x y 0,2 0,4 1,0000 0,9801 0,6 0,9211 0,8 0,8253 0,6967 Tính tích phân hàm sau theo cơng thức hình thang 0,8 ∫ I= f ( x)dx 15 Cho hàm f(x) dạng bảng sau: 0,0 0,5 1,0 1,5 y 1,50 0,75 0,50 0,75 2,0 2,5 3,0 3,5 4,0 1,50 2,75 4,50 6,75 10,00 C om x Tính tích phân hàm sau theo cơng thức hình thang ne I = ∫ f ( x )dx Zo 0,0 0,1 y 1,000 0,990 0,2 0,3 0,4 0,5 0,6 0,962 0,917 0,862 0,800 0,735 Vi x en 16 Cho hàm f(x) dạng bảng sau: 0,7 0,8 0,9 0,671 0,609 1,0 0,555 0,500 I = ∫ f ( x )dx Si nh Tính tích phân hàm sau theo cơng thức hình thang 17 Cho hàm f(x) dạng bảng sau: x 0,0 0,1 0,2 0,3 0,4 0,5 0,6 0,7 0,8 0,9 1,0 y 1/1,1 1/1,2 1/1,3 1/1,4 1/1,5 1/1,6 1/1,7 1/1,8 1/1,9 1/2 Tính tích phân hàm sau theo cơng thức hình thang I = ∫ f ( x )dx 18 Cho tích phân 53 SinhVienZone.com https://fb.com/sinhvienzonevn Mơn: Phương pháp tính sin x dx x I =∫ Hãy phân hoạch đoạn [0, 1] thành 10 đoạn tính gần tích Si nh Vi en Zo ne C om phân cho công thức hình thang 54 SinhVienZone.com https://fb.com/sinhvienzonevn Mơn: Phương pháp tính TÀI LIỆU THAM KHẢO [1] Ralston A, A first course in numberical analysis McGraw – Hill, NewYork, 2001, 576 pages [2] Đỗ Thị Tuyết Hoa, Giáo trình mơn Phương pháp tính, Đại học Đà Nẵng, 2007, 68 trang [3] Phan Văn Hạp, Hồng Đức Ngun, Lê Đình Thịnh, Phương pháp tính, NXB Khoa học Kỹ thuật, 1996 .C om [4] Phan Văn Hạp, Hồng Đức Ngun, Lê Đình Thịnh, Phương pháp tính (phần tập), NXB Khoa học Kỹ thuật, 1996, 204 trang [5] Tạ Văn Đĩnh, Phương pháp tính, NXB Khoa học Kỹ thuật, 2009, 118 trang [6] Đặng Quốc Lương, Phương pháp tính kỹ thuật, NXB Xây Dựng, 2001, 133 trang ne [7] Dương Thủy Vỹ, Giáo trình Phương pháp tính, NXB Khoa học Kỹ thuật, 2007, 180 trang Zo [8] Trần Văn Chính, Phương pháp tính với C++, NXB Đại học Quốc gia TPHCM, 2008 en [9] Nguyễn Hoài Sơn, Phương pháp tính ứng dụng tính tốn Kỹ thuật, NXB Đại học Quốc gia TPHCM, 2008, 260 trang Si nh Vi [10] Nguyễn Trọng Khiêm, Bài giảng Phương pháp tính, NXB Khoa học Kỹ thuật, 2009 55 SinhVienZone.com https://fb.com/sinhvienzonevn Môn: Phương pháp tính MỤC LỤC Trang CHƯƠNG.1 SAI SỐ .1 1.1 NHẬP MƠN PHƯƠNG PHÁP TÍNH 1.1.1 Giới thiệu môn phương pháp tính 1.1.2 Nhiệm vụ môn học 1.1.3 Trình tự giải tốn phương pháp tính C om 1.2 SAI SỐ .2 1.2.1 Khái niệm 1.2.2 Các loại sai số 1.2.3 Sai số tính tốn CHƯƠNG.2 GIẢI GẦN ĐÚNG PHƯƠNG TRÌNH ne 2.1 GIỚI THIỆU 2.2 TÁCH NGHIỆM Zo 2.3 TÁCH NGHIỆM CHO PHƯƠNG TRÌNH ĐẠI SỐ en 2.4 CHÍNH XÁC HÓA NGHIỆM 10 2.4.1 Phương pháp chia đôi 10 2.4.2 Phương pháp lặp 11 Vi 2.4.3 Phương pháp tiếp tuyến 13 nh 2.4.4 Phương pháp dây cung 15 CHƯƠNG.3 GIẢI HỆ PHƯƠNG TRÌNH ĐẠI SỐ TUYẾN TÍNH 20 Si 3.1 GIỚI THIỆU 20 3.2 PHƯƠNG PHÁP KRAME 20 3.3 PHƯƠNG PHÁP GAUSS .20 3.3.1 Nội dung phương pháp 20 3.3.2 Thuật toán 21 3.4 PHƯƠNG PHÁP LẶP GAUSS - SIEDEL (TỰ SỬA SAI) 22 3.4.1 Nội dung phương pháp 22 3.4.2 Thuật toán 24 3.5 PHƯƠNG PHÁP GIẢM DƯ 25 3.5.1 Nội dung phương pháp 25 3.5.2 Thuật toán 27 56 SinhVienZone.com https://fb.com/sinhvienzonevn Môn: Phương pháp tính CHƯƠNG.4 NỘI SUY VÀ PHƯƠNG PHÁP BÌNH PHƯƠNG BÉ NHẤT 32 4.1 GIỚI THIỆU 32 4.2 ĐA THỨC NỘI SUY LAGRANGE .33 4.2.1 Nội suy bậc (nội suy tuyến tính) 33 4.2.2 Nội suy bậc hai 33 4.2.3 Nội suy bậc ba 34 4.3 ĐA THỨC NỘI SUY LAGRANGE VỚI CÁC MỐI CÁCH ĐỀU 35 4.4 NỘI SUY NEWTON 37 4.4.1 Sai phân 37 C om 4.4.2 Công thức nội suy Newton 37 4.5 NỘI SUY TỔNG QUÁT 39 4.6 PHƯƠNG PHÁP BÌNH PHƯƠNG BÉ NHẤT 41 CHƯƠNG.5 TÍNH GẦN ĐÚNG ĐẠO HÀM VÀ TÍCH PHÂN XÁC ĐỊNH 47 ne 5.1 TÍNH GẦN ĐÚNG ĐẠO HÀM 47 Zo 5.2 TÍNH GẦN ĐÚNG TÍCH PHÂN XÁC ĐỊNH 47 5.2.1 Cơng thức hình thang 47 en 5.2.2 Công thức Parabol 48 Vi 5.2.3 Công thức Newton-Cotet 49 Si nh PHÒNG KHOA HỌC Nguyễn Viết Tuấn GV biên soạn Phạm Thị Ngọc Minh 57 SinhVienZone.com https://fb.com/sinhvienzonevn ... phương pháp - Biến đổi Ma trận A ma trận tam giác 20 SinhVienZone. com https://fb .com/ sinhvienzonevn Mơn: Phương pháp tính  a1 1 a A =  21    an1 a1 2 a2 2 a1 n a2 n an ann  a1 1 0 ⇒ A' ... a x + a x + + a x = a  21 22 2n n n +1    an1 x1 + an x2 + + ann xn = ann +1 Hệ phương trình cho ma trận: a1 2 a2 2 a1 n a2 n an ann a1 n +1  a2 n +1     ann +1  C om  a1 1 a Ann... cho phương trình đại số Có phương pháp xác h a nghiệm? Liệt kê phương pháp đó? Trình bày ý tưởng phương pháp chia đơi để xác h a nghiệm Trình bày thuật tốn phương pháp chia đơi để xác h a nghiệm