double
FemtoCellUrbanAreaChannelRealization::GetPathLoss (void)
{
	/*
	 * Path loss Models from sect. 5.2 in
	 * 3GPP TSG RAN WG4 R4-092042
	 *
	 * Alternative simplified model based on LTE-A evaluation methodology which avoids modeling any walls.
	 */

  double distance;
  double minimumCouplingLoss = 45; //[dB] - see 3GPP TSG RAN WG4 #42bis (R4-070456)
  double floorPenetration = 0.0;
  //18.3 n ((n+2)/(n+1)-0.46)

  if (GetSourceNode ()->GetNodeType () == NetworkNode::TYPE_UE
		  && ( GetDestinationNode ()->GetNodeType () == NetworkNode::TYPE_ENODEB || GetDestinationNode ()->GetNodeType () == NetworkNode::TYPE_HOME_BASE_STATION) )
    {
	  UserEquipment* ue = (UserEquipment*) GetSourceNode ();
	  ENodeB* enb = (ENodeB*) GetDestinationNode ();

	  if( enb->GetCell()->GetCellCenterPosition()->GetCoordinateZ() > 0
			  && ue->IsIndoor()
			     && enb->GetCell()->GetCellCenterPosition()->GetCoordinateZ() != ue->GetCell()->GetCellCenterPosition()->GetCoordinateZ())
	  {
		  int n = (int) abs( enb->GetCell()->GetCellCenterPosition()->GetCoordinateZ() - ue->GetCell()->GetCellCenterPosition()->GetCoordinateZ() );
		  floorPenetration = 18.3 * pow( n, ((n+2)/(n+1)-0.46));
	  }

	  distance =  ue->GetMobilityModel ()->GetAbsolutePosition ()->GetDistance (enb->GetMobilityModel ()->GetAbsolutePosition ());
    }

  else if (GetDestinationNode ()->GetNodeType () == NetworkNode::TYPE_UE
		  && ( GetSourceNode ()->GetNodeType () == NetworkNode::TYPE_ENODEB || GetSourceNode ()->GetNodeType () == NetworkNode::TYPE_HOME_BASE_STATION) )
  {
	  UserEquipment* ue = (UserEquipment*) GetDestinationNode ();
	  ENodeB* enb = (ENodeB*) GetSourceNode ();

	  if( enb->GetCell()->GetCellCenterPosition()->GetCoordinateZ() > 0
			  && ue->IsIndoor()
			  && enb->GetCell()->GetCellCenterPosition()->GetCoordinateZ() != ue->GetCell()->GetCellCenterPosition()->GetCoordinateZ())
	  {
		  int n = (int) abs( enb->GetCell()->GetCellCenterPosition()->GetCoordinateZ() - ue->GetCell()->GetCellCenterPosition()->GetCoordinateZ() );
		  floorPenetration = 18.3 * pow( n, ((n+2)/(n+1)-0.46));
	  }

	  distance =  ue->GetMobilityModel ()->GetAbsolutePosition ()->GetDistance (enb->GetMobilityModel ()->GetAbsolutePosition ());
  }




  m_pathLoss = max( minimumCouplingLoss, 127 + ( 30 * log10 (distance * 0.001) ) + floorPenetration);


  return m_pathLoss;
}
Example #2
0
void
Manhattan::UpdatePosition (double timestamp)
{
#ifdef MOBILITY_DEBUG
//	cout << "\t START MOBILITY MODEL for "<< GetNodeID () << endl;
#endif

  if (GetSpeed () == 0)
    {
#ifdef MOBILITY_DEBUG
	cout << "\t\t speed = 0 --> position has not been updated!"<< endl;
#endif
      return;
    }

  double timeInterval = timestamp - GetPositionLastUpdate ();
  double speedDirection;

  UserEquipment *thisNode = NetworkManager::Init ()->GetUserEquipmentByID (GetNodeID ());
  Cell *thisCell = thisNode->GetCell ();


  NetworkNode *targetNode = thisNode->GetTargetNode ();

#ifdef MOBILITY_DEBUG
	cout << "MOBILITY_DEBUG: User ID: " << GetNodeID ()
	    << "\n\t Cell ID " <<
					NetworkManager::Init()->GetCellIDFromPosition (GetAbsolutePosition()->GetCoordinateX(),
																   GetAbsolutePosition()->GetCoordinateY())
		<< "\n\t Initial Position (X): " << GetAbsolutePosition()->GetCoordinateX()
		<< "\n\t Initial Position (Y): " << GetAbsolutePosition()->GetCoordinateY()
		<< "\n\t Speed: " << GetSpeed()
		<< "\n\t Speed Direction: " << GetSpeedDirection()
		<< "\n\t Time Last Update: " << GetPositionLastUpdate()
		<< "\n\t Time Interval: " << timeInterval
		<< endl;
#endif
#ifdef MOBILITY_DEBUG_TAB
	cout <<  GetNodeID () << " "
		<< GetAbsolutePosition()->GetCoordinateX() << " "
		<< GetAbsolutePosition()->GetCoordinateY() << " "
		<< GetSpeed() << " "
		<< GetSpeedDirection() << " "
		<< GetPositionLastUpdate() << " "
		<< timeInterval << " "
		<< "-> ";
#endif
	const double pi = 3.14;

//  Init speedDirection if its not pointing into the direction of Manhattans Streets
	vector<double> directions;
	for(int i=0;i<=3;i++) {
	  directions.push_back(i * 1.57); // 1.57 = 90.0 * ((2.0*pi)/360.0))
	}

	double newdir;
	vector<double>::iterator it;
	it = find(directions.begin(), directions.end(), GetSpeedDirection());

	if(it == directions.end()) {
	  newdir = (double)(pi / 2) * round(GetSpeedDirection() * (double)(2 / pi));
	  if(newdir==6.28) newdir = 0;
	  SetSpeedDirection(newdir);
#ifdef MOBILITY_DEBUG_TAB
	  cout << "NEW SPEED-DIRECTION: " << newdir << " -> ";
#endif
	}

	double shift = timeInterval * (GetSpeed()*(1000.0/3600.0));

	CartesianCoordinates *newPosition = new CartesianCoordinates(GetAbsolutePosition()->GetCoordinateX(), GetAbsolutePosition()->GetCoordinateY());
	CartesianCoordinates *ENodeBPosition = targetNode->GetMobilityModel ()->GetAbsolutePosition ();

// Init Manhattan grid position
  if(fmod(GetAbsolutePosition()->GetCoordinateY(),100)!=0 && fmod(GetAbsolutePosition()->GetCoordinateX(),100)!=0){
	CartesianCoordinates *Correction = new CartesianCoordinates();
	double distfromEnB = newPosition->GetDistance (ENodeBPosition);
	double azim = newPosition->GetPolarAzimut (ENodeBPosition);

	//if it was randomly put outside the cell -> shift it inside
	if(distfromEnB > (thisCell->GetRadius()*1000)) {
	  Correction->SetCoordinates((newPosition->GetDistance (ENodeBPosition) - (thisCell->GetRadius()*1000)) * cos(azim),
			  (newPosition->GetDistance (ENodeBPosition) - (thisCell->GetRadius()*1000)) * sin(azim));
	  newPosition->SetCoordinates(newPosition->GetCoordinateX() - Correction->GetCoordinateX(),
			  newPosition->GetCoordinateY() - Correction->GetCoordinateY());
	}
	delete Correction;


	if(GetSpeedDirection()==0 || GetSpeedDirection()==3.14) {
			  if(newPosition->GetCoordinateY() < 0)
				  newPosition->SetCoordinateY( ceil( (double)(newPosition->GetCoordinateY() / 100) ) * 100);
			  else
				  newPosition->SetCoordinateY( floor( (double)(newPosition->GetCoordinateY() / 100) ) * 100);
		  }
	else {
			  if(newPosition->GetCoordinateX() < 0)
				  newPosition->SetCoordinateX( ceil( (double)(newPosition->GetCoordinateX() / 100)) * 100);
			  else
				  newPosition->SetCoordinateX( floor( (double)(newPosition->GetCoordinateX() / 100) ) * 100);
	}
  }

//Shift it
  double shift_x = shift * cos(GetSpeedDirection());
  double shift_y = shift * sin(GetSpeedDirection());
  if(GetSpeedDirection()==0 || GetSpeedDirection()==3.14) {
	  newPosition->SetCoordinateX(newPosition->GetCoordinateX()+shift_x);
  }
  else {
	  newPosition->SetCoordinateY(newPosition->GetCoordinateY()+shift_y);
  }

// if a node reaches a crossing, choose new speedDirection
  double old_x = abs( ((int)( GetAbsolutePosition()->GetCoordinateX() *1000))/1000.0 ); //cut after 3 decimal places
  double new_x = abs( ((int)( newPosition->GetCoordinateX() *1000))/1000.0 );
  double old_y = abs( ((int)( GetAbsolutePosition()->GetCoordinateY() *1000))/1000.0 );
  double new_y = abs( ((int)( newPosition->GetCoordinateY() *1000))/1000.0 );
  double rounded_x = abs( round(old_x/100)*100 );
  double rounded_y = abs( round(old_y/100)*100 );

  if( ((old_x<rounded_x && rounded_x<=new_x) || (old_x>rounded_x && rounded_x>=new_x)) ||
		  ((old_y<rounded_y && rounded_y<=new_y) || (old_y>rounded_y && rounded_y>=new_y)) ||
		  (rounded_x==0 && old_x<rounded_x && rounded_x<=new_x) || (rounded_x==0 && old_x>rounded_x && rounded_x>=new_x) ||
		  (rounded_y==0 && old_y<rounded_y && rounded_y<=new_y) || (rounded_y==0 && old_y>rounded_y && rounded_y>=new_y) )
  {
	  srand ( time(NULL) );
	  double prob_turn = (rand()%100)*0.01;
	  if(prob_turn<=0.25) {
		  speedDirection = GetSpeedDirection() + 1.57; //turn left;
		  newPosition->SetCoordinates(round(newPosition->GetCoordinateX()),round(newPosition->GetCoordinateY()));
#ifdef MOBILITY_DEBUG_TAB
		  cout << "TURN LEFT: " << speedDirection << " -> ";
#endif
	  }
	  if(prob_turn>0.25 && prob_turn<0.75) {
		  newPosition->SetCoordinates(round(newPosition->GetCoordinateX()),round(newPosition->GetCoordinateY()) );
		  speedDirection = GetSpeedDirection();
#ifdef MOBILITY_DEBUG_TAB
		  cout << "no TURN: straight -> ";
#endif
	  }
	  if(prob_turn>=0.75) {
		  newPosition->SetCoordinates(round(newPosition->GetCoordinateX()),round(newPosition->GetCoordinateY()) );
		  speedDirection = GetSpeedDirection() - 1.57;
#ifdef MOBILITY_DEBUG_TAB
		  cout << "TURN RIGHT: " << speedDirection << " -> ";
#endif
	  }

	  //Correction if speedDirection €! [0,2pi[
	  if(speedDirection>=2*pi) speedDirection -= 2*pi;
	  if(speedDirection<0) speedDirection += 2*pi;
	  SetSpeedDirection(speedDirection);
  }



//If node moves beyond the cell edge
  double azimut = newPosition->GetPolarAzimut (ENodeBPosition);
  double newDistanceFromTheENodeB = newPosition->GetDistance (ENodeBPosition);

  if (newDistanceFromTheENodeB >= (thisCell->GetRadius()*1000))
    {
	  if (GetHandover()== false)
		{
		  newPosition->SetCoordinateX(GetAbsolutePosition()->GetCoordinateX());
		  newPosition->SetCoordinateY(GetAbsolutePosition()->GetCoordinateY());

		  speedDirection = GetSpeedDirection() - pi;
#ifdef MOBILITY_DEBUG_TAB
		  cout << "Moved back to cell: SPEED DIRECTION CHANGE (HO): " << speedDirection << " -> ";
#endif
		  if(speedDirection<0) speedDirection = speedDirection + 2*pi;
		  SetSpeedDirection(speedDirection);
		}
	  else if (newPosition->GetDistance(0.0, 0.0) >= GetTopologyBorder ())
	      {
		  newPosition->SetCoordinateX(GetAbsolutePosition()->GetCoordinateX());
		  newPosition->SetCoordinateY(GetAbsolutePosition()->GetCoordinateY());

		  speedDirection = GetSpeedDirection() - 1.57;
#ifdef MOBILITY_DEBUG_TAB
		  cout << "Moved back to cell: SPEED DIRECTION CHANGE (TOPOLOGY-BORDER): " << speedDirection << " -> ";
#endif
		  if(speedDirection<0) speedDirection = speedDirection + 2*pi;
		  SetSpeedDirection(speedDirection);
	      }
    }

  SetAbsolutePosition(newPosition);
  SetPositionLastUpdate (timestamp);

#ifdef MOBILITY_DEBUG
  cout << "\n\t Final Position (X): " << GetAbsolutePosition()->GetCoordinateX()
			<< "\n\t Final Position (Y): " << GetAbsolutePosition()->GetCoordinateY()
			<< endl;
#endif
#ifdef MOBILITY_DEBUG_TAB
  cout << GetAbsolutePosition()->GetCoordinateX() << " "
			<< GetAbsolutePosition()->GetCoordinateY() << " " << GetSpeedDirection()
			<< endl;
#endif
  delete newPosition;
}