MIMO Fading Emulator Development with FPGA and Its Application to Performance Evaluation of Mobile Radio Systems

We present four new developments for a multiple-input multiple-output (MIMO) over-the-air measurement system based on our previous studies. The first two developments relate to the channel model for multipath environment generation. One is a further simplification of the circuit configuration without performance degradation by reducing the number of delay generation units, which dominate the performance limit when implementing the circuit on a field-programmable gate array (FPGA). The other is to realize spatial correlation characteristics among the input ports on the transmission side, whereas the previously proposed channel model did not consider this correlation. The third development involves the details of implementing the MIMO fading emulator on an FPGA as a two-stage scheme. The fourth is the demonstration of application examples of the developed system.


Introduction
In recent years, as the speed and capacity of information communication technology have grown, multiple-input multipleoutput (MIMO) transmission schemes have become a key technology for mobile radio systems. Effective evaluation of MIMO communication terminal performance in actual radio wave environments has been obtained in over-the-air (OTA) measurements, in which the device under test (DUT) is placed in a specifically constructed fading environment [1]. Various means of constructing MIMO-OTA measurement systems have been proposed and discussed in the 3rd Generation Partnership Project (3GPP) and other forums for the standardization of mobile communications [2].
Of these MIMO-OTA measurement systems, one type is a field emulation system, which is used to produce a multipath fading environment around the DUT by surrounding the terminal with multiple probe antennas [3][4][5][6][7][8][9]. Conventional OTA measurement systems, as typified in [5], have the benefit of generating an environment based on mature technologies used in commercially available fading simulators. However, the systems proposed so far result in large and complicated configurations for MIMO-OTA applications.
In this context, we have directed our research toward simplification of the circuitry to emulate the fading environment [10,11]. In [10], we proposed a new circuit configuration named "antenna-branch-controlled type." In [11], we used computer simulation of this proposed system to clarify the MIMO channel characteristics as a function of the number of probe antennas. In [12,13], we reported an implementation of the signal-processing unit on a field-programmable gate array (FPGA).
When measuring DUT performance using a MIMO-OTA system, it is necessary to construct a complicated experimental setup in a large radio anechoic chamber with a number of probe antennas. To simplify these measurements, a two-stage method has been proposed [14]. In the two-stage method, the DUT's antenna radiation pattern is measured in the first stage, and the actual performance evaluation is carried out by incorporating the measured antenna pattern data into the second-stage emulation. We used computer simulation to verify the functionality of this two-stage method in a configuration incorporating our proposed antenna-branchcontrolled scheme [15].
In this paper, we present four new developments for a MIMO-OTA measurement system based on our previous studies [10][11][12][13]15]. The first two relate to the channel model for multipath environment generation. One development is a further simplification of the circuit configuration without performance degradation by reducing the number of delay generation units, which dominate the performance limit when implementing the circuit on an FPGA. The other is to realize spatial correlation characteristics among input ports at the transmission side, whereas our previous channel model [11] did not consider this correlation. The third development is implementation of the proposed system on an FPGA incorporating the two-stage scheme, and the fourth is to show application examples of the developed system.

Basic Functions for OTA Measurements.
In OTA testing, a realistic propagation environment is generated around the receiving terminals for measurement of transmission performance characteristics. For OTA system development, the multipath propagation channel models proposed thus far by various researchers [16][17][18][19][20] require the following functions: (i) transmission (TX) side spatial correlation, (ii) Doppler spread, (iii) delay spread, and (iv) reception (RX) side spatial correlation. The angular spread (i and iv), Doppler spread (ii), and delay spread (iii) cause space-, time-, and frequencyselective fading, respectively, forming a so-called "triply selective" multipath environment. Figure 1 shows how these selectivities fit into the overall MIMO-OTA measurement system. Further details on MIMO-OTA systems and their standardization works can be found in [2]. In OTA test methods, the environment for MIMO terminal evaluation is generated by either a reverberation chamber (RC) or a fading emulator (FE) [2,21]. In the RC system, a chamber with metallic walls that effectively reflect radio waves is used to generate a rich multipath propagation environment [22][23][24]. In FE systems, a number of virtual scattering probe antennas encircle DUTs to generate a fading environment [3][4][5][6][7][8][9][25][26][27][28][29]. FE systems can provide flexible control of propagation parameters and various propagation environments, but system construction is complex and costly. RC systems lack flexibility in propagation parameter control but can readily provide a multipath-rich environment over a wide range of frequencies. To incorporate the benefits of both FE and RC, a hybrid OTA testing scheme has been proposed [30]. Although both methods have advantages and limitations, this paper focuses on FE OTA systems. Figure 2 shows a basic configuration of an FE MIMO-OTA measurement system. The overall system comprises three parts: the transmitting antenna ports (M ports), the probe antennas (L antennas), and the receiving antenna ports (N ports) of the terminal (DUT), together with the multipath channel generation part that links the input ports and the probe antennas.

