WifiSpectrumValueHelper::GetSpectrumModel (uint32_t centerFrequency, uint32_t channelWidth, double bandBandwidth)

WifiSpectrumValueHelper::GetSpectrumModel (uint32_t centerFrequency, uint32_t channelWidth, double bandBandwidth, uint32_t guardBandwidth)

  NS_LOG_FUNCTION (centerFrequency << channelWidth << bandBandwidth);

  NS_LOG_FUNCTION (centerFrequency << channelWidth << bandBandwidth << guardBandwidth);

      // Overall bandwidth will be channelWidth plus 10 MHz guards on each side

      double bandwidth = (channelWidth + (2 * guardBandwidth)) * 1e6;

      double bandwidth = (channelWidth + 20) * 1e6;

      uint32_t numBands = static_cast<uint32_t> (bandwidth / bandBandwidth + 0.5);

      uint32_t numBands = static_cast<uint32_t> ((bandwidth / bandBandwidth) + 0.5);

WifiSpectrumValueHelper::CreateDsssTxPowerSpectralDensity (uint32_t centerFrequency, double txPowerW)

WifiSpectrumValueHelper::CreateDsssTxPowerSpectralDensity (uint32_t centerFrequency, double txPowerW, uint32_t guardBandwidth)

  NS_LOG_FUNCTION (centerFrequency << txPowerW);

  NS_LOG_FUNCTION (centerFrequency << txPowerW << guardBandwidth);

  Ptr<SpectrumValue> c = Create<SpectrumValue> (GetSpectrumModel (centerFrequency, channelWidth, 312500));

  double bandBandwidth = 312500;

  Ptr<SpectrumValue> c = Create<SpectrumValue> (GetSpectrumModel (centerFrequency, channelWidth, bandBandwidth, guardBandwidth));

  double txPowerPerBand;

  uint32_t nGuardBands = static_cast<uint32_t>(((2 * guardBandwidth * 1e6) / bandBandwidth) + 0.5);

  // Evenly spread power across 22 MHz (71 bands)

  uint32_t nAllocatedBands = static_cast<uint32_t>(((channelWidth * 1e6) / bandBandwidth) + 0.5);

  NS_ASSERT (c->GetSpectrumModel ()->GetNumBands () == 135);

  NS_ASSERT (c->GetSpectrumModel ()->GetNumBands () == (nAllocatedBands + nGuardBands + 1));

  txPowerPerBand = txPowerW / 71;

  // Evenly spread power across 22 MHz

  double txPowerPerBand = txPowerW / nAllocatedBands;

      if (i >= 32 && i <= 102)

      if ((i >= (nGuardBands / 2)) && (i <= ((nGuardBands / 2) + nAllocatedBands - 1)))

WifiSpectrumValueHelper::CreateOfdmTxPowerSpectralDensity (uint32_t centerFrequency, uint32_t channelWidth, double txPowerW)

WifiSpectrumValueHelper::CreateOfdmTxPowerSpectralDensity (uint32_t centerFrequency, uint32_t channelWidth, double txPowerW, uint32_t guardBandwidth)

  NS_LOG_FUNCTION (centerFrequency << channelWidth << txPowerW);

  NS_LOG_FUNCTION (centerFrequency << channelWidth << txPowerW << guardBandwidth);

  Ptr<SpectrumValue> c = Create<SpectrumValue> (GetSpectrumModel (centerFrequency, channelWidth, 312500));

  double bandBandwidth = 312500;

  Values::iterator vit = c->ValuesBegin ();

  Ptr<SpectrumValue> c = Create<SpectrumValue> (GetSpectrumModel (centerFrequency, channelWidth, bandBandwidth, guardBandwidth));

  Bands::const_iterator bit = c->ConstBandsBegin ();

  uint32_t nGuardBands = static_cast<uint32_t>(((2 * guardBandwidth * 1e6) / bandBandwidth) + 0.5);

  double txPowerPerBand;

  uint32_t nAllocatedBands = static_cast<uint32_t>(((channelWidth * 1e6) / bandBandwidth) + 0.5);

  NS_ASSERT_MSG (c->GetSpectrumModel ()->GetNumBands () == (nAllocatedBands + nGuardBands + 1), "Unexpected number of bands");

  double txPowerPerBand = 0;

  uint32_t start1 = 0;

  uint32_t stop1 = 0;

  uint32_t start2 = 0;

  uint32_t stop2 = 0;

      // skip 38 subbands, then place power in 26 subbands, then

      // skip guard band and 6 subbands, then place power in 26 subbands, then

      // the final 38 subbands.  

      // the final 6 subbands and the guard band.

      NS_ASSERT_MSG (c->GetSpectrumModel ()->GetNumBands () == 129, "Unexpected number of bands");

      for (size_t i = 0; i < c->GetSpectrumModel ()->GetNumBands (); i++, vit++, bit++)

      start1 = (nGuardBands / 2) + 6;

        {

      stop1 = start1 + 26 - 1;

          if ((i >= 38 && i <= 63) || (i >= 65 && i <= 90))

      start2 = stop1 + 2;

            {

      stop2 = start2 + 26 - 1;

              *vit = txPowerPerBand / (bit->fh - bit->fl);

            }

          else

            {

              *vit = 0;

            }

        }

      NS_LOG_DEBUG ("Added signal power to subbands 38-63 and 65-90");

      NS_LOG_DEBUG ("Integrated power " << Integral (*c));

      NS_ASSERT_MSG (std::abs (txPowerW - Integral (*c)) < 1e-6, "Power allocation failed"); 

      // skip 34 subbands, then place power in 14 subbands, then

      // skip guard band and 2 subbands, then place power in 14 subbands, then

      // the final 34 subbands.  

      // the final 2 subbands and the guard band.

      NS_ASSERT_MSG (c->GetSpectrumModel ()->GetNumBands () == 97, "Unexpected number of bands");

      for (size_t i = 0; i < c->GetSpectrumModel ()->GetNumBands (); i++, vit++, bit++)

      start1 = (nGuardBands / 2) + 2;

        {

      stop1 = start1 + 14 - 1;

          if ((i >= 34 && i <= 47) || (i >= 49 && i <= 62))

      start2 = stop1 + 2;

            {

      stop2 = start2 + 14 - 1;

              *vit = txPowerPerBand / (bit->fh - bit->fl);

            }

          else

            {

              *vit = 0;

            }

        }

      NS_LOG_DEBUG ("Added signal power to subbands 34-47 and 49-62");

      NS_LOG_DEBUG ("Integrated power " << Integral (*c));

      NS_ASSERT_MSG (std::abs (txPowerW - Integral (*c)) < 1e-6, "Power allocation failed"); 

      // skip 34 subbands, then place power in 14 subbands, then

      // skip guard band and 2 subbands, then place power in 8 subbands, then

      // skip the center subband, then place power in 14 subbands, then skip

      // skip the center subband, then place power in 8 subbands, then skip

      // the final 34 subbands.  

      // the final 2 subbands and the guard.

      NS_ASSERT_MSG (c->GetSpectrumModel ()->GetNumBands () == 81, "Unexpected number of bands");

          if ((i >= 32 && i <= 39) || (i >= 41 && i <= 48))

          if ((i >= start1 && i <= stop1) || (i >= start2 && i <= stop2))

      NS_LOG_DEBUG ("Added signal power to subbands 32-39 and 41-48");

      NS_LOG_DEBUG ("Added signal power to subbands " << start1 << "-" << stop1 <<

      NS_LOG_DEBUG ("Integrated power " << Integral (*c));

                    " and " << start2 << "-" << stop2);

      NS_ASSERT_MSG (std::abs (txPowerW - Integral (*c)) < 1e-6, "Power allocation failed"); 

      break;

    case 40:

      // 112 subcarriers (104 data + 8 pilot) 

      // possible alternative:  114 subcarriers (108 data + 6 pilot)

      txPowerPerBand = txPowerW / 112;

      NS_LOG_DEBUG ("Power per band " << txPowerPerBand);

      for (size_t i = 0; i < c->GetSpectrumModel ()->GetNumBands (); i++, vit++, bit++)

        {

          if ((i >= start1 && i <= stop1) || (i >= start2 && i <= stop2) ||

              (i >= start3 && i <= stop3) || (i >= start4 && i <= stop4))

            {

              *vit = txPowerPerBand / (bit->fh - bit->fl);

            }

          else

            {

              *vit = 0;

            }

        }

      NS_LOG_DEBUG ("Added signal power to subbands " << start1 << "-" << stop1 <<

                    ", " << start2 << "-" << stop2 <<

                    ", " << start3 << "-" << stop3 <<

                    ", " << start4 << "-" << stop4);

      break;

    case 80:

      // 224 subcarriers (208 data + 16 pilot) 

      // possible alternative:  242 subcarriers (234 data + 8 pilot)

      txPowerPerBand = txPowerW / 224;

      NS_LOG_DEBUG ("Power per band " << txPowerPerBand);

      for (size_t i = 0; i < c->GetSpectrumModel ()->GetNumBands (); i++, vit++, bit++)

        {

          if ((i >= start1 && i <= stop1) || (i >= start2 && i <= stop2) ||

              (i >= start3 && i <= stop3) || (i >= start4 && i <= stop4) ||

              (i >= start5 && i <= stop5) || (i >= start6 && i <= stop6) ||

              (i >= start7 && i <= stop7) || (i >= start8 && i <= stop8))

            {

              *vit = txPowerPerBand / (bit->fh - bit->fl);

            }

          else

            {

              *vit = 0;

            }

        }

      NS_LOG_DEBUG ("Added signal power to subbands " << start1 << "-" << stop1 <<

                    ", " << start2 << "-" << stop2 <<

                    ", " << start3 << "-" << stop3 <<

                    ", " << start4 << "-" << stop4 <<

                    ", " << start5 << "-" << stop5 <<

                    ", " << start6 << "-" << stop6 <<

                    ", " << start7 << "-" << stop7 <<

                    ", " << start8 << "-" << stop8);

      break;

    case 160:

      // 448 subcarriers (416 data + 32 pilot)

      // possible alternative:  484 subcarriers (468 data + 16 pilot)

      txPowerPerBand = txPowerW / 448;

      NS_LOG_DEBUG ("Power per band " << txPowerPerBand);

      for (size_t i = 0; i < c->GetSpectrumModel ()->GetNumBands (); i++, vit++, bit++)

        {

          if ((i >= start1 && i <= stop1) || (i >= start2 && i <= stop2) ||

              (i >= start3 && i <= stop3) || (i >= start4 && i <= stop4) ||

              (i >= start5 && i <= stop5) || (i >= start6 && i <= stop6) ||

              (i >= start7 && i <= stop7) || (i >= start8 && i <= stop8) ||

              (i >= start9 && i <= stop9) || (i >= start10 && i <= stop10) ||

              (i >= start11 && i <= stop11) || (i >= start12 && i <= stop12) ||

              (i >= start13 && i <= stop13) || (i >= start14 && i <= stop14) ||

              (i >= start15 && i <= stop15) || (i >= start16 && i <= stop16))

            {

              *vit = txPowerPerBand / (bit->fh - bit->fl);

            }

          else

            {

              *vit = 0;

            }

        }

      NS_LOG_DEBUG ("Added signal power to subbands " << start1 << "-" << stop1 <<

                    ", " << start2 << "-" << stop2 <<

                    ", " << start3 << "-" << stop3 <<

                    ", " << start4 << "-" << stop4 <<

                    ", " << start5 << "-" << stop5 <<

                    ", " << start6 << "-" << stop6 <<

                    ", " << start7 << "-" << stop7 <<

                    ", " << start8 << "-" << stop8 <<

                    ", " << start9 << "-" << stop9 <<

                    ", " << start10 << "-" << stop10 <<

                    ", " << start11 << "-" << stop11 <<

                    ", " << start12 << "-" << stop12 <<

                    ", " << start13 << "-" << stop13 <<

                    ", " << start14 << "-" << stop14 <<

                    ", " << start15 << "-" << stop15 <<

                    ", " << start16 << "-" << stop16);

