CN102904015A - Short-wave small circular receiving antenna array - Google Patents

Short-wave small circular receiving antenna array Download PDF

Info

Publication number
CN102904015A
CN102904015A CN2012104426714A CN201210442671A CN102904015A CN 102904015 A CN102904015 A CN 102904015A CN 2012104426714 A CN2012104426714 A CN 2012104426714A CN 201210442671 A CN201210442671 A CN 201210442671A CN 102904015 A CN102904015 A CN 102904015A
Authority
CN
China
Prior art keywords
array
antenna
monopole antenna
directional diagram
radius
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
CN2012104426714A
Other languages
Chinese (zh)
Other versions
CN102904015B (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.)
Wuhan University WHU
Original Assignee
Wuhan University WHU
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 Wuhan University WHU filed Critical Wuhan University WHU
Priority to CN201210442671.4A priority Critical patent/CN102904015B/en
Publication of CN102904015A publication Critical patent/CN102904015A/en
Application granted granted Critical
Publication of CN102904015B publication Critical patent/CN102904015B/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Variable-Direction Aerials And Aerial Arrays (AREA)

Abstract

The invention relates to a short-wave small circular receiving antenna array. The short-wave small circular receiving antenna array comprises a small uniform circular array which consists of monopole antennae, a matching network, a radar radio frequency receiving channel, a weighted generator, a cross-coupling corrector and a weighting device, wherein the radius of the circular array is less than or equal to a quarter of wavelength; the quantity of antennae and the radius of the circular array are determined according to the characteristic of a directional diagram which is required to be formed; the matching network is in impedance matching with the receiving channel; the weighted generator generates a weight according to the distribution condition of the circular array and the characteristic of the directional diagram; the cross-coupling corrector forms correction data by measuring the cross-coupling condition of the antennae to correct the generated weight; and the weighting device weights an antenna unit. On the premise that signal to noise ratio and angular resolution are ensured, occupation area is greatly decreased, a portable high-frequency ground wave radar antenna array is realized, and sea surface kinetic parameters can be effectively detected.

Description