FE OTA Systems.
When designing an OTA system, the number of delayed paths gives the performance limit because they consume circuit resources significantly. Figure 3 shows two typical configurations in the delay generation part. Figure 3(a) shows the conventional configuration, the "path-controlled type," which uses MLK delay units. Almost all existing systems are of this type. Figure 3(b) shows our previously proposed "antenna-branch-controlled type" [11], which uses LK units. Although this configuration is extremely simple, the performance for creating a triply selective multipath environment is basically the same as for the path-controlled configuration. In the latter section in this paper, we propose a further simplified configuration, named "advanced antenna-branch-controlled   type," where the requirement for LK delay units is reduced to MK.

Two-Stage
Method. Construction of the measurement system shown in Figure 2 in a radio anechoic chamber is generally tedious work. Moreover, if the DUT is too large, for example, antennas installed in a vehicle [29], the area required to arrange the probe antenna array around the DUT is particularly large and may exceed the size of a conventional radio anechoic chamber. Even in these cases, it is possible to measure antenna array radiation patterns if the distance required for measurements is available in one direction by emitting radio waves from that direction only while rotating the DUT. These needs have recently led to the proposal of a two-stage method, with measurement of the antenna characteristics in the first stage followed by their input into an FE for performance evaluation in the second stage [14,30,31]. Figure 4 shows the two-stage measurement scheme. The basic configuration is shown in Figure 4(a), while an image of the environment generated in the second stage is shown in Figure 4(b). Since the two-stage scheme does not radiate any radio waves during fading emulation, it would be better to call this the "MIMO fading emulator (MIMO-FE)," rather than an OTA measurement scheme.
As stated above, a large benefit of the two-stage scheme is the simplicity of the system configuration and DUT performance evaluation. The largest difficulty is a requirement for connecting test cables directly to the antenna ports for the antenna pattern measurement during the first stage. Cabled connection to a very small terminal, such as a  smartphone, is not appropriate, since the existence of the cable will interfere with the device antenna pattern. For this purpose, a special data collection scheme without disturbing the true antenna patterns has to be incorporated into the DUT. Although there are some countermeasures for this problem [30,31], it remains generally difficult. Nevertheless, if numerically calculated antenna patterns are available, for example, the performance evaluation can be done very easily by incorporating them into the second-stage emulation.

FPGA-Based MIMO-FE.
When developing OTA fading emulators, it has been common to adopt well-developed fading simulator techniques [5], but this has led to large-scale, high-cost systems. A system configuration using an FPGA is promising for making the fading generation part small and light-weight. Table III of [32] summarizes the performance comparisons of various FPGA-based prototypes. As of 2010, the system proposed in [32] realized NMK < 160 and reached the best performance. The FPGA this system used was an Altera Stratix III, EP3SL150F.
In recent years, a MIMO configuration using generalpurpose software-defined radio transceivers made up of FPGAs has been developed [33,34]. These systems provide a 2 × 2 configuration, and the delay remains within 1 s. On the other hand, we had implemented a MIMO-OTA system with = 4, = 8, = 8, and a maximum delay of 200 s, based on the antenna-branch-controlled method in Figure 3(b) on an FPGA board [13]. The FPGA used was a Xilinx Virtex-6, LX 240T, which is the same FPGA used in Section 4.
When implementing more complicated circuits on an FPGA, devising the circuit configuration is the key for enhancing function and performance against the limited gate array resources.

