Tích phân Monte Carlo
Tích phân Monte Carlo ước tính một tích phân xác định bằng cách lấy trung bình của hàm dưới dấu tích phân trên các điểm mẫu ngẫu nhiên, chuyển đổi tích phân thành ước tính một kỳ vọng.
Definition
Tích phân Monte Carlo là phép xấp xỉ một tích phân bằng cách viết nó dưới dạng kỳ vọng của một hàm dưới một phân phối lấy mẫu và ước tính kỳ vọng đó bằng giá trị trung bình mẫu trên các lần lấy mẫu từ phân phối.
Scope
Chủ đề này bao gồm việc biểu diễn một tích phân dưới dạng kỳ vọng, ước lượng Monte Carlo thô (đơn giản) và tính không chệch của nó, tốc độ hội tụ căn bậc hai của n và tính độc lập với số chiều, ước tính sai số thông qua độ lệch chuẩn mẫu, và so sánh với phép cầu phương xác định. Các cải tiến giảm phương sai được coi là các phần mở rộng được đề cập ở nơi khác.
Core questions
- Làm thế nào để một tích phân tùy ý được biểu diễn dưới dạng kỳ vọng phù hợp để lấy mẫu?
- Tại sao ước lượng Monte Carlo thô lại không chệch và nhất quán?
- Điều gì chi phối tốc độ sai số căn bậc hai của n, và tại sao nó độc lập với số chiều?
- Khi nào tích phân Monte Carlo vượt trội hơn phép cầu phương xác định?
Key concepts
- Ước lượng Monte Carlo thô
- Tính không chệch
- Sai số chuẩn
- Tính độc lập với số chiều
- Mật độ lấy mẫu
Key theories
- Tích phân dưới dạng kỳ vọng
- Viết một tích phân dưới dạng kỳ vọng của hàm dưới dấu tích phân chia cho mật độ lấy mẫu biến tích phân thành ước tính một giá trị trung bình, mà giá trị trung bình mẫu ước tính không chệch.
- Tốc độ hội tụ và ước tính sai số
- Định lý giới hạn trung tâm đưa ra một sai số chuẩn tỷ lệ nghịch với căn bậc hai của kích thước mẫu, độc lập với số chiều của tích phân, và độ lệch chuẩn thực nghiệm của các số hạng cung cấp một ước tính sai số có thể sử dụng được.
Clinical relevance
Tích phân Monte Carlo tính toán các hằng số chuẩn hóa, kỳ vọng hậu nghiệm, xác suất biên và các kỳ vọng đa chiều phát sinh trong thống kê và khoa học vật lý; tốc độ sai số độc lập với số chiều của nó làm cho nó trở thành phương pháp được lựa chọn khi phép cầu phương dựa trên lưới trở nên không khả thi.
History
Ý tưởng ước tính tích phân bằng cách lấy mẫu có từ những tính toán của Los Alamos vào những năm 1940 và bài báo năm 1949 của Metropolis và Ulam; nó trở thành thực hành thường xuyên khi sức mạnh tính toán tăng lên và khi các nhà thống kê nhận ra lợi thế của nó so với phép cầu phương trong các không gian có số chiều cao.
Key figures
- Stanislaw Ulam
- Nicholas Metropolis
- Christian P. Robert
Related topics
Seminal works
- robert2004
- metropolis1949
Frequently asked questions
- Tích phân Monte Carlo chính xác đến mức nào?
- Sai số của nó giảm theo tỷ lệ nghịch với căn bậc hai của số lượng mẫu, vì vậy việc tăng gấp bốn lần kích thước mẫu sẽ giảm một nửa sai số. Ước lượng này cũng đi kèm với một ước tính sai số tích hợp từ độ lệch chuẩn mẫu của các giá trị hàm dưới dấu tích phân.
- Khi nào tôi nên ưu tiên Monte Carlo hơn phép cầu phương tiêu chuẩn?
- Đối với các tích phân trơn có số chiều thấp, phép cầu phương xác định thường hội tụ nhanh hơn. Monte Carlo thắng thế trong các không gian có số chiều cao, nơi chi phí của một lưới tăng theo cấp số nhân nhưng tốc độ sai số Monte Carlo vẫn giữ nguyên.