CN105119683A - Unmanned plane communication interference countermeasure method based on real-time embedded control system - Google Patents

Abstract

The invention relates to an unmanned plane communication interference countermeasure method based on a real-time embedded control system and belongs to the technical field of unmanned plane communication countermeasure. The method comprises steps of: S1, constructing an unmanned plane antenna model; S2: analyzing the characteristic of the antenna by using a time-domain finite difference method; S3, acquiring the field distribution of the antenna fixed in a coordinate system on a spherical surface at a certain distance in a far area; S4, transforming an acquired value into field distribution of the antenna in any attitude by using transformation of coordinates; S5, transforming the field distribution on the spherical surface into ground in order to obtain a relation between ground coordinates and effective interference power; and S6, adjusting the flying attitude of an unmanned plane in order to perform ground communication interference. The method may effectively and accurately adjust the flying attitude of the unmanned plane in order to maximize the interference effect of the unmanned plane on a ground communication system and effectively prevent defects of conventional ground communication countermeasure, and has high application value in the field of unmanned plane communication countermeasure.

Description

A kind of interference of the UAV Communication based on real-time embedded control system countercheck

Technical field

The present invention relates to a kind of UAV Communication based on real-time embedded control system interference countercheck, belong to UAV Communication countermeasure techniques field.

Background technology

Communication countermeasures is the important component part of electronic countermeasures.Radio interference of communications is the jamming countermeasure taked for the other side's radio communication.Due to the medium of radio communication--electromagnetic wave has a characteristic: when the electromagnetic wave co-operation of two or more same frequency, and receiving equipment will receive the superposition of this two or more signal, thus makes Received signal strength smudgy.Utilize this characteristic electromagnetic, launch the interference signal identical with the other side's electrical signal frequency by jamming device, the radio receiving equipment of the other side can be made malfunctioning.Communication Jamming generally includes to be covered true information and manufactures false information two aspects.Communication Jamming uses in communication countermeasures field the most extensively, acts on maximum electronic attack means.Implement effective communication countermeasures, the communication of enemy, command ability can be reduced.

Traditional communication countermeasures is mainly conceived to the electronic communication antagonism on ground, its limitation is more obvious in modern communications war, there is the contradictory problems that interference communicates between ourselves and the enemy in tradition ground communication countermeasures technology, and ground communication interference always, and the antenna of employing will regulate as required.Directional antenna is used when orientation is scouted.In short-wave band, antenna size is comparatively large, and rhombus and logarithm periodic antenna use comparatively general in fixed station; But when mobile use, be just difficult to realize, can only based on whip antenna.In ultrashort wave, microwave band, parabolic reflector antenna and logarithm periodic antenna are widely used.Use directional antenna, can increase the gain of antenna, improve reception, but have limitation on direction, sometimes can only use the omnidirectional antenna of low gain, therefore its antenna adjustments complex procedures, cannot reach the heightening gain of unmanned plane high-altitude radio wave propagation.

Summary of the invention

The invention provides a kind of UAV Communication based on real-time embedded control system interference countercheck, for the traditional ground communication countermeasures technology interference of solution between ourselves and the enemy, antenna adjustments complex procedures, and the problem that cannot reach high-altitude radio wave propagation gain.

The concrete steps of the interference of the UAV Communication based on real-time embedded control system countercheck of the present invention are as follows:

Step1, structure unmanned plane antenna model: utilize the large scale of unmanned plane flank itself to form conformal array antenna and be used as unmanned plane antenna model; Because the VNF potato masher antenna of routine is installed more difficult on unmanned plane body, therefore adopt conformal array antenna, namely utilize the large scale of unmanned plane flank itself to form conformal array antenna, effective Enhanced Radiation Reduced Blast area, improve antenna efficiency;

Step2, utilize FDTD analysis antenna performance:

Step3, ask for and be fixed on the Electric Field Distribution on a certain distance sphere of far field of antenna in coordinate system: for the earth axes S in a certain moment _gin be in the unmanned plane antenna of any attitude, it is at body axis system S _bin position be fixing, therefore in body axis system, calculate the Electric Field Distribution on its far field sphere under a certain operating frequency with Finite-Difference Time-Domain Method;

Step4, recycling coordinate transform, be converted to Electric Field Distribution when antenna is in any attitude above: wherein, prerequisite defines earth axes Ox by the value of trying to achieve _{g '}y _{g '}z _g(S _{g '}) and body axis system O _bx _by _bz _b(S _b), relation then between two coordinate systems is determined by unmanned plane three flight attitudes: yaw angle α, angle of pitch β and roll angle γ, its direction is all just to transfer to along corresponding coordinate dextrorotation, and the above primitive spin matrix corresponding to three attitude transformation of note is respectively L _gbz(α), L _gby(β) and L _gbx(γ), then from body axis system S _bto earth axes S _gtransformation matrix of coordinates Lgb (α, beta, gamma) be:

L _gb(α,β,γ)＝L _gbz(α)L _gby(β)L _gbx(γ)

According to coordinate transform formula above, the Electric Field Distribution of unmanned plane antenna when any attitude can be tried to achieve;

Step5, finally the Electric Field Distribution on sphere is transformed into ground, draws the corresponding relation of geographical coordinates and effective interference power: at earth axes S _gin, the coordinate of N point is (x _g, y _g, z _g), the intersection point N on note ON and ground ₀coordinate is (x _g0, y _g0, z _g0), then there is z _g0=-h, h are unmanned plane during flying height, obtain according to triangle geometric ratio relation:

$x_{g0}=\frac{-h}{z_g}{x}_g,y_{g0}=\frac{-h}{z_g}{y}_g$

Because the radio wave propagation mode between air-ground is considered as free space, far field electric field and propagation distance are inversely proportional to, and according to this relation, the electric field of ground N point are transformed into N ₀point, its component array is (E ₀) _ghave:

${\left(E_0\right)}_g=\frac{ON}{{\mathrm{ON}}_0}{\left(E\right)}_g=\frac{z_g}{-h}{\left(E\right)}_g$

By above-mentioned conversion, try to achieve earth axes S _gmiddle z is the corresponding relation of xg0 and yg0 and electric field E0 on ground, namely obtains ground field strength distribution;

Step6, adjustment unmanned plane during flying attitude carry out ground communication interference: according to the directivity of terrestrial receiving antenna, in conjunction with the ground field strength distribution calculated above, for regulating unmanned plane during flying attitude, carry out effective, maximized interference enemy communication system.

In described step Step2, the concrete steps of FDTD analysis antenna performance are utilized to be:

Step2.1, utilize conformal grid to sky line boundary Fdtd Method grid not intersection revise;

Step2.2, utilize thin wire Finite-Difference Time-Domain Method to loading thin thread part be divided into row relax;

Step2.3, driving source arrange the cutler feed method adopting additional incentive;

Step2.4, in finite region, carry out recurrence calculation, adopt Berenger completely permutation absorbing boundary condition;

Step2.5, calculating radiation pattern.

The invention has the beneficial effects as follows: the present invention adopts Finite-Difference Time-Domain Method, in conjunction with antenna asymmetry, analyze for the unmanned plane antenna model built, process and calculating, show that antenna greatest irradiation direction is long-armed one end, Electric Field Distribution when antenna being in any attitude according to coordinate transform is again transformed into ground, thus unmanned plane during flying attitude can be regulated effectively accurately, to realize the maximization of unmanned plane to Ground Communication System interference effect, effectively avoid the drawback of traditional ground communication countermeasures, higher using value is had in UAV Communication antagonism field.

Accompanying drawing explanation

Fig. 1 is flow chart of the present invention;

Fig. 2 is unmanned plane antenna model figure of the present invention;

Fig. 3 is unmanned plane coordinate system figure of the present invention.

Embodiment

Embodiment 1: as Figure 1-3, a kind of interference of the UAV Communication based on real-time embedded control system countercheck, the concrete steps of the described interference of the UAV Communication based on real-time embedded control system countercheck are as follows:

Step1, structure unmanned plane antenna model: utilize the large scale of unmanned plane flank itself to form conformal array antenna and be used as unmanned plane antenna model; Because the VNF potato masher antenna of routine is installed more difficult on unmanned plane body, therefore adopt conformal array antenna, namely utilize the large scale of unmanned plane flank itself to form conformal array antenna, effective Enhanced Radiation Reduced Blast area, improve antenna efficiency;

Step2, utilize FDTD analysis antenna performance:

Step3, ask for and be fixed on the Electric Field Distribution on a certain distance sphere of far field of antenna in coordinate system: for the earth axes S in a certain moment _gin be in the unmanned plane antenna of any attitude, it is at body axis system S _bin position be fixing, therefore in body axis system, calculate the Electric Field Distribution on its far field sphere under a certain operating frequency with Finite-Difference Time-Domain Method;

Step4, recycling coordinate transform, be converted to Electric Field Distribution when antenna is in any attitude above: wherein, prerequisite defines earth axes Ox by the value of trying to achieve _{g '}y _{g '}z _g(S _{g '}) and body axis system O _bx _by _bz _b(S _b), relation then between two coordinate systems is determined by unmanned plane three flight attitudes: yaw angle α, angle of pitch β and roll angle γ, its direction is all just to transfer to along corresponding coordinate dextrorotation, and the above primitive spin matrix corresponding to three attitude transformation of note is respectively L _gbz(α), L _gby(β) and L _gbx(γ), then from body axis system S _bto earth axes S _gtransformation matrix of coordinates Lgb (α, beta, gamma) be:

L _gb(α,β,γ)＝L _gbz(α)L _gby(β)L _gbx(γ)

According to coordinate transform formula above, the Electric Field Distribution of unmanned plane antenna when any attitude can be tried to achieve;

Step5, finally the Electric Field Distribution on sphere is transformed into ground, draws the corresponding relation of geographical coordinates and effective interference power: at earth axes S _gin, the coordinate of N point is (x _g, y _g, z _g), the intersection point N on note ON and ground ₀coordinate is (x _g0, y _g0, z _g0), then there is z _g0=-h, h are unmanned plane during flying height, obtain according to triangle geometric ratio relation:

$x_{g0}=\frac{-h}{z_g}{x}_g,y_{g0}=\frac{-h}{z_g}{y}_g$

Because the radio wave propagation mode between air-ground is considered as free space, far field electric field and propagation distance are inversely proportional to, and according to this relation, the electric field of ground N point are transformed into N ₀point, its component array is (E ₀) _ghave:

${\left(E_0\right)}_g=\frac{ON}{{\mathrm{ON}}_0}{\left(E\right)}_g=\frac{z_g}{-h}{\left(E\right)}_g$

By above-mentioned conversion, try to achieve earth axes S _gmiddle z is the corresponding relation of xg0 and yg0 and electric field E0 on ground, namely obtains ground field strength distribution;

Step6, adjustment unmanned plane during flying attitude carry out ground communication interference: according to the directivity of terrestrial receiving antenna, in conjunction with the ground field strength distribution calculated above, for regulating unmanned plane during flying attitude, carry out effective, maximized interference enemy communication system.

In described step Step2, the concrete steps of FDTD analysis antenna performance are utilized to be:

Step2.1, utilize conformal grid to sky line boundary Fdtd Method grid not intersection revise;

Step2.2, utilize thin wire Finite-Difference Time-Domain Method to loading thin thread part be divided into row relax;

Step2.3, driving source arrange the cutler feed method adopting additional incentive;

Step2.4, in finite region, carry out recurrence calculation, adopt Berenger completely permutation absorbing boundary condition;

Step2.5, calculating radiation pattern.

Embodiment 2: as Figure 1-3, a kind of interference of the UAV Communication based on real-time embedded control system countercheck, the concrete steps of the described interference of the UAV Communication based on real-time embedded control system countercheck are as follows:

Step1, structure unmanned plane antenna model: utilize the large scale of unmanned plane flank itself to form conformal array antenna and be used as unmanned plane antenna model; Because the VNF potato masher antenna of routine is installed more difficult on unmanned plane body, therefore adopt conformal array antenna, namely utilize the large scale of unmanned plane flank itself to form conformal array antenna, effective Enhanced Radiation Reduced Blast area, improve antenna efficiency;

Step2, utilize FDTD analysis antenna performance:

Step3, ask for and be fixed on the Electric Field Distribution on a certain distance sphere of far field of antenna in coordinate system: for the earth axes S in a certain moment _gin be in the unmanned plane antenna of any attitude, it is at body axis system S _bin position be fixing, therefore in body axis system, calculate the Electric Field Distribution on its far field sphere under a certain operating frequency with Finite-Difference Time-Domain Method;

Step4, recycling coordinate transform, be converted to Electric Field Distribution when antenna is in any attitude above: wherein, prerequisite defines earth axes Ox by the value of trying to achieve _{g '}y _{g '}z _g(S _{g '}) and body axis system O _bx _by _bz _b(S _b), relation then between two coordinate systems is determined by unmanned plane three flight attitudes: yaw angle α, angle of pitch β and roll angle γ, its direction is all just to transfer to along corresponding coordinate dextrorotation, and the above primitive spin matrix corresponding to three attitude transformation of note is respectively L _gbz(α), L _gby(β) and L _gbx(γ), then from body axis system S _bto earth axes S _gtransformation matrix of coordinates Lgb (α, beta, gamma) be:

L _gb(α,β,γ)＝L _gbz(α)L _gby(β)L _gbx(γ)

According to coordinate transform formula above, the Electric Field Distribution of unmanned plane antenna when any attitude can be tried to achieve;

Step5, finally the Electric Field Distribution on sphere is transformed into ground, draws the corresponding relation of geographical coordinates and effective interference power: at earth axes S _gin, the coordinate of N point is (x _g, y _g, z _g), the intersection point N on note ON and ground ₀coordinate is (x _g0, y _g0, z _g0), then there is z _g0=-h, h are unmanned plane during flying height, obtain according to triangle geometric ratio relation:

$x_{g0}=\frac{-h}{z_g}{x}_g,y_{g0}=\frac{-h}{z_g}{y}_g$

Because the radio wave propagation mode between air-ground is considered as free space, far field electric field and propagation distance are inversely proportional to, and according to this relation, the electric field of ground N point are transformed into N ₀point, its component array is (E ₀) _ghave:

${\left(E_0\right)}_g=\frac{ON}{{\mathrm{ON}}_0}{\left(E\right)}_g=\frac{z_g}{-h}{\left(E\right)}_g$

By above-mentioned conversion, try to achieve earth axes S _gmiddle z is the corresponding relation of xg0 and yg0 and electric field E0 on ground, namely obtains ground field strength distribution;

Step6, adjustment unmanned plane during flying attitude carry out ground communication interference: according to the directivity of terrestrial receiving antenna, in conjunction with the ground field strength distribution calculated above, for regulating unmanned plane during flying attitude, carry out effective, maximized interference enemy communication system.

In described step Step2, the concrete steps of FDTD analysis antenna performance are utilized to be:

Step2.1, utilize conformal grid to sky line boundary Fdtd Method grid not intersection revise;

Step2.2, utilize thin wire Finite-Difference Time-Domain Method to loading thin thread part be divided into row relax;

Step2.3, driving source arrange the cutler feed method adopting additional incentive;

Step2.4, in finite region, carry out recurrence calculation, adopt Berenger completely permutation absorbing boundary condition;

Step2.5, calculating radiation pattern.

The concrete implementation step of described method is:

Step s1: build unmanned plane antenna model.

Because the VNF potato masher antenna of routine is installed more difficult on unmanned plane body, therefore adopt conformal array antenna, namely utilize the large scale of unmanned plane flank itself to form conformal array antenna, effective Enhanced Radiation Reduced Blast area, improve antenna efficiency.As shown in Figure 2, it is the asymmetric vibrator of a plane to antenna model figure, and operating frequency is 20-100MHz.Unmanned plane flank is made up of wing centre section, middle outer wing and outer wing, forms the two poles of the earth of antenna at side wing centre section and middle outer wing lower surface according to certain shape copper foil.

Step s2: utilize FDTD analysis antenna performance.

1, utilize conformal grid to sky line boundary Fdtd Method grid not intersection revise, consider during calculating that wing metal framework and dielectric material are on the impact of antenna performance, thus reach the real operational environment of antenna;

2, thin wire Finite-Difference Time-Domain Method is utilized to be divided into row relax to loading thin thread part;

3, driving source arranges the cutler feed method adopting additional incentive.In order to the broadband character by once calculating antenna, pumping signal adopts differential Gaussian pulse, and its computing formula is:

$f\left(n\mathrm{\Δ}t\right)=\frac{n\mathrm{\Δ}t-t_0}{\mathrm{\τ}}\exp \[-\frac{4\mathrm{\π}{(n\mathrm{\Δ}t-t_0)}^2}{{\mathrm{\τ}}^2}\]$

Wherein n △ t is pulse duration; π is circumference ratio value, and t0, τ are constant, determines zero point and the width of pulse respectively.

4, due to the restriction of calculated capacity, in finite region, recurrence calculation is carried out; In order to the aerial radiation process of open domain can be simulated, need to arrange absorbing boundary condition at the cutoff boundary place of zoning.Namely Berenger completely permutation absorbing boundary condition is adopted;

5, radiation pattern is calculated, in order to save amount of calculation, first the output boundary transient electric field obtained under pulse excitation and magnetic-field component are carried out Fourier conversion and obtain frequency domain value, and then inquire into outside the venue in corresponding frequencies point enterprising line frequency territory distance and far zone field.

Step s3: ask for and be fixed on the Electric Field Distribution of antenna on a certain distance sphere of far field in coordinate system.

For the earth axes S in a certain moment _gin be in the unmanned plane antenna of any attitude, it is at body axis system S _bin position be fixing, therefore in body axis system, calculate the Electric Field Distribution on its far field sphere under a certain operating frequency with Finite-Difference Time-Domain Method, be wherein M is a bit on sphere, and θ is O _bm and y _bthe angle of axle, for O _bm is at z _bo _bx _bplane projection and z _bthe angle of axle.

Step s4: recycling coordinate transform, is converted to Electric Field Distribution when antenna is in any attitude above: precondition definition earth axes Ox by the value of trying to achieve _{g '}y _{g '}z _g(S _g) and body axis system O _bx _by _bz _b(S _b), then the relation between two coordinate systems can be determined by unmanned plane three flight attitudes: yaw angle α, angle of pitch β and roll angle γ, its direction is all just to transfer to along corresponding coordinate dextrorotation, and unmanned plane reference axis figure is as shown in Figure 3.The above primitive spin matrix corresponding to three attitude transformation of note is respectively L _gbz(α), L _gby(β) and L _gbx(γ), then from body axis system S _bto earth axes S _gtransformation matrix of coordinates L _gb(α, beta, gamma) is:

L _gb(α,β,γ)＝L _gbz(α)L _gby(β)L _gbx(γ)

According to coordinate transform formula above, the Electric Field Distribution of unmanned plane antenna when any attitude can be tried to achieve.

Step s5: finally the Electric Field Distribution on sphere is transformed into ground, draws the corresponding relation of geographical coordinates and effective interference power:

At earth axes S _gin, the coordinate of N point is (x _g, y _g, z _g), the intersection point N on note ON and ground ₀coordinate is (x _g0, y _g0, z _g0), then there is z _g0=-h (h is unmanned plane during flying height), obtains according to triangle geometric ratio relation:

$x_{g0}=\frac{-h}{z_g}{x}_g,y_{g0}=\frac{-h}{z_g}{y}_g$

Because the radio wave propagation mode between air-ground can be considered as free space, far field electric field and propagation distance are inversely proportional to, and the electric field of ground N point can be transformed into N according to this relation ₀point, its component array is (E ₀) _ghave:

${\left(E_0\right)}_g=\frac{ON}{{\mathrm{ON}}_0}{\left(E\right)}_g=\frac{z_g}{-h}{\left(E\right)}_g$

By above-mentioned conversion, try to achieve earth axes S _gmiddle z is x on ground _g0and y _g0with electric field E ₀corresponding relation, namely obtain ground field strength distribution.Due to the directivity of reception antenna, useful interference field strength is E ₀a certain component E ₀', so interference effective coverage is exactly E on ground ₀' be greater than the region of certain value.

Simulated experiment setting ground receiver is vertical polarized antenna, the conformal array antenna model that unmanned plane adopts Fig. 2 to design, its flying height is 2500m, antenna input power is 110W, operating frequency is 60MHz, unmanned plane with different attitude to immediately below ground receiver to carry out the data of Communication Jamming as shown in table 1:

Table 1 UAV Communication interference experiment result

Unmanned plane during flying attitude	Communication Jamming radius (Km)
		Put down and fly γ=-0 °	9
Left bank γ=-20 °	17
		Left bank γ=-30 °	15
Right bank γ=20 °	5
		Right bank γ=30 °	3

By above experiment and instance data analysis, can find out that its Communication Jamming radius is also different when unmanned plane during flying attitude is different, when Communication Jamming is carried out in left bank flight, the communication radius of interference is maximum, puts down and flies to take second place, and right bank flight is the poorest.

Wherein earth axes Ox _{g '}y _{g '}z _g(S _{g '}) be represent according to the letter representation in Fig. 3, wherein Ox _{g '}y _{g '}z _gthree-dimensional coordinate, (S _{g '}) representing earth axes code name, g represents ground, corresponding (S _b) be body axis system code name, b represents body; Component array is (E ₀) g, wherein E ₀electric field code name, (E ₀) g represents the component array of electric field intensity.

Earth axes Ox is defined in Fig. 3 _{g '}y _{g '}z _g(S _{g '}) and body axis system O _bx _by _bz _b(S _b), by Ox _{g '}, Oy _{g '}move to O _bx _g, O _by _b, thus obtain the earth axes O that involves mutually with body _bx _gy _gz _g(S _g).If O _bx _bat x _go _by _gplane be projected as O _bx _b', O _bz _bat x _bo _bz _gplane be projected as O _bz _b', then the relation between two coordinate systems can be determined by unmanned plane three flight attitudes: yaw angle α, angle of pitch β and roll angle γ, its direction all just to transfer to along corresponding coordinate dextrorotation, M ₀point is postulated point on the sphere of far field.

By reference to the accompanying drawings the specific embodiment of the present invention is explained in detail above, but the present invention is not limited to above-mentioned execution mode, in the ken that those of ordinary skill in the art possess, various change can also be made under the prerequisite not departing from present inventive concept.

Claims (2)

1. based on a UAV Communication interference countercheck for real-time embedded control system, it is characterized in that: the concrete steps of the described interference of the UAV Communication based on real-time embedded control system countercheck are as follows:

Step1, structure unmanned plane antenna model: utilize the large scale of unmanned plane flank itself to form conformal array antenna and be used as unmanned plane antenna model;

Step2, utilize FDTD analysis antenna performance:

Step3, ask for and be fixed on the Electric Field Distribution on a certain distance sphere of far field of antenna in coordinate system: for the earth axes S in a certain moment _gin be in the unmanned plane antenna of any attitude, it is at body axis system S _bin position be fixing, therefore in body axis system, calculate the Electric Field Distribution on its far field sphere under a certain operating frequency with Finite-Difference Time-Domain Method;

Step4, recycling coordinate transform, be converted to Electric Field Distribution when antenna is in any attitude above by the value of trying to achieve:

Wherein, prerequisite defines earth axes Ox _{g '}y _{g '}z _g(S _{g '}) and body axis system O _bx _by _bz _b(S _b), relation then between two coordinate systems is determined by unmanned plane three flight attitudes: yaw angle α, angle of pitch β and roll angle γ, its direction is all just to transfer to along corresponding coordinate dextrorotation, and the above primitive spin matrix corresponding to three attitude transformation of note is respectively L _gbz(α), L _gby(β) and L _gbx(γ), then from body axis system S _bto earth axes S _gtransformation matrix of coordinates Lgb (α, beta, gamma) be:

L _gb(α,β,γ)＝L _gbz(α)L _gby(β)L _gbx(γ)

According to coordinate transform formula above, the Electric Field Distribution of unmanned plane antenna when any attitude can be tried to achieve;

Step5, finally the Electric Field Distribution on sphere is transformed into ground, draws the corresponding relation of geographical coordinates and effective interference power: at earth axes S _gin, the coordinate of N point is (x _g, y _g, z _g), the intersection point N on note ON and ground ₀coordinate is (x _g0, y _g0, z _g0), then there is z _g0=-h, h are unmanned plane during flying height, obtain according to triangle geometric ratio relation:

$x_{g0}=\frac{-h}{z_g}{x}_g,y_{g0}=\frac{-h}{z_g}{y}_g$

Because the radio wave propagation mode between air-ground is considered as free space, far field electric field and propagation distance are inversely proportional to, and according to this relation, the electric field of ground N point are transformed into N ₀point, its component array is (E ₀) _ghave:

${\left(E_0\right)}_g=\frac{ON}{{\mathrm{ON}}_0}{\left(E\right)}_g=\frac{z_g}{-h}{\left(E\right)}_g$

By above-mentioned conversion, try to achieve earth axes S _gmiddle z is the corresponding relation of xg0 and yg0 and electric field E0 on ground, namely obtains ground field strength distribution;

Step6, adjustment unmanned plane during flying attitude carry out ground communication interference: according to the directivity of terrestrial receiving antenna, in conjunction with the ground field strength distribution calculated above, for regulating unmanned plane during flying attitude, carry out effective, maximized interference enemy communication system.

2. the UAV Communication based on real-time embedded control system according to claim 1 interference countercheck, is characterized in that: in described step Step2, utilizes the concrete steps of FDTD analysis antenna performance to be:

Step2.1, utilize conformal grid to sky line boundary Fdtd Method grid not intersection revise;

Step2.2, utilize thin wire Finite-Difference Time-Domain Method to loading thin thread part be divided into row relax;

Step2.3, driving source arrange the cutler feed method adopting additional incentive;

Step2.4, in finite region, carry out recurrence calculation, adopt Berenger completely permutation absorbing boundary condition;

Step2.5, calculating radiation pattern.