User Manual, Developers Guide and API Documentation

SNR2MI.hpp

Go to the documentation of this file.

/*******************************************************************************
 * This file is part of openWNS (open Wireless Network Simulator)
 * _____________________________________________________________________________
 *
 * Copyright (C) 2004-2007
 * Chair of Communication Networks (ComNets)
 * Kopernikusstr. 5, D-52074 Aachen, Germany
 * phone: ++49-241-80-27910,
 * fax: ++49-241-80-22242
 * email: info@openwns.org
 * www: http://www.openwns.org
 * _____________________________________________________________________________
 *
 * openWNS is free software; you can redistribute it and/or modify it under the
 * terms of the GNU Lesser General Public License version 2 as published by the
 * Free Software Foundation;
 *
 * openWNS 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 Lesser General Public License for more
 * details.
 *
 * You should have received a copy of the GNU Lesser General Public License
 * along with this program.  If not, see <http://www.gnu.org/licenses/>.
 *
 ******************************************************************************/

#ifndef RISE_PLMAPPING_SNR2MI_HPP
#define RISE_PLMAPPING_SNR2MI_HPP

#include <WNS/PowerRatio.hpp>
#include <WNS/container/Mapping.hpp>
#include <WNS/Enum.hpp>
#include <WNS/Singleton.hpp>
#include <WNS/service/phy/phymode/SNR2MIInterface.hpp>

namespace rise { namespace plmapping { 

    namespace tests {
        class SNR2MITest;
        class CoderFullTests;
    }

    const unsigned int MAX_MODULATIONS_FOR_TABLE = 6; // only up to QAM64

    class SNR2MI :
        virtual public wns::service::phy::phymode::SNR2MIInterface
    {
        friend class tests::SNR2MITest;
        friend class tests::CoderFullTests;
    public:
        virtual ~SNR2MI(){}
        virtual double convertSNR2MIB(const wns::Ratio& snr, const wns::service::phy::phymode::PhyModeInterface& phyMode) const;
        virtual double convertSNR2MI(const wns::Ratio& snr, const wns::service::phy::phymode::PhyModeInterface& phyMode) const;
        virtual wns::Ratio convertMIB2SNR(const double& mib, const wns::service::phy::phymode::Modulation& mod) const;

    private:
        virtual double mapSNR2MI(const wns::Ratio& snr, wns::service::phy::phymode::Modulation modTYPE) const = 0;
        double convertSNR2MIB(const wns::Ratio& snr, wns::service::phy::phymode::Modulation modTYPE) const;

        static double BER2MIB(double ber);
        static double MIB2BER(double mib);
    };


    class SNR2MIFormula :
        public SNR2MI
    {
    public:
        SNR2MIFormula(){};
        virtual ~SNR2MIFormula(){};

    private:
        virtual double mapSNR2MI(const wns::Ratio& snr, wns::service::phy::phymode::Modulation modTYPE) const;
    };


    class SNR2MIMap :
        public SNR2MI,
        public wns::container::Mapping<wns::service::phy::phymode::Modulation, double>
    {
    public:
        SNR2MIMap();
        virtual ~SNR2MIMap(){};

    private:
        virtual double mapSNR2MI(const wns::Ratio& snr, wns::service::phy::phymode::Modulation modTYPE) const;

        void initBPSK();
        void initQPSK();
        void initQAM8();
        void initQAM16();
        void initQAM32();
        void initQAM64();

    };

    inline double rootFindByBisection(double (*func)(double), double offset, double x1, double x2, double accuracy)
    {
        const unsigned int maxIterations = 100;
        unsigned int iteration;
        double f    = (*func)(x1) - offset;
        double fmid = (*func)(x2) - offset;
        assure(f*fmid<0.0, "rootFindByBisection(): Root must be bracketed within ["<<x1<<","<<x2<<"]: f="<<f<<", fmid="<<fmid);
        double rtb = (f<0.0) ? x1 : x2;
        double dx  = (f<0.0) ? (x2-x1) : (x1-x2);
        for(iteration=1; iteration<=maxIterations; iteration++) {
            double xmid=rtb+(dx *= 0.5);
            fmid = (*func)(xmid) - offset;
            if (fmid <= 0.0) rtb=xmid;
            if ((fabs(dx) < accuracy) || (fmid == 0.0)) return rtb;
        }
        assure(false,"rootFindByBisection(): too many iterations: "<<maxIterations);
        return rtb;
    }

}}
#endif // RISE_PLMAPPING_SNR2MI_HPP

---