Enhanced channel simulator for efficient antenna evaluation

0Associated
Cases 
0Associated
Defendants 
0Accused
Products 
1Forward
Citation 
0
Petitions 
2
Assignments
First Claim
1. A method of simulating a propagation channel, comprising:
 obtaining a channel power angular distribution such that one or more receiving clusters are identified at a receiving side of the propagation channel and one or more transmitting clusters are identified at a transmitting side of the propagation channel to map the power relationship between the receiving clusters and the transmitting clusters;
generating a timevarying channel realization sequence which represents how the response of the propagation channel varies with time by using the relationship between the receiving clusters and the transmitting clusters and one or more actual antenna radiation patterns as measured from one or more antennae;
adding an antenna element isolation effect to said timevarying channel realization sequence such that the response of the propagation channel includes mutual coupling impedance; and
applying baseband abstraction to the time varying channel realization sequence to generate one or more effects of the propagation channel including, but not limited to, one or more of, scattering, diffraction, and lineofsight propagation so that one or more system metrics are generated based on the time varying channel realization sequence.
2 Assignments
0 Petitions
Accused Products
Abstract
Method and apparatus for channel simulation is disclosed. The claimed invention provides method and apparatus 1200 to simulate a propagation channel, particularly a multipleinputmultipleinput (MIMO) channel. The claimed invention further provides a method and apparatus for efficient optimization of antenna by the enhanced channel simulation. The claimed invention takes both antenna characteristics and channel characteristics as inputs, and output timevarying channel realizations to generate the system metrics as the optimization target for antenna under optimization. The claimed invention advantageous provides enhanced channel simulation to meet the accuracy requirement of antenna evaluation.
17 Citations
View as Search Results
SYSTEMS AND METHODS OF TESTING ACTIVE DIGITAL RADIO ANTENNAS  
Patent #
US 20120086612A1
Filed 10/07/2011

Current Assignee
Andrew LLC

Sponsoring Entity
Andrew LLC

Wireless location using multiple mobile station location techniques  
Patent #
US 7,764,231 B1
Filed 09/08/1997

Current Assignee
TracBeam LLC

Sponsoring Entity
TracBeam LLC

Gateway and hybrid solutions for wireless location  
Patent #
US 7,525,484 B2
Filed 01/26/2001

Current Assignee
TracBeam LLC

Sponsoring Entity
TracBeam LLC

RADIATED PERFORMANCE OF A WIRELESS DEVICE  
Patent #
US 20080129615A1
Filed 09/07/2007

Current Assignee
Qualcomm Inc.

Sponsoring Entity


PILOT DESIGN FOR IMPROVED CHANNEL AND INTERFERENCE ESTIMATION  
Patent #
US 20080165891A1
Filed 03/26/2007

Current Assignee
Qualcomm Inc.

Sponsoring Entity
Qualcomm Inc.

Apparatus, and associated method, for testing a mobile terminal in test conditions that emulate an operating environment  
Patent #
US 20070019769A1
Filed 06/30/2004

Current Assignee
Nokia Corporation

Sponsoring Entity
Nokia Corporation

Spectrum sharing communications system for monitoring available spectrum  
Patent #
US 5,655,217 A
Filed 06/05/1995

Current Assignee
SBC Technology Resources Inc.

Sponsoring Entity
SBC Technology Resources Inc.

Smart antenna channel simulator and test system  
Patent #
US 6,236,363 B1
Filed 08/17/1999

Current Assignee
MICRONETICS WIRELESS

Sponsoring Entity
MICRONETICS WIRELESS

Spectrum sharing communications system and system for monitoring available spectrum  
Patent #
US 5,548,809 A
Filed 07/15/1992

Current Assignee
SBC Technology Resources Inc.

Sponsoring Entity
SOUTHWESTERN BELL TECHNOLOGY RESOURCES INC.

Radio channel simulation  
Patent #
US 7,054,781 B2
Filed 05/25/2004

Current Assignee
Keysight Technologies Singapore Holdings Pte. Limited.

Sponsoring Entity
ELEKTROBIT PRODUCTION SOLUTIONS OY

