CN103186697A - Method for AEWA (Airborne Early Warning Aircraft) radar detection target simulation - Google Patents

Method for AEWA (Airborne Early Warning Aircraft) radar detection target simulation Download PDF

Info

Publication number
CN103186697A
CN103186697A CN2013100982335A CN201310098233A CN103186697A CN 103186697 A CN103186697 A CN 103186697A CN 2013100982335 A CN2013100982335 A CN 2013100982335A CN 201310098233 A CN201310098233 A CN 201310098233A CN 103186697 A CN103186697 A CN 103186697A
Authority
CN
China
Prior art keywords
guided missile
rho
target
target guided
prime
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.)
Pending
Application number
CN2013100982335A
Other languages
Chinese (zh)
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.)
Xidian University
Original Assignee
Xidian University
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 Xidian University filed Critical Xidian University
Priority to CN2013100982335A priority Critical patent/CN103186697A/en
Publication of CN103186697A publication Critical patent/CN103186697A/en
Pending legal-status Critical Current

Links

Images

Landscapes

  • Radar Systems Or Details Thereof (AREA)

Abstract

The invention discloses a method for AEWA (Airborne Early Warning Aircraft) radar detection target simulation and mainly solves the problem that a radar cross section of a target guided missile in a half space cannot be calculated under high frequency in the prior art. The method is implemented by the following steps: establishing a model schematic diagram of the target guided missile by using commercial HFSS (High Frequency Structure Simulator) software, and performing triangulation on the model schematic diagram of the target guided missile by using commercial femap software; importing the model schematic diagram of the target guided missile after triangulation in the form of STL (Standard Template Library) text to acquire the profile data information of the target guided missile in the form of the STL text; introducing a green function in the half space into a physical optics method, and calculating the profile data information of the target guided in the form of STL text to acquire the radar cross section of the target guided missile in the half space; and calculating echo signals of the target guided missile received by radar according to the radar cross section of the target guided missile in the half space. According to the invention, the echo signals of the target guided missile can be accurately figured out. The method can be used for AEWA radar system stimulation, and provides technical support for maximally playing the power of a cruise missile weapon system.

Description

Early warning plane radar detection target simulation method
Technical field
The invention belongs to wireless communication field, relate to detection system and the detection method thereof of antenna resonant frequency, can be used for the early warning plane radar detection simulation.
Technical background
Early warning plane is a kind of perfect electronic equipments such as Reconnaissance, Long Range radar, data processing, enemy and we's identification and communication and navigation, commander's control, electronic countermeasure that are equipped with, integrate early warning, Command, Control, Communication And Intelligence, for searching for, monitor and follow the tracks of operational support aircraft aerial and that naval target is also commanded, guiding one's own side aircraft is carried out combat duty, play the effect of moving radar station and Air Operation Center.Early warning plane can reduce the influence of earth curvature and various atural objects, thereby can increase hunting zone and the detection range of radar greatly, increases pre-warning time.For the cruise missile's target of flight in the home radar chain blind area, early warning plane can and give a warning in remote discovery target, becomes the prominent anti-biggest threat of cruise missile's target.The core equipment of early warning plane is exactly its airborne radar, and therefore, simulation study is to carry out the early warning plane operation to use one of important content of research to the airborne early warn ing radar detection of a target.
Present early warning plane radar detection target simulator mainly comprises: early warning plane modeling, the modeling of missile-radars characteristic, clutter Modeling Research and guided missile detection obtain echoed signal.Conventional method is that above-mentioned four modules are calculated emulation respectively separately, in missile-radars characteristic MBM, a lot of RCS computer-aided analyses and optimizing Design Software bag have been emerged based on various high frequency algorithms at present, though the accurate emulation of free space TV university size objectives RCS can be simulated and can be finished in real time to these software packages, but up to the present, the RCS high frequency computing method about semispace target guided missile yet there are no report.
Summary of the invention
The purpose of this invention is to provide a kind of early warning plane radar detection target simulation method, can not under high frequency, calculate the problem of the RCS of semispace target guided missile to solve prior art, obtain target guided missile echoed signal.
For achieving the above object: technical scheme of the present invention comprises the steps:
(1) obtains the target guided missile shape data information of STL text formatting
1a) utilize HFSS software to set up the model synoptic diagram of target guided missile according to target guided missile physical size;
1b) from the discrete data point at the model synoptic diagram edge of target guided missile, take out 4 wherein adjacent points arbitrarily; Appoint the quadratic sum minimum criteria of the difference on both sides to carry out triangle division according to triangle to these 4, selected wherein 3 and connect;
1c) for remaining discrete data point, repeating step 1b), be joined together the target guided missile synoptic diagram after obtaining handling up to all discrete data points;
1d) the target guided missile synoptic diagram after will handling is with the output of STL text formatting, and extracts the shape data information of target guided missile from output file, for the target guided missile provides shape data accurately in the calculating of the RCS of semispace;
(2) calculate the target guided missile at the RCS σ of semispace
To introduce in the physical optics method at the Green function of semispace, the shape data of STL text formatting target guided missile will be calculated, obtain the target guided missile at the scattering cross-section σ of semispace:
σ = lim R → ∞ 2 π R E s · e ^ r E o exp ( jkR )
That is: σ = ( lim R → ∞ 2 π R E s · e ^ r E o exp ( jkR ) ) 2 ,
In the formula, R is the round trip distance from the early warning plane radar to the target Missile Center, E oBe the incident electric field,
Figure BDA00002961429900023
Be the unit vector of target guided missile electric polarization direction, k is propagation constant, E sFor target guided missile scattering electric field, under far field condition, arrive the distance of far field observation point much larger than the characteristics of target guided missile size according to the target guided missile, obtain E sFollowing expression:
E s = - jω ∫ ∫ ∫ v G ‾ A ( r , r ′ ) · J ( r ′ ) d r ′ + k · s ^ ω ∫ ∫ ∫ v [ G x qe ∂ ∂ x J ( r ′ ) + G y qe ∂ ∂ y J ( r ′ ) + G z qe ∂ ∂ z J ( r ′ ) ] d r s ′
In the formula, ω is the incident wave frequency, and J (r ') is current density, and k is propagation constant,
Figure BDA00002961429900025
Be the unit vector of scattering wave polarised direction,
Figure BDA00002961429900026
Be the Green function at semispace:
G ‾ A ( r , r ′ ) = G A xx 0 0 0 G A yy 0 G A zx G A zy G A zz
In the formula
G A xx = μ 4 π ( - j 2 ) ∫ - ∞ + ∞ k ρ k z H 0 ( 2 ) ( k ρ ρ ) [ e - j k z | z - z ′ | + A h e e - j k z z ] d k ρ
G A zx = - μ 8 π ∂ ∂ x ∫ - ∞ + ∞ k ρ k z H 0 ( 2 ) ( k ρ ρ ) [ k z k ρ 2 ( C h e - A h e ) e - j k z z ] d k ρ
G A zz = μ 4 π ( - j 2 ) ∫ - ∞ + ∞ k ρ k z H 0 ( 2 ) ( k ρ ρ ) [ e - j k z | z - z ′ | + A v e e - j k z z ] d k ρ
G A yy = G A xx G A zy = G A zx ;
Figure BDA00002961429900036
Be the scalar potential function:
G x qe = - j 8 πϵ ∫ - ∞ + ∞ k ρ k z H 0 ( 2 ) ( k ρ ρ ) [ e - j k z | z - z ′ | + k 2 A h e + k z 2 C h e k ρ 2 e j k z z ] d k ρ
G y qe = G x qe
G z qe = 1 4 πϵ ( - j 2 ) ∫ - ∞ + ∞ k ρ k z H 0 ( 2 ) ( k ρ ρ ) [ e - j k z | z - z ′ | + C v e e - j k z z ] d k ρ
Wherein:
A h e = R TE e - j k z z ′ , C h e = - R TM e - j k z z ′ , A v e = R TE e - j k z z ′ , C v e = - R TM e - j k z z ′
In the formula, R TE, R TMBe fresnel reflection coefficient, k x, k y, k zFor along the propagation constant on x, y, the z direction,
Figure BDA000029614299000314
Figure BDA000029614299000315
Be the second class zeroth order Hankel function, ε is specific inductive capacity;
The physical optics of current density, J (r ') is approximate:
J ( r ′ ) = 2 n ^ × H 0
In the formula,
Figure BDA000029614299000317
Be the method vector of target guided missile, H 0The incident wave magnetic field intensity;
(3) calculate target guided missile echoed signal
Based on radar equation, consider radar system loss, radar operation wavelength, pulse peak power, radar system loss parameter situation, utilize the target guided missile at the scattering cross-section σ of semispace, calculate radar and receive target guided missile echo signal amplitude For:
Figure BDA00002961429900042
In the formula: P tBe the radar emission signal peak power, λ is the radar operation wavelength, and σ is that the target guided missile is long-pending in the RCS of semispace, G (α T, β T) be the radar antenna voltage gain, L is that radar emission receives composite loss, R is that radar arrives target Missile Center distance.
The present invention has following advantage:
1. the present invention utilizes HFSS software to set up the model synoptic diagram of target guided missile, utilize commercial femap software that the model synoptic diagram of target guided missile is carried out triangle division, target guided missile synoptic diagram after the triangle division is exported with the STL text formatting, and from output file, extract the shape data information of target guided missile, for the target guided missile provides shape data accurately in the calculating of the RCS of semispace.
2. the present invention introduces the semispace Green function in the physical optical method, and calculating is quick, simple, has solved prior art and can not calculate the target guided missile in the problem of the RCS of semispace under high frequency.
3. the present invention can calculate the echoed signal of target guided missile accurately according to the RCS of the target guided missile that calculates at semispace under high frequency.
Description of drawings
Fig. 1 is realization flow figure of the present invention.
Specific implementation process
With reference to Fig. 1, performing step of the present invention is as follows:
Step 1 is obtained the target guided missile shape data information of STL text formatting.
In order to calculate the target guided missile in the RCS of semispace, must obtain the shape data of the target guided missile of STL text formatting, its obtaining step is as follows:
1a) utilize software to set up the model synoptic diagram of target guided missile according to target guided missile physical size, this software comprises commercial HFSS software and commercial CST software;
1b) the model synoptic diagram of target guided missile is opened in triangle division software, software begins to carry out automatically triangle division, after software demarcation finishes, obtains the target guided missile synoptic diagram after the triangle division, and this software comprises commercial femap software and commercial Ansys software.
1c) the target guided missile synoptic diagram after the triangle division is exported with the STL text formatting, and from output file, extracted the shape data information of target guided missile, for the target guided missile provides shape data accurately in the calculating of the RCS of semispace;
Step 2 is calculated the target guided missile at the RCS σ of semispace.
In order to calculate target guided missile echoed signal, must at first calculate the target guided missile at the RCS σ of semispace, its calculation procedure is as follows:
2a) the current density, J (r ') of calculating target guided missile
Physical optical method gives simple clear and definite value for the current density, J (r ') on target guided missile surface, and current density is the incident magnetic field H 0Twice, so the calculating formula of current density, J (r ') as:
Figure BDA00002961429900051
In the formula
Figure BDA00002961429900052
Be the method vector of target guided missile, H 0Be incident magnetic field;
2b) the scattered field E of calculating target guided missile s
In physical optical method, carry out the integration stack by the shape data to STL text formatting target guided missile and calculate, obtain the scattered field E of target guided missile sAs follows:
Figure BDA00002961429900053
In the formula, ω is the incident wave frequency, and k is propagation constant,
Figure BDA00002961429900054
Be the unit vector of scattering wave polarised direction,
Figure BDA00002961429900055
Be the Green function at semispace:
G ‾ A ( r , r ′ ) = G A xx 0 0 0 G A yy 0 G A zx G A zy G A zz
In the formula
G A xx = μ 4 π ( - j 2 ) ∫ - ∞ + ∞ k ρ k z H 0 ( 2 ) ( k ρ ρ ) [ e - j k z | z - z ′ | + A h e e - j k z z ] d k ρ
G A zx = - μ 8 π ∂ ∂ x ∫ - ∞ + ∞ k ρ k z H 0 ( 2 ) ( k ρ ρ ) [ k z k ρ 2 ( C h e - A h e ) e - j k z z ] d k ρ
G A zz = μ 4 π ( - j 2 ) ∫ - ∞ + ∞ k ρ k z H 0 ( 2 ) ( k ρ ρ ) [ e - j k z | z - z ′ | + A v e e - j k z z ] d k ρ
G A yy = G A xx G A zy = G A zx ;
Be the scalar potential function:
G x qe = - j 8 πϵ ∫ - ∞ + ∞ k ρ k z H 0 ( 2 ) ( k ρ ρ ) [ e - j k z | z - z ′ | + k 2 A h e + k z 2 C h e k ρ 2 e j k z z ] d k ρ
G y qe = G x qe
G z qe = 1 4 πϵ ( - j 2 ) ∫ - ∞ + ∞ k ρ k z H 0 ( 2 ) ( k ρ ρ ) [ e - j k z | z - z ′ | + C v e e - j k z z ] d k ρ
Wherein:
A h e = R TE e - j k z z ′ , C h e = - R TM e - j k z z ′ , A v e = R TE e - j k z z ′ , C v e = - R TM e - j k z z ′
In the formula, R TE, R TMBe fresnel reflection coefficient, k x, k y, k zFor along the propagation constant on x, y, the z direction,
Figure BDA000029614299000612
Figure BDA000029614299000613
Be the second class zeroth order Hankel function, ε is specific inductive capacity;
2c) calculate the target guided missile at the RCS σ of semispace
Scattered field E according to the target guided missile sCalculate the target guided missile at the scattering cross-section σ of semispace:
σ = lim R → ∞ 2 π R E s · e ^ r E o exp ( jkR )
That is: σ = ( lim R → ∞ 2 π R E s · e ^ r E o exp ( jkR ) ) 2 ,
In the formula, R is the round trip distance from the early warning plane radar to the target Missile Center, E oBe the incident electric field,
Figure BDA000029614299000616
Be the unit vector of target guided missile electric polarization direction, k is propagation constant.
Step 3 is calculated target guided missile echoed signal
Based on radar equation, consider radar system loss, radar operation wavelength, pulse peak power, radar system loss parameter situation, utilize target guided missile that step 2 calculates at the scattering cross-section σ of semispace, calculate radar and receive target guided missile echo signal amplitude
Figure BDA00002961429900071
For:
Figure BDA00002961429900072
In the formula: P tBe the radar emission signal peak power, λ is the radar operation wavelength, and σ is that the target guided missile is long-pending in the RCS of semispace, G (α T, β T) be the radar antenna voltage gain, L is that radar emission receives composite loss, R is that radar arrives target Missile Center distance.
Calculate target guided missile echoed signal, be the early warning plane radar and receive target guided missile echoed signal, finish early warning plane radar detection target simulation process.

Claims (2)

1. an early warning plane radar detection target simulation method comprises the steps:
(1) obtains the target guided missile shape data information of STL text formatting
1a) utilize HFSS software to set up the model synoptic diagram of target guided missile according to target guided missile physical size;
1b) from the discrete data point at the model synoptic diagram edge of target guided missile, take out 4 wherein adjacent points arbitrarily; Appoint the quadratic sum minimum criteria of the difference on both sides to carry out triangle division according to triangle to these 4, selected wherein 3 and connect;
1c) for remaining discrete data point, repeating step 1b), be joined together the target guided missile synoptic diagram after obtaining handling up to all discrete data points;
1d) the target guided missile synoptic diagram after will handling is with the output of STL text formatting, and extracts the shape data information of target guided missile from output file, for the target guided missile provides shape data accurately in the calculating of the RCS of semispace;
(2) calculate the target guided missile at the RCS σ of semispace
To introduce in the physical optics method at the Green function of semispace, the shape data of STL text formatting target guided missile will be calculated, obtain the target guided missile at the scattering cross-section σ of semispace:
σ = lim R → ∞ 2 π R E s · e ^ r E o exp ( jkR )
That is: σ = ( lim R → ∞ 2 π R E s · e ^ r E o exp ( jkR ) ) 2 ,
In the formula, R is the round trip distance from the early warning plane radar to the target Missile Center, E oBe the incident electric field, Be the unit vector of target guided missile electric polarization direction, k is propagation constant, E sFor target guided missile scattering electric field, under far field condition, arrive the distance of far field observation point much larger than the characteristics of target guided missile size according to the target guided missile, obtain E sFollowing expression:
E s = - jω ∫ ∫ ∫ v G ‾ A ( r , r ′ ) · J ( r ′ ) d r ′ + k · s ^ ω ∫ ∫ ∫ v [ G x qe ∂ ∂ x J ( r ′ ) + G y qe ∂ ∂ y J ( r ′ ) + G z qe ∂ ∂ z J ( r ′ ) ] dr ′ s
In the formula, ω is the incident wave frequency, and J (r ') is current density, and k is propagation constant, Be the unit vector of scattering wave polarised direction,
Figure FDA00002961429800016
Be the Green function at semispace:
G ‾ A ( r , r ′ ) = G A xx 0 0 0 G A yy 0 G A zx G A zy G A zz
In the formula
G A xx = μ 4 π ( - j 2 ) ∫ - ∞ + ∞ k ρ k z H 0 ( 2 ) ( k ρ ρ ) [ e - j k z | z - z ′ | + A h e e - j k z z ] d k ρ
G A zx = - μ 8 π ∂ ∂ x ∫ - ∞ + ∞ k ρ k z H 0 ( 2 ) ( k ρ ρ ) [ k z k ρ 2 ( C h e - A h e ) e - j k z z ] d k ρ
G A zz = μ 4 π ( - j 2 ) ∫ - ∞ + ∞ k ρ k z H 0 ( 2 ) ( k ρ ρ ) [ e - j k z | z - z ′ | + A v e e - j k z z ] d k ρ
G A yy = G A xx G A zy = G A zx ;
Figure FDA00002961429800027
Be the scalar potential function:
G x qe = - j 8 πϵ ∫ - ∞ + ∞ k ρ k z H 0 ( 2 ) ( k ρ ρ ) [ e - j k z | z - z ′ | + k 2 A h e + k z 2 C h e k ρ 2 e j k z z ] d k ρ
G y qe = G x qe
G z qe = 1 4 πϵ ( - j 2 ) ∫ - ∞ + ∞ k ρ k z H 0 ( 2 ) ( k ρ ρ ) [ e - j k z | z - z ′ | + C v e e - j k z z ] d k ρ
Wherein:
A h e = R TE e - j k z z ′ , C h e = - R TM e - j k z z ′ , A v e = R TE e - j k z z ′ , C v e = - R TM e - j k z z ′
In the formula, R TE, R TMBe fresnel reflection coefficient, k x, k y, k zFor along the propagation constant on x, y, the z direction,
Figure FDA000029614298000215
Figure FDA000029614298000216
Be the second class zeroth order Hankel function, ε is specific inductive capacity;
The physical optics of current density, J (r ') is approximate:
J ( r ′ ) = 2 n ^ × H 0
In the formula, Be the method vector of target guided missile, H 0The incident wave magnetic field intensity;
(3) calculate target guided missile echoed signal
Based on radar equation, consider radar system loss, radar operation wavelength, pulse peak power, radar system loss parameter situation, utilize the target guided missile at the scattering cross-section σ of semispace, calculate radar and receive target guided missile echo signal amplitude
Figure FDA000029614298000219
For:
In the formula: P tBe the radar emission signal peak power, λ is the radar operation wavelength, and σ is that the target guided missile is long-pending in the RCS of semispace, G (α T, β T) be the radar antenna voltage gain, L is that radar emission receives composite loss, R is that radar arrives target Missile Center distance.
2. early warning plane radar detection target simulation method according to claim 1, it is characterized in that: appointing the quadratic sum minimum criteria of the difference on both sides to carry out triangle division according to triangle to adjacent 4 described step 1b), is to utilize commercial femap software to carry out triangle division.
CN2013100982335A 2013-03-25 2013-03-25 Method for AEWA (Airborne Early Warning Aircraft) radar detection target simulation Pending CN103186697A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN2013100982335A CN103186697A (en) 2013-03-25 2013-03-25 Method for AEWA (Airborne Early Warning Aircraft) radar detection target simulation

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN2013100982335A CN103186697A (en) 2013-03-25 2013-03-25 Method for AEWA (Airborne Early Warning Aircraft) radar detection target simulation

Publications (1)

Publication Number Publication Date
CN103186697A true CN103186697A (en) 2013-07-03

Family

ID=48677864

Family Applications (1)

Application Number Title Priority Date Filing Date
CN2013100982335A Pending CN103186697A (en) 2013-03-25 2013-03-25 Method for AEWA (Airborne Early Warning Aircraft) radar detection target simulation

Country Status (1)

Country Link
CN (1) CN103186697A (en)

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP0579143A1 (en) * 1992-07-13 1994-01-19 Hughes Aircraft Company Method and apparatus for missile interface testing

Patent Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP0579143A1 (en) * 1992-07-13 1994-01-19 Hughes Aircraft Company Method and apparatus for missile interface testing

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
李晓峰等: "半空间电大涂敷目标散射的高频分析方法", 《物理学报》 *
李磊等: "一种探测低飞目标的PD 雷达仿真系统", 《系统工程与电子技术》 *
陈博韬: "低空目标与环境符合电磁散射特性研究", 《信息科技辑》 *

Similar Documents

Publication Publication Date Title
CN104077601B (en) A kind of method that based target integrated identification is carried out using different types of information
CN102981064B (en) Aircraft external radio frequency electromagnetic environment prediction method and prediction system
CN109901159A (en) A kind of metre wave radar target identification method based on deep learning
CN109613530B (en) Control method for multi-source information fusion of low-small slow air target
Arvas et al. A tutorial on the method of moments [testing ourselves]
CN104502911A (en) Wall parameter estimation method of through-wall imaging radar
Shen et al. Two‐step locating method for aircraft wake vortices based on Gabor filter and velocity range distribution
Ye et al. Gaf representation of millimeter wave drone rcs and drone classification method based on deep fusion network using resnet
CN106446733A (en) Label reception power prediction method of ultra high frequency (UHF) radio frequency identification (RFID) electronic toll collection (ETC) applications
Özdemir Radar cross section analysis of unmanned aerial vehicles using predics
CN103852648B (en) The method obtaining spatial electromagnetic intensity data
CN104101869B (en) A kind of ground wave radar moving target simulation modeling method under polar coordinate
Pandey Modeling and simulation of chaff cloud with random orientation and distribution
CN103675777B (en) Based on airborne radar clutter analogy method and the device of fitting process
CN104732050A (en) Method for estimating distribution of electromagnetism in flying targets made of carbon fiber materials under lightning pulses
CN103186697A (en) Method for AEWA (Airborne Early Warning Aircraft) radar detection target simulation
Ptak et al. Aircraft classification based on radar cross section of long-range trajectories
Mansukhani et al. Rcs based target classification using deep learning methods
Honda et al. An estimation algorithm of scattered powers caused by a moving aircraft
Zheng et al. Research on channel optimization of Ads‐B aviation target surveillance radar based on improved filtering algorithm
Weinmann Accurate prediction of EM scattering by wind turbines
Hashimov et al. ASSESSMENT OF DEAD ZONE OF JOINTLY OPERATING RADARS
Heidar et al. Simulation and reduction of radar cross section the unmanned aerial vehicle in X-band using Shaping technique
Oikawa et al. Analysis of lightning electromagnetic field on large‐scale terrain model using three‐dimensional MW‐FDTD parallel computation
Petchmak et al. Comparing Coverage Models of Radar: Evaluating Viewshed Against Ground Truth Detection

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C02 Deemed withdrawal of patent application after publication (patent law 2001)
WD01 Invention patent application deemed withdrawn after publication

Application publication date: 20130703