การคำนวณพหุนามลักษณะของเมทริกซ์กระจัดกระจายจริง


9

รับเมทริกซ์กระจัดกระจายทั่วไป ARn×nด้วยm << n (การแก้ไข:ม.«n2) องค์ประกอบที่ไม่ใช่ศูนย์ (โดยทั่วไป ม.O(n)) A เป็นเรื่องทั่วไปในแง่ที่ว่ามันไม่มีคุณสมบัติที่เฉพาะเจาะจง (เช่นความชัดเจนเชิงบวก) และไม่มีโครงสร้าง (เช่นความเป็นแถบสี)

อะไรคือวิธีการเชิงตัวเลขที่ดีในการคำนวณทั้งพหุนามลักษณะหรือพหุนามน้อยที่สุดของA?


3
ดูเหมือนว่าคุณต้องการคำนวณค่าลักษณะเฉพาะทั้งหมด ทำไมคุณถึงต้องการพหุนามและคุณต้องการให้มันแสดงออกมาอย่างไร? พื้นฐาน monomial นั้นมีสภาพที่ไม่ดีดังนั้นค่าสัมประสิทธิ์จึงไม่สามารถคำนวณได้อย่างเสถียรในคณิตศาสตร์ที่มีความแม่นยำแน่นอน
Jed Brown

@JedBrown ใคร่ครวญมากขึ้น ในคำตอบของคำถามนี้ฉันได้ให้วิธีพีชคณิตในการแปลงเมทริกซ์ซึ่งเป็นที่รู้จักกันดีในพีชคณิตคอมพิวเตอร์ ฉันต้องการทราบว่าฉันสามารถใช้เป็นเมทริกซ์เชิงตัวเลขได้หรือไม่ โปรดทราบว่าเพื่อจุดประสงค์ของคำถามนี้ฉันสนใจวิธีการเชิงตัวเลขสำหรับการค้นหาลักษณะ / พหุนามน้อยที่สุดแทนที่จะเป็นอินเวอร์ส

คำตอบ:


1

ถ้า O(n3)ความซับซ้อนไม่ได้เป็นตัวหยุดแล้วคุณอาจต้องการดูวิธี Danilevskii เป็นที่รู้จักกันดีในวรรณคดีรัสเซียเกี่ยวกับพีชคณิตเชิงเส้นเชิงตัวเลข แต่มีข้อมูลไม่มากนักในภาษาอังกฤษ คุณสามารถเริ่มต้นจากลิงค์นี้

แนวคิดนี้ค่อนข้างตรงไปตรงมา: เมทริกซ์จะค่อยๆลดลงเป็นรูปแบบปกติของ Frobeniusโดยการแปลงความคล้ายคลึงกันของ "Gaussian elimination-like" หากคุณไม่พบข้อมูลฉันสามารถทำให้อัลกอริทึมซับซ้อนยิ่งขึ้น


1

คุณสามารถใช้วิธีตัวเลขเช่น QR Factorization หรือ Power Method และ realtives (inverse power เป็นต้น) เพื่อคำนวณค่าลักษณะเฉพาะของเมทริกซ์ทั่วไปของคุณ หลังจากนั้นคุณสามารถคำนวณพหุนามลักษณะของคุณโดยแยกตัวประกอบเป็น: (λ-λ1) (λ-λ2) ... (λ-λn) = 0 โดยที่λiเป็นค่าลักษณะเฉพาะที่คำนวณ นี่คือการนำเสนอสั้น ๆ เกี่ยวกับวิธีการใช้พลังงานและ QR:

QR-Power


0

โดยวิธีการ: คุณต้องการที่จะบอกว่าคุณมี ม.«O(n2)รายการ? ถ้าแน่นอน ม.«O(n) ดังนั้นส่วนใหญ่ของแถวและคอลัมน์จะว่างเปล่าอย่างสมบูรณ์และเป็นไปได้ว่าพหุนามลักษณะในความเป็นจริงไม่ใช่ระดับ n แต่ระดับ O(ม.).


Ops ไม่ฉันตั้งใจจะพูดม.«n2, กล่าวคือ ม.O(n). ขอโทษด้วยกับเรื่องนั้น.
โดยการใช้ไซต์ของเรา หมายความว่าคุณได้อ่านและทำความเข้าใจนโยบายคุกกี้และนโยบายความเป็นส่วนตัวของเราแล้ว
Licensed under cc by-sa 3.0 with attribution required.