CN103870654A - Electromagnetic scattering simulation method based on parallel moment method and physical optics mixing - Google Patents

Electromagnetic scattering simulation method based on parallel moment method and physical optics mixing Download PDF

Info

Publication number
CN103870654A
CN103870654A CN201410115110.2A CN201410115110A CN103870654A CN 103870654 A CN103870654 A CN 103870654A CN 201410115110 A CN201410115110 A CN 201410115110A CN 103870654 A CN103870654 A CN 103870654A
Authority
CN
China
Prior art keywords
radar target
uneven surface
composite model
field
electric current
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
CN201410115110.2A
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 CN201410115110.2A priority Critical patent/CN103870654A/en
Publication of CN103870654A publication Critical patent/CN103870654A/en
Pending legal-status Critical Current

Links

Images

Landscapes

  • Radar Systems Or Details Thereof (AREA)

Abstract

The invention discloses an electromagnetic scattering simulation method based on parallel moment method and physical optics mixing to mainly solve the problem that when electrically large models are processed in the prior art, the simulation efficiency of electromagnetic scattering and the precision are low. The electromagnetic scattering simulation method comprises the steps that a rough surface power spectrum and rough surface parameters are input, and a rough surface is simulated based on the Monte Carlo method; geometric modeling is carried out on a radar target through 3dsMax software, and the subdivided radar target is output; a combined model is formed by combining the rough surface and the radar target, the rough surface is divided into a physical optics area, and the radar target is divided into a parallel moment method area; induced currents on the rough surface and induced currents on the surface of the radar target are obtained; a far scattered field of the combined model is calculated based on the induced currents on the surface of the combined model, a radar cross section of the combined model is obtained, and simulation between the rough surface and the combined electromagnetic scattering of the radar target is achieved. According to the electromagnetic scattering simulation method, the memory consumption is little, the simulation efficiency is high, and the electromagnetic scattering simulation method can be used in the study on the rough surface and the radar target combined electromagnetic scattering character in the field of national defense and the civil using.

Description

The electromagnetic scattering emulation mode of mixing with physical optics based on Parallel MoM
Technical field
The invention belongs to radar Electromagnetic Simulation technical field, relate to a kind of electromagnetic scattering emulation mode, for obtaining the composite radar scattering cross-section of target and rough surface.
Background technology
Along with the fast development of Radar Technology, the composite electromagnetic scattering research of target and rough surface has great importance in theoretical analysis and practical application.In the time that bank base, airborne or spaceborne radar carry out electromagnetic surveying to the target in all kinds of environment such as ocean, dark woods, desert, in the echoed signal that radar receives, not only comprise target information, also comprise the electromagnetic scattering signal of rugged face background.By analyzing these echoed signals, and then can draw some characteristic of surveyed environment and target.Thereby the radar scattering characteristic research of target and rough surface all has significant learning value and application prospect widely at national defence and civil area.
In decades, many Electromagnetic Simulation technology are proposed the composite electromagnetic scattering problem in order to process sea, ground and target by scholar, are broadly divided into approximation method and numerical method in the past.The advantage of approximation method is that EMS memory occupation is few, analysis speed is fast, but approximation method is all approximate based on specific physics, and narrow application range and precision are lower.Than approximation method, numerical method can keep higher simulation accuracy, but its EMS memory occupation is many, and analysis speed is slow.For TV university size uneven surface and the compound single station of complex target Electromagnetic Simulation problem, the configuration that numerical method has been no longer ordinary individual's computing machine to the demand of internal memory can meet, and simulation velocity is too slow, therefore these class methods are difficult to be generalized in the composite electromagnetic scattering emulation of actual electrical large scale uneven surface and complex target.Mixed method can be in conjunction with the advantage of the high-level efficiency of approximation method and numerical method accuracy, when conventional hybrid method is carried out emulation to the composite electromagnetic scattering of uneven surface and target, often adopt plane wave as incident source, can form artificial reflection at time-limited uneven surface two ends like this, affect the accuracy of result, and conventional hybrid method still can not meet engineering application requirements completely for the simulation efficiency of TV university size uneven surface and complex target composite electromagnetic scattering.
Summary of the invention
The object of the invention is to the deficiency for above-mentioned prior art, propose a kind of Electromagnetic Simulation method of mixing with physical optics based on Parallel MoM, guaranteeing, under the prerequisite of computational solution precision, to reduce unknown quantity number and computing time, reduce the memory requirements of emulation, improve simulation efficiency.
The technical scheme that realizes the object of the invention, comprises the steps:
(1) power spectral density function and the uneven surface parameter of input uneven surface, obtains the position coordinates of N × M sampled point on uneven surface by Monte Carlo method:
r -n+1, -m+1=(x -n+1,y -m+1,z -n+1, -m+1),...,
r 0, 0=(x 0,y 0z 0, 0),...,r n, m=(x n,y m?z n, m),
Wherein n=N2, m=M2, N represents x direction hits, and M represents y direction hits, and subscript represents the sequence number of uneven surface up-sampling point;
(2) adopt triangular facet tuple subdivision uneven surface, obtain the tri patch of series of discrete;
(3) utilize 3dsMax software to carry out Geometric Modeling to radar target, and export the radar target model after subdivision;
(4) with uneven surface and radar target composition composite model, and the uneven surface in this composite model is divided into physical optics region, the radar target in this composite model is divided into Parallel MoM region;
(5) according to the coupling between two regions, calculate the surface induction electric current of composite model, basic step is as follows:
5.1) irradiate composite model with conical wave, and the initial induction current of establishing radar target surface is zero, uneven surface is encouraged by conical wave only, approximate according to physical optics, obtains the surface induction electric current of uneven surface
Figure BDA0000482513070000021
5.2) use uneven surface induction current
Figure BDA0000482513070000022
the scattered field exciting encourages radar target together with conical wave, according to boundary condition, sets up the improved Electric Field Integral Equation on radar target surface, obtains the surface induction electric current of radar target
Figure BDA0000482513070000023
5.3) with radar target surface induction electric current
Figure BDA0000482513070000024
the scattered field exciting encourages uneven surface together with conical wave, approximate according to physical optics, obtains the surface induction electric current of uneven surface
Figure BDA0000482513070000025
5.4) repeating step 5.2) and step 5.3) iterative process, until
Figure BDA0000482513070000026
when <err, the surface induction electric current of whole composite model reaches stable, and iterative process stops, and obtains the final induction current in surface of composite model, and wherein i+1 is iterative steps, and err is the iteration error of setting;
(6) according to the induction current of radar target and uneven surface and Huygens' principle, calculate the far field scattered field of composite model, finally obtain the RCS of composite model.
Compared with prior art, tool has the following advantages in the present invention:
1) the present invention adopts conical wave E ias the irradiation wave source of radar, avoid the artificial reflection of uneven surface sample area edge, improve the accuracy of this emulation mode;
2) the present invention, due to TV university size uneven surface is adopted to physical optics approximate treatment, has greatly improved counting yield, has reduced memory requirements, and the target with labyrinth is calculated with Parallel MoM, has guaranteed the accuracy of result of calculation;
3) the present invention has adopted parallel computing calculating conductor target surface induction current, has reduced the time of calculating, the electric size of having expanded calculating object, is applicable to solving of TV university size uneven surface and complex target composite electromagnetic scattering problem.
Accompanying drawing explanation
Fig. 1 is realization flow figure of the present invention;
Fig. 2 is uneven surface and the spheroid target combination scattering model instance figure that the present invention sets up;
Fig. 3 is the uneven surface set up of the present invention and section and the parameter schematic diagram of spheroid target combination scattering model;
Fig. 4 is the composite radar by the present invention and existing multilevel fast multipole algorithm emulation uneven surface and spheroid target
Scattering cross-section comparison diagram;
Fig. 5 is with the different wind speed Rough Sea Surfaces of emulation of the present invention and the composite radar scattering cross-section figure of spheroid target;
Fig. 6 is with the different wind territory Rough Sea Surfaces of emulation of the present invention and the composite radar scattering cross-section figure of spheroid target;
Fig. 7 is different size target and the Rough Sea Surfaces composite radar scattering cross-section figure with emulation of the present invention;
Fig. 8 is diverse location target and the Rough Sea Surfaces composite radar scattering cross-section figure with emulation of the present invention.
Embodiment
Below in conjunction with accompanying drawing, embodiments of the invention and effect are described in further detail.
With reference to Fig. 1, performing step of the present invention is as follows:
Step 1: generate the required combination scattering model of emulation.
(1.1) obtain by experiment uneven surface parameter and the power spectral density function S (k) thereof on emulation ground or the sea of wanting, for example: be S (k)=2 π δ exp (k for ground power spectral density function 2/ δ 2)/l 2, be S (k)=δ for sea power spectral density function 2l/[π (1+k 2l 2)], k is space wave number, and l is persistence length, and δ is root-mean-square height;
(1.2) produce N random number that meets Gaussian distribution, N is greater than 0 even number;
(1.3) frequency f of input incident wave, establishes sampling interval △ x=3.0 × 10 8/ (10f), get N discrete sampling point x -N/ 2+1..., x -1, x 0, x 1..., x n/ 2, make x 0=0, x -N/ 2+1=(N/2+1). △ x, x -1=-△ x, x 1=△ x ..., x n/ 2=N/2. △ x;
(1.4) according to the uneven surface parameter and the power spectral density function S (k) that obtain, the random number respectively step (1.2) being produced in various discrete sample point is revised, and obtains respectively N revised random number at various discrete point place;
(1.5) horizontal ordinate using N the discrete sampling point producing in step (1.3) as uneven surface, superposes various discrete sample point N revised random number respectively, obtains the ordinate z that each sampled point is corresponding -N/ 2+1..., z -1, z 0, z 1..., z n/ 2, the uneven surface of the generation emulation of wanting;
(1.6) utilize 3ds Max software to carry out Geometric Modeling to radar target, by triangle bin subdivision for radar target model, derive the radar target after subdivision;
(1.7) by the radar target composition combination scattering model of deriving in the uneven surface in step (1.5) and step (1.6), as shown in Figure 2, and the uneven surface in combination scattering model is divided into physical optics region, the radar target in combination scattering model is divided into method of moment region.
Step 2: obtain uneven surface surface and radar target surface current.
(2.1) irradiate composite model with conical wave, the electric field of conical wave
Figure BDA0000482513070000041
and magnetic field
Figure BDA0000482513070000042
be expressed as follows respectively:
Figure BDA0000482513070000044
Wherein,
Figure BDA0000482513070000045
the sampled point on composite model surface,
Figure BDA0000482513070000046
j is complex unit,
Figure BDA0000482513070000047
for space wave number,
Figure BDA0000482513070000048
it is the unit vector of space wave number; η is composite model space outerpace wave impedance,
Figure BDA00004825130700000410
μ afor the relative magnetic specific inductive capacity of composite model space outerpace, ε afor the relatively electric specific inductive capacity of composite model space outerpace, g is the conical wave factor, θ ifor conical wave incident angle,
Figure BDA00004825130700000411
Figure BDA00004825130700000418
for conical wave incident orientation angle;
(2.2) establish the initial induction current in radar target surface
Figure BDA00004825130700000412
be zero, uneven surface is encouraged by conical wave only, approximate according to physical optics, obtains uneven surface surface induction electric current
Figure BDA00004825130700000413
Figure BDA00004825130700000414
Wherein,
Figure BDA00004825130700000415
for the unit normal vector on uneven surface surface,
Figure BDA00004825130700000416
for the sampled point on uneven surface surface,
Figure BDA00004825130700000417
the magnetic field of the scattered field that expression radar target surface induction electric current excites:
Figure BDA0000482513070000051
In formula
Figure BDA00004825130700000523
represent radar target surface sample point coordinate to ask curl, S pMOMfor radar target induction current integral area,
Figure BDA0000482513070000052
green function in three dimensions:
Figure BDA0000482513070000053
for the position of putting on composite model,
Figure BDA0000482513070000055
for the position of far zone field point;
(2.3) with uneven surface surface induction electric current
Figure BDA0000482513070000056
the scattered field exciting and conical wave E iencourage together radar target, again calculate the surface induction electric current of radar target according to boundary condition
Figure BDA0000482513070000057
(2.4) repeating step (2.2)-(2.3), until
Figure BDA0000482513070000058
<err, whole region induction current tends towards stability, and iterative process stops, and obtains the final induction current in surface of composite model, and wherein i is iterative steps, and err is the iteration error of setting.
Step 3: the far field scattered field and the RCS that obtain composite model.
(3.1) calculate radar target surface induction electric current
Figure BDA0000482513070000059
far field scattered field:
Figure BDA00004825130700000510
Wherein, k is space wave number, and η is space wave impedance,
Figure BDA00004825130700000521
for unit dyad,
Figure BDA00004825130700000511
for dyad,
Figure BDA00004825130700000512
for far zone field
The unit vector of point, r is the distance of far zone field point from true origin,
Figure BDA00004825130700000513
for the position of putting on radar target,
S pMOMfor radar target surface induction current integration area;
(3.2) calculate uneven surface surface induction electric current far field scattered field:
Figure BDA00004825130700000515
Wherein, for the position of uneven surface surface point, S pOfor radar target surface induction current integration area;
(3.3), according to the result of step (3.1) and step (3.2), calculate the far field scattered field of composite model
Figure BDA00004825130700000517
(3.4) by the far field scattered field of composite model transform to cylindrical coordinate by rectangular coordinate system, obtain cylindrical coordinate θ direction coordinate components
Figure BDA00004825130700000519
with
Figure BDA00004825130700000522
direction coordinate components
Figure BDA00004825130700000520
and then the RCS that obtains composite model is:
Wherein,
Figure BDA0000482513070000062
for the electric field of conical wave, r is the distance of far zone field point from true origin.
Effect of the present invention can further illustrate by following test examples:
1. test simulation condition
The radar incident frequency f=300MHz using in emulation experiment, incident wave wavelength X=1.0m, conical wave factor g=25.6, uneven surface is selected the Random ocean-like rough surfaces based on JONSWAP spectrum, its power spectral density function S ( &omega; , &theta; ) = ag c 2 1 &omega; 5 exp [ - 5 4 ( &omega; 0 &omega; ) 4 ] . &gamma; exp [ - ( &omega; - &omega; 0 ) 2 2 &sigma; 2 &omega; 0 2 ] . G ( &omega; , &theta; ) , Wherein zero dimension constant
Figure BDA0000482513070000064
nondimensional fetch its span is 10 -1~10 5, x Wei Feng district, U 10for the wind speed of 10m At The Height, gc=9.81m/s 2, angular frequency=2pf, peak frequency
Figure BDA0000482513070000066
γ is the peak rising factor, and observed reading is 1.5 6, and mean value is 3.3; σ is peak shape parameter, ω≤ω 0time σ=0.07, ω > ω 0time
Figure BDA0000482513070000067
wind direction
Figure BDA00004825130700000610
direction of observation θ=30.Due to the randomness of uneven surface profile, thereby the mean value of getting 30 samples obtains stable monostatic radar scattering parameter.
2. test simulation example and interpretation of result
Emulation experiment 1:
Incident wave frequency f=300MHz, incidence angle θ i=0, incident orientation angle region, sea S=25.6m*25.6m, sea wind speed U=2m/s, wind district x=5km.Under the condition that is 128*128 at sea sampling number, generate 100 sea samples, after these 100 sea samples are averaged, obtained the needed sea of this emulation experiment.Fig. 3 has provided this sea and sphere centre coordinate is (0,0,6), the composite model of the conducting sphere of r=2m.
The RCS of composite model in Fig. 3 is calculated by the inventive method, and the RCS of the composite model calculating with the quick multistage submethod MLFMM of existing multilayer contrasts, its result is as Fig. 4.
As can be seen from Figure 4, the RCS figure of the composite model that the quick multistage submethod of the RCS figure of the composite model that the inventive method is calculated and multilayer calculates is basically identical in whole scattering region, thereby has proved the accuracy of the RCS of the inventive method calculating composite model.
The time and the memory source that during by the RCS of composite model in the inventive method calculating chart 3, consume, and the time and the memory source that consume while calculating the RCS of composite model with the quick multistage submethod of existing multilayer contrast, its result table 1.
Figure BDA0000482513070000071
As can be seen from Table 1, the memory source simultaneously consuming also reduces to some extent far less than the quick multistage submethod of existing multilayer the time that the inventive method consumes while calculating the RCS of composite model.
Emulation experiment 2, according to Fig. 3 model, carries out emulation to the RCS of different size target and Rough Sea Surfaces composite model with the present invention, and result is as Fig. 5.
As can be seen from Figure 5, the inventive method can be calculated the RCS of different size target and Rough Sea Surfaces composite model.
Emulation experiment 3, according to Fig. 3 model, carries out emulation to the RCS of different Rough Sea Surfaces and target composite model with the present invention, and result is as Fig. 6.
As can be seen from Figure 6, the present invention can calculate the RCS of different Rough Sea Surfaces and target composite model.
Emulation experiment 4, according to Fig. 3 model, carries out emulation to the RCS figure cross section of different size target and Rough Sea Surfaces composite model with the present invention, and result is as Fig. 7.
As can be seen from Figure 7, the inventive method can be calculated the RCS of different size target and Rough Sea Surfaces composite model.
Emulation experiment 5, according to Fig. 3 model, carries out emulation to the RCS of diverse location target and Rough Sea Surfaces composite model with the present invention, and result is as Fig. 8.
As can be seen from Figure 8, the present invention can calculate the RCS of diverse location target and Rough Sea Surfaces composite model.

Claims (6)

1. an electromagnetic scattering emulation mode of mixing with physical optics based on Parallel MoM, comprises the steps:
(1) power spectral density function and the uneven surface parameter of input uneven surface, obtains the position coordinates of N × M sampled point on uneven surface by Monte Carlo method:
r -n+1, -m+1=(x -n+1,y -m+1,z -n+1, -m+1),...,
r 0, 0=(x 0,y 0z 0, 0),...,r n, m=(x n,y m?z n, m),
Wherein n=N2, m=M2, N represents x direction hits, and M represents y direction hits, and subscript represents the sequence number of uneven surface up-sampling point;
(2) adopt triangular facet tuple subdivision uneven surface, obtain the tri patch of series of discrete;
(3) utilize 3dsMax software to carry out Geometric Modeling to radar target, and export the radar target after subdivision;
(4) with uneven surface and radar target composition composite model, and the uneven surface in composite model is divided into physical optics region, radar target is divided into Parallel MoM region;
(5) according to the coupling between two regions, calculate the surface induction electric current of composite model, basic step is as follows:
5.1) irradiate composite model with conical wave, the initial induction current of establishing radar target surface is zero, and uneven surface is encouraged by conical wave only, the approximate and boundary condition according to physical optics, the surface induction electric current of acquisition uneven surface
Figure FDA0000482513060000011
5.2) use uneven surface induction current
Figure FDA0000482513060000012
the scattered field exciting encourages radar target together with conical wave, sets up the improved Electric Field Integral Equation on radar target surface according to boundary condition, obtains the surface induction electric current of radar target
Figure FDA0000482513060000013
5.3) with radar target surface induction electric current
Figure FDA0000482513060000014
the scattered field exciting encourages uneven surface together with conical wave, the approximate and boundary condition according to physical optics, the surface induction electric current of acquisition uneven surface
5.4) repeating step 5.2) and step 5.3) iterative process, until
Figure FDA0000482513060000016
when <err, the surface induction electric current of whole composite model reaches stable, and iterative process stops, and obtains the final induction current in surface of composite model, wherein i, and i+1 is iterative steps, err is the iteration error of setting;
(6) according to the induction current of radar target and uneven surface and Huygens' principle, calculate the far field scattered field of composite model, finally obtain the RCS of composite model.
2. electromagnetic scattering emulation mode of mixing with physical optics based on Parallel MoM according to claim 1, step (1) described according to ground or sea parameter and power spectral density function, by Monte-carlo Simulation uneven surface, concrete steps are as follows:
(1a) obtain by experiment parameter and the power spectral density function on emulation rough earth or the sea of wanting, and choose independently Gauss's harmonic wave of a series of amplitudes;
(1b), according to uneven surface parameter and power spectral density function, selected Gauss's harmonic wave is carried out to Modulation and Amplitude Modulation, then the Gauss's harmonic wave after modulation is carried out to the uneven surface of Fourier transform generation emulation.
3. electromagnetic scattering emulation mode of mixing with physical optics based on Parallel MoM according to claim 1, step (2) adopts triangular facet tuple subdivision uneven surface, obtain the triangle bin of series of discrete, on uneven surface, search for along x direction, mutual three adjacent sampled points on uneven surface are connected into a tri patch, each dough sheet is pressed to line number, sequence number, summit, normal vector and the seamed edge information of output dough sheet.
4. electromagnetic scattering emulation mode of mixing with physical optics based on Parallel MoM according to claim 1, step 5.1) the described approximate and boundary condition according to physical optics, obtain uneven surface surface induction electric current
Figure FDA0000482513060000021
carry out as follows:
(5.1a) irradiate composite model with conical wave, the electric field of conical wave
Figure FDA0000482513060000022
and magnetic field be expressed as follows respectively:
Figure FDA0000482513060000024
Figure FDA0000482513060000025
Wherein,
Figure FDA0000482513060000026
the sampled point on composite model surface,
Figure FDA0000482513060000027
j is complex unit,
Figure FDA0000482513060000028
for space wave number,
Figure FDA0000482513060000029
it is the unit vector of space wave number;
Figure FDA00004825130600000210
η is composite model space outerpace wave impedance,
Figure FDA00004825130600000211
μ afor the relative magnetic specific inductive capacity of composite model space outerpace, ε afor the relatively electric specific inductive capacity of composite model space outerpace; K is space wave number, and g is the conical wave factor, θ ifor conical wave incident angle,
Figure FDA00004825130600000212
Figure FDA00004825130600000213
for conical wave incident orientation angle;
(5.1b) establish the initial induction current in radar target surface
Figure FDA0000482513060000031
be zero, uneven surface is encouraged by conical wave only, and according to physical optics, approximate and boundary condition, obtains uneven surface surface induction electric current
Figure FDA0000482513060000032
Figure FDA0000482513060000033
Wherein,
Figure FDA0000482513060000034
for the unit normal vector on uneven surface surface,
Figure FDA0000482513060000035
for the sampled point on uneven surface surface,
Figure FDA0000482513060000036
the magnetic field of the scattered field that expression radar target surface induction electric current excites:
Figure FDA0000482513060000037
Wherein,
Figure FDA00004825130600000322
represent radar target surface sample point coordinate to ask curl, for the surface induction electric current of radar target, S pMOMfor radar target induction current integral area.
Figure FDA0000482513060000039
green function in three dimensions:
Figure FDA00004825130600000310
Wherein,
Figure FDA00004825130600000311
for the position of putting on composite model,
Figure FDA00004825130600000312
for the position of far zone field point.
5. electromagnetic scattering emulation mode of mixing with physical optics based on Parallel MoM according to claim 1, wherein step 5.2) the described improved Electric Field Integral Equation of setting up radar target surface according to boundary condition, obtain the surface induction electric current of radar target carry out as follows:
(5.2a) set up the improved Electric Field Integral Equation on radar target surface according to boundary condition:
Figure FDA00004825130600000314
Wherein,
Figure FDA00004825130600000315
for the electric field of conical wave, S pOfor uneven surface induction current integral area, subscript tan represents the tangential direction of electric field on radar target surface,
Figure FDA00004825130600000316
for uneven surface surface induction electric current
Figure FDA00004825130600000317
the electric field of the scattered field exciting at radar target:
Wherein,
Figure FDA00004825130600000321
represent uneven surface surface sample point coordinate to differentiate;
(5.2b) select one group of complete RWG basis function, by the surface current of radar target
Figure FDA00004825130600000319
launch:
Figure FDA00004825130600000320
In formula
Figure FDA0000482513060000041
for RWG basis function, β kfor expansion coefficient;
(5.2c) by the surface current of the radar target after launching
Figure FDA0000482513060000042
radar target surface field integral equation in substitution (5.2a):
Figure FDA0000482513060000043
(5.2d) select one group of weight function
Figure FDA0000482513060000044
and do inner product with (5.2c) middle improved Electric Field Integral Equation and carry out inspection by sampling, obtain discrete system of linear equations, adopt parallel method of conjugate gradient to solve system of linear equations, obtain radar target surface induction electric current expansion coefficient β k, finally obtain the surface induction electric current of radar target
Figure FDA0000482513060000046
6. electromagnetic scattering emulation mode of mixing with physical optics based on Parallel MoM according to claim 1, in wherein said step (6), calculate the far field scattered field of composite model, obtain the RCS of composite model, carry out as follows:
(6.1) calculate radar target surface induction electric current
Figure FDA0000482513060000047
far field scattered field:
Figure FDA0000482513060000048
In formula, k is wave number, and I is unit dyad, for dyad, for far zone field point unit vector, r is the distance of far zone field point from true origin;
(6.2) calculate uneven surface surface induction electric current
Figure FDA00004825130600000411
far field scattered field:
Figure FDA00004825130600000412
(6.3), according to the result of step (6.1) and step (6.2), calculate the far field scattered field of composite model
Figure FDA00004825130600000413
(6.4) by the far field scattered field of composite model
Figure FDA00004825130600000414
transform to cylindrical coordinate by rectangular coordinate system, obtain cylindrical coordinate θ durection component with durection component
Figure FDA00004825130600000416
and then the RCS that obtains composite model is:
Figure FDA00004825130600000417
CN201410115110.2A 2014-03-26 2014-03-26 Electromagnetic scattering simulation method based on parallel moment method and physical optics mixing Pending CN103870654A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201410115110.2A CN103870654A (en) 2014-03-26 2014-03-26 Electromagnetic scattering simulation method based on parallel moment method and physical optics mixing

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201410115110.2A CN103870654A (en) 2014-03-26 2014-03-26 Electromagnetic scattering simulation method based on parallel moment method and physical optics mixing

Publications (1)

Publication Number Publication Date
CN103870654A true CN103870654A (en) 2014-06-18

Family

ID=50909181

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201410115110.2A Pending CN103870654A (en) 2014-03-26 2014-03-26 Electromagnetic scattering simulation method based on parallel moment method and physical optics mixing

Country Status (1)

Country Link
CN (1) CN103870654A (en)

Cited By (28)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104317984A (en) * 2014-09-09 2015-01-28 中国舰船研究设计中心 Ship electromagnetic scattering prediction method and system based on sub-domain modeling
CN104992035A (en) * 2015-07-28 2015-10-21 中国人民解放军国防科学技术大学 Quick calculation method for terahertz band surface rough target electromagnetic scattering
CN106156394A (en) * 2015-04-22 2016-11-23 南京理工大学 Electromagnetic property extracting method based on explicit difference scheme
CN106202595A (en) * 2015-05-07 2016-12-07 南京理工大学 Analyze the time domain non-conformal grid method of dielectric object Transient Em Fields Scattering characteristic
CN106547950A (en) * 2016-10-11 2017-03-29 西安电子科技大学 A kind of method for reducing geometric optics and physical optics hybrid algorithm error
CN106599331A (en) * 2015-10-16 2017-04-26 南京理工大学 Antenna optimization method capable of combining moment method with confidence region space mapping algorithm
CN106844919A (en) * 2017-01-11 2017-06-13 三峡大学 A kind of method for solving of transmission line of electricity scattered field
CN107292014A (en) * 2017-06-15 2017-10-24 北京航空航天大学 A kind of method of the qualitative assessment radar station guidance capability based on Monte Carlo simulation
CN107391880A (en) * 2017-08-18 2017-11-24 上海无线电设备研究所 A kind of accurate efficient radome electrical property energy analysis method
CN107562981A (en) * 2017-07-14 2018-01-09 西安电子科技大学 Electric Field Calculation method and device in conductor scattering from rough surface problem
CN107967377A (en) * 2017-11-02 2018-04-27 上海无线电设备研究所 Inhomogeneous plasma electromagnetic property emulation mode based on Monte Carlo method
CN108614251A (en) * 2016-12-11 2018-10-02 南京理工大学 The radar cross section Visual calculation method of computing of multi-coated targets with RAM
CN109100692A (en) * 2018-06-19 2018-12-28 西安电子科技大学 Rough surface and multiple target combination scattering emulation modes based on iterative physical optics
CN109188384A (en) * 2018-08-31 2019-01-11 西安电子科技大学 The Electromagnetic Simulation method of extraterrestrial target echo dynamic observation
CN109215072A (en) * 2018-07-24 2019-01-15 杭州电子科技大学 Chaff cloud RCS acquisition methods based on tanh-sinh integration method
CN109212498A (en) * 2017-07-03 2019-01-15 上海东峻信息科技有限公司 A kind of fast algorithm of air formation radar cross section
CN109783829A (en) * 2017-11-13 2019-05-21 西安电子科技大学 A kind of electromagnetic field prediction technique of three-dimensional FEM mixing two dimension FMM
CN110309524A (en) * 2019-03-21 2019-10-08 杭州电子科技大学 A kind of method of magnetostatic coupling field intensity between calculating magnetic-particle
CN110377974A (en) * 2019-06-28 2019-10-25 南京理工大学 Thunderstorm activity composite plate members structural optimization method
CN110502782A (en) * 2019-07-04 2019-11-26 西安电子科技大学 The big scale of construction chaff cloud electromagnetic scattering measuring method transported based on outline vector
CN111123225A (en) * 2019-12-25 2020-05-08 西安电子科技大学 Sea background foil strip cloud scattering method based on vector transport theory
CN111274703A (en) * 2020-01-20 2020-06-12 中国人民解放军国防科技大学 Electromagnetic characteristic analysis method and system for large-size target with complex structure
CN111368398A (en) * 2020-02-06 2020-07-03 北京环境特性研究所 Method and device for analyzing electromagnetic scattering characteristics of large electric target with uncertain structure
CN111832183A (en) * 2020-07-22 2020-10-27 西安电子科技大学 Electromagnetic simulation method based on octree structure MLFMA region decomposition
CN112329308A (en) * 2020-11-09 2021-02-05 西安电子科技大学 Ship model electromagnetic scattering simulation method based on iterative multilayer fast multipole
CN112507553A (en) * 2020-12-08 2021-03-16 中国人民解放军空军工程大学 Complex regional environment modeling method
CN113376597A (en) * 2021-06-18 2021-09-10 西安电子科技大学 Complex terrain electromagnetic scattering rapid simulation method based on digital elevation map and GPU
CN114741880A (en) * 2022-04-15 2022-07-12 北京理工大学 Method and device for calculating electromagnetic characteristics of group targets

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103279600A (en) * 2013-05-15 2013-09-04 西安电子科技大学 Medium rough surface finite element electromagnetic simulation method based on integral boundary

Patent Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103279600A (en) * 2013-05-15 2013-09-04 西安电子科技大学 Medium rough surface finite element electromagnetic simulation method based on integral boundary

Non-Patent Citations (8)

* Cited by examiner, † Cited by third party
Title
J.M.TABOADA 等: "IMPROVEMENT OF THE HYBRID", 《MICROWAVE AND OPTICAL TECHNOLOGY LETTERS》 *
常伟 等: "《2013全国天线年会论文集(下册)》", 30 November 2013 *
李杰 等: "复合电磁散射分析中的精确几何建模", 《西安电子科技大学学报(自然科学版)》 *
王蕊: "粗糙面及其与目标复合电磁散射中的相关问题研究", 《中国博士学位论文全文数据库 信息科技辑(月刊)》 *
秦三团: "随机粗糙面与目标复合电磁散射的快速算法研究", 《中国博士学位论文全文数据库 信息科技辑(月刊)》 *
聂小春 等: "计算开槽电大目标电磁散射的IPO -MOM混合法", 《西安电子科技大学学报》 *
郭立新 等: "MPI并行矩量法计算二维粗糙面波束电磁散射", 《***工程与电子技术》 *
郭立新 等: "二维带限分形粗糙面电磁散射的基尔霍夫近似", 《电子学报》 *

