CN111123250A - Pulse Doppler radar based on pattern search algorithm and beam forming method - Google Patents

Pulse Doppler radar based on pattern search algorithm and beam forming method Download PDF

Info

Publication number
CN111123250A
CN111123250A CN201911389423.6A CN201911389423A CN111123250A CN 111123250 A CN111123250 A CN 111123250A CN 201911389423 A CN201911389423 A CN 201911389423A CN 111123250 A CN111123250 A CN 111123250A
Authority
CN
China
Prior art keywords
search
objective function
different
sum beam
coefficient
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
CN201911389423.6A
Other languages
Chinese (zh)
Other versions
CN111123250B (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.)
Chengdu Huirong Guoke Microsystem Technology Co ltd
Original Assignee
Chengdu Huirong Guoke Microsystem 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 Chengdu Huirong Guoke Microsystem Technology Co ltd filed Critical Chengdu Huirong Guoke Microsystem Technology Co ltd
Priority to CN201911389423.6A priority Critical patent/CN111123250B/en
Publication of CN111123250A publication Critical patent/CN111123250A/en
Application granted granted Critical
Publication of CN111123250B publication Critical patent/CN111123250B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
    • G01S13/02—Systems using reflection of radio waves, e.g. primary radar systems; Analogous systems
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/02—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
    • G01S7/28—Details of pulse systems
    • G01S7/2813—Means providing a modification of the radiation pattern for cancelling noise, clutter or interfering signals, e.g. side lobe suppression, side lobe blanking, null-steering arrays
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/02—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
    • G01S7/28—Details of pulse systems
    • G01S7/285—Receivers
    • G01S7/292—Extracting wanted echo-signals
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
    • G01S13/02—Systems using reflection of radio waves, e.g. primary radar systems; Analogous systems
    • G01S2013/0236—Special technical features
    • G01S2013/0245—Radar with phased array antenna

