Lý thuyết nền tảng hữu ích
Một thực tế nhỏ mà bạn có thể sử dụng để giúp hiểu liệu một đạo hàm số có được tính toán chính xác hay không là phần còn lại Cauchy của bản mở rộng Taylor. Đó là,
ξ∈[x,x+h]f( x + h ) = f( x ) + h f'( x ) + h22f′ ′( ξ)ξ∈ [ x , x + h ]
Điều này rất hữu ích, vì có lẽ bạn đã xấp xỉ đạo hàm đầu tiên của mình bằng
f'( x ) ≈ f( x + h ) - f( x - h )2 h
với một số nhỏ (tôi thường sử dụng , nhưng tôi chắc chắn một ngày nào đó tôi sẽ gặp phải trường hợp không phù hợp).10 - 4h10- 4
Sau một chút đại số, chúng ta có thể sử dụng phần còn lại Cauchy để thấy rằng về mặt lý thuyết gần đúng số của chúng ta phải nằm trong của .f ′ ( x )h f′ ′( ξ) , Ξ∈ [ x - h , x + h ]f'( x )
Trên thực tế, bạn thực sự có thể ràng buộc nó bằng , trong đó và ... tương đương với , .ξ 1 ∈ [ x - h , x ] ξ 2 ∈ [ x , x + h ] h 2 f ‴ ( ξ ) ξ ∈ [ x - h , x + h ]h ( f′ ′( ξ1) - f′ ′( ξ2) )ξ1∈ [ x - h , x ]ξ2∈ [ x , x + h ]h2f′( ξ)ξ∈ [ x - h , x + h ]
Vấn đề trong thực tế
Được rồi, chúng tôi có lý thuyết tốt đẹp ràng buộc lỗi của đạo hàm số. Nhưng có hai lỗ hổng khi trực tiếp cố gắng sử dụng các kết quả đó:
1.) Chúng tôi không biết (và có lẽ không muốn dành thời gian xấp xỉ nó)f′( x )
2.) vì , bị mất ổn định sốf ( x + h ) - f ( x - h )h → 0f( x + h ) - f( x - h )2 h
Vì vậy, bằng cách sử dụng những gì chúng ta biết từ trước đó, cách tôi kiểm tra các đạo hàm phân tích của mình (có thể không phải là cách tốt nhất) là tôi viết hàm đạo hàm số dưới dạng hàm của . Nếu tôi không thể biết được sự khác biệt giữa các đạo hàm số và phân tích là do lỗi mã hóa hay chỉ là xấp xỉ số, tôi có thể giảm và xem liệu đạo hàm số của mình có tiếp cận đạo hàm phân tích của tôi hay không trước khi mất ổn định số (khi điều này xảy ra, các xấp xỉ số của bạn sẽ trở nên ít nhất quán hơn khi trở nên nhỏ hơn). Lưu ý rằng thuật ngữ sẽ biến mất theo phương trình bậc hai, vì vậy nếu lỗi của tôi là khoảng vớih h f ‴ ( ξ ) 0,01 h = 10 - 4 0,0001 h = 10 - 5hhhf′( ξ)0,01h = 10- 4 , nó sẽ ở khoảng với giả sử không ổn định số chưa được kích hoạt .0,0001h = 10- 5
Thật không may, không có hướng dẫn cứng và nhanh để luôn xác định những điều này; nó phụ thuộc rất nhiều vào mức độ ổn định của hàm (và ý tôi là cả về độ ổn định số và đạo hàm cao hơn). Nhưng theo kinh nghiệm của tôi, tôi chưa bao giờ thấy trường hợp lỗi từ không dứt khoát về 0 (tức là sử dụng đã đưa ra câu trả lời gần giống như ) tại thời điểm mất ổn định số từ .h = 10 - 4 h = 10 - 5 h → 0h2f′( ξ)h = 10- 4h = 10- 5h → 0