From jfurman-dsl@sbcglobal.net Sun Apr 27 03:17:41 2008
Date: Sun, 27 Apr 2008 03:14:00 -0700
From: Jeffrey Furman <jfurman-dsl@sbcglobal.net>
To: jeff.furman@gmail.com
Subject: smplot.C


 smplot.C


    [ Part 2: "Attached Text" ]

/* generate smith chart data to plot */


#define gcc 0
#if gcc
#include <stdio.h>
#include <math.h>
#else
#include "..\..\include\math.h"
#include "..\..\include\stdio.h"
#endif

#define p 20
#define f1 (3e6)            

#define f2 f1
#define fr (1.02)

#define BASE 2.0
#define NUMBER_OF_LOOPS 5

#define plot 1

#define sgn(x) ((x == 0)? 1 :((x > 0)?1 :-1)))

#define twopi (2*3.1415926535) 

#define DELTAT 0.01

/*
double r_list[] = {0.1,0.2,0.3,0.4,0.5,0.6,0.7,0.8,0.9,
                   1.0,1.2,1.4,1.6,1.8,
                   2.0,2.5,3.0,4.0,5.0,6.0,7.0,8.0,9.0,
                   10.0,15.0,20.0,50.0,};
                   */
double r_list[] = {0.2,0.4,0.6,0.8,
                   1.0,1.5,
                   2.0,3.0,5.0,7.0,9.0,
                   10.0,15.0,20.0,50.0,};

#define rcount (sizeof(r_list)/sizeof(double))

raxis(r)
float r;
{
float t;
printf("\t \ta-%g\n",r);
for(t = 0;t <= 2*r; t = t+DELTAT)
printf("%.10g\t  %.10g\n",t,0);

for(t = 2*r;t >= 0; t = t-DELTAT)
printf("%.10g\t  %.10g\n",t,0);
}

rcircle(r)
float r;
{
float t;
printf("\t \tr-%g\n",r);
for(t = 0;t <= 2*r; t = t+0.02)
printf("%.10g\t  %.10g\n",t,sqrt(t*((2.0*r)-t)));

for(t = 2*r;t >= 0; t = t-DELTAT)
printf("%.10g\t  %.10g\n",t,-sqrt(t*((2.0*r)-t)));
}

xcircle(x)
float x;
{
float t,tmax;
printf("\t \tx-%g\n",x);
/*tmax = 2.0/(x*x+1.0);*/
tmax = x*(1.0+((1.0-x*x)/(x*x+1.0)));

for(t = 0;t <= tmax; t = t+DELTAT)
printf("%.10g\t  %.10g\n",sqrt(t*((2.0*x)-t)),t);

for(t = tmax;t >= 0; t = t-DELTAT)
printf("%.10g\t  %.10g\n",sqrt(t*((2.0*x)-t)),t);


for(t = 0;t <= tmax; t = t+DELTAT)
printf("%.10g\t  %.10g\n",sqrt(t*((2.0*x)-t)),-t);

for(t = tmax;t >= 0; t = t-DELTAT)
printf("%.10g\t  %.10g\n",sqrt(t*((2.0*x)-t)),-t);


}

smithchart()

{
double ratio,val,g1,g2;
int i,j;
ratio = exp(log(BASE)/NUMBER_OF_LOOPS);
printf("\t \tratio= %g\n",ratio);
rcircle(1.0);
raxis(1.0);
for(i=0;i<rcount;i++){
g1 =(r_list[i] -1.0)/(r_list[i] + 1.0);
g2 = (1.0 - g1)/2.0;
printf("\t \tg1 = %g, g2 = %g\n",g1,g2);
rcircle(g2);
/*rcircle(0.5/val);*/
/*rcircle(val);*/
xcircle((1.0-g1)/(1.0+g1));
/*xcircle(val);*/
/*
rcircle(val/2.0);
rcircle(2.0/val);
*/
};
}



float l(turns)
float turns;
{
#define L0  (1.0e-7)
#define L1  (2.653e-8)
#define L2  (1.313e-7)
#define T0  (5.75)

return (L0 + (L1*turns*turns) + (L2*turns));

}

float c(turns)
float turns;
{
#define C0  (15.e-12)
#define C1  ((3500.e-12)-C0)

return (C0 + (C1*(1.0+(((2*(turns-floor(turns)))-1)*(sgn(sin(twopi*(turns-floor(turns))))))));

}

float z1r(r,xs,xp)
float r,xs,xp;
{
return (float)((xp*xp*r)/((r*r)+((xs+xp)*(xs+xp))));
}

float z1i(r,xs,xp)
float r,xs,xp;
{
return (float)((xp*((r*r)+(xs*(xs+xp))))/((r*r)+((xs+xp)*(xs+xp))));
}

float z2r(r,xp,xs)
float r,xs,xp;
{
return (float)(xp*xp*r/((r*r)+(xp*xp)));
}

float z2i(r,xp,xs)
float r,xs,xp;
{
return (float)(xs+(xp*r*r/((r*r)+(xp*xp))));
}



float xl(f,t)
float f,t;
{
float tmp;
tmp = l(t);
return (float)(twopi*f*tmp);
}


float xc(f,t)
float f,t;
{
float tmp;
tmp = c(t);
return ((float)((-1.0)/(float)((float)twopi*f*(float)tmp)));
/*return (1.0/(float)c(t));*/
}

#define r0 (50.0)
#define f0 (1.0e6)

main()
{
int i,j; float t,f,r,x,xt,yt,d;




smithchart();


}