Romberg Integration EP1: จาก Trapezoid สู่การกำจัด error

ภาพรวมของบทเรียน

Romberg Integration เริ่มจาก Trapezoidal Rule ที่แบ่งช่วงละเอียดขึ้นเป็นลำดับ แล้วใช้ Richardson Extrapolation กำจัดพจน์ error ทีละระดับ

เป้าหมาย: ให้ผู้อ่านเห็นทั้งภาพรวมทางเรขาคณิต ที่มาของสูตร และเหตุผลว่าค่าประมาณดีขึ้นได้อย่างไร

แนวคิดสำคัญ

สร้าง T₁, T₂, T₄
มองรูปแบบ error ของ Trapezoid
สร้างค่า R₂,₂ จากสองระดับแรก
อธิบายว่าทำไมค่าผสมจึงแม่นกว่าเดิม

Romberg Integration เริ่มจาก Trapezoidal Rule หลายระดับ เช่น \(n=1,2,4,8\) แล้วนำค่าที่ได้มาผสมกันเพื่อลบพจน์ error หลัก แนวคิดนี้เรียกว่า Richardson Extrapolation

ภาพใหญ่คือ Trapezoid เมื่อแบ่งละเอียดขึ้นจะเข้าใกล้ค่าจริงอย่างมีรูปแบบ ถ้าเรารู้ว่า error หลักลดลงประมาณ 4 เท่าเมื่อ \(h\) ลดครึ่งหนึ่ง ก็สามารถผสมค่าหยาบและค่าละเอียดเพื่อหัก error ส่วนนั้นออกได้

แก่น EP1: อย่ามอง \(T_1,T_2,T_4\) เป็นคำตอบแยกกัน แต่มองเป็นข้อมูลหลายระดับที่ใช้เดารูปแบบ error ได้

สูตรที่ใช้

\[R_{k,j}=R_{k,j-1}+\frac{R_{k,j-1}-R_{k-1,j-1}}{4^{j-1}-1}\]

ก่อนใช้สูตร ควรอธิบายความหมายของตัวแปร เงื่อนไขจำนวนช่วง และสิ่งที่สูตรกำลังประมาณ ไม่ควรวางสูตรโดยไม่มีบริบท

ในตาราง Romberg ให้ \(R_{k,1}\) เป็นค่า Trapezoidal Rule ที่ใช้จำนวนช่วง \(2^{k-1}\) จากนั้นคอลัมน์ถัดไปสร้างจากค่าสองระดับก่อนหน้า

\[R_{k,j}=R_{k,j-1}+\frac{R_{k,j-1}-R_{k-1,j-1}}{4^{j-1}-1}\]

ตัวประกอบ \(4^{j-1}\) มาจาก error ของ Trapezoid ที่มีพจน์เป็นกำลังคู่ของ \(h\) เมื่อเลื่อนไปคอลัมน์ถัดไป เรากำลังกำจัด error ชั้นถัดไปทีละระดับ

ตัวอย่างทีละขั้น

ใช้ ∫₀¹ 1/(1+x²) dx และเทียบกับ π/4 เพื่อเห็นความแม่นเพิ่มขึ้น

กำหนดฟังก์ชันและช่วงที่ต้องการประมาณ
เลือกจำนวนช่วงหรือระดับการคำนวณ
หาค่าฟังก์ชันที่จุดสำคัญ
แทนค่าลงในสูตรโดยแสดงทุกบรรทัด
เปรียบเทียบกับค่าจริงหรือวิธีอื่นถ้ามี

ใช้ \(\int_0^1 \dfrac{1}{1+x^2}\,dx=\dfrac{\pi}{4}\approx0.785398\) ค่า Trapezoid หลายระดับคือ \(R_{1,1}=0.750000\), \(R_{2,1}=0.775000\), \(R_{3,1}=0.782794\)

\[R_{2,2}=0.775000+\frac{0.775000-0.750000}{3}=0.783333\]

\[R_{3,2}=0.782794+\frac{0.782794-0.775000}{3}=0.785392\]

ระดับ	\(R_{k,1}\)	\(R_{k,2}\)
1	0.750000	-
2	0.775000	0.783333
3	0.782794	0.785392

จะเห็นว่าเพียงผสมค่า Trapezoid สองระดับ ค่า \(R_{3,2}\) เข้าใกล้ \(\pi/4\) มากกว่า Trapezoid ดิบอย่างชัดเจน

สิ่งที่ควรสังเกตจากภาพ

ส่วนนี้ใช้สำหรับแทรกภาพนิ่งจาก Manim หรือภาพเปรียบเทียบสำคัญ เช่น จุดตัวอย่าง พื้นที่ที่ประมาณเกิน–ขาด ตารางค่า หรือเส้นทางของ error

IMAGE SLOT: เพิ่มภาพจากคลิปภายหลัง พร้อม alt text และคำอธิบายใต้ภาพ

ข้อควรระวัง

อย่าใช้สูตรโดยไม่ตรวจเงื่อนไขจำนวนช่วง
อย่าปัดเศษเร็วเกินไปในขั้นกลาง
ต้องแยก “ค่าประมาณ” ออกจาก “ค่าจริง” ให้ชัด
ควรอธิบายว่า error ที่กล่าวถึงวัดจากอะไร

ต้องมีลำดับระดับ: Romberg ต้องการค่า Trapezoid ที่ลด \(h\) อย่างเป็นระบบ เช่นจำนวนช่วง 1,2,4,8 ไม่ใช่เลือก \(n\) สะเปะสะปะ

หากฟังก์ชันไม่เรียบหรือมี singularity ในช่วง รูปแบบ error ที่ Romberg คาดไว้จะเสีย ทำให้ extrapolation ไม่เสถียร ควรตรวจโดเมนและไม่เร่งสร้างคอลัมน์ลึกเกินจากข้อมูลน้อยเกินไป

สรุป

เมื่อมีคลิปแล้ว ให้วางคลิปด้านบนเพื่อใช้ภาพเคลื่อนไหวอธิบาย intuition ส่วนหน้านี้เก็บที่มา สูตร การคำนวณ และรายละเอียดที่คลิปสั้นใส่ไม่หมด