GEO 1145 ANNALS OF GEOPHYSICS, SUPPLEMENT TO VOL. 47, N. 2/3, 2004 13 A way to increase the bit rate in ionospheric radio links CLENCY PERRINE (1), YVON M. ERHEL (2), DOMINIQUE LEMUR (1), LOUIS BERTEL (1) and ALAIN BOURDILLON (1) (1) Institut d’Electronique et Télécommunications de Rennes, Université de Rennes 1, France (2) Institut d’Electronique et Télécommunications de Rennes, Centre de Recherches des Ecoles de Coëtquidan, Guer, France The aim of this paper is to present a high data rate transmission system through the ionospheric channel in the HF-band (3-30 MHz). The applications expected in this study are image transmitting and real-time videoconferencing. Very high rates are required compared to the standard modems. Therefore, an array processing is performed with a set of antennas whose spatial response differs from one another arranged in a circular array or in a collocated sensor. Synchronization (Zero Cross- ing Detector) and source separation (LMS algorithm) resort to classical well-tested techniques in- volving training sequences. Experimental results are presented for both antenna configurations. These techniques improve data rate, reaching 20 kbits/s within the 6 kHz bandwidth (QAM 64) with- out coding or interleaving. 13.1. INTRODUCTION HF-waves (3-30 MHz), when propagated through the ionosphere, can achieve very long dis- tance transmissions with a minimal infrastructure compared to satellite links, for example. This multipath multimode channel strongly degrades transmissions (mainly fading and frequency se- lectivity). For this reason, the data rate of standard scalar HF modems does not exceed 4.8 kbps in 3 kHz bandwidth. The originality of this study is to use a heterogeneous array to improve the HF transmission. Indeed, it has been shown that such a device could achieve the direction finding regarding the in- coming polarization as a decorrelation factor and, consequently, that the diversity of the spatial responses could, to some extent, replace space diversity (Erhel et al., 1998). The array processing is performed on a set of four active collocated antennas or a circular ar- ray and the corresponding coherent receiving channels. The spatio-temporal equalization resorts to the LMS (Least Mean Square) algorithm as the synchronisation is based on a «Zero Crossing Detector». In order to introduce the heterogeneous array and to bring its advantages, it is necessary to start with a short description of ionospheric propagation and the effects on transmission. Then it is shown that a heterogeneous array can improve the HF transmission and the real transmission system is described. Finally, experimental results are provided to illustrate the theory by practi- cal examples. 1146 Clency Perrine, Yvon M. Erhel, Dominique Lemur, Louis Bertel and Alain Bourdillon 13.2. PROPAGATION OF HF-WAVES IN THE IONOSPHERE 13.2.1. Multipath and multimode Two main phenomena must be considered in ionospheric propagation. The first is that waves can follow several distinct paths by reflecting themselves against different layers of the ionosphere. Then the transmitted signal gives birth to as many replicas that are affected in different ways by the chan- nel: their differential time delays can be large, and their Doppler shifts and angles of arrival can be spaced out too. The second phenomenon appears because of the Earth’s magnetic field that makes the medium anisotropic. Within each path, the wave splits into two complementary propagation modes called O (ordinary) and X (extraordinary). Both are very close and therefore have similar time delays, Doppler shifts and DOA’s (Direction Of Arrival). On the other hand, their polarization properties are very different. In the plane orthogonal to the direction of propagation, their electromagnetic fields de- scribe orthogonal ellipses clockwise and counterclockwise (Budden, 1985). Multipath and multi- mode propagation are illustrated in fig. 13.1. 13.2.2. Effects on transmission All the transmitted signal replicas created by multipath and multimode propagation interfere, de- grading transmission quality. Again, two cases can be distinguished that do not have the same con- sequences, although they are of the same nature. If two signals are delayed by an interval smaller than 1/B (where B is their bandwidth), they interfere in a way that the resultant signal fades, making the Signal to Noise Ratio (SNR) vary in time (an example is given in fig. 13.2). In the opposite case, the SNR remains constant but frequency selectivity occurs. Fading is due to multimode propagation. Only the O- and X-modes of the same path can be close enough to cause it. It is the most important problem for transmitting data since the signal is suscep- tible to completely disappear for a while. Concerning frequency selectivity, it is caused by modes coming from distinct paths. When several modes are present, both phenomena may be observed. The time-frequency analysis of the signal displayed in fig. 13.2 is plotted in fig. 13.3. It represents the mapping of its power in the time-frequency plane, and is obtained directly through the observation of the broadband modulated signal (Bisiaux, 2001). Two sets of interference patterns are clearly visi- ble, generating both fading and frequency selectivity. They denote the presence of at least 3 modes, among which 2 come from the same path (the line spectrum that appears every 4 s is not a computa- tion aberration, it is due to a periodic pattern inserted into the transmitted signal). Fig. 13.1. Multipath and multimode propagation in the ionosphere. 1147 A way to increase the bit rate in ionospheric radio links 13.2.3. Improvement of HF transmissions resorting to arrays of receiving antennas 13.2.3.1. Classical arrays Most of the time, arrays are made of identical elements set distant from one another (they are ho- mogeneous arrays). Space diversity allows them to steer their main lobe toward the DOA of a desired signal, or to create notches in the direction of undesired ones. However, such arrays would not be a very good solution for HF communications applications. The reason for that is illustrated in fig. 13.4: they would be efficient to combat multipath effects, but would not be able to separate the com- plementary modes of a single path because they usually leave the ionosphere at very close angles. The spatially filtered signal would then still be affected by fading. Furthermore, given that the wave- length ranges from 10 to 100 m, such arrays would spread over a large area, making them rather im- practical. For these reasons, homogeneous arrays are rarely used for transmission purposes in the HF-band. Yet it is possible to suppress fading by separating the O- and X-modes, but it requires in- cluding some kind of polarization treatment in the whole process, which is not possible with homo- geneous arrays. Fig. 13.2. Example of a signal fading (QAM-16, 3 kHz bandwidth, F~ 8 MHz). Fig. 13.3. Time-frequency analysis of an HF signal at Intermediate Frequency IF = 2.25 kHz. 13.313.2 13.4 13.5 Fig. 13.4. Source separation with a classical array. Fig. 13.5. Source separation with a heterogeneous array. 1148 Clency Perrine, Yvon M. Erhel, Dominique Lemur, Louis Bertel and Alain Bourdillon 13.2.3.2. Heterogeneous arrays With antenna diversity, it is possible to join the ability of an array to alter its radiation pattern on the one hand, and polarization filtering on the other hand. This leads to the realization of a hetero- geneous array, that is to say, made of antennas with various shapes or orientations, as opposed to classic homogeneous arrays. Its application to HF transmissions can be illustrated in a simplified way by the drawing in fig. 13.5: the radiation pattern of the array is not the same for all the waves because it depends on their polarization. Since O- and X-modes coming out of the ionosphere have different polarization properties, two different radiation patterns must be considered instead of one, as is usually the case. One regards O-modes, and is completely unrelated to the other one that con- cerns X-modes. For this reason it is theoretically possible to separate two complementary modes, even if they come from the same direction, and still being able to create notches in the two radiation patterns. Such a principle is ideal for combating both multipath and multimode propagation. As a consequence, the steering vector of a heterogeneous array is expressed as a F e F e F e / ( ) ( ) ( ) hO X j j NC j 1 2 2 2 2 / / / O X O X O X $ $ $ f= { i { i { i R T S S S S V X W W W W where the ϕκ (θ) are the geometrical phases of a classical array and FnO/X is the spatial response of an- tenna n for the two expected polarizations. 13.2.3.3. Expressions of the spatial correlation For an array of isotropic sensors the spatial correlation of two incident sources with DOA θ1 and θ2 is expressed as 1 NC e a a a a H j NC 1 2 1 2 n n2 1 $ = =t i i i i - = { i { i n ! ^ ^ ^ ^ ] ] h h h h g g8 B Using a heterogeneous array with the same geometry, the expression of the spatial correlation be- comes in a general case (whatever is the incident polarization) a a a a F a F a F a F a h h H h H 1 2 1 2 1 1 2 2 1 1 2 2 $ 7 $ 7 7 7 = =t i i i i i i i i i i i i ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ h h h h h h h h h h h h8 B where F(θ) is the vector of the spatial responses for the DOA θ. This quantity can be written under the following form: 1 1 1 e F F F F h n NC n NC n n j NC 1 2 2 2 1 2 n n2 1 $ $ =t i i i i) = = - = { i { i n n n ! ! ! ^ ^ ^ ^ ] ] h h h h g g8 B The modulus of the spatial correlation is modified for this type of array: ρh≠ρ. 1149 A way to increase the bit rate in ionospheric radio links 13.2.3.4. Improvement due to the heterogeneous array A comparison of the two quantities ρh and ρ does not provide a general result: the hetero- geneity of the array does not systematically decrease the spatial correlation. However, it can be men- tioned that current situations of reception are characterized by very close angles of arrival though the spatial responses are significantly different for the DOA θ1 and θ2: it is the case with diversely po- larized electromagnetic waves impinging from the same direction. The geometrical phases are nearly equal in these conditions (ϕn(θ1) ≅ ϕn(θ2) ) and the modulus of the spatial correlation for the homogeneous array is close to 1: ρ ≈1. On the other hand, considering the correlation for the heterogeneous array, the modulus ρh ap- proaches the value 1 1 1 F F F F h n NC n NC n n NC 1 2 2 2 1 2 $ $ ,t i i i i) = = = n n n ! ! ! ^ ^ ^ ^ h h h h . This ratio is less or equal to one, according to the Schwartz inequality: the spatial correlation de- creases in this case. 13.2.4. Numerical simulation In the following simulation, HF radio waves are assumed to be received after propagation through the ionosphere. It involves two uniform circular arrays made up with 8 loop antennas (figs. 13.6 and 13.7). This geometry has a practical advantage since the accuracy of the angular estima- tion is not dependant on the actual azimuth of arrival in this case. The first array is made up with 8 identical sensors (homogeneous array) and its spatial correlation is ρ. In the second one, the anten- nas are submitted to a rotation of 30° around a vertical axis every two positions within the array, so Fig. 13.6. Homogeneous circular array. Fig. 13.7. Heterogeneous circular array. 13.6 13.7 1150 Clency Perrine, Yvon M. Erhel, Dominique Lemur, Louis Bertel and Alain Bourdillon that it becomes heterogeneous with a spatial correlation ρh. This angular value realizes a compro- mise between a significant diversity of the spatial responses and the existence of redundancy if the rotation exceeds 180°. The spatial responses of the antennas are computed according to the model described in Section 13.2. The numerical parameters are: radius R = 0.5 λ, λ denoting the wave- length, and carrier frequency fo =15 MHz. Two incident waves are impinging with a common elevation and orthogonal polarizations. The first signal has a fixed azimuth of arrival (Az1=180°), and the second azimuth Az2 varies in the range of 0° to 360°. The simulation compares the spatial correlation of this couple of incident waves (con- sidered as a function of the variable azimuth Az2) for the homogeneous and the heterogeneous array. Whatever the method of direction finding, the less is the modulus of the spatial correlation, the easi- er is the separation of the signals. The critical case corresponds to a modulus close or equal to one for which the separation fails. Figure 13.8 plots the spatial correlations for both arrays (modulus of ρ and ρh) as function of the variable azimuth Az2. The benefit provided by the heterogeneous array appears in an interval of Az2 containing the value of Az1=180°. In that area, the modulus ρh is significantly smaller than ρ. In particular, in the critical case of two superimposed incident waves (Az2 = Az1=180°), the mod- ulus ρh is equal to 0.82 and ρ to 1: the heterogeneous array can separate the incident waves (with orthogonal polarizations) while the homogeneous array cannot. ρh is not less than ρwhat- ever the azimuth, but when it exceeds this reference value, the corresponding values are significant- ly smaller than 1, so that the DOA estimation works correctly. 13.3. THE TRANSMISSION SYSTEM The link used for the experiments has been established between Poitiers (latitude: 46°35′N, longitude: 0°25′E) and Monterfil (latitude: 48°03′N, longitude: 2°00′O) near Rennes (both in Fig. 13.8. Comparison of the spatial correlation of signals. 1151 A way to increase the bit rate in ionospheric radio links France). It is about 300 km long. Figure 13.9 gives a schematic view of the link and of the devices employed. This is a short link for an HF transmission that is normally used to cover several thou- sands of kilometers. 13.3.1. The transmitting system The transmitting system is located in Poitiers. In order to provide a good flexibility, it was decid- ed that the whole digital signal processing would be done with a computer: digital modulation and filtering, possibly interleaving and coding, etc. This way it is possible to control very easily the type of modulation used (PSK, QAM, OFDM, …), the sampling frequency and the signal’s bandwidth. The system is composed of a personal computer, a complex modulator, an HF power amplifier con- nected to a delta antenna. The complex modulator’s architecture is shown in fig. 13.9. Its role is simply to convert digital data sent in real-time by the PC into analog signals, through the means of DACs (Digital to Analog Converters), and then to transpose them into the HF-band. The operation is complex in that both real and imaginary parts of the complex baseband signal (I and Q channels respectively) are sent to the modulator, which also generates a complex carrier with a DDS (Direct Digital Synthesiser). Let x(t) = a(t) + jb(t) be the complex baseband signal. The synthesized signal can be expressed as ( ) ( ) ( ) ( ) ( )cos sins t a t t b t t= -~ ~ . which represent the real part of the equivalent complex signal (called analytical signal) ( ) ( )s t x t e j t 0= ~ # -. The carrier frequency is programmable through the PC, from 0 to 13 MHz by 0.01 Hz steps. The phase shift between channels I and Q of the DDS is adjustable too, which can be used to compensate for phase errors due to electronic components. The sampling frequency can be set between 180 kHz and 1 MHz with 20 ns steps in period. It should be set high enough in order to make the rejection of the mirror spectrum easy, while enabling the emission of a wide-band signal. However, it remains limited by the throughput of the PC’s bus, which must send data in real-time. The DACs are followed by low-pass filters with a 100 kHz cut-off frequency. Thus the maximum bandwidth of the signal is 200 kHz. Fig. 13.9. Schematic of the Poitiers-Monterfil link. 1152 Clency Perrine, Yvon M. Erhel, Dominique Lemur, Louis Bertel and Alain Bourdillon 13.3.2. The receiving system The system used as a receiver is located in Monterfil. It was at first designed as a direction find- ing system and sounder, and then adapted to the current application. Therefore it is not as flexible as the transmitter. The array is a set of 8 collocated active antennas, as described by Marie et al. (2000) or a hetero- geneous circular array, as depicted in fig. 13.7. The collocated array is followed by 8 identical receivers which transpose the incoming HF signals onto a 30 kHz intermediate frequency. Then they are filtered using a quartz filter with either with a 3 kHz or a 12 kHz bandwidth. The sampling is performed at a programmable rate and the data is stored on the PC disk. The most restrictive features of the receiver are its bandwidth (12 kHz maximum), which might become a concern in future experiments, and the access time of the PC’s hard drive. Indeed, a few seconds of data result in files of several megabytes which must be saved in real-time and too long ac- quisitions often make the acquisition board overflow. 13.3.3. Spatio-temporal equalization As NS sources are incident on the array, the signal at the output of one antenna has the following form: 1 ( ) ( ) ( )x t s t n tFi i k k i NS = +i =k ! ^ h . A source separation method tends to estimate a set of weighting factors wik so that the filtered signal ( )s w x tk ik i NC 1 = !t is the closest estimate of sk(t). In the application, the classical Least Mean Square (LMS) algorithm is processed. The goal is to minimize the power of the error between the estimated signal ( )s tkt and the reference signal sk (Widrow and McColl, 1976) resorting to a descent along the gradient. It is assumed that a training se- quence is transmitted to make the estimation of the weighting factors converge. Then, the effective data transmission begins and their real-time computation involves the estimated received data instead of the training data. 13.3.4. Synopsis of the signal processing performed at reception In the complete receiving system which has been developed, the array processing is only one among a set of functions including: – matched filtering (to the transmitted wave form); – phase and frequency synchronisation; – time synchronisation (optimal placing of the sampling times); – interpolation. Figure 13.10 summarizes all these operations. The input signals xi(t) are the received signals transposed in base band. Their phases are corrected by the phase synchronisation loop. The synchro- nisation techniques are based on an approach of the maximum of likelihood running with the spa- tially filtered signal. Each symbol is over sampled with a ratio of 10: T/ Te =10 if T denotes the sym- bol duration and Te the sampling period. 1153 A way to increase the bit rate in ionospheric radio links A good temporal resolution in the placing of the sample times needs an interpolation: a parabolic interpolator with 4 coefficients has been used. The loop yields an estimate kxt of the time delay at the kth iteration; its recursive expression integrates an error signal generated by a «zero crossing detector» 1+ e= +x ck k kxt t with 1 1- -- /c c kT Te 2k k kk 0= - + x ) ) st t t t_ bi j( 2 . The quantities with a star denote the conjugates of the estimated complex data and the sample of the filtered signal is picked up with a rate twice faster than the symbol rate (Mengali and D’Andrea, 1997). 13.4. RESULTS 13.4.1. Characteristics of the signals received using the 4 collocated antennas The results presented in this section are based on an example that reflects the behavior of the sys- tem in the majority of cases. Only 4 antennas were used out of the 8 available ones because the oth- er 4 proved to be unnecessary: they did not help to improve the transmission performances due a low level of signals for the tested link and, moreover, they induced a larger computational time. The 4 el- ements of the array, as pictured in fig. 13.11, are: – a vertical square loop contained in the East-West plane; Fig. 13.10. Synopsis of the receiver. 1154 Clency Perrine, Yvon M. Erhel, Dominique Lemur, Louis Bertel and Alain Bourdillon – a vertical square loop contained in the North-South plane; – a horizontal square loop; – an XYZ dipole. The transmitted signal has the following properties: – 16-QAM modulation with rolloff factor of 0.2; – 3 kHz bandwidth; – 10 kbps (2500 baud); – carrier frequency of 8 MHz. The received signals are plotted in fig. 13.12. They are all affected by fading, which denotes that propagation is mainly provided by 2 complementary modes of a principal path. In such a case, these modes come from very close directions, however the array’s antenna diversity is enough to introduce a significant difference between the 4 observations. This remark is reinforced by the time-frequency representations of the first two signals (figs. 13.13 and 13.14). Two interference patterns appear (fig. 13.3); in addition to fading (the most visible one), some frequency selectivity shows that a second Fig. 13.13. Time-frequency analysis of the second channel at IF (VNS). Fig. 13.14. Time-frequency analysis of the first channel at IF (VEW). 13.13 13.14 13.11 13.12 Fig. 13.11. Schematic of the collocated array actually used (4 antennas). Fig. 13.12. Signals received by the array (VEW, VNS, H, XYZ). 1155 A way to increase the bit rate in ionospheric radio links path also exists. But for these two interference patterns, the signal fades at different instants and fre- quencies from one antenna to the other, what underlines the decorrelation provided by the diversity. 13.4.2. The heterogeneous array: a way to improve HF transmissions The adaptive method based on the LMS algorithm has been applied successively to 2, 3 and 4 an- tennas in order to evaluate the influence of the array’s size. – 2 antennas: the propagation conditions are too constraining during the first seconds of the trans- mission and the algorithm is unable to keep track. However, it starts to converge again after 6 s and then the BER remains below 10-2 independently of the fading (fig. 13.15). The average BER in this case is 2 × 10-2, not much higher than in the best case with a simple polarization filter. Obviously, two antennas are not enough to realize a reliable filter. – 3 antennas: this time the algorithm converges immediately (fig. 13.16) and the BER always re- mains below 3 × 10-2 with a mean value of 3.7 × 10-3. – 4 antennas: since the previous case led to very good results, the fourth antenna does not bring as much change as the third one (fig. 13.17), but still the BER is improved with an average of 2.8 × 10-3. The SNR obtained in the latter case (after adaptive filtering) is reproduced in fig. 13.18, along with the raw SNR received on the first antenna (the same pattern can be observed on the 3 others with Fig. 13.15. BER obtained using 2 antennas. Fig. 13.16. BER obtained using 3 antennas. 13.15 13.16 Fig. 13.17. BER obtained using 4 antennas. Fig. 13.18. Comparison of SNR on antenna VEW with SNR after filtering. 13.17 13.18 1156 Clency Perrine, Yvon M. Erhel, Dominique Lemur, Louis Bertel and Alain Bourdillon fading occurring at different times). This figure shows the advantage of antenna diversity: whatever the method employed to demodulate the signal, the quality of the transmission can only be improved when the SNR remains almost constant and is not subject to deep fading. 13.4.3. Experimental results with the array of 4 collocated antennas The following experiments are obtained with 4 collocated antennas as shown in fig. 13.11. Sev- eral transmissions have been made at the end of year 2000 with carrier frequencies between 6 and 9 MHz (between Poitiers and Monterfil). Acquisitions were done over 20 s. No coding, interleaving or soft decision was employed. 13.4.3.1. Use of a 3 kHz bandwidth The waveform is a QAM-16 or QAM-64 modulation. Its envelope has a roll-off factor equal to 0.2 and the symbol duration is equal to 0.4 ms (2500 bauds). 13.4.3.1.1. First example: QAM-16 The results presented here have been obtained for a typical case where deep fading appears, due to interferences between two close complementary modes. Figure 13.19 shows the set of 4 signals at the array output. Their minima and maxima do not occur at the same time, which proves that their association can bring diversity other than space diversity. The complex samples which are collected at the output of the spatial filtering (before the hard de- cision) are plotted in fig. 13.20. However, the most energetic symbols are not exactly centered at the corners of the map: that indicates that the power stage of the transmitter is not linear for correspon- ding signals. The corresponding Bit Error Rate (BER) is plotted in fig. 13.21. The mean value is 2.7×10-3 and the data rate is 9700 bps. Fig. 13.19. Four channel acquisition. Fig. 13.20. Constellation at the output of the spatial filter. 13.19 13.20 1157 A way to increase the bit rate in ionospheric radio links 13.4.3.1.2. Second example: QAM-64 In the following experiment, the waveform was a QAM-64 modulation and the bit transfer rate equal to 15 kbps. The average signal to noise ratio (before the spatial filtering) equals 20 dB. The cor- responding bit error rate is plotted in fig. 13.22 and the complex samples collected at the output of the spatial filtering are plotted in fig. 13.23. In this case the mean BER value is 7.2 × 10-3. 13.4.3.2. Use of a 6 kHz bandwidth The waveform is a QAM-16 modulation. Its envelope has a roll-off factor equal to 0.2 and the sym- bol duration is equal to 0.2 ms (5000 bauds) involving a bandwidth of 6 kHz. If the channel parameters are still the same, the SNR is decreased by 3 dB. Figure 13.24 shows the set of 4 signals delivered at the array output. The mean SNR is 18 dB. The variations of the BER is plotted in fig. 13.25 and the mean val- ue is 2.3 × 10-4. The data rate is 19.4 kbps and the corresponding constellation is plotted in fig. 13.26. Fig. 13.21. Variations of BER. Fig. 13.22. Temporal variations of BER. 13.21 13.22 Fig. 13.23. Constellation with a QAM-64. Fig. 13.24. Signals received on the 4 collocated antennas. 13.23 13.24 1158 Clency Perrine, Yvon M. Erhel, Dominique Lemur, Louis Bertel and Alain Bourdillon 13.4.4. Experimental results: heterogeneous circular array of 4 loop antennas The following experiments involved a heterogeneous circular array of 4 loop antennas and the trans- mitted data are images. The goal is to design a still image transmission system with a good reception quality in presence of transmissions errors (Souhard et al., 2003). That work was developed at the IR- COM-SIC laboratory of Poitiers (France) jointly with the IETR in Rennes (France). The aim is to find an efficient transmission system, in terms of visual quality for the received image. Usual source coding schemes are inappropriate when errors occur during the transmission. Many errors correction systems exist but they significantly decrease the compression ratios for a low errors resistance. To solve this problem, a Vector Quantization (VQ) algorithm based on the Kohonen topological maps (SOFM: Self Organization Feature Map) is used. The image used for testing is shown in fig. 13.29. Fig. 13.25. Temporal variations of the BER. Fig. 13.26. Constellation at the output of the spatial filter. 13.25 13.26 Fig. 13.27. Constellation at the output of the spatial filter. Fig. 13.28. State constellation diagram (6.7 MHz). Fig. 13.29. Original transmitted image. 13.27 13.28 13.29 1159 A way to increase the bit rate in ionospheric radio links 13.4.4.1. Results Several acquisitions were made on 20th February 2003 at frequencies between 6 and 9 MHz. For this experiment, the signal bandwidth was limited to 3 kHz. The modulation is a QAM-16 and the symbol’s duration is 0.4 ms. The transmitted files duration varies between 4.9 s and 5.8 s. They are sent continuously over 20 s. The number of acquisition files is 72. For many reasons, some of them are unusable: too weak signal due to bad propagation conditions, jammer or important noise pulse, etc. The number of useable files is 57 and table 13.I presents a list of the performances. The results obtained at the first carrier frequency (6.99 MHz) are rather good: the constellation points are clean (fig. 13.27) and the BER is small. When using the second and third carrier frequencies (7.97 MHz and 8.91 MHz), the results are not flawless: the constellation points are not as clean as in the previous example (fig. 13.28) and the BER is larger. The image in the fig. 13.30d contains twice as much errors than others images. However, the cod- ing gain is 7.19 dB in relation to SPIHT (fig. 13.30c), 8.13 dB in relation to JPEG2000 (fig. 13.30b) coding and 18.14 dB in relation to JPEG (fig. 13.30a). For the same compression rate, the SOFM as- sociated with QAM modulation fits very well for transmission through the ionospheric channel. Table 13.I. Number of data files recorded at 3 frequencies and corresponding BER. Carrier frequency BER interval Number of files 0 ≤ BER ≤ 10-4 21 6.994 MHz 10-4 < BER ≤ 10-3 0 (27 files) 10-3 < BER ≤ 10-2 0 Unusable 6 0 ≤ BER ≤ 10-4 2 7.967 MHz 10-4 < BER ≤ 10-3 16 (22 files) 10-3 < BER ≤ 10-2 3 Unusable 1 0 ≤ BER ≤ 1.10-4 0 8.917 MHz 10-4 < BER ≤ 10-3 5 (8 files) 10-3 < BER ≤ 10-2 3 Unusable 0 Fig. 13.30a-d. Images received using different number of antennas. a b c d 1160 Clency Perrine, Yvon M. Erhel, Dominique Lemur, Louis Bertel and Alain Bourdillon 13.5. CONCLUSIONS AND OUTLOOKS This paper presents a way to achieve a significant improvement of the data transfer rate in a HF radio link by the use of a heterogeneous array. The actual rate of the system reaches up to 20 kbits/s (in 6 kHz bandwidth) whereas the current standard for scalar modems is 4.8 kbits/s in a 3 kHz band- width. Actually, only one link has been tested with a relatively short range and limited number of paths. Moreover, no coding, interleaving or soft decision techniques were employed to reach these performances. The results obtained in terms of BER (about 10-3) are considered to be acceptable for HF modems. With a 4 circular antenna array, a QAM modulation and a robust source coding, the SOFM method fits very well for transmission through the ionospheric channel. Future work will consider longer range links. In this context, the number of received signals will increase in relation to the number of ionospheric paths so that the demodulation processing will have to cope with more severe conditions of reception. In addition, investigations will be performed to adjust the transmitted polarization in order to reduce the number of modes propagating through the ionosphere and this will reinforce the spatial filtering. Another development will be the implementation of a blind algorithm for spatio-temporal equal- ization in the receiving system, replacing the LMS algorithm. The benefit will be to avoid the trans- mission of training sequences and, consequently, to increase the effective data transfer rate. REFERENCES BISIAUX, A. (2001): Définition, conception et réalisation d’un modem vectoriel large bande utilisable dans la gamme Hautes Fréquences, Ph.D. Thesis (Université de Rennes 1). BISIAUX, A. and L. BERTEL (2000): Antenna diversity as a means of improving HF transmissions, in Proceedings of the IEE «HF Radio Systems and Techniques», Guildford, U.K. BUDDEN, K.G. (1985): The Magneto-Ionic Theory and its Applications to the Ionosphere (Cambridge University Press). ERHEL, Y.M., L. BERTEL and F. MARIE (1998): A method of direction finding operating on an array of collocated antennas, in IEEE-AP International Symposium, Atlanta. MARIE, F., Y.M. ERHEL, L. BERTEL and D. LEMUR (2000): Design of a HF compact direction finding system based on collocated antennas, in Proceedings of the IEE, HF Radio Systems and Tech- niques, Guilford, U.K. MENGALI, U. and A. D’ANDREA (1997): Synchronization Techniques for Digital Receivers (Kluwer Academic/Plenum Publishers). SOUHARD, B., C. CHATELIER and C. OLIVIER (2003): On the robustness of joint source/channel coding for transmission through an ionospheric channel, in Proceedings of the 9th International Confer- ence on HF Radio Systems and Techniques, June 2003, Bath, U.K. WIDROW, B. and J.M. MCCOLL (1976): A comparison of adaptative algorithms based on the methods of steepest descent and random search, IEEE Trans. Ant. Prop., 24, 615-637.