CN105842546A - Method of solving inverse matrix of electromagnetic wave propagation matrix based on antenna pattern information - Google Patents

Abstract

The invention discloses a method of solving the inverse matrix of an electromagnetic wave propagation matrix based on antenna pattern information of a wireless terminal. A microwave anechoic chamber comprises n test antennas inside, a wireless terminal comprises m receiving antennas, and the m receiving antennas and selected m test antennas carry out data transmission in the form of electromagnetic wave, wherein m>=2, and n>=m. The method of solving the inverse matrix of an electromagnetic wave propagation matrix comprises the following steps: getting the pattern information of the m receiving antennas; calculating out coordinate points suitable for solving the inverse matrix of an electromagnetic wave propagation matrix according to the pattern information; and solving the inverse matrix of a corresponding propagation matrix at the suitable coordinate points. Through the solving method, an inverse matrix meeting requirements can be solved quickly during calculation based on a radiation two-step approach (one of the international standards for MIMO test), no human intervention is needed, and time and effort are saved.

Description

Electromagnetic wave propagation inverse of a matrix matrix solving method based on antenna radiation pattern information

Technical field

The invention belongs to wireless device performance technical field of measurement and test, particularly relate to a kind of based on wireless terminal antenna pattern information The method for solving of electromagnetic wave propagation inverse of a matrix matrix, and the method tests state at LTE MIMO OTA (Over The Air) One of border standard radiation two-step method is applied.

Technical background

The purpose of MIMO OTA test is to ensure that the test result at laboratory can truly reflect that wireless terminal is in various complexity Wireless performance under practical service environment and user's use state.

At present, international wireless highest standard tissue, 3GPP (3rd Generation Partnership Project) provides Main testing standard has Multiple probe technique, radiation two-step method, reverberation method.Wherein, in the related, disclose a kind of based on Radiation two-step method, quick, accurate, economic MIMO OTA test Solution, as it is shown in figure 1, terminal test environment Include base station emulator, channel simulator, shielding anechoic chamber, turntable, antenna, it is possible to realize completing in efficient Measurement to the performance of LTE terminal.

Radiation two-step method is that conduction two-step method develops.The realization of conduction two-step method is by calculate in a computer should Arrive the signal of each receiver, be then delivered directly on this receiver by call wire, it is achieved that and Multi probe equivalence Measurement process.But, the OTA that disturbs mutually and do not meet that can not comprise reception antenna introducing owing to conducting two-step method surveys Examination theory, the most progressively to radiation two-step method development.Radiation two-step method is on the basis of conduction two-step method, removes call wire, Allow the free-space propagation that signal is simulated in darkroom to measured piece, then utilize corresponding method to demodulate this spatial transmission matrix Inverse matrix, so, it is possible to realize the effect of two virtual call wires in space, make signal on virtual call wire It is transmitted directly to reach receiver.Therefore, one of the difficulty that radiation two-step method the faces space matrix to electromagnetic wave propagation Solving of inverse matrix.

The space matrix of electromagnetic wave propagation is not only relevant with dark room conditions, and with the aerial position of selected transmitting signal, Polarizing and relevant with the attitude of measured piece, the combination of these variablees can produce very many electromagnetic wave propagation matrixes, due to Whether one matrix exists inverse matrix is conditional on, therefore, in these propogator matrixs, it is understood that there may be much can not carry out Inverse matrix solves, if which propogator matrix cannot be predicted be suitable for solving inverse matrix, unique way be exactly in order wheel sequence or Whether each possible forward-propagating inverse of a matrix matrix of person's random challenge meets requirement.When such way is extremely to expend Between and poor efficiency, often occur all not having training in rotation to arrive satisfied inverse matrix even more for a long time several hours.It is the most fatal, Tester must test by all-the-way tracking, and the performance of each solution inverse matrix out is made the most satisfied subjective requirement Judge.Therefore, development and the popularization of radiation two-step method are seriously hampered.

Summary of the invention

It is contemplated that one of technical problem solved the most to a certain extent in correlation technique.To this end, the present invention needs to carry Going out the method for solving of a kind of electromagnetic wave propagation matrix inversion matrix based on wireless terminal antenna pattern information, these computational methods exist Carry out radiating inverse matrix in two-step method and the inverse matrix meeting requirement when calculating, can be solved rapidly, it is not necessary to human intervention, save time Laborsaving.

In order to solve the problems referred to above, the embodiment of the present invention proposes a kind of electromagnetic wave propagation matrix inversion based on antenna radiation pattern information The method for solving of matrix, wherein, includes turntable and n test antenna in microwave dark room, wherein, wireless terminal is placed on institute Stating on turntable, described wireless terminal includes m reception antenna, described m reception antenna and m selected test antenna with Electromagnetic wave carries out data transmission, and wherein, m >=2, n >=m, described method for solving comprises the following steps: obtain described m The pattern information of individual reception antenna；Pattern information according to described m reception antenna calculates and is suitable for solving electromagnetic wave propagation The coordinate points of inverse of a matrix matrix, wherein, described coordinate points includes the parameter information affecting described electromagnetic wave propagation matrix；With And in the described coordinate points being suitable for solving electromagnetic wave propagation inverse of a matrix matrix, to corresponding electromagnetic wave propagation inverse of a matrix matrix Solve.

