CN105119683B - A kind of UAV Communication interference countercheck based on real-time embedded control system - Google Patents

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

Info

Publication number
CN105119683B
CN105119683B CN201510508339.7A CN201510508339A CN105119683B CN 105119683 B CN105119683 B CN 105119683B CN 201510508339 A CN201510508339 A CN 201510508339A CN 105119683 B CN105119683 B CN 105119683B
Authority
CN
China
Prior art keywords
msub
mrow
electric field
antenna
unmanned plane
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.)
Active
Application number
CN201510508339.7A
Other languages
Chinese (zh)
Other versions
CN105119683A (en
Inventor
张晶
薛冷
范洪博
容会
崔毅
李润鑫
汤守国
潘晟旻
孙俊
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Kunming University of Science and Technology
Original Assignee
Kunming University of Science and Technology
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 Kunming University of Science and Technology filed Critical Kunming University of Science and Technology
Priority to CN201510508339.7A priority Critical patent/CN105119683B/en
Publication of CN105119683A publication Critical patent/CN105119683A/en
Application granted granted Critical
Publication of CN105119683B publication Critical patent/CN105119683B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Variable-Direction Aerials And Aerial Arrays (AREA)
  • Control Of Position, Course, Altitude, Or Attitude Of Moving Bodies (AREA)

Abstract

The present invention relates to a kind of UAV Communication based on real-time embedded control system to disturb countercheck, belongs to UAV Communication countermeasure techniques field.The present invention includes step：S1：Build unmanned plane antenna model；S2：Utilize FDTD analysis antenna performance；S3：Antenna is asked for being fixed in coordinate system in a certain Electric Field Distribution on sphere in far field；S4：Coordinate transform is recycled, the value tried to achieve above is converted into Electric Field Distribution when antenna is in any attitude；S5：The Electric Field Distribution on sphere is finally transformed into ground, draws the corresponding relation of geographical coordinates and effective jamming power；S6：Adjust unmanned plane during flying posture and carry out ground communication interference.The present invention effectively can accurately adjust unmanned plane during flying posture, to realize maximization of the unmanned plane to GCS Ground Communication System interference effect, the drawbacks of effectively avoiding traditional ground communication countermeasure, have higher application value in UAV Communication confrontation field.

Description

A kind of UAV Communication interference countercheck based on real-time embedded control system
Technical field
The present invention relates to a kind of UAV Communication based on real-time embedded control system to disturb countercheck, belongs to nobody Machine communication countermeasure technical field.
Background technology
Communication countermeasure is the important component of electronic countermeasure.Radio interference of communications be for other side's radio communication and The jamming countermeasure taken.Due to the medium of radio communication -- electromagnetic wave has a characteristic：When two or more identical frequencies Electromagnetic wave cooperation when, receiving device will receive the superposition of the two or multiple signals, so that reception signal is fuzzy not Clearly.Using this characteristic of electromagnetic wave, believed by jamming device transmitting and the interference of other side's electrical signal frequency identical Number, the radio receiving equipment failure of other side can be made.Communication Jamming generally includes to cover true information and manufacture false information two Aspect.Communication Jamming is using the electronic attack means most extensive, effect is maximum in communication countermeasure field.Implement effective communication Confrontation, communication, the command ability of enemy can be reduced.
Traditional communication countermeasure is mainly focused on the electronic communication confrontation on ground, and its limitation is cured in modern communicationses war Invention is aobvious, the contradictory problems that traditional ground communication countermeasure technology communicates between ourselves and the enemy in the presence of interference always, and ground communication is done Disturb, the antenna of use will be adjusted on demand.Directional aerial is used when orienting and scouting.In short-wave band, antenna size is larger, rhombus With log-periodic antenna in fixed station using more universal；But, can only be with whip antenna in movement in use, being difficult to realize Based on.In ultrashort wave, microwave band, parabolic reflector antenna and log-periodic antenna are widely used., can using directional aerial To increase the gain of antenna, improve reception, but have limitation on direction, the omnidirectional antennas of low gain can only be used sometimes Line, therefore its antenna adjustments complex procedures, it is unable to reach the heightening gain of unmanned plane high-altitude radio wave propagation.
The content of the invention
The invention provides a kind of UAV Communication based on real-time embedded control system to disturb countercheck, for Solve traditional ground communication countermeasure technology interference between ourselves and the enemy, antenna adjustments complex procedures, and be unable to reach high-altitude electric wave and pass The problem of broadcasting gain.
The specific steps of UAV Communication interference countercheck of the present invention based on real-time embedded control system are such as Under：
Step1, structure unmanned plane antenna model：Conformal array antenna is formed to make using the large scale of unmanned plane flank in itself For unmanned plane antenna model；Due to the VNF potato masher antennas of routine installed on unmanned plane body it is relatively difficult, therefore using conformal Array antenna, i.e., conformal array antenna is formed using the large scale of unmanned plane flank in itself, effective Enhanced Radiation Reduced Blast area, improve antenna Efficiency；
Step2, utilize FDTD analysis antenna performance：
Step3, antenna is asked for being fixed in coordinate system in a certain Electric Field Distribution on sphere in far field：For certain for the moment The earth axes S mutually involved with body carvedgIn be in any attitude unmanned plane antenna, it is in body axis system SbIn Position is fixed, therefore it is remote under a certain working frequency in body axis system with Finite-Difference Time-Domain Method it to be calculated Electric Field Distribution on area's sphere；
Step4, coordinate transform is recycled, the value tried to achieve above is converted into the electric field point when antenna is in any attitude Cloth：
Wherein, on condition that defining earth axes Oxg’yg’zg(Sg’) and body axis system Obxbybzb(Sb), by OXg′、 OYg' move to ObXg、ObYgIt is O to obtain with the earth axes that body mutually involvesbXgYgZg(Sg), then ObXgYgZg(Sg) with Obxbybzb(Sb) relation between two coordinate systems determines by three flight attitudes of unmanned plane：Yaw angle α, angle of pitch β and rolling Corner γ, its direction along corresponding coordinate dextrorotation all to switch to just, the primitive spin moment corresponding to note three above attitude transformation Battle array is respectively Lgbz(α)、Lgby(β) and Lgbx(γ), then from body axis system SbTo the earth axes S mutually involved with bodyg's Transformation matrix of coordinates Lgb (α, beta, gamma) is：
Lgb(α, β, γ)=Lgbz(α)Lgby(β)Lgbx(γ)
According to coordinate transform formula above, Electric Field Distribution of the unmanned plane antenna in any attitude can be tried to achieve；
Step5, the Electric Field Distribution on sphere is finally transformed into ground, draws pair of geographical coordinates and effective jamming power It should be related to：In earth axes SgIn, N points are a bit on the sphere electric field of far field, and the coordinate of N points is (xg,yg,zg), remember ObN with The intersection point N on ground0Coordinate is (xg0,yg0,zg0), then there is zg0=-h, h are unmanned plane during flying height, compare relation according to triangle etc. ：
Because the radio wave propagation mode between air-ground is considered as free-space propagation, far field electric field is inversely proportional with propagation distance, The electric field of N points is transformed into by N according to this relation0Point, its component array are (E0)gHave：
By above-mentioned conversion, the earth axes S mutually involved with body is tried to achievegMiddle z is x on groundg0And yg0With electric field E0 Corresponding relation, that is, obtained the field strength distribution on ground；
Wherein：Earth axes are OXg′Yg′Zg(Sg'), body axis system ObXbYbZb(Sb), by OXg′、OYg' translation To ObXg、ObYgIt is O to obtain with the earth axes that body mutually involvesbXgYgZg(Sg)；
E0It is that N points are mapped to N0The electric field of point；E is the electric field at origin O；(E)gFor electricity corresponding to the electric field E at origin O The component array of field vector；
Step6, regulation unmanned plane during flying posture carry out ground communication interference：According to the directionality of terrestrial receiving antenna, knot The field strength distribution on ground that the face of closing is calculated, for adjusting unmanned plane during flying posture, effectively, maximumlly disturb Enemy's communication system.
In the step Step2, concretely comprising the following steps for FDTD analysis antenna performance is utilized：
Step2.1, using conformal grid the day misaligned part of line boundary Fdtd Method grid is modified；
Step2.2, using thin wire Finite-Difference Time-Domain Method loading fine rule part is handled；
Step2.3, driving source set the cutler feed method using additional incentive；
Step2.4, recurrence calculation is carried out in finite region, using Berenger completely permutation absorbing boundary conditions；
Step2.5, calculate radiation pattern.
The beneficial effects of the invention are as follows：The present invention uses Finite-Difference Time-Domain Method, with reference to antenna asymmetry, for structure Unmanned plane antenna model analyzed, handled and calculated, it is long-armed one end to draw antenna greatest irradiation direction, further according to coordinate Become Electric Field Distribution when antenna of changing commanders is in any attitude and be transformed into ground, so as to effectively accurately adjust unmanned plane during flying Posture, to realize maximization of the unmanned plane to GCS Ground Communication System interference effect, effectively avoid traditional ground communication countermeasure The drawbacks of, there is higher application value in UAV Communication confrontation field.
Brief description of the drawings
Fig. 1 is the flow chart of the present invention；
Fig. 2 is the unmanned plane antenna model figure of the present invention；
Fig. 3 is the unmanned plane coordinate system figure of the present invention.
Embodiment
Embodiment 1：As Figure 1-3, a kind of UAV Communication interference confrontation side based on real-time embedded control system Method, the UAV Communication interference countercheck based on real-time embedded control system comprise the following steps that：
Step1, structure unmanned plane antenna model：Conformal array antenna is formed to make using the large scale of unmanned plane flank in itself For unmanned plane antenna model；Due to the VNF potato masher antennas of routine installed on unmanned plane body it is relatively difficult, therefore using conformal Array antenna, i.e., conformal array antenna is formed using the large scale of unmanned plane flank in itself, effective Enhanced Radiation Reduced Blast area, improve antenna Efficiency；
Step2, utilize FDTD analysis antenna performance：
Step3, antenna is asked for being fixed in coordinate system in a certain Electric Field Distribution on sphere in far field：For certain for the moment The earth axes S mutually involved with body carvedgIn be in any attitude unmanned plane antenna, it is in body axis system SbIn Position is fixed, therefore it is remote under a certain working frequency in body axis system with Finite-Difference Time-Domain Method it to be calculated Electric Field Distribution on area's sphere；
Step4, coordinate transform is recycled, the value tried to achieve above is converted into the electric field point when antenna is in any attitude Cloth：
Wherein, on condition that defining earth axes Oxg’yg’zg(Sg’) and body axis system Obxbybzb(Sb), by OXg′、 OYg' move to ObXg、ObYgIt is O to obtain with the earth axes that body mutually involvesbXgYgZg(Sg), then ObXgYgZg(Sg) with Obxbybzb(Sb) relation between two coordinate systems determines by three flight attitudes of unmanned plane：Yaw angle α, angle of pitch β and rolling Corner γ, its direction along corresponding coordinate dextrorotation all to switch to just, the primitive spin moment corresponding to note three above attitude transformation Battle array is respectively Lgbz(α)、Lgby(β) and Lgbx(γ), then from body axis system SbTo the earth axes S mutually involved with bodyg's Transformation matrix of coordinates Lgb (α, beta, gamma) is：
Lgb(α, β, γ)=Lgbz(α)Lgby(β)Lgbx(γ)
According to coordinate transform formula above, Electric Field Distribution of the unmanned plane antenna in any attitude can be tried to achieve；
Step5, the Electric Field Distribution on sphere is finally transformed into ground, draws pair of geographical coordinates and effective jamming power It should be related to：In earth axes SgIn, N points are a bit on the sphere electric field of far field, and the coordinate of N points is (xg,yg,zg), remember ObN with The intersection point N on ground0Coordinate is (xg0,yg0,zg0), then there is zg0=-h, h are unmanned plane during flying height, compare relation according to triangle etc. ：
Because the radio wave propagation mode between air-ground is considered as free-space propagation, far field electric field is inversely proportional with propagation distance, The electric field of N points is transformed into by N according to this relation0Point, its component array are (E0)gHave：
By above-mentioned conversion, the earth axes S mutually involved with body is tried to achievegMiddle z is x on groundg0And yg0With electric field E0 Corresponding relation, that is, obtained the field strength distribution on ground；
Wherein：Earth axes are OXg′Yg′Zg(Sg'), body axis system ObXbYbZb(Sb), by OXg′、OYg' translation To ObXg、ObYgIt is O to obtain with the earth axes that body mutually involvesbXgYgZg(Sg)；
E0It is that N points are mapped to N0The electric field of point；E is the electric field at origin O；(E)gFor electricity corresponding to the electric field E at origin O The component array of field vector；
Step6, regulation unmanned plane during flying posture carry out ground communication interference：According to the directionality of terrestrial receiving antenna, knot The field strength distribution on ground that the face of closing is calculated, for adjusting unmanned plane during flying posture, effectively, maximumlly disturb Enemy's communication system.
In the step Step2, concretely comprising the following steps for FDTD analysis antenna performance is utilized：
Step2.1, using conformal grid the day misaligned part of line boundary Fdtd Method grid is modified；
Step2.2, using thin wire Finite-Difference Time-Domain Method loading fine rule part is handled；
Step2.3, driving source set the cutler feed method using additional incentive；
Step2.4, recurrence calculation is carried out in finite region, using Berenger completely permutation absorbing boundary conditions；
Step2.5, calculate radiation pattern.
Embodiment 2：As Figure 1-3, a kind of UAV Communication interference confrontation side based on real-time embedded control system Method, the UAV Communication interference countercheck based on real-time embedded control system comprise the following steps that：
Step1, structure unmanned plane antenna model：Conformal array antenna is formed to make using the large scale of unmanned plane flank in itself For unmanned plane antenna model；Due to the VNF potato masher antennas of routine installed on unmanned plane body it is relatively difficult, therefore using conformal Array antenna, i.e., conformal array antenna is formed using the large scale of unmanned plane flank in itself, effective Enhanced Radiation Reduced Blast area, improve antenna Efficiency；
Step2, utilize FDTD analysis antenna performance：
Step3, antenna is asked for being fixed in coordinate system in a certain Electric Field Distribution on sphere in far field：For certain for the moment The earth axes S at quartergIn be in any attitude unmanned plane antenna, it is in body axis system SbIn position be it is fixed, because The electric field on its far field sphere under a certain working frequency is calculated with Finite-Difference Time-Domain Method in body axis system in this Distribution；
Step4, coordinate transform is recycled, the value tried to achieve above is converted into the electric field point when antenna is in any attitude Cloth：
Wherein, on condition that defining earth axes Oxg’yg’zg(Sg’) and body axis system Obxbybzb(Sb), by OXg′、 OYg' move to ObXg、ObYgIt is O to obtain with the earth axes that body mutually involvesbXgYgZg(Sg), then ObXgYgZg(Sg) with Obxbybzb(Sb) relation between two coordinate systems determines by three flight attitudes of unmanned plane：Yaw angle α, angle of pitch β and rolling Corner γ, its direction along corresponding coordinate dextrorotation all to switch to just, the primitive spin moment corresponding to note three above attitude transformation Battle array is respectively Lgbz(α)、Lgby(β) and Lgbx(γ), then from body axis system SbTo the earth axes S mutually involved with bodyg's Transformation matrix of coordinates Lgb (α, beta, gamma) is：
Lgb(α, β, γ)=Lgbz(α)Lgby(β)Lgbx(γ)
According to coordinate transform formula above, Electric Field Distribution of the unmanned plane antenna in any attitude can be tried to achieve；
Step5, the Electric Field Distribution on sphere is finally transformed into ground, draws pair of geographical coordinates and effective jamming power It should be related to：In earth axes SgIn, N points are a bit on the sphere electric field of far field, and the coordinate of N points is (xg,yg,zg), remember ObN with The intersection point N on ground0Coordinate is (xg0,yg0,zg0), then there is zg0=-h, h are unmanned plane during flying height, compare relation according to triangle etc. ：
Because the radio wave propagation mode between air-ground is considered as free-space propagation, far field electric field is inversely proportional with propagation distance, The electric field of N points is transformed into by N according to this relation0Point, its component array are (E0)gHave：
By above-mentioned conversion, the earth axes S mutually involved with body is tried to achievegMiddle z is x on groundg0And yg0With electric field E0 Corresponding relation, that is, obtained the field strength distribution on ground；
Wherein：Earth axes are OXg′Yg′Zg(Sg'), body axis system ObXbYbZb(Sb), by OXg′、OYg' translation To ObXg、ObYgIt is O to obtain with the earth axes that body mutually involvesbXgYgZg(Sg)；
E0It is that N points are mapped to N0The electric field of point；E is the electric field at origin O；(E)gFor electricity corresponding to the electric field E at origin O The component array of field vector；
Step6, regulation unmanned plane during flying posture carry out ground communication interference：According to the directionality of terrestrial receiving antenna, knot The field strength distribution on ground that the face of closing is calculated, for adjusting unmanned plane during flying posture, effectively, maximumlly disturb Enemy's communication system.
In the step Step2, concretely comprising the following steps for FDTD analysis antenna performance is utilized：
Step2.1, using conformal grid the day misaligned part of line boundary Fdtd Method grid is modified；
Step2.2, using thin wire Finite-Difference Time-Domain Method loading fine rule part is handled；
Step2.3, driving source set the cutler feed method using additional incentive；
Step2.4, recurrence calculation is carried out in finite region, using Berenger completely permutation absorbing boundary conditions；
Step2.5, calculate radiation pattern.
The specific implementation step of methods described is：
Step s1：Build unmanned plane antenna model.
Because the VNF potato masher antennas of routine are installed relatively difficult on unmanned plane body, therefore conformal array antenna is used, i.e., Conformal array antenna is formed using the large scale of unmanned plane flank in itself, effective Enhanced Radiation Reduced Blast area, improves antenna efficiency.Antenna Illustraton of model is as shown in Fig. 2 it is the asymmetric vibrator of a plane, working frequency 20-100MHz.Unmanned plane flank is in The wing, middle outer wing and outer wing composition, in side, wing centre section and middle outer wing lower surface form antenna according to certain shape copper foil The two poles of the earth.
Step s2：Utilize FDTD analysis antenna performance.
1st, the day misaligned part of line boundary Fdtd Method grid is modified using conformal grid, during calculating The influence of wing metal framework and dielectric material to antenna performance is considered, so as to reach the real working environment of antenna；
2nd, loading fine rule part is handled using thin wire Finite-Difference Time-Domain Method；
3rd, driving source sets the cutler feed method using additional incentive.In order to the broadband by the way that antenna is once calculated Characteristic, pumping signal use differential Gaussian pulse, and its calculation formula is：
Wherein n △ t are pulse width；π is pi value, and t0, τ are constant, determine the zero point and width of pulse respectively.
4th, due to the limitation of calculating capacity, recurrence calculation is carried out in finite region；In order to simulate the antenna spoke of open domain Process is penetrated, it is necessary to set absorbing boundary condition at the cutoff boundary of zoning.Inhaled using Berenger completely permutations Receive boundary condition；
5th, radiation pattern is calculated, in order to save amount of calculation, the output boundary wink that will be obtained first under pulse excitation State electric field and magnetic-field component carry out Fourier and convert to obtain frequency domain value, then again in the enterprising line frequency domain distance field of corresponding frequencies point Far zone field is obtained in extrapolation.
Step s3：Antenna is asked for being fixed in coordinate system in a certain Electric Field Distribution on sphere in far field.
For the earth axes S at a certain momentgIn be in any attitude unmanned plane antenna, it is in body axis system Sb In position be fixed, therefore it is calculated under a certain working frequency with Finite-Difference Time-Domain Method in body axis system Far field sphere on Electric Field Distribution, beWherein M be sphere on a bit, θ ObM and ybThe angle of axle,For ObM In zbObxbPlane projection and zbThe angle of axle.
Step s4：Coordinate transform is recycled, the value tried to achieve above is converted into the electric field point when antenna is in any attitude Cloth：
Wherein, on condition that defining earth axes Oxg’yg’zg(Sg’) and body axis system Obxbybzb(Sb), by OXg′、 OYg' move to ObXg、ObYgIt is O to obtain with the earth axes that body mutually involvesbXgYgZg(Sg), then ObXgYgZg(Sg) with Obxbybzb(Sb) relation between two coordinate systems determines by three flight attitudes of unmanned plane：Yaw angle α, angle of pitch β and rolling Corner γ, all to switch to just along corresponding coordinate dextrorotation, unmanned plane reference axis figure is as shown in Figure 3 in its direction.Remember three above posture Primitive spin matrix corresponding to angular transformation is respectively Lgbz(α)、Lgby(β) and Lgbx(γ), then from body axis system SbArrive and machine The earth axes S that body phase involvesgTransformation matrix of coordinates Lgb(α, beta, gamma) is：
Lgb(α, β, γ)=Lgbz(α)Lgby(β)Lgbx(γ)
According to coordinate transform formula above, Electric Field Distribution of the unmanned plane antenna in any attitude can be tried to achieve.
Step s5：The Electric Field Distribution on sphere is finally transformed into ground, draws geographical coordinates and effective jamming power Corresponding relation：
In earth axes SgIn, N points are a bit on the sphere electric field of far field, and the coordinate of N points is (xg,yg,zg), remember ON With the intersection point N on ground0Coordinate is (xg0,yg0,zg0), then there is zg0=-h (h is unmanned plane during flying height), compares according to triangle etc. Relation obtains：
Because the radio wave propagation mode between air-ground is considered as free-space propagation, far field electric field is inversely proportional with propagation distance, The electric field of N points is transformed into by N according to this relation0Point, its component array are (E0)gHave：
By above-mentioned conversion, the earth axes S mutually involved with body is tried to achievegMiddle z is x on groundg0And yg0With electric field E0 Corresponding relation, that is, obtained the field strength distribution on ground.Due to the directionality of reception antenna, useful interference field strength is E0 A certain component E0', so interference effective coverage is exactly E on ground0' it is more than the region of certain value.
Simulated experiment sets ground receiver as vertical polarized antenna, and unmanned plane uses the conformal array antenna mould of Fig. 2 designs Type, its flying height are 2500m, antenna input power 110W, working frequency 60MHz, and unmanned plane is aligned with different postures The data that the ground receiver of lower section carries out Communication Jamming are as shown in table 1：
The UAV Communication interference experiment result of table 1
Experiment and instance data more than are analyzed, it can be seen that when unmanned plane during flying posture difference, its Communication Jamming Radius is also different, and when left bank, which is flown, carries out Communication Jamming, the communication radius of interference is maximum, flat to fly to take second place, right bank flight It is worst.
Wherein earth axes Oxg’yg’zg(Sg’) be in Fig. 3 letter represent be indicated, wherein Oxg’yg’zgIt is Three-dimensional coordinate, (Sg’) earth axes code name is represented, g represents ground, corresponding (Sb) it is body axis system code name, b represents machine Body；Component array is (E0) g, wherein E0It is electric field code name, (E0) g represent electric field intensity component array.
Earth axes Ox defined in Fig. 3g’yg’zg(Sg’) and body axis system Obxbybzb(Sb), by Oxg’、Oyg’Translation To Obxg、Obyb, so as to obtain the earth axes O mutually involved with bodybxgygzg(Sg).If ObxbIn xgObygThe projection of plane For Obxb', ObzbIn xgObzgPlane is projected as Obzb', ObXgYgZg(Sg) and Obxbybzb(Sb) pass between two coordinate systems System can be determined by three flight attitudes of unmanned plane：Yaw angle α, angle of pitch β and roll angle γ, its direction along corresponding all to sit Mark dextrorotation switchs to just, M0Point is postulated point on the sphere of far field.
Above in conjunction with accompanying drawing to the present invention embodiment be explained in detail, but the present invention be not limited to it is above-mentioned Embodiment, can also be before present inventive concept not be departed from those of ordinary skill in the art's possessed knowledge Put that various changes can be made.

Claims (2)

  1. A kind of 1. UAV Communication interference countercheck based on real-time embedded control system, it is characterised in that：It is described to be based on The UAV Communication interference countercheck of real-time embedded control system comprises the following steps that：
    Step1, structure unmanned plane antenna model：Conformal array antenna, which is formed, using the large scale of unmanned plane flank in itself is used as nothing Man-machine antenna model；
    Step2, utilize FDTD analysis antenna performance：
    Step3, antenna is asked for being fixed in coordinate system in a certain Electric Field Distribution on sphere in far field：For a certain moment The earth axes S mutually involved with bodygIn be in any attitude unmanned plane antenna, it is in body axis system SbIn position It is fixed, therefore its far field ball under a certain working frequency is calculated with Finite-Difference Time-Domain Method in body axis system Electric Field Distribution on face；
    Step4, coordinate transform is recycled, the value tried to achieve above is converted into Electric Field Distribution when antenna is in any attitude：
    Wherein, on condition that defining earth axes Oxg’yg’zg(Sg’) and body axis system Obxbybzb(Sb), by OXg′、OYg' flat Move to ObXg、ObYgIt is O to obtain with the earth axes that body mutually involvesbXgYgZg(Sg), then ObXgYgZg(Sg) and Obxbybzb (Sb) relation between two coordinate systems determines by three flight attitudes of unmanned plane：Yaw angle α, angle of pitch β and roll angle γ, All to switch to just along corresponding coordinate dextrorotation, the primitive spin matrix corresponding to note three above attitude transformation is respectively in its direction Lgbz(α)、Lgby(β) and Lgbx(γ), then from body axis system SbTo the earth axes S mutually involved with bodygCoordinate transform Matrix L gb (α, beta, gamma) is：
    Lgb(α, β, γ)=Lgbz(α)Lgby(β)Lgbx(γ)
    According to coordinate transform formula above, Electric Field Distribution of the unmanned plane antenna in any attitude can be tried to achieve；
    Step5, the Electric Field Distribution on sphere is finally transformed into ground, show that geographical coordinates is corresponding with effective jamming power and close System：In earth axes SgIn, N points are a bit on the sphere electric field of far field, and the coordinate of N points is (xg,yg,zg), remember ObN and ground Intersection point N0Coordinate is (xg0,yg0,zg0), then there is zg0=-h, h are unmanned plane during flying height, are obtained according to triangle etc. than relation：
    <mrow> <msub> <mi>x</mi> <mrow> <mi>g</mi> <mn>0</mn> </mrow> </msub> <mo>=</mo> <mfrac> <mrow> <mo>-</mo> <mi>h</mi> </mrow> <msub> <mi>z</mi> <mi>g</mi> </msub> </mfrac> <msub> <mi>x</mi> <mi>g</mi> </msub> <mo>,</mo> <msub> <mi>y</mi> <mrow> <mi>g</mi> <mn>0</mn> </mrow> </msub> <mo>=</mo> <mfrac> <mrow> <mo>-</mo> <mi>h</mi> </mrow> <msub> <mi>z</mi> <mi>g</mi> </msub> </mfrac> <msub> <mi>y</mi> <mi>g</mi> </msub> </mrow>
    Because the radio wave propagation mode between air-ground is considered as free-space propagation, far field electric field is inversely proportional with propagation distance, according to The electric field of N points is transformed into N by this relation0Point, its component array are (E0)gHave：
    <mrow> <msub> <mrow> <mo>(</mo> <msub> <mi>E</mi> <mn>0</mn> </msub> <mo>)</mo> </mrow> <mi>g</mi> </msub> <mo>=</mo> <mfrac> <mrow> <mi>O</mi> <mi>N</mi> </mrow> <mrow> <msub> <mi>O</mi> <mi>b</mi> </msub> <msub> <mi>N</mi> <mn>0</mn> </msub> </mrow> </mfrac> <msub> <mrow> <mo>(</mo> <mi>E</mi> <mo>)</mo> </mrow> <mi>g</mi> </msub> <mo>=</mo> <mfrac> <msub> <mi>z</mi> <mi>g</mi> </msub> <mrow> <mo>-</mo> <mi>h</mi> </mrow> </mfrac> <msub> <mrow> <mo>(</mo> <mi>E</mi> <mo>)</mo> </mrow> <mi>g</mi> </msub> </mrow>
    By above-mentioned conversion, the earth axes S mutually involved with body is tried to achievegMiddle z is x on groundg0And yg0With electric field E0Pair It should be related to, that is, obtain the field strength distribution on ground；
    Wherein：Earth axes are OXg′Yg′Zg(Sg'), body axis system ObXbYbZb(Sb), by OXg′、OYg' move to ObXg、ObYgIt is O to obtain with the earth axes that body mutually involvesbXgYgZg(Sg)；
    E0It is that N points are mapped to N0The electric field of point；E is the electric field at origin O；(E)gSweared for electric field corresponding to the electric field E at origin O The component array of amount；
    Step6, regulation unmanned plane during flying posture carry out ground communication interference：According to the directionality of terrestrial receiving antenna, with reference to upper The field strength distribution on ground that face is calculated, for adjusting unmanned plane during flying posture, progress is effective, maximumlly disturbs enemy Communication system.
  2. 2. the UAV Communication interference countercheck according to claim 1 based on real-time embedded control system, it is special Sign is：In the step Step2, concretely comprising the following steps for FDTD analysis antenna performance is utilized：
    Step2.1, using conformal grid the day misaligned part of line boundary Fdtd Method grid is modified；
    Step2.2, using thin wire Finite-Difference Time-Domain Method loading fine rule part is handled；
    Step2.3, driving source set the cutler feed method using additional incentive；
    Step2.4, recurrence calculation is carried out in finite region, using Berenger completely permutation absorbing boundary conditions；
    Step2.5, calculate radiation pattern.
CN201510508339.7A 2015-08-18 2015-08-18 A kind of UAV Communication interference countercheck based on real-time embedded control system Active CN105119683B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201510508339.7A CN105119683B (en) 2015-08-18 2015-08-18 A kind of UAV Communication interference countercheck based on real-time embedded control system

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201510508339.7A CN105119683B (en) 2015-08-18 2015-08-18 A kind of UAV Communication interference countercheck based on real-time embedded control system

Publications (2)

Publication Number Publication Date
CN105119683A CN105119683A (en) 2015-12-02
CN105119683B true CN105119683B (en) 2018-03-06

Family

ID=54667591

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201510508339.7A Active CN105119683B (en) 2015-08-18 2015-08-18 A kind of UAV Communication interference countercheck based on real-time embedded control system

Country Status (1)

Country Link
CN (1) CN105119683B (en)

Families Citing this family (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US11032148B2 (en) * 2016-04-07 2021-06-08 Qualcomm Incorporated Managing network communication of an unmanned autonomous vehicle
CN105842717B (en) * 2016-05-19 2018-10-16 湖南矩阵电子科技有限公司 A kind of method and system of anti-unmanned plane
CN107037886A (en) * 2017-05-27 2017-08-11 成都索微通讯技术有限公司 A kind of system and its method of work extracted for human action
CN108494521B (en) * 2018-03-23 2020-07-17 青岛中邦防务智能装备有限公司 Enemy unmanned ship capturing method and system based on unmanned aerial vehicle interference
CN108573629B (en) * 2018-04-13 2020-06-09 扬州宇安电子科技有限公司 Unmanned aerial vehicle airborne radar confrontation simulation training system and training method
CN110311816B (en) * 2019-06-28 2021-08-13 上海交通大学 VNF (virtual network function) placement method for VNF co-location interference sensing in NFV (network virtual function) environment
CN111614418A (en) * 2020-05-28 2020-09-01 中国人民解放军陆军工程大学 Unmanned aerial vehicle data chain electromagnetic environment effect prediction method and device

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102426458A (en) * 2011-11-28 2012-04-25 北京航空航天大学 Ground control system applicable to rotor-wing unmanned aerial vehicle
WO2014209465A2 (en) * 2013-04-09 2014-12-31 Raytheon Company Methods and apparatuses for allocating electromagnetic-spectrum jamming assets
CN104573376A (en) * 2015-01-22 2015-04-29 北京航空航天大学 Method of calculating far extrapolation of transient field of electromagnetic scattering through finite difference time domain
CN104656664A (en) * 2015-02-16 2015-05-27 南京航空航天大学 Vehicle-mounted multi-rotor unmanned helicopter landing guidance control system and guidance control method

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2005020445A2 (en) * 2003-08-05 2005-03-03 University Of Hawai'i Microwave self-phasing antenna arrays for secure data transmission & satellite network crosslinks

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102426458A (en) * 2011-11-28 2012-04-25 北京航空航天大学 Ground control system applicable to rotor-wing unmanned aerial vehicle
WO2014209465A2 (en) * 2013-04-09 2014-12-31 Raytheon Company Methods and apparatuses for allocating electromagnetic-spectrum jamming assets
CN104573376A (en) * 2015-01-22 2015-04-29 北京航空航天大学 Method of calculating far extrapolation of transient field of electromagnetic scattering through finite difference time domain
CN104656664A (en) * 2015-02-16 2015-05-27 南京航空航天大学 Vehicle-mounted multi-rotor unmanned helicopter landing guidance control system and guidance control method

Also Published As

Publication number Publication date
CN105119683A (en) 2015-12-02

Similar Documents

Publication Publication Date Title
CN105119683B (en) A kind of UAV Communication interference countercheck based on real-time embedded control system
CN106407723B (en) The determination method of sparse arrangement array antenna exciting current amplitude towards Sidelobe
CN105511296A (en) Inner field radiation type complex electromagnetic environment building method and simulation system
CN104485513B (en) Real-time integrating method for wideband antenna array
CN104063544B (en) Near-field effect error analysis method
CN104980236B (en) The method of testing of low orbit satellite Ka frequency range data receiving system tracking performances
CN107888241A (en) A kind of conformal circular polarisation phase array antenna beam composition algorithm of curved surface
CN104733850B (en) Artificial tensor impedance skin antenna and its implementation based on holographic modulation
CN107657070A (en) Decoupling method and antenna based on the conformal array beamses of genetic algorithm frequency diversity
CN104243070A (en) Wireless channel realtime simulator and simulation method oriented to deep space communication
CN106469850A (en) A kind of Thickness Design Method of antenna house
CN202817189U (en) S-band broadband wide beam aircraft antenna
CN103872463A (en) Small S-band passive phased array antenna with high radiation efficiency
CN104536291A (en) Method for simulating influence of elastic vibration on seeker measuring signal based on RF system
CN107017468A (en) A kind of cell layout of conformal array antenna of fuselage and Electromagnetic Desigu Method
CN108173578A (en) Array antenna simulates multi-beam shaping method
CN107333290A (en) A kind of antenna for base station unit for electrical property parameters monitoring method
CN107026331A (en) A kind of stepped Aperture distribution design method transmitted for microwave wireless energy
CN109286078A (en) Null tone domain gradient Meta Materials and its design method
CN104076336B (en) Based on radar folding face battle array rule numeral submatrix and the forming method of/difference beam
CN204103045U (en) A kind of quadrangular pyramid trumpet array antenna for 2.65GHz ~ 2.91GHz frequency range
Semenikhin et al. Active cancellation of radar cross section of large aircraft using conformai 2-port magnetic T-shaped antenna
CN113218379B (en) Lithium battery unmanned aerial vehicle electromagnetic interference suppression method coupled with flight trajectory
Litinskaya et al. A Ku-band low-profile wide-angle scanning antenna array with combined beam steering
CN108828508B (en) A kind of method for analyzing performance of the direct location model of over the horizon radiation source

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant