diff options
Diffstat (limited to '25/CH2')
| -rwxr-xr-x | 25/CH2/EX2.1/2_1.sce | 16 | ||||
| -rwxr-xr-x | 25/CH2/EX2.2/2_2.sce | 18 | ||||
| -rwxr-xr-x | 25/CH2/EX2.3/2_3.sce | 23 | ||||
| -rwxr-xr-x | 25/CH2/EX2.4/2_4.sce | 21 | ||||
| -rwxr-xr-x | 25/CH2/EX2.5/2_5.sce | 10 | ||||
| -rwxr-xr-x | 25/CH2/EX2.6/2_6.sce | 18 | ||||
| -rwxr-xr-x | 25/CH2/EX2.7/2_7.sce | 27 | ||||
| -rwxr-xr-x | 25/CH2/EX2.7/ex2_7.jpg | bin | 0 -> 24265 bytes |
8 files changed, 133 insertions, 0 deletions
diff --git a/25/CH2/EX2.1/2_1.sce b/25/CH2/EX2.1/2_1.sce new file mode 100755 index 000000000..84802b310 --- /dev/null +++ b/25/CH2/EX2.1/2_1.sce @@ -0,0 +1,16 @@ +// example 2.1,page no.-24.
+// program to calculate wavelength,phase velocity and wave impedence.
+f=3*10^9;
+mur=3;
+muo=4*%pi*10^-7;
+eipsilao=8.854*10^-12;
+eipsilar=7;
+mue=muo*mur;
+eipsila=eipsilao*eipsilar;
+Vp=sqrt(1/(mue*eipsila));
+lamda=Vp/f;
+eta=sqrt(mue/eipsila);
+//Result
+disp(Vp,'phase velocity in meter per second=') // phase velocity.
+disp(lamda,'wavelength in meter=') // wavelength.
+disp(eta,'wave impedence in ohm=') // wave impedence.
\ No newline at end of file diff --git a/25/CH2/EX2.2/2_2.sce b/25/CH2/EX2.2/2_2.sce new file mode 100755 index 000000000..14dd788f5 --- /dev/null +++ b/25/CH2/EX2.2/2_2.sce @@ -0,0 +1,18 @@ +// example:-2.2.page no.-26.
+// progarm to find out skin depth of aluminium,copper,gold and silver at frequency 10GHZ.
+f=10*10^9;
+muo=4*%pi*10^-7; // permeability in free space.
+omega=2*%pi*f;
+sigma_aluminium=3.816*10^7;
+sigma_copper=5.813*10^7;
+sigma_gold=4.098*10^7;
+sigma_silver=6.173*10^7;
+delta1=sqrt(2/(omega*muo*sigma_aluminium));
+delta2=sqrt(2/(omega*muo*sigma_copper));
+delta3=sqrt(2/(omega*muo*sigma_gold));
+delta4=sqrt(2/(omega*muo*sigma_silver));
+//result
+disp(delta1,'skin depth of aluminium in meter=') // skin depth of aluminium.
+disp(delta2,'skin depth of copper in meter=') //skin depth of copper.
+disp(delta3,'skin depth of gold in meter=') //skin depth of gold.
+disp(delta4,'skin depth of silver in meter=') //skin depth of silver.
\ No newline at end of file diff --git a/25/CH2/EX2.3/2_3.sce b/25/CH2/EX2.3/2_3.sce new file mode 100755 index 000000000..a949b68fb --- /dev/null +++ b/25/CH2/EX2.3/2_3.sce @@ -0,0 +1,23 @@ +// example:-2.3,page no.-31.
+//program to find the resulting fields by assumibg plane waves on either side of the current sheet and enforcing the boundary conditions.
+syms E x E1 E2 H1 H2 z Jo A B c N n d ko y;
+sym('n*(E2-E1)=0'); //boundary condition to be satisfied at z=0
+sym('z*(E2-E1)=0'); // " " " "
+sym('n*(H2-H1)=Jo'); // " " " "
+sym('z*(H2-H1)=Jo'); // " " " "
+E1=A*N*exp(%i*ko*z)*x; // x component of electric field (region z<0).
+H1=A*N*exp(%i*ko*z)*(-y); // -y component of magnetic field (region z<0).
+E2=B*N*exp(-%i*ko*z)*x; // x component of electric field (region z>0).
+H2=B*N*exp(-%i*ko*z)*y; // y component of electric field (region z>0).
+disp(E1,'for z<0, E1=')
+disp(H1,'for z<0, H1=')
+disp(E2,'for z>0, E2=')
+disp(H2,'for z>0, H2=')
+//from boundary conditions imposed.we get:-
+c=[-1 -1;1 -1];
+d=[A;B];
+c*d==[Jo;0];
+d=inv(c)*[Jo;0];
+//result
+// A=-Jo/2; B=-Jo/2.
+disp(d)
\ No newline at end of file diff --git a/25/CH2/EX2.4/2_4.sce b/25/CH2/EX2.4/2_4.sce new file mode 100755 index 000000000..eba66e05c --- /dev/null +++ b/25/CH2/EX2.4/2_4.sce @@ -0,0 +1,21 @@ +// example:-2.4,page no.-34.
+// program to show that a circularly polarized plane wave can be decomposed in to RHCP and LHCP.
+A=sym('A');
+B=sym('B');
+Eo=sym('Eo');
+x=sym('x');
+y=sym('y');
+Ko=sym('Ko');
+z=sym('z');
+E=Eo*(x+2*y)*exp(-%i*Ko*z); // given
+// can be written as:=>E=A*(x-y)*exp(-%i*Ko*z)+B*(x+y)*exp(-%i*Ko*z),so
+p=[1 1;-%i/2 %i/2];
+q=[A;B];
+r=[1;1];
+p*q==Eo*r;
+q=inv(p)*Eo*r;
+//result
+disp('value of A and B will be=')
+disp(q)
+disp(q(1,1)*(x-y)*exp(-%i*Ko*z)+q(2,1)*(x+y)*exp(-%i*Ko*z),'E=')
+//conclusion:-any linearly polarized wave can be decomposed in to two circularly polarized waves.
\ No newline at end of file diff --git a/25/CH2/EX2.5/2_5.sce b/25/CH2/EX2.5/2_5.sce new file mode 100755 index 000000000..ddef6cb5b --- /dev/null +++ b/25/CH2/EX2.5/2_5.sce @@ -0,0 +1,10 @@ +// example:-2.5,page no.-36.
+// program to compute the poynting vector for the plane wave field.
+syms E Eo H k s n N x r;
+E=Eo*exp(-%i*k*r); // electric field.
+H=(E/N)*n; //N is intrinsic impedence,n is unit vector.
+H1=conj(H) // conjugate of magnetic field.
+s=E*H1;
+//result
+disp(s,'poynting vector is(meter square)=')
+disp('which shows that power density is flowing in the direction of propagation.')
\ No newline at end of file diff --git a/25/CH2/EX2.6/2_6.sce b/25/CH2/EX2.6/2_6.sce new file mode 100755 index 000000000..97d875848 --- /dev/null +++ b/25/CH2/EX2.6/2_6.sce @@ -0,0 +1,18 @@ +// example:-2.6,page no.-46.
+// program to compute propagation constan,impedence,skin depth,reflection and transmission coefficient.
+f=1*10^9;
+omega=2*%pi*f;
+sigma=5.813*10^7; // for copper.
+mue=4*%pi*10^-7; // permeability in free space.
+delta=sqrt(2/(mue*sigma*omega)); // skin depth.
+gama=((1+%i)/delta); //propagation constant.
+eta=gama/sigma; // impedence
+etao=377; //intrinsic impedence in free space.
+tao=((eta-etao)/(eta+etao)); // reflection coefficient.
+t=(2*eta)/(eta+etao); //transmission coefficient.
+// result
+disp(delta,'skin depth in meter=')
+disp(gama,'propagation constant=')
+disp(eta,'intrinsic impedence in ohm=')
+disp(tao,'reflection coefficient=')
+disp(t,'transmission coefficient=')
\ No newline at end of file diff --git a/25/CH2/EX2.7/2_7.sce b/25/CH2/EX2.7/2_7.sce new file mode 100755 index 000000000..3e358c8e6 --- /dev/null +++ b/25/CH2/EX2.7/2_7.sce @@ -0,0 +1,27 @@ +// example:-2.7.page no.-50.
+// program to plot the reflection coefficients for parallel and perpendicular polarized plane waves incident from free space on to a dielectric region with Er=2.55,versus incidence angle.
+Er=2.55; // relaitve permittivity of dielectric medium.
+N1=377; // intrinsic impedence
+N2=N1/sqrt(Er); // intrinsic impedence of dielectric medium.
+xb=asin(sqrt(1/(1+1/2.55))); // brewster angle valid only in case of parallel polarization.
+xt=acos(sqrt(1-(1/Er)^2*sin(xb))); // angle of transmission.
+xi=[0:0.01:%pi/2]; // incidence angle.
+// for parallel polarization
+N2=N2*cos(xt);
+N1=N1*cos(xi);
+Tpar=(N2-N1)./(N2+N1);
+w=abs(Tpar);
+// result
+subplot(1,2,1)
+xtitle("parallel polarization","xi(incidence angle)","Tpar(reflection coefficient)")
+plot2d(xi,w,style=3,rect=[0,0,%pi/2,1])
+// for perpendicular polarization. //NOTE:- in case of this polarization.there is no brewster angle.
+xt=acos(sqrt(1-(1/Er)^2*sin(xi)));
+n1=377.*cos(xt);
+n2=(377/sqrt(Er)).*cos(xi);
+Tper=(n2-n1)./(n1+n2);
+z=abs(Tper);
+//result
+subplot(1,2,2)
+xtitle("perpendicular polarization","xi(inxidence angle)","Tper(reflection coefficient)")
+plot2d(xi,z,style=2,rect=[0,0,%pi/2,1])
\ No newline at end of file diff --git a/25/CH2/EX2.7/ex2_7.jpg b/25/CH2/EX2.7/ex2_7.jpg Binary files differnew file mode 100755 index 000000000..f7a8c0272 --- /dev/null +++ b/25/CH2/EX2.7/ex2_7.jpg |
