diff options
author | prashantsinalkar | 2017-10-10 12:27:19 +0530 |
---|---|---|
committer | prashantsinalkar | 2017-10-10 12:27:19 +0530 |
commit | 7f60ea012dd2524dae921a2a35adbf7ef21f2bb6 (patch) | |
tree | dbb9e3ddb5fc829e7c5c7e6be99b2c4ba356132c /3773/CH7/EX7.4/Ex7_4.sce | |
parent | b1f5c3f8d6671b4331cef1dcebdf63b7a43a3a2b (diff) | |
download | Scilab-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/CH7/EX7.4/Ex7_4.sce')
-rw-r--r-- | 3773/CH7/EX7.4/Ex7_4.sce | 37 |
1 files changed, 37 insertions, 0 deletions
diff --git a/3773/CH7/EX7.4/Ex7_4.sce b/3773/CH7/EX7.4/Ex7_4.sce new file mode 100644 index 000000000..ce861f911 --- /dev/null +++ b/3773/CH7/EX7.4/Ex7_4.sce @@ -0,0 +1,37 @@ +//Chapter 7: Loop, Slot and Horn Antennas +//Example 7-11.2 +clc; + +//Variable Initialization +n = 10 //Number of turns (unitless) +dia = 1e-3 //Diameter of copper wire (m) +dia_rod = 1e-2 //Diameter of ferrite rod (m) +len_rod = 10e-2 //Length of ferrite rod (m) +mu_r = 250 - 2.5*%i //Relative permeability (unitless) +mu_er = 50 //Effective relative permeability (unitless) +f = 1e6 //Frequency (Hz) +c = 3e8 //Speed of light (m/s) +mu_0 = %pi*4e-7 //Absolute permeability (H/m) + +//Calculations +wave_lt = c/f //Wavelength (m) +radius = dia_rod/2 +C_l = (2*%pi*radius)/(wave_lt) //Circumference of loop (m) +Rr = 197*(mu_er**2)*(n**2)*(C_l**4) //Radiation resistance (ohm) +Rf = 2*%pi*f*mu_er*(imag(mu_r)/real(mu_r))*mu_0*(n**2)*(%pi*radius**2)/len_rod //Loss resistance(ohm) +conduc = 1/((7e-5**2)*f*%pi*mu_er) //Conductivity (S/m) +delta = 1/(sqrt(f*%pi*mu_er*conduc)) //Depth of penetration(m) + +RL = n*(C_l/dia)*sqrt((f*mu_0)/(%pi*conduc)) //Ohmic resistance (ohm) +k = Rr/(RL+abs(Rf)) //Radiation efficiency (unitless) + +L = mu_er*(n**2)*(radius**2)*mu_0/len_rod //Inductance (H) +Q = 2*%pi*f*L/(abs(Rf) + Rr + RL) //Ratio of energy stored to energy lost per cycle (unitless) + +fHP = f/Q //Bandwidth at half power (Hz) + + +//Results +mprintf("The radiation efficiency is %.2e",k) +mprintf("\nThe value of Q is %.3f",Q) +mprintf("\nThe half-power bandwidth is %d Hz",fHP) |