Background on the Novelty of the Developed System.
The early stage of FE-type MIMO-OTA developments was mainly focused on control of the direction of arrival (DOA) of radio waves around the terminal by using many probe antennas with 7 to 31 elements [3,4,6]. This makes it possible to create Rayleigh fading environments with Doppler frequency variations and RX-side spatial correlation characteristics. However, at this stage, incorporation of delayed wave generation units capable of controlling the delay profile exceeding 1 s was not achieved. On the other hand, in the fading simulator for evaluating the transmission performance of single-input single-output (SISO) systems, a multipath fading link having necessary functions, such as Rayleigh and/or Nakagami-Rice amplitude fluctuations with variable Doppler power spectrum and delay profile, has been available at that time. By applying the technology cultivated in this fading simulator development, an extremely sophisticated FE-type MIMO-OTA measurement system was proposed in [5] based on the path-controlled configuration shown in Figure 3(a) composed of ML links. The developed system based on this scheme has played an important role in standardization of DUT performance measurement schemes [2]. However, since the fading generation function must be created using ML links individually, the system configuration results in large scale and costly. For this reason, more simplification of the system configuration by enhancing the circuit efficiency of multipath generation part in L × M matrix is required from a practical viewpoint.
In addition, the transmission side spatial correlation characteristics, namely, the distribution of the path directions in the base station side as shown in Figure 1, have not been incorporated explicitly in the L × M matrix of developed systems available so far.
Our group, on the other hand, has proposed an antennabranch-controlled method as shown in Figure 3(b) that can realize a MIMO-OTA measurement system with a simple configuration retaining the necessary functions [10,11]. Particularly in [11], we provide detailed configuration and operation principle of the antenna-branch-controlled scheme and evaluate the accuracy of generated propagation environments by computer simulation. The paper [13] deals with a hardware implementation of the method of [11] using an FPGA and realizing the probe antennas of L = 8. On the other hand, [14] proposes a configuration of the antenna-branch-control method with the two-stage method and confirms the validity through computer simulation.
The novelty of this paper based on the above-mentioned previous works is summarized as follows: (i) Further simplification of the configuration of the signal-processing unit, that is, the "advanced antennabranch-controlled method": this scheme can drastically reduce the number of delay units in the L × M connection matrix, which achieves L = 32 probe antennas with 8 delay paths varying up to 200 s (ii) Incorporation of the TX-side spatial correlation property to realize the propagation characteristics around the base station site (iii) Implementation of two-stage MIMO-FE on an FPGA incorporating the functions of (i) and (ii) above (iv) Evaluation and demonstration of effectiveness in a wireless local area network (WLAN) and terrestrial digital broadcasting reception.

Basic Configuration.
We pursued a multipurpose propagation channel model with a triply selective multipath environment where parameters such as direction of arrival, delay profile, and maximum Doppler frequency can be controlled flexibly. In that way, existing MIMO propagation channel models, such as 3GPP SCM, COST 259/273, and Winner [35], can be generated by setting the parameter values as requested.
In the simplified antenna-branch-controlled system [11], the multiple required functions of the multipath channel generator shown in Figure 3(b) are apportioned and connected in a cascade configuration. The spatially uncorrelated feature of the transmission signals is produced using Walsh-Hadamard code sets; then multipath delayed waves are produced for each signal, and finally a different Doppler shift is added for each probe antenna. Rayleigh fading is thus obtained by spatial synthesis of the signals emitted from the probe antennas in a configuration for amplitude and delay control for each path in the connection matrix that is far simpler than conventional path-control methods.
× MIMO transmission and reception signals of transmitter antennas and receiver antennas can be expressed as follows: where s is the transmission signal vector (M dimensions) given as a function of time t, n is the noise vector (N dimensions), r is the reception signal vector (N dimensions), and H FE represents the impulse response matrix ( × dimensions) of channels described as a function of time and delay time . In the equation, the symbol * represents a convolutional integral. The antenna-branch-control method [11] provides a simpler configuration for channel generation by conducting signal processing in branch units of probe antennas that create a multipath environment. H FE is given as follows: where W TX is a TX-side connection matrix with × dimensions consisting of a set of Walsh-Hadamard code vectors, with the matrix producing independent Rayleigh fluctuations against input signals, and A RX is the × channel matrix connecting the probe antenna array with the reception antenna array. For conventional OTA testing systems, A RX is automatically determined by the arrangement of both probe antennas and the DUT receiver array without taking receiver array patterns into account. In the case of two-stage systems, all necessary data for A RX are obtained during the first-stage measurement of the antenna patterns in a radio anechoic chamber. In (2), A Doppler represents an × diagonal matrix with diagonal elements of exp( 2 ) ( = 1, 2, . . . , ) that provide the Doppler shifts corresponding to the angular positions of probe antennas. Each Doppler shift is given by where V is the vehicular velocity and is the wavelength of the radio wave. The small angular offset Δ gives spectral variety [11].

Advanced Antenna-Branch-Controlled Scheme.
When implementing the circuit on an FPGA, increasing the number of delay generation units mainly decreases performance because they consume gate resources significantly. As stated in the previous section, the antenna-branchcontrolled scheme in Figure 3(b) realizes fairly large circuitscale reduction compared with the path-controlled scheme in Figure 3(a). Now, we propose the "advanced antenna-branchcontrolled" type shown in Figure 5, which uses only MK delay units. Considering values of = 4 and = 32, the reduction ratio for the number of necessary units for the new configuration is 1/32 for the scheme in Figure 3(a) and 1/8 for the scheme in Figure 3(b). Here we present the channel model for the proposed configuration.

TX Connection Matrix.
The TX connection matrix W TX with × dimensions provides specific spatial correlation characteristics for the signals from the input ports. Since the determination of this matrix is one of the key developments reported in this paper, a detailed description is provided in Section 3.3. The input signal vector s( ) and the associated output signal vector s ( ) are linked by

Delayed Wave Generation.
To generate delayed signals for each of the delayed waves from the above outputs, we construct an × delayed wave matrix: where ( = 1, 2, . . . , ) is the delay of the kth delay path and is the delta function.

Weight Matrix for Delayed
where ⊗ denotes the Kronecker product, which is used to construct a vector extended to the dimension of KM.
The vector w ( ) delay is the kth element of the weight matrix W delay ( × matrix) composed of Walsh-Hadamard code vectors. Accordingly, the matrix of W delay ⊗ W TX is dimensionally extended to × . The input signal for the respective probe antenna comprises a summation of all the delayed waves. Taking x to be the input vector for the probe antenna array having as its element, we then have where vector b represents the real-valued amplitudes in the angular direction of each probe antenna for the delayed wave .
The proposed configuration is shown in Figure 6, and the overall channel model is represented by the following equation:

TX Connection Matrix Incorporating TX-Side Spatial
Correlation. For the angular profile of path directions in mobile communication systems, a Laplacian distribution with the center in the direction of the mobile terminal can be assumed at the base station side [36]. The angular profile Ω( ) and spatial correlation coefficient (Δ ) can be written as follows: where is the wavenumber of the radio waves, is the arrival angle referenced to the baseline (direction of the array arrangement), 0 is the central direction, and is the standard deviation of the path direction angles. Equation (11) is an approximation conditioned on the angular spread not being overly large [ ≪ 1 (rad)] and is not valid when 0 is approximately equal to 0 [37]. The spatial correlation matrix R TX between ports of the TX-side array antenna that comprises element antennas is represented by the following equation for a linear array equally spaced in intervals of in the direction of = 0: International Journal of Antennas and Propagation  Owing to the eigenvalue properties for the above correlation matrix, the × connection matrix W TX that achieves spatial correlation for the TX-side is given by where TX_ ( = 1, 2, . . . , ) is the mth eigenvalue of the correlation matrix R TX , e TX_ is the eigenvector corresponding to the eigenvalue TX_ , and the superscript T represents the transpose operation. Based on the scheme explained above, a function of the Kronecker model [16] for TX-side spatial correlation is implicitly incorporated.
To create independent fluctuations against all input signals, we recommend using the Walsh-Hadamard code matrix introduced in [11].

