Dream是一个数字收音机Digital Radio Mondiale (DRM) 的软件实现。

/******************************************************************************\ * Technische Universitaet Darmstadt, Institut fuer Nachrichtentechnik * Copyright (c) 2001 * * Author(s): * Volker Fischer * * Description: * Multi-level-channel (de)coder (MLC) * ****************************************************************************** * * This program is free software; you can redistribute it and/or modify it under * the terms of the GNU General Public License as published by the Free Software * Foundation; either version 2 of the License, or (at your option) any later * version. * * This program is distributed in the hope that it will be useful, but WITHOUT * ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or FITNESS * FOR A PARTICULAR PURPOSE. See the GNU General Public License for more * details. * * You should have received a copy of the GNU General Public License along with * this program; if not, write to the Free Software Foundation, Inc., * 59 Temple Place, Suite 330, Boston, MA 02111-1307 USA * \******************************************************************************/ #include "MLC.h" /* Implementation *************************************************************/ /******************************************************************************\ * MLC-encoder * \******************************************************************************/ void CMLCEncoder::ProcessDataInternal(CParameter& Parameter) { int i, j; int iNoBitsOutDec; int iElementCounter; /* Energy dispersal ----------------------------------------------------- */ /* VSPP is treated as a separate part for energy dispersal */ EnergyDisp.ProcessData(pvecInputData); /* Partitioning of input-stream ----------------------------------------- */ iElementCounter = 0; if (iL[2] == 0) { /* Standard departitioning */ /* Protection level A */ for (j = 0; j < iLevels; j++) { /* Bits */ for (i = 0; i < iM[j][0]; i++) { vecbiEncInBuffer[j][i] = (*pvecInputData)[iElementCounter]; iElementCounter++; } } /* Protection level B */ for (j = 0; j < iLevels; j++) { /* Bits */ for (i = 0; i < iM[j][1]; i++) { vecbiEncInBuffer[j][iM[j][0] + i] = (*pvecInputData)[iElementCounter]; iElementCounter++; } } } else { /* Special departitioning with hierarchical modulation. First set hierarchical bits at the beginning, then append the rest */ /* Hierarchical frame (always "iM[0][1]"). "iM[0][0]" is always "0" in this case */ for (i = 0; i < iM[0][1]; i++) { vecbiEncInBuffer[0][i] = (*pvecInputData)[iElementCounter]; iElementCounter++; } /* Protection level A (higher protected part) */ for (j = 1; j < iLevels; j++) { /* Bits */ for (i = 0; i < iM[j][0]; i++) { vecbiEncInBuffer[j][i] = (*pvecInputData)[iElementCounter]; iElementCounter++; } } /* Protection level B (lower protected part) */ for (j = 1; j < iLevels; j++) { /* Bits */ for (i = 0; i < iM[j][1]; i++) { vecbiEncInBuffer[j][iM[j][0] + i] = (*pvecInputData)[iElementCounter]; iElementCounter++; } } } /* Convolutional encoder ------------------------------------------------ */ for (j = 0; j < iLevels; j++) { iNoBitsOutDec = ConvEncoder[j].Encode(vecbiEncInBuffer[j], vecbiEncOutBuffer[j]); #ifdef _DEBUG_ if (iNoBitsOutDec != iNoEncBits) { DebugError("MLC decoded bits test", "No bits decoder out", iNoBitsOutDec, "No encoded bits", iNoEncBits); } #endif } /* Bit interleaver ------------------------------------------------------ */ for (j = 0; j < iLevels; j++) if (piInterlSequ[j] != -1) BitInterleaver[piInterlSequ[j]].Interleave(vecbiEncOutBuffer[j]); /* QAM mapping ---------------------------------------------------------- */ QAMMapping.Map(vecbiEncOutBuffer[0], vecbiEncOutBuffer[1], vecbiEncOutBuffer[2], vecbiEncOutBuffer[3], vecbiEncOutBuffer[4], vecbiEncOutBuffer[5], pvecOutputData); } void CMLCEncoder::InitInternal(CParameter& TransmParam) { int i; CalculateParam(TransmParam, eChannelType); iNoInBits = iL[0] + iL[1] + iL[2]; /* Init modules --------------------------------------------------------- */ /* Energy dispersal */ EnergyDisp.Init(iNoInBits, iL[2]); /* Encoder */ for (i = 0; i < iLevels; i++) ConvEncoder[i].Init(eCodingScheme, eChannelType, iN[0], iN[1], iM[i][0], iM[i][1], iCodeRate[i][0], iCodeRate[i][1], i); /* Bit interleaver */ /* First init all possible interleaver (According table "TableMLC.h" -> "Interleaver sequence") */ if (eCodingScheme == CParameter::CS_3_HMMIX) { BitInterleaver[0].Init(iN[0], iN[1], 13); BitInterleaver[1].Init(iN[0], iN[1], 21); } else { BitInterleaver[0].Init(2 * iN[0], 2 * iN[1], 13); BitInterleaver[1].Init(2 * iN[0], 2 * iN[1], 21); } /* QAM-mapping */ QAMMapping.Init(iN_mux, eCodingScheme); /* Allocate memory for internal bit-buffers ----------------------------- */ for (i = 0; i < iLevels; i++) { /* Buffers for each encoder on all different levels */ /* Add bits from higher protected and lower protected part */ vecbiEncInBuffer[i].Init(iM[i][0] + iM[i][1]); /* Encoder output buffers for all levels. Must have the same length */ vecbiEncOutBuffer[i].Init(iNoEncBits); } /* Define block-size for input and output */ iInputBlockSize = iNoInBits; iOutputBlockSize = iN_mux; } /******************************************************************************\ * MLC-decoder * \******************************************************************************/ void CMLCDecoder::ProcessDataInternal(CParameter& ReceiverParam) { int i, j, k; int iElementCounter; _BOOLEAN bIteration; /* Save input signal for signal constellation. We cannot use the copy operator of vector because the input vector is not of the same size as our intermediate buffer, therefore the "for"-loop */ for (i = 0; i < iInputBlockSize; i++) vecSigSpacBuf[i] = (*pvecInputData)[i].cSig; #if 0 // TEST static FILE* pFile = fopen("test/constellation.dat", "w"); if (eChannelType == CParameter::CT_MSC) { for (i = 0; i < iInputBlockSize; i++) fprintf(pFile, "%e %e\n", vecSigSpacBuf[i].real(), vecSigSpacBuf[i].imag()); fflush(pFile); } // close all;load constellation.dat;constellation=complex(constellation(:,1),constellation(:,2));plot(constellation,'.') #endif /* Iteration loop */ for (k = 0; k < iNoIterations + 1; k++) { for (j = 0; j < iLevels; j++) { /* Metric ------------------------------------------------------- */ if (k > 0) bIteration = TRUE; else bIteration = FALSE; MLCMetric.CalculateMetric(pvecInputData, vecMetric, vecbiSubsetDef[0], vecbiSubsetDef[1], vecbiSubsetDef[2], vecbiSubsetDef[3], vecbiSubsetDef[4], vecbiSubsetDef[5], j, bIteration); /* Bit deinterleaver -------------------------------------------- */ if (piInterlSequ[j] != -1) BitDeinterleaver[piInterlSequ[j]].Deinterleave(vecMetric); /* Viterbi decoder ---------------------------------------------- */ rAccMetric = ViterbiDecoder[j].Decode(vecMetric, vecbiDecOutBits[j]); /* The last branch of encoding and interleaving must not be used at the very last loop */ /* "iLevels - 1" for iLevels = 1, 2, 3 "iLevels - 2" for iLevels = 6 */ if ((k < iNoIterations) || ((k == iNoIterations) && !(j >= iIndexLastBranch))) { /* Convolutional encoder ------------------------------------ */ ConvEncoder[j].Encode(vecbiDecOutBits[j], vecbiSubsetDef[j]); /* Bit interleaver ------------------------------------------ */ if (piInterlSequ[j] != -1) BitInterleaver[piInterlSequ[j]]. Interleave(vecbiSubsetDef[j]); } } } /* De-partitioning of input-stream ---------------------------------