Transmission path simulation method and transmission path simulator  
Patent #
US 20060148429A1
Filed 12/24/2003

Current Assignee
Matsushita Electric Industrial Company Limited

Sponsoring Entity
Matsushita Electric Industrial Company Limited

Systems, methods and apparatus for determining a radiated performance of a wireless device  
Patent #
US 20060194553A1
Filed 10/24/2005

Current Assignee
Qualcomm Inc.

Sponsoring Entity
Qualcomm Inc.

System and method for emulating a multiple input, multiple output transmission channel  
Patent #
US 7,154,959 B2
Filed 08/29/2001

Current Assignee
Intel Corporation

Sponsoring Entity
Intel Corporation

Method and device for simulating radio channel  
Patent #
US 20040088628A1
Filed 08/20/2003

Current Assignee
Keysight Technologies Singapore Holdings Pte. Limited.

Sponsoring Entity
Keysight Technologies Singapore Holdings Pte. Limited.

Broadband monopole/ dipole antenna with parallel inductorresistor load circuits and matching networks  
Patent #
US 20030103011A1
Filed 07/29/2002

Current Assignee
Clemson University

Sponsoring Entity


Wireless location using hybrid techniques  
Patent #
US 20030222820A1
Filed 09/30/2002

Current Assignee
TracBeam LLC

Sponsoring Entity
TracBeam LLC

TRANSMISSION LINE COUPLING AND COMBINING NETWORK FOR HIGH FREQUENCY ANTENNA ARRAY  
Patent #
US 3,824,500 A
Filed 04/19/1973

Current Assignee
Sperry Rand Corporation

Sponsoring Entity
Sperry Rand Corporation

11 Claims
 1. A method of simulating a propagation channel, comprising:
obtaining a channel power angular distribution such that one or more receiving clusters are identified at a receiving side of the propagation channel and one or more transmitting clusters are identified at a transmitting side of the propagation channel to map the power relationship between the receiving clusters and the transmitting clusters; generating a timevarying channel realization sequence which represents how the response of the propagation channel varies with time by using the relationship between the receiving clusters and the transmitting clusters and one or more actual antenna radiation patterns as measured from one or more antennae; adding an antenna element isolation effect to said timevarying channel realization sequence such that the response of the propagation channel includes mutual coupling impedance; and applying baseband abstraction to the time varying channel realization sequence to generate one or more effects of the propagation channel including, but not limited to, one or more of, scattering, diffraction, and lineofsight propagation so that one or more system metrics are generated based on the time varying channel realization sequence.  View Dependent Claims (2, 3, 4, 5)
 6. A method of evaluating and optimizing antenna design, comprising:
defining one or more antenna structure and material parameters as a population of genetic algorithm; generating one or more antenna characteristics by one or more electromagnetic computational methods; computing one or more system metrics by a channel simulator through which one or more receiving clusters are identified at a receiving side of the propagation channel and one or more transmitting clusters are identified at a transmitting side of the propagation channel to map the power relationship between the receiving clusters and the transmitting clusters; choosing a fitness function as the one or more system metrics to generate a new population; and obtaining one or more optimized antenna structure and material parameters.  View Dependent Claims (7, 8, 9, 10, 11)
1 Specification
There are no related applications.
The claimed invention relates generally to a channel simulation method and apparatus, more particularly to a MIMO (MultipleInput and MultipleOutput) channel simulation and a MIMO channel simulator. The claimed invention further relates to an antenna evaluation and optimization method, more particularly to a method of improving antenna designs by MIMO channel simulation.
For antenna design and evaluation according to existing and established methods, four groups of people are generally involved. The first group is the antenna designers, who design antenna module based on the application, carrier frequency, antenna gain and other various requirements. For terminal MIMO antenna design, the characteristics such as antenna radiation pattern and antenna element isolation are the key parameters. The fabricated antenna module is evaluated by another group called the measurement group, who performs a series of field tests so that an antenna is tested in real channel environments. The antenna designers, based on the feedback from the measurement group, improve their designs so that the field test results meet the requirements.
The designed antenna module is delivered to the hardware system design group, who designs specific wireless hardware systems for the given application. For hardware system design, the characteristics such as capacity, BER (biterror rate), and FER (frameerror rate) are the target metrics to optimize. The final wireless hardware system with antenna mounted is traditionally measured and tested in real radio environments by the measurement group. Specific measurement platform is required to perform the measurement on a series of measurement sites. The choices of the sites should cover the target channel environments of the wireless hardware system. On each site, the wireless hardware system is tested with different configurations and the antenna module is tested with different orientations. And the results will be feedback to the system designers as well as the antenna designers to improve wireless system designs and antenna designs.
The fourth group is CAD (Computer Aided Design) tool designers. To enable a fast design cycle, CAD tools, especially certain antenna electromagnetic (EM) simulators, are generally used in the antenna design. The EM simulators generate the antenna characteristics based on antenna structure and antenna materials. The antenna characteristics generated include antenna radiation pattern, antenna transmission characteristics or reflection characteristics such as MIMO antenna element isolation, antenna gain, insertion loss, return loss, and VSWR (voltage standing wave ratio). The antenna designers use the EM simulators to find suitable antenna physical dimensions and antenna materials to achieve satisfactory antenna characteristics.
The claimed invention is on an enhanced channel simulator for antenna evaluation which provides the antenna virtual measurement capability in a simulation environment.
The aforesaid antenna design and evaluation process is advantageous in that the antenna is tested in real channel environments with all possible effects in reality reflected in the results.
However, field measurements generally take long time. The measurement is not repeatable even on the same site because of environment variables. The measurement platform uses a hardware antenna module, so it is not possible to change the antenna parameters on site. Thus if better performance is desired, the antenna needs to be redesigned, fabricated and mounted on the measurement platform before testing continues. The antenna may still have great potential of optimization which may not be completed due to time constraints.
The cost of the field measurement is demanding. Apart from slow evaluation speed and associated costs, antenna performance evaluation is directly limited by the accessibility of the desired testing locations. Iterative cycles can result further compounding time to market as well as underlying design costs.
The major drawbacks of the current antenna design and evaluation process can be summarized as timeconsuming, nonrepeatable, lack of flexibility and high cost.
Consequently, it is desired to have an antenna evaluation method which is efficient, repeatable, flexible and low cost. More critically, new methods should be able to use the actual 3D antenna radiation pattern, generate the timevarying channel realization sequence, and simulate the real radio channel precisely. Thus a preferred embodiment of the new antenna evaluation method is an MIMO channel simulation method. The MIMO channel simulation method ideally uses channel data obtained from either measurement or from channel probing and importing antenna radiation pattern and antenna element isolation, to robustly and accurately simulate the radio channel. Results should be consistent compared with the measurement results.
Given these requirements and constraints, an enhanced MIMO channel simulator is hereby disclosed and claimed. The MIMO channel simulator and its associated methods have clear distinction with the known methods. The previous approaches include the channel modeling methods for baseband simulation, for antenna evaluation and for network planning.
Baseband MIMO channel simulation approaches can be categorized into correlation based channel modeling and doubledirectional channel modeling. With respect to correlation based channel modeling, signal correlation at the transmitter and receiver are characterized as correlation matrices with a simple antenna radiation pattern assumed for the computation of the correlation. Doubledirectional channel modeling extracts all channel effects and models them as scatterers. The radiation pattern is idealistic, and the rays propagating through the scatterer have simple assumptions. For doubledirectional channel modeling, a distinct pair of transmitting and receiving cluster is assumed. These baseband modeling methods generally have simple channel assumptions, thus these method may not be accurately enough to match with the practical measurement results. These baseband modeling methods can be implemented in the channel emulators for hardware system designers to test their wireless hardware systems.
On antenna evaluation, previous known approaches are only limited to utilizing antenna radiation pattern and receiving signal power distribution, and computing antenna correlation, and received power or effective gain. Previous approaches do not generate timevarying channel coefficients and they only use the channel power distribution on receiving side without actually modeling the channel characteristics between transmitter and receiver. Thus these methods are inadequate for channel simulation which characterizes the propagating channel characteristics.
On network planning, raytracing channel modeling has been used. But such modeling often requires very detailed real world floor layout information. Without the ability to rely on robust predictive modeling, raytracing channel modeling cannot be equally scaled to other implementations.
For the claimed simulation method, it differs from previous methods in the following aspects: Compared with baseband channel simulation methods, the disclosed method utilizes actual 3D antenna radiation pattern measured in chamber, and has an enhancement on modeling the properties of clusters which improves its performance in predicting the field measurement results. Additionally, the simulated results can be consistent with the field measurement results. Compared with the antenna evaluation methods using antenna radiation pattern, the claimed method models the physical channel more precisely—including transmitting signal distribution, receiving signal distribution and the relationship of transmitting and receiving clusters. Consequently, the timevarying channel realization sequence is generated to the compute system metrics. The system metrics is finally used to optimize the antenna'"'"'s parameters. Compared with the networking planning methods, the claimed method does not require the detailed real world floor layout information. The timevarying channel realization sequence is the output other than the received signal power in most network planning situations.
In contrast to the available channel modeling methods, the claimed invention provides an enhanced MIMO channel simulation method. The enhanced MIMO channel simulation method simulates the propagation channel between any transmitters and any receivers based upon a more robust and predictive model. As a result, the enhanced channel simulation provides a virtual measurement platform to evaluate the antenna so that the effort of field tests can be reduced. This saves time and makes the evaluation of antenna design parameters more flexible. At the same time, the enhanced MIMO channel simulation still provides accurate evaluation results that are comparable to the field measurements.
The claimed channel simulation method includes obtaining channel power angular distribution, generating timevarying channel realization sequence, adding antenna element isolation effect and generating system metrics. The channel power angular distribution can be generated from the predefined channel scenario, or estimated for channel probing. In generating the timevarying channel realization sequence, the claimed simulation method is enhanced by proposing a cluster mapping matrix defining the power mapping relationship of the transmittingreceiving cluster pairs, which is a more accurate modeling of channel characteristics. The generated channel realization sequence is used to generate the system metrics as capacity, BER and FER indicator.
The channel simulation method is further packaged into a channel simulator with specific inputs and outputs. The inputs of the simulator include the antenna radiation pattern data, antenna element isolation data, and other antenna characteristics if needed. The antenna radiation pattern can be obtained from antenna testing or electromagnetic simulation, and the antenna element isolation can be obtained from the network analyzer measurement. The outputs of the simulator include the timevarying channel realization sequence and the generated system metrics.
The claimed invention also provides a method of antenna evaluation and optimization by using the said enhanced MIMO channel simulator. The system metrics such as capacity, BER and FER indicator, are used for antenna optimization. The optimization process in one aspect of the claimed invention is thus different from the traditional ones which aim at optimizing antenna metrics. The optimization process provided in the claimed invention optimizes the system performance directly. A possible way of optimize the antenna design parameter is to define system metrics as the fitness function in the genetic algorithm optimization process of the antenna parameters. The resulting provides a handy antenna evaluation tool for antenna design under given set of channel scenarios.
Other aspects of the claimed invention are also disclosed.
These and other objects, aspects and embodiments of this invention will be described hereinafter in more details with reference to the following drawings, in which:
Mainly two categories of information are needed by channel simulation. The first category is channel characteristics which in particular are power angular distributions as acquired in the obtaining step 110, and the other category is antenna characteristics which include antenna radiation pattern and antenna element isolation. The antenna radiation pattern is to be input in the channel realization step 120. The antenna element isolation is to be input in the antenna element isolation effect step 130.
Generally there are two ways to obtain the channel characteristics. One way is by generation as shown in
Antenna characteristics are generally provided by either measuring the real antenna module using chamber and network analyzer, or by performing electromagnetic simulation of the antenna design. For terminal MIMO antenna module, the major antenna characteristics are the antenna radiation pattern, and the antenna element isolation of the MIMO antenna.
Secondly, Azimuth parameters of the selected channel scenario are loaded in loading step 220. The Azimuth distributions of indoor and open space in one embodiment are Laplacian distributions, and for outdoor scenarios are Gaussian distributions in another embodiment. The Azimuth parameters include: 1) Number of ray clusters at both transmitting and receiving sides; 2) Angle Spread (AS) at transmitting side and the mean Angle of Departure (AoD) for each cluster; 3) AS at receiving side and the mean Angle of Arrival (AoA) for each cluster; and 4) Power Delay Profile (PDP) for each cluster. The Azimuth Power Angular Distribution (Azimuth PAD) can then be generated from the Azimuth parameters of the selected channel scenario for both transmitting and receiving sides.
Thirdly, Elevation parameters of the selected channel scenario are input in a providing step 230. The Elevation distributions of all scenarios in one embodiment are doubleexponential distributions. The Elevation parameters include: peak angle of doubleexponential distributions, standard deviation of exponential distribution above the peak angle, and standard deviation of exponential distribution under the peak angle. The Elevation Power Angular Distribution (Elevation PAD) is subsequently generated using the Elevation parameters of the selected channel scenario for both transmitting and receiving sides.
Finally, the 3D Power Angular Distribution (3D PAD) is computed in a computing step 240 using the Azimuth PAD and the Elevation PAD for at both transmitting and receiving sides. The Azimuth PAD and the Elevation PAD are independent in a preferred embodiment.
The purpose of channel probing experiments is to find the channel power distribution of the receiving side and transmitting side in the channel measurements. In channel probing, different antennas other than the aforesaid antenna modules for channel measurement experiments are used. These antennas are specificpurpose antennas for channel probing only. Two types of antennas are often used in one embodiment: omnidirectional antennas, and directional antennas. To probe the receiving side, the omnidirectional antenna is put at transmitting side, and the directional antenna is put at receiving side. The highly directional antenna is then rotated for N directions, and the received power values of the directional antenna at the receiving side are recorded.
To probe the transmitting side, the omnidirectional antenna is put at the receiving side, and the directional antenna is put at transmitting side. And the highly directional antenna is rotated for N directions, and the received power values of the directional antenna at the transmitting side are recorded.
Secondly, the estimation algorithm is performed to obtain power angular distribution at both transmitting side and receiving side in estimation step 221. The general computation of the estimation algorithm is shown as follows, which is to be applied to the received power values at the transmitting side as well as the received power values at the receiving side respectively:
The received signal power has N values, and each value is the result of rotation angle 2π/N. In general, the received signal power for either the transmitting side or the receiving side is
where G_{v}(q,f) is the vertical polarized component of the antenna radiation pattern gain at angle (q,f), and G_{h}(q,f) is the horizontal polarized component. P_{v}(q,f) is the vertical polarized component of the channel power at angle (q,f), and P_{h}(q,f) is the horizontal polarized component. For 2D summation approximation, it becomes
Given P_{v,n}=P_{h,n}=P_{n}, G_{v,n}+G_{h,n}=G_{n}, and considering the rotation measurement experiment setup, the channel gains are shifted in one embodiment by one element compared with the previous measurement, the following N linear equations are generated:
R_{1 }to R_{N }are measured by channel probing, and G_{1 }to G_{N }can also be measured in chamber. To estimate P_{1 }to P_{N}, LS (leastsquare) solution can be applied as:
P=(G^{H}G)^{−1}G^{H}R
where P is an N×1 matrix representing the power angular distribution, G is an N×N matrix representing the antenna gain and R is an N×1 matrix representing the received signal power.
The above general computation is also applicable for 3D channel parameter estimation. For wideband system, the above general computation should be solved for each channel tap.
As an embodiment of the claimed invention, MIMO channel simulation is adopted. The principle of the MIMO channel simulation is illustrated in
In general, a doubledirectional channel may have n taps which constitute to rays with n possible delay time. For each tap, with a wireless system setup of P transmitting antennas and Q receiving antennas, the channel realization between the pth transmitting antenna and qth receiving antenna for one tap is computed as:
 where M_{0 }is the number of rays with the assumption that each ray has equal power, thus the ray density represents the power angular distribution. The angle for each ray is (θ_{m},φ_{m});
 R_{T,p}^{(v)}(m) is the antenna complex response for vertical polarized component of the pth transmitting antenna at the direction of (θ_{m},φ_{m});
 R_{R,q}^{(v)}(m) is the antenna complex response for vertical polarized component of the qth receiving antenna at the direction of (θ_{m},φ_{m});
 R_{T,p}^{(h)}(m) is the antenna complex response for horizontal polarized component of the pth transmitting antenna at the direction of (θ_{m},φ_{m});
 R_{R,q}^{(h)}(m) is the antenna complex response for horizontal polarized component of the qth receiving antenna at the direction of (θ_{m},φ_{m});
 Ω_{T,p}(m,t) is the phase shift at transmitting side for the ray at the direction of (θ_{m},φ_{m}), which is Ω_{T,p}(m,t)=exp(jα_{m,T,p})·exp (jβ_{m,T}(t)), where α_{m,T,p }is the phase shift contributed by antenna structure, and β_{m,T}(t) is the phase shift contributed by Doppler. For linear array in one embodiment, the two values are computed as
