IMath.Msvd

Data. Входной параметр, двумерный массив размерности n×m;

U. Выходной параметр, двумерный массив размерности n×n;

S. Выходной параметр, одномерный массив размерности n;

V. Выходной параметр, двумерный массив размерности m×n.

Примечание. Параметр S является вектором и содержит элементы главной диагонали матрицы S, получаемой через сингулярное разложение.

Описание

Метод Msvd возвращает результат сингулярного разложения вещественной матрицы.

Комментарии

Сингулярное разложение (англ. singular value decomposition, SVD) - разложение прямоугольной вещественной матрицы на ортогональные матрицы U и V и диагональную матрицу S такие, что произведение U·S·V* дает исходную матрицу, где V* - транспонированная матрица к V.

Примечание. В «Форсайт. Аналитическая платформа» реализован экономный вариант SDV-разложения, из-за чего может не выполняться условие ортогональности матриц U и V: матрицы U и V могут получаться прямоугольными, а не квадратными, и произведения U·U* и V·V* могут не давать единичную матрицу. Единичную матрицу в таком случае будут давать только произведения U*·U и V*·V.

В общем случае соотношение Data = U·S·V* не имеет единственного решения. Метод Msvd рассчитывает один из вариантов.

Пример

Для выполнения примера добавьте ссылку на системную сборку MathFin.

Sub UserProc;
Var
    dimcol, dimrow: Integer;
    matr: Array Of Double;
    i, j : Integer;
    U, S, V : Array Of Double;
    str : String;
Begin
    dimrow := 3;
    dimcol := 5;
    // Входная матрица:
    matr := New Double[dimrow, dimcol];
    matr[0,0] := 1;     matr[0,1] := 2;     matr[0,2] := 0.2;   matr[0,3] := 3;     matr[0,4] := 8;
    matr[1,0] := 0.5;   matr[1,1] := -1;    matr[1,2] := 2.2;   matr[1,3] := 1.4;   matr[1,4] := 2.1;
    matr[2,0] := 0.7;   matr[2,1] := 12;    matr[2,2] := 10;    matr[2,3] := 17;    matr[2,4] := 0;
    // Сингулярное разложение матрицы:
    math.Msvd(matr, U, S, V);
    Debug.WriteLine("U = ");
    For i:=0 To U.GetUpperBound(1) Do
        str := "";
        For j:=0 To U.GetUpperBound(2) Do
            str := str + U[i,j].ToString + "; ";
        End For;
        Debug.WriteLine(str);
    End For;
    Debug.WriteLine("S = " );
    str := "";
    For i:=0 To S.GetUpperBound(1) Do
        str := str + S[i].ToString + "; ";
    End For;
    Debug.WriteLine(str);
    Debug.WriteLine("V = " );
    For i:=0 To V.GetUpperBound(1) Do
        str := "";
        For j:=0 To V.GetUpperBound(2) Do
            str := str + V[i,j].ToString + "; ";
        End For;
        Debug.WriteLine(str);
    End For;
End Sub UserProc;

После выполнения примера в окно консоли будет выведено SVD-разложение исходной матрицы.

IMath.Msvd

Синтаксис

Параметры

Описание

Комментарии

Пример