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 bool DecodeMessage(int[] codeword, int[,] parityMatrix, float rank, out int[] originalMessage)
{
int[] finalCodeWord = new int[codeword.Length];
bool isSindromConfirmed = false;
int[,] parityMatrixTransposed = MatrixOperation.Transponse(parityMatrix);
int[] syndrome = MatrixOperation.LogicMultiplicateVectorMatrix(codeword, parityMatrixTransposed);
for (int i = 0; i < syndrome.Length; i++)
{
if (syndrome[i] == 1)
{
isSindromConfirmed = true;
break;
}
}
Console.WriteLine("Синдром = " + isSindromConfirmed);
if (!isSindromConfirmed)
{
originalMessage = new int[(int)(codeword.Length * rank)];
for (int i = 0; i < originalMessage.Length; i++)
{
originalMessage[i] = codeword[i];
}
return(isSindromConfirmed);
}
int[] F = MatrixOperation.MultiplicateVectorMatrix(syndrome, parityMatrix);
int threshold = 0;
int[] codewordCopy = (int[])codeword.Clone();
for (int i = 0; i < F.Length; i++)
{
if (threshold < F[i])
{
threshold = F[i];
}
}
Console.WriteLine("Порог " + threshold);
for (int i = 0; i < F.Length; i++)
{
if (F[i] >= threshold)
{
Console.WriteLine("Исправление в {0} элементе", i);
codewordCopy[i] = (codewordCopy[i] + 1) % 2;
break;
}
}
originalMessage = codewordCopy;
Console.WriteLine("Новое кодовое слово");
for (int i = 0; i < codewordCopy.Length; i++)
{
Console.Write(codeword[i] + " ");
}
Console.WriteLine();
return(isSindromConfirmed);
}