Cited By (43)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104317984B (en) * 2014-09-09 2017-06-27 中国舰船研究设计中心 Forecasting Methodology and system are scattered based on a point ship electromagnetic for domain modeling
CN104317984A (en) * 2014-09-09 2015-01-28 中国舰船研究设计中心 Ship electromagnetic scattering prediction method and system based on sub-domain modeling
CN106156394A (en) * 2015-04-22 2016-11-23 南京理工大学 Electromagnetic property extracting method based on explicit difference scheme
CN106202595A (en) * 2015-05-07 2016-12-07 南京理工大学 Analyze the time domain non-conformal grid method of dielectric object Transient Em Fields Scattering characteristic
CN104992035B (en) * 2015-07-28 2017-03-15 中国人民解放军国防科学技术大学 A kind of quick calculation method of Terahertz frequency range rough surface Electromagnetic Scattering of Target
CN104992035A (en) * 2015-07-28 2015-10-21 中国人民解放军国防科学技术大学 Quick calculation method for terahertz band surface rough target electromagnetic scattering
CN106599331B (en) * 2015-10-16 2019-11-29 南京理工大学 The antenna optimization method of moment method combination confidence region space mapping algorithm
CN106599331A (en) * 2015-10-16 2017-04-26 南京理工大学 Antenna optimization method capable of combining moment method with confidence region space mapping algorithm
CN106547950A (en) * 2016-10-11 2017-03-29 西安电子科技大学 A kind of method for reducing geometric optics and physical optics hybrid algorithm error
CN106547950B (en) * 2016-10-11 2019-11-01 西安电子科技大学 A method of reducing geometric optics and physical optics hybrid algorithm error
CN108614251A (en) * 2016-12-11 2018-10-02 南京理工大学 The radar cross section Visual calculation method of computing of multi-coated targets with RAM
CN106844919A (en) * 2017-01-11 2017-06-13 三峡大学 A kind of method for solving of transmission line of electricity scattered field
CN106844919B (en) * 2017-01-11 2020-09-01 三峡大学 Solution method for scattered field of power transmission line
CN107292014A (en) * 2017-06-15 2017-10-24 北京航空航天大学 A kind of method of the qualitative assessment radar station guidance capability based on Monte Carlo simulation
CN107292014B (en) * 2017-06-15 2018-03-30 北京航空航天大学 A kind of method of the qualitative assessment radar station guidance capability based on Monte Carlo simulation
CN109212498A (en) * 2017-07-03 2019-01-15 上海东峻信息科技有限公司 A kind of fast algorithm of air formation radar cross section
CN107562981B (en) * 2017-07-14 2020-01-10 西安电子科技大学 Electric field calculation method and device in conductor rough surface scattering problem
CN107562981A (en) * 2017-07-14 2018-01-09 西安电子科技大学 Electric Field Calculation method and device in conductor scattering from rough surface problem
CN107391880A (en) * 2017-08-18 2017-11-24 上海无线电设备研究所 A kind of accurate efficient radome electrical property energy analysis method
CN107391880B (en) * 2017-08-18 2020-07-31 上海无线电设备研究所 Accurate and efficient antenna housing electrical property analysis method
CN107967377B (en) * 2017-11-02 2021-01-19 上海无线电设备研究所 Monte Carlo method-based non-uniform plasma electromagnetic characteristic simulation method
CN107967377A (en) * 2017-11-02 2018-04-27 上海无线电设备研究所 Inhomogeneous plasma electromagnetic property emulation mode based on Monte Carlo method
CN109783829A (en) * 2017-11-13 2019-05-21 西安电子科技大学 A kind of electromagnetic field prediction technique of three-dimensional FEM mixing two dimension FMM
CN109783829B (en) * 2017-11-13 2023-02-17 西安电子科技大学 Electromagnetic field prediction method for three-dimensional FEM (field emission modeling) and two-dimensional FMM (frequency modulation modeling)
CN109100692A (en) * 2018-06-19 2018-12-28 西安电子科技大学 Rough surface and multiple target combination scattering emulation modes based on iterative physical optics
CN109215072A (en) * 2018-07-24 2019-01-15 杭州电子科技大学 Chaff cloud RCS acquisition methods based on tanh-sinh integration method
CN109215072B (en) * 2018-07-24 2021-07-20 杭州电子科技大学 Foil strip cloud RCS acquisition method
CN109188384B (en) * 2018-08-31 2023-03-21 西安电子科技大学 Electromagnetic simulation method for space target echo dynamic observation
CN109188384A (en) * 2018-08-31 2019-01-11 西安电子科技大学 The Electromagnetic Simulation method of extraterrestrial target echo dynamic observation
CN110309524B (en) * 2019-03-21 2022-11-11 杭州电子科技大学 Method for calculating strength of magnetostatic interaction field between magnetic particles
CN110309524A (en) * 2019-03-21 2019-10-08 杭州电子科技大学 A kind of method of magnetostatic coupling field intensity between calculating magnetic-particle
CN110377974A (en) * 2019-06-28 2019-10-25 南京理工大学 Thunderstorm activity composite plate members structural optimization method
CN110502782A (en) * 2019-07-04 2019-11-26 西安电子科技大学 The big scale of construction chaff cloud electromagnetic scattering measuring method transported based on outline vector
CN111123225B (en) * 2019-12-25 2022-01-04 西安电子科技大学 Sea background foil strip cloud scattering method based on vector transport theory
CN111123225A (en) * 2019-12-25 2020-05-08 西安电子科技大学 Sea background foil strip cloud scattering method based on vector transport theory
CN111274703A (en) * 2020-01-20 2020-06-12 中国人民解放军国防科技大学 Electromagnetic characteristic analysis method and system for large-size target with complex structure
CN111368398A (en) * 2020-02-06 2020-07-03 北京环境特性研究所 Method and device for analyzing electromagnetic scattering characteristics of large electric target with uncertain structure
CN111368398B (en) * 2020-02-06 2023-12-29 北京环境特性研究所 Electromagnetic scattering characteristic analysis method and device for electrically large target with uncertain structure
CN111832183A (en) * 2020-07-22 2020-10-27 西安电子科技大学 Electromagnetic simulation method based on octree structure MLFMA region decomposition
CN112329308A (en) * 2020-11-09 2021-02-05 西安电子科技大学 Ship model electromagnetic scattering simulation method based on iterative multilayer fast multipole
CN112507553A (en) * 2020-12-08 2021-03-16 中国人民解放军空军工程大学 Complex regional environment modeling method
CN113376597A (en) * 2021-06-18 2021-09-10 西安电子科技大学 Complex terrain electromagnetic scattering rapid simulation method based on digital elevation map and GPU
CN114741880A (en) * 2022-04-15 2022-07-12 北京理工大学 Method and device for calculating electromagnetic characteristics of group targets

