ถ้าฉันสร้างเมทริกซ์สมมาตรแบบสุ่มโอกาสที่จะเป็นบวกแน่นอนคืออะไร


32

ฉันมีคำถามแปลก ๆ เมื่อฉันทดลองการเพิ่มประสิทธิภาพของนูน คำถามคือ:

สมมติว่าฉันสุ่ม (พูดการแจกแจงแบบปกติมาตรฐาน) สร้างเมทริกซ์สมมาตร (ตัวอย่างเช่นฉันสร้างเมทริกซ์สามเหลี่ยมด้านบนและเติมครึ่งล่างเพื่อให้แน่ใจว่ามันสมมาตร) โอกาสที่จะเป็นบวกแน่นอน เมทริกซ์? อย่างไรก็ตามมีการคำนวณความน่าจะเป็นหรือไม่?N×N


1
ลองจำลอง ...
kjetil b halvorsen

1
@kjetilbhalvorsen ขอบคุณ แต่ฉันสงสัยว่าโอกาสของค่าลักษณะเฉพาะทั้งหมดมีค่ามากกว่า 0 หรือเราอาจจะทำการวิเคราะห์ด้วยซ้ำ
Haitao Du

6
คำตอบขึ้นอยู่กับวิธีที่คุณสร้างเมทริกซ์ ตัวอย่างเช่นวิธีหนึ่งสร้างค่าลักษณะเฉพาะจริงตามการแจกแจงบางส่วนแล้วผันเมทริกซ์ทแยงมุมนั้นโดยเมทริกซ์มุมฉากแบบสุ่ม ผลที่ได้จะเป็นบวกแน่นอนถ้าหากค่าลักษณะเฉพาะทั้งหมดนั้นเป็นค่าบวก หากคุณกำลังจะสร้างลักษณะเฉพาะอิสระตามการกระจายสมมาตรเกี่ยวกับศูนย์แล้วว่าโอกาสที่เห็นได้ชัดคือที่มากที่สุดn} เพื่อสร้าง PD matrix จากนั้นเลือกค่าเฉพาะของคุณ! (สำหรับการทำงานที่รวดเร็วฉันสร้างเมทริกซ์เช่นโคn2n
วาเรียส

11
ไม่ใช่คำตอบสำหรับคำถามที่ถาม แต่โปรดทราบว่าหากจำลองเมทริกซ์เป็นครั้งแรกกับแต่ละรายการ iid ปกติและขนาดเดียวกันของดังนั้นเป็นสมมาตรและแน่นอนแน่นอนที่มีความน่าจะเป็น 1LNN=LLT
Cliff AB

คำตอบ:


41

หากเมทริกซ์ของคุณถูกดึงมาจากรายการ iid มาตรฐานปกติความน่าจะเป็นที่จะเป็นบวกแน่นอนคือประมาณดังนั้นตัวอย่างเช่นถ้าโอกาสคือ 1/1000 และลงไปอย่างรวดเร็วหลังจากนั้น คุณสามารถค้นหาขยายการสนทนาของคำถามนี้ที่นี่pN3N2/4N=5

คุณค่อนข้างสามารถตอบคำถามนี้ได้ด้วยการยอมรับว่าการกระจายค่าลักษณะเฉพาะของเมทริกซ์ของคุณจะเป็นครึ่งวงกลม Wignerซึ่งมีความสมมาตรประมาณศูนย์ หากค่าลักษณะเฉพาะเป็นอิสระทั้งหมดคุณจะมีโอกาสในเชิงบวกโดยตรรกะนี้ ในความเป็นจริงคุณจะได้รับพฤติกรรมทั้งเนื่องจากความสัมพันธ์ระหว่างค่าลักษณะเฉพาะและกฎหมายที่ควบคุมค่าเบี่ยงเบนขนาดใหญ่ของค่าลักษณะเฉพาะโดยเฉพาะขนาดที่เล็กที่สุดและใหญ่ที่สุด โดยเฉพาะค่าลักษณะเฉพาะแบบสุ่มคล้ายมากที่อนุภาคที่มีประจุและไม่ชอบที่จะได้อยู่ใกล้ชิดกับแต่ละอื่น ๆ พวกเขาจึงขับไล่แต่ละอื่น ๆ (แปลกพอมีข้อมูลที่อาจเกิดขึ้นเช่นเดียวกับอนุภาคที่มีประจุ,ที่(1/2)NN21/rrคือระยะห่างระหว่างค่าลักษณะเฉพาะที่อยู่ติดกัน) ขอให้พวกเขาทั้งหมดเป็นบวกจึงจะขอสูงมาก

นอกจากนี้เนื่องจากกฎความเป็นสากลในทฤษฎีเมทริกซ์แบบสุ่มฉันสงสัยอย่างยิ่งน่าจะเป็นที่สูงกว่าน่าจะเหมือนกันสำหรับเมทริกซ์แบบสุ่มที่ "สมเหตุสมผล" ใด ๆ โดยมีรายการ iid ที่มีค่าเฉลี่ยและส่วนเบี่ยงเบนมาตรฐานpN


5
มันดีที่รู้ว่ามันต่ำมาก ดังนั้นฉันจะไม่ใช้การสุ่มตัวอย่างการปฏิเสธเพื่อสร้างเมทริกซ์ SPD ในอนาคต
Haitao Du

5
@ hxd1011: หากคุณพยายามที่จะลองตัวอย่างเมทริกซ์ SPD ฉันขอแนะนำวิธีที่ฉัน desribed ในความคิดเห็นข้างต้น นอกจากนี้อาจเป็นประโยชน์ในการอ่านข้อมูลเกี่ยวกับการสลายตัวของ Cholesky
Cliff AB

@CliffAB ขอบคุณ ฉันมักจะสร้างเมทริกซ์ความแปรปรวนร่วมแบบเมทริกซ์ SPD ของข้อมูลบางส่วนหรือจากคล้ายกับที่คุณแนะนำ ฉันมีเวลาของการพยายามใส่ตัวเลขบางส่วนลงในเมทริกซ์ขนาดเล็กด้วยตนเองว่าและหวังว่ามันจะเป็นเมทริกซ์ PD AA2×2
Haitao Du
โดยการใช้ไซต์ของเรา หมายความว่าคุณได้อ่านและทำความเข้าใจนโยบายคุกกี้และนโยบายความเป็นส่วนตัวของเราแล้ว
Licensed under cc by-sa 3.0 with attribution required.