{
"metadata": {
"name": "",
"signature": "sha256:ffa36c37285243e82f8b7c3abbc469b9a273b2bd68fa586a9dc529ef122af349"
},
"nbformat": 3,
"nbformat_minor": 0,
"worksheets": [
{
"cells": [
{
"cell_type": "heading",
"level": 1,
"metadata": {},
"source": [
"Ch-4: Microwave resonators & Waveguide components"
]
},
{
"cell_type": "heading",
"level": 2,
"metadata": {},
"source": [
"Page Number: 193 Example 4.1"
]
},
{
"cell_type": "code",
"collapsed": false,
"input": [
"from __future__ import division\n",
"from math import sqrt \n",
"#Given\n",
"a=5 #cm\n",
"a1=a/100 #m\n",
"b=2 #cm\n",
"b1=b/100 #m\n",
"c=15 #cm\n",
"c1=c/100 #m\n",
"\n",
"#(i) Air filled cavity\n",
"m=1 \n",
"n=0 \n",
"p=1 \n",
"c=3e8 #for air\n",
"fr=(1/2)*c*sqrt((m/a1)**2+(n/b1)**2+(p/c1)**2) #hz\n",
"print 'Resonant frequency for an air filled cavity: %0.3f'%(fr/10**9),'Ghz'\n",
"\n",
"#(ii) Dielctric filled cavity\n",
"er=2.56 \n",
"fr1=(1/2)*(c/sqrt(er))*sqrt((m/a1)**2+(n/b1)**2+(p/c1)**2) #hz\n",
"print 'Resonant frequency for dielectric cavity: %0.3f'%(fr1/10**9),'Ghz'"
],
"language": "python",
"metadata": {},
"outputs": [
{
"output_type": "stream",
"stream": "stdout",
"text": [
"Resonant frequency for an air filled cavity: 3.162 Ghz\n",
"Resonant frequency for dielectric cavity: 1.976 Ghz\n"
]
}
],
"prompt_number": 30
},
{
"cell_type": "heading",
"level": 2,
"metadata": {},
"source": [
"Page Number: 193 Example 4.2"
]
},
{
"cell_type": "code",
"collapsed": false,
"input": [
" \n",
"#Given\n",
"a=0.38 #cm\n",
"a1=a/100 #m\n",
"b=0.76 #cm\n",
"b1=b/100 #m\n",
"f=50e9 \n",
"c=3e8 \n",
"\n",
"#Length for TE102\n",
"m=1 \n",
"n=0 \n",
"p=2 \n",
"l=1/sqrt((f/c)**2-(1/(4*b1**2))) #m\n",
"print 'Length c: %0.3f'%(l*100),'cm'"
],
"language": "python",
"metadata": {},
"outputs": [
{
"output_type": "stream",
"stream": "stdout",
"text": [
"Length c: 0.653 cm\n"
]
}
],
"prompt_number": 31
},
{
"cell_type": "heading",
"level": 2,
"metadata": {},
"source": [
"Page Number: 194 Example 4.3"
]
},
{
"cell_type": "code",
"collapsed": false,
"input": [
"from math import pi \n",
"#Given\n",
"c=3e8 #m/s\n",
"a=2.286 #cm\n",
"a1=a/100 #m\n",
"b=1.024 #cm\n",
"b1=b/100 #m\n",
"f=10e9 #hz\n",
"sig=6e7 \n",
"u=4e-7*pi \n",
"w=2*pi*f \n",
"eet=377 \n",
"\n",
"#Shortest cavity length\n",
"lamc=2*a1 #m\n",
"fc=c/lamc #hz\n",
"lam=c/f #m\n",
"lamg=lam/sqrt(1-(fc/f)**2) #m\n",
"sc=lamg/2 #m\n",
"print 'Shortest cavity length: %0.3f'%(sc*100), 'cm'\n",
"\n",
"#Qw of the resonator operating in TE101 mode\n",
"rs=sqrt((w*u)/(2*sig)) #ohm\n",
"lamr=c/f \n",
"x=(((a1*b1)/(sc**2))+((sc**2+a1**2)/(2*sc*a1))+(b1*sc/a1**2)) \n",
"qw=(2*pi*eet*a1*b1*sc)/(rs*(lamr**3)*x) \n",
"print 'Qw of the resonator operating in TE101 mode %0.3f'%qw"
],
"language": "python",
"metadata": {},
"outputs": [
{
"output_type": "stream",
"stream": "stdout",
"text": [
"Shortest cavity length: 1.988 cm\n",
"Qw of the resonator operating in TE101 mode 7990.324\n"
]
}
],
"prompt_number": 32
},
{
"cell_type": "heading",
"level": 2,
"metadata": {},
"source": [
"Page Number: 195 Example 4.4"
]
},
{
"cell_type": "code",
"collapsed": false,
"input": [
" \n",
"#Given\n",
"c=3e8 #m/s\n",
"a=4.8 #cm\n",
"a1=a/100 #m\n",
"b=2.2 #cm\n",
"b1=b/100 #m\n",
"f=5e9 #hz\n",
"er=2.25 \n",
"tandel=4e-4 \n",
"sig=5.813e7 \n",
"oneby=3e8 \n",
"u=4e-7*pi \n",
"w=2*pi*f \n",
"eet=377 \n",
"\n",
"#Length at p=1\n",
"m=1 \n",
"n=0 \n",
"p=1 \n",
"z=(f*2*sqrt(er))/c \n",
"cp1=p/sqrt((z**2)-((m/a1)**2)-((n/b1)**2)) \n",
"print 'Length of resonator at p=1: %0.3f'%(cp1*100), 'cm'\n",
"\n",
"#At p=2\n",
"cp2=cp1*2 \n",
"print 'Length of resonator at p=2: %0.3f'%(cp2*100), 'cm'\n",
"\n",
"#Qw\n",
"rs=sqrt((w*u)/(2*sig)) #ohm\n",
"lamr=c/(f*sqrt(er)) \n",
"x=(((a1*b1)/(cp1**2))+((cp1**2+a1**2)/(2*cp1*a1))+(b1*cp1/a1**2)) \n",
"qw=(2*pi*(eet/sqrt(er))*a1*b1*cp1)/(rs*(lamr**3)*x) \n",
"qd=1/tandel \n",
"q=(qw*qd)/(qw+qd) \n",
"print 'Q for TE101 mode: %0.3f'%q"
],
"language": "python",
"metadata": {},
"outputs": [
{
"output_type": "stream",
"stream": "stdout",
"text": [
"Length of resonator at p=1: 2.200 cm\n",
"Length of resonator at p=2: 4.400 cm\n",
"Q for TE101 mode: 1925.612\n"
]
}
],
"prompt_number": 33
},
{
"cell_type": "heading",
"level": 2,
"metadata": {},
"source": [
"Page Number: 196 Example 4.5"
]
},
{
"cell_type": "code",
"collapsed": false,
"input": [
" \n",
"#Given\n",
"c=3e8 #m/s\n",
"a=2 #cm\n",
"a1=a/100 #m\n",
"b=2.5 #cm\n",
"b1=b/100 #m\n",
"\n",
"print 'TE modes' \n",
"h01=3.832 \n",
"fr=(c/(2*pi))*sqrt((h01/a1)**2+(pi/b1)**2) #hz\n",
"print 'Resonant frequency for mode TE010: %0.3f'%(fr/10**9),'Ghz'\n",
"\n",
"h11=1.841 \n",
"fr1=(c/(2*pi))*sqrt((h11/a1)**2+(pi/b1)**2) #hz\n",
"print 'Resonant frequency for mode TE111: %0.3f'%(fr1/10**9),'Ghz'\n",
"\n",
"h21=3.054 \n",
"fr2=(c/(2*pi))*sqrt((h21/a1)**2+(pi/b1)**2) #hz\n",
"print 'Resonant frequency for mode TE211: %0.3f'%(fr2/10**9),'Ghz'\n",
"\n",
"print 'TM modes:' \n",
"l1=0 \n",
"h011=2.405 \n",
"fr3=(c/(2*pi))*sqrt((h011/a1)**2+(pi*l1/b1)**2) #hz\n",
"print 'Resonant frequency for mode TM010 %0.3f'%(fr3/10**9),'Ghz'\n",
"\n",
"l2=1 \n",
"fr4=(c/(2*pi))*sqrt((h011/a1)**2+(pi*l2/b1)**2) #hz\n",
"print 'resonant frequency for mode TM011: %0.3f'%(fr4/10**9),'Ghz'\n",
"\n",
"l3=1 \n",
"h111=3.832 \n",
"fr5=(c/(2*pi))*sqrt((h111/a1)**2+(pi*l3/b1)**2) #hz\n",
"print 'Resonant frequency for mode TM111: %0.3f'%(fr5/10**9),'Ghz'"
],
"language": "python",
"metadata": {},
"outputs": [
{
"output_type": "stream",
"stream": "stdout",
"text": [
"TE modes\n",
"Resonant frequency for mode TE010: 10.940 Ghz\n",
"Resonant frequency for mode TE111: 7.438 Ghz\n",
"Resonant frequency for mode TE211: 9.442 Ghz\n",
"TM modes:\n",
"Resonant frequency for mode TM010 5.742 Ghz\n",
"resonant frequency for mode TM011: 8.305 Ghz\n",
"Resonant frequency for mode TM111: 10.940 Ghz\n"
]
}
],
"prompt_number": 34
},
{
"cell_type": "heading",
"level": 2,
"metadata": {},
"source": [
"Page Number: 196 Example 4.6"
]
},
{
"cell_type": "code",
"collapsed": false,
"input": [
" \n",
"#Given\n",
"QTM010=1.202 \n",
"QTE101=1.11 \n",
"\n",
"r=QTM010/QTE101 \n",
"print 'Ratio of Qs of cylindrical and rectangular resonators: %0.3f'%r"
],
"language": "python",
"metadata": {},
"outputs": [
{
"output_type": "stream",
"stream": "stdout",
"text": [
"Ratio of Qs of cylindrical and rectangular resonators: 1.083\n"
]
}
],
"prompt_number": 35
},
{
"cell_type": "heading",
"level": 2,
"metadata": {},
"source": [
"Page Number: 197 Example 4.7"
]
},
{
"cell_type": "code",
"collapsed": false,
"input": [
"#Given\n",
"f=7.07e9 #hz\n",
"a=3 #cm\n",
"a1=a/100 #m\n",
"sig=5.8e7 \n",
"er=2.25 \n",
"tandel=4e-4 \n",
"ur=1 \n",
"n=377 \n",
"w=2*pi*f \n",
"u=4e-7*pi \n",
"\n",
"#Q of resonantor\n",
"rs=sqrt(w*u/(2*sig)) #ohm\n",
"qw=(0.7419*n)/(rs*sqrt(2.25)) \n",
"qd=1/tandel \n",
"q=(qw*qd)/(qw+qd) \n",
"print 'Q of resonator: %0.3f'%q"
],
"language": "python",
"metadata": {},
"outputs": [
{
"output_type": "stream",
"stream": "stdout",
"text": [
"Q of resonator: 1931.819\n"
]
}
],
"prompt_number": 36
},
{
"cell_type": "heading",
"level": 2,
"metadata": {},
"source": [
"Page Number: 198 Example 4.8"
]
},
{
"cell_type": "code",
"collapsed": false,
"input": [
" \n",
"#Given\n",
"a=5 #cm\n",
"a1=a/100 #m\n",
"b=4 #cm\n",
"b1=b/100 #m\n",
"c=10 #cm\n",
"c1=c/100 #m\n",
"sig=5.8e7 \n",
"u0=4e-7*pi \n",
"er=3 \n",
"eet=377 \n",
"\n",
"ur=1 \n",
"spl=3e8 \n",
"tandel=2.5e-4 \n",
"\n",
"#TE101 mode\n",
"m=1 \n",
"n=0 \n",
"p=1 \n",
"fr=(spl/(2*sqrt(er*ur)))*sqrt((m/a1)**2+(n/b1)**2+(p/c1)**2) #hz\n",
"print 'Resonant frequency: %0.3f'%(fr/10**9), 'Ghz'\n",
"\n",
"w=2*pi*fr \n",
"rs=sqrt((w*u0)/(2*sig)) #ohm\n",
"lamr=spl/(fr*sqrt(er)) \n",
"x=(((a1*b1)/(c1**2))+((c1**2+a1**2)/(2*c1*a1))+((b1*c1)/a1**2)) \n",
"qw=(2*pi*(eet/sqrt(er))*a1*b1*c1)/(rs*(lamr**3)*x) \n",
"print 'Q for TE101 mode: %0.3f'%qw\n",
"\n",
"qd=1/tandel \n",
"q=(qw*qd)/(qw+qd) \n",
"print 'Q for lossy dielectric: %0.3f'%q\n",
"\n",
"#Value of qw is calculated wrong in book as lamr comes to be 0.08 not 0.89 m"
],
"language": "python",
"metadata": {},
"outputs": [
{
"output_type": "stream",
"stream": "stdout",
"text": [
"Resonant frequency: 1.936 Ghz\n",
"Q for TE101 mode: 10916.466\n",
"Q for lossy dielectric: 2927.360\n"
]
}
],
"prompt_number": 40
},
{
"cell_type": "heading",
"level": 2,
"metadata": {},
"source": [
"Page Number: 198 Example 4.9"
]
},
{
"cell_type": "code",
"collapsed": false,
"input": [
" \n",
"#Given\n",
"c=3e8 #m/s\n",
"a=2.286 #cm\n",
"a1=a/100 #m\n",
"b=1.106 #cm\n",
"b1=b/100 #m\n",
"\n",
"#For fr1=9.3e9 \n",
"fr1=9.3e9 #hz\n",
"lamr1=c/fr1 #m\n",
"c1=(2*a1)/sqrt((((2*a1)/lamr1)**2)-1) \n",
"\n",
"#For fr2=10.2e9 \n",
"fr2=10.2e9 #hz\n",
"lamr2=c/fr2 #m\n",
"c2=(2*a1)/sqrt((((2*a1)/lamr2)**2)-1) \n",
"\n",
"r=c1-c2 \n",
"print 'Range of piston movement: %0.3f'%(r*100), 'cm'"
],
"language": "python",
"metadata": {},
"outputs": [
{
"output_type": "stream",
"stream": "stdout",
"text": [
"Range of piston movement: 0.710 cm\n"
]
}
],
"prompt_number": 41
},
{
"cell_type": "heading",
"level": 2,
"metadata": {},
"source": [
"Page Number: 199 Example 4.10"
]
},
{
"cell_type": "code",
"collapsed": false,
"input": [
" \n",
"#Given\n",
"a=3 #cm\n",
"a1=a/100 #m\n",
"d=10 #cm\n",
"d1=d/100 #m\n",
"df=2.5e6 \n",
"er=2.25 \n",
"p11=1.841 \n",
"c=3e8 #m/s\n",
"\n",
"#Resonant frequency\n",
"fr=(c/2)*(sqrt((p11/a1)**2+(pi/d1)**2)) #hz\n",
"print 'Resonant frequency: %0.3f'%(fr/10**9),'Ghz'\n",
"\n",
"#Q without dielectric\n",
"q0=fr/df \n",
"print 'Q wirhout dielectric constant: %0.3f'%q0\n",
"\n",
"# Q with dielectric\n",
"fr1=fr/sqrt(er) \n",
"qd=1e3 \n",
"q=(q0*qd)/(q0+qd) \n",
"print 'Q with dielectric constant: %0.3f'%q"
],
"language": "python",
"metadata": {},
"outputs": [
{
"output_type": "stream",
"stream": "stdout",
"text": [
"Resonant frequency: 10.341 Ghz\n",
"Q wirhout dielectric constant: 4136.446\n",
"Q with dielectric constant: 805.313\n"
]
}
],
"prompt_number": 44
},
{
"cell_type": "heading",
"level": 2,
"metadata": {},
"source": [
"Page Number: 200 Example 4.11"
]
},
{
"cell_type": "code",
"collapsed": false,
"input": [
" \n",
"#Given\n",
"f=9.375e9 #hz\n",
"sig=5.8e7 \n",
"eet=377 \n",
"c=3e8 #m/s\n",
"w=2*pi*f \n",
"r=1.5 \n",
"u=4e-7*pi \n",
"\n",
"#Radius\n",
"a=c/(f*2.62) #m\n",
"print 'Radius of resonantor %0.3f'%(a*100), 'cm'\n",
"\n",
"#O\n",
"rs=sqrt((w*u)/(2*sig)) #ohm\n",
"x=1.202*eet \n",
"y=rs*(1+(1/r)) \n",
"q=x/y \n",
"print 'Q of the resonator: %0.3f'%q\n",
"\n",
"#Answer for Q is calculated as 10875 in book but it is 10763.303"
],
"language": "python",
"metadata": {},
"outputs": [
{
"output_type": "stream",
"stream": "stdout",
"text": [
"Radius of resonantor 1.221 cm\n",
"Q of the resonator: 10763.303\n"
]
}
],
"prompt_number": 45
},
{
"cell_type": "heading",
"level": 2,
"metadata": {},
"source": [
"Page Number: 215 Example 4.12"
]
},
{
"cell_type": "code",
"collapsed": false,
"input": [
" \n",
"#Given\n",
"f=5e9 #hz\n",
"sig=5.813e7 \n",
"er=2.25 \n",
"tandel=4e-4 \n",
"c=3e8 #m/s\n",
"h01=3.832 \n",
"u=4e-7*pi \n",
"\n",
"#Length of resonator\n",
"lamr=c/(f*sqrt(er)) \n",
"d=sqrt((((((2*3.832)**2)+(pi*pi))*(lamr*lamr))/(2*2*pi*pi)) )\n",
"print 'Length of resonator: %0.3f'%(d*100),'cm'\n",
"\n",
"#Q of resonator\n",
"n=(120*pi)/sqrt(er) \n",
"Rs=sqrt((f*u)/sig) \n",
"a=d/2 \n",
"Qw1=n*(((h01/a)**2+(pi/d)**2)**(3/2))\n",
"Qw2=2*Rs*(((h01*h01)/(a*a*a))+((2*pi*pi)/(d*d*d))) \n",
"Qw=Qw1/Qw2 \n",
"Qd=1/tandel \n",
"Q=(Qw*Qd)/(Qw+Qd) \n",
"print 'Q of resonator: %0.3f'%Q\n",
"\n",
"#Value of Qw is calculated wrong in the book, it should be 50057.91 instead of 53473.8\n",
"#Hence the value of Q also differs"
],
"language": "python",
"metadata": {},
"outputs": [
{
"output_type": "stream",
"stream": "stdout",
"text": [
"Length of resonator: 5.273 cm\n",
"Q of resonator: 2381.084\n"
]
}
],
"prompt_number": 46
},
{
"cell_type": "heading",
"level": 2,
"metadata": {},
"source": [
"Page Number: 215 Example 4.13"
]
},