Computer Simulation for Channel Model Verification.
A computer simulation was performed to verify the functionality of the proposed model represented by (9) and (13). First, the channel characteristics of the TX-side spatial correlation and the impact on telecommunications performance were examined in the case of a single delay wave for a narrowband Rayleigh fading environment. The specific settings were = = 4, = 1 with amplitude vector b 1 in (8b) for setting all 1 values, and L = 32. If 1 is set to 0, then H FE in (9) can be represented by the following equation using the narrow-band channel response matrix A FE : The time-dependent correlation matrix R FE ( ) of the channel is represented by where superscript H represents the complex conjugate transpose.
The four eigenvalues of this channel correlation matrix are set as 1 ( ) through 4 ( ) in descending order. The element interval for both the TX-side and the RX-side arrays is set to one wavelength. Simulation results using the Kronecker model [16] are used as a reference to see whether the probability distribution of channel eigenvalues for the multipath environment generated by (14) is correctly performed.
The cumulative distributions of the four eigenvalues are shown in Figure 7, where the angular profile of the TX-side paths is a Laplacian distribution with 0 = 90 ∘ and = 5 ∘ . A case of no correlation between transmission ports, that is, the eigenvalues for an i.i.d. environment, is also shown. Comparison of channel capacities for 0 = 90 ∘ , = 5 ∘ and 0 = 30 ∘ , = 5 ∘ are shown in Figure 8 for a signal-to-noise power ratio (SNR) of 10 dB. As shown in Figures 7 and 8, the TX-side spatial correlation imposes large channel capacity degradations compared with the i.i.d.   (1, 1, . . . , 1), b 2 = 0.8b 1 , b 3 = 0.6b 1 , and b 4 = 0.4b 1 . We examined the reception signal ( , , ) ( ), which was set as the port output signal for delay wave from the port input signal. For this signal, we performed the following correlation operation without amplitude normalization: An example of the expected values in the environment established with the above parameters, for = and = , would be 1, 0.64 (0.8 2 ), 0.36 (0.6 2 ), and 0.16 (0.4 2 ) for k = 1, 2, 3, and 4, respectively, and for ̸ = would be 0. In Table 1, a comparison is made for | ( , ) | between these given settings (theoretically expected values) and values with generated signals (simulated values) as "generated/expected." The table shows that the values match up almost perfectly. Moreover, for conditions of ̸ = , where the value | ( , ) | should be 0, all combinations were confirmed to be less than 0.001. Thus, the delay wave environment was clearly generated as expected. From this validation, we can conclude that the multipath delay environment generated by (9) provides a wideband MIMO propagation environment.

System Implementation on an FPGA
In this section, we describe the FPGA implementation of the MIMO-FE incorporating the two-stage functions described in Section 3. Multipath fading for MIMO with any combination of ≤ 4 and ≤ 4 can be created. The FPGA used in this system is a Xilinx Virtex-6, LX240T. The prototype system specifications are given in Table 2. The program language is VHDL, and signals in the FPGA are created in parallel, as shown in Figure 6. The resource consumption ratios in the FPGA are about 90% for block RAM (RAMB36E1), which creates all delay paths, and less than 20% for multipliers (DSP48E1), which create the logic circuits. The interface connecting the FPGA and a personal computer (PC) is USB 2.0, and its protocol is UART (universal asynchronous receiver transmitter). Fading emulation can be carried out in real-time operation, while the parameter setting is done offline before starting the emulation.
The operation clock is synchronized at 160 MHz with the input, output, and signal-processing units. The input and output signals are specified by an intermediate frequency (IF) band of IF = 40 MHz, a signal bandwidth of W = 40 MHz, and an accordingly fractional bandwidth / IF of 1. For distortion-free signal transmission at this bandwidth, we use the digital filter configuration shown in Figure 9(a), in which the Hilbert transformation output ( ) exhibits pass characteristics, as shown in Figure 9 Figure 10 shows an external view of the FE as implemented on the FPGA for all functions shown in Figure 6. The overall system is highly compact, as shown here by its size when mounted in a case measuring 28 × 22 × 5 cm. A functional schematic of the created multipath environment in the FPGA is given in Figure 4(b) for = 32, = 8, and = = 4. For functional verification, we first set the parameters to their most simple terms, namely, = 1, = 1, and = 1, and we inputted an unmodulated signal into input port 1. Figure 11 shows the power spectrum of port 1 output at = 200 Hz and b 1 = 1 when received by an omnidirectional antenna (see Figure 13, Case I). As discussed in [11], in order to exclude the same Doppler shift when setting the probe antennas symmetrically, we adopted small angular offsets to each probe antenna direction. Through the addition of a   Although the average power radiating from each probe antenna is uniform, the averages' spectral density profile becomes approximately U-shaped because the spectral density in the frequency axis becomes higher towards the frequency of ± . This condition corresponds to the Jakes model. Moreover, the system has a function to realize various channel models having an arbitrary angular power profile by adjusting the amplitude of radiated signals from the probe antennas [factors of b in (8b)].
With the four receiving antennas ( = 4) in a cardioid pattern, we rotated the beam direction in increments of 90 ∘ for each antenna and entered the characteristics into the receiving connection matrix A RX (see Figure 13, Case III). Figure 12 shows the received signal spectra of two output ports ( = 1 and = 2). The figure indicates, as shown by the shape of the cardioid pattern, that each of the antennas exhibited strong reception for as many as 16 waves in the beam direction and that the spectra corresponded to antenna direction (the directivity rotation direction).
The generated channel characteristics for input port m, output port n, and delayed wave are denoted by ( ) ( ). We statistically tested amplitude distribution, eigenvalue characteristics, and correlation characteristics and confirmed them as expected. We also found that, for all combinations of delayed waves, ( ) and ( ) were perfectly uncorrelated except for = .
Regarding the combination of and M, there is no particular restriction from the system performance viewpoint in principle. For example, when = 1 and = 4, the effectiveness is demonstrated in the evaluation of maximalratio combining (MRC) diversity (see Section 5.1). Moreover, as shown in Table 1, performance is also good when = 4 and = 1. Therefore, an imbalance between and does not limit the performance of the circuit configuration. However, in FPGA implementation, performance is limited by the gate array resources in the circuit. As the values of and increase, for example, in systems such as massive MIMO, the number of probe antennas must also be increased in order to realize accurate eigenvalue characteristics. With such a scale, in commercially available FPGAs, the number of gates in the circuit becomes insufficient. Therefore, the scale size with arbitrary combination of = = 4 or less is a practical limit to satisfy the specifications in this prototype.

Application to Wireless System Performance Evaluations
In this section, we will demonstrate the effectiveness of our MIMO-FE system with two application examples.

Digital Broadcasting Signal Reception with MRC Diversity in Mobile
Environments. Terrestrial digital broadcasting systems in Japan operate in the ultra-high-frequency band and conform to the Integrated Services Digital Broadcasting-Terrestrial (ISDB-T) standard (bandwidth of 5.6 MHz, orthogonal frequency-division multiplexing [OFDM] modulation using 5,617 subcarriers). In this study, we assessed the reception characteristics of a television receiver (Panasonic TH-L17F1) equipped with a commercial four-antenna diversity function with MRC on an individual subcarrier base. Specifically, we generated two types of multipath environments, a "fast fading" environment, and a "large delay difference" environment using the FE and compared a single antenna and four-antenna array.
We assessed diversity reception performance for the following four cases.
Case II. Four omnidirectional antennas are placed one wavelength from the center in a square array ( = 4).
Case III. Four cardioid-shape directional antennas are placed as in Case II with beams radiating in directions opposed to the center (see Figure 12). Case IV. Four cosine-shape (half-side only) directional antennas are placed as in Case III.
Considering the omnidirectional antenna pattern of I ( ) = II ( ) = 1 for the whole angular range in Cases I and II, we set III for Case III to III ( ) = III (1+cos( −( −1) /2)) and IV for Case IV at IV ( ) = IV cos( − ( − 1) /2), where III = √2/3, while IV = 2 for | | ≤ /2 and IV = 0 for | | > /2. The coefficient is the gain-adjustment factor so that the electric power integral in the horizontal plane is the same as for the omnidirectional antenna. Figure 13 shows the measurement system with the antenna positioning described above. An ISDB-T signal with a center frequency of 40 MHz is input to a port ( = 1) and the output of each antenna was obtained with = 1 for Case I and = 1, 2, 3, and 4 for Cases II, III, and IV. The outputs of each were up-converted from 40 MHz to 557 MHz to convert to the television frequency band. The up-converted signals were connected directly to the ports of the television receiver after disconnecting the existing antennas, and a noise signal was added to the input to adjust the carrier-to-noise power ratio (CNR). For MRC evaluation performance in fast fading conditions, we set = 3 with b 1 = b 2 = b 3 = (1, 1, . . . , 1), 1 = 0, 2 = 1 s, and 2 = 2 s. Under these conditions, reception was investigated while the maximum Doppler frequency was changed. Figure 14 shows the equivalent average CNR values relative to the normalized maximum Doppler frequency OFDM , where OFDM of 1.134 ms is the symbol period, including the guard interval (GI) in the OFDM signal. The results in this figure demonstrate that MRC diversity applications (Cases II, III, and IV) show a large improvement in characteristics. They also show that an array of higher directional antennas (IV > III > II) has better tolerance against Doppler spreading in a multipath environment.
To demonstrate a related application, we evaluated the differential delay tolerance of MRC reception for ISDB-T signals. We set = 4, where b 1 = b 2 = (1, 1, . . . , 1) and b 3 = b 4 = 0.562b 1 (Δ = −5 dB) or 0.316b 1 (Δ = −10 dB) with 1 = 0, 2 = 1 s, 3 = Δ , and 4 = Δ + 1 s. Under these conditions, the delay difference Δ is varied. Figure 15 shows the equivalent average CNR values as a function of delay difference Δ for Cases I and III. Considering that the GI period is 126 s and that the apparent degradation can be seen for Δ exceeding the GI period, the results obtained are reasonable.