A kind of shortwave miniature circular receiving antenna array
Technical field
The invention belongs to antenna technical field, particularly the high-frequency ground wave radar receiving antenna array.
Background technology
High-frequency ground wave radar is a kind of important detecting devices that can carry out the continuous large-area marine environmental monitoring, main drive marine mathematic(al) parameter and the sea low-altitude low-speed moving targets such as ocean surface wind, wave, stream, tide surveyed are of great significance national economy and national defense construction.In order to obtain narrower antenna main lobe width, namely in order to obtain higher angular resolution, the bore of receiving antenna array is often very large, array length reaches tens meters even upper km, therefore need to take very large " Gold Coast " area, greatly increase construction difficulty and construction cost, seriously hindered the promotion and application of high-frequency ground wave radar.
Summary of the invention
For the large problem of traditional high-frequency ground wave radar antenna aperture, the present invention proposes a kind of shortwave miniature circular receiving antenna array, this antenna array can obtain directivity factor and the main lobe width identical with traditional large aperture antenna battle array in the situation that guarantee certain signal to noise ratio, generate required directional diagram, and be easy to realize 360 ° of omnidirectional's scannings.
Technical scheme of the present invention is a kind of shortwave miniature circular receiving antenna array, comprises the little spacing monopole antenna of circle battle array, matching network, radar radio frequency reception channel, weighting generator, mutual coupling calibration device and a weighter;
Circular little spacing monopole antenna battle array comprises a plurality of monopole antennas, and each monopole antenna is evenly arranged on the circumference, and the radius of circular little spacing monopole antenna battle array is less than or equal to 1/4th of radar operation wavelength;
The output of each monopole antenna connects the input of corresponding matching network, and the output of each matching network is through corresponding radar radio frequency reception channel difference weighted input device;
The weights that the weighting generator produces output to the mutual coupling calibration device; The mutual coupling calibration device forms correction data according to the mutual coupling situation between each monopole antenna, proofread and correct weights and weighted input device that the weighting generator produces, weighter is weighted each monopole antenna according to the input by weights and corresponding radar radio frequency reception channel after the correction of mutual coupling calibration device input, form specific directional diagram, support 360 ° of omnidirectional's scannings.
And matching network is five rank broadband matching networks.
And the radius of monopole antenna number and circular little spacing monopole antenna battle array is determined according to the characteristics of the default directional diagram that will form; The characteristics of described directional diagram comprise main lobe direction, main lobe width and sidelobe level inhibiting value.
And, the weighting generator is according to the characteristics of the radius of monopole antenna number, circular little spacing monopole antenna battle array and the directional diagram that will form, utilize maximization directivity factor method under the secondary lobe constraint to calculate the weighted value of each antenna element, computational process comprises following substep
Step 1, input parameter comprises radius r, the operating frequency f of monopole antenna number N, circular little spacing monopole antenna battle array 0, main lobe direction
Figure BDA00002368068300021
With secondary lobe constraints Wherein, θ sWith
Figure BDA00002368068300023
Represent respectively the angle of pitch and azimuth, secondary lobe constraints
Figure BDA00002368068300024
Represent the amplitude that i secondary lobe will retrain;
Step 2 is according to array number N, circle battle array radius r, operating frequency f 0And main lobe direction
Figure BDA00002368068300025
Obtain one group of weight w with maximization directivity factor method D
Step 3 is with weight w DWeighing vector calculated direction figure as array;
Step 4, whether the sidelobe level of directional diagram all is less than or equal to the sidelobe level inhibiting value in the directional diagram characteristics in the determining step 3
Figure BDA00002368068300026
Then to organize weights output, process ends, otherwise enter step 5;
Step 5, find out backlog demand secondary lobe the position and add secondary lobe constraints;
Step 6 is obtained weights with maximization directivity factor method under current secondary lobe constraints, obtain new w D, return step 3 and calculate, until sidelobe level all meets the demands, export current weight w D
And the mutual coupling calibration device forms correction data according to the mutual coupling situation between each monopole antenna, and implementation is proofreaied and correct the weights that the weighting generator produces as correction data with mutual resistance matrix for generate the mutual resistance matrix of array with the S parameter; Described S parameter be antenna with adopt before matching network is connected the two-port network analyzer directly in actual working environment, to measure.
The invention has the beneficial effects as follows, under the prerequisite that guarantees certain signal to noise ratio and angular resolution, greatly reduced floor space, realized the portable high frequency groundwave radar antenna array, can effectively survey the ocean surface kinetic parameter.
Description of drawings
Fig. 1 is the structural representation of the embodiment of the invention.
Fig. 2 is the matching network structural representation of the embodiment of the invention.
Fig. 3 is 8 yuan of embodiment of the invention circle battle array main lobes array element weighting schematic diagrames when pointing to 15 °.
Fig. 4 is array element weighting schematic diagram when main lobe points to 60 ° after the counterclockwise by turns weighting of the round battle array of 8 yuan of the embodiment of the invention.
Fig. 5 is the lower maximization directivity factor of the sidelobe level constraint synthesis calculation flow chart of the embodiment of the invention.
Fig. 6 be the embodiment of the invention measure S parameter schematic diagram between two antennas with network analyzer.
8 yuan of radius 2.5m circle battle array main lobes pointed to 90 ° directional diagram when Fig. 7 was the frequency 10MHz of the embodiment of the invention.
8 yuan of radius 2.5m circle battle array Sidelobe Suppression values were that 10dB, main lobe point to 90 ° directional diagram when Fig. 8 was the frequency 10MHz of the embodiment of the invention.
Embodiment
Describe technical solution of the present invention in detail below in conjunction with drawings and Examples.
Referring to Fig. 1, the shortwave miniature circular receiving antenna array that the embodiment of the invention provides comprises the little spacing monopole antenna of circle battle array, matching network, radar radio frequency reception channel, weighting generator, mutual coupling calibration device and a weighter.Circular little spacing monopole antenna battle array comprises a plurality of monopole antennas, and each monopole antenna is evenly arranged on the circumference, and the radius of circular little spacing monopole antenna battle array is less than or equal to 1/4th of radar operation wavelength.For ease of outlines device, below shortwave miniature circular receiving antenna array is called for short the circle battle array.Matching network is used for and the radar radio frequency reception channel is carried out impedance matching; The weighting generator produces weights according to circle battle array distribution situation and directional diagram characteristics, and the weights that produce output to the mutual coupling calibration device; The mutual coupling calibration device forms correction data according to the mutual coupling situation between each monopole antenna, proofreaies and correct weights and weighted input device that the weighting generator produces; Weighter is weighted each monopole antenna according to the input by weights and corresponding radar radio frequency reception channel after the correction of mutual coupling calibration device input, forms specific directional diagram, supports 360 ° of omnidirectional's scannings.
Each monopole antenna is established total N monopole antenna D1, D2, D3, D4 as an antenna element ... DN is arranged on the circumference equably, and circle battle array radius is less than or equal to 1/4th of radar operation wavelength, and is propped up by support.Circle battle array radius is less than or equal to the spacing requirement that 1/4th of wavelength belongs to the Superdirective antenna array, therefore can obtain the directivity factor identical with the wide aperture array antenna.The height of monopole antenna and diameter can require to select according to gain size, weight requirement and wind resistance.
Antenna number and a circle battle array radius determined by the characteristics of the directional diagram that will form, for example main lobe direction, main lobe width and sidelobe level inhibiting value, and those skilled in the art can preset the characteristics of the directional diagram that will form according to demand.Can require to determine antenna main lobe direction and Sidelobe Suppression value according to radar performance first, then in the situation that can reach required main lobe width, optimize fate number and circle battle array radius with existing maximization directivity factor method with Computer Simulation, maximization directivity factor method is in following derivation and the step that detailed description is arranged.Namely elder generation is according to the requirement of antenna array miniaturization, select the antenna of right quantity and suitable circle battle array radius, and circle battle array radius is less than or equal to 1/4th of radar operation wavelength, in computer, obtain the directional diagram of antenna with maximization directivity factor method, see that whether the main lobe width of directional diagram is less than required main lobe width, if less than namely meeting the demands, otherwise need to increase the antenna number or increase circle battle array radius (but can not surpass radar operation wavelength 1/4th).Theoretical proof suitably increases the antenna number and can reduce main lobe width.Until main lobe width meets the demands.
Among the embodiment, the circle battle array is as the receiving antenna array of high-frequency ground wave radar, and operating frequency is 7.5MHz-25MHz, the height of each antenna element is 2m, adopts broadband omni directional monopoles sub antenna, in 7.5MHz-25MHz frequency band, its standing-wave ratio is less than 2, and variation is mild.Before a burst of type of circle was determined, embodiment optimized number of antennas and circle battle array radius according to the characteristics of directional diagram.In high-frequency ground wave radar, when operating frequency during at 7.5MHz-10MHz, require main lobe width on each scanning direction less than 35 °.Find when comprehensive, circle battle array radius is during less than 3m, can allow the main lobe width be 8 less than 35 ° of required minimum array numbers, so adopt 8 antennas, a circle battle array radius can be according to the degree of miniaturization and the variation comprehensive selection of main lobe width.
Because at high band, external noise is greater than the internal system noise, therefore when aerial array dwindles, antenna reception to signal and noise reduce in proportion, namely signal to noise ratio is constant, this is the theoretical foundation of high frequency antenna battle array miniaturization.The weighting scheme of circle battle array is fairly simple, because each antenna element is evenly distributed, the directional diagram that points to all directions is period profile, and therefore a demand goes out the weighted value of required antenna pattern in the one-period, then array element is rotated weighting, can form the directional diagram that points to other directions.For this reason, embodiment provides weighting generator and mutual coupling calibration device.
Described weighting generator namely according to antenna number, circle battle array radius and directional diagram characteristics, utilizes the superdirectivity synthesis to calculate the weighted value of each antenna element.Embodiment is according to formation distribution, main lobe direction and sidelobe level inhibiting value, and the maximization directivity factor synthesis under the constraint of employing secondary lobe is obtained the weighted value of each antenna.During implementation, can adopt computer software technology to realize the automatic operation of maximization directivity factor synthesis, namely the weighting generator can adopt computer or chip microcontroller, perhaps adopts the realizations such as software solidification module, integrated circuit.
Described mutual coupling calibration device, generate the mutual resistance matrix of array with the S parameter, weighted value is proofreaied and correct as correction data with mutual resistance matrix, weighter is then with the output multiplication of the weighted value after proofreading and correct to the respective radio-frequency receive path, realization is weighted each antenna element, by by turns weighting, realize 360 ° of omnidirectional's scannings.During implementation, can adopt computer software technology to generate mutual resistance matrix, the weighted calculation of array, namely the weighting generator can adopt computer or chip microcontroller, perhaps adopts the realizations such as software solidification module, integrated circuit.Weighter can adopt logical integrated circuit, for example adopts adder and multiplier to realize weighting.Described S parameter, can antenna with adopt before matching network is connected the two-port network analyzer directly in actual working environment, to measure.Because each antenna distance is less, the mutual coupling degree between antenna is very large, therefore must carry out mutual coupling calibration, to eliminate the interference between antenna.After antenna number and circle battle array radius are determined, be after the formation distribution is determined, the circle battle array is placed on the work-yard, measure S parameter between each antenna with network analyzer, again the S parameter is converted to mutual resistance matrix, according to the input impedance value of radar radio frequency reception channel, mutual resistance matrix is carried out normalization, proofread and correct the weights that the weighting generator generates with the normalization mutual resistance matrix that obtains.
Matching network is used for and receive path carries out impedance matching.The output of each monopole antenna connects the input of corresponding matching network, and the output of each matching network is connected with corresponding radar radio frequency reception channel; Each radar radio frequency reception channel is distinguished the weighted input device again.Described radar radio frequency reception channel can be various radar receivers, has expressed the annexation of signal, and its internal structure is not at the row of content of the present invention.The output of each antenna is matching connection network 1,2 respectively ... N is used for and the radar radio frequency reception channel input impedance of rear end 50 Ω is complementary.N receive path 1,2 ... N respectively N output with the mutual coupling calibration device is corresponding, and the mutual coupling calibration device links to each other with the weighting generator.Matching network can have multiple design, and scheme shown in Figure 2 is five rank matching networks, is made of capacitor C 1=52.9pF wherein, C2=97.1pF, C3=51.6pF, inductance L 1=689.5nH, L2=1.7uH three electric capacity and two inductance.Series capacitance C1 and inductance L 1 between port P1, the P2; Capacitor C 2 one ends are connected between capacitor C 1 and the inductance L 1, other end ground connection; Capacitor C 3 and inductance L 2 homogeneous end connectivity port P2, other end ground connection.Capacitance and inductance value are only in this scheme.Port P1 connects the feedback point of monopole antenna, and port P2 connects receive path.During implementation, the invention technician can arrange matching network voluntarily as the case may be.
For the purpose of efficient is provided, can based on asking for the weights in the scan period, obtain 360 ° weights during implementation of the present invention.In 8 yuan of circle battle arrays shown in Figure 3, the central angle alpha of two adjacent array elements is 45 °, therefore the scan period is 45 °, therefore only need calculate the weights that main lobe points to all assigned directions in 0 ° to 44 °, when then needing to point to other directions, only get final product with the weighted value of rotating certain number of times.When namely needing orientation angle β, calculate first γ wherein 1For
Figure BDA00002368068300052
The merchant, γ 2For
Figure BDA00002368068300053
Remainder, and γ 2<α.Then in 0 to α, find out and point to γ 2The time corresponding weights, again with the counterclockwise γ by turns successively of the weights on each array element 1Inferior, main lobe can orientation angle β.So both can by by turns weighting, realize 360 ° of omnidirectional's scannings.When pointing to 60 ° such as needs, calculate first 60 divided by 45, the merchant is 1, and remainder is 15.Weights corresponding to each array element when then finding out main lobe and pointing to 15 °, as shown in Figure 3, antenna D1, D2, D3, D4, D5, D6, the corresponding weighted value of D7, D8 are respectively W1, W2, W3, W4, W5, W6, W7, W8.All weights are rotated 1 time counterclockwise successively again, namely W1 is assigned to D2 to D8 successively to W7, W8 is assigned to D1, as shown in Figure 4, then main lobe points to 60 °, and the characteristics of directional diagram are constant.In high-frequency ground wave radar, main lobe is every 15 ° of run-downs, then only needs to calculate the weighted value that main lobe points to 0 °, 15 ° and 30 °, only need rotate weighted value by above-mentioned rule when needing other to point to, and is simple and convenient.
Maximization directivity factor method is prior art, for ease of implementing reference, provides related description as follows:
By the theory of Array Signal Processing as can be known, the directional diagram of N element array is:
Figure BDA00002368068300054
Wherein, w=[w 1, w 2..., w N] TBe the weighing vector of array, Be the steering vector of array, subscript H and T represent respectively conjugate transpose and transposition, θ and
Figure BDA00002368068300056
Represent respectively the angle of pitch and azimuth.
For isotropism array element, have:
Figure BDA00002368068300057
Wherein, f 0Be the incident wave frequency,
Figure BDA00002368068300058
Arrive the time delay of each array element for the wavefront of incident plane wave, r is circle battle array radius, and c is electromagnetic wave propagation speed in the air, and N is element number of array.Exp (j2 π f 0τ k) (k=1,2 ..., N) expression is by corresponding phase term of the time of this delay, and j is imaginary unit.
The directivity factor of array is defined as:
Figure BDA00002368068300062
Be the greatest irradiation direction of array, i.e. main lobe direction.
Directivity factor is as follows with the quotient representation of broad sense Rayleigh:
D = w H Bw w H Aw - - - ( 4 )
Wherein,
Figure BDA00002368068300064
Figure BDA00002368068300065
Figure BDA00002368068300066
That main lobe points to.
The method of available generalized eigenvalue decomposition obtains the corresponding weighing vector of maximum D, but is not easy to add constraints.The method travel direction of using numerical value is comprehensive, namely is converted into following problem:
min w w H Aw
(5)
Figure BDA00002368068300069
Wherein Re represents to get real part, subscript θ and
Figure BDA000023680683000611
Expression is respectively to its differentiate.Then this problem equivalence is:
min w w H Aw (6)
st.C Hw=f
Wherein, C is the row non-singular matrix of N * M, and f is the column vector of M * 1, and M is the number of constraints.
Constrained optimization problem has a variety of methods, and such as quadratic programming, method of Lagrange multipliers etc., it is as follows that embodiment draws weight coefficient with method of Lagrange multipliers:
w D=A -1C(C HA -1C) -1f (7)
w DBe the corresponding weighing vector of array maximum directivity coefficient.
For to anti-interference, need to suppress sidelobe level in practice.Make the arithmetic number function
Figure BDA000023680683000613
Be the amplitude function of expectation secondary lobe in the secondary lobe to be optimized zone, the weighing vector of establishing the front array of each iteration is w c, its corresponding directional diagram is
Figure BDA000023680683000614
Then the weighing vector that obtains of each iteration is the solution of following problems.
min w w H Aw
Figure BDA00002368068300072
Figure BDA00002368068300073
Figure BDA00002368068300074
Figure BDA00002368068300075
Last formula is illustrated in m the Linear Constraints that newly increases on the basis of M constraints in the following formula.f iBe the complex representation of Sidelobe Suppression value, use in this median as computing.This m Linear Constraints is used for controlling current directional diagram
Figure BDA00002368068300076
In the amplitude at m side lobe peak point place to set-point
Figure BDA00002368068300077
F then iValue be:
Figure BDA00002368068300078
Wherein,
Figure BDA00002368068300079
Currency for this m some place.
This problem also can equivalence be
min w w H Aw (10)
st.C 1 Hw=g
C wherein 1Be the complex matrix of N * (M+m), g is the column vector of M+m dimension.The side lobe peak point of current directional diagram in searching target zone during each iteration
Figure BDA000023680683000711
Be total to m sIndividual.Because the number of Linear Constraints can not surpass N, so can add m at most in the iteration Max=N-M secondary lobe constraints.If m s>m Max, then choose wherein m sIndividual large secondary lobe is as constraint object, so the number m=min (m of secondary lobe constraint s, m Max).
Still use method of Lagrange multipliers, draw weighing vector and be:
w D=A -1C 1(C 1 HA -1C 1) -1g (11)
By
Figure BDA000023680683000712
Can form main lobe points to
Figure BDA000023680683000713
Wave beam.
The calculation process explanation of the maximization directivity factor method under the secondary lobe constraint of embodiment is provided for the sake of ease of implementation.As shown in Figure 5, in an embodiment step is as follows:
Step 1, input parameter comprises that array number N(is the monopole antenna number), a circle battle array radius r (being the radius of circular little spacing monopole antenna battle array), operating frequency f 0, main lobe direction
Figure BDA000023680683000714
sWith
Figure BDA000023680683000715
Represent respectively the angle of pitch and azimuth), secondary lobe constraints
Figure BDA000023680683000716
(representing the amplitude that i secondary lobe will retrain);
Step 2 is obtained one group of weights with maximization directivity factor method, and (1) to (7) formula has provided the derivation of maximization directivity factor method, will can obtain one group of weight w in other parameter substitutions (7) except secondary lobe constraints s in the step 1 D, namely according to array number N, circle battle array radius r, operating frequency f 0, main lobe direction
Figure BDA00002368068300081
Obtain one group of weight w with maximization directivity factor method D, the formula of asking for employing is w D=A -1C (C HA -1C) -1F, wherein
Figure BDA00002368068300082
C is the row non-singular matrix of N * M;
Step 3 is with the weight w of obtaining DCalculate directional diagram, the weight w of namely obtaining according to step 2 in the substitution (1) DMake up the weighing vector w=[w of array 1, w 2..., w N] T, calculate F is the column vector of M * 1, and M is the number of secondary lobe constraints;
Step 4, whether the sidelobe level of directional diagram all is less than or equal to the sidelobe level inhibiting value in the directional diagram characteristics in the determining step 3
Figure BDA00002368068300084
If it is should organize weights output, process ends, otherwise enter step 5;
Step 5, find out backlog demand (namely greater than the sidelobe level inhibiting value in the directional diagram characteristics) secondary lobe the position and add secondary lobe constraints, soon do not satisfy in the location point substitution (9) of secondary lobe requirement, calculate median f iAfter, in the again substitution (8), and finally be reduced to (10) formula;
Step 6 is obtained weights with maximization directivity factor method under current secondary lobe constraints, namely use (11) formula w D=A -1C 1(C 1 HA -1C 1) -1G obtains the solution w of (10) formula DReturn step 3 and calculate, until sidelobe level all meets the demands.
For the sake of ease of implementation, provide the method for measurement of S parameter for reference: as shown in Figure 6, when adopting the two-port network analyzer to measure S parameter (being the microwave scattering parameter) between two antennas, if two antennas are designated as respectively array element 1 and array element 2, in the environment of circle battle array work, instrument is calibrated coupling first, then port one is connected the feedback point that connects respectively array element 1 and array element 2 with port, network analyzer is set to the single port mode of operation, can measure respectively S11 and S22 parameter under each Frequency point.Network analyzer is set to the dual-port transmission mode again, can measure respectively S12 and S21 parameter under each Frequency point.Owing to can not make the characteristic of two monopoles just the same in the manufacturing process, therefore the numerical value of S12 and S21 is understood some difference, be unequal.According to microwave network theory, the conversion relation of S parameter and mutual resistance matrix is as follows:
Z=Z 0(S+I)(I-S) -1 (12)
Wherein, Z 0Being characteristic impedance, is 50 Ω at this.
Z = Z 11 Z 12 Z 21 Z 22 Be mutual resistance matrix, I = 1 0 0 1 Be the second order unit matrix; S = S 11 S 12 S 21 S 22 Be the S parameter matrix.
In mutual resistance matrix, Z 11And Z 22The self-impedance that represents respectively array element 1 and array element 2, Z 12The mutual impedance of 2 pairs of array elements 1 of expression array element, Z 21The mutual impedance of 1 pair of array element 2 of expression array element.In the S parameter matrix, S 11During expression port 2 coupling, the reflection coefficient of port one; S 22During expression port one coupling, the reflection coefficient of port 2; S 12During expression port one coupling, port 2 is to the reverse transfer coefficient of port one; S 21During expression port 2 coupling, port one is to the forward transmission coefficient of port 2.
Measure the S parameter of all per two antennas according to above method, and be converted into mutual resistance matrix, finally form the mutual resistance matrix of array:
Z = Z 11 Z 12 . . . Z 1 N Z 21 Z 22 . . . Z 2 N . . . . . . . . . Z N 1 Z N 2 . . . Z NN - - - ( 13 )
Then the antenna array that proposes according to Gupta receives equivalent network model, draws normalized mutual resistance matrix and is:
Z u = 1 + Z 11 Z L Z 12 Z L . . . Z 1 N Z L Z 21 Z L 1 + Z 22 Z L . . . Z 2 N Z L . . . . . . . . . Z N 1 Z L Z N 2 Z L . . . 1 + Z NN Z L - - - ( 14 )
Wherein, Z LBe load impedance, be the input impedance of receiver at this, be i.e. 50 Ω.Z uBe final correction data.Weighing vector after then proofreading and correct is w L=w D HZ uUse w LEach antenna is weighted, can obtains required directional diagram.
Fig. 7 and shown in Figure 8 be that said method is used in 8 yuan of radius 2.5m circle battle arrays, be 10MHz, main lobe when pointing to 90 ° in operating frequency, secondary lobe is in the situation of 10dB without constraint and constraint level, the directional diagram that obtains respectively.Among Fig. 7, main lobe width is 32.49 °, and directivity factor is 12.21dB, and the highest sidelobe level is-8.2dB; Main lobe width is 31.34 ° among Fig. 8, and directivity factor is 11.71dB, and the highest sidelobe level is-10dB.Can find out, the superdirective antenna battle array is because antenna distance is very little, and mutual coupling is very large on the impact of directional diagram, must be proofreaied and correct, and the invention provides the effective scheme of mutual coupling calibration.
Specific embodiment described herein only is to the explanation for example of the present invention's spirit.Those skilled in the art can make various modifications or replenish or adopt similar mode to substitute described specific embodiment, but can't depart from spirit of the present invention or surmount the defined scope of appended claims.