The method for solving of the electromagnetic wave propagation matrix inversion matrix based on mobile terminal antenna pattern information according to the present invention, can The most satisfied condition properly solving inverse matrix of each coordinate is judged, more with the pattern information of direct basis reception antenna Calculate rapidly and be suitable for solving the coordinate points of electromagnetic wave propagation inverse of a matrix matrix, in applicable coordinate points, inverse matrix is solved, And can intelligent realize, it is not necessary to human intervention, time saving and energy saving.

Further, it is suitable for solving electromagnetic wave propagation inverse of a matrix square according to the pattern information calculating of described m reception antenna The coordinate points of battle array specifically includes: extract correspondenceEach reception antenna of coordinate points relative to Gain information and any two reception antenna of each test antenna test the phase place receiving electromagnetic wave of antenna relative to each Difference information, wherein, (Ant₁,Ant₂,...,Ant_m) it is the information of the selected test antenna of described coordinate points,For described The rotation angle information of the turntable of wireless terminal is placed in microwave dark room；According to described gain information and described phase information, The m bar void being loaded with between m reception antenna and m the test antenna realized after corresponding inverse matrix is obtained in computing respectively Intend the isolation of wire；Isolation according to described m bar virtual wires judges corresponding described Whether coordinate points is the qualified coordinate points solving described inverse matrix.

Further, judge that correspondence is described according to described m bar isolationCoordinate points is No is that the qualified coordinate points solving described inverse matrix specifically includes: if corresponding describedSit The sign symbol of the isolation of the m bar virtual wires of punctuate is identical and absolute value is all higher than predetermined threshold value, then judge this coordinate points For solving the qualified coordinate points of described inverse matrix.

Said method also includes: traversalAll coordinate points of affiliated scope, it is thus achieved that solve The all of qualified coordinate points of described inverse matrix；And determine in described all of qualified coordinate points and described most preferably to solve electromagnetism The coordinate points of the inverse matrix of ripple propogator matrix.

Further, described all of qualified coordinate points determines the described electromagnetic wave propagation inverse of a matrix matrix that most preferably solves Coordinate points specifically includes: take absolute value the isolation of the m bar virtual wires of coordinate points qualified each described, and will be definitely Isolation angle value minimum in value comprehensively isolates angle value as respective coordinates point；Select in described all of qualified coordinate points comprehensive The coordinate points of isolation maximum is the coordinate points most preferably solving electromagnetic wave propagation inverse of a matrix matrix.

Or, calculate the virtual path loss of described m bar virtual wires respectively；And the m bar to coordinate points qualified each described The virtual path loss of virtual wires takes absolute value, and using comprehensive empty as respective coordinates point of virtual circuit loss value maximum in absolute value Intend circuit loss value, select integrated virtual circuit loss value in described all of qualified coordinate points minimum coordinate points most preferably ask as described Solve the coordinate points of electromagnetic wave propagation inverse of a matrix matrix.

Or, the isolation of the m bar virtual wires of coordinate points qualified each described is taken absolute value, and by absolute value Little isolation angle value comprehensively isolates angle value as respective coordinates point；Calculate the virtual path loss of described m bar virtual wires respectively, The virtual path loss of the m bar virtual wires of coordinate points qualified each described is taken absolute value, and by maximum virtual in absolute value Circuit loss value is as the integrated virtual circuit loss value of respective coordinates point；Comprehensive isolation according to described all of qualified coordinate points is with comprehensive Close virtual circuit loss value and determine the described coordinate points most preferably solving electromagnetic wave propagation matrix inversion matrix.

Accompanying drawing explanation

Fig. 1 is the part-structure schematic diagram of wireless terminal test environment according to an embodiment of the invention；

Fig. 2 be the conduction two-step method in correlation technique realize schematic diagram；

Fig. 3 is to radiate two-step method in correlation technique to realize system schematic；

Fig. 4 is that in radiation two-step method according to an embodiment of the invention, spacing wave propagates schematic diagram；

Fig. 5 be in two receiver practical work process according to an embodiment of the invention signal propagate schematic diagram and etc. Effect schematic diagram；

Fig. 6 is based on wireless terminal antenna pattern information according to one embodiment of present invention electromagnetic wave propagation matrix inversion square The flow chart of the method for solving of battle array；

Fig. 7 is the realization configuration system schematic for method for solving described in Fig. 6 of a specific embodiment according to the present invention；

(1) and (2) in Fig. 8 is the directional diagram letter of two reception antennas of another specific embodiment according to the present invention The data form of breath.

Detailed description of the invention

Embodiments of the invention are described below in detail, and the example of described embodiment is shown in the drawings, the most identical Or similar label represents same or similar element or has the element of same or like function.Retouch below with reference to accompanying drawing The embodiment stated is exemplary, it is intended to is used for explaining the present invention, and is not considered as limiting the invention.

The inverse matrix asked for state the computational methods of the embodiment of the present invention can use, first in radiation two-step method The concept of the several definitionals related in radiation two-step method and the embodiment of the present invention is done simple introduction.

As a example by two receivers, the realization of conduction two-step method can be as in figure 2 it is shown, by channel model and the side of antenna Run in instrument to the computing of figure information, then operation result is directly inputted into by call wire the reception of measured piece Machine, the effect i.e. reached can be expressed as:

[R₁ R₂]=[T₁ T₂] (1)

Wherein, R₁And R₂Represent the signal that two receivers receive, T respectively₁And T₂It is expressed as channel simulator Two-way output signal.

As it is shown on figure 3, for radiation two-step method, propagate in darkroom due to signal and be equivalent at free-space propagation, It is assumed that spatial transmission matrix is:

$H=\left[\begin{array}{cc}{h}_{11}& h_{12}\\ h_{21}& h_{22}\end{array}\right]---\left(2\right)$

Then launch signal and reception signal relation be expressed as:

$\left[\begin{array}{cc}{R}_1^{,}& R_2^{,}\end{array}\right]=\left[\begin{array}{cc}{T}_1& T_2\end{array}\right]*\left[\begin{array}{cc}{h}_{11}& h_{12}\\ h_{21}& h_{22}\end{array}\right]---\left(3\right)$

If reaching the effect as conduction two-step method, before can launching launching signal, enter in instrument The following computing of row, is multiplied by a new matrix:

$C=\left[\begin{array}{cc}{c}_{11}& c_{12}\\ c_{21}& c_{22}\end{array}\right]---\left(4\right)$

The newest relation launching signal and reception signal is:

$\left[\begin{array}{cc}{R}_1^{,,}& R_2^{,,}\end{array}\right]=\left[\begin{array}{cc}{T}_1& T_2\end{array}\right]*\left[\begin{array}{cc}{h}_{11}& h_{12}\\ h_{21}& h_{22}\end{array}\right]*\left[\begin{array}{cc}{c}_{11}& c_{12}\\ c_{21}& c_{22}\end{array}\right]---\left(5\right)$

Then when meeting condition:

$C*H=\left[\begin{array}{cc}x& 0\\ 0& x\end{array}\right]---\left(6\right)$

Wherein, x is not zero, it is possible to realizes radiation two-step method and the equivalence of conduction two-step method, i.e. reaches formula (1) Effect, i.e. equivalence can think space exists two " wire ", be referred to as the most in an embodiment of the present invention " virtual wires ", can be directly inputted to receiver by signal, as shown in Figure 3 by these two virtual wires.With Upper formulae express is for equally applicable during M (M >=3) individual receiver.

Forward electromagnetic wave propagation matrix and selected test antenna, the measured piece anglec of rotation propagated in darkroom are relevant.Connect with 2 As a example by 22 forward matrixes of receipts machine, define coordinateWherein, Ant₁First survey selected by expression Examination antenna, Ant₂Second test antenna selected by expression, θ represents the angle that measured piece rotates on θ axle in such as Fig. 1 Degree,Represent the angle that measured piece rotates on φ axle in such as Fig. 1, then the corresponding forward electromagnetic wave of each coordinate passes Broadcast matrix.

But, there is inverse matrix or be not each in the forward electromagnetic wave propagation matrix that not each coordinate is corresponding Inverse matrix can meet actual demand.It is directed to this, proposes the finger of " isolation " in an embodiment of the present invention Mark.In case of 2 receivers, after adding inverse matrix, real work situation such as Fig. 5 of 2 receivers (1) shown in, if inverse matrix C, C writes in instrument, is expressed as above-mentioned formula (4) form, and in darkroom, space passes Broadcast matrix table to be shown as H, H and be expressed as:

$H=\left[\begin{array}{cc}{h}_{11}& h_{12}\\ h_{21}& h_{22}\end{array}\right],$ Wherein, when meeting H*C equal to formula (6), in Fig. 5 (1), situation can wait Imitate in as shown in Fig. 5 (2).Actually, receiver 1 (R1) not only can receive the signal that transmitter 1 (T1) sends, If above-mentioned watt level is expressed as P₁₁, unit is dBm, and receiver 1 (R1) also can receive the interference of transmitter 2 (T2), If the above-mentioned watt level that sets is as P₁₂, then it is assumed that the jamproof isolation size that virtual wires 1 is capable of is:

Iso₁=P₁₁-P₁₂ (7)

Similarly, the signal making the transmitter 2 (T2) that receiver 2 (R2) receives is P₂₂, receiver 2 (R2) receives The interference of transmitter 1 (T1) be P₂₁, then it is assumed that the jamproof isolation size that virtual wires 2 is capable of is:

