การแปลงแชงค์
ในการวิเคราะห์เชิงตัวเลขการแปลง Shanksเป็น วิธี การเร่งความเร็วอนุกรมแบบไม่เชิงเส้น เพื่อเพิ่มอัตราการล convergenceของลำดับวิธีนี้ตั้งชื่อตามDaniel Shanksผู้ค้นพบการแปลงลำดับนี้อีกครั้งในปี 1955 โดย R. Schmidt เป็นผู้คิดค้นและเผยแพร่เป็นครั้งแรกในปี 1941 [ 1 ]
เราสามารถคำนวณได้เพียงไม่กี่พจน์ของการขยายอนุกรมการรบกวนโดยปกติไม่เกินสองหรือสามพจน์ และแทบจะไม่เกินเจ็ดพจน์เลย อนุกรมที่ได้มักจะลู่เข้าอย่างช้าๆ หรือแม้กระทั่งลู่ออก อย่างไรก็ตาม พจน์เพียงไม่กี่พจน์เหล่านั้นกลับมีข้อมูลที่น่าทึ่งมากมาย ซึ่งผู้วิจัยควรพยายามอย่างเต็มที่ที่จะดึงข้อมูลเหล่านั้นออกมา มุมมองนี้ได้รับการนำเสนออย่างน่าเชื่อถือในบทความที่ยอดเยี่ยมโดย Shanks (1955) ซึ่งแสดงตัวอย่างที่น่าทึ่งมากมาย รวมถึงตัวอย่างจากกลศาสตร์ของไหลด้วย
สูตร
สำหรับลำดับอนุกรม
จะต้องพิจารณาต่อไป ก่อนอื่นเรากำหนดนิยามของผลรวมย่อยดังนี้:
และสร้างลำดับใหม่หากอนุกรมลู่เข้าก็จะเข้าใกล้ลิมิตเช่นกัน การแปลง Shanks ของลำดับคือลำดับใหม่ที่กำหนดโดย[ 2 ] [ 3 ]
โดยลำดับนี้มักจะลู่เข้าได้เร็วกว่าลำดับ อาจเพิ่มความเร็วได้อีกโดยการใช้การแปลง Shanks ซ้ำๆ โดยการคำนวณฯลฯ
โปรดทราบว่าการแปลงแบบไม่เชิงเส้นที่ใช้ในการแปลง Shanks นั้นโดยพื้นฐานแล้วเหมือนกับที่ใช้ในกระบวนการเดลต้ากำลังสองของ Aitkenดังนั้นเช่นเดียวกับวิธีของ Aitken นิพจน์ทางขวาสุดในนิยามของ (เช่น) จึงมีเสถียรภาพเชิงตัวเลขมากกว่านิพจน์ทางซ้าย (เช่น) ทั้งวิธีของ Aitken และการแปลง Shanks ทำงานกับลำดับ แต่ลำดับที่การแปลง Shanks ทำงานนั้นมักถูกมองว่าเป็นลำดับของผลรวมย่อย แม้ว่าลำดับใดๆ ก็สามารถมองได้ว่าเป็นลำดับของผลรวมย่อยเช่นกัน
ตัวอย่าง

ตัวอย่างเช่น พิจารณาอนุกรมลู่เข้าช้า[ 3 ]
ซึ่งมีผลรวมที่แน่นอน คือ π ≈ 3.14159265 ผลรวมย่อยมีความแม่นยำเพียงหลักเดียว ในขณะที่ความแม่นยำหกหลักต้องใช้การบวกประมาณ 400,000 พจน์
ในตารางด้านล่างนี้ แสดงผลรวมย่อยการแปลง Shanks บนผลรวมย่อยเหล่านั้น รวมถึงการแปลง Shanks ซ้ำๆสำหรับค่าสูงสุดถึง 12 รูปทางด้านขวาแสดงข้อผิดพลาดสัมบูรณ์สำหรับผลรวมย่อยและผลลัพธ์การแปลง Shanks ซึ่งแสดงให้เห็นถึงความแม่นยำและอัตราการล convergence ที่ดีขึ้นอย่างชัดเจน
| 0 | 4.00000000 | — | — | — |
| 1 | 2.66666667 | 3.16666667 | — | — |
| 2 | 3.46666667 | 3.13333333 | 3.14210526 | — |
| 3 | 2.89523810 | 3.14523810 | 3.14145022 | 3.14159936 |
| 4 | 3.33968254 | 3.13968254 | 3.14164332 | 3.14159086 |
| 5 | 2.97604618 | 3.14271284 | 3.14157129 | 3.14159323 |
| 6 | 3.28373848 | 3.14088134 | 3.14160284 | 3.14159244 |
| 7 | 3.01707182 | 3.14207182 | 3.14158732 | 3.14159274 |
| 8 | 3.25236593 | 3.14125482 | 3.14159566 | 3.14159261 |
| 9 | 3.04183962 | 3.14183962 | 3.14159086 | 3.14159267 |
| 10 | 3.23231581 | 3.14140672 | 3.14159377 | 3.14159264 |
| 11 | 3.05840277 | 3.14173610 | 3.14159192 | 3.14159266 |
| 12 | 3.21840277 | 3.14147969 | 3.14159314 | 3.14159265 |
การแปลง Shanks มีความแม่นยำสองหลักอยู่แล้ว ในขณะที่ผลรวมย่อยดั้งเดิมให้ความแม่นยำเท่ากันที่ ที่น่าทึ่งคือมีความแม่นยำหกหลัก ซึ่งได้มาจากการแปลง Shanks ซ้ำๆ ที่ใช้กับเจ็ดพจน์แรกดังที่กล่าวไว้ก่อนหน้านี้จะได้ความแม่นยำหกหลักก็ต่อเมื่อบวกพจน์ประมาณ 400,000 พจน์เท่านั้น
แรงจูงใจ
การแปลง Shanks ได้รับแรงบันดาลใจจากการสังเกตว่า — สำหรับค่าที่มากขึ้น— ผลรวมย่อยมักจะมีลักษณะโดยประมาณเป็น[ 2 ]
โดยที่ลำดับจะลู่เข้าสู่ผลลัพธ์ของอนุกรมชั่วคราวสำหรับ ดังนั้นสำหรับและผลรวมย่อยที่เกี่ยวข้องคือ:
สมการทั้งสามนี้มีตัวแปรที่ไม่ทราบค่าสามตัว: และแก้หาค่าจะได้[ 2 ]
ในกรณีพิเศษที่ตัวส่วนเท่ากับศูนย์: แล้วสำหรับทุก ๆ
การแปลง Shanks ทั่วไป
การแปลง Shanks อันดับkทั่วไป จะแสดงเป็นอัตราส่วนของ ดีเทอร์มิแนนต์ : [ 4 ]
โดยเป็นวิธีการแก้ปัญหาของแบบจำลองสำหรับพฤติกรรมการลู่เข้าของผลรวมย่อยที่มีการเปลี่ยนแปลงชั่วคราวที่แตกต่างกัน:
แบบจำลองพฤติกรรมการลู่เข้าประกอบด้วยตัวแปรที่ไม่ทราบค่า โดยการประเมินสมการข้างต้นที่องค์ประกอบและแก้หาค่าการแสดงออกข้างต้นสำหรับ การแปลง Shanks อันดับที่ kจะได้ การแปลง Shanks ทั่วไปอันดับแรกเท่ากับการแปลง Shanks แบบธรรมดา:
การแปลง Shanks ทั่วไปมีความเกี่ยวข้องอย่างใกล้ชิดกับตัวประมาณ PadéและตารางPadé [ 4 ]
หมายเหตุ: การคำนวณดีเทอร์มิแนนต์ต้องใช้การดำเนินการทางคณิตศาสตร์หลายอย่าง อย่างไรก็ตามปีเตอร์ วินน์ค้นพบขั้นตอนการประเมินแบบเรียกซ้ำที่เรียกว่าอัลกอริธึมเอปซิลอน ซึ่งหลีกเลี่ยงการคำนวณดีเทอร์มิแนนต์[ 5 ] [ 6 ]
ดูเพิ่มเติม
- กระบวนการเดลต้ากำลังสองของไอท์เคน
- การเร่งความเร็วของแอนเดอร์สัน
- อัตราการบรรจบกัน
- การคาดการณ์ของริชาร์ดสัน
- การแปลงลำดับ