diff options
Diffstat (limited to '2219/CH9')
-rwxr-xr-x | 2219/CH9/EX9.1/Ex9_1.sce | 15 | ||||
-rwxr-xr-x | 2219/CH9/EX9.10/Ex9_10.sce | 18 | ||||
-rwxr-xr-x | 2219/CH9/EX9.11/Ex9_11.sce | 16 | ||||
-rwxr-xr-x | 2219/CH9/EX9.12/Ex9_12.sce | 20 | ||||
-rwxr-xr-x | 2219/CH9/EX9.13/Ex9_13.sce | 22 | ||||
-rwxr-xr-x | 2219/CH9/EX9.14/Ex9_14.sce | 22 | ||||
-rwxr-xr-x | 2219/CH9/EX9.15/Ex9_15.sce | 19 | ||||
-rwxr-xr-x | 2219/CH9/EX9.16/Ex9_16.sce | 24 | ||||
-rwxr-xr-x | 2219/CH9/EX9.17/Ex9_17.sce | 11 | ||||
-rwxr-xr-x | 2219/CH9/EX9.2/Ex9_2.sce | 19 | ||||
-rwxr-xr-x | 2219/CH9/EX9.3/Ex9_3.sce | 25 | ||||
-rwxr-xr-x | 2219/CH9/EX9.4/Ex9_4.sce | 15 | ||||
-rwxr-xr-x | 2219/CH9/EX9.5/Ex9_5.sce | 15 | ||||
-rwxr-xr-x | 2219/CH9/EX9.6/Ex9_6.sce | 24 | ||||
-rwxr-xr-x | 2219/CH9/EX9.7/Ex9_7.sce | 10 | ||||
-rwxr-xr-x | 2219/CH9/EX9.8/Ex9_8.sce | 18 | ||||
-rwxr-xr-x | 2219/CH9/EX9.9/Ex9_9.sce | 16 |
17 files changed, 309 insertions, 0 deletions
diff --git a/2219/CH9/EX9.1/Ex9_1.sce b/2219/CH9/EX9.1/Ex9_1.sce new file mode 100755 index 000000000..19acedf0b --- /dev/null +++ b/2219/CH9/EX9.1/Ex9_1.sce @@ -0,0 +1,15 @@ +// Chapter 9 example 1
+//------------------------------------------------------------------------------
+clc;
+clear;
+// Given Data
+PRF = 1000; // Pulse repetitive frequency in Hz
+t = 0.15*10^-3; // Round propagation time in s
+c = 3*10^8; // velocity of EM waves in m/s
+// calculations
+R = (c*t)/2; // Range
+Runamb = c/(2*PRF) // Max unambiguous range
+
+// Output
+mprintf('Target Range = %3.1f Km\n Maximum Unambiguous range = %d Km',R/1000,Runamb/1000);
+//------------------------------------------------------------------------------
diff --git a/2219/CH9/EX9.10/Ex9_10.sce b/2219/CH9/EX9.10/Ex9_10.sce new file mode 100755 index 000000000..4fcb94fc2 --- /dev/null +++ b/2219/CH9/EX9.10/Ex9_10.sce @@ -0,0 +1,18 @@ +// Chapter 9 example 10
+//------------------------------------------------------------------------------
+clc;
+clear;
+// Given Data
+f = 10^9; // CW radar waveform freq.
+fm = 100; // modulation freq. in Hz
+MaxfD = 500; // max freq deviation in Hz
+c = 3*10^8; // vel. of EM waves in m/s
+
+// Calculations
+Mf = MaxfD/fm // Modulation index
+BW = 2*(Mf + 1)*fm // Bandwidth
+RR = c/(2*BW); // Range Resolution in m
+
+// Output
+mprintf('Bandwidth = %d Hz\n Range Resolution = %d Km',BW,RR/1000);
+//------------------------------------------------------------------------------
diff --git a/2219/CH9/EX9.11/Ex9_11.sce b/2219/CH9/EX9.11/Ex9_11.sce new file mode 100755 index 000000000..175eec14b --- /dev/null +++ b/2219/CH9/EX9.11/Ex9_11.sce @@ -0,0 +1,16 @@ +// Chapter 9 example 11
+//------------------------------------------------------------------------------
+clc;
+clear;
+// Given Data
+f = 10^9; // Centre freq. of spectrum
+t = 13 // pulse width in us
+N = 13; // N-bit Barker code
+
+// Calculations
+Sub_PW = t/N; // sub pulsewidth
+match_BW= 1/Sub_PW; // Matched bandwidth in Mhz
+
+// Output
+mprintf('Matched Bandwidth = %d Mhz\n Center Frequency of the spectrum = %d Ghz',match_BW,f/10^9 );
+//-----------------------------------------------------------------------------
diff --git a/2219/CH9/EX9.12/Ex9_12.sce b/2219/CH9/EX9.12/Ex9_12.sce new file mode 100755 index 000000000..2d7cba099 --- /dev/null +++ b/2219/CH9/EX9.12/Ex9_12.sce @@ -0,0 +1,20 @@ +// Chapter 9 example 12
+//------------------------------------------------------------------------------
+clc;
+clear;
+// Given Data
+PW = 10^-6; // Pulse width in sec
+Pp = 100*10^3; // Peak power in watts
+PRF = 1000; // pulse rep.rate
+N_target= 20; // no of target hits in 1 dwell period
+
+// Calculations
+PE = Pp*PW; // Pulse energy in Joule
+LE = N_target *PE; // look energy
+DC = PW*PRF // Duty cycle
+Pav = Pp*DC; // Average power
+Pavg = 10*log10(Pav); // Avg power in dB
+
+// Output
+mprintf('Average power = %d dB\n Look Energy = %3.0f Joules',Pavg,LE);
+//-----------------------------------------------------------------------------
diff --git a/2219/CH9/EX9.13/Ex9_13.sce b/2219/CH9/EX9.13/Ex9_13.sce new file mode 100755 index 000000000..066470771 --- /dev/null +++ b/2219/CH9/EX9.13/Ex9_13.sce @@ -0,0 +1,22 @@ +// Chapter 9 example 13
+// Data taken from Ex 12
+//------------------------------------------------------------------------------
+clc;
+clear;
+// Given Data
+PW = 10^-6; // Pulse width in sec
+Pp = 100*10^3; // Peak power in watts
+PRF = 1000; // pulse rep.rate
+N_target= 20; // no of target hits in 1 dwell period
+
+// Calculations
+PE = Pp*PW; // Pulse energy in Joule
+LE = N_target *PE; // look energy
+DC = PW*PRF // Duty cycle
+Pav = Pp*DC; // Average power
+Pavg = 10*log10(Pav); // Avg power in dB
+Pp_dB = 10*log10(Pp); // Peak power in dB
+DCCF = Pp_dB - Pavg // Duty cycle correction factor
+// Output
+mprintf('Duty cycle correction factor = %d dB',DCCF);
+//-----------------------------------------------------------------------------
diff --git a/2219/CH9/EX9.14/Ex9_14.sce b/2219/CH9/EX9.14/Ex9_14.sce new file mode 100755 index 000000000..9c73a26d0 --- /dev/null +++ b/2219/CH9/EX9.14/Ex9_14.sce @@ -0,0 +1,22 @@ +// Chapter 9 example 14
+//------------------------------------------------------------------------------
+clc;
+clear;
+// Given Data
+G_rx = 97; // Rx gain in dB
+Bn = 5*10^6; // Bandwidth in Hz
+To = 300; // temperature in kelvin
+K = 1.38*10^-23; // Boltzmann constant in J/k
+n = -3 // o/p noise power in dBm
+
+// calculations
+Pn_dB = n-G_rx // input noise power
+Pn = 10^(Pn_dB/10)*10^-3 // converting from dBm to watts
+// Pn = KToBnF;
+F = Pn/(K*To*Bn) // Noise Factor
+T = To*[F - 1] // Equivalent Noise Temperature
+
+// Output
+mprintf('Equivalent Noise Temperature = %d°K',T );
+//------------------------------------------------------------------------------
+
diff --git a/2219/CH9/EX9.15/Ex9_15.sce b/2219/CH9/EX9.15/Ex9_15.sce new file mode 100755 index 000000000..4b3f69cab --- /dev/null +++ b/2219/CH9/EX9.15/Ex9_15.sce @@ -0,0 +1,19 @@ +// Chapter 9 example 15
+//------------------------------------------------------------------------------
+clc;
+clear;
+// Given Data
+Gx = 60; // gain of Rx 'X' in dB
+Gy = 70; // gain of Rx 'Y' in dB
+Fx = 3; // Noise factor of 'X'
+Fy = 2; // Noise factor of 'Y'
+
+// calculations
+Gx_W = 10^(Gx/10) // gain in watts
+Gy_W = 10^(Gy/10) // gain in watts
+// k = Pnx/Pny; // Ratio of noise power levels produced at the o/p's of Rx 'X' and 'Y'
+k = (Fx*Gx_W)/(Fy*Gy_W); // Ratio of noise power levels produced at the o/p's of Rx 'X' and 'Y'
+
+// output
+mprintf('Ratio of noise power levels produced at the outputs of Rx X and Y = %3.2f',k);
+//------------------------------------------------------------------------------
diff --git a/2219/CH9/EX9.16/Ex9_16.sce b/2219/CH9/EX9.16/Ex9_16.sce new file mode 100755 index 000000000..34a4a0d85 --- /dev/null +++ b/2219/CH9/EX9.16/Ex9_16.sce @@ -0,0 +1,24 @@ +// Chapter 9 example 16
+//------------------------------------------------------------------------------
+clc;
+clear;
+// Given Data
+Pn = -70; // Noise power in dBm
+fl = 10^6; // lower cut-off freq in Hz
+fh = 11*10^6; // upper cut-off freq in Hz
+BP_fl = 13*10^6; // Bandpass filter lower cutoff freq
+BP_fh = 14*10^6; // Bandpass filter lower cutoff freq
+
+// Calculations
+Pn_W = 10^(Pn/10)*10^-3; // coversion from dBm to Watts
+BW = fh - fl
+PSD = Pn_W/BW // Noise power spectral density
+// Since white noise has the same spectral power density through the frequency spectrum,
+// therefore Noise power in second case
+B = BP_fh - BP_fl
+Pn_2 = PSD*B; // Noise power in second case
+
+// Output
+mprintf('Noise power for BandPass filter having Cutoff frequencies 13Mhz and 14Mhz = %3.0e W',Pn_2);
+//--------------------------------------------------------------------------------
+
diff --git a/2219/CH9/EX9.17/Ex9_17.sce b/2219/CH9/EX9.17/Ex9_17.sce new file mode 100755 index 000000000..1c980595b --- /dev/null +++ b/2219/CH9/EX9.17/Ex9_17.sce @@ -0,0 +1,11 @@ +// Chapter 9 example 16
+//------------------------------------------------------------------------------
+clc;
+clear;
+// Given Data from Figure triagle OAB
+OA = 100 // in Km
+OB = OA*cos(60*%pi/180); // Range of Target 2
+
+// Output
+mprintf('Range of Target-2 = %d Km\n Azimuth angle of target-1 = 60°\n Azimuth angle of Target-2 = 120°',OB);
+//------------------------------------------------------------------------------
diff --git a/2219/CH9/EX9.2/Ex9_2.sce b/2219/CH9/EX9.2/Ex9_2.sce new file mode 100755 index 000000000..0b2a8a1b5 --- /dev/null +++ b/2219/CH9/EX9.2/Ex9_2.sce @@ -0,0 +1,19 @@ +// Chapter 9 example 2
+//------------------------------------------------------------------------------
+clc;
+clear;
+// Given Data
+f = 10*10^9; // radar Tx frequency
+c = 3*10^8; // velocity of EM waves in m/s
+V = 108; // vel of car in kmph
+
+// Calculations
+lamda = c/f; // wavelength in m
+Vr = V*(5/18); // vel of car in m/s
+fd = (2*Vr)/lamda // Doppler shift in Hz
+fr = f + fd // received freq
+fr_away = f-fd // Rx frequency if the car is moving away from radar
+
+// Output
+mprintf('Doppler Shift = %d Khz\n Frequency of Received signal = %3.6f Ghz\n Received Frequency if car is moving away from radar = %3.6f Ghz',fd/1000,fr/10^9,fr_away/10^9);
+//------------------------------------------------------------------------------
diff --git a/2219/CH9/EX9.3/Ex9_3.sce b/2219/CH9/EX9.3/Ex9_3.sce new file mode 100755 index 000000000..47a0c422f --- /dev/null +++ b/2219/CH9/EX9.3/Ex9_3.sce @@ -0,0 +1,25 @@ +// Chapter 9 example 3
+//------------------------------------------------------------------------------
+clc;
+clear;
+// Given Data
+f = 10*10^9; // radar Tx frequency
+PRF = 2000; // Pulse repetitive frequency in Hz
+Vr = 0.5; // radial vel in Mach
+c = 3*10^8; // velocity of EM waves in m/s
+vs = 330; // velocity of sound in m/s
+
+// Calculations
+lamda = c/f; // wavelength in m
+max_unamb_fd = PRF/2; // maximum unambiguous doppler shift
+Vrunamb = (lamda*max_unamb_fd)/2; // doppler shift
+Vaircraft = 0.5*vs; // Converting from Mach to m/s
+fd_desired = (2*Vaircraft)/lamda;
+PRF_desired = 2*fd_desired; // desired PRF
+
+// Output
+if Vrunamb < Vaircraft then
+ mprintf('The radar is not capable of determining unambiguously the velocity of the approaching aircraft\n');
+end
+mprintf(' Desired Pulse Repetition Rate = %d Khz',PRF_desired/1000);
+//------------------------------------------------------------------------------
diff --git a/2219/CH9/EX9.4/Ex9_4.sce b/2219/CH9/EX9.4/Ex9_4.sce new file mode 100755 index 000000000..1b1967445 --- /dev/null +++ b/2219/CH9/EX9.4/Ex9_4.sce @@ -0,0 +1,15 @@ +// Chapter 9 example 4
+//------------------------------------------------------------------------------
+clc;
+clear;
+// Given Data
+PW_tx = 10^-6; // Transmitted pulse width
+Rx_PW = 10^-6; // Received pulse width
+c = 3*10^8; // velocity of EM waves in m/s
+
+// Calculations
+RR = (c*Rx_PW)/2; // Range Resolution in m
+
+// output
+mprintf('This Radar can resolve upto an inter target separation in range of %d m\n Therefore,given radar will be able to resolve the targets',RR);
+//------------------------------------------------------------------------------
diff --git a/2219/CH9/EX9.5/Ex9_5.sce b/2219/CH9/EX9.5/Ex9_5.sce new file mode 100755 index 000000000..89f8b36d9 --- /dev/null +++ b/2219/CH9/EX9.5/Ex9_5.sce @@ -0,0 +1,15 @@ +// Chapter 9 example 5
+//------------------------------------------------------------------------------
+clc;
+clear;
+// Given Data
+CRR = 100; // Cross range resolution in m
+R = 3000; // radial range
+
+// Calculations
+// CRR = (R*theta3)*(%pi/180);
+theta3 = (180*CRR)/(%pi*R) // 3 dB beamwidth
+
+// Output
+mprintf('3 dB beamwidth = %3.2f°',theta3);
+//------------------------------------------------------------------------------
diff --git a/2219/CH9/EX9.6/Ex9_6.sce b/2219/CH9/EX9.6/Ex9_6.sce new file mode 100755 index 000000000..04ed7d169 --- /dev/null +++ b/2219/CH9/EX9.6/Ex9_6.sce @@ -0,0 +1,24 @@ +// Chapter 9 example 6
+//------------------------------------------------------------------------------
+clc;
+clear;
+// Given Data
+Vs = 330; // velocity of sound in m/s
+NM = 1.85*(5/18) // 1NM equivalent in m/s
+V1 = 0.5; // velocity of first aircraft in mach
+V2 = 400; // velocity of second aircraft in NM/hr
+theta = 30; // angle with radial axis in degrees
+lamda = 3*10^-2; // wavelength in m
+
+// Calculations
+v1 = V1*Vs // velocity of first aircraft in m/s
+fd1 = (2*v1)/lamda // doppler freq.
+v2 = V2*NM*cos(30*(%pi/180)) // velocity of second aircraft in m/s
+fd2 = (2*v2)/lamda // doppler freq
+dd = fd2 - fd1 // doppler difference
+Tl = 1/dd // look time in s
+
+// Output
+mprintf('Required minimum look time = %3.2f ms',Tl/10^-3);
+mprintf('\n Note: Cos(30) value is taken as 0.5 in textbook');
+//------------------------------------------------------------------------------
diff --git a/2219/CH9/EX9.7/Ex9_7.sce b/2219/CH9/EX9.7/Ex9_7.sce new file mode 100755 index 000000000..973587884 --- /dev/null +++ b/2219/CH9/EX9.7/Ex9_7.sce @@ -0,0 +1,10 @@ +// Chapter 9 example 7
+//------------------------------------------------------------------------------
+clc;
+clear;
+// Given Data
+// Rmax = [1000000/(12.4*PRF)]NM
+// = [1000000*t/12.4]NM
+mprintf('The Numerator represents round trip propagation time in us\n');
+mprintf(' Therefore, number 12.4 represents the units microseconds per nautical miles\n');
+mprintf(' In other words, this means that the round propagation time for one nautical mile is 12.4 us which is equivalent to 6.66us for 1km range')
diff --git a/2219/CH9/EX9.8/Ex9_8.sce b/2219/CH9/EX9.8/Ex9_8.sce new file mode 100755 index 000000000..7232cd01c --- /dev/null +++ b/2219/CH9/EX9.8/Ex9_8.sce @@ -0,0 +1,18 @@ +// Chapter 9 example 8
+//------------------------------------------------------------------------------
+clc;
+clear;
+// Given Data
+PW = 10*10^-6; // pulse width in sec
+f = 10*10^9; // frequency in Hz
+fm = 1000; // modulating frequency
+
+// calculations
+BW_M = 1/PW // matched bandwidth
+cf1 = f+fm; // closest freq.
+cf2 = f-fm; // closest freq.
+fo = f; // centre freq.
+
+// Output
+mprintf('Centre of frequency spectrum = %d Khz\n The two closet frequencies to the center of the spectrum are %d Khz and %d Khz',fo/10^3,cf1/10^3,cf2/10^3);
+//------------------------------------------------------------------------------
diff --git a/2219/CH9/EX9.9/Ex9_9.sce b/2219/CH9/EX9.9/Ex9_9.sce new file mode 100755 index 000000000..3fe39385e --- /dev/null +++ b/2219/CH9/EX9.9/Ex9_9.sce @@ -0,0 +1,16 @@ +// Chapter 9 example 9
+//------------------------------------------------------------------------------
+clc;
+clear;
+// Given Data
+fc1 = 495; // freq in Mhz
+fc2 = 505; // freq in Mhz
+
+// Calculations
+fo = (fc1 + fc2)/2; // Center of spectrum in Mhz
+BW = fc2 - fc1; // Bandwidth in Mhz
+PW = 1/BW; // compressed pulse width in us
+
+// Output
+mprintf('Center of spectrum = %d Mhz\n Matched Bandwidth = %d Mhz\n Compressed Pulse width = %3.1fus',fo,BW,PW);
+//------------------------------------------------------------------------------
|