private static int[,] GetGaussMatrix(int[,] parityMatrix, int gap)
{
int[,] resultMatrix = new int[parityMatrix.GetLength(0), parityMatrix.GetLength(0)];
int[,] EMatrix = new int[gap, parityMatrix.GetLength(0) - gap];
int[,] TMatrix = new int[parityMatrix.GetLength(0) - gap, parityMatrix.GetLength(0) - gap];
int matrixRow = 0;
for (int row = 0; row < parityMatrix.GetLength(0); row++)
{
int matrixColumn = 0;
for (int column = parityMatrix.GetLength(1) - TMatrix.GetLength(1); column < parityMatrix.GetLength(1); column++)
{
if (row < parityMatrix.GetLength(0) - gap)
{
TMatrix[row, matrixColumn] = parityMatrix[row, column];
}
else
{
EMatrix[matrixRow, matrixColumn] = parityMatrix[row, column];
}
matrixColumn++;
}
if (row >= parityMatrix.GetLength(0) - gap)
{
matrixRow++;
}
}
int[,] ETMatrix = MatrixOperation.LogicMultiplicateMatrixes(EMatrix, MatrixOperation.LogicGetInverseMatrix(TMatrix));
matrixRow = 0;
for (int row = 0; row < resultMatrix.GetLength(0); row++)
{
for (int column = 0; column < resultMatrix.GetLength(1); column++)
{
if (row < (resultMatrix.GetLength(0) - gap))
{
if (row == column)
{
resultMatrix[row, column] = 1;
}
}
else
{
if (column < resultMatrix.GetLength(0) - gap)
{
resultMatrix[row, column] = ETMatrix[matrixRow, column];
}
else
{
if (row == column)
{
resultMatrix[row, column] = 1;
}
}
}
}
if (row >= (resultMatrix.GetLength(0) - gap))
{
matrixRow++;
}
}
return(resultMatrix);
}
private int[] EncodeMessage(int[,] parityMatrix, int[] message, int gap, int oldGap)
{
int[,] p1 = null;
int[,] messageArray = new int[1, message.Length];
for (int i = 0; i < message.Length; i++)
{
messageArray[0, i] = message[i];
}
if (gap != 0)
{
int startRow = parityMatrix.GetLength(0) - gap;
int[,] DMatrix = new int[gap, gap];
int[,] CMatrix = new int[gap, parityMatrix.GetLength(1) - DMatrix.GetLength(1) - parityMatrix.GetLength(0) + gap];
for (int i = startRow; i < parityMatrix.GetLength(0); i++)
{
for (int j = 0; j < parityMatrix.GetLength(1) - parityMatrix.GetLength(0) + gap; j++)
{
if (j < CMatrix.GetLength(1))
{
CMatrix[i % startRow, j] = parityMatrix[i, j];
}
else
{
DMatrix[i % startRow, j % CMatrix.GetLength(1)] = parityMatrix[i, j];
}
}
}
int[,] DCMatrix = MatrixOperation.LogicMultiplicateMatrixes(MatrixOperation.LogicGetInverseMatrix(DMatrix), CMatrix);
p1 = MatrixOperation.Transponse(MatrixOperation.LogicMultiplicateMatrixes(DCMatrix, MatrixOperation.Transponse(messageArray)));
}
else
{
p1 = new int[1, oldGap];
}
int[,] BMatrix = new int[parityMatrix.GetLength(0) - gap, oldGap];
int[,] TMatrix = new int[parityMatrix.GetLength(0) - gap, parityMatrix.GetLength(0) - gap];
int[,] AMatrix = new int[BMatrix.GetLength(0), parityMatrix.GetLength(1) - BMatrix.GetLength(1) - TMatrix.GetLength(1)];
for (int i = 0; i < AMatrix.GetLength(0); i++)
{
for (int j = 0; j < parityMatrix.GetLength(1); j++)
{
if (j < AMatrix.GetLength(1))
{
AMatrix[i, j] = parityMatrix[i, j];
}
else if (j < AMatrix.GetLength(1) + BMatrix.GetLength(1))
{
BMatrix[i, j % AMatrix.GetLength(1)] = parityMatrix[i, j];
}
else
{
TMatrix[i, j % (AMatrix.GetLength(1) + BMatrix.GetLength(1))] = parityMatrix[i, j];
}
}
}
int[,] AuMatrix = MatrixOperation.LogicMultiplicateMatrixes(AMatrix, MatrixOperation.Transponse(messageArray));
int[,] addedMatrixes = null;
int[,] Bp1Matrix = MatrixOperation.LogicMultiplicateMatrixes(BMatrix, MatrixOperation.Transponse(p1));
addedMatrixes = MatrixOperation.LogicAddUpMatrixes(AuMatrix, Bp1Matrix);
int[,] p2 = MatrixOperation.Transponse(MatrixOperation.LogicMultiplicateMatrixes(MatrixOperation.LogicGetInverseMatrix(TMatrix), addedMatrixes));
int[] codeword = new int[message.Length + p1.GetLength(1) + p2.GetLength(1)];
for (int i = 0; i < codeword.Length; i++)
{
if (i < message.Length)
{
codeword[i] = message[i];
}
else if (i < message.Length + p1.GetLength(1))
{
codeword[i] = p1[0, i % message.Length];
}
else
{
codeword[i] = p2[0, i % (message.Length + p1.GetLength(1))];
}
}
return(codeword);
}
private static int[,] GetEncodeMatrix(int[,] parityMatrix, out int gap)
{
int[,] parityMatrixCopy = (int[, ])parityMatrix.Clone();
int permutations = 0;
for (int i = 0; i < parityMatrixCopy.GetLength(0) - permutations; i++)
{
if (parityMatrixCopy[i, parityMatrixCopy.GetLength(1) - 1] == 1)
{
MatrixOperation.PermuteRow(parityMatrixCopy, i, parityMatrixCopy.GetLength(0) - 1);
permutations++;
i = 0;
}
}
gap = permutations - 1;
// номер столбца, с которым мы работаем
int p = parityMatrixCopy.GetLength(1) - 2;
int j = 2;
for (; ;)
{
int diagonalRow = parityMatrixCopy.GetLength(0) - (gap + j);
if (diagonalRow < 0)
{
break;
}
int columnWithMinOnes = 0;
int minOnesInColumn = Int32.MaxValue;
for (int column = 0; column <= p; column++)
{
int onesInColumn = 0;
for (int row = 0; row <= diagonalRow; row++)
{
if (parityMatrixCopy[row, column] == 1)
{
onesInColumn++;
}
}
if (onesInColumn < minOnesInColumn && onesInColumn > 0)
{
minOnesInColumn = onesInColumn;
columnWithMinOnes = column;
}
}
MatrixOperation.PermuteColumn(parityMatrixCopy, columnWithMinOnes, p);
int rowWithOne = 0;
for (int row = 0; row <= diagonalRow; row++)
{
if (parityMatrixCopy[row, p] == 1)
{
rowWithOne = row;
break;
}
}
MatrixOperation.PermuteRow(parityMatrixCopy, rowWithOne, diagonalRow);
if (minOnesInColumn > 1)
{
for (int i = 0; i < minOnesInColumn - 1; i++)
{
rowWithOne = 0;
for (int row = 0; row < diagonalRow; row++)
{
if (parityMatrixCopy[row, p] == 1)
{
rowWithOne = row;
break;
}
}
MatrixOperation.PermuteRow(parityMatrixCopy, rowWithOne, parityMatrixCopy.GetLength(0) - 1);
gap++;
diagonalRow--;
}
}
j++;
p--;
}
int[,] gaussMatrix = GetGaussMatrix(parityMatrixCopy, gap);
parityMatrixCopy = MatrixOperation.LogicMultiplicateMatrixes(gaussMatrix, parityMatrixCopy);
Console.WriteLine("Проверочная");
MatrixOperation.ShowMatrix(parityMatrixCopy);
return(parityMatrixCopy);
}
