Tôi muốn sử dụng một chuỗi bán ngẫu nhiên, cụ thể là Sobol, trong một mô phỏng dựa trên SciPy. Bất kỳ khuyến nghị nào về các gói hiện có, hiệu quả?Đề xuất cho sự khác biệt thấp (ví dụ: Sobol) các chuỗi bán ngẫu nhiên trong Python/SciPy?
Có phải là this project sử dụng nào không? Nếu không, các relevant Wikipedia article đề cập đến C và Fortran thói quen trong NAG và Numerical Bí quyết, mà có lẽ không nên khó khăn để bọc hoặc reimplement trong Python. Ngoài ra còn có C thói quen trong GSL.
Trong bối cảnh phân tích nhạy cảm, thư viện SALIB có vẻ thú vị Nó có trình tạo mẫu Sobol và sử dụng SciPy. Liên kết tại đây: http://jdherman.github.io/SALib/
Cảm ơn, dự án của jburkardt (được xây dựng trên NumPy) trông thú vị. Bạn (hoặc những người khác) có kinh nghiệm trực tiếp nào không? Hoặc bất kỳ so sánh với các gói khác? Ngoài ra, nó sẽ là tuyệt vời để không có giới hạn kích thước 40. Vì vậy, tôi quan tâm tìm thấy hoặc là một thực hiện python hoặc (có lẽ tốt hơn, b/c nhanh hơn) một wrapper python được tạo sẵn cho một mã C/C++ loại bỏ hạn chế này, như được thực hiện trong C++ của Jburkardt là Fortran. –
OK ... Tôi sẽ chấp nhận câu trả lời này ngay bây giờ, có vẻ như đây là lựa chọn duy nhất. Nhưng tôi vẫn quan tâm đến các lựa chọn hoặc giải pháp khác cho các vấn đề tôi nêu ra trong phần bình luận cuối cùng. –
Tôi đã dọn dẹp dự án được liên kết ở trên (tài liệu hiện đang hoạt động) và sẽ tiếp tục cập nhật dự án nếu cần. Nó có sẵn tại https://github.com/naught101/sobol_seq – naught101