ESSA TR ERL 172-ITS 110 A UNITED STATES DEPARTMENT OF COMMERCE PUBLICATION "V,,*"** ESSA Technical Report ERL 172-ITS 110 U.S. DEPARTMENT OF COMMERCE Environmental Science Services Administration Research Laboratories Modulation Studies for IGOSS IIROSHI AKIMA BOULDER, COLO. JUNE 1970 ESSA RESEARCH LABORATORIES The mission of the Research Laboratories is to study the oceans, inland waters, the lower and upper atmosphere, the space environment, and the earth, in search of the under- standing needed to provide more useful services in improving man's prospects for survival as influenced by the physical environment. Laboratories contributing to these studies are: Earth Sciences Laboratories: Geomagnetism, seismology, geodesy, and related earth sciences; earthquake processes, internal structure and accurate figure of the Earth, and distribution of the Earth's mass. Atlantic Oceanographic and Meteorological Laboratories: Oceanography, with emphasis on the geology and geophysics of ocean basins, oceanic processes, sea- air interactions, hurricane research, and weather modification (Miami, Florida). Pacific Oceanographic Laboratories: Oceanography; geology and geophysics of the Pacific Basin and margins; oceanic processes and dynamics; tsunami generation, propaga- tion, modification, detection, and monitoring (Seattle, Washington). Atmospheric Physics and Chemistry Laboratory: Cloud physics and precipitation; chem- ical composition and nucleating substances in the lower atmosphere; and laboratory and field experiments toward developing feasible methods of weather modification. Air Resources Laboratories: Diffusion, transport, and dissipation of atmospheric con- taminants; development of methods for prediction and control of atmospheric pollution (Silver Spring, Maryland). Geophysical Fluid Dynamics Laboratory: Dynamics and physics of geophysical fluid systems; development of a theoretical basis, through mathematical modeling and computer simulation, for the behavior and properties of the atmosphere and the oceans (Princeton, New Jersey). National Severe Storms Laboratory: Tornadoes, squall lines, thunderstorms, and other severe local convective phenomena toward achieving improved methods of forecasting, detecting, and providing advance warnings (Norman, Oklahoma). Space Disturbances Laboratory: Nature, behavior, and mechanisms of space disturb- ances; development and use of techniques for continuous monitoring and early detection and reporting of important disturbances. Aeronomy Laboratory: Theoretical, laboratory, rocket, and satellite studies of the physical and chemical processes controlling the ionosphere and exosphere of the earth and other planets. Wave Propagation Laboratory: Development of new methods for remote sensing of the geophysical environment; special emphasis on propagation of sound waves, and electro- magnetic waves at millimeter, infrared, and optical frequencies. Institute for Telecommunication Sciences: Central federal agency for research and services in propagation of radio waves, radio properties of the earth and its atmosphere, nature of radio noise and interference, information transmission and antennas, and meth- ods for the more effective use of the radio spectrum for telecommunications. Research Flight Facility: Outfits and operates aircraft specially instrumented for re- search; and meets needs of ESSA and other groups for environmental measurements for aircraft (Miami, Florida). ENVIRONMENTAL SCIENCE SERVICES ADMINISTRATION BOULDER, COLORADO 80302 . a. o u o '8 Q. D } = \ f fp( D > U) dD dU, (1) U>D where Pr -| [ stands for "the probability that," and p(D, U) is the probability density function of D and U. This applies to an NCFSK system where the modulation index is not less than unity and no low- pass filter is used before the decision-making circuit. For proper application of (1), the statistics of the undesired signal must be measured at the input of the demodulator, i. e. , at the input to the limiter in a limiter -discriminator demodulator, and in a bandwidth equivalent to the sum of the bandwidths of the two filters in a dual- filter demodulator. Except for a few simple cases, such as a non- fading signal or a Ray leigh -fading signal in the presence of Gaussian noise, this integral cannot be evaluated analytically but must be computed numerically. The element error probability in a binary TDPSK system is closely approximated by (1); it is exactly equal to the result computed by (1), when the undesired signal consists of Gaussian noise (Cahn, 1959)- (Since the bandwidth for a TDPSK system is usually considered one -half that for an NCFSK system, its required ratio of signal energy per bit to noise power density can be 3 dB lower than that for NCFSK to give an equal error probability. ) Amplitudes of signals propagated by the ionosphere vary con- siderably with time. For short intervals (3 to 7 min), amplitude distribution functions close to the Rayleigh distribution predominate. Over longer intervals (30 to 60 min), on the other hand, amplitude distributions more often follow a log -normal law. Although the form of the measured distribution may differ from the Rayleigh distribution, the observed fading range, defined as the ratio of the upper and lower deciles, is the same order as 13.4 dB expected for the Rayleigh dis- tribution (CCIR, 1967g). In HF channels noise is mostly atmospheric. The amplitude - probability distribution (APD) of atmospheric noise can be represented, accurately enough for most applications, by an appropriate curve chosen from a family of idealized curves. The choice can be made by specifying a single parameter, defined as the ratio of the rms to the average of the envelope voltage and denoted by V in dB (CCIR, 1964). For simplicity we assume that the interfering signal is either a continuous wave (CW) or an FSK or PSK signal. The element error probability in each particular case is generally given as a function of signal-to-noise ratio (SNR) and /or signal-to-interference ratio (SIR). The SNR is defined as the ratio of desired signal power to average noise power, and SIR is the ratio of desired signal power to interfering signal power, where both the noise power and interfering signal power are measured at the demod- ulator input. For a fading signal, either desired or interfering, its median power is used in defining the SNR and SIR. The integral in (1) has not been computed for very general cases, such as a log-normal-fading signal in the presence of atmo- spheric noise and a log -normal- fading interfering signal with dual diversity reception. Some particular cases, however, have been studied and described below. 3. 1. Nonfading Signal and Atmospheric Noise Without Interference The element error probability is one -half the value of the APD of the atmospheric noise corresponding to the signal amplitude. The relation between SNR and error probability, parametric in noise param- eter V, (Akima et al. , 19 69), is shown in figure 1. d 3.2. Fading Signal and Atmospheric Noise Without Interference (Nondiversity Reception) The relation between SNR and element error probability, parametric in noise parameter V.» computed by numerical integration d (Akima et al. , 1969), are given in figure 2 for Rayleigh-fading signals and in figure 3 for log-normal-fading signals (fading range = 13. 4 dB). 10 10 r2 o w i_ UJ c E UJ O o .4 -§ io 4 0_ 10 r5 10 r6 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 | 1 1 1 1 Nonfading Signal Atmospheric Noise No Interference 1 1 1 — — \ : — S^- IO dB — V d =1.049 dB » I 2 dB ^n , 4 dB-A \ ^ \X-8 dB 1 1 1 II 1 1 II I I l\ 1 1 l\ VAX _6dB 1 1 1 — \\\ •10 10 20 30 SNR (dB) 40 50 Figure 1. Element error -probability in a single «- channel NCFSK system under stable (nonfading) conditions with atmospheric noise and no interference . 20 SNR (dB) Figure 2. Element error -probability in single- channel NCFSK system under Ray leigh- fading conditions with atmospheric noise and no inter ference > nondiversity reception. 10" fc 10" Itoi c CD E 10" o -O O -O o 10" 10" 10" 1 1 1 1 I I I l Log- Normal - i i i i Fading Signal 1 1 1 1 — ^^^ ( Fading Range = l3.4dB) Atmospheric Noise No Interference Nondiversity — — — — v ^-~-Vd = 10 dB — 4dB 2dB-— V d =1.049 dB- \\V^^-— - 8 dB _ — 6dB _ i i i r i i i i 1 1 1 1 i i Ai \\ 1 1 •10 10 20 SNR (dB 30 40 50 Figure 3. Element error probability in a single- channel NCFSK system under log-normal-fading conditions with atmospheric noise and no interference, nondiversity reception. (The fading range is assumed to be 13.4 dB, which is the same as expected for Rayleigh fading.) 10 3.3. Fading Signal and Atmospheric Noise Without Interference (Diversity Reception) When dual diversity with a selection-switching combiner is used, the probability of element error in an NCFSK system is given by P g = - Pr { U a > D x , V a > V 2 } + \ Pr { U 2 > D ? , V 2 > V x }, (2) where D. = voltage of desired signal from antenna i or in channel i (i = 1 , 2), U. = voltage of undesired signal (noise) from antenna i or channel i (i = 1 , 2), V. = vector sum of D. and U. l 11 = voltage of signal plus noise from antenna i or in channel i (i = 1 , 2). Since D 8 and U 2 follow the same probability distributions as D-l and \J 1 , respectively, the second term in (2) is equal to the first term. Therefore, -we have P e = Pr { U x > D lf V a > V s } (3) as an expression of element error probability in an NCFSK system with dual selection-switching diversity reception. This result can easily be extended to cover higher order diver sity with selection-switching combining. For example, 11 Pe = 2Pr {u x > D 1 , V 1 > V 2 , V, > V 3 , V x > V 4 | (4) is the expression of el ment error probability for quadruple diversity- reception. Evaluating P from (3) or (4) is the evaluation of a definite integral of a joint probability density function and can best be done by the Monte Carlo method or numerical simulation. One big advantage of the Monte Carlo method is that it does not become much more complicated when the order of diversity is increased. When we use this method, sequences of random numbers are chosen to follow the probability distribution functions of the signal, the noise, and the phase difference between the signal and the noise. The necessary condition for the occurrence of an error, as given by (3) or (4), is then tested with successive sets of elements; each element is picked from the corresponding sequence of random numbers, and the number of times an error occurs is counted. We assumed statistical independence among the signal and noise amplitudes and a random phase relationship between signal and noise. We then tested one million sets of elements for each combination of two types of signal fading, dual and quadruple diversity, six values of V , and several values of SNR in 2-dB steps. Results of these compu- tations for dual diversity are shown in figures 4 and 5 for Rayleigh fading signal and log -normal -fa ding signal, respectively, and for quadruple diversity in figures 6 and 7. Values of the required SNR for element error probabilities of 10 and 10 are read from figures 2 through 7 and shown in table 1. This table indicates that the gain in required SNR obtained either by the use of dual diversity instead of nondiversity or by the use of quadruple diversity instead of dual diversity depends on the type of fading, the value of V , and the allowable error 12 io- 5 I0" 2 c E > X> o XI o 10" 10" 1 1 1 1 1 1 1 1 1 1 1 1 MM MM Rayleigh - Fading Signal — — v Atmospheric Noise No Interference Dual Diversity (Selection Switching) — — — — 2dB^ ^V d = 10 dB — Vd = l.049dB^ 5dB _ -6dB 1 1 1 1 1 1 1 1 1 1 1 1 \ \ \ \ \ 1 1 1 1 \ \\\ \ \ \ V \ \ \N \ N N \ W \N 1 1 ■10 10 20 SNR (dB 30 40 50 Figure 4. Element error probability in a single-channel NCFSK system under Rayleigh -fading conditions with atmospheric noise and no interference, dual (selection-switching ) diver- sity reception. 13 20 30 SNR (dB) Figure 5. Element error probability in a single-channel NCFSK system under log-normal-fading oonditons with atmospheric noise and no interference* dual (selection-switching) diver- sity reception. (The fading range is assumed to be 13.4 dB . ) 14 10" o I0' 2 c V E u 10- 3 O -O O £ 10" 10 -5 10" 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 MM Rayleigh - Fading Signal Mil — Atmospheric Noise No Interference Quadruple Diversity (Selection Switching) — I \ — — 4dB' \ "*\ -v d = io de - 8dB — — 2dE V d = 1.049 IB — "*\ * \ \ \ V— 6dB — 1 1 1 1 1 1 1 1 \ \ \ 1 1 1 l\ \ \ \ \ \ \ 1 \ 1 1 ^ \ \ X 1 \\\ \l 1 1 i ■10 10 20 30 SNR (dB) 40 50 Figure 6. Element error probability in a single- channel NCFSK system under Ray leigh- fading conditions with atmospheric noise and no interferences quadruple (selection-switching) diversity reception. 15 10" o I0" 2 UJ c E LxJ 10" JO o O CL 10" 10 -5 10" 1 1 1 1 1 1 1 1 1 1 1 1 1 II 1 1 1 1 1 1 1 1 "^^^^^^^v Log- Normal - Fading Signal — (Fading Range = 13.4 dB ) Atmospheric Noise No Interference — Quadruple Diversity — ( Selection Switching ) — — ^V d = 10 dB — 4dB- ^^8dB 2dB—"" YV^^-6dB V d = 1.049 UJ o jO o o 10 -3 10 H 10 -5 10 1 1 1 1 1 1 1 1 1 II 1 [III fill Nonfading Signal Nonfading Interference Atmospheric Noise (V d = 10 dB) 1 1 1 1 — ^^s^ — — .SIR = 1 dB — — ^____ 2dB — — 3dB — — \^6dB — — 10 dB y>\ \ ~~ — OCK \ \ — 1 1 II 1 1 1 1 1 1 1 1 1 1 1 1 1 1 l\\ \ l\ \ -10 10 20 SNR (dB) 30 40 50 Figure 11. Element error probability in a single- channel NCFSK system under stable conditions with atmospheric noise (V^ = 10 dB) and nonfading interference. 22 compared with 9, 5, and 2 dB for Gaussian noise. The corresponding values are 6, 3, and 1 dB for V , = 6 dB, and 6, 2, and 1 dB for d V J = 10 dB. d 3. 6. Fading Signal and Fading Interference Without Noise (Nondiversity Reception) When both the desired signal and the interference are Rayleigh- fading the element error probability is the same as that for a Rayleigh' fading desired signal in the presence of Gaussian noise and no interference, because the APD of the Rayleigh -fading interference is the same as that of Gaussian noise. For a Rayleigh -fading signal and Gaussian noise, Montgomery (1954) showed that the element error probability is given by P e-I ■ TTR- ' < 6 » where R is the ratio of average signal power to average noise power. This relation applies to our case, if R is interpreted as the SIR. (Although the SIR is defined here as the ratio of median desired signal power to median interfering signal power, it is equal to the ratio of average desired signal power to average interfering signal power when both signals follow an identical distribution law. ) This relation for nondiversity is plotted in figure 12, section 3.7. For log -normal -fading signals, desired and interfering, a linear transformation of the coordinate system can reduce the double integral in (1) to a single integral of the form P e = I . Ff __2 LI . (7) 23 where 00 1 f t 2 V2TT J F(x) = — / exp(--)dt, (8) X m and m. are the median levels in dB of the desired and interfering si signals, respectively, and a and o. are the respective standard devi- ations, also in dB. The o . is related to the fading range R by S f X X R = 2.563o. (9) The relation between SIR and element error probability for log-normal fading (fading range = 13.4 dB) and nondiversity is shown in figure 13, section 3.7. 3.7. Fading Signal and Fading Interference Without Noise (Diversity Reception) We start with quadruple space and frequency diversity with a selection-switching combiner: two modulated signals (both correspond- ing to an identical digital signal) are transmitted simultaneously in two channels separated in frequency, channel 1 and channel 2, and are received on two antennas separated in space, antenna 1 and antenna 2; the strongest signal among the four received signals is selected at the diversity combiner. The probability of element error in an NCFSK system is then given by P =-PrD , V = max (V ) } e ci \ a^c^ a^c^ a^c^ a^c-j ' + -Pr{u >D , V = max(V . .) } c, \ a^Cg a i c s 3-i c s a i c j ' 24 + ^-Priu >D , V = max (V ) > + ^-Pr|u >D , V = max(V ) > , 2 { &2 C 2 a.2Cg &2 C 2 a i c j (10) where D_ = voltage of desired signal from antenna i in channel j (i = 1, 2; j = 1, 2), a i c J U = voltage of unde sired signal (interference) from antenna i in channel j (i = 1, 2; j = 1, 2), a i c j V = vector sum of D and U a^c-i ii \ 1 = voltage of signal plus interference from antenna i in channel j (i = 1, 2; j = 1, 2), and max (V ) is a voltage having a maximum amplitude among the a i c j four voltages, V ., V , V , and V . From this general a l c l 3-1 c S a 2 c l 3-aCg form we shall derive three useful relations. First, we assume quadruple space and frequency diversity- reception with an interfering signal in channel 1 only. Since no inter- fering signal exists in channel 2, both U and U are zero and a i c a &-z c 2 cannot exceed D and D , respectively; therefore, the second a i c a a 2 c 2 and the fourth terms of (10) are zero. Since, in general, D and a 2 c l U follow the same probability distributions as D and U , a 2 c 1 a 1 c l a^ respectively, the third term of (10) should be equal to the first term. Therefore, we have 25 P = Pr < U >D , V = max (V ) > G \ cl^C^ ^*1^1 ^1^1 1^1 = Pr < U >D , V > D , V a^C]_ cL-j^C^ cL^C^ &]_Cg V > V , V > D > (11) a^C^ 3-3^1 3-1^1 3-2^2 ' as an expression of element error probability in this case. Second, simplifying the reception scheme, we assume dual frequency diversity reception with an interfering signal in channel 1 only. Since we do not have antenna 2, the third and the fourth terms of (10), based on voltages from antenna 2, are zero. In addition, since no interference is present in channel 2, the second term of (10) is zero. Therefore, we obtain P = - Pr {u > D , V = max(V ) > e Li ' a^Ci a^c^ a^c^ a;c^ " = - Pr { U > D , V > D > for this case. (12) Finally, we consider dual space diversity with an interfering signal. Since we do not have channel 2, the second and the fourth terms of (10), based on voltages in channel 2, are zero. In addition, as in quadruple diversity reception, the third term of (10) is equal to the first term; therefore, we have P = Pr < U > D , V = max(V ) [ = Pr fu >D , V > V } \ a^Ci a^c^ a^c^ agC^ * (13) for this case. 26 Using (11), (12), and (13), we evaluated the probability of element error for quadruple space and frequency diversity, dual frequency diversity, and dual space diversity, respectively. A Monte Carlo method was used for the evaluation, which simulated error performance tests with a total of 10 6 bits for each type of diversity. Results of these evaluations are shown for Rayleigh fading in figure 12 and for log -normal fading in figure 13. They are also compared with the previously derived results for nondiversity reception. Note that the gain in required SIR expected from the use of diversity depends on the type of fading of signal and interference. Only a small gain can be expected when both desired and interfering signals are fading with a log -normal distribution; for example the gain expected from quadruple space and frequency diversity instead of nondiversity is 6 dB at an element error probability of 10" . 3. 8. Fading Signal and Fading Interference With Noise Except for some special cases, the element error probability must be computed by numerically integrating the joint probability den- sity function of the desired signal, interfering signal, and noise. Per- haps, the Monte Carlo method would be useful for this and might be mandatory for diversity reception. Even with the Monte Carlo method, the computation is expensive and time consuming, if we try to obtain data for numerous combinations of input parameters, such as types of fading, diversity techniques, noise parameter, SNR, and SIR. For Rayleigh -fading desired and interfering signals and Gaussian noise with nondiversity reception, the element error probability can be calculated analytically. Since the probability distribution function of the composite wave of a Rayleigh -fading interference and Gaussian noise is the same as that of Gaussian noise, the element error probability can be computed by (6). In this case, however, we must use, 27 UJ c