sat-markov-fading-trace-example.cc

Go to the documentation of this file.

/* -*- Mode:C++; c-file-style:"gnu"; indent-tabs-mode:nil; -*- */
/*
 * Copyright (c) 2014 Magister Solutions
 *
 * This program is free software; you can redistribute it and/or modify
 * it under the terms of the GNU General Public License version 2 as
 * published by the Free Software Foundation;
 *
 * 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
 *
 * Author: Frans Laakso <frans.laakso@magister.fi>
 *
 */
 
#include "ns3/core-module.h"
#include "ns3/gnuplot.h"
#include "ns3/satellite-module.h"
 
#include <fstream>
#include <vector>
 
NS_LOG_COMPONENT_DEFINE("sat-markov-trace-example");
 
namespace ns3
{

 
class SatMarkovFadingExamplePlot : public Object
{
  public:
    SatMarkovFadingExamplePlot();
    static TypeId GetTypeId();
    void Run();
 
  private:
    Gnuplot2dDataset GetGnuplotDataset(std::string title);
    Gnuplot GetGnuplot(std::string outputName, std::string title);
    void FadingTraceCb(std::string context,
                       double time,
                       SatEnums::ChannelType_t chType,
                       double fadingValue);
    double GetElevation();
    double GetVelocity();
 
    double m_elevation;
    double m_velocity;
    std::vector<std::pair<double, double>> m_fadingValues;
};
 
NS_OBJECT_ENSURE_REGISTERED(SatMarkovFadingExamplePlot);
 
TypeId
SatMarkovFadingExamplePlot::GetTypeId()
{
    static TypeId tid = TypeId("ns3::SatMarkovFadingExamplePlot")
                            .SetParent<Object>()
                            .AddConstructor<SatMarkovFadingExamplePlot>();
    return tid;
}
 
SatMarkovFadingExamplePlot::SatMarkovFadingExamplePlot()
{
    m_elevation = 45;
    m_velocity = 0;
}
 
void
SatMarkovFadingExamplePlot::FadingTraceCb(std::string context,
                                          double time,
                                          SatEnums::ChannelType_t chType,
                                          double fadingValue)
{
    std::cout << time << " " << chType << " " << 10 * log10(fadingValue) << std::endl;
    m_fadingValues.push_back(std::make_pair(time, 10 * log10(fadingValue)));
}
 
double
SatMarkovFadingExamplePlot::GetElevation()
{
    return m_elevation;
}
 
double
SatMarkovFadingExamplePlot::GetVelocity()
{
    return m_velocity;
}
 
Gnuplot2dDataset
SatMarkovFadingExamplePlot::GetGnuplotDataset(std::string title)
{
    Gnuplot2dDataset ret;
    ret.SetTitle(title);
    ret.SetStyle(Gnuplot2dDataset::LINES);
 
    for (uint32_t i = 0; i < m_fadingValues.size(); i++)
    {
        ret.Add(m_fadingValues[i].first, m_fadingValues[i].second);
    }
    return ret;
}
 
Gnuplot
SatMarkovFadingExamplePlot::GetGnuplot(std::string outputName, std::string title)
{
    Gnuplot ret(outputName + ".png");
    ret.SetTitle(title);
    ret.SetTerminal("png");
    ret.SetLegend("Time (s)", "Fading (dB)");
    ret.AppendExtra("set key top right");
    ret.AppendExtra("set grid xtics mxtics ytics");
    return ret;
}
 
void
SatMarkovFadingExamplePlot::Run()
{
    Config::SetDefault("ns3::SatEnvVariables::SimulationCampaignName",
                       StringValue("example-markov-fading-trace"));
    Config::SetDefault("ns3::SatEnvVariables::SimulationTag", StringValue(""));
    Config::SetDefault("ns3::SatEnvVariables::EnableSimulationOutputOverwrite", BooleanValue(true));

    Ptr<SatMarkovConf> markovConf = CreateObject<SatMarkovConf>();
 
    SatBaseFading::ElevationCallback elevationCb =
        MakeCallback(&SatMarkovFadingExamplePlot::GetElevation, this);
    SatBaseFading::VelocityCallback velocityCb =
        MakeCallback(&SatMarkovFadingExamplePlot::GetVelocity, this);

    Ptr<SatMarkovContainer> markovContainer =
        CreateObject<SatMarkovContainer>(markovConf, elevationCb, velocityCb);
 
    markovContainer->TraceConnect("FadingTrace",
                                  "The trace for fading values",
                                  MakeCallback(&SatMarkovFadingExamplePlot::FadingTraceCb, this));
 
    Address macAddress;

    for (uint32_t i = 0; i < 100000; i++)
    {
        Simulator::Schedule(MilliSeconds(1 * i),
                            &SatMarkovContainer::DoGetFading,
                            markovContainer,
                            macAddress,
                            SatEnums::FORWARD_USER_CH);
    }
 
    Simulator::Schedule(MilliSeconds(0),
                        &SatMarkovContainer::LockToSetAndState,
                        markovContainer,
                        0,
                        0);
    Simulator::Schedule(MilliSeconds(20000),
                        &SatMarkovContainer::LockToSetAndState,
                        markovContainer,
                        0,
                        1);
    Simulator::Schedule(MilliSeconds(40000),
                        &SatMarkovContainer::LockToSetAndState,
                        markovContainer,
                        0,
                        2);
    Simulator::Schedule(MilliSeconds(60000),
                        &SatMarkovContainer::LockToSetAndState,
                        markovContainer,
                        0,
                        0);
    Simulator::Schedule(MilliSeconds(80000),
                        &SatMarkovContainer::LockToSetAndState,
                        markovContainer,
                        0,
                        1);
 
    Simulator::Run();
 
    Gnuplot2dDataset dataset = GetGnuplotDataset("Markov Fading Trace");
    Gnuplot plot = GetGnuplot("markov_fading_trace", "Markov Fading Trace");
    plot.AddDataset(dataset);
 
    std::string plotFileName = "markov_fading_trace.plt";
    std::ofstream plotFile(plotFileName.c_str());
    plot.GenerateOutput(plotFile);
    plotFile.close();
 
    std::cout << "Output file written: " << plotFileName << std::endl;
 
    int result = system("gnuplot markov_fading_trace.plt");
 
    if (result < 0)
    {
        std::cout
            << "Unable to open shell process for Gnuplot file conversion, conversion not done!"
            << std::endl;
    }
    else
    {
        std::cout << "Output file converted to: markov_fading_trace.png" << std::endl;
    }
 
    Simulator::Destroy();
}
 
} // namespace ns3
 
int
main(int argc, char* argv[])
{
    ns3::Ptr<ns3::SatMarkovFadingExamplePlot> stub;
    stub = ns3::CreateObject<ns3::SatMarkovFadingExamplePlot>();
    ns3::Config::RegisterRootNamespaceObject(stub);
    stub->Run();
 
    return 0;
}