CN105842546B - Electromagnetic Wave Propagation inverse of a matrix matrix solving method based on antenna radiation pattern information - Google Patents
Electromagnetic Wave Propagation inverse of a matrix matrix solving method based on antenna radiation pattern information Download PDFInfo
- Publication number
- CN105842546B CN105842546B CN201510531984.0A CN201510531984A CN105842546B CN 105842546 B CN105842546 B CN 105842546B CN 201510531984 A CN201510531984 A CN 201510531984A CN 105842546 B CN105842546 B CN 105842546B
- Authority
- CN
- China
- Prior art keywords
- matrix
- coordinate points
- electromagnetic wave
- inverse
- antenna
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Landscapes
- Mobile Radio Communication Systems (AREA)
- Variable-Direction Aerials And Aerial Arrays (AREA)
Abstract
The method for solving for the Electromagnetic Wave Propagation matrix inverse matrix based on wireless terminal antenna pattern information that the invention discloses a kind of, wherein, it include n test antenna in microwave dark room, wireless terminal includes m receiving antenna, m receiving antenna is carried out data transmission with m selected test antenna with electromagnetic wave, m >=2, n >=m.This method is to Electromagnetic Wave Propagation matrix inversion Matrix Solving the following steps are included: obtaining the pattern information of m receiving antenna；Being calculated according to above-mentioned pattern information is suitble to the coordinate points for solving Electromagnetic Wave Propagation inverse of a matrix matrix to solve in above-mentioned suitable coordinate points to the inverse matrix of corresponding propogator matrix.The method for solving can rapidly solve the inverse matrix met the requirements when the radiation two-step method for carrying out one of MIMO test international standard calculates, be not necessarily to human intervention, time saving and energy saving.
Description
Technical field
The invention belongs to wireless device performance the field of test technology, more particularly to one kind to be based on wireless terminal antenna directional diagram
The method for solving and this method of the Electromagnetic Wave Propagation inverse of a matrix matrix of information are surveyed at LTE MIMO OTA (Over The Air)
It is applied in one of examination international standard radiation two-step method.
Technical background
The purpose of MIMO OTA test is to guarantee really reflect wireless terminal various in the test result in laboratory
Wireless performance under complicated practical service environment and user's use state.
Currently, international wireless highest standard tissue, 3GPP (3rd Generation Partnership Project) is mentioned
The main testing standard supplied has Multiple probe technique, radiation two-step method, reverberation method.Wherein, in the related art, one kind is disclosed to be based on
Two-step method, quick, accurate, economic MIMO OTA test Solution is radiated, as shown in Figure 1, terminal test environment includes
Base station emulator, channel simulator, shielding anechoic chamber, turntable, antenna can be realized and complete to LTE end in efficient
The measurement of the performance at end.
Radiation two-step method is that conduction two-step method develops.The realization for conducting two-step method is by calculating in a computer
The signal of each receiver should be reached, is then delivered directly on the receiver by call wire, is realized and Multi probe
Equivalent measurement process.But since what conduction two-step method can not be introduced comprising receiving antenna mutually interferes and does not meet OTA
Test theory, therefore gradually to radiation two-step method development.Radiating two-step method is to remove conduction on the basis of conducting two-step method
Then line, the free-space propagation for allowing signal to simulate in darkroom demodulate the spatial using corresponding method to measured piece
Inverse of a matrix matrix, so, so that it may realize the effect of two virtual call wires in space, make signal in virtual call wire
On be transmitted directly to up to receiver.Therefore, one of the difficulty that radiation two-step method faces is exactly the space matrix to Electromagnetic Wave Propagation
Inverse matrix solution.
The space matrix of Electromagnetic Wave Propagation is not only related with dark room conditions, but also the antenna position with selected transmitting signal
It sets, polarize and related with the posture of measured piece, the combination of these variables can produce very more Electromagnetic Wave Propagation matrixes, by
It is conditional in a matrix with the presence or absence of inverse matrix, therefore, in these propogator matrixes, it is understood that there may be not can be carried out very much
Inverse matrix solve, if can not predict which propogator matrix be suitble to solve inverse matrix, unique way be exactly in sequence take turns sequence or
Whether each possible forward-propagating inverse of a matrix matrix of person's random challenge meets the requirements.Such way is when extremely expending
Between and it is inefficient, often will appear several hours even more long all without training in rotation to satisfaction inverse matrix.It is more fatal, it surveys
Examination personnel must all-the-way tracking test, to the performance of each inverse matrix for solving out make it is subjective whether meet the requirements sentence
It is disconnected.Therefore, the development and popularization of radiation two-step method are seriously hampered.
Summary of the invention
The present invention is directed to solve at least some of the technical problems in related technologies.For this purpose, the present invention needs
It is proposed a kind of method for solving of Electromagnetic Wave Propagation matrix inverse matrix based on wireless terminal antenna pattern information, the calculating side
Method can rapidly solve the inverse matrix met the requirements when inverse matrix calculates in carrying out radiation two-step method, be not necessarily to human intervention,
It is time saving and energy saving.
To solve the above-mentioned problems, the embodiment of the present invention proposes a kind of Electromagnetic Wave Propagation square based on antenna radiation pattern information
The method for solving of battle array inverse matrix, wherein include turntable and n test antenna in microwave dark room, wherein wireless terminal is placed on institute
It states on turntable, the wireless terminal includes m receiving antenna, and the m receiving antenna and m selected test antenna are with electromagnetism
Waveshape carries out data transmission, wherein m >=2, n >=m, the method for solving is the following steps are included: obtain m reception day
The pattern information of line；It is calculated according to the pattern information of the m receiving antenna and is suitble to solve Electromagnetic Wave Propagation inverse of a matrix
The coordinate points of matrix, wherein the coordinate points include the parameter information for influencing the Electromagnetic Wave Propagation matrix；And described suitable
The coordinate points for solving Electromagnetic Wave Propagation inverse of a matrix matrix are closed, corresponding Electromagnetic Wave Propagation inverse of a matrix matrix is solved.
The solution side of Electromagnetic Wave Propagation matrix inverse matrix according to the present invention based on mobile terminal antenna pattern information
Method, can judge whether each coordinate meets the suitable item for solving inverse matrix with the pattern information of direct basis receiving antenna
Part is more quickly computed the coordinate points for being suitble to solve Electromagnetic Wave Propagation inverse of a matrix matrix, in suitable coordinate points to inverse square
Battle array solves, and intelligent can realize, is not necessarily to human intervention, time saving and energy saving.
Further, it is calculated according to the pattern information of the m receiving antenna and is suitble to solve Electromagnetic Wave Propagation matrix
The coordinate points of inverse matrix specifically include: extracting and correspond toEach receiving antenna phase of coordinate points
Gain information and any two receiving antenna for each test antenna test the reception electromagnetic wave of antenna relative to each
Phase information, wherein (Ant1,Ant2,...,Antm) it is the information for testing antenna that the coordinate points are selected,For
The rotation angle information of the turntable of wireless terminal is placed in the microwave dark room；Believed according to the gain information and the phase difference
Breath, operation finds out the m item being loaded between the m receiving antenna and m test antenna that corresponding inverse matrix is realized later respectively
The isolation of virtual wires；According to the isolation of m virtual wires judgement correspondenceWhether coordinate points are the qualified coordinate points for solving the inverse matrix.
Further, according to m isolation judgement correspondenceCoordinate points
It whether is to solve the qualified coordinate points of the inverse matrix to specifically include: if described in corresponding
The sign symbol of the isolation of m virtual wires of coordinate points is identical and absolute value is all larger than preset threshold, then judges the coordinate
Point is the qualified coordinate points for solving the inverse matrix.
The above method further include: traversalAll coordinate points of affiliated range obtain
Solve all qualified coordinate points of the inverse matrix；And the best solution is determined in all qualified coordinate points
The coordinate points of Electromagnetic Wave Propagation inverse of a matrix matrix.
Further, the best solution Electromagnetic Wave Propagation inverse of a matrix square is determined in all qualified coordinate points
The coordinate points of battle array specifically include: taking absolute value the isolation of m virtual wires of each qualified coordinate points, and will be exhausted
Comprehensive isolation angle value to isolation angle value the smallest in value as respective coordinates point；It selects comprehensive in all qualified coordinate points
Closing the maximum coordinate points of isolation is the best coordinate points for solving Electromagnetic Wave Propagation inverse of a matrix matrix.
Alternatively, calculating separately the virtual path loss of the m virtual wires；And the m each qualified coordinate points
The virtual path loss of virtual wires takes absolute value, and using virtual circuit loss value maximum in absolute value as the synthesis of respective coordinates point
Virtual circuit loss value selects the coordinate points that integrated virtual circuit loss value is the smallest in all qualified coordinate points as described best
Solve the coordinate points of Electromagnetic Wave Propagation inverse of a matrix matrix.
Alternatively, taking absolute value the isolation of m virtual wires of each qualified coordinate points, and will be in absolute value
Comprehensive isolation angle value of the smallest isolation angle value as respective coordinates point；The virtual path loss of the m virtual wires is calculated separately,
It takes absolute value the virtual path loss of m virtual wires of each qualified coordinate points, and will be maximum virtual in absolute value
Integrated virtual circuit loss value of the circuit loss value as respective coordinates point；According to the synthesis isolation of all qualified coordinate points and comprehensive
It closes virtual circuit loss value and determines the best coordinate points for solving Electromagnetic Wave Propagation matrix inverse matrix.
Detailed description of the invention
Fig. 1 is the partial structure diagram of wireless terminal test environment according to an embodiment of the invention；
Fig. 2 is the realization schematic diagram of conduction two-step method in the related technology；
Fig. 3 is that radiation two-step method realizes system schematic in the related technology；
Fig. 4 is that spacing wave propagates schematic diagram in radiation two-step method according to an embodiment of the invention；
Fig. 5 be in two receiver practical work process according to an embodiment of the invention signal propagate schematic diagram and
Its equivalent schematic；
Fig. 6 is Electromagnetic Wave Propagation matrix according to one embodiment of present invention based on wireless terminal antenna pattern information
The flow chart of the method for solving of inverse matrix；
Fig. 7 is that the realization of a specific embodiment according to the present invention is shown for the configuration system of method for solving described in Fig. 6
It is intended to；
(1) and (2) in Fig. 8 is the directional diagram letter of two receiving antennas of another specific embodiment according to the present invention
The data form of breath.
Specific embodiment
The embodiment of the present invention is described below in detail, examples of the embodiments are shown in the accompanying drawings, wherein from beginning to end
Same or similar label indicates same or similar element or element with the same or similar functions.Below with reference to attached
The embodiment of figure description is exemplary, it is intended to is used to explain the present invention, and is not considered as limiting the invention.
It can be used in radiation two-step method to state the inverse matrix that the calculation method of the embodiment of the present invention is sought, first
Simple introduction is done to the concept of several definitionals involved in radiation two-step method and the embodiment of the present invention.
By taking two receivers as an example, the realization for conducting two-step method can be as shown in Fig. 2, by the direction of channel model and antenna
The operation of figure information is run in instrument, then operation result is directly inputted into the receiver of measured piece by call wire, i.e.,
The effect reached can indicate are as follows:
[R1 R2]=[T1 T2] (1)
Wherein, R1And R2Respectively indicate the signal that two receivers receive, T1And T2It is expressed as channel simulator
Two-way output signal.
As shown in figure 3, being propagated and being equivalent in free-space propagation, such as in darkroom due to signal for radiating two-step method
Shown in Fig. 4, it is assumed that spatial matrix are as follows:
Then emitting signal and receiving signal relation indicates are as follows:
If reaching effect same as two-step method is conducted, can transmitting signal launch before, in instrument into
The following operation of row, multiplied by a new matrix:
The then relationship of new transmitting signal and reception signal are as follows:
Then when meeting condition:
Wherein, x is not zero, so that it may realize radiation two-step method and conduct the equivalent of two-step method, that is, reach the effect of formula (1)
Fruit, it can it is equivalent to think that there are two " conducting wires " in space, such as it is referred to as " virtual wires " in an embodiment of the present invention,
Signal can be directly inputted to receiver by this two virtual wires, as shown in Figure 3.Above formula statement for M (M >=
3) equally applicable the case where a receiver.
The positive Electromagnetic Wave Propagation matrix propagated in darkroom is related with selected test antenna, measured piece rotation angle.With 2
For the positive matrix of 2 ╳ 2 of receiver, coordinate is definedWherein, Ant1Indicate first selected test
Antenna, Ant2Indicate that selected second test antenna, θ indicate the angle that measured piece rotates on θ axis as shown in figure 1,Table
Show the angle that measured piece rotates on φ axis as shown in figure 1, then each coordinate pair answers a positive Electromagnetic Wave Propagation matrix.
But not the corresponding positive Electromagnetic Wave Propagation matrix of each coordinate there are inverse matrix or is not each
Inverse matrix can meet actual demand.In view of this, the index of one " isolation " is proposed in an embodiment of the present invention.With 2
In case where a receiver, after inverse matrix is added, shown in real work situation such as Fig. 5 (1) of 2 receivers, if inverse
Matrix C, C are written in instrument, are expressed as above-mentioned formula (4) form, and space propogator matrix is expressed as H in darkroom, and H is indicated are as follows:
Wherein, when meeting H*C equal to formula (6), in Fig. 5 (1) situation can be equivalent to as
Shown in Fig. 5 (2).Actually, receiver 1 (R1) can not only receive the signal of transmitter 1 (T1) sending, if above-mentioned power is big
It is small to be expressed as P11, unit dBm, receiver 1 (R1) can also receive the interference of transmitter 2 (T2), if above-mentioned set watt level
For P12, then it is assumed that the jamproof isolation size that virtual wires 1 can be realized are as follows:
Iso1=P11-P12 (7)
Similarly, the signal for the transmitter 2 (T2) for enabling receiver 2 (R2) receive is P22, receiver 2 (R2) receives
The interference of transmitter 1 (T1) is P21, then it is assumed that the jamproof isolation size that virtual wires 2 can be realized are as follows:
Iso2=P22-P21 (8)
When two isolation sizes meet certain requirements, it is believed that the interference signal from another way is to receiver work
The influence of work can be ignored, that is, think that same test result can be obtained with conduction measurement in the case.It should be noted that
Although isolation is not the sole indicator of virtual wires, isolation is to solve for whether virtual wires meet radiation two-step method
It is required that necessary index.In an embodiment of the present invention, isolation is still, right in practical applications with the statement of dB data mode
Data format does not require, and real number expression can also be with.
On the basis of virtual wires meet insulated degree requirement, " virtual path loss " is also defined 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 issues is respectively PT1And PT2, with
DBm data format indicates that the watt level that two receivers of R1 and R2 receive is respectively PR1And PR2, then virtual wires 1 are defined
It is respectively as follows: with the virtual path loss of virtual wires 2
Lp1=PT1-PR1 (9)
Lp2=PT2-PR2 (10)
The meaning of virtual path loss can assess loss of the virtual wires to power it can be seen from above-mentioned two formula.
In the above-described embodiments, coordinate pointsNumber it is very more.But due to that can not solve in advance
Each coordinate pair answers the isolation of its corresponding virtual wires after the inverse matrix of propogator matrix, can not just judge the coordinate points
If appropriate for solution inverse matrix.Radiation two-step method mainly takes traversal or random manner to the selection of coordinate, that is, passes through object
Reason rotation (and radio frequency switching) traverses all coordinate, then successively solve inverse matrix and verify corresponding virtual wires every
From degree, until finding satisfactory i.e. stopping.This mode very elapsed time and cannot be guaranteed find so that virtual wires
Isolation reaches maximum coordinate, thus greatly reduces the efficiency and anti-interference ability of radiation two-step method test.
In view of the above-mentioned problems, the embodiment of the present invention proposes a kind of new inverse matrix calculation method, can quickly, intelligently look for
To satisfactory coordinate points, and then to the corresponding forward-propagating inverse of a matrix Matrix Solving of the coordinate points.
Asking to the Electromagnetic Wave Propagation matrix inverse matrix of the wireless terminal antenna pattern information of the embodiment of the present invention below
Solution method is illustrated.Wherein, data transmission carries out in microwave dark room, includes turntable and n test antenna in microwave dark room,
Wireless terminal is placed on turntable, and the wireless terminal includes m receiving antenna, m >=2, the n test days in microwave dark room
M test antenna, n >=m are selected in line, the m receiving antenna carries out data biography with m test antenna with electromagnetic wave
Defeated, i.e., m receiving antenna receives electromagnetic wave, m test aerial radiation electromagnetic wave, to realize that data are transmitted, therefore micro- at this time
Electromagnetic Wave Propagation in wave darkroom between m receiving antenna and m test antenna amplitude and phase variation can with a m ×
The math matrix of m indicates that the method for solving of the embodiment of the present invention solves the inverse matrix of the math matrix.
As shown in fig. 6, the inverse matrix method for solving of the embodiment of the present invention the following steps are included:
S1 obtains the pattern information of m receiving antenna.
Wherein, antenna radiation pattern is one of the performance of antenna.Specifically, as shown in Figure 1, tested wireless terminal (DUT) is put
It sets at the center of a turntable, tests the distance between antenna and tested wireless terminal and meet standard regulation.Rotating and radio terminal
When, test antenna measurement obtains tested wireless terminal in space all directions antenna gain and any two receiving antenna
Phase difference relative to each test antenna.The measurement of the step can be returned by Radio Terminal Function such as mobile phone
It arrives.
S2 calculates the seat for being suitble to solve Electromagnetic Wave Propagation inverse of a matrix matrix according to the pattern information of m receiving antenna
Punctuate.
Wherein, coordinate points include the parameter information for influencing Electromagnetic Wave Propagation matrix, that is, include selected test antenna and nothing
Position of line terminal etc. can change all information of propogator matrix, for example, the location information of m test antenna and wireless terminal
Rotation angle information.Specifically, coordinate is definedWherein, (Ant1,Ant2,...,Antm)
For test antenna location information,For the rotation angle information of wireless terminal, m is the wireless of progress MIMO OTA test
The number of the receiving antenna of terminal, Ant1,Ant2...,AntmIt indicates the m chosen dark indoor test antenna, includes in darkroom
At least m test antenna, wherein m >=2, θ are rotation angle of the wireless terminal as shown in Figure 1 on θ axis, the step of rotation
Into size, slewing area is according to depending on darkroom turret design situation.For rotation of the wireless terminal as shown in Figure 1 on φ axis
Gyration, the step size of rotation, slewing area are related to depending on situation according to darkroom turntable.
Specifically, it extracts and corresponds toThe each receiving antenna of coordinate points is relative to each survey
Phase information, that is, any two for trying the reception electromagnetic wave between the gain information and any two receiving antenna of antenna receive
Phase difference of the antenna relative to the reception signal of each test antenna.For example, extraction wireless terminal is in turntableAngle
When, relative to the antenna radiation pattern information of each test antenna, antenna Ant is tested relative to any onei(0 < i≤m) is extracted
The pattern information of receiving antenna of wireless terminal include the yield value of m receiving antenna and the phase of any two receiving antenna
Potential difference.By taking Fig. 4 as an example, wherein T1,T2,...,TmFor selected m different test antennas, R1,R2,...,RmIt is connect for m
Antenna is received, the amplitude and phase difference information extracted can be expressed as:
Wherein, αxyIt indicates gain (real number representation),Indicate phase difference, above is a kind of representation, and unlimited
The representation of gain and any two receiving antenna phase difference is not required in this unique representation, and herein.
Further, according to the gain information of said extracted and phase information, operation, which is found out, respectively is loaded with accordingly
The isolation for the m virtual wires between m receiving antenna and m test antenna realized after inverse matrix.
The calculation formula of the isolation of above-mentioned m virtual wires and the value of m are in relation to (its definition may refer to formula (7)
With formula (8)), when m difference, the isolation calculation formula of virtual wires is different, and the application focuses on when providing m=2 (i.e.
For the method for solving of 2 × 2 propogator matrix inverse matrixs) 2 virtual wires isolation solution formula it is as follows.
As m=2, virtual wires isolation solves explanation:
Specifically, it extracts and corresponds toIncreasing of two receiving antennas of coordinate points relative to each test antenna
The phase information of reception electromagnetic wave between beneficial information and above-mentioned two receiving antenna.For example, extracting wireless terminal in turntable
ForWhen angle, relative to the antenna radiation pattern information of each test antenna, antenna Ant is tested relative to any onei(0
< i≤m=2), the pattern information of the receiving antenna of the wireless terminal of extraction includes that the yield value of 2 receiving antennas and two connect
Receive the phase information of antenna.Wherein it is possible to which the amplitude and phase difference information extracted is expressed as:
Wherein, αxyIt indicates gain (real number representation),Indicate phase difference.It is then defined by formula (7) and (8) virtual
The specific formula for calculation of the size of conducting wire isolation is as follows: the jamproof isolation size that virtual wires 1 can be realized be (every
Indicated from degree with dB digital form):
The jamproof isolation size that virtual wires 2 can be realized is (isolation is indicated with dB digital form):
Wherein, αxyFor gain information, δxyFor phase information, x=1,2, y=1,2.
Further, judge to correspond to according to the isolation of m virtual wiresCoordinate points
It whether is the qualified coordinate points for solving inverse matrix.Specifically, if it is correspondingThe m of coordinate points
The sign symbol of the isolation of virtual wires is identical and absolute value is all larger than preset threshold such as 6dB, then judges the coordinate points
For solve inverse matrix qualified coordinate points, such as mark the coordinate points be " qualification " coordinate points.
It should be noted that needing to be traversed forAll coordinate points of affiliated range, are obtained
All qualified coordinate points of inverse matrix must be solved, and then determine best solution Electromagnetic Wave Propagation in all qualified coordinate points
The coordinate points of inverse of a matrix matrix find Best Point and solve inverse matrix.
It is not particularly limited herein for the judgment criteria of Best Point, in one embodiment of the invention, to each
A coordinate points define the concept of comprehensive isolation, and value is equal to the exhausted of the isolation of all virtual wires in a coordinate points
To the minimum value of value.Specifically, it takes absolute value to the isolation of m virtual wires of each qualified coordinate points, and will be absolute
Comprehensive isolation angle value of the smallest isolation angle value as respective coordinates point in value, and then can choose in all qualified coordinate points
The comprehensive maximum coordinate points of isolation are the coordinate points for being most suitable for solving Electromagnetic Wave Propagation inverse of a matrix matrix, i.e. solution inverse matrix
Best Point.
In addition, each coordinate points are defined with the concept of integrated virtual path loss, value is equal to all void in coordinate points
The maximum value of the absolute value of the virtual path loss of quasi- conducting wire.Specifically, the virtual path loss for calculating separately m virtual wires, for void
The calculating of quasi- path loss is referred to formula (9) and (10), repeats no more herein, in turn, to the m item of each qualified coordinate points
The virtual path loss of virtual wires takes absolute value, and using virtual circuit loss value maximum in absolute value as the comprehensive empty of respective coordinates point
Quasi- circuit loss value selects the smallest coordinate points of integrated virtual path loss in all qualified coordinate points to pass as solution electromagnetic wave is most suitable for
Broadcast the coordinate points of inverse of a matrix matrix.
Alternatively, comprehensively considering comprehensive isolation and integrated virtual path loss to determine best coordinates, specifically, to each conjunction
The isolation of m virtual wires of lattice coordinate points takes absolute value, and using isolation angle value the smallest in absolute value as respective coordinates
The comprehensive isolation angle value of point, and, the virtual path loss of m virtual wires is calculated separately, the m of each qualified coordinate points
The virtual path loss of virtual wires takes absolute value, and using virtual circuit loss value maximum in absolute value as the synthesis of respective coordinates point
Virtual circuit loss value, and then best solution electromagnetism is determined according to isolation comprehensive in all qualified coordinate points and integrated virtual path loss
The coordinate points of the inverse matrix of wave propogator matrix, such as the integrated value of the two is determined into Best Point coordinate as criterion.
S3 is asked in the inverse matrix for being suitble to the coordinate points for solving Electromagnetic Wave Propagation inverse of a matrix matrix to correspond to propogator matrix
Solution.
In summary, wireless terminal is in microwave dark room or free space, receiving antenna and darkroom test antenna into
When row data are transmitted, the amplitude and phase variation of electromagnetic transmission can be indicated with one group of propogator matrix, the propogator matrix
It is related relative to the test position of antenna with the selection of test antenna and wireless terminal.The method for solving of the embodiment of the present invention, base
In the pattern information of the receiving antenna of wireless terminal, all propogator matrixes and inverse matrix phase can quickly, be intelligently calculated
The isolation and virtual path loss of the virtual wires formed after multiplying, and then according to the isolation size of virtual wires, it can plan
It is suitble to do the combination of the test antenna of inverse matrix and wireless terminal rotation angle out.
Below by taking the MIMO OTA performance test that two hairs two are received as an example to the inverse matrix method for solving of the embodiment of the present invention into
Row explanation.
In practical applications, as shown in fig. 7, being obtained according to system configuration as shown in the figure, system shown in Figure 7 configuration is said
Bright as follows: test instrumentation is Keysight E7515A (Chinese is referred to as UXM), is integrated with anolog base station and mimo channel is imitative
True instrument.Microwave dark room includes: cabinet, suction wave cotton, test antenna, turntable, RF switch, controller etc..Wherein, antenna position is tested
At position 1 to the position 11 marked in Fig. 7.There are two cross-polarized test antenna, communication antenna positions at each position
At top of the box position 0, turntable can only be along an axis rotation (such as Z axis in Fig. 7), and the angle rotated corresponds to
φ axis in Fig. 1, the Aux2 mouth on 1800 controllers can only connect to the horizontal polarization test antenna of position 6, and In/Out mouthfuls can
Be connected to position 1 to position 11 any one test antenna, Aux1 connection darkroom communication antenna, Aux3 link position 6
Vertical polarization tests antenna.As seen from Figure 7, a downlink of UXM is connected to Aux2 mouthfuls, i.e. level at link position 6
Polarization test antenna, another downlink are connected to In/Out mouthfuls, can choose any one test of position 1 to position 11
Antenna tests antenna including the horizontal polarization at position 6.
When testing, the pattern information for obtaining wireless terminal first, since in the present embodiment, wireless terminal has
Two receiving antennas are for example named as Main antenna and Sub antenna, therefore pattern information includes the H polarized gain of Main antenna
With V polarized gain, the H polarized gain of Sub antenna and the phase information of V polarized gain and two antennas.For example, such as Fig. 8
(1) and (2) are respectively the Main antenna of a terminal and the 2D pattern information of Sub antenna, and wherein test angle is in the direction θ
It is divided into 30 ° between 0 to 360 °, 90 ° are fixed as on the direction φ.
Wherein, in order to indicate Main antenna and Sub antenna phase difference, it is specified that Main antenna phase be 0 °, Sub antenna
Phase is the phase difference of two receiving antennas.
Secondly, calculating the coordinate for being suitble to solve inverse matrix according to the pattern information of the two of above-mentioned acquisition receiving antennas
Point, and inverse matrix is solved in the coordinate points.Specifically, coordinate is definedWherein, Ant1,Ant2It indicates to choose
Fig. 7 in test antenna,For the angle of turntable rotation.The coordinate points for traversing all aforementioned definitions, perform the following operation: for
More fully hereinafter illustrate, below with coordinate (6H,6V, 90) for, wherein 6HThe H polarization test antenna of expression position 6,6VTable
Show the V polarization test antenna of position 6,90 indicate that turntable rotates 90 ° of
MIMO wireless terminal is extracted in the amplitude and phase information of the coordinate, is indicated are as follows:
It finds out to be loaded with after corresponding inverse matrix according to acquisition amplitude and phase difference information operation and can be realized two void
The isolation of quasi- conducting wire is for example are as follows: ISO1=33dB, ISO2=23dB, and find out the virtual path loss of each virtual wires for example
Are as follows: Lp1=2dB, Lp2=12dB, it can be seen that the absolute value of the isolation of all virtual wires of the coordinate is both greater than 6dB,
Marking the coordinate is " qualification " coordinate points.
After traversing coordinate points defined above, in " qualification " coordinate points of all confirmations, according to institute on the coordinate
There are virtual wires isolation or virtual path loss or both integrated value to determine Best Point, and solves inverse matrix in Best Point.
The inverse matrix method for solving of the embodiment of the present invention can be received it can be seen from specific embodiments above with direct basis
The pattern information of antenna judges whether each coordinate meets the suitable condition for solving inverse matrix and institute's finding the inverse matrix is corresponding
Virtual wires maximum isolation degree, that is, anti-interference ability.The method of the embodiment of the present invention can solve cannot in radiation two-step method
Intelligently, the problem of the inverse matrix of automatic calculation space matrix.
In the explanation of this specification, any process described otherwise above or method description can in flow chart or herein
To be understood to, indicate to include the steps that one or more for realizing the executable instruction of specific logical function or process
Module, segment or the part of code, and the range of the preferred embodiment of the present invention includes other realization, wherein can not
By sequence shown or discussed, including according to related function by it is basic simultaneously in the way of or in the opposite order, to hold
Row function, this should be understood by the embodiment of the present invention person of ordinary skill in the field.
Expression or logic and/or step described otherwise above herein in flow charts, for example, being considered use
In the order list for the executable instruction for realizing logic function, may be embodied in any computer-readable medium, for
Instruction execution system, device or equipment (such as computer based system, including the system of processor or other can be held from instruction
The instruction fetch of row system, device or equipment and the system executed instruction) it uses, or combine these instruction execution systems, device or set
It is standby and use.For the purpose of this specification, " computer-readable medium ", which can be, any may include, stores, communicates, propagates or pass
Defeated program is for instruction execution system, device or equipment or the dress used in conjunction with these instruction execution systems, device or equipment
It sets.The more specific example (non-exhaustive list) of computer-readable medium include the following: there is the electricity of one or more wirings
Interconnecting piece (electronic device), portable computer diskette box (magnetic device), random access memory (RAM), read-only memory
(ROM), erasable edit read-only storage (EPROM or flash memory), fiber device and portable optic disk is read-only deposits
Reservoir (CDROM).In addition, computer-readable medium can even is that the paper that can print described program on it or other are suitable
Medium, because can then be edited, be interpreted or when necessary with it for example by carrying out optical scanner to paper or other media
His suitable method is handled electronically to obtain described program, is then stored in computer storage.
It should be appreciated that each section of the invention can be realized with hardware, software, firmware or their combination.Above-mentioned
In embodiment, software that multiple steps or method can be executed in memory and by suitable instruction execution system with storage
Or firmware is realized.It, and in another embodiment, can be under well known in the art for example, if realized with hardware
Any one of column technology or their combination are realized: having a logic gates for realizing logic function to data-signal
Discrete logic, with suitable combinational logic gate circuit specific integrated circuit, programmable gate array (PGA), scene
Programmable gate array (FPGA) etc..
Those skilled in the art are understood that realize all or part of step that above-described embodiment method carries
It suddenly is that relevant hardware can be instructed to complete by program, the program can store in a kind of computer-readable storage medium
In matter, which when being executed, includes the steps that one or a combination set of embodiment of the method.
In addition, term " first ", " second " are used for descriptive purposes only and cannot be understood as indicating or suggesting relative importance
Or implicitly indicate the quantity of indicated technical characteristic.Define " first " as a result, the feature of " second " can be expressed or
Implicitly include at least one this feature.In the description of the present invention, the meaning of " plurality " is at least two, such as two, three
It is a etc., unless otherwise specifically defined.
In the description of this specification, reference term " one embodiment ", " some embodiments ", " example ", " specifically show
The description of example " or " some examples " etc. means specific features, structure, material or spy described in conjunction with this embodiment or example
Point is included at least one embodiment or example of the invention.In the present specification, schematic expression of the above terms are not
It must be directed to identical embodiment or example.Moreover, particular features, structures, materials, or characteristics described can be in office
It can be combined in any suitable manner in one or more embodiment or examples.In addition, without conflicting with each other, the skill of this field
Art personnel can tie the feature of different embodiments or examples described in this specification and different embodiments or examples
It closes and combines.
Although the embodiments of the present invention has been shown and described above, it is to be understood that above-described embodiment is example
Property, it is not considered as limiting the invention, those skilled in the art within the scope of the invention can be to above-mentioned
Embodiment is changed, modifies, replacement and variant.
Claims (7)
1. a kind of method for solving of the Electromagnetic Wave Propagation matrix inverse matrix based on wireless terminal antenna pattern information, feature exist
In including n test antenna and turntable in microwave dark room, wireless terminal is placed on the turntable, and the wireless terminal includes m
A receiving antenna, the m receiving antenna are carried out data transmission with m selected test antenna with electromagnetic wave, wherein m
>=2, n >=m, the method for solving the following steps are included:
Obtain the pattern information of the m receiving antenna；
The coordinate for being suitble to solve Electromagnetic Wave Propagation inverse of a matrix matrix is calculated according to the pattern information of the m receiving antenna
Point, wherein the coordinate points include the parameter information for influencing the Electromagnetic Wave Propagation matrix；And
In the suitable coordinate points for solving Electromagnetic Wave Propagation inverse of a matrix matrix, to corresponding Electromagnetic Wave Propagation inverse of a matrix square
Battle array is solved.
2. the method as described in claim 1, which is characterized in that calculated according to the pattern information of the m receiving antenna suitable
The coordinate points for solving Electromagnetic Wave Propagation inverse of a matrix matrix are closed, are specifically included:
It extracts and corresponds toAt coordinate points, each receiving antenna is relative to each test antenna
The reception electromagnetic wave phase differences information of gain information and any two receiving antenna relative to each test antenna, wherein
(Ant1,Ant2,...,Antm) it is the information for testing antenna that the coordinate points are selected, θ is rotation of the wireless terminal on θ axis
Angle,Exist for wireless terminalRotation angle on axis,For the turntable of placement wireless terminal in the microwave dark room
Rotation angle information；
According to the gain information and the phase information, operation, which is found out, respectively is loaded with what corresponding inverse matrix was realized later
The isolation of m virtual wires between m receiving antenna and m test antenna；
According to the isolation of m virtual wires judgement correspondenceWhether coordinate points
For the qualified coordinate points for solving the inverse matrix.
3. method according to claim 2, which is characterized in that judge corresponding institute according to the isolation of the m virtual wires
It statesWhether coordinate points are to solve the qualified coordinate points of the inverse matrix to specifically include:
If described in correspondingThe symbol of the isolation of m virtual wires of coordinate points is identical
And absolute value is all larger than preset threshold, then judges the coordinate points for the qualified coordinate points of the solution inverse matrix.
4. method as claimed in claim 3, which is characterized in that further include:
TraversalAll coordinate points of affiliated range obtain and solve all of the inverse matrix
Qualified coordinate points；And
The best coordinate points for solving Electromagnetic Wave Propagation inverse of a matrix matrix are determined in all qualified coordinate points.
5. method as claimed in claim 4, which is characterized in that determination is described in all qualified coordinate points most preferably asks
The coordinate points for solving Electromagnetic Wave Propagation inverse of a matrix matrix, specifically include:
Take absolute value the isolation of m virtual wires of each qualified coordinate points, and by absolute value it is the smallest every
Comprehensive isolation angle value from angle value as respective coordinates point；And
Selecting the comprehensive maximum coordinate points of isolation in all qualified coordinate points is the best solution Electromagnetic Wave Propagation
The coordinate points of inverse of a matrix matrix.
6. method as claimed in claim 4, which is characterized in that determination is described in all qualified coordinate points is suitble to ask
The coordinate points for solving Electromagnetic Wave Propagation inverse of a matrix matrix, specifically include:
Calculate separately the virtual path loss of the m virtual wires；And
It takes absolute value the virtual path loss of m virtual wires of each qualified coordinate points, and will be maximum in absolute value
Virtual integrated virtual circuit loss value of the circuit loss value as respective coordinates point；
Select the coordinate points that integrated virtual circuit loss value is the smallest in all qualified coordinate points as the best solution electricity
The coordinate points of electromagnetic wave propagation inverse of a matrix matrix.
7. method as claimed in claim 4, which is characterized in that determination is described in all qualified coordinate points most preferably asks
The coordinate points for solving Electromagnetic Wave Propagation inverse of a matrix matrix, specifically include:
Take absolute value the isolation of m virtual wires of each qualified coordinate points, and by absolute value it is the smallest every
Comprehensive isolation angle value from angle value as respective coordinates point；
The virtual path loss for calculating separately the m virtual wires, the void of m virtual wires of each qualified coordinate points
Quasi- path loss takes absolute value, and using virtual circuit loss value maximum in absolute value as the integrated virtual circuit loss value of respective coordinates point；
The best solution electromagnetism is determined according to the synthesis isolation of all qualified coordinate points and integrated virtual circuit loss value
The coordinate points of the inverse matrix of wave propogator matrix.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510531984.0A CN105842546B (en) | 2015-08-26 | 2015-08-26 | Electromagnetic Wave Propagation inverse of a matrix matrix solving method based on antenna radiation pattern information |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510531984.0A CN105842546B (en) | 2015-08-26 | 2015-08-26 | Electromagnetic Wave Propagation inverse of a matrix matrix solving method based on antenna radiation pattern information |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN105842546A CN105842546A (en) | 2016-08-10 |
| CN105842546B true CN105842546B (en) | 2019-02-26 |
Family
ID=56580847
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201510531984.0A Active CN105842546B (en) | 2015-08-26 | 2015-08-26 | Electromagnetic Wave Propagation inverse of a matrix matrix solving method based on antenna radiation pattern information |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN105842546B (en) |
Families Citing this family (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107819527B (en) * | 2016-09-12 | 2020-10-13 | 中国移动通信有限公司研究院 | Testing device and testing method for large-scale antenna base station equipment |
| CN108574539B (en) * | 2017-03-08 | 2021-01-15 | 深圳市通用测试系统有限公司 | Signal generation method and device based on MIMO wireless terminal test |
| CN111372273A (en) * | 2018-12-26 | 2020-07-03 | 是德科技股份有限公司 | Test system for testing antenna of device under test and method for testing user equipment |
| CN111385036B (en) * | 2018-12-28 | 2021-08-24 | 深圳市通用测试系统有限公司 | Radio frequency performance test method and device of wireless equipment and tester |
| CN112188537B (en) * | 2019-07-05 | 2021-09-28 | 中国信息通信研究院 | Near-field wireless channel simulation measurement method and system based on forward optimization |
| CN113111524B (en) * | 2021-04-21 | 2021-09-21 | 北京航空航天大学 | Antenna isolation degree prediction method based on far field test, storage medium and device |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102130725A (en) * | 2010-01-20 | 2011-07-20 | 中兴通讯股份有限公司 | Test method and system for over-the-air (OTA) performance of multi-aerial system |
| US8331869B2 (en) * | 2006-07-24 | 2012-12-11 | Ets Lindgren, L.P. | Systems and methods for over the air performance testing of wireless devices with multiple antennas |
| JP5340344B2 (en) * | 2011-05-12 | 2013-11-13 | パナソニック株式会社 | Communication apparatus and communication method |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP5724496B2 (en) * | 2011-03-17 | 2015-05-27 | 富士通株式会社 | Wireless communication system and wireless communication method |
-
2015
- 2015-08-26 CN CN201510531984.0A patent/CN105842546B/en active Active
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US8331869B2 (en) * | 2006-07-24 | 2012-12-11 | Ets Lindgren, L.P. | Systems and methods for over the air performance testing of wireless devices with multiple antennas |
| CN102130725A (en) * | 2010-01-20 | 2011-07-20 | 中兴通讯股份有限公司 | Test method and system for over-the-air (OTA) performance of multi-aerial system |
| JP5340344B2 (en) * | 2011-05-12 | 2013-11-13 | パナソニック株式会社 | Communication apparatus and communication method |
Non-Patent Citations (1)
| Title |
|---|
| 一种改进的MIMO-OFDM系统Steiner估计算法;路晓菊 等;《信息化研究》;20100531;全文 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN105842546A (en) | 2016-08-10 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN105842546B (en) | Electromagnetic Wave Propagation inverse of a matrix matrix solving method based on antenna radiation pattern information | |
| EP2533572B1 (en) | Method and system for testing over the air performances in multi-antenna system | |
| Zhang et al. | A survey of testing for 5G: Solutions, opportunities, and challenges | |
| CN109617623B (en) | Multi-probe anechoic chamber (MPAC) Over The Air (OTA) test system and method | |
| CN103856272B (en) | The wireless performance method of testing of MIMO wireless terminals | |
| CN103562737B (en) | For the method and apparatus measuring the improvement of the performance of antenna, mobile phone and other wireless terminals | |
| EP2521282B1 (en) | Method and system for spatial radio-frequency performance testing based on multiple-antenna system | |
| US20120282863A1 (en) | Antenna testing system and antenna testing method | |
| CN109617638B (en) | Wireless channel emulator with dynamically variable channel model | |
| Arai | Measurement of mobile antenna systems | |
| CN106788791B (en) | Darkroom multi-wave-surface controller test system, method and device | |
| KR102436859B1 (en) | Method and system for testing wireless performance of wireless terminals | |
| CN102130725A (en) | Test method and system for over-the-air (OTA) performance of multi-aerial system | |
| CN102148649A (en) | Method and system for testing space radio frequency performance of multi-antenna device | |
| CN101667873A (en) | Method and system for testing radio-frequency performance of receiver in multi-antenna channel environment | |
| WO2017215020A1 (en) | Method and device for testing large-scale mimo system base station | |
| CN108574539A (en) | Signal creating method and device based on the test of MIMO wireless terminals | |
| Moss et al. | A simple model of the UHF cross-polar terrestrial channel for DVBNGH | |
| CN111314002B (en) | Automatic test system for omnibearing coverage performance of wireless equipment | |
| CN113573349B (en) | Performance test system and method of millimeter wave communication equipment | |
| WO2011097900A1 (en) | Method and system for spatial radio frequency performance test based on multi-antenna system | |
| Li et al. | Total isotropic sensitivity measurement in switched beam antenna systems | |
| CN111382587A (en) | Radio frequency reader-writer, test signal selection method and storage medium | |
| CN102148885A (en) | Method and system for testing multi-antenna terminal | |
| CN106851710A (en) | Base station testing system and method |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |