---

transformer.cc

Go to the documentation of this file.

// transformer.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:  Default is now to shadow device::Z0
// 9/27/99:   Added verbosity control for warning about impedance
// 11/11/98:  Changed table access to new syntax
// 9/18/98:   made constructors more efficient.
// 7/28/98:   added <math.h> to includes, since SIScmplx.h doesn't
// 12/29/97:  modified by FR for new sdata (with znorm)
// 12/19/97:  modified by FR for matmath classes

#include "transformer.h"
#include "global.h"
#include "error.h"
#include <math.h>

//
// Default constructor
//
transformer::transformer() : 
  nport(2), Z1(&device::Z0), Z2(&device::Z0), is_verbose(false)
{ info.noise = info.active = info.source = false; }

void transformer::recalc()
{
  double znorm = device::Z0 ;
  Complex zz1 = Z1 ;
  Complex zz2 = Z2 ;
  double rr1 = real(zz1)/znorm ;
  double xx1 = imag(zz1)/znorm ;
  if(rr1 < Tiny) {
    if(is_verbose)
      error::warning(
      "Very small or negative input resistance for transformer") ;
    rr1 = Tiny;
  }

  double rr2 = real(zz2)/znorm ;
  double xx2 = imag(zz2)/znorm ;
  if(rr2 < Tiny) {
    if(is_verbose)
      error::warning(
      "Very small or negative output resistance for transformer") ;
    rr2 = Tiny;
  }

  double sr = sqrt(rr1/rr2) ;
  double zeta = 0.5*(sr - 1./sr) ;
  double psi = 0.5*(xx2*sr + xx1/sr) ;
  complex denom = sqrt(1.+zeta*zeta) + I*psi;

  complex S11 = (zeta + I*psi)/denom ;
  complex S22 = (-zeta + I*psi)/denom ;
  complex S12 = 1./denom ;  // = S21 from reciprocity

  data.set_znorm(znorm);  // which is Z0.  added 12/29/97.
  data.S[1][1] = S11;
  data.S[2][2] = S22;
  data.S[1][2] = S12;
  data.S[2][1] = S12;

  // Assume that data.C is already a zero matrix.
}

---

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