Iso₂=P₂₂-P₂₁ (8)

When two isolation sizes meet certain requirements, it is believed that from the interference signal on another road to receiver work The impact made can be ignored, and i.e. thinks the test result as obtaining in the case with conduction measurement.Need explanation , although isolation is not the sole indicator of virtual wires, but, isolation is to solve for whether virtual wires meets The necessary index that radiation two-step method requires.In an embodiment of the present invention, isolation is stated with dB data mode, but, Data form does not do requirement in actual applications, and real number is expressed can also.

On the basis of virtual wires meets insulated degree requirement, also define " virtual path loss " in an embodiment of the present invention Concept.Specifically, as shown in Fig. 5 (2), it is assumed that the watt level that T1 and T2 two paths of signals sends is respectively P_T1 And P_T2, all representing with dBm data form, the watt level that two receivers of R1 and R2 receive is respectively P_R1And P_R2, Then the virtual path loss of definition virtual wires 1 and virtual wires 2 is respectively as follows:

L_p1=P_T1-P_R1 (9)

L_p2=P_T2-P_R2 (10)

By above-mentioned two formula it can be seen that the meaning of virtual path loss can assess the virtual wires loss to power.

In the above-described embodiments, coordinate pointsNumber the most.But, owing to asking in advance Solve the isolation of the virtual wires of its correspondence after the inverse matrix of each coordinate correspondence propogator matrix, just cannot judge this Coordinate points is if appropriate for solving inverse matrix.The selection of coordinate is mainly taked traversal or random manner by radiation two-step method, I.e. travel through all of coordinate by physical location (and radio frequency switching), solve inverse matrix the most successively and verify corresponding The isolation of virtual wires, satisfactory i.e. stop until finding.This mode very elapsed time and can not protect Card find so that virtual wires isolation reach maximum coordinate, thus greatly reduce radiation two-step method test efficiency with And capacity of resisting disturbance.

For the problems referred to above, the embodiment of the present invention proposes a kind of new inverse matrix computational methods, can quickly, intelligently Find satisfactory coordinate points, and then the forward-propagating inverse of a matrix Matrix Solving corresponding to this coordinate points.

The side of solving to the electromagnetic wave propagation matrix inversion matrix of the wireless terminal antenna pattern information of the embodiment of the present invention below Method illustrates.Wherein, data transmission is carried out in microwave dark room, includes turntable and n test days in microwave dark room Line, wireless terminal is placed on turntable, and described wireless terminal includes m reception antenna, m >=2, the n in microwave dark room Selected m test antenna, n >=m in individual test antenna, described m reception antenna and m test antenna are with electromagnetic wave Carrying out data transmission, i.e. m reception antenna receives electromagnetic wave, m test days beta radiation electromagnetic wave, thus realizes data and pass Defeated, the most now in microwave dark room, amplitude and the phase place of electromagnetic wave propagation between m reception antenna and m test antenna become Change can represent with the math matrix of a m × m, and the inverse matrix of this math matrix is i.e. entered by the method for solving of the embodiment of the present invention Row solves.

As shown in Figure 6, the inverse matrix method for solving of the embodiment of the present invention comprises the following steps:

S1, obtains the pattern information of m reception antenna.

Wherein, antenna radiation pattern is one of performance of antenna.Specifically, as it is shown in figure 1, tested wireless terminal (DUT) is put Putting the center at a turntable, the distance between test antenna and tested wireless terminal meets standard regulation.Rotating and radio terminal When, test antenna measurement obtains tested wireless terminal at space all directions antenna gain and any two reception antenna Phase contrast relative to each test antenna.The measurement of this step can be obtained by the return of Radio Terminal Function such as mobile phone.

S2, calculates according to the pattern information of m reception antenna and is suitable for solving the coordinate points of electromagnetic wave propagation inverse of a matrix matrix.

Wherein, coordinate points includes the parameter information affecting electromagnetic wave propagation matrix, i.e. includes the test antenna selected and wireless end The position etc. of end can change all information of propogator matrix, such as, the positional information of m test antenna and the rotation of wireless terminal Gyration information.Specifically, definition coordinateWherein, (Ant₁,Ant₂,...,Ant_m) For testing the positional information of antenna,For the rotation angle information of wireless terminal, m is the nothing carrying out MIMO OTA test The number of the reception antenna of line terminal, Ant₁,Ant₂...,Ant_mRepresent the test antenna in m the darkroom chosen, darkroom Nei Bao Including at least m test antenna, wherein, m >=2, θ is the wireless terminal as shown in Figure 1 anglec of rotation on θ axle, turns Dynamic step size, depending on slewing area is according to darkroom turret design situation.For wireless terminal as shown in Figure 1 at φ axle On the anglec of rotation, the step size of rotation, depending on slewing area relates to situation according to darkroom turntable.