Landscapes

  • Engineering & Computer Science (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Radar Systems Or Details Thereof (AREA)

Abstract

The invention provides a pulse Doppler radar and a beam forming method based on a pattern search algorithm, which comprises the following steps: s1: measuring radar receiving gain of a front array surface of the radar phased array according to different azimuth angles and pitching angles to obtain different array element receiving antenna directional diagrams; s2: initializing numbers and beam generation coefficients, designing different search step lengths, determining base vectors in different search directions, and initializing convergence conditions; s3: determining a sum beam generation coefficient objective function according to the sum beam distribution characteristics; s4: solving a target function value corresponding to the sum beam generation coefficient, and judging whether the target function reaches a convergence condition; s5: and obtaining different sum beam generation coefficients and corresponding objective function values by using different step lengths and basis vectors, and selecting the sum beam generation coefficient which enables the objective function value to be minimum as a final result for outputting. The invention synthesizes ideal and channel antenna directional patterns by a mode search optimization method.

Description

Pulse Doppler radar based on pattern search algorithm and beam forming method
Technical Field
The invention belongs to the technical field of radars, and particularly relates to a pulse Doppler radar based on a pattern search algorithm and a beam forming method.
Background
Radar, i.e. radio detection and ranging, as the name suggests, one of the basic functions of radar is target detection. When the radar completes the digital beam synthesis function, it is generally assumed that the front end of the radar phased array has completed amplitude-phase consistency calibration of each array element, that is, there is no difference between amplitude and phase between each array element, but in practical engineering application, because the amplitude and phase of the echo of each array element are different due to inconsistency of the manufacturing elements or manufacturing processes of each array element, there will be a certain difference between amplitude and phase, and under the condition of the amplitude-phase difference, if a traditional "sum beam" forming method is used, the problem that the synthesized beam may have asymmetry of the main lobe and attenuation of the main lobe gain will be caused. The attenuation of the main lobe gain can cause the attenuation of radar detection power, and the asymmetry of the main lobe can cause the reduction of subsequent radar angle measurement precision, so that how to synthesize ideal 'sum beams' through the echo data of each array element with amplitude-phase difference has important engineering application value.
Disclosure of Invention
The invention provides a pulse Doppler radar and beam forming method based on a pattern search algorithm based on the existing problems, and aims to solve the problem that when amplitude differences exist among array elements at the front end, an effective sum beam generation coefficient is obtained through a pattern search optimization algorithm.
S1: according to a front-end array surface of the radar phased array, measuring radar receiving gains of the array surface in different directions according to different azimuth angles and pitching angles, and obtaining different array element receiving antenna directional diagrams;
s2: initializing numbers and beam generation coefficients, designing different search step lengths, determining base vectors in different search directions, and initializing convergence conditions;
s3: determining a sum beam generation coefficient objective function according to the sum beam distribution characteristics;
s4: solving a target function value corresponding to the sum beam generation coefficient, judging whether the target function reaches a convergence condition, if not, determining the direction according to the basis vector, taking the step length as stepping to obtain the next sum beam generation coefficient, solving the corresponding target function, if the convergence condition is reached, ending the search, and generating the current sum beam into the coefficient;
s5: and (4) performing the operation of the step (S4) for multiple times according to the different step lengths and the basis vectors determined in the step (S2) to obtain different sum beam generation coefficients and corresponding objective function values, and selecting the sum beam generation coefficient which enables the objective function value to be minimum as a final result to be output.
In another embodiment of the above pulse doppler radar and beam forming method based on the pattern search algorithm of the present invention, the obtaining of the directional patterns of the different array element receiving antennas according to the radar receiving gains of the radar phased array front end array surfaces in different directions and according to different azimuth angles and elevation angles includes:
according to the azimuth angle range of the front array surface of the radar phased array
Figure BDA0002344522170000021
In angular steps
Figure BDA0002344522170000022
Obtaining a set of azimuthal measurements for a step
Figure BDA0002344522170000023
According to the pitch angle range [ theta ] of the front-end array surface of the radar phased arraymax，θmin]Taking the angle step delta theta as a stepping to obtain a pitch angle measurement set theta;
to obtainSet of angles
Figure BDA0002344522170000024
Echo power of each receiving array element at front end of phased array for measuring each angle coordinate in angle set
Figure BDA0002344522170000025
And m x n is the area array dimension of the front array surface of the phased array.
In another embodiment of the above pulse doppler radar and beam forming method based on the pattern search algorithm of the present invention, the initialization numbers and the beam generating coefficients are designed to have different search step lengths and to determine basis vectors in different search directions, and the initialization convergence condition includes:
initializing numbers and beamforming coefficients
Figure BDA0002344522170000026
Wherein 0 represents an initial time;
initial search step λ ═ λ1,λ2,...,λk]Designing a plurality of different search step lengths, and selecting a plurality of search step lengths to compare optimization results;
determining search basis vectors V ═ V in different directions1,v2,...,vk]Determining search coordinates at different moments in cooperation with the search step length;
initializing a convergence condition, wherein the search condition is used for judging whether the convergence result reaches a set index, and the convergence condition is as follows: and (4) giving a search threshold delta, judging that the search is finished when the objective function value is smaller than the search threshold, and returning an optimal value result.
In another embodiment of the above pulse doppler radar and beam forming method based on the pattern search algorithm of the present invention, the determining a sum beam generation coefficient objective function from the sum beam distribution characteristics includes:
the echo power of each receiving array element at the front end of the phased array
Figure BDA0002344522170000031
With initial sum beam-forming coefficients
Figure BDA0002344522170000032
The initial sum beam is calculated in combination with the following equation:
Figure BDA0002344522170000033
determining a corresponding objective function from the sum beam main lobe peak power
Figure BDA0002344522170000034
In the formula (f)1(w) determining a corresponding objective function, P, for the peak power of the sum beam mainlobepeakIs the peak power, the position L of the peak powerpeakIs the peak power PpeakThe location of the location;
from the position L of the peak power of the main lobepeakDetermining a corresponding objective function to the width of the main lobe
Figure BDA0002344522170000035
Figure BDA0002344522170000036
Ppeak-f2(w)≤ΔgmwIn the formula,. DELTA.gmwIs a main lobe gain attenuation value, AmwIs the main lobe width;
obtaining a corresponding objective function according to the included angle between the peak value of the first minor lobe and the main lobe and the width of the minor lobe
Figure BDA0002344522170000037
Figure BDA0002344522170000038
Figure BDA0002344522170000039
Ppeak-f3(w)≥Δgsw；
In the formula,. DELTA.gswThe ratio of the main lobe to the auxiliary lobe is,
Figure BDA00023445221700000310
is the included angle between the peak value of the first minor lobe and the main lobe, AswThe minor lobe width.
Obtaining a final objective function f (w) ═ k1f1(w)+k2f2(w)+k3f3(w) in the formula, k1,k2,k3Representing the weight coefficients of the respective objective functions.
In another embodiment of the above pulse doppler radar and beam forming method based on the pattern search algorithm of the present invention, the solving and beam forming coefficients correspond to objective function values, and determine whether the objective function meets a convergence condition, if not, determine a direction according to a basis vector, obtain a next sum beam forming coefficient by taking a step length as a step, solve the corresponding objective function, if the convergence condition is met, end the search, and include the current sum beam forming coefficient:
from the initial weight coefficient
Figure BDA00023445221700000311
And the objective function obtains an initial objective function value
Figure BDA00023445221700000312
Judging objective function value
Figure BDA0002344522170000041
Whether the value is less than a search termination threshold delta;
if less than, output the optimal weight coefficient
Figure BDA0002344522170000042
Otherwise, calculating by the search base vector and the search step length
Figure BDA0002344522170000043
Determining a next generation weight coefficient;
until the corresponding objective function meets the search termination condition, outputting the corresponding weight coefficient w1opt。
In another embodiment of the above pulse doppler radar and beam forming method based on the pattern search algorithm according to the present invention, the performing the operation of step S4 for a plurality of times according to the different step sizes and basis vectors determined in step S2 to obtain different sum beam generating coefficients and corresponding objective function values, and selecting the sum beam generating coefficient that minimizes the objective function value as a final result for outputting includes:
different iteration step lengths are selected to obtain corresponding weight coefficients [ w ]2opt,w3opt,...,wkopt]；
Finding out optimal weight coefficient set woptThe optimal weight coefficient w for minimizing the objective functionoptAs a final output result;
calculation of the resulting optimal weight coefficients
Figure BDA0002344522170000044
A sum beam antenna pattern is obtained.
According to the pulse Doppler radar and the beam forming method based on the mode search algorithm, when phase differences exist among array elements of a radar receiver, the problems that detection probability is reduced and angle measurement errors are large in partial angles due to the fact that gain of a directional diagram of the receiver is asymmetric when the radar receiver is synthesized and channels are synthesized by a traditional method are solved.
Drawings
Fig. 1 is a flowchart of an embodiment of a pulse doppler radar and a beam forming method based on a pattern search algorithm according to the present invention.
Fig. 2 is a measured data sum beam azimuth elevation antenna pattern obtained based on the method of the present invention.
Detailed Description
In order to make the objects, technical solutions and advantages of the embodiments of the present invention clearer, the technical solutions in the embodiments of the present invention will be clearly and completely described below with reference to the drawings in the embodiments of the present invention, and it is obvious that the described embodiments are only a part of the embodiments of the present invention, and not all of the embodiments. All other embodiments, which can be derived by a person skilled in the art from the embodiments given herein without making any creative effort, shall fall within the protection scope of the present invention.
The following describes a pulse doppler radar and a beam forming method based on a pattern search algorithm in more detail with reference to the accompanying drawings and embodiments.
As shown in fig. 1, the pulse doppler radar and beamforming method based on the pattern search algorithm includes:
s1: according to a front-end array surface of the radar phased array, measuring radar receiving gains of the array surface in different directions according to different azimuth angles and pitching angles, and obtaining different array element receiving antenna directional diagrams;
s2: initializing numbers and beam generation coefficients, designing different search step lengths, determining base vectors in different search directions, and initializing convergence conditions;
s3: determining a sum beam generation coefficient objective function according to the sum beam distribution characteristics;
s4: solving a target function value corresponding to the sum beam generation coefficient, judging whether the target function reaches a convergence condition, if not, determining the direction according to the basis vector, taking the step length as stepping to obtain the next sum beam generation coefficient, solving the corresponding target function, if the convergence condition is reached, ending the search, and generating the current sum beam into the coefficient;
s5: and (4) performing the operation of the step (S4) for multiple times according to the different step lengths and the basis vectors determined in the step (S2) to obtain different sum beam generation coefficients and corresponding objective function values, and selecting the sum beam generation coefficient which enables the objective function value to be minimum as a final result to be output.
In another embodiment of the above pulse doppler radar and beam forming method based on the pattern search algorithm of the present invention, the obtaining of the directional patterns of the different array element receiving antennas according to the radar receiving gains of the radar phased array front end array surfaces in different directions and according to different azimuth angles and elevation angles includes:
according to the azimuth angle range of the front array surface of the radar phased array
Figure BDA0002344522170000051
In angular steps
Figure BDA0002344522170000052
Obtaining a set of azimuthal measurements for a step
Figure BDA0002344522170000053
According to the pitch angle range [ theta ] of the front-end array surface of the radar phased arraymax，θmin]Taking the angle step delta theta as a stepping to obtain a pitch angle measurement set theta;
obtain a set of angles
Figure BDA0002344522170000054
Echo power of each receiving array element at front end of phased array for measuring each angle coordinate in angle set
Figure BDA00023445221700000610
And m x n is the area array dimension of the front array surface of the phased array.
In another embodiment of the above pulse doppler radar and beam forming method based on the pattern search algorithm of the present invention, the initialization numbers and the beam generating coefficients are designed to have different search step lengths and to determine basis vectors in different search directions, and the initialization convergence condition includes:
initializing numbers and beamforming coefficients
Figure BDA0002344522170000061
Wherein 0 represents an initial time;
initial search step λ ═ λ1,λ2,...,λk]Designing a plurality of different search step lengths, and selecting a plurality of search step lengths to compare optimization results;
determining search basis vectors V ═ V in different directions1,v2,...,vk]Determining search coordinates at different moments in cooperation with the search step length;
initializing a convergence condition, wherein the search condition is used for judging whether the convergence result reaches a set index, and the convergence condition is as follows: and (4) giving a search threshold delta, judging that the search is finished when the objective function value is smaller than the search threshold, and returning an optimal value result.
In another embodiment of the above pulse doppler radar and beam forming method based on the pattern search algorithm of the present invention, the determining a sum beam generation coefficient objective function from the sum beam distribution characteristics includes:
the echo power of each receiving array element at the front end of the phased array
Figure BDA0002344522170000062
With initial sum beam-forming coefficients
Figure BDA0002344522170000063
The initial sum beam is calculated in combination with the following equation:
Figure BDA0002344522170000064
determining a corresponding objective function from the sum beam main lobe peak power
Figure BDA0002344522170000065
In the formula (f)1(w) determining a corresponding objective function, P, for the peak power of the sum beam mainlobepeakIs the peak power, the position L of the peak powerpeakIs the peak power PpeakThe location of the location;
from the position L of the peak power of the main lobepeakDetermining a corresponding objective function to the width of the main lobe
Figure BDA0002344522170000066
Figure BDA0002344522170000067
Ppeak-f2(w)≤ΔgmwIn the formula,. DELTA.gmwIs a main lobe gain attenuation value, AmwIs the main lobe width;
obtaining a corresponding objective function according to the included angle between the peak value of the first minor lobe and the main lobe and the width of the minor lobe
Figure BDA0002344522170000068
Figure BDA0002344522170000069
Figure BDA0002344522170000071
Ppeak-f3(w)≥Δgsw；
In the formula,. DELTA.gswThe ratio of the main lobe to the auxiliary lobe is,
Figure BDA0002344522170000072
is the included angle between the peak value of the first minor lobe and the main lobe, AswThe minor lobe width.
Obtaining a final objective function f (w) ═ k1f1(w)+k2f2(w)+k3f3(w) in the formula, k1,k2,k3Representing the weight coefficients of the respective objective functions.
In another embodiment of the above pulse doppler radar and beam forming method based on the pattern search algorithm of the present invention, the solving and beam forming coefficients correspond to objective function values, and determine whether the objective function meets a convergence condition, if not, determine a direction according to a basis vector, obtain a next sum beam forming coefficient by taking a step length as a step, solve the corresponding objective function, if the convergence condition is met, end the search, and include the current sum beam forming coefficient:
from the initial weight coefficient
Figure BDA0002344522170000073
And the objective function obtains an initial objective function value
Figure BDA0002344522170000074
Judging objective function value
Figure BDA0002344522170000075
Whether the value is less than a search termination threshold delta;
if less than, output the optimal weight coefficient
Figure BDA0002344522170000076
Otherwise, calculating by the search base vector and the search step length
Figure BDA0002344522170000077
Determining a next generation weight coefficient;
until the corresponding objective function meets the search termination condition, outputting the corresponding weight coefficient w1opt。
In another embodiment of the above pulse doppler radar and beam forming method based on the pattern search algorithm according to the present invention, the performing the operation of step S4 for a plurality of times according to the different step sizes and basis vectors determined in step S2 to obtain different sum beam generating coefficients and corresponding objective function values, and selecting the sum beam generating coefficient that minimizes the objective function value as a final result for outputting includes:
different iteration step lengths are selected to obtain corresponding weight coefficients [ w ]2opt,w3opt,...,wkopt]；
Finding out optimal weight coefficient set woptThe optimal weight coefficient w for minimizing the objective functionoptAs a final output result;
calculation of the resulting optimal weight coefficients
Figure BDA0002344522170000078
A sum beam antenna pattern is obtained.
FIG. 2 is a diagram of a sum beam azimuth elevation antenna pattern obtained based on measured data according to the method of the present invention, as shown in FIG. 2, and in FIG. 2, a main lobe gain attenuation value Δ g is setmw2, main-to-side lobe ratio Δ g sw25; the result analysis shows that the method of the invention can obtain a more ideal antenna directional pattern.
It will be evident to those skilled in the art that the embodiments of the present invention are not limited to the details of the foregoing illustrative embodiments, and that the embodiments of the present invention are capable of being embodied in other specific forms without departing from the spirit or essential attributes thereof. The present embodiments are therefore to be considered in all respects as illustrative and not restrictive, the scope of the embodiments being indicated by the appended claims rather than by the foregoing description, and all changes which come within the meaning and range of equivalency of the claims are therefore intended to be embraced therein. Any reference sign in a claim should not be construed as limiting the claim concerned. Furthermore, it is obvious that the word "comprising" does not exclude other elements or steps, and the singular does not exclude the plural. Several units, modules or means recited in the system, apparatus or terminal claims may also be implemented by one and the same unit, module or means in software or hardware. The terms first, second, etc. are used to denote names, but not any particular order.
Finally, it should be noted that the above embodiments are only used for illustrating the technical solutions of the embodiments of the present invention and not for limiting, and although the embodiments of the present invention are described in detail with reference to the above preferred embodiments, it should be understood by those skilled in the art that modifications or equivalent substitutions can be made on the technical solutions of the embodiments of the present invention without departing from the spirit and scope of the technical solutions of the embodiments of the present invention.

Claims (6)

1. A pulse Doppler radar and a beam forming method based on a pattern search algorithm are characterized by comprising the following steps:
s1: according to a front-end array surface of the radar phased array, measuring radar receiving gains of the array surface in different directions according to different azimuth angles and pitching angles, and obtaining different array element receiving antenna directional diagrams;
s2: initializing numbers and beam generation coefficients, designing different search step lengths, determining base vectors in different search directions, and initializing convergence conditions;
s3: determining a sum beam generation coefficient objective function according to the sum beam distribution characteristics;
s4: solving a target function value corresponding to the sum beam generation coefficient, judging whether the target function reaches a convergence condition, if not, determining the direction according to the basis vector, taking the step length as stepping to obtain the next sum beam generation coefficient, solving the corresponding target function, if the convergence condition is reached, ending the search, and generating the current sum beam into the coefficient;
s5: and (4) performing the operation of the step (S4) for multiple times according to the different step lengths and the basis vectors determined in the step (S2) to obtain different sum beam generation coefficients and corresponding objective function values, and selecting the sum beam generation coefficient which enables the objective function value to be minimum as a final result to be output.
2. The pulse doppler radar and beam forming method based on the pattern search algorithm according to claim 1, wherein the obtaining of the antenna pattern of different array elements according to the radar receiving gain of the radar phased array front end array in different directions according to the measurement front of different azimuth angles and elevation angles comprises:
according to the azimuth angle range of the front array surface of the radar phased array
Figure FDA0002344522160000011
In angular steps
Figure FDA0002344522160000012
Obtaining a set of azimuthal measurements for a step
Figure FDA0002344522160000013
According to the pitch angle range [ theta ] of the front-end array surface of the radar phased arraymax，θmin]Taking the angle step delta theta as a stepping to obtain a pitch angle measurement set theta;
obtain a set of angles
Figure FDA0002344522160000014
Echo power of each receiving array element at front end of phased array for measuring each angle coordinate in angle set
Figure FDA0002344522160000015
And m x n is the area array dimension of the front array surface of the phased array.
3. The pulse doppler radar and beamforming method according to claim 2, wherein the initialization numbers and the beamforming coefficients are designed for different search steps and determining basis vectors for different search directions, and the initialization of convergence conditions includes:
initializing numbers and beamforming coefficients
Figure FDA0002344522160000016
Wherein 0 represents an initial time;
initial search step λ ═ λ1,λ2,...,λk]Designing a plurality of different search step lengths, and selecting a plurality of search step lengths to compare optimization results;
determining search basis vectors V ═ V in different directions1,v2,...,vk]Determining search coordinates at different moments in cooperation with the search step length;
initializing a convergence condition, wherein the search condition is used for judging whether the convergence result reaches a set index, and the convergence condition is as follows: and (4) giving a search threshold delta, judging that the search is finished when the objective function value is smaller than the search threshold, and returning an optimal value result.
4. The pulse doppler radar and beamforming method based on a pattern search algorithm according to claim 3, wherein the determining of the sum beam generation coefficient objective function from the sum beam profile feature comprises:
the echo power of each receiving array element at the front end of the phased array
Figure FDA0002344522160000021
With initial sum beam-forming coefficients
Figure FDA0002344522160000022
The initial sum beam is calculated in combination with the following equation:
Figure FDA0002344522160000023
determining a corresponding objective function from the sum beam main lobe peak power
Figure FDA0002344522160000024
In the formula (f)1(w) determining a corresponding objective function, P, for the peak power of the sum beam mainlobepeakIs the peak power, the position L of the peak powerpeakIs the peak power PpeakThe location of the location;
from the position L of the peak power of the main lobepeakDetermining a corresponding objective function to the width of the main lobe
Figure FDA0002344522160000025
Figure FDA0002344522160000026
Ppeak-f2(w)≤ΔgmwIn the formula,. DELTA.gmwIs a main lobe gain attenuation value, AmwIs the main lobe width;
obtaining a corresponding objective function according to the included angle between the peak value of the first minor lobe and the main lobe and the width of the minor lobe
Figure FDA0002344522160000027
Figure FDA0002344522160000028
Figure FDA0002344522160000029
Ppeak-f3(w)≥Δgsw；
In the formula,. DELTA.gswThe ratio of the main lobe to the auxiliary lobe is,
Figure FDA00023445221600000210
is the included angle between the peak value of the first minor lobe and the main lobe, AswThe minor lobe width.
Obtaining a final objective function f (w) ═ k1f1(w)+k2f2(w)+k3f3(w) in the formula, k1,k2,k3Representing the weight coefficients of the respective objective functions.
5. The pulse doppler radar and beam forming method according to claim 4, wherein the solving and beam forming coefficients correspond to objective function values, determine whether the objective function reaches a convergence condition, if not, determine a direction according to the basis vectors, obtain a next sum beam forming coefficient by taking a step length as a step, solve the corresponding objective function, if the convergence condition is reached, end the search, and include the current sum beam forming coefficient:
from the initial weight coefficient
Figure FDA0002344522160000031
And the objective function obtains an initial objective function value
Figure FDA0002344522160000032
Judging objective function value
Figure FDA0002344522160000033
Whether the value is less than a search termination threshold delta;
if less than, output the optimal weight coefficient
Figure FDA0002344522160000034
Otherwise, calculating by the search base vector and the search step length
Figure FDA0002344522160000035
Determining a next generation weight coefficient;
until the corresponding objective function meets the search termination condition, outputting the corresponding weight coefficient w1opt。
6. The pulse doppler radar and beam forming method based on pattern search algorithm of claim 5, wherein the step S4 is performed a plurality of times according to the different step sizes and basis vectors determined in the step S2 to obtain different sum beam generating coefficients and corresponding objective function values, and the selecting the sum beam generating coefficient that minimizes the objective function value as a final result comprises:
different iteration step lengths are selected to obtain corresponding weight coefficients [ w ]2opt,w3opt,...,wkopt]；
Finding out optimal weight coefficient set woptThe optimal weight coefficient w for minimizing the objective functionoptAs a final output result;
calculation of the resulting optimal weight coefficients
Figure FDA0002344522160000036
A sum beam antenna pattern is obtained.
CN201911389423.6A 2019-12-30 2019-12-30 Pulse Doppler radar based on pattern search algorithm and beam forming method Active CN111123250B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201911389423.6A CN111123250B (en) 2019-12-30 2019-12-30 Pulse Doppler radar based on pattern search algorithm and beam forming method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201911389423.6A CN111123250B (en) 2019-12-30 2019-12-30 Pulse Doppler radar based on pattern search algorithm and beam forming method

Publications (2)

Publication Number Publication Date
CN111123250A true CN111123250A (en) 2020-05-08
CN111123250B CN111123250B (en) 2021-09-28

Family

ID=70504519

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201911389423.6A Active CN111123250B (en) 2019-12-30 2019-12-30 Pulse Doppler radar based on pattern search algorithm and beam forming method

Country Status (1)

Country Link
CN (1) CN111123250B (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113221337A (en) * 2021-04-23 2021-08-06 电子科技大学 Near-field sparse array design method based on pattern search

Citations (22)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20030020646A1 (en) * 2001-06-15 2003-01-30 Kai-Bor Yu Adaptive digital sub-array beamforming and deterministic sum and difference beamforming, with jamming cancellation and monopulse ratio preservation
US20030117314A1 (en) * 2001-08-16 2003-06-26 Lockheed Martin Corporation Pri-staggered post-doppler adaptive monopulse processing for detection and location of a moving target in ground clutter
US6781541B1 (en) * 2003-07-30 2004-08-24 Raytheon Company Estimation and correction of phase for focusing search mode SAR images formed by range migration algorithm
CN1578424A (en) * 2003-07-30 2005-02-09 三星电子株式会社 Method and apparatus for receiving digital television signals using space diversity and beam-forming
US20050190427A1 (en) * 2004-02-21 2005-09-01 Frank Steinsiek Method and apparatus for transmitting energy via a laser beam
US7671789B1 (en) * 2008-10-03 2010-03-02 Lockheed Martin Corporation Method and system for target detection and angle estimation based on a radar signal
CN101995566A (en) * 2010-10-15 2011-03-30 西安电子科技大学 System and method for forming digital wave beams of two-dimensional digital array radar
CN102164374A (en) * 2011-05-10 2011-08-24 华为技术有限公司 Method, device and system for processing beam searching
US20120268314A1 (en) * 2011-02-11 2012-10-25 Honda Elesys Co., Ltd. Multibeam radar apparatus for vehicle, multibeam radar method, and multibeam radar program
CN103389493A (en) * 2013-06-25 2013-11-13 西安电子科技大学 Multi-beam single-pulse angle measuring method based on beam selection method
CN103616679A (en) * 2013-11-19 2014-03-05 北京航空航天大学 PD radar ranging and angle measuring method based on difference beam modulation and waveform analysis
CN104144036A (en) * 2014-07-24 2014-11-12 西安交通大学 PU2RC-based MU-MIMO user selection method
CN104391276A (en) * 2014-10-08 2015-03-04 西安电子工程研究所 Transmit-receive split planar array phased-array radar antenna array and beam former
CN105182292A (en) * 2015-08-24 2015-12-23 电子科技大学 Multi-waveform phase coding method based on mode search algorithm
CN106443663A (en) * 2016-09-30 2017-02-22 西安电子科技大学 Method for dimensional-reduction four-channel sum-difference beam angle measurement of phased array radar
CN106772260A (en) * 2017-03-31 2017-05-31 西安电子科技大学 Radar array and difference beam directional diagram optimization method based on convex optimized algorithm
CN107064880A (en) * 2017-04-10 2017-08-18 南京航空航天大学 Distributed many base radar transmit-receive beam Synchronizations and the accuracy method of wave beam control
CN107390039A (en) * 2017-08-22 2017-11-24 湖南卫导信息科技有限公司 A kind of method of quick test Antenna Array Pattern
CN107390197A (en) * 2017-08-04 2017-11-24 西安电子科技大学 The radar self-adaption and difference beam angle-measuring method in feature based space
CN109541557A (en) * 2018-12-30 2019-03-29 成都汇蓉国科微系统技术有限公司 A kind of radar antenna protection channel generation method
CN109946664A (en) * 2019-03-06 2019-06-28 西安电子科技大学 A kind of array radar target seeker Monopulse estimation method under major lobe suppression
CN110488255A (en) * 2019-09-09 2019-11-22 上海无线电设备研究所 A kind of phased-array radar pulse high-resolution angle measuring system and method

Patent Citations (22)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20030020646A1 (en) * 2001-06-15 2003-01-30 Kai-Bor Yu Adaptive digital sub-array beamforming and deterministic sum and difference beamforming, with jamming cancellation and monopulse ratio preservation
US20030117314A1 (en) * 2001-08-16 2003-06-26 Lockheed Martin Corporation Pri-staggered post-doppler adaptive monopulse processing for detection and location of a moving target in ground clutter
US6781541B1 (en) * 2003-07-30 2004-08-24 Raytheon Company Estimation and correction of phase for focusing search mode SAR images formed by range migration algorithm
CN1578424A (en) * 2003-07-30 2005-02-09 三星电子株式会社 Method and apparatus for receiving digital television signals using space diversity and beam-forming
US20050190427A1 (en) * 2004-02-21 2005-09-01 Frank Steinsiek Method and apparatus for transmitting energy via a laser beam
US7671789B1 (en) * 2008-10-03 2010-03-02 Lockheed Martin Corporation Method and system for target detection and angle estimation based on a radar signal
CN101995566A (en) * 2010-10-15 2011-03-30 西安电子科技大学 System and method for forming digital wave beams of two-dimensional digital array radar
US20120268314A1 (en) * 2011-02-11 2012-10-25 Honda Elesys Co., Ltd. Multibeam radar apparatus for vehicle, multibeam radar method, and multibeam radar program
CN102164374A (en) * 2011-05-10 2011-08-24 华为技术有限公司 Method, device and system for processing beam searching
CN103389493A (en) * 2013-06-25 2013-11-13 西安电子科技大学 Multi-beam single-pulse angle measuring method based on beam selection method
CN103616679A (en) * 2013-11-19 2014-03-05 北京航空航天大学 PD radar ranging and angle measuring method based on difference beam modulation and waveform analysis
CN104144036A (en) * 2014-07-24 2014-11-12 西安交通大学 PU2RC-based MU-MIMO user selection method
CN104391276A (en) * 2014-10-08 2015-03-04 西安电子工程研究所 Transmit-receive split planar array phased-array radar antenna array and beam former
CN105182292A (en) * 2015-08-24 2015-12-23 电子科技大学 Multi-waveform phase coding method based on mode search algorithm
CN106443663A (en) * 2016-09-30 2017-02-22 西安电子科技大学 Method for dimensional-reduction four-channel sum-difference beam angle measurement of phased array radar
CN106772260A (en) * 2017-03-31 2017-05-31 西安电子科技大学 Radar array and difference beam directional diagram optimization method based on convex optimized algorithm
CN107064880A (en) * 2017-04-10 2017-08-18 南京航空航天大学 Distributed many base radar transmit-receive beam Synchronizations and the accuracy method of wave beam control
CN107390197A (en) * 2017-08-04 2017-11-24 西安电子科技大学 The radar self-adaption and difference beam angle-measuring method in feature based space
CN107390039A (en) * 2017-08-22 2017-11-24 湖南卫导信息科技有限公司 A kind of method of quick test Antenna Array Pattern
CN109541557A (en) * 2018-12-30 2019-03-29 成都汇蓉国科微系统技术有限公司 A kind of radar antenna protection channel generation method
CN109946664A (en) * 2019-03-06 2019-06-28 西安电子科技大学 A kind of array radar target seeker Monopulse estimation method under major lobe suppression
CN110488255A (en) * 2019-09-09 2019-11-22 上海无线电设备研究所 A kind of phased-array radar pulse high-resolution angle measuring system and method

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
崔卫东: "基于模式搜索算法的仅相位加权相控阵天线波束赋形", 《火控雷达技术》 *
徐海; 崔连虎: "射频仿真目标回波的天线方向图调制方法", 《国外电子测量技术》 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113221337A (en) * 2021-04-23 2021-08-06 电子科技大学 Near-field sparse array design method based on pattern search
CN113221337B (en) * 2021-04-23 2022-09-09 电子科技大学 Near-field sparse array design method based on pattern search

Also Published As

Publication number Publication date
CN111123250B (en) 2021-09-28

Similar Documents

Publication Publication Date Title
CN104237844B (en) Distributed meter wave planar array radar azimuth angle measuring method based on phase compensation
CN107919535B (en) three-dimensional array antenna based on directional double circular arrays and construction method thereof
CN110058204B (en) Satellite-borne antenna beam center calibration method based on directional diagram matching
CN111046591B (en) Joint estimation method for sensor amplitude-phase error and target arrival angle
CN113835063B (en) Unmanned aerial vehicle array amplitude and phase error and signal DOA joint estimation method
CN104375129B (en) A kind of distributive array coherent synthesis radar emission coherent parameter calibrating method
CN112596022B (en) Wave arrival angle estimation method of low-orbit satellite-borne multi-beam regular hexagonal phased array antenna
CN117310707B (en) On-orbit azimuth direction diagram extraction method for azimuth multichannel satellite-borne SAR
CN111123250B (en) Pulse Doppler radar based on pattern search algorithm and beam forming method
CN113671485B (en) ADMM-based two-dimensional DOA estimation method for meter wave area array radar
KR101239165B1 (en) Method and apparatus for estimating target direction
CN106990385B (en) Non-homogeneous face battle array design method based on minimum-redundancy linear arrays
CN109783960A (en) A kind of Wave arrival direction estimating method based on meshing refinement
CN117741597A (en) MIMO radar antenna rapid calibration method based on matching metric
CN116148784B (en) Automatic phase calibration system and method for single-station lightning positioning system
CN109613474B (en) Angle measurement compensation method suitable for short-distance vehicle-mounted radar
CN109541557B (en) Radar antenna protection channel generation method
CN107144813B (en) Method and device for constructing four-array-element three-dimensional array
CN114063005B (en) Maximum posterior direction of arrival estimation method based on fusion center feedback information
CN109061594A (en) A kind of angle-measuring method based on the processing of submatrix space-time adaptive under strong clutter environment
CN114488142A (en) Radar two-dimensional angle imaging method and system based on difference-sum beam
CN111965602B (en) Phased array radar amplitude-phase consistency detection method and system
RU2686113C1 (en) Method of amplitude two-dimensional direction-finding
CN114234904B (en) Angle measuring method, system and medium based on lens multi-beam antenna
AU2019204348B2 (en) Method and system for determining an angle of arrival of a radioelectric signal

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