Phương pháp Halley là gì?
Phương pháp Halley là một kỹ thuật số học lặp dùng để giải phương trình dạng \(f(x) = 0\). Đây là "người anh em" bậc ba (hội tụ bậc ba) của phương pháp Newton-Raphson: trong khi Newton chỉ dùng hàm số và đạo hàm bậc nhất, thì Halley còn tận dụng thêm đạo hàm bậc hai, nhờ đó thường đạt được độ chính xác mong muốn chỉ sau ít vòng lặp hơn. Đây là một công cụ toán học / giải tích số mang tính phổ quát, áp dụng được ở mọi nơi; các góc bên trong hàm lượng giác đều được hiểu theo đơn vị radian.
Cách sử dụng
Nhập hàm số mà bạn muốn tìm nghiệm vào ô \(f(x)\), sau đó cung cấp đạo hàm bậc nhất \(f'(x)\) và đạo hàm bậc hai \(f''(x)\) dưới dạng các biểu thức theo \(x\). Cú pháp được hỗ trợ gồm +, -, *, /, ^ cho lũy thừa, dấu ngoặc đơn, và các hàm như sin, cos, tan, asin, acos, atan, sinh, cosh, tanh, exp, log/ln, log10, sqrt và abs, cùng các hằng số pi và e. Hãy chọn một giá trị khởi đầu \(x_0\) nằm gần nghiệm bạn cần, đặt số vòng lặp tối đa \(n\), và chọn số chữ số có nghĩa muốn hiển thị.
Giải thích công thức
Mỗi bước lặp được tính theo công thức $$x_{n+1} = x_n - \frac{2\,f(x_n)\,f'(x_n)}{2\,[f'(x_n)]^2 - f(x_n)\,f''(x_n)}$$ Tử số chính là số hiệu chỉnh Newton quen thuộc (đã được nhân hệ số), còn số hạng phụ \(- f(x_n) f''(x_n)\) ở mẫu số bù trừ cho độ cong của \(f\), giúp tăng tốc độ hội tụ. Vòng lặp sẽ dừng khi mức thay đổi của \(x\) hoặc giá trị dư \(f(x)\) giảm xuống dưới một sai số rất nhỏ, hoặc khi đạt tới giới hạn số vòng lặp. Nếu mẫu số bằng 0 thì phương pháp sẽ "đổ vỡ" và bạn nên thử một giá trị \(x_0\) khác.
Ví dụ minh họa
Với \(f(x) = x - \cos(x)\), \(f'(x) = 1 + \sin(x)\), \(f''(x) = \cos(x)\), bắt đầu tại \(x_0 = 1\): bước đầu tiên cho $$x_1 = 1 - \frac{2 \times 0.4596977 \times 1.8414710}{2 \times 1.8414710^2 - 0.4596977 \times 0.5403023} = 1 - \frac{1.6930504}{6.5336550} = 0.7408769$$ Phép lặp nhanh chóng ổn định tại số Dottie \(x = 0.7390851332151607\), nghiệm duy nhất của phương trình \(x = \cos(x)\).
Câu hỏi thường gặp
Phương pháp Halley khác phương pháp Newton ở điểm nào? Phương pháp Newton bỏ qua độ cong, còn Halley bổ sung thêm số hạng đạo hàm bậc hai, cho hội tụ bậc ba thay vì bậc hai, nên thường cần ít vòng lặp hơn để đạt mỗi chữ số chính xác.
Tại sao tôi phải nhập các đạo hàm? Máy tính này tin tưởng hoàn toàn vào các đạo hàm bạn cung cấp. Nếu chúng sai, quá trình hội tụ sẽ kém hoặc thất bại. Hãy tính \(f'(x)\) và \(f''(x)\) bằng cách lấy đạo hàm của \(f(x)\) một cách cẩn thận.
Nếu nó không hội tụ thì sao? Phương pháp này tìm nghiệm gần \(x_0\) nhất. Một điểm khởi đầu không tốt có thể khiến phép lặp phân kỳ hoặc hội tụ về một nghiệm khác, còn mẫu số bằng 0 sẽ khiến nó dừng hẳn. Hãy thử đổi \(x_0\) hoặc kiểm tra lại các biểu thức đạo hàm của bạn.