Claims (5)

1. a shortwave miniature circular receiving antenna array is characterized in that: comprise the little spacing monopole antenna of circle battle array, matching network, radar radio frequency reception channel, weighting generator, mutual coupling calibration device and a weighter;
Circular little spacing monopole antenna battle array comprises a plurality of monopole antennas, and each monopole antenna is evenly arranged on the circumference, and the radius of circular little spacing monopole antenna battle array is less than or equal to 1/4th of radar operation wavelength;
The output of each monopole antenna connects the input of corresponding matching network, and the output of each matching network is through corresponding radar radio frequency reception channel difference weighted input device;
The weights that the weighting generator produces output to the mutual coupling calibration device; The mutual coupling calibration device forms correction data according to the mutual coupling situation between each monopole antenna, proofread and correct weights and weighted input device that the weighting generator produces, weighter is weighted each monopole antenna according to the input by weights and corresponding radar radio frequency reception channel after the correction of mutual coupling calibration device input, form specific directional diagram, support 360 ° of omnidirectional's scannings.
2. shortwave miniature circular receiving antenna array according to claim 1, it is characterized in that: matching network is five rank broadband matching networks.
3. shortwave miniature circular receiving antenna array according to claim 1 is characterized in that: the radius of monopole antenna number and circular little spacing monopole antenna battle array is definite according to the characteristics of the default directional diagram that will form; The characteristics of described directional diagram comprise main lobe direction, main lobe width and sidelobe level inhibiting value.
4. shortwave miniature circular receiving antenna array according to claim 3, it is characterized in that: the weighting generator is according to the characteristics of the radius of monopole antenna number, circular little spacing monopole antenna battle array and the directional diagram that will form, utilize maximization directivity factor method under the secondary lobe constraint to calculate the weighted value of each antenna element, computational process comprises following substep
Step 1, input parameter comprises radius r, the operating frequency f of monopole antenna number N, circular little spacing monopole antenna battle array 0, main lobe direction
Figure FDA00002368068200011
With secondary lobe constraints
Figure FDA00002368068200012
Wherein, θ sWith
Figure FDA00002368068200013
Represent respectively the angle of pitch and azimuth, secondary lobe constraints
Figure FDA00002368068200014
Represent the amplitude that i secondary lobe will retrain;
Step 2 is according to array number N, circle battle array radius r, operating frequency f 0And main lobe direction
Figure FDA00002368068200015
Obtain one group of weight w with maximization directivity factor method D
Step 3 is with weight w DWeighing vector calculated direction figure as array;
Step 4, whether the sidelobe level of directional diagram all is less than or equal to the sidelobe level inhibiting value in the directional diagram characteristics in the determining step 3
Figure FDA00002368068200016
Then to organize weights output, process ends, otherwise enter step 5;
Step 5, find out backlog demand secondary lobe the position and add secondary lobe constraints;
Step 6 is obtained weights with maximization directivity factor method under current secondary lobe constraints, obtain new w D, return step 3 and calculate, until sidelobe level all meets the demands, export current weight w D
5. according to claim 1 and 2 or 3 or 4 described shortwave miniature circular receiving antenna arrays, it is characterized in that: the mutual coupling calibration device forms correction data according to the mutual coupling situation between each monopole antenna, implementation is proofreaied and correct the weights that the weighting generator produces as correction data with mutual resistance matrix for generate the mutual resistance matrix of array with the S parameter; Described S parameter be antenna with adopt before matching network is connected the two-port network analyzer directly in actual working environment, to measure.
CN201210442671.4A 2012-11-07 2012-11-07 Short-wave small circular receiving antenna array Expired - Fee Related CN102904015B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201210442671.4A CN102904015B (en) 2012-11-07 2012-11-07 Short-wave small circular receiving antenna array

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201210442671.4A CN102904015B (en) 2012-11-07 2012-11-07 Short-wave small circular receiving antenna array