Similar Documents

Publication Publication Date Title
CN103870654A (en) Electromagnetic scattering simulation method based on parallel moment method and physical optics mixing
CN103593510B (en) Matsurface and target composite electromagnetic scattering emulation mode based on reciprocity principle
CN102401898A (en) Quantified simulation method for forest remote sensing data of synthetic aperture radar
CN103777186B (en) Decompose and moving target near field radar return characteristic computing method under local irradiation based on wave beam
CN103400004A (en) Medium rough surface electromagnetic scattering simulation method based on multi-zone model moment method
CN105445711A (en) Sea level essential factor SAR original data simulation method based on inverse Omega-K algorithm
CN103279600A (en) Medium rough surface finite element electromagnetic simulation method based on integral boundary
CN103699810A (en) Modeling method of rough surface microwave band bidirectional reflectance distribution function
CN104318593A (en) Simulation method and system of radar sea clusters
CN107271977A (en) High-precision SAR echo simulation methods based on mobile driving source FDTD algorithms
CN103425816A (en) Matrix extraction method for rapidly acquiring electromagnetic scattering characteristics of metal bodies of revolution
CN103954969B (en) Based on the dynamic conical target attitude angle of Laser One-dimensional Range Profile inverting and physical dimension
CN104101869B (en) A kind of ground wave radar moving target simulation modeling method under polar coordinate
CN106649900B (en) Time domain analysis method for electromagnetic characteristics of non-uniform rotation symmetric body
CN104915326A (en) Domain decomposition order stepping time domain integration method based on equivalence principle
CN106326579A (en) Method for generating visualizing of land and sea environment based on electromagnetic scattering technique
CN103675777B (en) Based on airborne radar clutter analogy method and the device of fitting process
Tang et al. Simulation of the Doppler echoes from wind turbine based on scattered electric field calculation
CN104346488A (en) Hybrid modeling and electromagnetic scattering high-speed simulation method of electrically large metal object with complicated shape
CN105760343A (en) Time-domain integration method for analyzing electromagnetic scattering characteristics of multiple complicated metal targets
CN105095546A (en) Mixed-order Nystrom method for analyzing electromagnetic scattering characteristics of multi-scale conductive object
CN106156431A (en) Target conductor electromagnetic scattering emulation mode based on nurbs surface modeling
He et al. Analysis of the wind turbine RCS and micro-Doppler feature based on FEKO
CN114966870B (en) Multi-component joint detection transient electromagnetic method for arbitrary position of ground loop
CN104699870B (en) TV university complexity Lossy Dielectric Object electromagnetic scattering parabola rapid simulation method

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
WD01 Invention patent application deemed withdrawn after publication
WD01 Invention patent application deemed withdrawn after publication

Application publication date: 20140618