papropagation.cc
上传用户:rrhhcc
上传日期:2015-12-11
资源大小:54129k
文件大小:4k
源码类别:

通讯编程

开发平台:

Visual C++

  1. /* -*- Mode:C++; c-basic-offset:8; tab-width:8; indent-tabs-mode:t -*- */
  2. /*
  3.  * Copyright (c) 1997 Regents of the University of California.
  4.  * All rights reserved.
  5.  *
  6.  * Redistribution and use in source and binary forms, with or without
  7.  * modification, are permitted provided that the following conditions
  8.  * are met:
  9.  * 1. Redistributions of source code must retain the above copyright
  10.  *    notice, this list of conditions and the following disclaimer.
  11.  * 2. Redistributions in binary form must reproduce the above copyright
  12.  *    notice, this list of conditions and the following disclaimer in the
  13.  *    documentation and/or other materials provided with the distribution.
  14.  * 3. All advertising materials mentioning features or use of this software
  15.  *    must display the following acknowledgement:
  16.  * This product includes software developed by the Computer Systems
  17.  * Engineering Group at Lawrence Berkeley Laboratory.
  18.  * 4. Neither the name of the University nor of the Laboratory may be used
  19.  *    to endorse or promote products derived from this software without
  20.  *    specific prior written permission.
  21.  *
  22.  * THIS SOFTWARE IS PROVIDED BY THE REGENTS AND CONTRIBUTORS ``AS IS'' AND
  23.  * ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE
  24.  * IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE
  25.  * ARE DISCLAIMED.  IN NO EVENT SHALL THE REGENTS OR CONTRIBUTORS BE LIABLE
  26.  * FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL
  27.  * DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS
  28.  * OR SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS INTERRUPTION)
  29.  * HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT
  30.  * LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY
  31.  * OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF
  32.  * SUCH DAMAGE.
  33.  * 
  34.  *
  35.  * Ported 2006 from ns-2.29 
  36.  * by Federico Maguolo, Nicola Baldo and Simone Merlin 
  37.  * (SIGNET lab, University of Padova, Department of Information Engineering)
  38.  * 
  39.  */
  40. #include "papropagation.h"
  41. // methods for free space model
  42. static class PAFreeSpaceClass: public TclClass {
  43. public:
  44. PAFreeSpaceClass() : TclClass("Propagation/FreeSpace/PowerAware") {}
  45. TclObject* create(int, const char*const*) {
  46. return (new PAFreeSpace);
  47. }
  48. } class_pafreespace;
  49. double PAFreeSpace::Pr(Packet* , PacketStamp *t, PacketStamp *r,   double L, double lambda)
  50. {
  51.   //double L = ifp->getL(); // system loss
  52.   //double lambda = ifp->getLambda();   // wavelength
  53. double Xt, Yt, Zt; // location of transmitter
  54. double Xr, Yr, Zr; // location of receiver
  55. t->getNode()->getLoc(&Xt, &Yt, &Zt);
  56. r->getNode()->getLoc(&Xr, &Yr, &Zr);
  57. // Is antenna position relative to node position?
  58. Xr += r->getAntenna()->getX();
  59. Yr += r->getAntenna()->getY();
  60. Zr += r->getAntenna()->getZ();
  61. Xt += t->getAntenna()->getX();
  62. Yt += t->getAntenna()->getY();
  63. Zt += t->getAntenna()->getZ();
  64. double dX = Xr - Xt;
  65. double dY = Yr - Yt;
  66. double dZ = Zr - Zt;
  67. double d = sqrt(dX * dX + dY * dY + dZ * dZ);
  68. //printf("DISTANZA %f %f %f %f %f%n", d, Xr, Yr, Xt, Yt);  
  69. // get antenna gain
  70. double Gt = t->getAntenna()->getTxGain(dX, dY, dZ, lambda);
  71. double Gr = r->getAntenna()->getRxGain(dX, dY, dZ, lambda);
  72. // calculate receiving power at distance
  73. double Pr = Friis(t->getTxPr(), Gt, Gr, lambda, L, d);
  74. //printf("%lf: d: %lf, Pr: %en", Scheduler::instance().clock(), d, Pr);
  75. return Pr;
  76. }
  77. double
  78. PAFreeSpace::getDist(double Pr, double Pt, double Gt, double Gr, double hr, double ht, double L, double lambda)
  79. {
  80.         return sqrt((Pt * Gt * Gr * lambda * lambda) / (L * Pr)) /
  81.                 (4 * PI);
  82. }