Specifically, correspondence is extractedThe each reception antenna of coordinate points is relative to each test days The phase information i.e. any two reception antenna phase receiving electromagnetic wave between gain information and any two reception antenna of line Each test antenna is received to the phase contrast of signal.Such as, extracting wireless terminal at turntable isDuring angle, phase For the antenna radiation pattern information of each test antenna, test antenna Ant relative to any one_i(0 < i≤m), extraction The pattern information of the reception antenna of wireless terminal includes yield value and the phase place of any two reception antenna of m reception antenna Difference.As a example by Fig. 4, wherein, T₁,T₂,...,T_mFor selected m different test antenna, R₁,R₂,...,R_mFor m Reception antenna, can be expressed as the amplitude extracted and phase information:

Wherein, α_xyRepresent gain (real number representation),Represent phase contrast, more than a kind of representation, and do not limit In this unique representation, and the representation of gain and any two reception antenna phase contrast is not done requirement herein.

Further, according to gain information and the phase information of said extracted, computing is obtained and is loaded with corresponding inverse square respectively The isolation of the m bar virtual wires between m the reception antenna and m the test antenna that realize after Zhen.

The computing formula of the isolation of above-mentioned m bar virtual wires relevant with the value of m (its definition may refer to formula (7) and Formula (8)), during m difference, the isolation computing formula of virtual wires is different, the application focus on when providing m=2 ( Method for solving for 2 × 2 propogator matrix inverse matrixs) the isolation solution formula of 2 virtual wires as the most hereafter.

As m=2, virtual wires isolation solves explanation:

Specifically, correspondence is extractedTwo reception antennas of coordinate points are believed relative to the gain of each test antenna The phase information receiving electromagnetic wave between breath and above-mentioned two reception antenna.Such as, extracting wireless terminal at turntable isDuring angle, relative to the antenna radiation pattern information of each test antenna, test antenna relative to any one Ant_i(0 < i≤m=2), the pattern information of the reception antenna of the wireless terminal of extraction includes the gain of 2 reception antennas Value and the phase information of two reception antennas.Wherein it is possible to the amplitude extracted and phase information are expressed as:

Wherein, α_xyRepresent gain (real number representation),Represent phase contrast.The void then defined by formula (7) and (8) The specific formula for calculation of the size intending wire isolation is as follows: the jamproof isolation size that virtual wires 1 is capable of For (isolation represents with dB digital form):

${\mathrm{Iso}}_1=10*\lg (\frac{{{\mathrm{\&alpha;}}_{11}}^2}{{{\mathrm{\&alpha;}}_{21}}^2}+\frac{{{\mathrm{\&alpha;}}_{12}}^2}{{{\mathrm{\&alpha;}}_{22}}^2}-\frac{2*{\mathrm{\&alpha;}}_{11}*{\mathrm{\&alpha;}}_{12}*cos({\mathrm{\&delta;}}_{21}-{\mathrm{\&delta;}}_{22})}{{\mathrm{\&alpha;}}_{21}*{\mathrm{\&alpha;}}_{22}})+20.---\left(13\right)$

The jamproof isolation size that virtual wires 2 is capable of is (isolation represents with dB digital form):

${\mathrm{Iso}}_2=10*\lg (\frac{{{\mathrm{\&alpha;}}_{22}}^2}{{{\mathrm{\&alpha;}}_{12}}^2}+\frac{{{\mathrm{\&alpha;}}_{21}}^2}{{{\mathrm{\&alpha;}}_{11}}^2}-\frac{2*{\mathrm{\&alpha;}}_{21}*{\mathrm{\&alpha;}}_{22}*\cos ({\mathrm{\&delta;}}_{21}-{\mathrm{\&delta;}}_{22})}{{\mathrm{\&alpha;}}_{11}*{\mathrm{\&alpha;}}_{12}})+20---\left(14\right)$

Wherein, α_xyFor gain information, δ_xyFor phase information, x=1,2, y=1,2.

Further, correspondence is judged according to the isolation of m bar virtual wiresCoordinate points is The no qualified coordinate points being to solve inverse matrix.Specifically, if correspondingThe m of coordinate points The sign symbol of the isolation of bar virtual wires is identical and absolute value is all higher than predetermined threshold value such as 6dB, then judge this coordinate points For solving the qualified coordinate points of inverse matrix, such as this coordinate points of labelling is " qualified " coordinate points.

It should be noted that need traversalAll coordinate points of affiliated scope, it is thus achieved that ask Solve all of qualified coordinate points of inverse matrix, and then determine in all of qualified coordinate points and most preferably solve electromagnetic wave propagation matrix The coordinate points of inverse matrix, i.e. find optimum to solve inverse matrix.

Criterion herein for optimum is not particularly limited, and in one embodiment of the invention, sits each Punctuate defines the concept of comprehensive isolation, and its value is equal to the absolute value of the isolation of all virtual wires in a coordinate points Minima.Specifically, the isolation of the m bar virtual wires of each qualified coordinate points is taken absolute value, and by absolute value The isolation angle value of middle minimum comprehensively isolates angle value as respective coordinates point, and then can select to combine in all of qualified coordinate points The coordinate points closing isolation maximum is the coordinate points being best suitable for solving electromagnetic wave propagation inverse of a matrix matrix, i.e. solves inverse matrix Optimum.

