मैट्रिक्स इनवर्स का उपयोग MathNet.Numerics

Question

मैं कुछ सी # कोड लिख रहा हूं जो MathNet.Numerics.LinearAlgebra का उपयोग करता है, और पाठ्यपुस्तक उदाहरण के साथ परिणामों से मिलान करने का प्रयास करता है। कोड का एक हिस्सा एक और सरणी "Zbus" में एक complex32 सरणी "Ybus", और यह स्टोर के एक व्युत्क्रम करता है:

Matrix<Complex32> Ybus = Matrix<Complex32>.Build.Dense(numBuses, numBuses); 

Matrix<Complex32> Zbus = Matrix<Complex32>.Build.Dense(numBuses, numBuses); 

मेरे Ybus किताब में बिल्कुल उदाहरण मेल खाता है।

Ybus = j[ -13  5  4  0 
      5 -13.5 2.5  2 
      4  2.5  -9 2.5 
      0  2 2.5 -4.5] 

लेकिन जब मैं एक व्युत्क्रम करना

Zbus = Ybus.Inverse(); 

Zbus के परिणाम सभी NaN

जबकि किताब से सही परिणाम इस तरह दिखता है:

Zbus = j[ .15  .09 .12 .11 
      .09  .15 .12 .13 
      .12  .12 .25 .19 
      .11  .13 .19 .39] 

किसी को भी क्या कोई विचार है कि समस्या क्या हो सकती है? शायद एक जटिल मैट्रिक्स के उलटा कुछ मुद्दों है?

सबक सीखा:

: इतना है कि वे 0 के या उलटा को उड़ाने होगा :) ... की पंक्तियों अपने सरणियों बहुत बड़ा नहीं बनाते हैं यहाँ सही जवाब है

Answer 1

जैसा कि जेसन ने उल्लेख किया है, ऐसा लगता है कि यह ठीक काम करता है। उदाहरण के लिए:

var y = Complex32.ImaginaryOne * CreateMatrix.Dense(4, 4, new Complex32[] {-13f,5f,4f,0f,5f,-13.5f,2.5f,2f,4f,2.5f,-9f,2.5f,0f,2f,2.5f,-4.5f}); 
y.ToString("F3"); 
y.Inverse().ToString("F3"); 

(पुस्तक में बुरा राउंडिंग को छोड़ कर) अपनी पुस्तक परिणाम मिलान, निम्न उत्पादन प्रदान करता है:

DenseMatrix 4x4-Complex32 
(0.000, -13.000) (0.000, 5.000) (0.000, 4.000) (0.000, 0.000) 
    (0.000, 5.000) (0.000, -13.500) (0.000, 2.500) (0.000, 2.000) 
    (0.000, 4.000) (0.000, 2.500) (0.000, -9.000) (0.000, 2.500) 
    (0.000, 0.000) (0.000, 2.000) (0.000, 2.500) (0.000, -4.500) 

DenseMatrix 4x4-Complex32 
(0.000, 0.153) (0.000, 0.097) (0.000, 0.126) (0.000, 0.113) 
(0.000, 0.097) (0.000, 0.153) (0.000, 0.124) (0.000, 0.137) 
(0.000, 0.126) (0.000, 0.124) (0.000, 0.256) (0.000, 0.197) 
(0.000, 0.113) (0.000, 0.137) (0.000, 0.197) (0.000, 0.393)