บทนำ
ในการคำนวณรากของอนุกรม Hermite สามารถใช้เมธอด hermite.hermroots() จาก NumPy เมธอดนี้จะคืนค่าเป็นอาเรย์ของรากของอนุกรม หากรากทั้งหมดเป็นจำนวนจริง ผลลัพธ์จะเป็นจำนวนจริง มิฉะนั้นจะเป็นจำนวนเชิงซ้อน พารามิเตอร์ c คืออาเรย์ 1 มิติของสัมประสิทธิ์
การประมาณค่ารากได้มาจากค่าเฉพาะของเมทริกซ์คู่หู รากที่อยู่ไกลจากจุดกำเนิดของระนาบเชิงซ้อนอาจมีข้อผิดพลาดมาก เนื่องจากความไม่เสถียรทางตัวเลข สำหรับรากที่มีความหลายมากกว่า 1 ก็จะมีข้อผิดพลาดมากขึ้นเช่นกัน รากที่แยกจากกันและอยู่ใกล้จุดกำเนิดสามารถปรับปรุงได้ด้วยการวนซ้ำวิธีของนิวตัน ไปสักไม่กี่รอบ
ขั้นตอนการทำงาน
1. นำเข้าไลบรารี
from numpy.polynomial import hermite as H
2. คำนวณรากของอนุกรม Hermite
print("Result...\n", H.hermroots((-1, 0, 1)))
3. ตรวจสอบชนิดข้อมูล
print("\nType...\n", H.hermroots((-1, 0, 1)).dtype)
4. ตรวจสอบรูปร่าง
print("\nShape...\n", H.hermroots((-1, 0, 1)).shape)
ตัวอย่างโค้ดสมบูรณ์
from numpy.polynomial import hermite as H
# คำนวณรากของอนุกรม Hermite
print("Result...\n", H.hermroots((-1, 0, 1)))
# ตรวจสอบชนิดข้อมูล
print("\nType...\n", H.hermroots((-1, 0, 1)).dtype)
# ตรวจสอบรูปร่าง
print("\nShape...\n", H.hermroots((-1, 0, 1)).shape)
ผลลัพธ์ที่คาดหวัง
Result... [-0.8660254 0.8660254] Type... float64 Shape... (2,)
