Romberg Integration EP2: สร้างตารางทีละช่อง

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

จากคอลัมน์แรกที่เป็นค่า Trapezoid เราจะสร้างคอลัมน์ถัดไปด้วยตัวหาร 3, 15 และ 63 พร้อมติดตามว่าพจน์ error ชั้นใดถูกกำจัดออก

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

แนวคิดสำคัญ

คอลัมน์แรกมาจาก T₁,T₂,T₄,T₈
คอลัมน์สองกำจัด h²
คอลัมน์สามกำจัด h⁴
อธิบายกรณี x⁴ ที่ได้คำตอบตรงพอดี

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

ทิศทางการคำนวณสำคัญมาก: จะเติมช่อง \(R_{k,j}\) ได้ ต้องมีช่องซ้าย \(R_{k,j-1}\) และช่องเฉียงซ้ายบน \(R_{k-1,j-1}\) ก่อน

ภาพจำ: แต่ละช่องใช้ “ค่าละเอียดกว่า” ในแถวเดียวกัน ลบกับ “ค่าหยาบกว่า” จากแถวบน แล้วหารด้วยตัวประกอบ \(4^{j-1}-1\)

สูตรที่ใช้

\[R_{3,3}=R_{3,2}+\frac{R_{3,2}-R_{2,2}}{15}\]

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

สูตรเดิมยังใช้เหมือน EP1 แต่ใน EP2 เราเน้นความหมายของดัชนี \(k\) และ \(j\): \(k\) คือระดับการแบ่งช่วงของ Trapezoid ส่วน \(j\) คือระดับการ extrapolate

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

รูปนี้เท่ากับสูตรบวกส่วนแก้ไข แต่บางคนอ่านง่ายกว่า เพราะเห็นว่าให้ weight กับค่าละเอียด \(R_{k,j-1}\) มากกว่าค่าหยาบ \(R_{k-1,j-1}\)

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

ใช้ ∫₀¹ x⁴ dx=1/5 และสร้างตารางครบถึง R₄,₄

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

ใช้ตัวอย่างเดิม \(\int_0^1 \dfrac{1}{1+x^2}\,dx\) ตารางถึง \(n=8\) ได้ดังนี้

\(k\)	\(n\)	\(R_{k,1}\)	\(R_{k,2}\)	\(R_{k,3}\)	\(R_{k,4}\)
1	1	0.750000	-	-	-
2	2	0.775000	0.783333	-	-
3	4	0.782794	0.785392	0.785529	-
4	8	0.784747	0.785398	0.785399	0.785396

ค่าจริงคือ \(0.785398\) จะเห็นว่าคอลัมน์กลาง ๆ เข้าใกล้มาก แต่คอลัมน์ลึกสุดไม่ได้ต้องดีขึ้นทุกครั้งในเลขทศนิยมจริง เพราะการปัดเศษและข้อมูลตั้งต้นมีผล

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

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

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

ข้อควรระวัง

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

อย่าไล่คอลัมน์ลึกอย่างเดียว: ควรดูการลู่เข้าทั้งตามแถวและตามคอลัมน์ ถ้าค่าเริ่มแกว่งหรือไม่เสถียร ให้หยุดและตรวจฟังก์ชัน

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

สรุป

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