6 files changed, 109 insertions, 0 deletions
diff --git a/3773/CH4/EX4.1/Ex4_1.sce b/3773/CH4/EX4.1/Ex4_1.sce
new file mode 100644
index 000000000..aaf28c96d
--- /dev/null
+++ b/3773/CH4/EX4.1/Ex4_1.sce
@@ -0,0 +1,17 @@
+//Chapter 4: Radiation
+//Example 4-4.1
+clc;
+
+//Variable Initialization
+theta = 30              //Angle of radiation (degrees)
+epsilon_0 = 8.854e-12   //Permittivity of free space (F/m)
+I_dl = 10               //Current in length dl (A-m)
+r = 100e3               //Distance of point from origin (m)
+
+//Calculation
+E_mag = (I_dl*sin(theta*%pi/180))/(4*%pi*epsilon_0)                        //Magnitude of Electric field vector (V/m)
+H_mag = (I_dl*sin(theta*%pi/180))/(4)                        //Magnitude of Magnetic field vector (T)
+
+//Result
+disp(E_mag,"The magnitude of E vector in V/m ")
+mprintf("\nThe magnitude of H vector is %.3f /pi T", H_mag)
diff --git a/3773/CH4/EX4.2/Ex4_2.sce b/3773/CH4/EX4.2/Ex4_2.sce
new file mode 100644
index 000000000..e69146989
--- /dev/null
+++ b/3773/CH4/EX4.2/Ex4_2.sce
@@ -0,0 +1,14 @@
+//Chapter 4: Radiation
+//Example 4-4.2
+clc;
+
+//Variable Initialization
+v = 3e8         //Speed of light(m/s)
+f = 10e6        //Frequency (Hz)
+
+//Calculation
+w = 2*%pi*f     //Angular frequency(rad/s)
+r = v/w             //Distance (m)
+
+//Result
+mprintf("The distance for the specified condition is %.2f m",r)
diff --git a/3773/CH4/EX4.3/Ex4_3.sce b/3773/CH4/EX4.3/Ex4_3.sce
new file mode 100644
index 000000000..2cba15737
--- /dev/null
+++ b/3773/CH4/EX4.3/Ex4_3.sce
@@ -0,0 +1,15 @@
+//Chapter 4: Radiation
+//Example 4-4.3
+clc;
+
+//Variable Initialization
+c = 3e8             //Speed of light (m/s)
+f = 3e9             //Frequency (Hz)
+
+//Calculation
+v = 0.6*c           //60% of velocity of light (m/s)
+w = 2*%pi*f     //Angular frequency (rad/s)
+r = v/w             //Distance (m)
+
+//Result
+mprintf("The distance for the specified condition is %.6f m", r)
diff --git a/3773/CH4/EX4.4/Ex4_4.sce b/3773/CH4/EX4.4/Ex4_4.sce
new file mode 100644
index 000000000..e38f69daa
--- /dev/null
+++ b/3773/CH4/EX4.4/Ex4_4.sce
@@ -0,0 +1,18 @@
+//Chapter 4: Radiation
+//Example 4-5.1
+clc;
+
+//Variable Initialization
+dl = 1e-2       //Length of radiating element (m)
+I_eff = 0.5     //Effective current (A)
+f = 3e9         //Frequency (Hz)
+c = 3e8         //Velocity of light (m/s)
+
+//Calculation
+w = 2*%pi*f     //Angular Frequency (rad/s)
+P = 20*(w**2)*(I_eff**2)*(dl**2)/(c**2)     //Radiated power (W)
+
+//Result
+mprintf("The radiated power is %.2f W", P)
+
+//The answer obtained is varying compared with the textbook answer because of a calculation error
diff --git a/3773/CH4/EX4.5/Ex4_5.sce b/3773/CH4/EX4.5/Ex4_5.sce
new file mode 100644
index 000000000..fd3200079
--- /dev/null
+++ b/3773/CH4/EX4.5/Ex4_5.sce
@@ -0,0 +1,28 @@
+//Chapter 4: Radiation
+//Example 4-5.2
+clc;
+
+//Variable Initialization
+L = 5         //Length of radiating element (m)
+f1 = 30e3     //Frequency (Hz)  
+f2 = 30e6     //Frequency (Hz)  
+f3 = 15e6     //Frequency (Hz)
+c = 3e8       //Velocity of light (m/s)  
+
+//Calculation
+wave_lt1 = c/f1                 //Wavelength (m)
+wave_lt1 = wave_lt1 /10
+R_r1 = 800*(L/wave_lt1)**2      //Radiation resistance (ohm)
+
+wave_lt2 = c/f2                 //Wavelength (m)
+L = wave_lt2/2                  //Effective length (m)
+R_r2 = 200*(L/wave_lt2)**2      //Radiation resistance (ohm)
+
+wave_lt3 = c/f3                 //Wavelength (m)
+L = wave_lt3/4                  //Effective length (m)
+R_r3 = 400*(L/wave_lt3)**2      //Radiation resistance (ohm)
+
+//Result
+mprintf("The radiation resistance for f1 is %.2f ohms", R_r1)
+mprintf("\nThe radiation resistance for f2 is %d ohms",R_r2)
+mprintf("\nThe radiation resistance for f3 is %d ohms",R_r3)
diff --git a/3773/CH4/EX4.6/Ex4_6.sce b/3773/CH4/EX4.6/Ex4_6.sce
new file mode 100644
index 000000000..e018f715c
--- /dev/null
+++ b/3773/CH4/EX4.6/Ex4_6.sce
@@ -0,0 +1,17 @@
+//Chapter 4: Radiation
+//Example 4-6.1
+clc;
+
+//Variable Initialization
+Im = 5              //Maximum current (A)
+r = 1e3             //Distance (km)
+eta = 120*%pi   //Intrinsic impedance (ohm)
+theta = 60*%pi/180          //Angle of radiation (radians)
+
+//Calculation
+sin2 = sin(theta)**2       //Sine squared theta (unitless)
+P_av = (eta*(Im**2))/(8*(%pi**2)*(r**2))
+P_av = P_av*(cos(%pi/2*cos(theta))**2)/(sin2)                        //Average power (W)
+
+//Result
+mprintf("The average power available at 1km distance is %e W",P_av)