initial commit / add all books

7 files changed, 114 insertions, 0 deletions

diff --git a/3773/CH5/EX5.1/Ex5_1.sce b/3773/CH5/EX5.1/Ex5_1.sce
new file mode 100644
index 000000000..e0ea98851
--- /dev/null
+++ b/3773/CH5/EX5.1/Ex5_1.sce
@@ -0,0 +1,17 @@
+//Chapter 5: Point Source and Their Arrays
+//Example 5-6.1
+clc;
+clear;
+
+//Variable Initialization
+Um=1 //Maximum radiation intensity (unitless)
+
+//Calculation
+deff('z=f(x,y)','z=cos(x)*sin(x)')//Integrand(Unitless)
+X=[0 0;%pi/2 %pi/2;%pi/2 0];
+Y=[0 0;2*%pi 2*%pi;0 2*%pi];
+[I,err]=int2d(X,Y,f)//Total power radiated (relative to Um)
+
+D=(4*%pi)/I //Directivity (unitless)
+
+mprintf('The directivity is %f',D)
diff --git a/3773/CH5/EX5.2/Ex5_2.sce b/3773/CH5/EX5.2/Ex5_2.sce
new file mode 100644
index 000000000..44fea4da9
--- /dev/null
+++ b/3773/CH5/EX5.2/Ex5_2.sce
@@ -0,0 +1,15 @@
+//Chapter 5: Point Source and Their Arrays
+//Example 5-6.2
+clc;
+clear;
+
+//Variable Initialization
+Um=1 //Maximum radiation intensity (unitless)
+deff('z=f(x,y)','z=cos(x)*sin(x)')//Integrand(Unitless)
+X=[0 0;%pi/2 %pi/2;%pi/2 0];
+Y=[0 0;2*%pi 2*%pi;0 2*%pi];
+[I,err]=int2d(X,Y,f)//Total power radiated (relative to Um)
+
+D=(4*%pi)/(2*I) //Directivity (unitless)
+
+mprintf('The directivity is %f',D)
diff --git a/3773/CH5/EX5.3/Ex5_3.sce b/3773/CH5/EX5.3/Ex5_3.sce
new file mode 100644
index 000000000..f63f8b2e2
--- /dev/null
+++ b/3773/CH5/EX5.3/Ex5_3.sce
@@ -0,0 +1,15 @@
+//Chapter 5: Point Source and Their Arrays
+//Example 5-6.3
+clc;
+clear;
+
+//Variable Initialization
+Um=1 //Maximum radiation intensity (unitless)
+deff('z=f(x,y)','z=sin(x)**2')//Integrand(Unitless)
+X=[0 0;%pi %pi;%pi 0];
+Y=[0 0;2*%pi 2*%pi;0 2*%pi];
+[I,err]=int2d(X,Y,f)//Total power radiated (relative to Um)
+
+D=(4*%pi)/I //Directivity (unitless)
+
+mprintf('The directivity is %.3f',D)
diff --git a/3773/CH5/EX5.4/Ex5_4.sce b/3773/CH5/EX5.4/Ex5_4.sce
new file mode 100644
index 000000000..14fe3468c
--- /dev/null
+++ b/3773/CH5/EX5.4/Ex5_4.sce
@@ -0,0 +1,15 @@
+//Chapter 5: Point Source and Their Arrays
+//Example 5-6.4
+clc;
+clear;
+
+//Variable Initialization
+Um=1 //Maximum radiation intensity (unitless)
+deff('z=f(x,y)','z=sin(x)**3')//Integrand(Unitless)
+X=[0 0;%pi %pi;%pi 0];
+Y=[0 0;2*%pi 2*%pi;0 2*%pi];
+[I,err]=int2d(X,Y,f)//Total power radiated (relative to Um)
+
+D=(4*%pi)/I //Directivity (unitless)
+
+mprintf('The directivity is %.1f',D)
diff --git a/3773/CH5/EX5.5/Ex5_5.sce b/3773/CH5/EX5.5/Ex5_5.sce
new file mode 100644
index 000000000..b07b90df0
--- /dev/null
+++ b/3773/CH5/EX5.5/Ex5_5.sce
@@ -0,0 +1,15 @@
+//Chapter 5: Point Source and Their Arrays
+//Example 5-6.5
+clc;
+clear;
+
+//Variable Initialization
+Um=1 //Maximum radiation intensity (unitless)
+deff('z=f(x,y)','z=sin(x)*cos(x)**2')//Integrand(Unitless)
+X=[0 0;%pi/2 %pi/2;%pi/2 0];
+Y=[0 0;2*%pi 2*%pi;0 2*%pi];
+[I,err]=int2d(X,Y,f)//Total power radiated (relative to Um)
+
+D=(4*%pi)/I //Directivity (unitless)
+
+mprintf('The directivity is %.1f',D)
diff --git a/3773/CH5/EX5.6/Ex5_6.sce b/3773/CH5/EX5.6/Ex5_6.sce
new file mode 100644
index 000000000..7a6a493fc
--- /dev/null
+++ b/3773/CH5/EX5.6/Ex5_6.sce
@@ -0,0 +1,21 @@
+//Chapter 5: Point Source and Their Arrays
+//Example 5-6.6
+clc;
+
+//Variable Initialization
+lobes = [0.25,0.37,0.46,0.12,0.07]  //Normalized power of lobes (unitless)
+
+//Calculation
+ohm_a = 0                   //Beam area (sr)
+sum_lobes = 0               //Sum of all lobes (unitless)
+for i=lobes
+    ohm_a =ohm_a + 2*%pi*(%pi/36)*(i)
+    sum_lobes =sum_lobes + i
+end
+D = 4*%pi/ohm_a         //Directivity (unitless)
+D_db = 10*log10(D)     //Directivity (in dBi)
+e_m = lobes(1)/sum_lobes    //Beam efficiency (unitless)
+
+//Result
+mprintf("The directivity is %d or %.1f dBi",D,D_db)
+mprintf("\nThe beam efficiency is %.1f",e_m)
diff --git a/3773/CH5/EX5.7/Ex5_7.sce b/3773/CH5/EX5.7/Ex5_7.sce
new file mode 100644
index 000000000..3648e722f
--- /dev/null
+++ b/3773/CH5/EX5.7/Ex5_7.sce
@@ -0,0 +1,16 @@
+//Chapter 5: Point Source and Their Arrays
+//Example 5-21.1
+clc;
+
+//Variable Initialization
+a = 25              //Height of vertical conducting wall (m)
+r = 100             //Distance to the receiver (m)
+wave_lt = 10e-2     //Transmitter dimension (m)
+
+//Calculation
+k = sqrt(2/(r*wave_lt))    //Constant (unitless)
+S_av = (r*wave_lt)/(4*(%pi**2)*(a**2))  //Relative signal level (unitless)
+S_av_db = 10*log10(S_av)   //Signal level (in db)
+
+//Result
+mprintf("The signal level at the receiver is %.5f or %.0f dB",S_av,S_av_db)