It addition, each coordinate points defines the concept of integrated virtual path loss, its value is led equal to all virtual in coordinate points The maximum of the absolute value of the virtual path loss of line.Specifically, the virtual path loss of m bar virtual wires is calculated respectively, for virtual The calculating of path loss is referred to formula (9) and (10), here repeats no more, and then, to each qualified coordinate points The virtual path loss of m bar virtual wires takes absolute value, and using virtual circuit loss value maximum in absolute value combining as respective coordinates point Close virtual circuit loss value, select the coordinate points that in all of qualified coordinate points, integrated virtual path loss is minimum as being best suitable for solving electromagnetism The coordinate points of the inverse matrix of ripple propogator matrix.

Or, consider comprehensive isolation and integrated virtual path loss to determine best coordinates, specifically, qualified to each The isolation of the m bar virtual wires of coordinate points takes absolute value, and using isolation angle value minimum in absolute value as respective coordinates point Comprehensively isolate angle value, and, calculate the virtual path loss of m bar virtual wires respectively, to coordinate points qualified each described The virtual path loss of m bar virtual wires takes absolute value, and using virtual circuit loss value maximum in absolute value combining as respective coordinates point Close virtual circuit loss value, and then determine according to isolation comprehensive in all of qualified coordinate points and integrated virtual path loss and most preferably solve electricity The coordinate points of electromagnetic wave propagation inverse of a matrix matrix, such as using both integrated value as criterion to determine optimum coordinate.

S3, solves in the inverse matrix being suitable for solving the coordinate points correspondence propogator matrix of electromagnetic wave propagation inverse of a matrix matrix.

Putting it briefly, wireless terminal is in microwave dark room or free space, and antenna number is tested in its reception antenna and darkroom During according to transmission, the amplitude of electromagnetic transmission and phase place change can be indicated with one group of propogator matrix, this propogator matrix and survey The selection of examination antenna is relevant relative to the position of test antenna with wireless terminal.The method for solving of the embodiment of the present invention, based on nothing The pattern information of the reception antenna of line terminal, can calculate all propogator matrixs quickly, intelligently and be multiplied it with inverse matrix The isolation of the virtual wires of rear formation and virtual path loss, and then according to the isolation size of virtual wires, can cook up suitable Close test antenna and the combination of the wireless terminal anglec of rotation doing inverse matrix.

As a example by two the two MIMO OTA performance tests received, the inverse matrix method for solving of the embodiment of the present invention is said below Bright.

In actual applications, as it is shown in fig. 7, draw according to the configuration of system as shown in the figure, system shown in Figure 7 configuration is said Bright as follows: test instrumentation is Keysight E7515A (Chinese referred to as UXM), and it is integrated with anolog base station and mimo channel Emulating instrument.Microwave dark room includes: casing, suction ripple are cotton, test antenna, turntable, radio-frequency (RF) switch, controller etc..Wherein, Test antenna is positioned in Fig. 7 at 1 to the position, position 11 of mark.There are two cross-polarized test antennas each position, Communication antenna is positioned at casing tip position 0, and turntable is only along an axis rotation (Z axis in such as Fig. 7), its turn Dynamic angle is corresponding to the φ axle in Fig. 1, and the Aux2 mouth on 1800 controllers can only connect to the horizontal polarization test of position 6 Antenna, In/Out mouth may be coupled to any one test antenna of 1 to position, position 11, and Aux1 connects communication sky, darkroom Line, the vertical polarization test antenna of Aux3 link position 6.As seen from Figure 7, a downlink of UXM is connected to Aux2 Mouthful, i.e. horizontal polarization test antenna at link position 6, another downlink is connected to In/Out mouth, can select position Put any one test antenna of 1 to position 11, test antenna including the horizontal polarization at position 6.

When testing, first obtaining the pattern information of wireless terminal, due in the present embodiment, wireless terminal has two Individual reception antenna the most named Main antenna and Sub antenna, therefore pattern information includes that the H polarization of Main antenna increases Benefit and V polarized gain, the H polarized gain of Sub antenna and V polarized gain and the phase information of two antennas.Such as, As Fig. 8 (1) and (2) are respectively Main antenna and the 2D pattern information of Sub antenna of a terminal, wherein test Angle is 0 to 360 ° in θ direction and is spaced apart 30 °, is fixed as 90 ° on φ direction.

Wherein, in order to represent Main antenna and the phase contrast of Sub antenna, it is stipulated that the phase place of Main antenna is 0 °, Sub days Phase of line is the phase contrast of two reception antennas.

