Máy tính quy tắc của Simpson

Ước lượng độ dài đường cong bằng Quy tắc Simpson khi tích phân chính xác không khả thi hoặc không cần thiết.

Quy tắc Simpson là gì?

Quy tắc Simpson là một phương pháp tích phân số dùng đa thức bậc hai để xấp xỉ giá trị của một tích phân xác định. Thay vì ước lượng diện tích dưới đường cong bằng hình chữ nhật hay hình thang, phương pháp này khớp một parabol qua mỗi nhóm ba điểm mẫu liên tiếp rồi tính diện tích chính xác dưới mỗi parabol. Tổng các diện tích parabol này cho ra phép xấp xỉ.

Phương pháp mang tên Thomas Simpson, nhà toán học người Anh thế kỷ 18, dù ý tưởng nền tảng đã xuất hiện sớm hơn trong công trình của người khác. Sức mạnh của nó đến từ một quan sát đơn giản: một parabol có thể bám theo đường cong trơn sát hơn nhiều so với đường thẳng. Khi hàm được lấy tích phân là hàm trơn, các mảnh parabol ôm sát đường cong thực, và phép xấp xỉ thu được chính xác đáng kinh ngạc ngay cả với số khoảng con vừa phải.

Để hiểu vì sao ba điểm lại quan trọng, hãy nhớ rằng đường thẳng được xác định hoàn toàn bởi hai điểm, trong khi parabol (đa thức bậc hai) cần ba điểm. Quy tắc Simpson khai thác điều này bằng cách chia khoảng tích phân thành số chẵn các khoảng con, nhóm chúng thành từng cặp, rồi khớp một parabol trên mỗi cặp. Mỗi cặp trải dài hai khoảng con và ba điểm: đầu mút trái, trung điểm và đầu mút phải.

Dạng cổ điển được gọi là Quy tắc Simpson 1/3 vì hệ số 1/3 xuất hiện trong công thức. Cũng có Quy tắc Simpson 3/8 dùng đa thức bậc ba trên các nhóm ba khoảng con, nhưng Quy tắc 1/3 là phương pháp chuẩn và là phiên bản được cài đặt trong hầu hết các máy tính, kể cả máy tính trên trang này.

Quy tắc Simpson thuộc họ công thức Newton-Cotes, xấp xỉ tích phân bằng cách tính hàm dưới dấu tích phân tại các điểm cách đều rồi kết hợp các giá trị đó với trọng số cố định. Điều làm Quy tắc Simpson khác biệt là mẫu trọng số: 1, 4, 2, 4, 2, ..., 4, 1. Các đầu mút nhận trọng số 1, các điểm trong ở vị trí lẻ nhận trọng số 4, và các điểm trong ở vị trí chẵn nhận trọng số 2. Mẫu xen kẽ này là dấu ấn của phương pháp và là nguồn gốc độ chính xác của nó.

Hiểu về độ dài cung trong toán học

Độ dài của một đoạn thẳng rất dễ tính: đo khoảng cách giữa hai đầu mút. Đường cong thì khác. Nó uốn lượn, nên không một phép đo thẳng đơn lẻ nào nắm bắt được độ trải dài thực sự của nó. Để tìm độ dài của một đường cong, các nhà toán học chia nó thành những mảnh nhỏ, xấp xỉ mỗi mảnh như một đoạn thẳng, rồi cộng các kết quả lại.

Hãy tưởng tượng bạn đi bộ dọc một con đường quanh co. Nếu bạn bước những bước rất nhỏ, mỗi bước gần như thẳng, dù toàn bộ con đường thì cong. Cộng tất cả những bước thẳng nhỏ ấy lại, bạn có được ước lượng tốt về tổng quãng đường đã đi. Làm cho các bước ngày càng nhỏ, ước lượng sẽ hội tụ về độ dài chính xác. Đó là ý tưởng cốt lõi đằng sau độ dài cung.