WifiSpectrumValueHelper::CreateHtOfdmTxPowerSpectralDensity (uint32_t centerFrequency, uint32_t channelWidth, double txPowerW)

WifiSpectrumValueHelper::CreateHeOfdmTxPowerSpectralDensity (uint32_t centerFrequency, uint32_t channelWidth, double txPowerW, uint32_t guardBandwidth)

  NS_LOG_FUNCTION (centerFrequency << channelWidth << txPowerW);

  NS_LOG_FUNCTION (centerFrequency << channelWidth << txPowerW << guardBandwidth);

  Ptr<SpectrumValue> c = Create<SpectrumValue> (GetSpectrumModel (centerFrequency, channelWidth, 312500));

  double bandBandwidth = 78125;

  Ptr<SpectrumValue> c = Create<SpectrumValue> (GetSpectrumModel (centerFrequency, channelWidth, bandBandwidth, guardBandwidth));

      // 56 subcarriers (52 data + 4 pilot) 

      // 242 subcarriers (234 data + 8 pilot)

      NS_ASSERT_MSG (c->GetSpectrumModel ()->GetNumBands () == 129, "Unexpected number of bands");

      txPowerPerBand = txPowerW / 242;

      // skip 32 subbands, then place power in 28 of the next 32 subbands, then

      NS_LOG_DEBUG ("Power per band " << txPowerPerBand);

      // skip the center subband, then place power in 28 of the next 32

      // skip the guard band and 11 subbands, then place power in 121 subbands, then

      // subbands, then skip the final 32 subbands.  

      // skip 3 DC, then place power in 121 subbands, then skip

      txPowerPerBand = txPowerW / 56;

      // the final 11 subbands and the guard band.

      start1 = (nGuardBands / 2) + 12;

      stop1 = start1 + 121 - 1;

      start2 = stop1 + 4;

      stop2 = start2 + 121 - 1;

          if ((i >= 36 && i <= 63) || (i >= 65 && i <= 92))

          if ((i >= start1 && i <= stop1) || (i >= start2 && i <= stop2))

            {

              *vit = txPowerPerBand / (bit->fh - bit->fl);

            }

        }

      NS_LOG_DEBUG ("Added signal power to subbands 36-63 and 65-92");

      NS_LOG_DEBUG ("Integrated power " << Integral (*c));

      NS_ASSERT_MSG (std::abs (txPowerW - Integral (*c)) < 1e-6, "Power allocation failed"); 

      break;

    case 40:

      // 112 subcarriers (104 data + 8 pilot) 

      // possible alternative:  114 subcarriers (108 data + 6 pilot)

      NS_ASSERT_MSG (c->GetSpectrumModel ()->GetNumBands () == 193, "Unexpected number of bands");

      txPowerPerBand = txPowerW / 112;

      for (size_t i = 0; i < c->GetSpectrumModel ()->GetNumBands (); i++, vit++, bit++)

        {

          if ((i >= 36 && i <= 63) || (i >= 65 && i <= 92) ||

              (i >= 100 && i<= 127) || (i >= 129 && i <= 156))

      NS_LOG_DEBUG ("Added signal power to subbands 36-63, 65-92, 100-127, and 129-156");

      NS_LOG_DEBUG ("Added signal power to subbands " << start1 << "-" << stop1 <<

      NS_LOG_DEBUG ("Integrated power " << Integral (*c));

                    " and " << start2 << "-" << stop2);

      NS_ASSERT_MSG (std::abs (txPowerW - Integral (*c)) < 1e-6, "Power allocation failed"); 

    case 80:

    case 40:

      // 224 subcarriers (208 data + 16 pilot) 

      // 484 subcarriers (468 data + 16 pilot)

      // possible alternative:  242 subcarriers (234 data + 8 pilot)

      txPowerPerBand = txPowerW / 484;

      NS_ASSERT_MSG (c->GetSpectrumModel ()->GetNumBands () == 321, "Unexpected number of bands");

      NS_LOG_DEBUG ("Power per band " << txPowerPerBand);

      txPowerPerBand = txPowerW / 224;

      // skip the guard band and 11 subbands, then place power in 242 subbands, then

      // skip 5 DC, then place power in 242 subbands, then skip

      // the final 11 subbands and the guard band.

      start1 = (nGuardBands / 2) + 12;

      stop1 = start1 + 242 - 1;

      start2 = stop1 + 6;

      stop2 = start2 + 242 - 1;

          if ((i >= 36 && i <= 63) || (i >= 65 && i <= 92) || 

          if ((i >= start1 && i <= stop1) || (i >= start2 && i <= stop2))

              (i >= 100 && i <= 127) || (i >= 129 && i <= 156) ||

              (i >= 164 && i <= 191) || (i >= 193 && i <= 220) ||

              (i >= 228 && i <= 255) || (i >= 257 && i <= 284))

      NS_LOG_DEBUG ("Added signal power to subbands 36-63, 65-92, 100-127, 129-156, 164-191, 193-220, 228-255, 257-284");

      NS_LOG_DEBUG ("Added signal power to subbands " << start1 << "-" << stop1 <<

      NS_LOG_DEBUG ("Integrated power " << Integral (*c));

                    " and " << start2 << "-" << stop2);

      NS_ASSERT_MSG (std::abs (txPowerW - Integral (*c)) < 1e-6, "Power allocation failed"); 

    case 160:

    case 80:

      // 448 subcarriers (416 data + 32 pilot)

      // 996 subcarriers (980 data + 16 pilot)

      // possible alternative:  484 subcarriers (468 data + 16 pilot)

      txPowerPerBand = txPowerW / 996;

      NS_ASSERT_MSG (c->GetSpectrumModel ()->GetNumBands () == 577, "Unexpected number of bands");

      NS_LOG_DEBUG ("Power per band " << txPowerPerBand);

      txPowerPerBand = txPowerW / 448;

      // skip the guard band and 11 subbands, then place power in 498 subbands, then

      // skip 5 DC, then place power in 498 subbands, then skip

      // the final 11 subbands and the guard band.

      start1 = (nGuardBands / 2) + 12;

      stop1 = start1 + 498 - 1;

      start2 = stop1 + 6;

      stop2 = start2 + 498 - 1;

          if ((i >= 36 && i <= 63) || (i >= 65 && i <= 92) || 

          if ((i >= start1 && i <= stop1) || (i >= start2 && i <= stop2))

              (i >= 100 && i <= 127) || (i >= 129 && i <= 156) ||

              (i >= 164 && i <= 191) || (i >= 193 && i <= 220) ||

              (i >= 228 && i <= 255) || (i >= 257 && i <= 284) ||

              (i >= 292 && i <= 319) || (i >= 321 && i <= 348) ||

              (i >= 356 && i <= 383) || (i >= 385 && i <= 412) ||

              (i >= 420 && i <= 447) || (i >= 449 && i <= 476) ||

              (i >= 484 && i <= 511) || (i >= 513 && i <= 540))

      NS_LOG_DEBUG ("Added signal power to subbands 36-63, 65-92, 100-127, 129-156, 164-191, 193-220, 228-255, 257-284, 292-319, 321-348, 356-383, 385-412, 420-447, 449-476, 484-511, and 513-540");

      NS_LOG_DEBUG ("Added signal power to subbands " << start1 << "-" << stop1 <<

      NS_LOG_DEBUG ("Integrated power " << Integral (*c));

                    " and " << start2 << "-" << stop2);

      NS_ASSERT_MSG (std::abs (txPowerW - Integral (*c)) < 1e-6, "Power allocation failed"); 

      break;

    case 160:

      // 2 x 996 subcarriers (2 x 80 MHZ bands)

      txPowerPerBand = txPowerW / (2 * 996);

      NS_LOG_DEBUG ("Power per band " << txPowerPerBand);

      start1 = (nGuardBands / 2) + 12;

      stop1 = start1 + 498 - 1;

      start2 = stop1 + 6;

      stop2 = start2 + 498 - 1;

      start3 = stop2 + (2 * 12);

      stop3 = start3 + 498 - 1;

      start4 = stop3 + 6;

      stop4 = start4 + 498 - 1;

      for (size_t i = 0; i < c->GetSpectrumModel ()->GetNumBands (); i++, vit++, bit++)

        {

          if ((i >= start1 && i <= stop1) || (i >= start2 && i <= stop2) ||

              (i >= start3 && i <= stop3) || (i >= start4 && i <= stop4))

            {

              *vit = txPowerPerBand / (bit->fh - bit->fl);

            }

          else

            {

              *vit = 0;

            }

        }

      NS_LOG_DEBUG ("Added signal power to subbands " << start1 << "-" << stop1 <<

                    ", " << start2 << "-" << stop2 <<

                    ", " << start3 << "-" << stop3 <<

                    ", " << start4 << "-" << stop4);

WifiSpectrumValueHelper::CreateHeOfdmTxPowerSpectralDensity (uint32_t centerFrequency, uint32_t channelWidth, double txPowerW)

WifiSpectrumValueHelper::CreateNoisePowerSpectralDensity (uint32_t centerFrequency, uint32_t channelWidth, double bandBandwidth, double noiseFigure, uint32_t guardBandwidth)

  NS_LOG_FUNCTION (centerFrequency << channelWidth << txPowerW);

  Ptr<SpectrumModel> model = GetSpectrumModel (centerFrequency, channelWidth, bandBandwidth, guardBandwidth);

  Ptr<SpectrumValue> c = Create<SpectrumValue> (GetSpectrumModel (centerFrequency, channelWidth, 78125));

  Values::iterator vit = c->ValuesBegin ();

  Bands::const_iterator bit = c->ConstBandsBegin ();

  double txPowerPerBand;

  switch (channelWidth)

    {

    case 20:

      // 242 subcarriers (234 data + 8 pilot)

      NS_ASSERT_MSG (c->GetSpectrumModel ()->GetNumBands () == 513, "Unexpected number of bands");

      // skip 133 subbands, then place power in 121 subbands, then

      // skip 3 DC, then place power in 121 subbands, then skip

      // the final 133 subbands.

      txPowerPerBand = txPowerW / 242;

      for (size_t i = 0; i < c->GetSpectrumModel ()->GetNumBands (); i++, vit++, bit++)

        {

          if ((i >= 134 && i <= 254) || (i >= 258 && i <= 378))

            {

              *vit = txPowerPerBand / (bit->fh - bit->fl);

            }

        }

      NS_LOG_DEBUG ("Added signal power to subbands 134-254 and 258-378");

      NS_LOG_DEBUG ("Integrated power " << Integral (*c));

      NS_ASSERT_MSG (std::abs (txPowerW - Integral (*c)) < 1e-6, "Power allocation failed"); 

      break;

    case 40:

      // 484 subcarriers (468 data + 16 pilot)

      NS_ASSERT_MSG (c->GetSpectrumModel ()->GetNumBands () == 769, "Unexpected number of bands");

      // skip 139 subbands, then place power in 242 subbands, then

      // skip 5 DC, then place power in 242 subbands, then skip

      // the final 139 subbands.

      txPowerPerBand = txPowerW / 484;

      for (size_t i = 0; i < c->GetSpectrumModel ()->GetNumBands (); i++, vit++, bit++)

        {

          if ((i >= 140 && i <= 381) || (i >= 387 && i <= 628))

            {

              *vit = txPowerPerBand / (bit->fh - bit->fl);

            }

        }

      NS_LOG_DEBUG ("Added signal power to subbands 140-381 and 387-628");

      NS_LOG_DEBUG ("Integrated power " << Integral (*c));

      NS_ASSERT_MSG (std::abs (txPowerW - Integral (*c)) < 1e-6, "Power allocation failed"); 

      break;

    case 80:

      // 996 subcarriers (980 data + 16 pilot)

      NS_ASSERT_MSG (c->GetSpectrumModel ()->GetNumBands () == 1281, "Unexpected number of bands");

      // skip 139 subbands, then place power in 498 subbands, then

      // skip 5 DC, then place power in 498 subbands, then skip

      // the final 139 subbands.

      txPowerPerBand = txPowerW / 996;

      for (size_t i = 0; i < c->GetSpectrumModel ()->GetNumBands (); i++, vit++, bit++)

        {

          if ((i >= 140 && i <= 637) || (i >= 643 && i <= 1140))

            {

              *vit = txPowerPerBand / (bit->fh - bit->fl);

            }

        }

      NS_LOG_DEBUG ("Added signal power to subbands 140-637 and 643-1140");

      NS_LOG_DEBUG ("Integrated power " << Integral (*c));

      NS_ASSERT_MSG (std::abs (txPowerW - Integral (*c)) < 1e-6, "Power allocation failed"); 

      break;

    case 160:

      // 2 x 996 subcarriers (2 x 80 MHZ bands)

      NS_ASSERT_MSG (c->GetSpectrumModel ()->GetNumBands () == 2305, "Unexpected number of bands");

      txPowerPerBand = txPowerW / (2 * 996);

      for (size_t i = 0; i < c->GetSpectrumModel ()->GetNumBands (); i++, vit++, bit++)

        {

          if ((i >= 140 && i <= 637) || (i >= 643 && i <= 1140) ||

              (i >= 1164 && i <= 1661) || (i >= 1667 && i <= 2164))

            {

              *vit = txPowerPerBand / (bit->fh - bit->fl);

            }

        }

      NS_LOG_DEBUG ("Added signal power to subbands 140-637, 643-1140, 1164-1661 and 1667-2164");

      NS_LOG_DEBUG ("Integrated power " << Integral (*c));

      NS_ASSERT_MSG (std::abs (txPowerW - Integral (*c)) < 1e-6, "Power allocation failed"); 

      break;

    default:

      NS_FATAL_ERROR ("ChannelWidth " << channelWidth << " unsupported");

      break;

    }

  return c;

}