Publications (2)

Publication Number Publication Date
CN102904015A true CN102904015A (en) 2013-01-30
CN102904015B CN102904015B (en) 2014-09-17

Family

ID=47576148

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201210442671.4A Expired - Fee Related CN102904015B (en) 2012-11-07 2012-11-07 Short-wave small circular receiving antenna array

Country Status (1)

Country Link
CN (1) CN102904015B (en)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107037283A (en) * 2017-03-06 2017-08-11 中国人民解放军军械工程学院 Radiosensitivity method of testing and device under the conditions of reverberation chamber
CN113193359A (en) * 2021-04-22 2021-07-30 中国人民解放军海军工程大学 Short wave phased array antenna mutual coupling suppression method based on passive parasitic array elements
US11183774B2 (en) 2019-05-31 2021-11-23 The Mitre Corporation High frequency system using a circular array

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1418028A (en) * 2001-11-06 2003-05-14 深圳市中兴通讯股份有限公司上海第二研究所 Correction method and device of intelligent antenna subsystem
CN1417965A (en) * 2001-11-06 2003-05-14 深圳市中兴通讯股份有限公司上海第二研究所 Fully adaptive weight value updating method and device for array antenna
CN202330704U (en) * 2011-11-25 2012-07-11 武汉大学 Small-sized broadband reception array for multi-frequency high-frequency ground wave radar

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1418028A (en) * 2001-11-06 2003-05-14 深圳市中兴通讯股份有限公司上海第二研究所 Correction method and device of intelligent antenna subsystem
CN1417965A (en) * 2001-11-06 2003-05-14 深圳市中兴通讯股份有限公司上海第二研究所 Fully adaptive weight value updating method and device for array antenna
CN202330704U (en) * 2011-11-25 2012-07-11 武汉大学 Small-sized broadband reception array for multi-frequency high-frequency ground wave radar

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107037283A (en) * 2017-03-06 2017-08-11 中国人民解放军军械工程学院 Radiosensitivity method of testing and device under the conditions of reverberation chamber
CN107037283B (en) * 2017-03-06 2019-08-20 中国人民解放军军械工程学院 Radiosensitivity test method and device under the conditions of reverberation chamber
US11183774B2 (en) 2019-05-31 2021-11-23 The Mitre Corporation High frequency system using a circular array
CN113193359A (en) * 2021-04-22 2021-07-30 中国人民解放军海军工程大学 Short wave phased array antenna mutual coupling suppression method based on passive parasitic array elements
CN113193359B (en) * 2021-04-22 2023-03-31 中国人民解放军海军工程大学 Short wave phased array antenna mutual coupling suppression method based on passive parasitic array elements

Also Published As

Publication number Publication date
CN102904015B (en) 2014-09-17

Similar Documents

Publication Publication Date Title
CN104615854A (en) Beam broadening and sidelobe suppression method based on sparse constraint
CN108462521A (en) The anti-interference realization method of adaptive array antenna
CN105259550A (en) Compressed sensing-based multi-input and multi-output radar two-dimensional angle estimation method
CN103178881B (en) Main lobe interference suppression method and device
CN105929384B (en) The joint low sidelobe Multibeam synthesis method of ring array reception is penetrated based on cross paroxysm
CN101931124A (en) Embattling method of logarithmic spiral array antennas
CN104600428B (en) A kind of system for improving continuous wave radar isolation between transmitting and receiving antenna
CN103000996B (en) Uniform circular array direction-finder antenna receiving mutual impedance test and mutual coupling compensation system
WO2023279737A1 (en) Multi-constraint millimeter-wave vehicle-mounted mimo radar antenna array design method, and apparatus
CN109490820A (en) A kind of arrival direction estimation method based on parallel nested battle array
CN104515909A (en) Large antenna directional pattern measuring method based on correlation method
CN107276658A (en) The Beamforming Method reconstructed under coloured noise based on covariance matrix
CN102721965B (en) Omni-directional side lobe suppression method for synthetic aperture radar imaging processing
CN102904015B (en) Short-wave small circular receiving antenna array
CN110501675A (en) One kind being based on MIMO radar low sidelobe transmitting pattern design method
CN202871978U (en) Short wave miniature circular receiving antenna array
CN104020459A (en) Waveform optimization method for improving MIMO-STAP detection performance
CN103852749A (en) Robust waveform optimization method for improving MIMO-STAP detection performance
CN101483280B (en) Weight solving method for stable wave beam synthesizer
CN104102775A (en) Side-lobe suppression based beam optimization method for electrical tilt parasitic antennas
CN105024166A (en) Planar array antenna pattern synthesis method based on sub arrays
CN113540791B (en) Method for optimizing aperture-level transmit-receive simultaneous array
CN104362437A (en) S-frequency-band monopulse self-tracking antenna system
CN111812607B (en) Meter wave MIMO radar low elevation angle estimation method based on beam space
WO2020114303A1 (en) Radio-frequency energy collection system and method using blind adaptive beamforming algorithm

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20140917

Termination date: 20191107

CF01 Termination of patent right due to non-payment of annual fee