✨Tích phân Monte-Carlo

nhỏ|Một minh họa về tích hợp Monte-Carlo. Tích phân Monte Carlo là một phương pháp tìm giá trị số của tích phân, đặc biệt là các tích phân đa chiều có dạng:

:$I\; =\; \backslash int\_\{V\}\; f(\backslash vec\; x)\backslash ,\; d\backslash vec\; x$

trên một miền không gian đa chiều V sử dụng một số hữu hạn các lần gọi hàm f.

Các phương pháp tích phân Monte-Carlo bao gồm phương pháp cơ bản, phương pháp lấy mẫu có trọng tâm,... Các phương pháp này cũng cho biết ước lượng sai số thống kê của phép tính, tuy rằng ước lượng này có thể không chính xác do việc khảo sát ngẫu nhiên hàm số trên miền không gian đa chiều có thể không cho thấy hết mọi biểu hiện của hàm.

Tích phân Monte Carlo cơ bản

Tích phân một chiều

Ở dạng cơ bản nhất, giá trị của tích phân một chiều: :$I\; =\; \backslash int\_\{a\}^\{b\}\; f(x)\backslash ,\; dx$

được dự đoán là tổng: :$E(I)\; =\; V=\; \backslash frac\{V\}\{N\}\; \backslash sum\_\{i=1\}^N\; f(x\_i)$

trong đó :V là thể tích mở rộng của miền tích phân ::$V\; =\; \backslash int\_\{a\}^\{b\}\backslash ,\; dx$ :x_i là các giá trị lấy ngẫu nhiên đều trong khoảng [a, b]. :N là tổng số lần lấy mẫu x_i

Sai số của dự đoán được tính bằng căn của phương sai của giá trị trung bình:

Khi số lần lấy mẫu, N, tăng, phương sai giảm theo 1/N, tức là sai số của phép tính giảm theo $\backslash frac\{1\}\{\backslash sqrt\{N$.

Tích phân đa chiều

Phương pháp trên được mở rộng cho tích phân đa chiều: :

Lấy mẫu có trọng tâm

Tích phân một chiều

Nếu biết hàm cần tích phân f(x) cư xử như nào, ta có thể chọn được một hàm g(x) có giá trị biến đổi gần giống |f(x)| trên miền cần tích phân, ta có thể biến đổi tích phân thành: :$I\; =\; \backslash int${a}^{b} f(x)\, dx = \int{a}^{b} \frac{f(x)}{g(x)} g(x)\, dx = \int_{a}^{b} \frac{f(x)}{g(x)} \, dG(x) với: :$g(x)\; =\backslash frac\{dG(x)\}\{dx\}$ và g(x) thỏa mãn điều kiện chuẩn hóa: :$\backslash int\_\{a\}^\{b\}\; g(x)\backslash ,\; dx\; =\; 1$ Lúc này có thể lấy các điểm x_i ngẫu nhiên trong khoảng [a, b] theo phân bố xác suất g(x') để tìm giá trị tích phân: :$I\; \backslash approx\; \backslash frac\{1\}\{N\}\; \backslash sum\_\{i=1\}^N\; \backslash frac\{f(x\_i)\}\{g(x\_i)\}$

Hàm g(x) càng giống f(x) thì phương sai của f(x)/g(x) càng nhỏ và sai số của phép tính càng nhỏ.

Một bất lợi của phương pháp này là sai số có thể lớn nếu hàm g(x) được chọn gần bằng 0 tại những điểm mà f(x) khác 0. Lúc đó, phương sai của f(x)/g(x) có thể lớn đến vô cùng. Lỗi này có thể khó phát hiện khi miền giá trị tại đó g(x) bằng 0 là rất nhỏ.