summaryrefslogtreecommitdiff
path: root/3773/CH23/EX23.5/Ex23_5.sce
diff options
context:
space:
mode:
authorprashantsinalkar2017-10-10 12:27:19 +0530
committerprashantsinalkar2017-10-10 12:27:19 +0530
commit7f60ea012dd2524dae921a2a35adbf7ef21f2bb6 (patch)
treedbb9e3ddb5fc829e7c5c7e6be99b2c4ba356132c /3773/CH23/EX23.5/Ex23_5.sce
parentb1f5c3f8d6671b4331cef1dcebdf63b7a43a3a2b (diff)
downloadScilab-TBC-Uploads-7f60ea012dd2524dae921a2a35adbf7ef21f2bb6.tar.gz
Scilab-TBC-Uploads-7f60ea012dd2524dae921a2a35adbf7ef21f2bb6.tar.bz2
Scilab-TBC-Uploads-7f60ea012dd2524dae921a2a35adbf7ef21f2bb6.zip
initial commit / add all books
Diffstat (limited to '3773/CH23/EX23.5/Ex23_5.sce')
-rw-r--r--3773/CH23/EX23.5/Ex23_5.sce29
1 files changed, 29 insertions, 0 deletions
diff --git a/3773/CH23/EX23.5/Ex23_5.sce b/3773/CH23/EX23.5/Ex23_5.sce
new file mode 100644
index 000000000..5ca7f473e
--- /dev/null
+++ b/3773/CH23/EX23.5/Ex23_5.sce
@@ -0,0 +1,29 @@
+//Chapter 23: Ground Wave Propagation
+//Example 23-5.1
+clc;
+
+//Variable Initialization
+f1 = 5e3 //Frequency (Hz)
+f2 = 50e3 //Frequency (Hz)
+f3 = 500e3 //Frequency (Hz)
+sigma = 5e-5 //Standard deviation of surface irregularities (unitless)
+eps_r = 15.0 //Relative permittivity (unitless)
+mu = %pi*4e-7 //Absolute Permeability (H/m)
+
+//Calculations
+w1 = 2*%pi*f1 //Angular frequency (rad/s)
+w2 = 2*%pi*f2 //Angular frequency (rad/s)
+w3 = 2*%pi*f3 //Angular frequency (rad/s)
+
+
+Zs1 = sqrt((w1*mu)/sqrt(sigma**2 + (w1**2)*eps_r)) //Surface impedence for f1 (ohm)
+Zs2 = sqrt((w2*mu)/sqrt(sigma**2 + (w2**2)*eps_r)) //Surface impedence for f2 (ohm)
+Zs3 = sqrt((w3*mu)/sqrt(sigma**2 + (w3**2)*eps_r)) //Surface impedence for f3 (ohm)
+
+//Result
+mprintf( "The surface impedence for 5kHz is %.5f ohms",Zs1)
+mprintf( "\nThe surface impedence for 50kHz is %.5f ohms", Zs2)
+mprintf( "\nThe surface impedence for 500kHz is %.5f ohms", Zs3)
+
+//An error has been made in calculation/substitution of square root of
+//(sigma**2 + (w1**2)*eps_r) and in the second case, the mistake in the calculation of (w2*mu)/sqrt(sigma**2 + (w2**2)*eps_r)