summaryrefslogtreecommitdiff
path: root/25/CH4
diff options
context:
space:
mode:
authorpriyanka2015-06-24 15:03:17 +0530
committerpriyanka2015-06-24 15:03:17 +0530
commitb1f5c3f8d6671b4331cef1dcebdf63b7a43a3a2b (patch)
treeab291cffc65280e58ac82470ba63fbcca7805165 /25/CH4
downloadScilab-TBC-Uploads-b1f5c3f8d6671b4331cef1dcebdf63b7a43a3a2b.tar.gz
Scilab-TBC-Uploads-b1f5c3f8d6671b4331cef1dcebdf63b7a43a3a2b.tar.bz2
Scilab-TBC-Uploads-b1f5c3f8d6671b4331cef1dcebdf63b7a43a3a2b.zip
initial commit / add all books
Diffstat (limited to '25/CH4')
-rwxr-xr-x25/CH4/EX4.1/4_1.sce26
-rwxr-xr-x25/CH4/EX4.2/4_2.sce24
-rwxr-xr-x25/CH4/EX4.3/4_3.sce9
-rwxr-xr-x25/CH4/EX4.4/4_4.sce15
-rwxr-xr-x25/CH4/EX4.4/ex4_4.jpgbin0 -> 40661 bytes
-rwxr-xr-x25/CH4/EX4.5/4_5.sce22
-rwxr-xr-x25/CH4/EX4.7/4_7.sce14
-rwxr-xr-x25/CH4/EX4.9/4_9.sce15
8 files changed, 125 insertions, 0 deletions
diff --git a/25/CH4/EX4.1/4_1.sce b/25/CH4/EX4.1/4_1.sce
new file mode 100755
index 000000000..2552cc333
--- /dev/null
+++ b/25/CH4/EX4.1/4_1.sce
@@ -0,0 +1,26 @@
+// example:-4.1,page no.-148.
+// program to find the cut off frequency fo the first four propagating modes.
+a=0.02286;b=0.01016;f=10*10^9;k=209.44;sigma=5.8*10^7;mue=4*%pi*10^-7;
+c=3*10^8;
+m=0;n=1;
+fc=(c/(%pi*2))*sqrt(((%pi*m)/a)^2+((%pi*n)/b)^2);
+fc=fc/(10^9);
+disp(fc,'cut-off frequency for TE01 mode in GHZ=')
+m=1;n=0;
+fc=(c/(%pi*2))*sqrt(((%pi*m)/a)^2+((%pi*n)/b)^2);
+fc=fc/(10^9);
+disp(fc,'cut-off frequency for TE10 mode in GHZ=')
+m=2;n=0;
+fc=(c/(%pi*2))*sqrt(((%pi*m)/a)^2+((%pi*n)/b)^2);
+fc=fc/(10^9);
+disp(fc,'cut-off frequency for TE20 mode in GHZ=')
+m=1;n=1;
+fc=(c/(%pi*2))*sqrt(((%pi*m)/a)^2+((%pi*n)/b)^2);
+fc=fc/(10^9);
+disp(fc,'cut-off frequency for TE11 mode in GHZ=')
+B=sqrt(k^2-(%pi/a)^2) // for TE10 mode
+Rs=sqrt(((2*%pi*f)*mue)/(2*sigma)); // surface resistance.
+disp(Rs,'surface resistance in ohm=')
+ac=(Rs/(a^3*b*B*k*377))*((2*b*%pi^2)+(a^3*k^2)) //attenuation constant.
+ac=-20*(-ac)*log10(%e);
+disp(ac,'attenuation constant in dB/m=') \ No newline at end of file
diff --git a/25/CH4/EX4.2/4_2.sce b/25/CH4/EX4.2/4_2.sce
new file mode 100755
index 000000000..678822a5f
--- /dev/null
+++ b/25/CH4/EX4.2/4_2.sce
@@ -0,0 +1,24 @@
+//example:-4.2,page no.-160.
+//program to find the cut off frequency of two propagating modes of a circular waveguide.
+a=0.005;eipsilar=2.25;f=13*10^9;c=3*10^8;d=0.001;sigma=6.17*01^7;muo=4*%pi*10^-7;
+m=1;n=1;
+p11=1.841;p01=2.405;
+fc=(p11*c)/(2*%pi*a*sqrt(eipsilar));
+kc=p11/a;
+fc=fc/(10^9);
+disp(fc,'cut-off frequency for TE11 mode in GHZ')
+m=0;n=1;
+fc=(p01*c)/(2*%pi*a*sqrt(eipsilar));
+fc=fc/(10^9);
+disp(fc,'cut-off frequency for TE01 mode in GHZ')
+// so,TE01 can't be propagating mode.only TE11 will be.
+k=(2*%pi*f*sqrt(eipsilar))/c;
+disp(k,'k in m-1=')
+B=sqrt(k^2-kc^2);
+disp(B,'propagation constant of TE11 mode')
+ac=(k^2*d)/(2*B);
+Rs=sqrt((2*%pi*f*muo)/(2*sigma)); // surface resistance.
+acm=(Rs/(a*k*377*B))*(kc^2+((k^2)/(p11^2-1)));
+a=ac+acm;
+a=-20*(-0.547*0.5)*log10(%e);
+disp(a,'total attenuation factor in dB=') \ No newline at end of file
diff --git a/25/CH4/EX4.3/4_3.sce b/25/CH4/EX4.3/4_3.sce
new file mode 100755
index 000000000..828e03cad
--- /dev/null
+++ b/25/CH4/EX4.3/4_3.sce
@@ -0,0 +1,9 @@
+//example:-4.3,page no.-167.
+//program to find out the highest usable frequency.
+a=0.000889;b=0.0029464;eipsilar=2.2;c=3*10^8;
+// here (b/a)=3.3,so for this kc*a=0.47
+kc=0.47/a;
+fc=(c*kc)/(2*%pi*sqrt(eipsilar))
+fc=fc/(10^9);
+fmax=0.95*fc;
+disp(fmax,'maximum usable frequency in GHZ=') \ No newline at end of file
diff --git a/25/CH4/EX4.4/4_4.sce b/25/CH4/EX4.4/4_4.sce
new file mode 100755
index 000000000..a809edfb3
--- /dev/null
+++ b/25/CH4/EX4.4/4_4.sce
@@ -0,0 +1,15 @@
+//example:-4.4,page no.-175.
+// program to calculate and plot the propagation constant of first three propagating surface wave mode.
+eipsilar=2.55;c=3*10^8; // x=d/lamdao;
+x=0.001:0.01:1.2;
+for n=0:1:4
+y=sqrt(eipsilar-((n^2)./(4.*(x^2)*(eipsilar-1)))); // y=beta/lamdao;
+plot2d(x,y,style=2,rect=[0,0,1.2,1.6])
+end
+x=0.001:0.01:1.2;
+for n=1:1:4
+ y=sqrt(eipsilar-((((2.*n)-1)^2)./(16.*(x^2)*(eipsilar-1))))
+ plot2d(x,y,style=6,rect=[0,0,1.2,1.6])
+end
+xtitle("plot of propagation constant of first 4 mode","d/lamdao","beta/Ko");
+legend("TE MODE") \ No newline at end of file
diff --git a/25/CH4/EX4.4/ex4_4.jpg b/25/CH4/EX4.4/ex4_4.jpg
new file mode 100755
index 000000000..3056a4567
--- /dev/null
+++ b/25/CH4/EX4.4/ex4_4.jpg
Binary files differ
diff --git a/25/CH4/EX4.5/4_5.sce b/25/CH4/EX4.5/4_5.sce
new file mode 100755
index 000000000..2836fb01b
--- /dev/null
+++ b/25/CH4/EX4.5/4_5.sce
@@ -0,0 +1,22 @@
+// example:-4.5,page no.-180.
+// program to find width of a copper strip line conductor.
+eipsilar=2.20;Zo=50;b=0.0032;d=0.001,f=10^10;t=0.00001;
+c=3*10^8;Rs=0.026;A=4.74;
+x=(30*%pi)/(sqrt(eipsilar)*Zo);
+x=x-0.441;
+w=b*x;
+if ((sqrt(eipsilar)*Zo)<120)
+ disp("width of copper strip line conductor is 0.00266m")
+end
+K=(2*%pi*f*sqrt(eipsilar))/c;
+ad=(K*d)/2;
+ac=(2.7*(10^-3)*Rs*eipsilar*Zo*A)/(30*%pi*(b-t));
+a=ac+ad;
+a=20*a*log10(%e);
+lamda=c/(sqrt(eipsilar)*f);
+alamda=lamda*a;
+disp(K,'wave number=')
+disp(ad,'dielectric aattenuation=')
+disp(ac,'conductor attenuation=')
+disp(a,'total attenuation constant=')
+disp(alamda,'attenuation in dB/lamda=') \ No newline at end of file
diff --git a/25/CH4/EX4.7/4_7.sce b/25/CH4/EX4.7/4_7.sce
new file mode 100755
index 000000000..c08ce8f7a
--- /dev/null
+++ b/25/CH4/EX4.7/4_7.sce
@@ -0,0 +1,14 @@
+//example:-4.7,page no.-187.
+//program to calculate the width and length of microstrip line.
+eipsilae=1.87;//effective dielectric constant.
+Zo=50;q=%pi/2;c=3*10^8;
+f=2.5*10^9;
+ko=(2*%pi*f)/c;
+d=0.00127;
+eipsilar=2.20;
+// for w/d>2;
+B=7.985;
+w=3.081*d*100;
+disp(w,'width in centi meter=')
+l=(q*100)/(sqrt(eipsilae)*ko);
+disp(l,'length of microstrip in centi meter=') \ No newline at end of file
diff --git a/25/CH4/EX4.9/4_9.sce b/25/CH4/EX4.9/4_9.sce
new file mode 100755
index 000000000..944ed182f
--- /dev/null
+++ b/25/CH4/EX4.9/4_9.sce
@@ -0,0 +1,15 @@
+//example:-4.9,page no.-197.
+//program to calculate the group velocity.
+syms w c v;
+B=sym('B');
+ko=sym('ko');
+kc=sym('kc');
+ko=w/c;
+B=sqrt(ko^2-kc^2);
+v=diff(B,w);
+vg=v^(-1);
+vg=(c*B)/ko;
+vp=w/B;
+disp(vg,'group velocity=')
+disp(vp,'phase velocity=')
+disp('conclusion:-since B<ko,we have that vg<c<vp,which indicates that the phase velocity of a waveguide mode may be greater than the speed of light.but the group velocity will be lesser than the speed of light.') \ No newline at end of file