summaryrefslogtreecommitdiff
path: root/3773/CH12/EX12.2/Ex12_2.sce
diff options
context:
space:
mode:
authorprashantsinalkar2017-10-10 12:27:19 +0530
committerprashantsinalkar2017-10-10 12:27:19 +0530
commit7f60ea012dd2524dae921a2a35adbf7ef21f2bb6 (patch)
treedbb9e3ddb5fc829e7c5c7e6be99b2c4ba356132c /3773/CH12/EX12.2/Ex12_2.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/CH12/EX12.2/Ex12_2.sce')
-rw-r--r--3773/CH12/EX12.2/Ex12_2.sce32
1 files changed, 32 insertions, 0 deletions
diff --git a/3773/CH12/EX12.2/Ex12_2.sce b/3773/CH12/EX12.2/Ex12_2.sce
new file mode 100644
index 000000000..ab3113674
--- /dev/null
+++ b/3773/CH12/EX12.2/Ex12_2.sce
@@ -0,0 +1,32 @@
+//Chapter 12: The Cylindrical Antenna and the Moment Method
+//Example 12-12.2
+clc;
+
+//Variable Initialization
+z_load = 2.083 + 1605*%i //Conjugate matched load (ohm)
+e0 = 1.0 //Electric field magnitude (unitless)
+l = 1/10.0 //Length of dipole (lambda)
+ima = 0+1*%i //Imaginary number
+
+z11_exact = 0.4935 - 3454*%i //Exact impedance vector(ohm)
+z11_apprx = 0.4944 - 3426*%i //Approximate impedance vector(ohm)
+z12_exact = 0.4935 + 1753*%i //Exact impedance vector(ohm)
+z12_apprx = 0.4945 + 1576*%i //Approximate impedance vector(ohm)
+z13_exact = 0.4935 + 129.9*%i //Exact impedance vector(ohm)
+z13_apprx = 0.4885 + 132.2*%i //Approximate impedance vector(ohm)
+
+//Calculation
+d = l/4 //Length of each segment (lambda)
+vm = (2*e0/(2*%pi))*tan(2*%pi*d/2) //Voltage vector (VA)
+z22 = z11_exact + z_load //Impedance matrix for loaded dipole (VA)
+zmat_exact =([z11_exact+z13_exact,z12_exact;2*z12_exact,z22])//Z(impedance) matrix (unitless)
+vmat = ([vm;vm]) //Voltage matrix (unitless)
+[i1]= linsolve(zmat_exact,vmat) //Current matrix (unitless)
+i1=i1*-1
+i3 = i1(1) //Current vector (unitless)
+e_zn = (60*tan(2*%pi*d/2))*ima //Free space electric field (V/m)
+e_s = i1(1)*e_zn + i1(2)*e_zn + i3*e_zn //Scattered field (V/m)
+sigma = 4*%pi*(abs(e_s)**2)/(abs(e0)**2) //Radar Cross section (lambda**2)
+
+//Result
+mprintf("The radar cross section using exact values of Z matrix is %.4f lambda square",sigma(1))