Simpson’s 1/3 Rule: ใช้พาราโบลาแทนกราฟ

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

Simpson’s 1/3 Rule ใช้พาราโบลาประมาณกราฟทีละสองช่วง จึงมักแม่นกว่า Trapezoidal Rule อย่างมากเมื่อฟังก์ชันเรียบ

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

แนวคิดสำคัญ

จำนวนช่วงต้องเป็นเลขคู่
น้ำหนักสลับ 1,4,2,4,...,1
มองเป็นการประมาณด้วยพาราโบลา
เปรียบเทียบ error กับวิธีเส้นตรง

Simpson’s 1/3 Rule ขยับจากการใช้เส้นตรงมาใช้พาราโบลา โดยใช้จุด 3 จุดต่อหนึ่งชุด ได้แก่ ปลายซ้าย จุดกลาง และปลายขวา ถ้ากราฟในช่วงนั้นโค้งคล้ายพาราโบลา วิธีนี้จะตามรูปทรงได้ดีกว่า Trapezoidal Rule อย่างเห็นได้ชัด

ชื่อ 1/3 มาจากตัวคูณ \(\dfrac{h}{3}\) ในสูตร composite ไม่ได้แปลว่าใช้หนึ่งในสามของช่วง วิธีนี้ต้องใช้จำนวนช่วงย่อยเป็นเลขคู่ เพราะหนึ่งชุดของพาราโบลาต้องกิน 2 ช่วงย่อยเสมอ

ภาพจำ: ทุก 2 ช่วงย่อยใช้จุด 3 จุดสร้างพาราโบลา แล้วหาพื้นที่ใต้พาราโบลานั้นแทนพื้นที่ใต้กราฟจริง

สูตรที่ใช้

\[S_n=\frac{h}{3}\left[f(x_0)+4f(x_1)+2f(x_2)+\cdots+4f(x_{n-1})+f(x_n)\right]\]

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

สำหรับช่วงย่อยกว้างเท่ากัน \(h=\dfrac{b-a}{n}\) และ \(n\) เป็นเลขคู่ สูตร composite คือ

\[S_{1/3}=\frac{h}{3}\left[f(x_0)+4\sum_{i\text{ odd}}f(x_i)+2\sum_{i\text{ even},\,i\ne0,n}f(x_i)+f(x_n)\right]\]

น้ำหนัก 4 เกิดกับจุดกลางของแต่ละพาราโบลา เช่น \(x_1,x_3,x_5\) ส่วนน้ำหนัก 2 เกิดกับจุดภายในที่เป็นรอยต่อระหว่างพาราโบลาสองชุด เช่น \(x_2,x_4\) ส่วนปลายสุดสองด้านมีน้ำหนัก 1

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

ใช้ตัวอย่างเดียวกับ Trapezoidal Rule เพื่อให้เห็นความต่างของความแม่นอย่างเป็นธรรม

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

ประมาณ \(\int_0^2 x^4\,dx\) ด้วย \(n=4\) จะมี \(h=0.5\) และจุด \(0,0.5,1.0,1.5,2.0\) ค่าจริงคือ \(\dfrac{2^5}{5}=6.4\)

จุด	ค่า \(x\)	\(f(x)=x^4\)	น้ำหนัก
\(x_0\)	0	0	1
\(x_1\)	0.5	0.0625	4
\(x_2\)	1.0	1	2
\(x_3\)	1.5	5.0625	4
\(x_4\)	2.0	16	1

\[S_{1/3}=\frac{0.5}{3}\left[0+4(0.0625)+2(1)+4(5.0625)+16\right]=6.416667\]

error คือประมาณ \(0.016667\) แม้ฟังก์ชัน \(x^4\) ไม่ใช่พาราโบลา แต่ Simpson ยังทำได้ดีเพราะใช้ข้อมูลความโค้งผ่านจุดกลาง ถ้าใช้ฟังก์ชันพหุนามดีกรีไม่เกิน 3 วิธีนี้จะให้ค่าตรงพอดีในเลขคณิต exact

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

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

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

ข้อควรระวัง

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

เงื่อนไขสำคัญ: Composite Simpson’s 1/3 ต้องมีจำนวนช่วง \(n\) เป็นเลขคู่ และจุดต้องห่างเท่ากัน

ความผิดพลาดที่เจอบ่อยคือใส่น้ำหนัก 4 และ 2 สลับกัน หรือใช้ \(n\) เป็นจำนวนจุดแทนจำนวนช่วง ถ้ามีจุดข้อมูล 5 จุดที่ห่างเท่ากัน แปลว่ามี 4 ช่วง ไม่ใช่ 5 ช่วง

สำหรับข้อมูลที่ไม่เรียบหรือมีจุดกระโดด Simpson อาจไม่ได้ชนะ Trapezoid เสมอ เพราะพาราโบลาที่ fit ผ่านจุดอาจแกว่งผิดธรรมชาติ ควรดูกราฟหรือแบ่งช่วงบริเวณที่พฤติกรรมเปลี่ยนเร็วเป็นพิเศษ

สรุป

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