Với hàm y = f(x) trên khoảng [a, b], xét một bước ngang nhỏ dx. Trên bước đó, hàm tăng thêm khoảng dy = f'(x) dx, trong đó f'(x) là đạo hàm (độ dốc). Đoạn cong nhỏ có độ dài cho bởi định lý Pythagore: ds = √((dx)2 + (dy)2) = √(1 + (f'(x))2) dx. Cộng tất cả những mảnh vô cùng nhỏ này qua phép tích phân cho ra công thức độ dài cung:

L = ∫ab √(1 + (f'(x))2) dx

Mỗi phần của công thức này có một ý nghĩa rõ ràng. Ký hiệu ∫ab cho biết ta tích lũy các đóng góp từ x = a đến x = b. Biểu thức f'(x) là đạo hàm, đo độ dốc của đường cong tại mỗi điểm. Bình phương nó, cộng 1 rồi lấy căn bậc hai sẽ chuyển độ dốc thành hệ số giãn cục bộ: đường cong dài hơn đường nằm ngang phẳng tại vị trí đó bao nhiêu. Nơi đường cong phẳng, f'(x) = 0 và hệ số bằng 1. Nơi đường cong dốc, hệ số tăng lên.

Công thức này đòi hỏi hàm phải đủ trơn: f phải có đạo hàm liên tục trên khoảng, và hàm dưới dấu tích phân √(1 + (f'(x))2) phải khả tích. Với hầu hết các hàm gặp trong các khóa học giải tích và công việc kỹ thuật, các điều kiện này đều thỏa mãn. Vấn đề nảy sinh ở các góc nhọn, tiếp tuyến đứng hoặc điểm gián đoạn, nơi đạo hàm hoạt động bất thường.

Khó khăn là tích phân này thường không có nguyên hàm sơ cấp. Ngay cả với các hàm đơn giản như y = x2, nguyên hàm đã liên quan đến các hàm hypebolic ngược. Với nhiều đường cong thực tế, hoàn toàn không tồn tại nguyên hàm dạng đóng. Đây chính là lúc các phương pháp số như Quy tắc Simpson trở nên thiết yếu: chúng ước lượng trực tiếp tích phân mà không cần nguyên hàm.

Quy tắc Simpson tính độ dài cung như thế nào

Áp dụng Quy tắc Simpson cho độ dài cung gồm một chuỗi các bước rõ ràng. Ý tưởng then chốt là Quy tắc Simpson không được áp dụng cho hàm gốc f(x). Nó được áp dụng cho hàm dưới dấu tích phân của độ dài cung g(x) = √(1 + (f'(x))2). Sự phân biệt này rất quan trọng và là nguồn gốc phổ biến của các sai lầm.

Đây là toàn bộ quy trình:

  1. Bắt đầu với hàm số. Xác định y = f(x) và xác nhận nó khả vi trên khoảng của bạn.
  2. Tính đạo hàm. Tìm f'(x) bằng giải tích hoặc xác minh nó đã có sẵn.
  3. Xây dựng hàm dưới dấu tích phân. Lập g(x) = √(1 + (f'(x))2). Đây là hàm mà Quy tắc Simpson sẽ lấy tích phân.
  4. Xác định khoảng. Xác nhận cận dưới a và cận trên b của đoạn cong bạn muốn đo.
  5. Chọn số khoảng con. Chọn số chẵn n khoảng con. Số này phải chẵn vì Quy tắc Simpson làm việc theo cặp.
  6. Tính độ dài bước. Tính h = (b - a) / n, độ rộng của mỗi khoảng con.
  7. Tìm các điểm mẫu. Các điểm là x0 = a, x1 = a + h, x2 = a + 2h, cứ thế đến xn = b.
  8. Tính giá trị hàm dưới dấu tích phân. Tính g(x0), g(x1), ..., g(xn) tại mọi điểm mẫu.
  9. Áp trọng số. Nhân với mẫu Simpson: 1 cho các đầu mút, 4 cho các điểm trong có chỉ số lẻ, 2 cho các điểm trong có chỉ số chẵn.
  10. Kết hợp và co giãn. Cộng các giá trị đã nhân trọng số rồi nhân với h/3 để được ước lượng độ dài cung.

Công thức Quy tắc Simpson 1/3 tổng hợp cho độ dài cung là:

L ≈ (h/3)[g(x0) + 4g(x1) + 2g(x2) + 4g(x3) + ... + 4g(xn-1) + g(xn)]

trong đó h = (b - a)/n và g(x) = √(1 + (f'(x))2).

Yêu cầu n phải chẵn xuất phát từ cấu trúc cặp. Mỗi parabol trải dài hai khoảng con (ba điểm). Nếu n lẻ, khoảng con cuối cùng sẽ không có cặp, và phương pháp sẽ sụp đổ. Khi vì lý do nào đó bạn cần số panel lẻ, người ta dùng quy tắc hình thang hoặc một cách tiếp cận hỗn hợp thay thế.

Hãy chú ý kỹ mẫu trọng số. Giá trị đầu tiên và cuối cùng (g(x0) và g(xn)) mỗi giá trị xuất hiện một lần với trọng số 1. Các điểm trong có chỉ số lẻ (x1, x3, x5, ...) nhận trọng số 4 vì mỗi điểm đóng vai trò trung điểm của panel parabol của nó, nơi ảnh hưởng của parabol mạnh nhất. Các điểm trong có chỉ số chẵn (x2, x4, ...) nhận trọng số 2 vì mỗi điểm đóng vai trò đầu mút chung giữa hai panel kề nhau.

Giải thích công thức Quy tắc Simpson từng bước

Công thức Quy tắc Simpson tổng hợp có nhiều thành phần, mỗi thành phần đóng một vai trò cụ thể. Hiểu từng thành phần riêng lẻ khiến phương pháp bớt đáng sợ hơn nhiều.

Độ dài bước (h). Đây là độ rộng của mỗi khoảng con, tính bằng h = (b - a)/n. h càng nhỏ nghĩa là càng nhiều panel và thường độ chính xác càng tốt, với cái giá là nhiều lần tính hàm hơn. Với độ dài cung, h kiểm soát đường cong được lấy mẫu mịn đến mức nào.

Cận tích phân (a và b). Chúng xác định đoạn cong được đo. Việc chọn khoảng ảnh hưởng lớn đến kết quả: kéo dài b xa hơn dọc theo đường cong sẽ cộng thêm độ dài. Luôn kiểm tra kỹ rằng a < b và cả hai đều nằm trong miền xác định của hàm.

Điểm mẫu. n+1 điểm từ x0 đến xn cách đều nhau trên [a, b]. Quy tắc Simpson tính hàm dưới dấu tích phân tại từng điểm này. Không điểm nào bị bỏ qua.

Đóng góp của đầu mút. Các giá trị g(x0) và g(xn) mỗi giá trị tham gia tổng một lần, với trọng số 1. Chúng neo phép xấp xỉ tại các biên.

Đóng góp chỉ số lẻ. Các điểm x1, x3, x5, v.v. nhận trọng số 4. Đây là các trung điểm của mỗi panel parabol. Trọng số 4 phản ánh thực tế rằng trung điểm có ảnh hưởng lớn nhất đến diện tích dưới một parabol.

Đóng góp trong chỉ số chẵn. Các điểm x2, x4, x6, v.v. (không tính các đầu mút) nhận trọng số 2. Mỗi điểm này được chia sẻ làm đầu mút giữa hai panel kề nhau, nên tổng đóng góp của nó được nhân đôi.

Tổng có trọng số. Nhân mỗi giá trị hàm với trọng số của nó, cộng tất cả lại, rồi nhân với h/3. Hệ số h/3 đưa tổng có trọng số về đúng đơn vị và độ lớn.

Bảng sau cho thấy mẫu trọng số với n = 6:

ĐiểmGiá trị xTrọng sốVai trò
x0a1Đầu mút trái
x1a + h4Trung điểm panel 1
x2a + 2h2Đầu mút chung
x3a + 3h4Trung điểm panel 2
x4a + 4h2Đầu mút chung
x5a + 5h4Trung điểm panel 3
x6b1Đầu mút phải

Cách dùng máy tính Quy tắc Simpson

Máy tính trên trang này ước lượng độ dài cung bằng Quy tắc Simpson qua nhiều chế độ nhập liệu. Mỗi chế độ xử lý một cách mô tả đường cong khác nhau. Hiểu chế độ nào phù hợp với bài toán của bạn là bước đầu tiên.

Chế độ hàm số. Nhập một công thức như x^2, sin(x) hoặc exp(x), cùng với các cận a và b của khoảng. Máy tính lấy đạo hàm bằng số hoặc dùng cấu trúc bạn cung cấp, xây dựng hàm dưới dấu tích phân của độ dài cung, rồi áp dụng Quy tắc Simpson. Đây là chế độ phổ biến nhất cho các bài toán giải tích.

Chế độ bán kính và góc. Với các cung tròn, nhập bán kính và góc ở tâm. Độ dài cung suy ra trực tiếp từ hình học, và Quy tắc Simpson không cần thiết cho riêng trường hợp tròn, nhưng chế độ này hữu ích để so sánh.

Chế độ tham số. Nhập x(t) và y(t) dưới dạng hàm của tham số t, cùng với khoảng tham số. Máy tính dùng công thức độ dài cung tham số với Quy tắc Simpson áp dụng cho √((dx/dt)2 + (dy/dt)2).

Chế độ cực. Nhập r(θ) và khoảng góc. Quy tắc Simpson lấy tích phân √(r2 + (dr/dθ)2) trên miền góc.

Chế độ điểm. Nhập một danh sách các cặp tọa độ. Máy tính cộng các khoảng cách đường thẳng giữa các điểm liên tiếp. Đây là phép xấp xỉ đa giác chứ không phải Quy tắc Simpson, hữu ích khi bạn có dữ liệu đo đạc thay vì công thức.

Chế độ tham số 3D. Nhập x(t), y(t) và z(t) cho một đường cong không gian. Quy tắc Simpson lấy tích phân √((dx/dt)2 + (dy/dt)2 + (dz/dt)2).

Để chạy một phép tính độ dài cung bằng Quy tắc Simpson ở chế độ hàm số:

  1. Nhập hàm số. Gõ công thức bằng ký hiệu chuẩn, ví dụ x^2 hoặc sin(x).
  2. Đặt cận dưới. Nhập giá trị x bắt đầu a cho đoạn cong.
  3. Đặt cận trên. Nhập giá trị x kết thúc b. Xác nhận rằng a < b.
  4. Chọn số khoảng con. Chọn một số chẵn n. Bắt đầu với giá trị vừa phải như 20 hoặc 50.
  5. Xác nhận phương pháp. Trang này khóa phương pháp số ở Quy tắc Simpson, nên không cần chọn.
  6. Chạy phép tính. Công cụ tính hàm dưới dấu tích phân tại n+1 điểm và trả về ước lượng.
  7. Kiểm tra hội tụ. Chạy lại với n lớn hơn (ví dụ gấp đôi) và so sánh. Nếu kết quả hầu như không đổi, nó đáng tin cậy.

Thử một ví dụ thực tế: f(x) = x2 trên [0, 1] với n = 20, rồi n = 40. Hai kết quả nên khớp nhau đến vài chữ số thập phân, xác nhận ước lượng ổn định.

Ví dụ minh họa: Tính độ dài cung của y = x<sup>2</sup> bằng Quy tắc Simpson

Phần này trình bày một ví dụ đầy đủ với mọi bước trung gian được thể hiện. Hàm số là y = x2 trên khoảng [0, 1], dùng n = 4 khoảng con. Mọi giá trị đều được tính độc lập và có thể kiểm chứng bằng bất kỳ máy tính nào.

Bước 1: Hàm số. Ta bắt đầu với y = x2, một parabol đơn giản.

Bước 2: Đạo hàm. Lấy đạo hàm cho y' = 2x.

Bước 3: Hàm dưới dấu tích phân của độ dài cung. Thay vào công thức: g(x) = √(1 + (2x)2) = √(1 + 4x2).

Bước 4: Tích phân. Độ dài cung chính xác là L = ∫01 √(1 + 4x2) dx. Tích phân này có dạng đóng liên quan đến các hàm hypebolic ngược, và giá trị của nó xấp xỉ 1.4789428575. Ta sẽ xem Quy tắc Simpson tiến gần đến mức nào.

Bước 5: Độ dài bước. Với n = 4 trên [0, 1], h = (1 - 0)/4 = 0.25.

Bước 6: Điểm mẫu. Năm điểm cách đều là:

  • x0 = 0.00
  • x1 = 0.25
  • x2 = 0.50
  • x3 = 0.75
  • x4 = 1.00

Bước 7: Giá trị hàm. Tính g(x) = √(1 + 4x2) tại mỗi điểm:

  • g(0.00) = √(1 + 0) = 1.0000000000
  • g(0.25) = √(1 + 0.25) = √1.25 ≈ 1.1180339887
  • g(0.50) = √(1 + 1) = √2 ≈ 1.4142135624
  • g(0.75) = √(1 + 2.25) = √3.25 ≈ 1.8027756377
  • g(1.00) = √(1 + 4) = √5 ≈ 2.2360679775

Bước 8: Áp trọng số. Mẫu Simpson với n = 4 là 1, 4, 2, 4, 1:

  • 1 × g(x0) = 1 × 1.0000000000 = 1.0000000000
  • 4 × g(x1) = 4 × 1.1180339887 = 4.4721359550
  • 2 × g(x2) = 2 × 1.4142135624 = 2.8284271247
  • 4 × g(x3) = 4 × 1.8027756377 = 7.2111025508
  • 1 × g(x4) = 1 × 2.2360679775 = 2.2360679775

Bước 9: Tổng có trọng số. Cộng lại được 1.0000000000 + 4.4721359550 + 2.8284271247 + 7.2111025508 + 2.2360679775 = 17.7477336082.

Bước 10: Kết quả cuối cùng. Nhân với h/3 = 0.25/3 ≈ 0.0833333333: L ≈ 0.0833333333 × 17.7477336082 ≈ 1.4789778007.

So với giá trị chính xác 1.4789428575, sai số khoảng 0.0000349, tức xấp xỉ 3.5 × 10-5. Với chỉ 4 khoảng con (5 lần tính hàm), Quy tắc Simpson đã đạt độ chính xác 4 chữ số thập phân. Đây là điểm mạnh của phương pháp: hội tụ nhanh trên các hàm trơn.

Để so sánh, tăng gấp đôi lên n = 6 cho 1.4789405084 (sai số ≈ 2.4 × 10-6), và ngay cả n = 2 thô cũng cho 1.4821537045 (sai số ≈ 0.0032). Mỗi lần tăng gấp đôi n đều giảm sai số mạnh mẽ, đặc trưng của hội tụ bậc bốn của Simpson.

Phép xấp xỉ parabol của Quy tắc Simpson qua ba điểm mẫu trên một đường cong trơn

Hiểu về độ chính xác của Quy tắc Simpson

Mọi phương pháp số đều cho ra một phép xấp xỉ, không phải đáp án chính xác. Hiệu giữa phép xấp xỉ và giá trị thực được gọi là sai số. Hiểu cách sai số này hoạt động giúp bạn chọn thiết lập phù hợp và tin tưởng vào kết quả của mình.

Quy tắc Simpson là phương pháp bậc bốn, nghĩa là sai số co lại tỉ lệ với h4 (hay tương đương 1/n4) khi độ dài bước giảm. Trên thực tế, tăng gấp đôi số khoảng con thường giảm sai số khoảng 16 lần. Đó là lý do ví dụ minh họa trên cải thiện mạnh mẽ từ n = 2 lên n = 4 rồi n = 6.

Giới hạn sai số cổ điển cho Quy tắc Simpson tổng hợp là:

|Es| ≤ M(b - a)5 / (180n4)

Ở đây, M là một cận trên của giá trị tuyệt đối của đạo hàm cấp bốn của hàm dưới dấu tích phân g(x) trên khoảng [a, b]. Tức là, M ≥ |g(4)(x)| với mọi x trong [a, b]. Đạo hàm cấp bốn đo tốc độ thay đổi của chính độ cong của hàm dưới dấu tích phân. Khi g rất trơn (đạo hàm cấp bốn của nó nhỏ), giới hạn đảm bảo sai số nhỏ. Khi g dao động nhanh, M lớn và giới hạn kém thuyết phục hơn.

Có một số lưu ý quan trọng. Thứ nhất, giới hạn này đòi hỏi g có đạo hàm cấp bốn liên tục trên khoảng. Nếu hàm dưới dấu tích phân có điểm gãy, điểm nhọn hay gián đoạn ở đạo hàm cấp cao hơn, giới hạn không áp dụng. Thứ hai, trên thực tế việc tính M thường khó khăn, nên giới hạn hữu ích hơn cho hiểu biết lý thuyết so với ước lượng sai số hằng ngày. Thứ ba, giới hạn có thể rất bi quan: sai số thực tế thường nhỏ hơn nhiều so với công thức dự báo.

Trên thực tế, cách kiểm tra độ chính xác đáng tin cậy nhất là kiểm tra hội tụ. Tính độ dài cung với n khoảng con, rồi với 2n, rồi với 4n. Nếu các kết quả khớp nhau đến số chữ số thập phân bạn cần, ước lượng đáng tin cậy. Nếu chúng tiếp tục thay đổi đáng kể, hãy tăng n thêm nữa hoặc điều tra xem hàm dưới dấu tích phân có hành vi bất thường không (khúc cua gắt, gần điểm kỳ dị, dao động nhanh).

Độ trơn của hàm gốc f ảnh hưởng gián tiếp qua g. Nếu f có góc nhọn, f' nhảy gián đoạn, và g = √(1 + (f')2) thừa hưởng một điểm gãy. Quy tắc Simpson vẫn hoạt động, nhưng hội tụ chậm lại. Trong những trường hợp như vậy, chia khoảng tại điểm có vấn đề và áp dụng Quy tắc Simpson riêng cho từng mảnh trơn sẽ khôi phục hội tụ nhanh.

Mẫu trọng số 1 4 2 4 2 4 1 của Quy tắc Simpson tổng hợp trên các điểm cách đều

Quy tắc Simpson so với quy tắc hình thang

Quy tắc hình thang là họ hàng gần nhất và là lựa chọn thay thế phổ biến nhất của Quy tắc Simpson. Cả hai đều là công thức Newton-Cotes dùng các điểm cách đều, nhưng chúng khác nhau căn bản ở cách xấp xỉ hàm dưới dấu tích phân giữa các điểm mẫu.

Quy tắc hình thang nối các điểm mẫu liên tiếp bằng đường thẳng, tạo thành các hình thang. Quy tắc Simpson khớp các parabol qua các nhóm ba điểm. Vì parabol có thể cong theo hàm dưới dấu tích phân trong khi đường thẳng thì không, Quy tắc Simpson nắm bắt được độ cong mà quy tắc hình thang bỏ lỡ. Sự khác biệt hình học này chuyển trực tiếp thành độ chính xác: Quy tắc Simpson bậc bốn (sai số tỉ lệ với 1/n4) trong khi quy tắc hình thang bậc hai (sai số tỉ lệ với 1/n2).

Trên thực tế, điều này có nghĩa là Quy tắc Simpson thường đạt được với n = 20 điều mà quy tắc hình thang cần n = 200 trở lên mới sánh kịp, ít nhất với các hàm trơn. Cái giá là Quy tắc Simpson đòi hỏi n chẵn và việc ghi chép cẩn thận hơn một chút với các trọng số xen kẽ. Quy tắc hình thang hoạt động với mọi n và có mẫu trọng số đơn giản hơn (1, 2, 2, ..., 2, 1).

Không phương pháp nào vượt trội tuyệt đối. Quy tắc hình thang có một ưu điểm bất ngờ với các hàm tuần hoàn lấy tích phân trên cả chu kỳ, nơi nó có thể hội tụ nhanh theo cấp số nhân. Quy tắc Simpson có thể gặp khó với các hàm không đủ trơn để lý thuyết bậc bốn của nó áp dụng. Với dữ liệu thô hoặc nhiễu, quy tắc hình thang đơn giản hơn đôi khi vững chắc hơn.

Tiêu chíQuy tắc SimpsonQuy tắc hình thang
Dạng xấp xỉCác đoạn parabolCác đoạn thẳng
Bậc hội tụBậc bốn (1/n4)Bậc hai (1/n2)
Yêu cầu chia khoảngn phải chẵnMọi n dương
Mẫu trọng số1, 4, 2, 4, ..., 4, 11, 2, 2, ..., 2, 1
Số lần tính hàmn + 1n + 1
Phù hợp nhất choHàm trơn cần độ chính xác caoTrường hợp đơn giản, hàm tuần hoàn, dữ liệu thô
Độ trơn cần thiếtĐạo hàm cấp bốn liên tục cho lý thuyết đầy đủĐạo hàm cấp hai liên tục cho lý thuyết đầy đủ

Một quy trình làm việc hợp lý là dùng Quy tắc Simpson làm phương pháp chính và kiểm tra chéo bằng quy tắc hình thang ở độ phân giải cao hơn khi hành vi của hàm dưới dấu tích phân chưa chắc chắn. Nếu cả hai phương pháp đều đồng ý, độ tin cậy cao. Nếu chúng khác nhau đáng kể, hàm dưới dấu tích phân có lẽ cần nhiều khoảng con hơn hoặc chia khoảng.

So sánh phép xấp xỉ parabol của Simpson với phép xấp xỉ tuyến tính của quy tắc hình thang trên một đường cong trơn

Các cách tiếp cận toán học khác nhau cho độ dài cung

Các đường cong có thể được mô tả trong nhiều hệ tọa độ, và mỗi hệ có công thức độ dài cung riêng. Quy tắc Simpson áp dụng cho tất cả, vì mỗi công thức là một tích phân xác định của một hàm trơn.

Độ dài cung Descartes

Với y = f(x) trên [a, b], công thức là L = ∫ab √(1 + (dy/dx)2) dx. Đây là dạng chuẩn được dạy trong các khóa học giải tích. Quy tắc Simpson lấy tích phân của g(x) = √(1 + (f'(x))2) như đã mô tả xuyên suốt hướng dẫn này. Nó hoạt động tốt khi f khả vi với đạo hàm liên tục.

Độ dài cung tham số

Với đường cong mô tả bởi x = x(t) và y = y(t) với t từ a đến b, công thức là L = ∫ab √((dx/dt)2 + (dy/dt)2) dt. Dạng tham số xử lý các đường cong không qua được kiểm tra đường thẳng đứng, như đường tròn, elip và xycloit, những đường không thể viết dưới dạng y = f(x). Quy tắc Simpson lấy tích phân của hàm tốc độ tham số √((x'(t))2 + (y'(t))2), biểu thị tốc độ một điểm di chuyển dọc đường cong khi t tăng.

Độ dài cung cực

Với r = r(θ) cùng θ từ α đến β, công thức là L = ∫αβ √(r2 + (dr/dθ)2) dθ. Tọa độ cực mô tả đường xoắn ốc, đường cardioid và đường hoa hồng một cách tự nhiên. Hàm dưới dấu tích phân kết hợp khoảng cách xuyên tâm r với tốc độ thay đổi theo góc dr/dθ. Quy tắc Simpson áp dụng trực tiếp, lấy tích phân theo θ.

Độ dài cung ba chiều

Với đường cong không gian x = x(t), y = y(t), z = z(t), công thức mở rộng tự nhiên: L = ∫ab √((dx/dt)2 + (dy/dt)2 + (dz/dt)2) dt. Mỗi tọa độ đóng góp thành phần vận tốc bình phương của nó. Điều này thiết yếu cho các đường xoắn ốc, quỹ đạo 3D và các đường cong trong đồ họa máy tính. Quy tắc Simpson xử lý hàm ba chiều giống như các phiên bản 2D.

Trong mọi trường hợp, chiến lược đều giống nhau: xác định hàm dưới dấu tích phân, xác nhận nó trơn trên khoảng, chọn n chẵn, và áp dụng tổng có trọng số. Máy tính trên trang này hỗ trợ các chế độ hàm số, tham số, cực và tham số 3D, mỗi chế độ đều được định tuyến qua Quy tắc Simpson một cách phù hợp.

Hình học độ dài cung thể hiện vi phân đoạn ds với các thành phần dx và dy trên một đường cong

Ứng dụng thực tế của Quy tắc Simpson

Quy tắc Simpson cho độ dài cung xuất hiện trong nhiều lĩnh vực, nơi khoảng cách cong phải được ước lượng từ mô hình toán học hoặc dữ liệu đo đạc.

Kỹ thuật cơ khí. Biên dạng cam, đường cong răng bánh răng và ray dẫn hướng được thiết kế dưới dạng các đường cong toán học trơn. Ước lượng độ dài của chúng giúp xác định nhu cầu vật liệu, khoảng cách tiếp xúc và mẫu mòn. Quy tắc Simpson cung cấp các kiểm tra số nhanh trong chu kỳ lặp thiết kế.

Hình học đường bộ và giao thông. Các đường cong xa lộ, đoạn chuyển tiếp đường sắt và bố trí đường băng tuân theo các biên dạng hình học được thiết kế cẩn thận (đường clothoid, cung tròn, đường cong đứng parabol). Tính độ dài thực của các tuyến này quan trọng cho ước tính chi phí, đặt biển báo và lập kế hoạch thi công.

CAD và thiết kế có sự hỗ trợ của máy tính. Phần mềm thiết kế biểu diễn các đường cong dưới dạng spline và đoạn Bezier. Tích phân số ước lượng độ dài của các đường cong tổng hợp này cho sản xuất, nơi đường cắt và lượng vật liệu dùng phụ thuộc vào dữ liệu độ dài chính xác.

Robot và lập kế hoạch chuyển động. Cánh tay robot hay xe tự hành đi theo một quỹ đạo đã lập kế hoạch trong không gian. Tổng độ dài đường đi ảnh hưởng đến thời gian di chuyển, tiêu thụ năng lượng và lập kế hoạch nhiệm vụ. Quy tắc Simpson lấy tích phân biên dạng tốc độ dọc theo các quỹ đạo tham số hóa.

Tính toán khoa học. Các nhà nghiên cứu mô hình hóa hiện tượng vật lý, từ quỹ đạo hạt đến đường sức, thường xuyên cần độ dài đường cong. Khi các phương trình điều khiển không cho độ dài cung dạng đóng, tích phân số là cách tiếp cận chuẩn.

Kỹ thuật kết cấu. Cầu vòm, dầm cong và kết cấu vỏ có độ dài và số đo bề mặt suy ra từ định nghĩa hình học của chúng. Tính độ dài cung chính xác hỗ trợ phân tích tải trọng và bóc tách khối lượng vật liệu.

Mô hình toán học và nghiên cứu. Trong toán học thuần túy và ứng dụng, độ dài cung xuất hiện trong hình học vi phân, các bài toán biến phân và tối ưu hóa. Các ước lượng số cung cấp trực giác và kiểm chứng cho các kết quả lý thuyết.

Trong tất cả các bối cảnh này, Quy tắc Simpson được đánh giá cao vì mang lại độ chính xác cao với tương đối ít lần tính hàm, khiến nó hiệu quả cho cả công cụ tương tác và tính toán quy mô lớn.

Các lỗi thường gặp khi dùng Quy tắc Simpson

Có một số lỗi thường gặp khiến học sinh và người hành nghề vấp phải. Nhận biết các cạm bẫy này giúp ngăn ngừa hầu hết các kết quả tệ.

Dùng số khoảng con lẻ. Quy tắc Simpson 1/3 về cơ bản đòi hỏi n chẵn. Với n lẻ, cấu trúc cặp bị phá vỡ. Luôn xác minh n chẵn trước khi tính. Nếu bài toán của bạn gợi ý số lẻ, hãy làm tròn lên số chẵn tiếp theo.

Đạo hàm sai. Hàm dưới dấu tích phân dùng f'(x), không phải f(x). Lỗi thường gặp là lấy đạo hàm sai (ví dụ viết đạo hàm của x2 là x thay vì 2x) hoặc quên quy tắc dây chuyền với hàm hợp. Hãy kiểm tra kỹ mọi đạo hàm.

Áp Quy tắc Simpson cho f thay vì g. Liên quan đến điều trên, một số người dùng áp trọng số trực tiếp cho các giá trị của hàm gốc f(xi) thay vì hàm dưới dấu tích phân của độ dài cung g(xi) = √(1 + (f'(xi))2). Việc này tính ra một tích phân hoàn toàn sai.

Cận khoảng sai. Nhập a > b, hoặc các cận nằm ngoài miền xác định của hàm, sẽ cho kết quả vô nghĩa. Với các hàm có miền xác định bị hạn chế (như √(x) đòi hỏi x ≥ 0), hãy xác minh khoảng hợp lệ.

Quá ít khoảng con. Dù Quy tắc Simpson hội tụ nhanh, n = 2 hoặc n = 4 có thể không đủ với các hàm cong mạnh. Luôn thực hiện kiểm tra hội tụ với n lớn hơn.

Nhầm lẫn độ và radian. Khi f liên quan đến hàm lượng giác, hãy xác nhận máy tính dùng radian hay độ. Trộn lẫn chúng sẽ âm thầm làm hỏng đạo hàm và kết quả. Hầu hết các công cụ toán học dùng radian.

Cho rằng kết quả chính xác. Quy tắc Simpson chỉ chính xác với các đa thức bậc không quá 3. Với các hàm tổng quát, kết quả là một phép xấp xỉ. Coi nó là chính xác mà không có bằng chứng hội tụ là một sai lầm.

Bỏ qua điểm kỳ dị. Nếu f' bùng nổ trong [a, b] (ví dụ gần tiếp tuyến đứng), hàm dưới dấu tích phân g trở nên lớn và Quy tắc Simpson có thể hội tụ chậm hoặc không đáng tin cậy. Hãy chia khoảng hoặc dùng các phương pháp chuyên biệt.

Khi nào nên dùng Quy tắc Simpson?

Quy tắc Simpson là lựa chọn mặc định tuyệt vời cho việc tính độ dài cung bằng số, nhưng biết khi nào nó tỏa sáng và khi nào nên tìm hướng khác sẽ khiến bạn trở thành người giải quyết vấn đề hiệu quả hơn.

Dùng Quy tắc Simpson khi:

  • Hàm dưới dấu tích phân g(x) = √(1 + (f'(x))2) trơn trên khoảng, không có điểm gãy hay kỳ dị.
  • Không tồn tại nguyên hàm sơ cấp, hoặc việc tìm ra nó tốn công hơn mức một ước lượng số xứng đáng.
  • Bạn cần độ chính xác tốt với số lần tính hàm vừa phải.
  • Bạn đang tính độ dài cung cho bài tập, kiểm tra kỹ thuật hay thẩm định thiết kế.
  • Đường cong được mô tả tham số hoặc trong tọa độ cực và hàm thu được trơn.

Cân nhắc phương án khác hoặc thận trọng hơn khi:

  • Hàm dưới dấu tích phân có điểm gián đoạn hoặc góc nhọn trong khoảng. Hãy chia khoảng tại điểm có vấn đề trước.
  • Hàm số dao động nhanh. Bạn có thể cần n rất lớn, hoặc một phương pháp dành cho tích phân dao động.
  • Bạn chỉ có dữ liệu đo đạc nhiễu thay vì công thức trơn. Quy tắc hình thang hoặc các cách tiếp cận dựa trên spline có thể phù hợp hơn.
  • Cần độ chính xác cực cao. Phép cầu phương Gauss hoặc Quy tắc Simpson thích nghi (tự động làm mịn khoảng con ở nơi cần thiết) có thể tốt hơn.
  • Khoảng vô hạn hoặc bán vô hạn. Hãy biến đổi tích phân về khoảng hữu hạn trước, hoặc dùng phép cầu phương chuyên biệt.

Với hầu hết sinh viên giải tích và kỹ sư đang làm việc, Quy tắc Simpson kèm kiểm tra hội tụ (so sánh n, 2n và 4n) mang lại mọi độ chính xác cần thiết. Nó đơn giản để cài đặt, có sẵn rộng rãi trong phần mềm, và được hiểu rõ về mặt lý thuyết.

Câu hỏi thường gặp về Quy tắc của Simpson

Quy tắc Simpson dùng để làm gì? +

Quy tắc Simpson xấp xỉ các tích phân xác định bằng số. Với độ dài cung, nó ước lượng L = ∫<sub>a</sub><sup>b</sup> √(1 + (f'(x))<sup>2</sup>) dx bằng cách tính hàm dưới dấu tích phân tại các điểm cách đều rồi kết hợp chúng với mẫu trọng số 1, 4, 2, 4, ..., 4, 1. Nó được dùng rộng rãi khi không tồn tại nguyên hàm sơ cấp hoặc khi một ước lượng số nhanh là đủ.

Quy tắc Simpson tính độ dài cung như thế nào? +

Nó áp công thức trọng số Simpson cho hàm dưới dấu tích phân của độ dài cung g(x) = √(1 + (f'(x))<sup>2</sup>), chứ không phải hàm gốc. Khoảng [a, b] được chia thành số chẵn n khoảng con có độ rộng h = (b - a)/n. Ước lượng là L ≈ (h/3)[g(x<sub>0</sub>) + 4g(x<sub>1</sub>) + 2g(x<sub>2</sub>) + ... + 4g(x<sub>n-1</sub>) + g(x<sub>n</sub>)].

Vì sao Quy tắc Simpson đòi hỏi số khoảng con chẵn? +

Mỗi panel parabol trong Quy tắc Simpson trải dài hai khoảng con (ba điểm). Phương pháp nhóm các khoảng con thành từng cặp, khớp một parabol trên mỗi cặp. Nếu n lẻ, khoảng con cuối cùng sẽ không có cặp, và cấu trúc cặp sẽ sụp đổ. n chẵn đảm bảo mỗi khoảng con thuộc về đúng một panel parabol.

Quy tắc Simpson 1/3 và quy tắc hình thang khác nhau thế nào? +

Quy tắc Simpson 1/3 khớp các parabol qua các nhóm ba điểm, đạt độ chính xác bậc bốn (sai số tỉ lệ với 1/n<sup>4</sup>). Quy tắc hình thang nối các điểm bằng đường thẳng, chỉ đạt độ chính xác bậc hai (sai số tỉ lệ với 1/n<sup>2</sup>). Quy tắc Simpson thường cần ít khoảng con hơn nhiều cho cùng độ chính xác trên các hàm trơn, nhưng nó đòi hỏi n chẵn.

Cách tính độ dài bước trong Quy tắc Simpson? +

Độ dài bước là h = (b - a)/n, trong đó a và b là cận tích phân dưới và trên, còn n là số (chẵn) khoảng con. Ví dụ, trên [0, 1] với n = 4, h = 0.25. h càng nhỏ nghĩa là càng nhiều panel và thường độ chính xác càng tốt.

Quy tắc Simpson có chính xác hơn quy tắc hình thang không? +

Với các hàm trơn, đúng vậy, thường với biên độ lớn. Hội tụ bậc bốn của Simpson nghĩa là tăng gấp đôi n giảm sai số khoảng 16 lần, so với 4 lần của quy tắc hình thang. Tuy nhiên, ưu thế này giả định độ trơn đủ. Với các hàm thô, nhiễu hoặc tuần hoàn, quy tắc hình thang đôi khi hoạt động tương đương hoặc tốt hơn.

Tôi nên dùng bao nhiêu khoảng con? +

Bắt đầu với một số chẵn vừa phải như 20 hoặc 50, rồi tăng gấp đôi và so sánh kết quả. Nếu ước lượng hầu như không đổi, nó đã hội tụ và bạn có thể tin tưởng. Nếu nó thay đổi đáng kể, hãy tiếp tục tăng gấp đôi. Với các hàm trơn điển hình trong bài toán độ dài cung, n = 50 đến 200 thường là quá đủ.

Quy tắc Simpson có tính được độ dài đường cong tham số không? +

Có. Với x = x(t), y = y(t) trên [a, b], áp dụng Quy tắc Simpson cho hàm tốc độ tham số √((dx/dt)<sup>2</sup> + (dy/dt)<sup>2</sup>). Phương pháp hoàn toàn giống: chọn n chẵn, tính tại các giá trị t cách đều, áp trọng số 1-4-2-...-4-1, rồi nhân với h/3.

Quy tắc Simpson có dùng được cho đường cong ba chiều không? +

Có. Với đường cong không gian x(t), y(t), z(t), hàm dưới dấu tích phân của độ dài cung trở thành √((dx/dt)<sup>2</sup> + (dy/dt)<sup>2</sup> + (dz/dt)<sup>2</sup>). Quy tắc Simpson lấy tích phân của nó y hệt trường hợp 2D. Máy tính trên trang này có chế độ tham số 3D.

Công thức sai số của Quy tắc Simpson là gì? +

Giới hạn sai số của Quy tắc Simpson tổng hợp là |E<sub>s</sub>| ≤ M(b - a)<sup>5</sup>/(180n<sup>4</sup>), trong đó M là cận trên của |g<sup>(4)</sup>(x)|, giá trị tuyệt đối của đạo hàm cấp bốn của hàm dưới dấu tích phân, trên [a, b]. Điều này đòi hỏi g khả vi liên tục bốn lần. Trên thực tế, kiểm tra hội tụ (so sánh n, 2n, 4n) hữu ích hơn việc tính M.

Quy tắc Simpson có cho đáp án chính xác không? +

Chỉ với các hàm là đa thức bậc không quá 3. Với các hàm trơn tổng quát, nó cho một phép xấp xỉ có sai số co lại như 1/n<sup>4</sup>. Ví dụ minh họa trong hướng dẫn này (y = x<sup>2</sup> trên [0,1]) đạt sai số dưới 0.00004 chỉ với n = 4, cho thấy phép xấp xỉ trở nên gần như chính xác nhanh đến mức nào.

Làm sao kiểm tra kết quả số của tôi có đáng tin cậy không? +

Dùng kiểm tra hội tụ: tính với n, rồi 2n, rồi 4n khoảng con. Nếu các kết quả khớp nhau đến số chữ số thập phân bạn yêu cầu, ước lượng đáng tin cậy. Để kiểm tra thêm, hãy so sánh với quy tắc hình thang ở độ phân giải cao. Sự đồng thuận giữa hai phương pháp khác nhau là dấu hiệu mạnh của tính đúng đắn.