wssuschannel2.h

Go to the documentation of this file.

/*-*-c++-*-*****************************************************************
                          wssuschannel.h  -  MIMO WSSUS channel classes
                             -------------------
    begin                : May 2003
    copyright            : (C) 2003 by Christian Stimming
    email                : stimming@tuhh.de
 ***************************************************************************/

/***************************************************************************
 *                                                                         *
 *   This library is free software; you can redistribute it and/or         *
 *   modify it under the terms of the GNU Lesser General Public            *
 *   License as published by the Free Software Foundation; either          *
 *   version 2.1 of the License, or (at your option) any later version.    *
 *                                                                         *
 *   This library 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 library; if not, write to the Free Software   *
 *   Foundation, Inc., 59 Temple Place, Suite 330, Boston,                 *
 *   MA  02111-1307  USA                                                   *
 *                                                                         *
 ***************************************************************************/

#ifndef WSSUSCHANNEL2_H
#define WSSUSCHANNEL2_H

#include <iostream>
#include <fstream>
#include <complex>

#include <simthetic/phbib.h>
#include <simthetic/phrand.h>
#include <simthetic/signals.h>
#include <simthetic/basicdevice.h>
#include <simthetic/interfaces.h>

#include <simtheticlib/delayprofile.h>

namespace simthlib {
    class OfdmReceiver;
    class OfdmTransmitter;
}

using namespace simthlib;

#include "matrixseq.h"
#include "mimo-channel.h"

// MIMO channels are only available if the matrix data types are
// available

namespace mimo {

class WSSUSChannel2 : public MIMOChannel
{
   public:
      class PathParam {
     public:
        double phase;

        int delay;

        double ampl;
      };

      class ScattererParam 
      {
     public:
        double angle;

        double dopp_freq;

        double delta_phi;
      };
   private:

    double rice_factor;

  protected:
    unsigned num_paths;

    double max_delay;

    double max_doppler;

    unsigned L_T;
    unsigned L_R;
    unsigned paths_per_scatterer;

    typedef std::vector<ScattererParam> ScatVect;
    ScatVect scatterers_tx;
    CMatrix steering_tx;
    CMatrix doppler_tx;

    ScatVect scatterers_rx;
    CMatrix steering_rx;
    CMatrix doppler_rx;

    typedef std::vector<PathParam> PathVect;
    PathVect h_pathparams;
    CMatrixSeq::storage_type h_paths;

    AntennaArray *array_geometry_tx;
    AntennaArray *array_geometry_rx;

    double delta_t;

   protected:
    VectorTimeSignal::storage_type prevsym;
    bool firstSequence;

    mutable MatrixTimeSignal::storage_type imp_resp;

   private:
    enum Normalization_Mode {
    STD,
    TIME,
    FREQ,
    TIME_TOTAL,
    FREQ_TOTAL
    };
    Normalization_Mode normMode;
    double normFactor;
    void resetNormalizeFactor();
    double freqDist;  
    int numFreqs;
    unsigned normDuration;
    unsigned normDurationLeft;

    // Possibly being used for the amplitude
    RayleighRng rayleigh_rng;
    
    void adjustDeltaT(double newDeltaT);

    simth::UniformRng AoARng;
    simth::UniformRng AoDRng;

    simth::UniformRng coupling_phase_rng;

  protected:
    // sequence to store the (ideal) estimated transfer factors in
    // order avoid (unessecary) multiple estimations in
    // getTransfer(..)
    mutable MatrixFreqSignal::storage_type cacheTransfer;
    mutable int cacheTime;
    mutable bool isTransferCacheValid;
    bool timevariant;

    int current_time;

    //------------------------------

  protected:
    int roundTimeToSamples(double value) const;

    virtual double diceNewDelay() const = 0;

    virtual double diceNewPhase() const = 0;

    virtual double diceNewAmpl() const;

    virtual double diceNewAoA() const;

    virtual double diceNewAoD() const;

  public:
    WSSUSChannel2(int deviceID, 
         DeviceSystemIntf* system,
         const simth::PropertyList& pl,
         double maxDelay,
          double max_doppler);

    virtual ~WSSUSChannel2();

    virtual void refresh();

    void getTransfer(int relative_time, MatrixFreqSignal& trans) const;

    void transmit(const VectorTimeSignal &insym, VectorTimeSignal& outsym);

    void setRiceFactor(double rf);

    void setRiceFactor_dB(double rf);


    void setTimeNormalization(unsigned duration = 0);

    void setFreqNormalization(double freqDist, int numFreqs, 
                  unsigned duration = 0);

    virtual void startOfSimulation();

    virtual void print(std::ostream &os) const;

   protected:
    void getImpulseResponse(int relative_time, MatrixTimeSignal &t, 
                int increment_time = 0) const;

};


inline int WSSUSChannel2::roundTimeToSamples(double value) const
{
    return int(value / delta_t + 0.5);
}






class WSSUSChannel2JakesUniform : public WSSUSChannel2
{
  private:
    UniformRng uniformPhaseRng;
  protected:
    virtual double diceNewPhase() const;
  public:
    WSSUSChannel2JakesUniform(int deviceID, 
                 DeviceSystemIntf* system,
                 const simth::PropertyList& pl,
                 double max_doppler,
                             double maxDelay);
    virtual ~WSSUSChannel2JakesUniform();
};


class WSSUSChannel2ExpDelay : public WSSUSChannel2JakesUniform
{
  private:
    ExponentialDelayProfile *expDelay;
  protected:
    virtual double diceNewDelay() const;
    WSSUSChannel2ExpDelay(int deviceID, 
             DeviceSystemIntf* system,
             const simth::PropertyList& pl,
             double max_doppler, double delayAttenuation,
                         double max_delay);
  public:
    WSSUSChannel2ExpDelay(int deviceID, 
             DeviceSystemIntf* system,
             const simth::PropertyList& pl);

    virtual ~WSSUSChannel2ExpDelay();
};



class WSSUSChannel2Scattering : public WSSUSChannel2ExpDelay
{
  public:
    WSSUSChannel2Scattering(int deviceID, 
               DeviceSystemIntf* system,
               const simth::PropertyList& pl)
    : WSSUSChannel2ExpDelay(deviceID, system, pl,
                   pl.getProperty<double>("max_doppler"), 
                   (-log(DELAY_THRESHOLD)/pl.getProperty<double>("max_signal_delay")), 
                   pl.getProperty<double>("max_signal_delay")) {};
};



class WSSUSChannel2TU : public WSSUSChannel2ExpDelay
{
  public:
    WSSUSChannel2TU(int deviceID, 
           DeviceSystemIntf* system,
           const simth::PropertyList& pl)
        :WSSUSChannel2ExpDelay(deviceID, system, pl,
                  pl.getProperty<double>("max_doppler"), 
                  1e6,
                  7e-6) {};
};



} // namespace mimo




#endif /* WSSUSCHANNEL_H */

---