CN114966239A - Quasi-far field measuring method based on separable excitation coefficient variables - Google Patents

Quasi-far field measuring method based on separable excitation coefficient variables Download PDF

Info

Publication number
CN114966239A
CN114966239A CN202210902679.8A CN202210902679A CN114966239A CN 114966239 A CN114966239 A CN 114966239A CN 202210902679 A CN202210902679 A CN 202210902679A CN 114966239 A CN114966239 A CN 114966239A
Authority
CN
China
Prior art keywords
antenna
equation
far
quasi
field
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.)
Granted
Application number
CN202210902679.8A
Other languages
Chinese (zh)
Other versions
CN114966239B (en
Inventor
栗曦
张鹏
韩伟
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Shaanxi Shibei Communication Technology Co ltd
Original Assignee
Shaanxi Shibei Communication Technology Co ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Shaanxi Shibei Communication Technology Co ltd filed Critical Shaanxi Shibei Communication Technology Co ltd
Priority to CN202210902679.8A priority Critical patent/CN114966239B/en
Publication of CN114966239A publication Critical patent/CN114966239A/en
Application granted granted Critical
Publication of CN114966239B publication Critical patent/CN114966239B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R29/00—Arrangements for measuring or indicating electric quantities not covered by groups G01R19/00 - G01R27/00
    • G01R29/08—Measuring electromagnetic field characteristics
    • G01R29/10—Radiation diagrams of antennas
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R29/00—Arrangements for measuring or indicating electric quantities not covered by groups G01R19/00 - G01R27/00
    • G01R29/08—Measuring electromagnetic field characteristics
    • G01R29/0864—Measuring electromagnetic field characteristics characterised by constructional or functional features
    • G01R29/0892—Details related to signal analysis or treatment; presenting results, e.g. displays; measuring specific signal features other than field strength, e.g. polarisation, field modes, phase, envelope, maximum value
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10—Complex mathematical operations
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10—Complex mathematical operations
    • G06F17/14—Fourier, Walsh or analogous domain transformations, e.g. Laplace, Hilbert, Karhunen-Loeve, transforms
    • G06F17/141—Discrete Fourier transforms
    • G06F17/142—Fast Fourier transforms, e.g. using a Cooley-Tukey type algorithm
    • Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02D—CLIMATE CHANGE MITIGATION TECHNOLOGIES IN INFORMATION AND COMMUNICATION TECHNOLOGIES [ICT], I.E. INFORMATION AND COMMUNICATION TECHNOLOGIES AIMING AT THE REDUCTION OF THEIR OWN ENERGY USE
    • Y02D30/00—Reducing energy consumption in communication networks
    • Y02D30/70—Reducing energy consumption in communication networks in wireless communication networks

Landscapes

  • Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Mathematical Physics (AREA)
  • Data Mining & Analysis (AREA)
  • Theoretical Computer Science (AREA)
  • Computational Mathematics (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Pure & Applied Mathematics (AREA)
  • Algebra (AREA)
  • Databases & Information Systems (AREA)
  • Software Systems (AREA)
  • General Engineering & Computer Science (AREA)
  • Electromagnetism (AREA)
  • Discrete Mathematics (AREA)
  • Variable-Direction Aerials And Aerial Arrays (AREA)

Abstract

The invention discloses a quasi-far field measuring method based on separable excitation coefficient variables, which comprises the following steps: the method comprises the following steps: determining the quasi-far field distance of the antenna to be measured, and step two: and (3) extracting the amplitude and the phase of the quasi-far field position of the antenna, and performing the third step: obtaining a far field directional diagram function of the one-dimensional linear array, and the fourth step: obtaining a far field directional diagram function of the two-dimensional antenna array, and the fifth step: an amplitude and phase pattern and gain are obtained. The one-dimensional antenna array quasi-far field measurement method based on cylindrical wave expansion is suitable for measuring antennas with large one-dimensional electrical size and small another-dimensional electrical size, is suitable for measuring two-dimensional antenna arrays, increases the universality of a quasi-far field extrapolation algorithm, obtains an antenna far field directional diagram accurately, measures in a quasi-far field range of an antenna to be measured, and avoids the problems of large test distance and low near field measurement efficiency required by far field measurement, so that the antenna test efficiency is remarkably improved.

Description

Quasi-far field measuring method based on separable excitation coefficient variables
Technical Field
The invention belongs to the technical field of antenna measurement, and particularly relates to a quasi-far field measurement method based on separable excitation coefficient variables.
Background
The existing antenna far field measurement needs to meet the far field measurement distance, the conventional darkroom environment is difficult to meet the far field measurement requirement of a larger-caliber antenna, if a near field measurement method is adopted, the required time is long, the efficiency is low, in order to realize rapid measurement in a short distance, the prior art provides a one-dimensional antenna array quasi-far field measurement method based on cylindrical wave expansion, which is suitable for measuring the antenna with large one-dimensional electrical size and small other-dimensional electrical size, the method only needs to measure the distance to meet the far field condition of the antenna to be measured in the dimension with small size, therefore, the measurement distance is at the quasi-far-field position of the dimension with large size, and only information on one section can be obtained by single measurement, so the method is named as quasi-far-field measurement.
Disclosure of Invention
The invention aims to provide a quasi-far-field measuring method based on separable excitation coefficient variables, which solves the problems in the prior art.
In order to achieve the purpose, the invention provides the following technical scheme: the quasi-far field measurement method based on separable excitation coefficient variables comprises the following steps:
the method comprises the following steps: determining the quasi-far field distance of the antenna to be measured;
step two: extracting the amplitude and the phase at the quasi-far field position of the antenna;
step three: obtaining a far field directional diagram function of the one-dimensional linear array, which specifically comprises the following steps:
electric field of observation point located on plane under cylindrical coordinate system:
Figure 920370DEST_PATH_IMAGE001
equation 1
Wherein the content of the first and second substances,
Figure 371074DEST_PATH_IMAGE002
is shown at
Figure 769694DEST_PATH_IMAGE003
The electric field at the observation point on the plane,
Figure 291942DEST_PATH_IMAGE004
are all variables in a standard coordinate system under a spherical coordinate system,
Figure 944772DEST_PATH_IMAGE005
in order to be able to set the integer number,
Figure 831956DEST_PATH_IMAGE006
，
Figure 452293DEST_PATH_IMAGE007
the setting is 2-10, and the device is,
Figure 512653DEST_PATH_IMAGE008
to encompass the minimum cylinder radius of the antenna,jis the unit of an imaginary number,kas the number of free-space waves,
Figure 19989DEST_PATH_IMAGE009
in a standard coordinate system
Figure 999447DEST_PATH_IMAGE010
The direction unit vector of the direction unit vector,
Figure 869531DEST_PATH_IMAGE011
in a standard coordinate system
Figure 61478DEST_PATH_IMAGE012
The direction unit vector of the direction unit vector,
Figure 688899DEST_PATH_IMAGE013
and
Figure 370417DEST_PATH_IMAGE014
is the coefficient of expansion of the cylindrical wave,
Figure 184920DEST_PATH_IMAGE015
is composed of
Figure 180558DEST_PATH_IMAGE016
To the second type of Hankel function, at
Figure 193644DEST_PATH_IMAGE003
General form of the electric field of an observation point on a plane, its arbitrary linearly polarized electric field component
Figure 265637DEST_PATH_IMAGE017
Can be expressed as follows:
Figure 347862DEST_PATH_IMAGE018
equation 2
Wherein
Figure 960240DEST_PATH_IMAGE019
Represents the cylindrical wave expansion coefficient (of equation 1)A n OrB n ) Defined in the antenna under testρ = ρ 0 An electric field ofE m ，eTo the bottom of the exponent, the electric field at the observation point can be expressed as:
Figure 827833DEST_PATH_IMAGE020
equation 3
Then, inverting the above equation can result in:
Figure 601885DEST_PATH_IMAGE021
equation 4
Since the observation point is located in the far field, combining the properties of the Hankel function:
Figure 905828DEST_PATH_IMAGE022
equation 5
Bringing the above intoρ = ρ 0 In the expression of the electric field component, it can be found that:
Figure 259580DEST_PATH_IMAGE023
equation 6
Removing the constant irrelevant to the angle in the above formula to obtain the far field directional diagram function of the antenna to be measured
Figure 699788DEST_PATH_IMAGE024
Comprises the following steps:
Figure 441479DEST_PATH_IMAGE025
equation 7
Step four: obtaining a far-field directional pattern function of the two-dimensional antenna array, specifically:
the excitation coefficient of each cell can be expressed as:
Figure 717871DEST_PATH_IMAGE026
equation 8
Wherein the content of the first and second substances,Mis a value of the number of column-wise antenna elements,Nis a quantity value of a row-wise antenna element if the array antenna hasM×NA plurality of antenna units, each of which has a plurality of antenna elements,
Figure 265527DEST_PATH_IMAGE027
Figure 91401DEST_PATH_IMAGE028
Figure 941676DEST_PATH_IMAGE029
Figure 361156DEST_PATH_IMAGE030
the excitation coefficients are normalized for the columns of the two-dimensional antenna array,
Figure 571558DEST_PATH_IMAGE031
Figure 268249DEST_PATH_IMAGE032
Figure 679639DEST_PATH_IMAGE029
Figure 711049DEST_PATH_IMAGE033
for the row normalization excitation coefficient of the two-dimensional antenna array, a two-dimensional antenna array directional diagram can be obtained according to the directional diagram product theorem and the superposition theorem
Figure 600507DEST_PATH_IMAGE034
Comprises the following steps:
Figure 417285DEST_PATH_IMAGE035
equation 9
Wherein the content of the first and second substances,
Figure 468418DEST_PATH_IMAGE036
are variables in a standard coordinate system under spherical coordinates,
Figure 987124DEST_PATH_IMAGE037
the distance between the units in the direction isd x ，
Figure 290060DEST_PATH_IMAGE038
The distance between the direction units isd y ，I MN Is prepared from (a)M , N) The excitation coefficients of the location elements, m being the column-wise antenna element integer number, n being the row-wise antenna element integer number,
Figure 679453DEST_PATH_IMAGE039
for the excitation coefficients of the (m, n) -th element, the continued simplification can result in:
Figure 432645DEST_PATH_IMAGE040
equation 10
Then
Figure 189380DEST_PATH_IMAGE041
Equation 11
In the formula (I), the compound is shown in the specification,
Figure 420641DEST_PATH_IMAGE042
normalizing the excitation coefficient for the mth column of the two-dimensional antenna array;
Figure 930120DEST_PATH_IMAGE043
equation 12
In the formula (I), the compound is shown in the specification,
Figure 452282DEST_PATH_IMAGE044
normalizing the excitation coefficient for the nth row of the two-dimensional antenna array;
therefore, the temperature of the molten metal is controlled,
Figure 820946DEST_PATH_IMAGE045
can be expressed as
Figure 714953DEST_PATH_IMAGE046
Equation 13
Wherein the content of the first and second substances,
Figure 829671DEST_PATH_IMAGE047
and
Figure 924666DEST_PATH_IMAGE048
directional diagram functions of the column direction one-dimensional linear array and the row direction one-dimensional linear array are respectively calculated by using the formula in the step three;
step five: obtaining an amplitude directional diagram, a phase directional diagram and a gain, specifically:
amplitude directional diagram
Figure 639681DEST_PATH_IMAGE049
Comprises the following steps:
Figure 212744DEST_PATH_IMAGE050
equation 14
Wherein (A), (B), (C) and CdB) Is in amplitude units;
phase directional diagram
Figure 713127DEST_PATH_IMAGE051
Comprises the following steps:
Figure 838078DEST_PATH_IMAGE052
equation 15
Wherein angle () is a phase taking function inρ = ρ 0 An electric field ofE m The electric field can be expressed as:
Figure 650176DEST_PATH_IMAGE053
equation 16
At the same time, is located atρ = ρ 0 Electric field ofE m Can also be expressed as:
Figure 495772DEST_PATH_IMAGE054
equation 17
By the above two formulae, it is possible to obtain:
Figure 788345DEST_PATH_IMAGE055
equation 18
If the observation point is located in the far field region, the electric field can be written as:
Figure 553038DEST_PATH_IMAGE056
equation 19
The directional pattern function of the far field region
Figure 993378DEST_PATH_IMAGE057
Can be expressed as:
Figure 501720DEST_PATH_IMAGE058
equation 20
Gain compensation of antenna under test between quasi-far field and far field
Figure 711115DEST_PATH_IMAGE059
The calculation is as follows:
Figure 53235DEST_PATH_IMAGE060
equation 21
Wherein the content of the first and second substances,
Figure 230138DEST_PATH_IMAGE061
a far-field pattern function is represented,
Figure 558483DEST_PATH_IMAGE062
representing a quasi-far-field pattern function.
Preferably, in the first step, the quasi-far-field distance measurement satisfies the far-field distance of the antenna unit, the distance between the antenna to be measured and the test probe is calculated, and the antenna is placed on the turntable at the distance.
Preferably, in the second step, the control computer is used to control a test probe with a known characteristic and a working frequency consistent with that of the antenna to be tested, and the amplitude and phase on a certain surface at the quasi-far-field position of the antenna are corresponded, and the amplitude and phase information obtained by the test is stored in the test file of the control computer.
Compared with the prior art, the invention has the beneficial effects that:
the universality of a quasi far-field extrapolation algorithm is increased: the quasi-far field measurement method of the one-dimensional antenna array based on cylindrical wave expansion is suitable for measuring the antenna with large one-dimensional electrical size and small other-dimensional electrical size, is suitable for measuring the two-dimensional antenna array, increases the universality of a quasi-far field extrapolation algorithm, and finally obtains an antenna far field directional pattern accurately;
the antenna test efficiency is obviously improved: the antenna test system is used for measuring in the quasi-far field range of the antenna to be tested, so that the problems of large test distance and low near field measurement efficiency required by far field measurement are solved, and the antenna test efficiency is remarkably improved.
The measurement distance is significantly reduced: the distance between the antenna to be tested and the probe only needs to meet the quasi-far field condition, the testing distance is obviously reduced, and the cost of darkroom construction is greatly saved.
The reconstructed far-field directional diagram has high precision: by combining the quasi-far-field measurement method based on separable excitation coefficient variables, the finally obtained far-field directional pattern of the antenna to be measured is high in accuracy, and far-field data can be efficiently reconstructed through an algorithm.
Drawings
FIG. 1 is a schematic diagram of a quasi-far field antenna test of the present invention;
FIG. 2 is a schematic view of the (r, φ, z) coordinate system of the present invention;
FIG. 3 is a simplified model diagram of an array antenna of the present invention;
FIG. 4 is a simplified model diagram of an antenna under test according to the present invention;
FIG. 5 is a schematic flow chart of a quasi-far-field measurement method according to the present invention;
fig. 6 is a schematic diagram comparing a reconstructed amplitude pattern with a theoretically calculated amplitude pattern.
Detailed Description
The technical solution in the embodiments of the present invention will be clearly and completely described below with reference to the accompanying drawings in the embodiments of the present invention.
The invention discloses a quasi-far field measuring method based on separable excitation coefficient variables, which relates to an antenna to be tested arranged on a rotary table, a control computer, a test probe and a test instrument in an antenna measuring system in the measuring process, wherein the characteristics of all the test probes are known before testing, and the method comprises the following steps:
the method comprises the following steps: determining the quasi-far field distance of the antenna to be measured: the method is suitable for measuring a two-dimensional antenna array, the quasi-far-field distance can be met only by the far-field distance of the antenna unit, the distance between the antenna to be measured and the test probe is calculated, the antenna is placed on the turntable at the distance, referring to fig. 1, the distance between the antenna to be measured and the test probe placed on the turntable is the quasi-far-field distance of the antenna to be measured, and the antenna to be measured and the test probe are required to be consistent in working frequency, identical in height and matched in polarization before testing.
Step two: the amplitude and phase on a certain surface at the quasi-far field position of the antenna are extracted: and controlling a test probe with the same working frequency and known characteristics with the antenna to be tested by using a control computer, corresponding to the amplitude and the phase on the quasi-far-field surface of the antenna, and storing the amplitude and phase information in a test file of the control computer.
Step three: obtaining a far field directional diagram function of the one-dimensional linear array: according to the amplitude and the phase of a radiation field of an antenna to be detected, a one-dimensional antenna array quasi-far-field extrapolation algorithm based on cylindrical wave expansion firstly deduces the relation between a measured value and a coefficient on a certain tangent plane observation point, carries out fast Fourier transform calculation on the measured value to obtain the coefficient, combines the property of a Hankel function, and can extrapolate a far-field directional diagram function of a one-dimensional linear array by carrying out inverse Fourier transform on the coefficient, specifically:
referring to FIGS. 2 and 3, the cylindrical coordinate system is located at
Figure 12598DEST_PATH_IMAGE003
Electric field of observation point on plane:
Figure 650252DEST_PATH_IMAGE063
equation 1
Wherein the content of the first and second substances,
Figure 189818DEST_PATH_IMAGE002
representing the electric field at an observation point lying on a plane,
Figure 587433DEST_PATH_IMAGE004
are all variables in a standard coordinate system under a spherical coordinate system,
Figure 896054DEST_PATH_IMAGE005
in order to be able to set the integer number,
Figure 970189DEST_PATH_IMAGE006
，
Figure 731472DEST_PATH_IMAGE007
the setting is 2-10, and the device is,
Figure 932777DEST_PATH_IMAGE008
to encompass the minimum cylinder radius of the antenna,jis the unit of an imaginary number,kas the number of free-space waves,
Figure 954960DEST_PATH_IMAGE009
in a standard coordinate system
Figure 75363DEST_PATH_IMAGE010
The direction unit vector of the direction unit vector,
Figure 464887DEST_PATH_IMAGE011
in a standard coordinate system
Figure 63359DEST_PATH_IMAGE012
The direction unit vector of the direction unit vector,
Figure 940048DEST_PATH_IMAGE013
and
Figure 372297DEST_PATH_IMAGE014
is the coefficient of expansion of the cylindrical wave,
Figure 373751DEST_PATH_IMAGE015
is composed of
Figure 775914DEST_PATH_IMAGE016
To the second type of Hankel function, at
Figure 38268DEST_PATH_IMAGE003
General form of the electric field of an observation point on a plane, its arbitrary linearly polarized electric field component
Figure 766052DEST_PATH_IMAGE017
Can be expressed as follows:
Figure 599010DEST_PATH_IMAGE018
equation 2
Wherein
Figure 70443DEST_PATH_IMAGE019
Represents the cylindrical wave expansion coefficient (of equation 1)A n OrB n ) Defined in the antenna under testρ = ρ 0 An electric field ofE m ，eTo the bottom of the exponent, the electric field at the observation point can be expressed as:
Figure 187303DEST_PATH_IMAGE064
equation 3
By inverting the above equation, we can obtain:
Figure 820410DEST_PATH_IMAGE021
equation 4
Since the observation point is located in the far field, combining the properties of the Hankel function:
Figure 406243DEST_PATH_IMAGE022
equation 5
Bringing the above intoρ = ρ 0 In the expression of the electric field component, it can be found that:
Figure 415788DEST_PATH_IMAGE023
equation 6
Removing the constant irrelevant to the angle in the above formula to obtain the far field directional diagram function of the antenna to be measured
Figure 387155DEST_PATH_IMAGE024
Comprises the following steps:
Figure 456742DEST_PATH_IMAGE065
equation 7
Step four: obtaining a far field directional pattern function of the two-dimensional antenna array: according to the characteristic that the two-dimensional antenna array has separable variables of the excitation coefficient, a far field pattern function of the two-dimensional antenna array is obtained through calculation through a far field pattern function of the one-dimensional linear array, and specifically:
although the physical size of a common two-dimensional antenna array does not satisfy the conditions of large one-dimensional size and small one-dimensional size, the antenna has the characteristic of separable variables of excitation coefficients, namely the excitation coefficient of each unit can be expressed as:
Figure 264292DEST_PATH_IMAGE026
equation 8
Wherein the content of the first and second substances,Mis a value of the number of column-wise antenna elements,Nis a quantity value of a row-wise antenna element if the array antenna hasM×NAn antenna unit, wherein,I MN is prepared from (a)m , n) The excitation coefficient of the location unit is,
Figure 77527DEST_PATH_IMAGE027
Figure 168980DEST_PATH_IMAGE028
Figure 143889DEST_PATH_IMAGE029
Figure 704315DEST_PATH_IMAGE030
the excitation coefficients are normalized for the columns of the two-dimensional antenna array,
Figure 55662DEST_PATH_IMAGE031
Figure 1621DEST_PATH_IMAGE032
Figure 147431DEST_PATH_IMAGE029
Figure 941293DEST_PATH_IMAGE033
for the row normalization excitation coefficient of the two-dimensional antenna array, a two-dimensional antenna array directional diagram can be obtained according to the directional diagram product theorem and the superposition theorem
Figure 955385DEST_PATH_IMAGE034
Comprises the following steps:
Figure 896796DEST_PATH_IMAGE066
equation 9
Wherein the content of the first and second substances,
Figure 479087DEST_PATH_IMAGE036
are variables in a standard coordinate system under spherical coordinates,
Figure 748526DEST_PATH_IMAGE037
the distance between each unit in the direction isd x ，
Figure 441675DEST_PATH_IMAGE038
The distance between the direction units isd y ，I MN Is prepared from (a)M , N) The excitation coefficients of the location elements, m being the column-wise antenna element integer number, n being the row-wise antenna element integer number,
Figure 362227DEST_PATH_IMAGE039
for the excitation coefficients of the (m, n) -th element, the continued simplification can result in:
Figure 115419DEST_PATH_IMAGE040
equation 10
Then
Figure 606574DEST_PATH_IMAGE041
Equation 11
In the formula (I), the compound is shown in the specification,
Figure 103415DEST_PATH_IMAGE042
normalizing the excitation coefficient for the mth column of the two-dimensional antenna array;
Figure 612894DEST_PATH_IMAGE043
equation 12
In the formula (I), the compound is shown in the specification,
Figure 412353DEST_PATH_IMAGE044
normalizing the excitation coefficient for the nth row of the two-dimensional antenna array;
therefore, the temperature of the molten metal is controlled,
Figure 781018DEST_PATH_IMAGE045
can be expressed as
Figure 81549DEST_PATH_IMAGE046
Equation 13
Wherein the content of the first and second substances,
Figure 711114DEST_PATH_IMAGE047
and
Figure 806109DEST_PATH_IMAGE048
directional diagram functions of the column direction one-dimensional linear array and the row direction one-dimensional linear array are respectively calculated by using the formula in the step three;
step five: calculating to obtain an amplitude directional diagram, a phase directional diagram and gain of the antenna to be measured: according to the far field directional diagram function of the antenna to be measured, taking the absolute value of the far field directional diagram function of the antenna to be measured to obtain the amplitude directional diagram of the antenna to be measured, taking the angle of the far field directional diagram function of the antenna to be measured to obtain the phase directional diagram of the antenna to be measured, calculating to obtain the gain compensation of the antenna to be measured between a quasi far field and a far field, and completing the quasi far field measurement of the antenna, specifically:
amplitude directional diagram
Figure 271856DEST_PATH_IMAGE049
Comprises the following steps:
Figure 110499DEST_PATH_IMAGE050
equation 14
Wherein (A), (B), (C), (D), (C), (B), (C)dB) Is in amplitude units;
phase directional diagram
Figure 860149DEST_PATH_IMAGE051
Comprises the following steps:
Figure 126046DEST_PATH_IMAGE052
equation 15
Wherein angle () is a phase taking function inρ = ρ 0 Is under an electric field ofE m The electric field can be expressed as:
Figure 79089DEST_PATH_IMAGE067
equation 16
At the same time, is located atρ = ρ 0 Electric field ofE m Can also be expressed as:
Figure 455844DEST_PATH_IMAGE054
equation 17
By the above two formulae, it is possible to obtain:
Figure 60001DEST_PATH_IMAGE055
equation 18
If the observation point is located in the far field region, the electric field can be written as:
Figure 496798DEST_PATH_IMAGE056
equation 19
The directional pattern function of the far field region
Figure 937138DEST_PATH_IMAGE057
Can be expressed as:
Figure 117584DEST_PATH_IMAGE058
equation 20
Gain compensation of antenna to be tested between quasi-far field and far field
Figure 107405DEST_PATH_IMAGE059
The calculation is as follows:
Figure 715104DEST_PATH_IMAGE060
equation 21
Wherein the content of the first and second substances,
Figure 377161DEST_PATH_IMAGE061
a far-field pattern function is represented,
Figure 361297DEST_PATH_IMAGE062
representing a quasi-far-field pattern function.
The invention provides an overall technical scheme of a quasi-far-field measurement method based on separable excitation coefficient variables, which is suitable for measuring a two-dimensional antenna array and has the characteristic of separable excitation coefficient variables.
The quasi-far field measurement method based on separable excitation coefficient variables, which is provided by the invention, expands the antenna to be measured applicable to quasi-far field measurement from a one-dimensional antenna array to a two-dimensional antenna array, does not influence the precision and accuracy of reconstructing a far field pattern, and calculates the far field pattern of the antenna to be measured quickly, efficiently and accurately by a fast Fourier transform method.
The two-dimensional antenna array has the characteristic of separable excitation coefficient variables, is applied to the quasi-far-field measurement of an actual antenna, improves the universality of a quasi-far-field extrapolation algorithm, does not influence the precision of a far-field directional diagram, can obviously improve the testing efficiency of the quasi-far field of the antenna, and saves the cost of microwave darkroom construction.
The invention has the advantages that: the method improves the universality of the quasi-far-field extrapolation algorithm, is suitable for the quasi-far-field measurement of the two-dimensional antenna array, remarkably reduces the test distance of the antenna to be tested, does not influence the accuracy of a far-field directional pattern, and improves the efficiency of the quasi-far-field test.
The following is a simulation measurement based on the quasi-far-field measurement method with separable excitation coefficient variables disclosed by the invention:
simulation conditions and contents:
the antenna to be tested shown in fig. 4 is an array antenna, the frequency of the antenna to be tested is 3 GHz, the array antenna is arranged along two directions, each unit consists of two half-wave oscillators and consists of Ne antenna units, the distance between the y-axis antenna units is dy, the array amplitude distribution is Taylor distribution of-30 dB, and the detailed parameters of the antenna to be tested are shown in table 1.
Table 1: detailed parameters of the antenna to be measured
Parameter(s) Means of Value taking
Ne Number of array antenna elements 10*10
dy Unit pitch in y-axis direction 0.5λ
dx Cell pitch in x-axis direction 0.5λ
In a simulation experiment, the quasi-far-field test distance is selected to be 2 meters, the sampling interval is 1 degree, and the number of sampling points is 361 points.
Simulation results and analysis
Fig. 6 is a schematic diagram comparing the amplitude directional diagram reconstructed by the present invention with the amplitude directional diagram obtained by theoretical calculation, in which the abscissa represents angle change, the ordinate represents the amplitude value of the antenna at different angles, the solid line represents the amplitude directional diagram of the antenna obtained by theoretical calculation, and the dotted line represents the amplitude directional diagram reconstructed by the present invention, and refer to fig. 6 that the amplitude directional diagram reconstructed by the present invention matches with the amplitude directional diagram obtained by theoretical calculation.
In short, the quasi-far field measurement method based on separable excitation coefficient variables solves the problem of quasi-far field measurement of a two-dimensional antenna array, determines the test distance of an antenna to be tested according to the quasi-far field measurement condition of the antenna, extracts the amplitude and the phase on a certain surface at the quasi-far field position of the antenna, firstly calculates the far field directional diagram function of a one-dimensional linear array, obtains the far field directional diagram function of the two-dimensional antenna array based on the far field directional diagram function, realizes the quasi-far field measurement of the two-dimensional antenna array, provides the quasi-far field measurement mode of the two-dimensional antenna array by combining the characteristic that the two-dimensional antenna array has separable variables of the excitation coefficients, reconstructs the far field directional diagram of the antenna to be tested efficiently and accurately, and remarkably improves the test efficiency.