where D_{max }is the distance of the two antenna elements at both ends, and λ is the wavelength; and
where v is the moving speed and φ_{v }is the moving direction with the assumption of moving in Azimuth plane, t is the time instance of the current channel realization.
 Ω_{R,q}(m,t) is the phase shift at receiving side for each ray, and is defined as Ω_{R,q}(m,t)=exp(jα_{m,R,q})·exp(jβ_{m,R}(t)) representing the phase shift values for the receiving ray of the qth receiving antenna can be obtained similarly as the transmitting side.
 P_{SM }is the polarization scattering matrix which describes the polarization shift of each ray. It is defined as
where η is the portion of power that the polarization is converted, and ψ_{m}^{(x,y) }represents the phase shifts of each ray from x polarization to y polarization. Random phases can be applied on these values.
The resulting channel realization needs to be scaled by the square root of the current tap power. LOS component may need to be added to the first tap, and lognormal fading and path loss need to be added.
Channel simulation is performed on a plurality of channel taps, and a plurality of antenna pairs, and a series of time instances. Therefore, the generated channel realization sequence has the dimension of time instances, delay taps and antenna pairs.
In an embodiment of the claimed invention, channel simulation is enhanced by cluster mapping as shown in
One example of cluster mapping is to map the transmitting rays observed in the second observing direction 492 to the receiving rays observed in the first observing direction 491 despite of different combinations of contributory scatterers.
for M receiving clusters and N transmitting clusters
 where the notation a_{mn }represents the power from m^{th }receiving cluster to n^{th }transmitting cluster. The incorporation of the cluster mapping matrix can simulate the reality closely.
The introduction of cluster mapping matrix is necessary for keeping simulation results for antenna performance consistent with measurement results by taking channel characteristics into considerations. While for baseband simulation, which is not used for evaluating antenna performance, the incorporation of cluster mapping matrix may not be required.
To obtain the element values a_{mn }in the cluster mapping matrix, the following steps illustrate one of the possible embodiments: The first step is to perform channel probing to obtain the receiving power angular distribution and transmitting power angular distribution. From the transmitting and receiving power angular distributions, it is possible to identify a number of receiving clusters and transmitting clusters existing in the propagation channel, and also the range covered by each cluster. The second step is to perform another measurement to identify the power mapping relationship between each receivingtransmitting cluster pairs. A directional antenna is used at transmitting side so that the directional antenna points in a direction to the center of the first transmitting cluster. The halfpower beamwidth and firstnull beamwidth of the directional antenna are chosen to cover the range of the first transmitting cluster only. Channel probing is performed on the receiving side to obtain the corresponding receiving power angular distribution. Different power levels can be observed for different receiving clusters. Based on the power ratio of different receiving clusters, the values of the first column in the cluster mapping matrix, a_{11}, a_{12}, a_{13 }. . . a_{M1}, can be obtained. Similar operations can be carried out to obtain other values in the remaining columns in the cluster mapping matrix.
For the portions of randomly selected rays, channel realizations are generated in channel realizations step 540 by applying transmitting ray propagation processing, polarization shifter and receiving ray propagation processing. The transmitting ray propagation processing (not shown) uses the antenna radiation pattern at transmitter, Doppler at transmitter and transmitting antenna array structure to process the transmitting rays. The polarization shifter (not shown) applies the polarization scattering matrix (PSM), which reflects the polarization change by scattering, reflection and diffraction of the environment. The receiving ray propagation processing (not shown) uses the receiving rays, antenna radiation pattern at receiver, Doppler at receiver and receiving antenna array structure, to process the receiving rays. The iteration 515 through steps 520, 530 and 540 will come to an end if there is enough number of iterations to represent the real channel in checking step 550. The iteration 515 comes to an end in the checking step when all clusters at both transmitter and receiver have been selected even for those cluster pairs with low probability as reflected by their low power mapping relationship in the cluster mapping matrix, and the rays generated has distributed over the whole sphere throughout all possible directions. The specific number of iteration needs to be determined for specific cases given the number of transmitting and receiving clusters, the values in the cluster mapping matrix, and the number of rays each time generated.
All the channel realizations of different ray iterations from different cluster pairs are summed together in summing step 560. The resulting channel realization is subsequently scaled by the power of current path in scaling step 570 and added with effects like LOS (lineofsight) component, path loss and shadowing fading in effect adding step 580.
The embodiment of the claimed invention for MIMO simulation represents the broadest simulation method that covers singleantenna simulations, the claimed channel simulation method can also be applied to nonMIMO cases, such as SISO (singleinput singleoutput), MISO (multipleinput singleoutput), and SIMO (singleinput multipleoutput). To perform the simulation under these cases, the transmitting side and/or the receiving side needs to be set to the single antenna specification.
For transmitting side 810, the voltages V_{S1 }and V_{S2 }are signals transmitted to antenna. The voltages V_{T1 }and V_{T2 }are signals transmitted to the physical channel 820. V_{R1 }and V_{R2 }are the received signals from the physical channel 820 at receiving side 830. And V_{L1 }and V_{L2 }are signals finally received by the wireless system.
The channel realization described in previous embodiments is the channel observed from V_{T1}, V_{T2 }to V_{R1}, V_{R2}. In the following, a method of adding antenna element isolation effect will be described. The antenna element isolation effect is also referred to mutual coupling effect.
The antenna element isolation effect is modeled as a voltage excited by the other element. For example, the additional voltage introduced by the mutual coupling on the first transmitting antenna element 811, V_{12}^{(T)}, is excited by the current on the second transmitting antenna element 812, I_{2}^{(T)}. A mutual coupling impedance Z_{12}^{(T) }is used to describe the currentvoltage relationship. Similar conditions are for the receiving side 830. The channel realization with antenna element isolation effect will subsequently be given by the following equation:
H_{M}=Z_{L}[Z_{L}+Z^{(R)}]^{−1}·H_{R}·Z_{A}^{(T)}[Z_{S}+Z^{(T)}]^{−1 }
Where the load impedance

 the source impedance

 the antenna impedance

 the mutual coupling impedance matrices Z^{(R) }and Z^{(T) }are matrices having the mutual coupling impedance included as matrix elements. On the transmitting side 810,
and on the receiving side 830
channel realization H_{R }represents the channel generated from the doubledirectional channel modeling, and H_{M }is the channel realization with the mutual coupling effect added.
The mutual coupling impedances Z_{12}^{(T)}, Z_{21}^{(T)}, Z_{12}^{(R) }and Z_{21}^{(R) }can be measured by measuring S parameter using network analyzer 850. In
If the antenna elements are symmetric, and Z_{A1}=Z_{A2}=Z_{A}, Z_{12}=Z_{21}=Z_{C}, the relationship will be
Under most circumstances, the symmetric assumption cannot be held true, then the mutual coupling impedances Z_{12 }and Z_{21 }have to be transformed from S parameter using asymmetric antenna model.
The postdetection SINRs (Signal to Interference Noise Ratio) are obtained from the channel realization in step 1020. The spatial filtering matrix is firstly computed, and the postdetection SINR can be computed based on the filtered signal power and noise plus interference power. The symbol error is computed using the approximation functions (Q functions) based on the postdetection SINRs. The symbol error rate is subsequently converted to bit error rate.
In step 1030, the FER indicators are computed based the postdetection SINRs. A possible embodiment of computing the FER indicator is
where β is a parameter related to coding block length and the coding rate. And the range N of the SINRs to be computed depends on the system configuration.
Traditionally antenna designers use the EM simulator to help the design. The electromagnetic computational methods implemented in EM simulators include FDTD (finitedifferential timedomain), FEM (finite element method) and various approaches to solve the Maxwell equations under different conditions. The antenna structure and material parameters as geometrical layout information and material permittivity, are input to the EM simulator, and the antenna characteristics are output. The antenna characteristics generated include the antenna radiation pattern, transmission and reflection characteristics of the antenna. The transmission characteristics include Sparameters S_{21 }and S_{12}, gain, insertion loss and MIMO antenna element isolation. The reflection characteristic includes Sparameters S_{11 }and S_{22}, reflection coefficient, return loss, and VSWR.
One possible embodiment for antenna optimization method can be a genetic algorithmbased optimization. The genetic algorithm (GA) defines a population which is composed of chromosomes. A chromosome is a sequence of bits which determines the key characteristics of the object under study. A fitness function is defined to evaluate the quality of the chromosome. The GA works in the loop of selection, crossover, mutation and acception. In the selection step, two parent chromosomes are selected according to the fitness function. In the crossover step, new offspring is generated with a crossover probability. In the mutation step, the offspring generated in the crossover step is mutated with a mutation probability. In the acception step, the offspring mutated in the mutation step is accepted in the population, and the chromosome with lowest fitness function is deleted. The loop is terminated when an end condition is satisfied.
The key of GA is the fitness function. In current antenna optimizations, the fitness function is generally chosen as minimization or maximization of certain antenna characteristics. Some examples of fitness function include minimizing return loss, minimizing VSWR, or maximizing antenna gain at certain carrier frequency. The population is the antenna structure and material parameters.
In the claimed method, the fitness function is chosen as the system metrics generated by channel simulation. The GA works with the new fitness function to find the optimized antenna design.
In population defining step 1110, antenna structure and material parameters are defined as the population for GA optimization.
In EM computation step 1120, the antenna characteristics is generated by EM computational methods such as FDTD.
In system metric generating step 1130, the system metrics of the antenna characteristics is generated by channel simulator for the application channel conditions. In one embodiment, the system metrics include ergodic capacity, outage capacity, received signal power, postdetection SINR, raw symbol error rate, raw bit error rate, coded bit error rate, and frame error rate.
In obtaining new population step 1140, the new population is obtained by using the system metrics as the fitness function. The selection, crossover, mutation and accepting steps of GA are applied in this step.
The end condition 1150 refers to the variation of the fitness function is under a threshold. Then finally the optimized antenna design parameters are output in step 1160.
The GA is suitable for the computational infeasible problems and adapt the solutions in response to a changing environment. Thus GA is a suitable method for antenna optimization. GA works by discovering, emphasizing and recombining the good schemas (building blocks) of solutions in a highly parallel fashion. For antenna optimization, the good schemas of the antenna for the channel environment is kept, and those not suitable for the channel environment are weakened or destroyed by GA. The application of channel simulator in the antenna parameter optimization can be described by the above GA assisted process.
The antenna evaluation method can be applied to compact terminal antenna isolation and radiation pattern optimization, and outdoor basestation antenna radiation pattern optimization.
The description of preferred embodiments of this invention are not exhaustive and any update or modifications to them are obvious to those skilled in the art, and therefore reference is made to the appending claims for determining the scope of this invention.
The disclosed methods and related system have industrial applicability in evaluation and design of antenna software/hardware, in particular those with MIMO antennas, including MIMO antenna isolation optimization for compact terminal antennas as well as antenna pattern optimization for indoor/outdoor applications.