mirror of
https://gitlab.ibr.cs.tu-bs.de/tschuber/ns-3-leo.git
synced 2025-06-08 01:53:58 +02:00
103 lines
2.4 KiB
C++
103 lines
2.4 KiB
C++
/* -*- Mode:C++; c-file-style:"gnu"; indent-tabs-mode:nil; -*- */
|
|
|
|
#include "ns3/log.h"
|
|
#include "ns3/mobility-model.h"
|
|
#include "ns3/double.h"
|
|
#include "math.h"
|
|
|
|
#include "isl-propagation-loss-model.h"
|
|
|
|
namespace ns3 {
|
|
|
|
NS_LOG_COMPONENT_DEFINE ("IslPropagationLossModel");
|
|
|
|
NS_OBJECT_ENSURE_REGISTERED (IslPropagationLossModel);
|
|
|
|
TypeId
|
|
IslPropagationLossModel::GetTypeId (void)
|
|
{
|
|
static TypeId tid = TypeId ("ns3::IslPropagationLossModel")
|
|
.SetParent<PropagationLossModel> ()
|
|
.SetGroupName ("Leo")
|
|
.AddConstructor<IslPropagationLossModel> ()
|
|
;
|
|
return tid;
|
|
}
|
|
|
|
IslPropagationLossModel::IslPropagationLossModel ()
|
|
{
|
|
}
|
|
|
|
IslPropagationLossModel::~IslPropagationLossModel ()
|
|
{
|
|
}
|
|
|
|
bool
|
|
IslPropagationLossModel::GetLos (Ptr<MobilityModel> moda, Ptr<MobilityModel> modb)
|
|
{
|
|
// origin of LOS
|
|
Vector3D apos = moda->GetPosition ();
|
|
Vector3D bpos = modb->GetPosition ();
|
|
|
|
// select upper satellite as origin
|
|
Vector oc = apos.GetLength () > bpos.GetLength () ? apos : bpos;
|
|
Vector bp = apos.GetLength () > bpos.GetLength () ? bpos : apos;
|
|
|
|
Vector3D u = bp - oc;
|
|
// distance from s1 to s2
|
|
double s2 = u.GetLength ();
|
|
// direction unit vector
|
|
u = Vector3D (u.x / s2, u.y / s2, u.z / s2);
|
|
|
|
double a = u.x*u.x + u.y*u.y + u.z*u.z;
|
|
double b = 2.0 * (oc.x*u.x + oc.y*u.y + oc.z*u.z);
|
|
double c = (oc.x*oc.x + oc.y*oc.y + oc.z*oc.z) - (LEO_EARTH_RAD*LEO_EARTH_RAD);
|
|
double discriminant = b*b - 4*a*c;
|
|
|
|
NS_LOG_DEBUG ("a_pos="<<moda->GetPosition ()<<";b_pos"<<modb->GetPosition ()
|
|
<<";u="<<u
|
|
<<";a="<<a
|
|
<<";b="<<b
|
|
<<";c="<<c
|
|
<<";disc="<<discriminant);
|
|
|
|
if (discriminant < 0)
|
|
{
|
|
return true;
|
|
}
|
|
else
|
|
{
|
|
// distances from oc (first sat)
|
|
double t1 = (-b - sqrt (discriminant)) / (2.0 * a);
|
|
double t2 = (-b + sqrt (discriminant)) / (2.0 * a);
|
|
|
|
NS_LOG_DEBUG ("t1="<<t1
|
|
<<"t2="<<t2
|
|
<<"s2="<<s2);
|
|
|
|
// check if second sat is behind earth
|
|
return (s2 < abs (t1) && s2 < abs (t2));
|
|
}
|
|
}
|
|
|
|
double
|
|
IslPropagationLossModel::DoCalcRxPower (double txPowerDbm,
|
|
Ptr<MobilityModel> a,
|
|
Ptr<MobilityModel> b) const
|
|
{
|
|
NS_LOG_FUNCTION (this << txPowerDbm << a << b);
|
|
if (!GetLos (a, b))
|
|
{
|
|
return -1000.0;
|
|
}
|
|
|
|
return txPowerDbm;
|
|
}
|
|
|
|
int64_t
|
|
IslPropagationLossModel::DoAssignStreams (int64_t stream)
|
|
{
|
|
return 0;
|
|
}
|
|
|
|
};
|