Claims (3)

1. The quasi-far-field measurement method based on separable excitation coefficient variables is characterized by comprising the following steps of:
the method comprises the following steps: determining the quasi-far field distance of the antenna to be measured;
step two: extracting the amplitude and the phase at the quasi-far field position of the antenna;
step three: obtaining a far field directional diagram function of the one-dimensional linear array, which specifically comprises the following steps:
electric field of observation point located on plane under cylindrical coordinate system:
Figure 3877DEST_PATH_IMAGE001
equation 1
Wherein, the first and the second end of the pipe are connected with each other,
Figure 907111DEST_PATH_IMAGE002
is shown lying on a plane
Figure 446677DEST_PATH_IMAGE003
The electric field of the observation point of (a),
Figure 434089DEST_PATH_IMAGE004
are all variables in a standard coordinate system under a spherical coordinate system,
Figure 742711DEST_PATH_IMAGE005
in order to be able to set the integer number,
Figure 957792DEST_PATH_IMAGE006
，
Figure 109287DEST_PATH_IMAGE007
the setting is 2-10, and the device is,
Figure 169647DEST_PATH_IMAGE008
to encompass the minimum cylinder radius of the antenna,jis the unit of an imaginary number,kas the number of free-space waves,
Figure 208142DEST_PATH_IMAGE009
in a standard coordinate system
Figure 594123DEST_PATH_IMAGE010
The direction unit vector of the direction unit vector,
Figure 967336DEST_PATH_IMAGE011
in a standard coordinate system
Figure 831387DEST_PATH_IMAGE012
The direction unit vector of the direction unit vector,
Figure 255546DEST_PATH_IMAGE013
and
Figure 546850DEST_PATH_IMAGE014
is the coefficient of expansion of the cylindrical wave,
Figure 407359DEST_PATH_IMAGE015
is composed of
Figure 75100DEST_PATH_IMAGE016
To the second type of Hankel function, at
Figure 88187DEST_PATH_IMAGE003
General form of the electric field of an observation point on a plane, its arbitrary linearly polarized electric field component
Figure 550392DEST_PATH_IMAGE017
Can be expressed as follows:
Figure 39142DEST_PATH_IMAGE018
equation 2
Wherein
Figure 369630DEST_PATH_IMAGE019
Represents the cylindrical wave expansion coefficient (in equation 1)A n OrB n ) Defined in the antenna under testρ = ρ 0 An electric field ofE m ，eTo the bottom of the exponent, the electric field at the observation point can be expressed as:
Figure 361856DEST_PATH_IMAGE020
equation 3
Then, inverting the above equation can result in:
Figure 135908DEST_PATH_IMAGE021
equation 4
Since the observation point is located in the far field, combining the properties of the Hankel function:
Figure 846376DEST_PATH_IMAGE022
equation 5
Bringing the above intoρ = ρ 0 In the expression of the electric field component, it can be found that:
Figure 714974DEST_PATH_IMAGE023
equation 6
Removing the constant irrelevant to the angle in the above formula to obtain the far field directional diagram function of the antenna to be measured
Figure 827287DEST_PATH_IMAGE024
Comprises the following steps:
Figure 506661DEST_PATH_IMAGE025
equation 7
Step four: obtaining a far-field directional pattern function of the two-dimensional antenna array, specifically:
the excitation coefficient for each cell can be expressed as:
Figure 704424DEST_PATH_IMAGE026
equation 8
Wherein the content of the first and second substances,Mis a value of the number of column-wise antenna elements,Nis a quantity value of a row-wise antenna element if the array antenna hasM×NA plurality of antenna units, each of which has a plurality of antenna elements,
Figure 376714DEST_PATH_IMAGE027
Figure 609112DEST_PATH_IMAGE028
Figure 724967DEST_PATH_IMAGE029
Figure 144447DEST_PATH_IMAGE030
the excitation coefficients are normalized for the columns of the two-dimensional antenna array,
Figure 620428DEST_PATH_IMAGE031
Figure 441753DEST_PATH_IMAGE032
Figure 728509DEST_PATH_IMAGE029
Figure 900864DEST_PATH_IMAGE033
for the row normalization excitation coefficient of the two-dimensional antenna array, a two-dimensional antenna array directional diagram can be obtained according to the directional diagram product theorem and the superposition theorem
Figure 790323DEST_PATH_IMAGE034
Comprises the following steps:
Figure 856368DEST_PATH_IMAGE035
equation 9
Wherein the content of the first and second substances,
Figure 438659DEST_PATH_IMAGE036
are variables in a standard coordinate system under spherical coordinates,
Figure 708098DEST_PATH_IMAGE037
the distance between the units in the direction isd x ，
Figure 666826DEST_PATH_IMAGE038
The distance between the direction units isd y ，I MN Is prepared from (a)M , N) The excitation coefficients of the location elements, m being the column-wise antenna element integer number, n being the row-wise antenna element integer number,
Figure 321799DEST_PATH_IMAGE039
for the excitation coefficients of the (m, n) -th element, continued simplification can result in:
Figure 74991DEST_PATH_IMAGE040
equation 10
Then
Figure 831725DEST_PATH_IMAGE041
Equation 11
In the formula (I), the compound is shown in the specification,
Figure 328566DEST_PATH_IMAGE042
normalizing the excitation coefficient for the mth column of the two-dimensional antenna array;
Figure 103624DEST_PATH_IMAGE043
equation 12
In the formula (I), the compound is shown in the specification,
Figure 762138DEST_PATH_IMAGE044
normalizing the excitation coefficient for the nth row of the two-dimensional antenna array;
therefore, the temperature of the molten metal is controlled,
Figure 396382DEST_PATH_IMAGE045
can be expressed as
Figure 306700DEST_PATH_IMAGE046
Equation 13
Wherein the content of the first and second substances,
Figure 77210DEST_PATH_IMAGE047
and
Figure 31260DEST_PATH_IMAGE048
directional diagram functions of the column direction one-dimensional linear array and the row direction one-dimensional linear array are respectively calculated by using the formula in the step three;
step five: obtaining an amplitude directional diagram, a phase directional diagram and a gain, specifically:
amplitude directional diagram
Figure 887220DEST_PATH_IMAGE049
Comprises the following steps:
Figure 866809DEST_PATH_IMAGE050
equation 14
Wherein (A), (B), (C), (D), (C), (B), (C)dB) Is in amplitude units;
phase direction diagram function
Figure 226246DEST_PATH_IMAGE051
Comprises the following steps:
Figure 616776DEST_PATH_IMAGE052
equation 15
Wherein angle () is a phase taking function inρ = ρ 0 An electric field ofE m The electric field can be expressed as:
Figure 960033DEST_PATH_IMAGE053
equation 16
At the same time, is located atρ = ρ 0 Electric field ofE m Can also be expressed as:
Figure 602366DEST_PATH_IMAGE054
equation 17
By the above two formulae, it is possible to obtain:
Figure 957256DEST_PATH_IMAGE055
equation 18
If the observation point is located in the far field region, the electric field can be written as:
Figure 394053DEST_PATH_IMAGE056
equation 19
The directional pattern function of the far field region
Figure 83660DEST_PATH_IMAGE057
Can be expressed as:
Figure 529685DEST_PATH_IMAGE058
equation 20
Gain compensation of antenna under test between quasi-far field and far field
Figure 739081DEST_PATH_IMAGE059
The calculation is as follows:
Figure 346780DEST_PATH_IMAGE060
equation 21
Wherein the content of the first and second substances,
Figure 523683DEST_PATH_IMAGE061
a far-field pattern function is represented,
Figure 242240DEST_PATH_IMAGE062
representing a quasi-far-field pattern function.
2. The quasi-far-field measurement method based on separable excitation coefficient variables according to claim 1, characterized in that: in the first step, the quasi-far-field distance measurement meets the far-field distance of the antenna unit, the distance between the antenna to be measured and the test probe is calculated, and the antenna is placed on the turntable at the distance.
3. The quasi-far-field measurement method based on separable excitation coefficient variables according to claim 1, characterized in that: in the second step, a control computer is used for controlling a test probe which is consistent with the working frequency of the antenna to be tested and has known characteristics, the amplitude and the phase position of the test probe on a certain surface at the quasi-far-field position of the antenna are corresponded, and the amplitude and the phase position information obtained by the test are stored in a test file of the control computer.
CN202210902679.8A 2022-07-29 2022-07-29 Quasi-far field measuring method based on separable excitation coefficient variables Active CN114966239B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202210902679.8A CN114966239B (en) 2022-07-29 2022-07-29 Quasi-far field measuring method based on separable excitation coefficient variables

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202210902679.8A CN114966239B (en) 2022-07-29 2022-07-29 Quasi-far field measuring method based on separable excitation coefficient variables

Publications (2)

Publication Number Publication Date
CN114966239A true CN114966239A (en) 2022-08-30
CN114966239B CN114966239B (en) 2022-10-28

Family

ID=82969584

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202210902679.8A Active CN114966239B (en) 2022-07-29 2022-07-29 Quasi-far field measuring method based on separable excitation coefficient variables

Country Status (1)

Country Link
CN (1) CN114966239B (en)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115524541A (en) * 2022-11-28 2022-12-27 陕西拾贝通讯技术有限公司 Arbitrary polarization line antenna array directional diagram measuring method based on quasi-far-field mode filtering
CN116819431A (en) * 2023-08-31 2023-09-29 杭州岸达科技有限公司 Phase interferometer direction finding method based on anomalous phase mode excitation
WO2024045649A1 (en) * 2022-08-31 2024-03-07 东南大学 Near-field air interface rapid measurement system and method for pattern of active array antenna

Citations (19)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4270129A (en) * 1979-01-30 1981-05-26 Sperry Corporation Apparatus and method for realizing preselected free space antenna patterns
US6456235B1 (en) * 2001-03-29 2002-09-24 Raytheon Company Method of predicting the far field pattern of a slotted planar array at extreme angles using planar near field data
US6469672B1 (en) * 2001-03-15 2002-10-22 Agence Spatiale Europeenne (An Inter-Governmental Organization) Method and system for time domain antenna holography
JP2003028905A (en) * 2001-07-11 2003-01-29 Matsushita Electric Ind Co Ltd Analyzing device for antenna excitation coefficient
CN102590656A (en) * 2012-01-03 2012-07-18 西安电子科技大学 Antenna cover electric property forecasting method based on distant field
US20140348255A1 (en) * 2013-05-21 2014-11-27 Telefonica, S.A. Method and system for performing multiple access in wireless ofdm cellular systems considering both space and frequency domains
US9053268B1 (en) * 2007-12-19 2015-06-09 The United States Of America As Represented By The Secretary Of The Navy Analytic antenna design for a dipole antenna
CN106291129A (en) * 2015-06-01 2017-01-04 北京空间飞行器总体设计部 Phased array antenna far-field pattern measuring method
US20170155193A1 (en) * 2015-11-30 2017-06-01 Elwha Llc Beam pattern projection for metamaterial antennas
WO2018023929A1 (en) * 2016-08-01 2018-02-08 深圳市新益技术有限公司 Integrated antenna test system
CA2946011A1 (en) * 2016-10-21 2018-04-21 C-Com Satellite Systems Inc. Method and apparatus for phased antenna array calibration
CN108594026A (en) * 2018-01-31 2018-09-28 清华大学 A kind of near-field test method for rebuilding antenna array far-field pattern in radiading near-field region
US20200028598A1 (en) * 2018-05-28 2020-01-23 Keysight Technologies, Inc. Method and system for testing antenna array using middle field antenna pattern
CN110764068A (en) * 2019-10-25 2020-02-07 上海霍莱沃电子系统技术股份有限公司 Multi-probe quasi-far-field electromagnetic scattering cross section (RCS) extrapolation test system
US20200213016A1 (en) * 2018-12-26 2020-07-02 Keysight Technologies, Inc. System and method for obtaining far field radiated power with multiple radiated power measurements in middle field range
TWI727848B (en) * 2020-07-06 2021-05-11 川升股份有限公司 Quasi-far-field measurement system quasi-far-field measurement method
CN113533867A (en) * 2021-07-14 2021-10-22 西安电子科技大学 Fourier interpolation-based far field pattern rapid measurement method
EP3952023A1 (en) * 2020-08-03 2022-02-09 Rohde & Schwarz GmbH & Co. KG Method for manufacturing an antenna array, antenna array, and test system
WO2022088213A1 (en) * 2020-11-02 2022-05-05 网络通信与安全紫金山实验室 Phase difference measurement method and system and phase compensation method and system for antenna array

