1
using System;  
using System.Collections.Generic;  
using System.ComponentModel;  
using System.Data;  
using System.Drawing  ;  
using System.Linq;  
using System.Text;  
using System.Windows.Forms;  
using Meta.Numerics.Matrices;  

 double[,] ptsImgOne = new double[3, 5];  
 double[,] ptsImgTwo = new double[3, 5];  


 ptsImgOne = new double[,] { { 0.41726, 0.94478, 0.33771, 0.11120, 0.24169 }, 
{ 0.04965, 0.49086, 0.90005, 0.78025, 0.40391}, 
{ 0.90271, 0.48925, 0.36924, 0.38973, 0.09645 } };  
ptsImgTwo = new double[,] { { 0.13197, 0.57520, 0.35315, 0.04302, 0.73172 },
 { 0.94205, 0.05977, 0.82119, 0.16899, 0.64774 },
 { 0.95613, 0.23477, 0.01540, 0.64911, 0.45092 } };


            RectangularMatrix q1 = new RectangularMatrix(ptsImgOne);
            RectangularMatrix q2 = new RectangularMatrix(ptsImgTwo);


            RectangularMatrix Qmultiply = new RectangularMatrix(5, 9);
            for (int i = 0; i < q2.ColumnCount; i++)
            {
                Qmultiply[i, 0] = q1[0, i] * q2[0, i]; 
                Qmultiply[i, 1] = q1[1, i] * q2[0, i]; 
                Qmultiply[i, 2] = q1[2, i] * q2[0, i]; 

                Qmultiply[i, 3] = q1[0, i] * q2[1, i];
                Qmultiply[i, 4] = q1[1, i] * q2[1, i];
                Qmultiply[i, 5] = q1[2, i] * q2[1, i];

                Qmultiply[i, 6] = q1[0, i] * q2[2, i];
                Qmultiply[i, 7] = q1[1, i] * q2[2, i];
                Qmultiply[i, 8] = q1[2, i] * q2[2, i]; 
            }


SingularValueDecomposition singValDec=    Qmultiply.SingularValueDecomposition();




        }
    }
}

Meta.Numerics.Matricesを使用して行列を操作し、SVDを計算して
いますが、コンパイラがコードの最後の行に到達するたびに、
「SingularValueDecomposition singValDec = Qmultiply.SingularValueDecomposition();」インデックスが配列の境界を超えて
いるというエラーが表示されます。Qmultiply配列と残りの境界を再確認しました。デバッグしようとしましたが、エラーが解決しません。 誰かが私がこの問題を修正するのを手伝ってもらえますか?

4

1 に答える 1

0

MetaパッケージはSVDの計算にいくつかのバグがあるようです。

同じコードでMathNetを使用しましたが、非常にうまく機能しました。

安全のために、MathNet(Math.Numeric)を使用してください。

于 2012-10-22T07:44:23.350 に答える