
{
"cell_type": "code",
"collapsed": false,
"input": [
"#Given\n",
"p=100 #mW\n",
"#As 2 and 3 are matched terminals\n",
"x=1/2 \n",
"y=1/sqrt(2) \n",
"s=[x,-x, y, -x, 0, y, y, y, 0] \n",
"\n",
"#Power delivered\n",
"#Port 1\n",
"p1=p*(1-s[1]**2) \n",
"print 'Power at port 1:',p1, 'mW'\n",
"\n",
"#Port2\n",
"p2=p*s[2]**2 \n",
"print 'Power at port 2:',p2, 'mW'\n",
"\n",
"#Port 3\n",
"p3=p*s[3]**2 \n",
"print 'Power at port 3:' ,p3,'mW'"
],
"language": "python",
"metadata": {},
"outputs": [
{
"output_type": "stream",
"stream": "stdout",
"text": [
"Power at port 1: 75.0 mW\n",
"Power at port 2: 50.0 mW\n",
"Power at port 3: 25.0 mW\n"
]
}
],
"prompt_number": 14
},
{
"cell_type": "heading",
"level": 2,
"metadata": {},
"source": [
"Page Number: 216 Example 4.14"
]
},
{
"cell_type": "code",
"collapsed": false,
"input": [
" \n",
"#Given\n",
"p=40 #mW\n",
"#Since port 3 is matched\n",
"x=sqrt(2) \n",
"s=[1, 1, x, 1, 1, -x, x, -x, 0] \n",
"r1=40 #ohm\n",
"r2=60 #ohm\n",
"w=50 #ohm\n",
"\n",
"#Reflection coefficients\n",
"T1=(w-r1)/(w+r1) \n",
"T2=(r2-w)/(r2+w) \n",
"\n",
"#As power is fed into 1 and 2 equally\n",
"pd=p/2 \n",
"\n",
"#Power delivered\n",
"#Port 1\n",
"p1=pd*(1-T1**2) \n",
"print 'Power at port 1: %0.3f'%p1,'mW'\n",
"\n",
"#Port2\n",
"p2=pd*(1-T2**2) \n",
"print 'Power at port 2: %0.3f'%p2,'mW'"
],
"language": "python",
"metadata": {},
"outputs": [
{
"output_type": "stream",
"stream": "stdout",
"text": [
"Power at port 1: 19.753 mW\n",
"Power at port 2: 19.835 mW\n"
]
}
],
"prompt_number": 48
},
{
"cell_type": "heading",
"level": 2,
"metadata": {},
"source": [
"Page Number: 216 Example 4.15"
]
},
{
"cell_type": "code",
"collapsed": false,
"input": [
"from __future__ import division \n",
"#Given\n",
"T1=1/2 \n",
"T2=3/5 \n",
"T3=0 \n",
"T4=4/5 \n",
"p=500e-3 #W\n",
"#S matrix for magic Tee\n",
"x=1/sqrt(2) \n",
"s=[0, 0, x, x, 0, 0, x, -x, x ,x, 0, 0, x, -x, 0, 0] \n",
"#Using the input output relation\n",
"#[b]=[s]*[a]\n",
"b=[0.6565, 0.7576, 0.5536, 0.0892] \n",
"\n",
"#(i) Power transmitted through ports\n",
"#Port 1\n",
"p1=(1/2)*b[0]**2*(1-T1**2) \n",
"print 'Power at port 1: %0.3f'%p1,'W'\n",
"\n",
"#Port2\n",
"p2=(1/2)*(b[1]**2)*(1-(T2**2)) \n",
"print 'Power at port 2: %0.3f'%p2, 'W'\n",
"\n",
"#Port 4\n",
"p4=(1/2)*b[3]**2*(1-T4**2) \n",
"print 'Power at port 4: %0.3f'%p4,'W'\n",
"\n",
"#(ii) Power reflected at port 3\n",
"#Port 3\n",
"p3=p*b[2]**2 \n",
"print 'Power at port 3: %0.3f'%p3,'W'\n",
"\n",
"#(iii) Power absorbed\n",
"pabs=p-(p1+p2+p3+p4) \n",
"print 'Power absorbed: %0.3f'%pabs,'W'\n",
"\n",
"#Answer for power absorbed is calculated wrong in book"
],
"language": "python",
"metadata": {},
"outputs": [
{
"output_type": "stream",
"stream": "stdout",
"text": [
"Power at port 1: 0.162 W\n",
"Power at port 2: 0.184 W\n",
"Power at port 4: 0.001 W\n",
"Power at port 3: 0.153 W\n",
"Power absorbed: 0.000 W\n"
]
}
],
"prompt_number": 53
},
{
"cell_type": "heading",
"level": 2,
"metadata": {},
"source": [
"Page Number: 236 Example 4.18"
]
},
{
"cell_type": "code",
"collapsed": false,
"input": [
"from numpy import mat, set_printoptions\n",
"#Given\n",
"C=10 #dB\n",
"D=30 #dB\n",
"\n",
"#Parameters\n",
"bet=10**(-C/20) \n",
"x=bet/(10**(D/20)) \n",
"a=sqrt(1-(bet*bet)) \n",
"#Scattering matrix\n",
"#Assuming symmetery\n",
"s=mat([[0, a ,x ,(bet*1J)],[ a, 0, (bet*1J), x],[ x ,(bet*1J), 0 ,a],[ (bet*1J) ,x ,a, 0] ])\n",
"set_printoptions(precision=3)\n",
"print 'Scattering matrix:\\n',s "
],
"language": "python",
"metadata": {},
"outputs": [
{
"output_type": "stream",
"stream": "stdout",
"text": [
"Scattering matrix:\n",
"[[ 0.000+0.j 0.949+0.j 0.010+0.j 0.000+0.316j]\n",
" [ 0.949+0.j 0.000+0.j 0.000+0.316j 0.010+0.j ]\n",
" [ 0.010+0.j 0.000+0.316j 0.000+0.j 0.949+0.j ]\n",
" [ 0.000+0.316j 0.010+0.j 0.949+0.j 0.000+0.j ]]\n"
]
}
],
"prompt_number": 57
},
{
"cell_type": "heading",
"level": 2,
"metadata": {},
"source": [
"Page Number: 238 Example 4.20"
]
},
{
"cell_type": "code",
"collapsed": false,
"input": [
" \n",
"#Given\n",
"vswr=2 \n",
"D1=8 #mW\n",
"D2=2 #mW\n",
"\n",
"#Reflection coefficient at arm 4\n",
"T=(vswr-1)/(vswr+1) \n",
"#Powwe delivered to D1\n",
"P=(D1*100)/(1-T**2) \n",
"P1=0.99*P \n",
"#Power reflected at D1\n",
"W1=(P/100)*T*T \n",
"#Power reflected at load\n",
"W2=D2-W1 \n",
"Tt=sqrt((W2*100)/(P1)) \n",
"pt=(1+Tt)/(1-Tt) \n",
"print 'VSWR:%0.3f'%pt\n",
"Pl=P1*(1-(Tt*Tt)) \n",
"print 'Power delivered:' ,Pl,'mW'\n",
"\n",
"#Answer for P1 should be 792 but it is given as 800"
],
"language": "python",
"metadata": {},
"outputs": [
{
"output_type": "stream",
"stream": "stdout",
"text": [
"VSWR:2.008\n",
"Power delivered: 791.0 mW\n"
]
}
],
"prompt_number": 59
},
{
"cell_type": "heading",
"level": 2,
"metadata": {},
"source": [
"Page Number: 239 Example 4.21"
]
},
{
"cell_type": "code",
"collapsed": false,
"input": [
"from numpy import mat, set_printoptions\n",
"#Given\n",
"I=30 #dB\n",
"Il=0.4 #dB\n",
"\n",
"S12=10**(I/-20) \n",
"S21=10**(Il/-20) \n",
"s=mat([[0, S12],[ S21, 0] ])\n",
"set_printoptions(precision=3)\n",
"print 'Scattering matrix:\\n' ,s"
],
"language": "python",
"metadata": {},
"outputs": [
{
"output_type": "stream",
"stream": "stdout",
"text": [
"Scattering matrix:\n",
"[[ 0. 0.032]\n",
" [ 0.955 0. ]]\n"
]
}
],
"prompt_number": 60
},
{
"cell_type": "heading",
"level": 2,
"metadata": {},
"source": [
"Page Number: 240 Example 4.22"
]
},
{
"cell_type": "code",
"collapsed": false,
"input": [
"from numpy import mat, set_printoptions\n",
"#Given\n",
"I=30 #dB\n",
"Il=2 #dB\n",
"p=1.3 \n",
"\n",
"#Elelments\n",
"T=(p-1)/(p+1) \n",
"S11=T \n",
"S22=T \n",
"S33=T \n",
"S12=10**(-Il/20) \n",
"S13=10**(-I/20) \n",
"S21=S13 \n",
"S32=S13 \n",
"S23=S12 \n",
"S31=S23 \n",
"s=mat([[S11, S21, S31] ,[S12, S22, S32], [S13, S23 ,S33] ])\n",
"set_printoptions(precision=3)\n",
"print 'Scattering matrix:\\n' ,s"
],
"language": "python",
"metadata": {},
"outputs": [
{
"output_type": "stream",
"stream": "stdout",
"text": [
"Scattering matrix:\n",
"[[ 0.13 0.032 0.794]\n",
" [ 0.794 0.13 0.032]\n",
" [ 0.032 0.794 0.13 ]]\n"
]
}
],
"prompt_number": 61
},
{
"cell_type": "heading",
"level": 2,
"metadata": {},
"source": [
"Page Number: 249 Example 4.23"
]
},
{
"cell_type": "code",
"collapsed": false,
"input": [
" \n",
"#Given\n",
"f=10e9 #Hz\n",
"u=4e-7*pi \n",
"c=3e8 #m/s\n",
"a=2.29 #cm\n",
"a1=a/100 \n",
"b=1.02 #cm\n",
"b1=b/100 \n",
"\n",
"#E/H\n",
"w=2*pi*f \n",
"EbyH=(w*u)/sqrt(((w/c)**2)+((pi/a1)**2)) \n",
"lam=c/f \n",
"lamc=2*a1 \n",
"d=(1/4)*(lam/sqrt(1-((lam/lamc)**2))) \n",
"print 'Position: %0.3f'%(d*100), 'cm'\n",
"\n",
"#Answer for positon is calculated wrong in book"
],
"language": "python",
"metadata": {},
"outputs": [
{
"output_type": "stream",
"stream": "stdout",
"text": [
"Position: 0.993 cm\n"
]
}
],
"prompt_number": 62
},
{
"cell_type": "heading",
"level": 2,
"metadata": {},
"source": [
"Page Number: 250 Example 4.24"
]
},
{
"cell_type": "code",
"collapsed": false,
"input": [
"from numpy import mat, set_printoptions\n",
"#Given\n",
"#As it is perfectly matched\n",
"S12=1/sqrt(2) \n",
"S21=S12 \n",
"s=mat([[0 ,S12] ,[S21, 0] ])\n",
"set_printoptions(precision=3)\n",
"print 'Scattering matrix:\\n' , s"
],
"language": "python",
"metadata": {},
"outputs": [
{
"output_type": "stream",
"stream": "stdout",
"text": [
"Scattering matrix:\n",
"[[ 0. 0.707]\n",
" [ 0.707 0. ]]\n"
]
}
],
"prompt_number": 65
}
],
"metadata": {}
}
]
}