Wireless Local Area
Networks. Evaluation of transmission characteristics is a relatively simple matter in cases of oneway transmission, as described above. In cases where a transmission signal is controlled by reception acknowledgments (ACK/NAK), the measurement procedure is more complex because of the need to set up a feedback loop. The following demonstration provides an example assessment procedure for an IEEE 802.11n-compliant WLAN.
The WLAN base unit (i.e., the access point, AP) used in this experiment was a Buffalo WZR-AMPG300NH router, designed to IEEE 802.11a/b/g/n specifications, with two antenna ports for transmission and reception and one for reception only. In our measurement, the 802.11n mode with 40-MHz bandwidth in the 5-GHz band is used for data transmission. Its theoretical maximum speed is 300 Mbps. In this experiment, the original-equipment antennas were replaced with directly connected signal cables. The user equipment (UE) used in this experiment was an Intel Centrino Advanced-N 6200 adapter, which has the same specifications as the AP unit. It was mounted on a board to  Figure 16, included 2 × 2 MIMO communication channels with fading added by the emulator to the downlink (AP→UE). Since the fading addition was performed in the 40-MHz band, down-and up-converters were introduced before and after the emulator, respectively. For upand down-converters, we made circuits using commercially available radiofrequency components, such as mixers, lowpass filters, and circulators with the circuit configuration in Figure 16. The uplink (AP←UE) was for ACK/NAK signal transmission only and was therefore directly connected by a cable without frequency conversion. The transmission and reception signals were separated by a circulator, and unwanted signals were suppressed using attenuators to avoid leakage of the downlink signal via the uplink channel. The transmission signal from the PC was continuously sent to the PC receiving end without interruption, and the throughput was assessed with our evaluation software. Figure 17 shows an example of this evaluation where the throughput in terms of bits per second is depicted as a function of the maximum Doppler frequency and delay difference between two equally powered delayed waves. Only a relative trend for the throughput behavior against the delay and Doppler spread can be presented here, since there is no reference channel for making an absolute comparison. As can be seen from the figure, in Rayleigh fading environments, statistically stable characteristics are obtained within a range of = 100 Hz and a delay difference of about 400 ns. For IEEE 802.11n, GI = 800 ns is the standard, but when the delay spread is smaller, there is also a mode that operates with GI = 400 ns. Therefore, in the vicinity of the delay of 400 ns, the operation tends to become unstable because it is in the threshold position of this switching. A Doppler shift of = 100 Hz corresponds to the normalized maximum Doppler frequency OFDM ( OFDM = 4 s) of approximately 0.0004, which never falls into the category of fast fading but is degraded nonetheless. This may be attributed to the fact that WLAN is designed under the assumption of slow-moving indoor conditions; that is, fast fading tolerance is not considered as a design objective. No further details concerning the WLAN performance are given here, since the objective of the evaluation was to demonstrate the effectiveness of the emulator and not to assess commercial WLAN.

Conclusions
In this study, we achieved implementation of a two-stage MIMO-FE circuit on an FPGA with an extremely simple configuration that is able to generate a variety of fading environments. The input of propagation parameters and receiving antenna characteristics from a PC allows it to produce a 4 × 4 MIMO fading environment in real-time. The results are summarized as follows.
(i) We introduced a more simplified configuration of a multipath fading generation unit, that is, the "advanced antenna-branch-controlled" method. The scheme can drastically reduce the number of delay units in the × connection matrix, which achieves L = 32 probe antennas with 8 delay paths varying up to 200 s.
(ii) In order to realize the propagation characteristics around the base station site, the TX-side spatial correlation property was incorporated, and its validity was confirmed through computer simulation.
(iii) We successfully demonstrated two evaluations of transmission characteristics: application to WLAN and terrestrial digital broadcasting diversity reception.
The next goal is to apply our method to the performance evaluation of next-generation wireless terminals.