Ptr<SpectrumValue>

WifiSpectrumValueHelper::CreateNoisePowerSpectralDensity (uint32_t centerFrequency, uint32_t channelWidth, double bandBandwidth, double noiseFigure)

{

  Ptr<SpectrumModel> model = GetSpectrumModel (centerFrequency, channelWidth, bandBandwidth);

WifiSpectrumValueHelper::CreateRfFilter (uint32_t centerFrequency, uint32_t channelWidth, double bandGranularity)

WifiSpectrumValueHelper::CreateRfFilter (uint32_t centerFrequency, uint32_t channelWidth, double bandGranularity, uint32_t guardBandwidth)

  NS_LOG_FUNCTION (centerFrequency << channelWidth << bandGranularity);

  NS_LOG_FUNCTION (centerFrequency << channelWidth << bandGranularity << guardBandwidth);

  Ptr<SpectrumValue> c = Create <SpectrumValue> (GetSpectrumModel (centerFrequency, channelWidth, bandGranularity));

  Ptr<SpectrumValue> c = Create <SpectrumValue> (GetSpectrumModel (centerFrequency, channelWidth, bandGranularity, guardBandwidth));

  static Ptr<SpectrumModel> GetSpectrumModel (uint32_t centerFrequency, uint32_t channelWidth, double bandBandwidth);

  static Ptr<SpectrumModel> GetSpectrumModel (uint32_t centerFrequency, uint32_t channelWidth, double bandBandwidth, uint32_t guardBandwidth);

  /**

   * Create a transmit power spectral density corresponding to DSSS 

   *

   * The center frequency typically corresponds to 802.11b channel 

   * center frequencies but is not restricted to those frequencies.

   *

   * \param centerFrequency center frequency (MHz)

   * \param txPowerW  transmit power (W) to allocate

   * \param guardBandwidth width of the guard band (MHz)

   */

  static Ptr<SpectrumValue> CreateDsssTxPowerSpectralDensity (uint32_t centerFrequency, double txPowerW, uint32_t guardBandwidth);

   * High Efficiency (HE) (802.11ax).  Channel width may vary between 

   * (802.11a/g).  Channel width may vary between 20, 10, and 5 MHz.

   * 20, 40, 80, and 160 MHz.

  static Ptr<SpectrumValue> CreateHeOfdmTxPowerSpectralDensity (uint32_t centerFrequency, uint32_t channelWidth, double txPowerW);

  static Ptr<SpectrumValue> CreateOfdmTxPowerSpectralDensity (uint32_t centerFrequency, uint32_t channelWidth, double txPowerW, uint32_t guardBandwidth);

  static Ptr<SpectrumValue> CreateHtOfdmTxPowerSpectralDensity (uint32_t centerFrequency, uint32_t channelWidth, double txPowerW);

  static Ptr<SpectrumValue> CreateHtOfdmTxPowerSpectralDensity (uint32_t centerFrequency, uint32_t channelWidth, double txPowerW, uint32_t guardBandwidth);

   * (802.11a/g).  Channel width may vary between 20, 10, and 5 MHz.

   * High Efficiency (HE) (802.11ax).  Channel width may vary between 

   * 20, 40, 80, and 160 MHz.

  static Ptr<SpectrumValue> CreateOfdmTxPowerSpectralDensity (uint32_t centerFrequency, uint32_t channelWidth, double txPowerW);

  static Ptr<SpectrumValue> CreateHeOfdmTxPowerSpectralDensity (uint32_t centerFrequency, uint32_t channelWidth, double txPowerW, uint32_t guardBandwidth);

  /**

   * Create a transmit power spectral density corresponding to DSSS 

   *

   * The center frequency typically corresponds to 802.11b channel 

   * center frequencies but is not restricted to those frequencies.

   *

   * \param centerFrequency center frequency (MHz)

   * \param txPowerW  transmit power (W) to allocate

   */

  static Ptr<SpectrumValue> CreateDsssTxPowerSpectralDensity (uint32_t centerFrequency, double txPowerW);

  static Ptr<SpectrumValue> CreateNoisePowerSpectralDensity (uint32_t centerFrequency, uint32_t channelWidth, double bandBandwidth, double noiseFigure);

  static Ptr<SpectrumValue> CreateNoisePowerSpectralDensity (uint32_t centerFrequency, uint32_t channelWidth, double bandBandwidth, double noiseFigure, uint32_t guardBandwidth);

  static Ptr<SpectrumValue> CreateRfFilter (uint32_t centerFrequency, uint32_t channelWidth, double bandGranularity);

  static Ptr<SpectrumValue> CreateRfFilter (uint32_t centerFrequency, uint32_t channelWidth, double bandGranularity, uint32_t guardBandwidth);

          m_rxSpectrumModel = WifiSpectrumValueHelper::GetSpectrumModel (GetFrequency (), GetChannelWidth (), GetBandBandwidth ());

          m_rxSpectrumModel = WifiSpectrumValueHelper::GetSpectrumModel (GetFrequency (), GetChannelWidth (), GetBandBandwidth (), GetGuardBandwidth ());

  Ptr<SpectrumValue> filter = WifiSpectrumValueHelper::CreateRfFilter (GetFrequency (), GetChannelWidth (), GetBandBandwidth ());

  Ptr<SpectrumValue> filter = WifiSpectrumValueHelper::CreateRfFilter (GetFrequency (), GetChannelWidth (), GetBandBandwidth (), GetGuardBandwidth ());

      v = WifiSpectrumValueHelper::CreateOfdmTxPowerSpectralDensity (centerFrequency, channelWidth, txPowerW);

      v = WifiSpectrumValueHelper::CreateOfdmTxPowerSpectralDensity (centerFrequency, channelWidth, txPowerW, GetGuardBandwidth ());

      v = WifiSpectrumValueHelper::CreateDsssTxPowerSpectralDensity (centerFrequency, txPowerW);

      v = WifiSpectrumValueHelper::CreateDsssTxPowerSpectralDensity (centerFrequency, txPowerW, GetGuardBandwidth ());

      v = WifiSpectrumValueHelper::CreateHtOfdmTxPowerSpectralDensity (centerFrequency, channelWidth, txPowerW);

      v = WifiSpectrumValueHelper::CreateHtOfdmTxPowerSpectralDensity (centerFrequency, channelWidth, txPowerW, GetGuardBandwidth ());

      v = WifiSpectrumValueHelper::CreateHeOfdmTxPowerSpectralDensity (centerFrequency, channelWidth, txPowerW);

      v = WifiSpectrumValueHelper::CreateHeOfdmTxPowerSpectralDensity (centerFrequency, channelWidth, txPowerW, GetGuardBandwidth ());

  return bandBandwidth ;

  return bandBandwidth;

}

