วิธี 1/t ที่ปรับปรุงแล้วสำหรับการหาปริพันธ์เชิงตัวเลขของฟังก์ชันที่คำนวณได้ยาก

ปัญหาที่งานวิจัยแก้ไข

การสุ่มมอนติคาร์โลทั่วไปอาจไม่พบพื้นที่แคบที่ฟังก์ชันมียอดสูงมาก จึงประเมินปริพันธ์ที่คำนวณยากต่ำกว่าความเป็นจริงอย่างมาก แม้อัลกอริทึม 1/t แบบเดิมสำรวจพื้นที่ดังกล่าวได้ดีขึ้น แต่การใช้ค่าคงที่แทนแต่ละช่วงทำให้เกิดความคลาดเคลื่อนจากความกว้างช่วงและเกิดภาวะความคลาดเคลื่อนอิ่มตัว

ผู้ได้รับประโยชน์:

• นักวิทยาศาสตร์เชิงคำนวณที่ประเมินปริพันธ์หลายมิติซึ่งคำนวณได้ยาก
• นักวิจัยด้านมอนติคาร์โลและความหนาแน่นของสถานะ
• นักคณิตศาสตร์ประยุกต์ที่ศึกษาการลู่เข้าเชิงตัวเลข
• นักวิเคราะห์ฟิสิกส์และวิศวกรรมที่ทำงานกับฟังก์ชันยอดแหลม

วิธีการ

ตัวประมาณที่เสนอคงการเดินสุ่มแบบ 1/t ไว้เพื่อประมาณความหนาแน่นของสถานะ g(y) แต่แทนค่าคงที่ประจำแต่ละช่วงในปริภูมิ y ด้วยค่าเฉลี่ยของค่าฟังก์ชันที่สุ่มพบจริงภายในช่วงนั้นและปรับค่าอย่างต่อเนื่อง ค่าปริพันธ์จึงประมาณจากผลรวมของ g(y) คูณค่าเฉลี่ยภายในช่วง งานวิจัยทดสอบกับปริพันธ์เกาส์เซียนยอดแหลมหนึ่ง สอง และสามมิติ โดยกำหนดค่าเฉลี่ยเป็นศูนย์ ส่วนเบี่ยงเบนมาตรฐานเท่ากับรากที่สองของ 0.4 และโดเมน [-10, 10] ในแต่ละมิติ จากนั้นเปรียบเทียบการสุ่มอย่างง่าย อัลกอริทึม 1/t แบบเดิม และวิธีที่เสนอภายใต้ความกว้างช่วงสี่ระดับ การทดลอง 10^10 ครั้งต่อรอบ และการจำลองอิสระ 50 ชุด โดยใช้ค่าปริพันธ์จริงและความคลาดเคลื่อนเชิงสัดส่วน

ผลการศึกษา

การสุ่มอย่างง่ายประเมินปริพันธ์ยอดแหลมต่ำกว่าความเป็นจริงอย่างมาก ขณะที่ค่าประมาณจาก 1/t แบบเดิมยังไวต่อความกว้างช่วงและเกิดภาวะความคลาดเคลื่อนอิ่มตัว เมื่อใช้ความกว้างช่วง 0.0005 วิธีที่เสนอให้ค่าประมาณกรณีหนึ่ง สอง และสามมิติเท่ากับ 1.585321, 2.513090 และ 3.984088 ซึ่งใกล้ค่าจริง 1.585331, 2.513274 และ 3.984371 ตามลำดับ ความคลาดเคลื่อนลดลงโดยประมาณตาม N^-1/2 โดยไม่อิ่มตัว และเมื่อความแม่นยำไม่ถูกกำหนดด้วยความกว้างช่วงแล้ว จึงใช้ช่วงที่กว้างเพื่อเข้าถึงทุกสถานะได้เร็วขึ้นและลดภาระการคำนวณ