---

sources.cc

Go to the documentation of this file.

// sources.cc
// SuperMix version 1.0  C++ source file
//
// Copyright (c) 1999 California Institute of Technology.
// All rights reserved.
//
// Redistribution and use in source and binary forms for noncommercial
// purposes are permitted provided that the above copyright notice and
// this paragraph are duplicated in all such forms and that any
// documentation and other materials related to such distribution and
// use acknowledge that the software was developed by California
// Institute of Technology. Redistribution and/or use in source or
// binary forms is not permitted for any commercial purpose. Use of
// this software does not include a permitted use of the Institute's
// name or trademark for any purpose.
//
// DISCLAIMER:
// THIS SOFTWARE AND/OR RELATED MATERIALS ARE PROVIDED "AS-IS" WITHOUT
// WARRANTY OF ANY KIND INCLUDING ANY WARRANTIES OF PERFORMANCE OR
// MERCHANTABILITY OR FITNESS FOR A PARTICULAR USE OR PURPOSE (AS SET
// FORTH IN UCC 23212-2313) OR FOR ANY PURPOSE WHATSOEVER, FOR THE
// LICENSED PRODUCT, HOWEVER USED.  IN NO EVENT SHALL CALTECH/JPL BE
// LIABLE FOR ANY DAMAGES AND/OR COSTS, INCLUDING BUT NOT LIMITED TO
// INCIDENTAL OR CONSEQUENTIAL DAMAGES OF ANY KIND, INCLUDING ECONOMIC
// DAMAGE OR INJURY TO PROPERTY AND LOST PROFITS, REGARDLESS OF
// WHETHER CALTECH/JPL SHALL BE ADVISED, HAVE REASON TO KNOW, OR IN
// FACT SHALL KNOW OF THE POSSIBILITY.  THE USER BEARS ALL RISK
// RELATING TO QUALITY AND PERFORMANCE OF THE SOFTWARE AND/OR RELATED
// MATERIALS.
//
// Change history:
// 10/29/99: Generator R now shadows device::Z0 by default
// 11/11/98: Changed table access to new syntax
// 11/10/98: Removed limits on voltage and current in
//           voltage_source and current_sink
// 9/18/98:  Made constructors more efficient.
// 9/15/98:  Fixed a nasty bug: wrong argument order in
//           passive noise calculations.
// 12/29/97: modified by FR for new sdata (with znorm)
// 12/19/97: modified by FR for matmath classes

#include "sources.h"
#include <math.h>
#include "error.h"


// generator:

generator::generator() :
  nport(1), 
  R(&device::Z0), Temp(&device::T), 
  source_f(0.0), source_width(0.0), source_power(0.0), source_phase(0.0)
{
  R.set_min(0.0);
  source_f.set_min(0.0);
  source_width.set_min(0.0);
  source_power.set_min(0.0);
}

const nport::data_info & generator::get_data_info()
{ 
  info.active = (Temp != device::T && R != 0.0);
  info.source = source_power != 0.0;
  return info;
}

void generator::recalc_S()
{
  // Set normalization impedance. added 12/29/97
  data.set_znorm(Z0);

  // Calculate S11
  data.S[1][1] = (R - Z0)/(R + Z0);

  // Convert the source power to an amplitude wave and put in B1
  if(fabs(f - source_f) <= (source_width/2.0))
    data.B[1] = polar(2.0*sqrt(Z0*R*source_power)/(Z0+R), source_phase);
  else
    data.B[1] = 0.0;
}

void generator::recalc()
{
  recalc_S();

  // Calculate C11
  compute_passive_noise(data, f, Temp);
}


// voltage_source:

voltage_source::voltage_source() :
  nport(1),
  R(0.0), Temp(&T),
  source_f(0.0), source_width(0.0), source_voltage(0.0), source_phase(0.0)
{
  R.set_min(0.0);
  source_f.set_min(0.0);
  source_width.set_min(0.0);
}

const nport::data_info & voltage_source::get_data_info()
{ 
  info.active = (Temp != device::T && R != 0.0);
  info.source = source_voltage != 0.0;
  return info;
}

void voltage_source::recalc_S()
{
  // Set normalization impedance:
  data.set_znorm(Z0);

  // Calculate S11
  data.S[1][1] = (R - Z0)/(R + Z0);

  // Convert the source voltage to an amplitude wave and put in B1
  if(fabs(f - source_f) <= (source_width/2.0))
    data.B[1] = polar(sqrt(Z0)*source_voltage/(Z0+R), source_phase);
  else
    data.B[1] = 0.0;
}

void voltage_source::recalc()
{
  recalc_S();

  // Calculate C11
  compute_passive_noise(data, f, Temp);
}


// current_sink:

current_sink::current_sink() :
  nport(1),
  Y(0.0), Temp(&T),
  source_f(0.0), source_width(0.0), sink_current(0.0), sink_phase(0.0)
{
  Y.set_min(0.0);
  source_f.set_min(0.0);
  source_width.set_min(0.0);
}

const nport::data_info & current_sink::get_data_info()
{ 
  info.active = (Temp != device::T && Y != 0.0);
  info.source = sink_current != 0.0;
  return info;
}

void current_sink::recalc_S()
{
  // Set normalization impedance:
  data.set_znorm(Z0);

  // Calculate S11
  data.S[1][1] = (1 - Y*Z0)/(1 + Y*Z0);

  // Calculate C11
  compute_passive_noise(data, f, Temp);

  // Convert the sink current to an amplitude wave and put in B1
  if(fabs(f - source_f) <= (source_width/2.0))
    data.B[1] = - polar(sqrt(Z0)*sink_current/(1 + Y*Z0), sink_phase);
  else
    data.B[1] = 0.0;
}

void current_sink::recalc()
{
  recalc_S();

  // Calculate C11
  compute_passive_noise(data, f, Temp);
}

---

Please direct comments and corrections to supermix@submm.caltech.edu
Go to the supermix home page
Generated by 1.2.7