uint32_t

SpectrumWifiPhy::GetGuardBandwidth (void) const

{

  double guardBandwidth = 0;

  switch (GetStandard ())

    {

    case WIFI_PHY_STANDARD_80211a:

    case WIFI_PHY_STANDARD_80211g:

    case WIFI_PHY_STANDARD_holland:

    case WIFI_PHY_STANDARD_80211_10MHZ:

    case WIFI_PHY_STANDARD_80211_5MHZ:

    case WIFI_PHY_STANDARD_80211b:

    case WIFI_PHY_STANDARD_80211n_2_4GHZ:

    case WIFI_PHY_STANDARD_80211n_5GHZ:

    case WIFI_PHY_STANDARD_80211ac:

    case WIFI_PHY_STANDARD_80211ax_2_4GHZ:

    case WIFI_PHY_STANDARD_80211ax_5GHZ:

      // Use 10 MHZ

      guardBandwidth = 10;

      break;

    default:

      NS_FATAL_ERROR ("Standard unknown: " << GetStandard ());

      break;

    }

  return guardBandwidth;

  Ptr<SpectrumValue> txPowerSpectrum = WifiSpectrumValueHelper::CreateOfdmTxPowerSpectralDensity (FREQUENCY, CHANNEL_WIDTH, txPowerWatts);

  Ptr<SpectrumValue> txPowerSpectrum = WifiSpectrumValueHelper::CreateOfdmTxPowerSpectralDensity (FREQUENCY, CHANNEL_WIDTH, txPowerWatts, GUARD_WIDTH);