dadada/ns-3-leo
dadada/ns-3-leo
1
0
Fork 0
mirror of https://gitlab.ibr.cs.tu-bs.de/tschuber/ns-3-leo.git synced 2025-06-08 10:03:58 +02:00
Code Issues Releases Wiki Activity

Fix bugs and write data directly to file

Browse source
This commit is contained in:
Tim Schubert 2020-08-17 13:30:39 +02:00
parent f8804533cc
commit bf59e66526
3 changed files with 34 additions and 21 deletions
Download patch file Download diff file Expand all files Collapse all files

13
examples/leo-circular-orbit-tracing-example.cc
View file

@ -1,5 +1,7 @@
/* -*- Mode:C++; c-file-style:"gnu"; indent-tabs-mode:nil; -*- */
#include <fstream>
#include "ns3/core-module.h"
#include "ns3/mobility-module.h"
#include "ns3/leo-module.h"
@ -8,15 +10,19 @@ using namespace ns3;
NS_LOG_COMPONENT_DEFINE ("LeoCircularOrbitTracingExample");
ofstream outfile ("leo-circular-orbit-tracing-example.csv");
void CourseChange (std::string context, Ptr<const MobilityModel> position)
{
Vector pos = position->GetPosition ();
Ptr<const Node> node = position->GetObject<Node> ();
std::cout << Simulator::Now () << "," << node->GetId () << "," << pos.x << "," << pos.y << "," << pos.z << "," << position->GetVelocity ().GetLength() << std::endl;
outfile << Simulator::Now () << "," << node->GetId () << "," << pos.x << "," << pos.y << "," << pos.z << "," << position->GetVelocity ().GetLength() << std::endl;
}
int main(int argc, char *argv[])
{
for (double incl: { 70.0, 100.0 })
{
NodeContainer c;
c.Create (600);
@ -26,14 +32,15 @@ int main(int argc, char *argv[])
"NumSatellites", IntegerValue (100));
mobility.SetMobilityModel ("ns3::LeoCircularOrbitMobilityModel",
"Altitude", DoubleValue (1200.0),
"Inclination", DoubleValue (80.0),
"Inclination", DoubleValue (incl),
"Precision", TimeValue (Minutes (1)));
mobility.Install (c);
}
Config::Connect ("/NodeList/*/$ns3::MobilityModel/CourseChange",
MakeCallback (&CourseChange));
std::cout << "Time,Satellite,x,y,z,Speed" << std::endl;
outfile << "Time,Satellite,x,y,z,Speed" << std::endl;
Simulator::Stop (Minutes (60.0));
Simulator::Run ();

14
model/leo-circular-orbit-mobility-model.cc
View file

@ -76,6 +76,7 @@ DotProduct (const Vector3D &l, const Vector3D &r)
double
LeoCircularOrbitMobilityModel::GetSpeed () const
{
// force non-retrograd movement
return sqrt (LEO_EARTH_GM_KM_E10 / m_orbitHeight) * 1e5;
}
@ -90,9 +91,14 @@ LeoCircularOrbitMobilityModel::DoGetVelocity () const
Vector3D
LeoCircularOrbitMobilityModel::PlaneNorm () const
{
return Vector3D (sin (-m_inclination) * cos (m_latitude),
sin (-m_inclination) * sin (m_latitude),
cos (m_inclination));
int sign = 1;
if (m_inclination < M_PI/4)
{
sign = -1;
}
return Vector3D (sign * sin (-m_inclination) * cos (m_latitude),
sign * sin (-m_inclination) * sin (m_latitude),
sign * cos (m_inclination));
}
double
@ -173,7 +179,7 @@ void LeoCircularOrbitMobilityModel::SetAltitude (double h)
double LeoCircularOrbitMobilityModel::GetInclination () const
{
return (m_inclination / M_PI) * 180.0;
return (m_inclination / (M_PI)) * 180.0;
}
void LeoCircularOrbitMobilityModel::SetInclination (double incl)

6
utils/plot-satellites.gnuplot
View file

@ -7,13 +7,13 @@ unset xtics
unset ytics
unset ztics
unset border
set hidden3d
unset key
set view equal xyz
n=600
# number of nodes per time slot
n=1200
do for [j=0:100] {
set title 'time '.j
splot 'somefile' using 3:4:5:2 every ::(j*n)::((j+1)*n)
splot 'leo-circular-orbit-tracing-example.csv' using 3:4:5:2 every ::(j*n)::((j+1)*n)
}