Patent Citations (19)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4270129A (en) * 1979-01-30 1981-05-26 Sperry Corporation Apparatus and method for realizing preselected free space antenna patterns
US6469672B1 (en) * 2001-03-15 2002-10-22 Agence Spatiale Europeenne (An Inter-Governmental Organization) Method and system for time domain antenna holography
US6456235B1 (en) * 2001-03-29 2002-09-24 Raytheon Company Method of predicting the far field pattern of a slotted planar array at extreme angles using planar near field data
JP2003028905A (en) * 2001-07-11 2003-01-29 Matsushita Electric Ind Co Ltd Analyzing device for antenna excitation coefficient
US9053268B1 (en) * 2007-12-19 2015-06-09 The United States Of America As Represented By The Secretary Of The Navy Analytic antenna design for a dipole antenna
CN102590656A (en) * 2012-01-03 2012-07-18 西安电子科技大学 Antenna cover electric property forecasting method based on distant field
US20140348255A1 (en) * 2013-05-21 2014-11-27 Telefonica, S.A. Method and system for performing multiple access in wireless ofdm cellular systems considering both space and frequency domains
CN106291129A (en) * 2015-06-01 2017-01-04 北京空间飞行器总体设计部 Phased array antenna far-field pattern measuring method
US20170155193A1 (en) * 2015-11-30 2017-06-01 Elwha Llc Beam pattern projection for metamaterial antennas
WO2018023929A1 (en) * 2016-08-01 2018-02-08 深圳市新益技术有限公司 Integrated antenna test system
CA2946011A1 (en) * 2016-10-21 2018-04-21 C-Com Satellite Systems Inc. Method and apparatus for phased antenna array calibration
CN108594026A (en) * 2018-01-31 2018-09-28 清华大学 A kind of near-field test method for rebuilding antenna array far-field pattern in radiading near-field region
US20200028598A1 (en) * 2018-05-28 2020-01-23 Keysight Technologies, Inc. Method and system for testing antenna array using middle field antenna pattern
US20200213016A1 (en) * 2018-12-26 2020-07-02 Keysight Technologies, Inc. System and method for obtaining far field radiated power with multiple radiated power measurements in middle field range
CN110764068A (en) * 2019-10-25 2020-02-07 上海霍莱沃电子系统技术股份有限公司 Multi-probe quasi-far-field electromagnetic scattering cross section (RCS) extrapolation test system
TWI727848B (en) * 2020-07-06 2021-05-11 川升股份有限公司 Quasi-far-field measurement system quasi-far-field measurement method
EP3952023A1 (en) * 2020-08-03 2022-02-09 Rohde & Schwarz GmbH & Co. KG Method for manufacturing an antenna array, antenna array, and test system
WO2022088213A1 (en) * 2020-11-02 2022-05-05 网络通信与安全紫金山实验室 Phase difference measurement method and system and phase compensation method and system for antenna array
CN113533867A (en) * 2021-07-14 2021-10-22 西安电子科技大学 Fourier interpolation-based far field pattern rapid measurement method

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
XI LI: "Fast Determination of Single-Cut Far-Field Pattern of Base Station Antenna at a Quasi-Far-Field Distance", 《IEEE TRANSACTIONS ON ANTENNAS AND PROPAGATION》 *
可阳: "基于傅里叶级数展开的柱面近场测量算法", 《2021年全国天线年会论文集》 *
邝浩欣: "准远场天线测量修正方法研究", 《宇航计测技术》 *

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2024045649A1 (en) * 2022-08-31 2024-03-07 东南大学 Near-field air interface rapid measurement system and method for pattern of active array antenna
CN115524541A (en) * 2022-11-28 2022-12-27 陕西拾贝通讯技术有限公司 Arbitrary polarization line antenna array directional diagram measuring method based on quasi-far-field mode filtering
CN115524541B (en) * 2022-11-28 2023-03-10 陕西拾贝通讯技术有限公司 Arbitrary polarization line antenna array directional diagram measuring method based on quasi-far-field mode filtering
CN116819431A (en) * 2023-08-31 2023-09-29 杭州岸达科技有限公司 Phase interferometer direction finding method based on anomalous phase mode excitation
CN116819431B (en) * 2023-08-31 2023-12-08 杭州岸达科技有限公司 Phase interferometer direction finding method based on anomalous phase mode excitation