Secondly, calculate according to the pattern information of two reception antennas of above-mentioned acquisition and be suitable for solving the coordinate points of inverse matrix, and Inverse matrix is solved in this coordinate points.Specifically, definition coordinateWherein, Ant₁,Ant₂Expression is chosen Test antenna in Fig. 7,The angle rotated for turntable.Travel through the coordinate points of all aforementioned definitions, carry out following operation: In order to more fully hereinafter illustrate, below with coordinate (6_H,6_V, 90) as a example by, wherein, 6_HRepresent the H polarization test antenna of position 6, 6_VRepresenting the V polarization test antenna of position 6,90 represent that turntable rotates 90 °.

Extract MIMO wireless terminal in the amplitude of this coordinate and phase information, be expressed as:

Obtain after being loaded with corresponding inverse matrix, be capable of two virtual wires according to acquisition amplitude and phase information computing Isolation for example: I_SO1=33dB, I_SO2=23dB, and obtain the virtual path loss of each virtual wires for example: L_p1=2dB, L_p2=12dB, it can be seen that the absolute value of the isolation of all virtual wires of this coordinate is both greater than 6dB, I.e. this coordinate of labelling is " qualified " coordinate points.

After traveling through coordinate points defined above, in " qualified " coordinate points of all confirmations, according to all on this coordinate Virtual wires isolation or virtual path loss or both integrated value determine optimum, and solve inverse matrix in optimum.

By specific embodiments above it can be seen that the inverse matrix method for solving of the embodiment of the present invention can be with direct basis reception antenna Pattern information judge that whether each coordinate meets the condition properly solving inverse matrix, and the void that institute's finding the inverse matrix is corresponding Intend the maximum isolation degree i.e. capacity of resisting disturbance of wire.The method of the embodiment of the present invention can solve to radiate can not intelligence in two-step method Ground, the problem of inverse matrix of automatic calculation space matrix.

In the explanation of this specification, describe can be by for any process described otherwise above or method in flow chart or at this It is interpreted as, represents and include one or more code for the executable instruction of the step that realizes specific logical function or process Module, fragment or part, and the scope of the preferred embodiment of the present invention includes other realization, wherein can not press Order that is shown or that discuss, including according to involved function by basic mode simultaneously or in the opposite order, performs Function, this should be understood by embodiments of the invention person of ordinary skill in the field.

Represent in flow charts or the logic described otherwise above at this and/or step, for example, it is possible to be considered as reality The sequencing list of the executable instruction of existing logic function, may be embodied in any computer-readable medium, for instruction Execution system, device or equipment (system such as computer based system, including processor or other can perform from instruction System, device or equipment instruction fetch also perform the system instructed) use, or combine these instruction execution systems, device or set Standby and use.For the purpose of this specification, " computer-readable medium " can be any can to comprise, store, communicate, propagate Or transmission procedure for instruction execution system, device or equipment or combines these instruction execution systems, device or equipment and uses Device.The more specifically example (non-exhaustive list) of computer-readable medium includes following: have one or more cloth The electrical connection section (electronic installation) of line, portable computer diskette box (magnetic device), random access memory (RAM), read-only Memorizer (ROM), erasable read only memory (EPROM or flash memory) of editing, fiber device, and portable Compact disc read-only memory (CDROM).It addition, computer-readable medium can even is that the paper that can print described program thereon Or other suitable media, because then can carry out editing, solving such as by paper or other media are carried out optical scanning Translate or be processed to electronically obtain described program with other suitable methods if desired, be then stored in computer In memorizer.

Should be appreciated that each several part of the present invention can realize by hardware, software, firmware or combinations thereof.In above-mentioned reality Execute in mode, software that multiple steps or method in memory and can be performed by suitable instruction execution system with storage or Firmware realizes.Such as, if realized with hardware, with the most the same, available well known in the art under Any one or their combination in row technology realize: have the logic gates for data signal realizes logic function Discrete logic, there is the special IC of suitable combination logic gate circuit, programmable gate array (PGA), existing Field programmable gate array (FPGA) etc..

Those skilled in the art are appreciated that realizing all or part of step that above-described embodiment method carries is can Completing instructing relevant hardware by program, described program can be stored in a kind of computer-readable recording medium, This program upon execution, including one or a combination set of the step of embodiment of the method.

Additionally, term " first ", " second " are only used for describing purpose, and it is not intended that instruction or hint relative importance Or the implicit quantity indicating indicated technical characteristic.Thus, define " first ", the feature of " second " can be expressed Or implicitly include at least one this feature.In describing the invention, " multiple " are meant that at least two, such as two Individual, three etc., unless otherwise expressly limited specifically.

In the description of this specification, reference term " embodiment ", " some embodiments ", " example ", " specifically show Example " or the description of " some examples " etc. means to combine this embodiment or example describes specific features, structure, material or Feature is contained at least one embodiment or the example of the present invention.In this manual, the schematic representation to above-mentioned term Necessarily it is directed to identical embodiment or example.And, the specific features of description, structure, material or feature are permissible One or more embodiment in office or example combine in an appropriate manner.Additionally, in the case of the most conflicting, ability The feature of the different embodiments described in this specification or example and different embodiment or example can be entered by the technical staff in territory Row combines and combination.

Although above it has been shown and described that embodiments of the invention, it is to be understood that above-described embodiment is exemplary, Being not considered as limiting the invention, those of ordinary skill in the art within the scope of the invention can be to above-described embodiment It is changed, revises, replaces and modification.

Claims (7)

1. the method for solving of an electromagnetic wave propagation matrix inversion matrix based on wireless terminal antenna pattern information, it is characterized in that, n test antenna and turntable is included in microwave dark room, wireless terminal is placed on described turntable, and described wireless terminal includes that m reception antenna, described m reception antenna carry out data transmission with electromagnetic wave with m selected test antenna, wherein, m >=2, n >=m, described method for solving comprises the following steps:

Obtain the pattern information of described m reception antenna；

Pattern information according to described m reception antenna calculates and is suitable for solving the coordinate points of electromagnetic wave propagation inverse of a matrix matrix, and wherein, described coordinate points includes the parameter information affecting described electromagnetic wave propagation matrix；And

In the described coordinate points being suitable for solving electromagnetic wave propagation inverse of a matrix matrix, corresponding electromagnetic wave propagation inverse of a matrix matrix is solved.

2. the method for claim 1, it is characterised in that calculate according to the pattern information of described m reception antenna and be suitable for solving the coordinate points of electromagnetic wave propagation inverse of a matrix matrix, specifically include:

Extract correspondenceAt coordinate points, each reception antenna tests the reception electromagnetic wave phase differences information of antenna, wherein, (Ant relative to gain information and any two reception antenna of each test antenna relative to each₁,Ant₂,...,Ant_m) it is the information of the selected test antenna of described coordinate points,For placing the rotation angle information of the turntable of wireless terminal in described microwave dark room；

According to described gain information and described phase information, the isolation of the m bar virtual wires being loaded with between m reception antenna and m the test antenna realized after corresponding inverse matrix is obtained in computing respectively；

Isolation according to described m bar virtual wires judges corresponding describedWhether coordinate points is the qualified coordinate points solving described inverse matrix.

3. method as claimed in claim 2, it is characterised in that judge corresponding described according to the isolation of described m bar virtual wiresWhether coordinate points is that the qualified coordinate points solving described inverse matrix specifically includes:

If it is corresponding describedThe symbol of the isolation of the m bar virtual wires of coordinate points is identical and absolute value is all higher than predetermined threshold value, then judge that this coordinate points is the qualified coordinate points solving described inverse matrix.

4. method as claimed in claim 3, it is characterised in that also include:

TraversalAll coordinate points of affiliated scope, it is thus achieved that solve all of qualified coordinate points of described inverse matrix；And

The coordinate points most preferably solving electromagnetic wave propagation inverse of a matrix matrix is determined in described all of qualified coordinate points.

5. method as claimed in claim 4, it is characterised in that determine the described coordinate points most preferably solving electromagnetic wave propagation inverse of a matrix matrix in described all of qualified coordinate points, specifically include:

The isolation of the m bar virtual wires of coordinate points qualified each described is taken absolute value, and isolation angle value minimum in absolute value is comprehensively isolated angle value as respective coordinates point；And

Selecting the coordinate points that in described all of qualified coordinate points, comprehensive isolation is maximum is the described coordinate points most preferably solving electromagnetic wave propagation inverse of a matrix matrix.

6. method as claimed in claim 4, it is characterised in that determine the described coordinate points being suitable for solving electromagnetic wave propagation inverse of a matrix matrix in described all of qualified coordinate points, specifically include:

Calculate the virtual path loss of described m bar virtual wires respectively；And

The virtual path loss of the m bar virtual wires of coordinate points qualified each described is taken absolute value, and using the virtual circuit loss value of maximum in absolute value as the integrated virtual circuit loss value of respective coordinates point；

Select integrated virtual circuit loss value in described all of qualified coordinate points minimum coordinate points as the described coordinate points most preferably solving electromagnetic wave propagation inverse of a matrix matrix.

7. method as claimed in claim 4, it is characterised in that determine the described coordinate points most preferably solving electromagnetic wave propagation inverse of a matrix matrix in described all of qualified coordinate points, specifically include:

The isolation of the m bar virtual wires of coordinate points qualified each described is taken absolute value, and isolation angle value minimum in absolute value is comprehensively isolated angle value as respective coordinates point；

Calculate the virtual path loss of described m bar virtual wires respectively, the virtual path loss of the m bar virtual wires of coordinate points qualified each described is taken absolute value, and using the virtual circuit loss value of maximum in absolute value as the integrated virtual circuit loss value of respective coordinates point；

Comprehensive isolation and integrated virtual circuit loss value according to described all of qualified coordinate points determine the described coordinate points most preferably solving electromagnetic wave propagation inverse of a matrix matrix.