Also Published As

Publication number Publication date
CN114966239B (en) 2022-10-28

Similar Documents

Publication Publication Date Title
CN114966239B (en) Quasi-far field measuring method based on separable excitation coefficient variables
Persson et al. Reconstruction of equivalent currents using a near-field data transformation-with radome applications
Joy et al. Applications of probe-compensated near-field measurements
Bucci et al. Accurate diagnosis of conformal arrays from near-field data using the matrix method
CN106291130B (en) A kind of near field antenna measurements method of arbitrary surface scanning
CN113533867B (en) Fourier interpolation-based far field pattern rapid measurement method
CN107783086B (en) Method for diagnosing distorted position of antenna array aperture amplitude phase field
CN105445566A (en) Spherical multi-probe antenna test data processing method
CN111553046B (en) Antenna radiation calculation method based on spherical wave expansion and source reconstruction
CN108959806B (en) Equivalent radiation modeling method based on spherical surface near-field measurement and spherical mode source
Varela et al. Fast spherical near-field to far-field transformation for offset-mounted antenna measurements
CN109286069A (en) Array antenna field source reconstructing method and device based on dipole
Hirose et al. Equivalence between Kim method and Omi method in single-cut near-to-far-field transformations for antenna near-field measurements
CN115524541B (en) Arbitrary polarization line antenna array directional diagram measuring method based on quasi-far-field mode filtering
CN111460900B (en) Complex electromagnetic environment construction equivalence quantitative evaluation method
Khashimov et al. Dynamic patterns technique for circular phased array diagnostics
Varela et al. Single-cut phaseless near-field measurements using specialized probes
Liu et al. An automatic antenna near-field measurement system for narrow beam array antennas
Liu et al. A full-polarization scattering matrix measurement method based on antenna spatial polarization characteristic analysis
Liu et al. Extraction of equivalent dipole-moment model based on improve mapped matrix
Omi et al. Subdivision technique in near-field far-field transformation for 2-D radar cross section measurement for large objects
Miloserdov et al. A novel approach to measurement and calibration of space-borne SAR active phased array
Wang et al. Partial discharge ultrasonic localization simulation based on double helix ultrasonic sensor array
WO2023050407A1 (en) Electromagnetic field near-field test method and system, and readable storage medium and computer device
Di Massa et al. Large Conformal Array Diagnostics by Fast Matrix Method

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant