CN107066693A - The spaceborne AIS reconnaissance signals simulation system of multi-channel multi-target - Google Patents
The spaceborne AIS reconnaissance signals simulation system of multi-channel multi-target Download PDFInfo
- Publication number
- CN107066693A CN107066693A CN201710148410.4A CN201710148410A CN107066693A CN 107066693 A CN107066693 A CN 107066693A CN 201710148410 A CN201710148410 A CN 201710148410A CN 107066693 A CN107066693 A CN 107066693A
- Authority
- CN
- China
- Prior art keywords
- mrow
- msub
- signal
- ais
- ship
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S19/00—Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
- G01S19/38—Determining a navigation solution using signals transmitted by a satellite radio beacon positioning system
- G01S19/39—Determining a navigation solution using signals transmitted by a satellite radio beacon positioning system the satellite radio beacon positioning system transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
- G01S19/42—Determining position
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S19/00—Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
- G01S19/38—Determining a navigation solution using signals transmitted by a satellite radio beacon positioning system
- G01S19/39—Determining a navigation solution using signals transmitted by a satellite radio beacon positioning system the satellite radio beacon positioning system transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
- G01S19/52—Determining velocity
Landscapes
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Networks & Wireless Communication (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- Position Fixing By Use Of Radio Waves (AREA)
Abstract
The present invention relates to the spaceborne AIS reconnaissance signals simulation system of multi-channel multi-target.Be it is a kind of to spaceborne AIS multichannels reconnaissance system in orbit when receive the system simulated of dynamic ais signal of multiple target aliasing.It is digital calculating that present system, which produces two paths of signals, all filtering, the frequency conversion of four passage ais signals, and interchannel uniformity is good, and the spuious small, harmonic component of signal is small；Clock sampling rate is high, and the temporal resolution for generating ais signal is high, and the control accuracy of time delay can be of about 8us magnitudes；It can produce that AIS target numbers are more, the continuous simulation time is long（More than more than five minutes）, any scene of 0~10000 ship in same visual field can be emulated.The AIS radiofrequency signals of the dynamic change multiple target aliasing with Doppler frequency shift, amplitude scintillation etc. caused by relative motion are received when can realize the spaceborne scouting of generation true to nature.
Description
Technical field
The invention belongs to signal imitation and simulation technical field, it is related to signal imitation system, it is particularly a kind of that satellite is detectd
Examine the multi-channel multi-target simulation transposition of ship automatic identification system (AIS) signal of the various ship transmittings in sea.
Background technology
Ship automatic identification system (Automatic Identification System, AIS) is a kind of in very high frequency(VHF) sea
Upper cellular band uses self-organizing time division multiple acess access way Auto broadcast and receives the information such as Ship dynamic situation, static state to realize
The system of identification, monitoring and communication.Realized extensively on all kinds of big-and-middle crafts at present and short disappear is received and dispatched between ship-to-shore, ship
The functions such as breath, interaction text, navigation collision prevention.Ship effective monitoring is managed to realize, spaceborne ship is known automatically in the recent decade
Other system is arisen at the historic moment, and naval vessel progress all-weather effectively to global marine site, the monitoring of big visual field can be realized to China
The information battle state display and early warning on surrounding sea areas or even global marine site naval vessel.
Because AIS systems employ ITU-R M.1371-5 self-organizing time division multiple acess (SOTDMA) specified in technical specification
Transmission signal causes satellite reception end signal aliasing phenomenon occur to ship between communication system, different districts simultaneously.In addition,
Spaceborne AIS is received receives also following what time different from onboard AIS：One is due to that satellite causes to mix with ship versus high motion
There is Doppler frequency shift between folded signal source；Two be due to that different distance power attenuation causes between the aliasing signal source that receives
There is difference power；Three be due to that ship and reception antenna position difference cause there is time delay between the aliasing signal source received.
Therefore the time domain and the multiple aliasing of frequency domain, Doppler frequency shift of multiple naval vessel cell signals will be produced during satellite reconaissance, amplitude rises
Phenomenon is lied prostrate, the performance to spaceborne AIS reconnaissance systems causes to have a strong impact on.Current most spaceborne AIS systems such as Norway is sent out for 2010
Two penetrated have the SDI of spaceborne AIS reception processings function：AISSat-1 and NORAIS etc. is with one or more
Many receiving channels.After being gone up to the sky in orbit due to satellite launch, it is impossible to find some standard scenes to each of spaceborne AIS systems
Item technical indicator is verified and demarcated, and in order to carry out assessment that is accurate, quantifying to the performance of spaceborne AIS systems, having must
Develop the spaceborne AIS reconnaissance signals simulator of special multi-channel multi-target.
Chinese patent application 201610786782.5 proposes the AIS digital simulation sides in a kind of Marine Electronics application system
Method, but this method is mainly used in emulating the data of single AIS targets, it is impossible to being received when producing spaceborne scouting has relative motion
The dynamic change AIS radiofrequency signals of caused Doppler frequency shift, amplitude scintillation etc., can not more emulate the spaceborne of generation multiple target
AIS scouts the aliasing signal of multichannel receiver coherent, and these exactly influence the most critical of spaceborne AIS reconnaissance systems performance
Several factors.
The content of the invention
The technical problem to be solved in the present invention is received when being and proposing one kind to spaceborne AIS multichannels reconnaissance system in orbit
The system that the dynamic ais signal of multiple target aliasing is simulated.
The technical scheme is that, a kind of spaceborne AIS reconnaissance signals simulation system of multi-channel multi-target, including one
CPCI industrial computers (or common industrial computer, all-purpose computer), signal produce card, two pieces of DAC subcards,
The CPCI industrial computers include computer motherboard, cabinet, power supply, keyboard, mouse, display, and operation thereon has scene number
Software is loaded according to software, data are produced；
The contextual data generation software is an application software run on the CPCI industrial computers, and it sets according to user
Scenario simulation generates the multiple target multichannel AIS baseband signal data files of satellite reception.For many ships hair on sea
Ais signal scene is penetrated, because each ship is all launched according to SOTDMA mechanism, therefore can mould in accordance with the following steps
Intend：The first step, read user's selection includes satellite position, the satellite ephemeris file of speed, on the visual visual field sea of satellite
Randomly generate position, speed, operational configuration of N ship etc.；Second step, first ship is according to ITU-RM.1371-5 technical specifications
According to its headway and Determines its transmission interval Tr, a time slot moment progress is randomly choosed in its transmission interval Tr
Transmitting, then the transmission time slot according to transmission interval and first, determines the slot transmission moment all behind this ship successively；3rd
Step, to follow-up n-th ship, it is first determined whether having ship in its visual field cell range；4th step, kth ship random selection is first
Secondary transmission time slot simultaneously determines following slots according to transmission interval；5th step, judges all transmission time slots and its visual field of kth ship
Ship has distributed whether time slot conflicts in cell；If not conflicting, next ship generation is carried out, if conflict, comes back to the 4th
Step ..., until not conflicting；6th step, judges whether the time slot distribution of all ships finishes, i.e. whether k is more than N, if it is not, then repeating
3rd step, if so, the timeslot number then recorded the emission time of every ship and taken.
According to timeslot number and emission time, the AIS data of the time slot are generated, detailed process is as follows：
According to k-th of radiation source, i-th of timeslot number corresponding time slot moment tk,iAnd ITU-R M.1371-4 technical specifications, it is raw
Into sample rate fs1The AIS baseband signals of=GMSK modulation for being 38.4kHzWhereinExpression is rounded
Computing.The signal data of multiple time slots is multiplied by after boat-carrying transmitter antenna gain (dBi), spaceborne receiving antenna gain and Doppler frequency shift
Time domain superposition is carried out according to the corresponding moment of each time slot, you can obtain spaceborne AIS baseband signal datas：
Above formula C is P in the amplitude normalization factor, formulark,iFor the received signal power of each time slot, fdk,iFor each time slot
Doppler frequency shift, τk,iFor the time delay of each time slot, φk,iFor the phase difference of the interchannel of each time slot, Ik
For slot index, n is superposed signal quantity, and j is superposed signal sequence number；
The calculation respectively measured is as follows：
A) calculating of power (amplitude)
The P in (1) formulark,iFor received signal power space attenuation amount, according to reconnaissance equation, its calculating power is as follows
(2) P in formulaTFor known onboard AIS transmission power, λ is ais signal wavelength, Gtk,iExist for onboard AIS transmitting antenna
tk,iGain on moment satellite pointing direction, the empirical formula of typical onboard AIS antenna gain is as follows
Wherein θt,iFor the local elevation angle of satellite place ship, its calculation formula is
Wherein xsk,i=[xsk,i,ysk,i,zsk,i]TAnd xtk,i=[xtk,i,ytk,i,ztk,i]TRespectively satellite and kth ship is in tk,i
Three-dimensional position of the moment ground under solid (ECF) coordinate system of heart, " | | | | " represent 2 norms, rk,iFor the distance of satellite and ship：
rk,i=| | xsk,i-xtk,i|| (5)
(6) G in formulark,iIt is spaceborne reception antenna in tk,iGain on moment kth ship pointing direction, its typical experience
Calculation formula is
Wherein uak,iThe sensing unit vector for being satellite antenna in ECF coordinate systems, typically according to the attitude of satellite and install can be with
Previously known and hypothesis is fixed in celestial body coordinate system；
B) calculating of Doppler frequency shift
(7) f in formuladk,iFor Doppler frequency shift, its calculation formula is as follows
For satellite and Relative motion of ship speed, λk,iFor transmitted signal wavelengths；
C) calculating of binary channels phase difference
Because the line direction of dual-channel antenna is consistent with satellite motion direction, the phase difference that can be calculated between two antennas is
WhereinD is distance, v between passagesk,iFor speed horizontal component, vsi,iIt is vertical for speed
Component.
D) calculating of signal time delay
According to each target reach satellite distance, can calculate obtain each time slot signal signal time delay be
Wherein c ' is the light velocity；
The parameter that each time slot is calculated in (1) formula can be generated according to (2)~(9) formula, many mesh can be calculated by (1) formula
The signal waveform data of superposition is marked, signal waveform data binary file is written as to disk file, you can obtain baseband signal
Data file.
The data loading software is operate in another application software on CPCI industrial computers, and its major function is selection
Scene generates the data file that software is produced with baseband signal, is downloaded to by CPCI interfaces on signal generation motherboard, for broadcasting
Put generated data file waveform.Because the baseband signal data file that the contextual data generates software generation is stored
On disk, therefore the same or different pieces of information file of broadcasting can be iteratively repeated by data loading software, for examining repeatedly
Test and check the performance of spaceborne AIS reconnaissance systems.
It is one piece of CPCI plug-in card (or PCI plug-in cards) that the signal, which produces card, signal produce plate be based on 66MHz, 64
The CPCI plug-in cards of pci bus, for connecting DAC subcards and completing Some digital signal processing function.Include clock electricity on support plate
Road, power module, signal processing module, memory, and cpci bus control module etc., and pass through FPGA and the company of expansion slot
Connect to external expansion multiple interfaces.Bottom plate receives the signal sampling point sent from daughter board or data, the carry out portion in FPGA and dsp chip
The real-time processing of point parameter, or store and stay post processing of doing things onto the DDR memory on plate.It will be located by the cpci bus of 64
Reason result or initial data pass out to display or processing in industrial control host.Processing module is a piece of ADI companies on signal generation card
ADSP-TS201S dsp chips and two panels Xilinx companies XC5VSX95T fpga chips.There is number between DSP and FPGA
Connected according to bus, while also having the high speed interconnection of definable between two panels FPGA, realize direct data between FPGA
Interaction；FPGA is connected with four external data expansion interfaces on bottom plate, and the friendship with external data is realized by piggy-back daughter board mode
Change；There is a DSP expansion interface, FPGA, DSP and memory are increased by piggy-back daughter board mode, so that the money of expansion board clamping
Source and disposal ability；There is self-defined bus connection between FPGA and CPCI expansion slot J5, for realizing between other CPCI boards
Self-defined connection.All self-defined bus may be defined as LVDS differential signal lines or LVCMOS single-ended signal lines, depend on
In FPGA design；Clock integration module on plate can carry out parameter setting by DSP, clock frequency and phase required for obtaining
Position relation.
The DAC subcards be DAC subcards to be connected in the form of high speed connector interface on data processing support plate, including one
Individual DAC chip, the data that can produce FPGA pass through digital to analog conversion (D/A) directly output AIS radio frequency analog signals.DA is changed
Device is 16 DAC chips of DAC5687 models.The chip operation clock is 500MHz, the dynamic range with more than 80dB.Carry
Plate provides power supply and clock by high speed connector for daughter board, and signal sampling point is transferred into DAC module realizes conversion.
In view of AIS different channels centre frequency be respectively 161.975MHz/162.025MHz (87B/88B channels) and
156.775MHz/156.825MHz (75/76 channel), signal bandwidth is all 25kHz, in 161.975MHz/162.025MHz letters
Differed between number and 50kHz is also differed between 50kHz, 156.775MHz/156.825MHz, contextual data generation of the present invention
Software produces the sample rate of baseband signal data, and representative value selection is that sample rate is 122.07kHz, can be with undistorted realization pair
Imitated while two channels of ais signal 161.975MHz/162.025MHz or two channels of 156.775MHz/156.825MHz
True generation.
In order to be simulated while the signal produces and realizes above-mentioned AIS channel signals on card, signal can be used
The technology of waveform multistage zero padding interpolation and three-level Vector Modulation.I.e. described contextual data generation software is same producing data file
When produce the base band of two channels of 161.975MHz/162.025MHz and two channels of person 156.775MHz/156.825MHz
Data, contextual data generation software exports fs=122.07kHz baseband signal data, then loads on signal generation card card
FPGA according to multistage zero padding interpolation, 4096 times of interpolations are obtained altogether and carry out up-conversion again, the direct radio frequency to ais signal is realized
Simulation.
In order to reduce the demand to FPGA resource, the data of 87B/88B channels are created within sample rate for fs=
In 122.07kHz same I/Q data, by amounting to 8 × 8=64 times of filtering interpolation after two-stage 8 times of interpolation, LPFs
Afterwards, by a complex carrier signal complex modulation to 2.6MHz digital intermediate frequency；The data of 75/76 channel are created within fs=
In 122.07kHz another I/Q data, amounted to by 8 times of interpolation of two-stage, LPF after 64 times of filtering interpolations by one
Complex carrier signal complex modulation is on -2.6MHz digital intermediate frequency, and it is 7.8125MHz's that two complex carriers synthesize a sample rate again
In complex carrier, then 8 times, 4 times, 2 times of interpolation, LPF 8 × 4 × 2=64 times of filtering interpolations of formation, generation by three-level
The zero intermediate frequency i/q signal of 500MHz samplings, then on carrier wave by Vector Modulation to 159.4MHz, formed two-way, 4 channels/
The data signal of the 500MHz samplings on road, is controlled the D/A chips on two-way DAC subcards to carry out digital to analog conversion by FPGA, produced
Examine the dynamic spaceborne ais signal of binary channels needed for satellite AIS reconnaissance equipments.
Following beneficial effect can be reached using the present invention：
1st, received when present system can realize the spaceborne scouting of generation true to nature have Doppler frequency shift caused by relative motion,
The AIS radiofrequency signals of the dynamic change multiple target aliasing of amplitude scintillation etc.；
2nd, it is digital calculating that present system, which produces two paths of signals, all filtering, the frequency conversion of four passage ais signals, is led to
Duplication between channels is good, and the spuious small, harmonic component of signal is small；
3rd, clock sampling rate of the invention is high, and the temporal resolution for generating ais signal is high, and the control accuracy of time delay is reachable
About 8us magnitudes；
4th, the present invention can produce that AIS target numbers are more, continuous simulation time length (being more than more than five minutes), can emulate same
Any scene of 0~10000 ship in visual field.
Brief description of the drawings
Fig. 1：The spaceborne AIS reconnaissance signals simulation system composition frame chart of multi-channel multi-target of the present invention；
Fig. 2：The spaceborne AIS reconnaissance signals simulator workflow diagram of multi-channel multi-target；
Fig. 3：The composition connection block diagram of signal generation card hardware module；
Fig. 4：The hardware module composition connection block diagram of DAC subcards；
Fig. 5：Single channel signal waveform multistage zero padding interpolation and three-level Vector Modulation functional block diagram；
Fig. 6：The spaceborne AIS of typical case scouts simulating scenes two-dimension projection；
Fig. 7：The AIS baseband signal waveform figures of multiple target aliasing.
Embodiment
Below, with reference to the drawings and specific embodiments, the invention will be further described：
As shown in figure 1, the spaceborne AIS reconnaissance signals simulator of the multi-channel multi-target proposed, its by a CPCI industrial computer (or
The common industrial computer of person, all-purpose computer), signal produce the hardware such as card, two pieces of DAC subcards and the field run on CPCI industrial computers
Scape data produce the software compositions such as software, data loading software.Wherein contextual data produces software according to user's setting generation base
Baseband signal data hair is simultaneously passed through cpci bus by band signal data file, data loading software loading baseband signal data file
Being sent on signal generation card card, signal generation card card, there is FPGA, DDR memory, DSP etc. to calculate memory device, by base band
Signal generates the AIS rf digital signal data of two-way 500MHz samplings through multistage zero padding interpolation and three-level Vector Modulation,
Signal generation card card is connected with two pieces of DAC subcards by high speed connector and AIS rf digital signal data is sent into DAC
Card, D/A devices on DAC subcards are by AIS rf digital signals data conversion into analog signal output.
Clicked on as shown in Fig. 2 starting scene after user's start and being generated with baseband signal data after software, scene set parameter
Start emulation, software is to start to calculate visual range according to the satellite motion scene of setting, and is generated in visual range multiple
Ship Target, and SOTDMA time slot emulation is carried out, ais signal then is generated according to time slot as defined in each target, to multiple target
Signal is overlapped, and the baseband signal data of generation multiple target superposition is simultaneously stored on disk, then added by signal loading software
Carry corresponding signal data file output.The baseband signal data file of any history trendline also may be selected in signal loading software,
Play out or play back repeatedly.
As shown in figure 3, signal generation card is the CPCI plug-in cards based on 66MHz, 64 pci bus, for connecting DAC subcards
And complete the digital signal processing functions such as multistage zero padding interpolation, filtering, three-level Vector Modulation.Include clock circuit, electricity on support plate
Source module, signal processing module, memory, and cpci bus control module etc., and pass through the connection pair of FPGA and expansion slot
External expansion multiple interfaces.Bottom plate receives the signal sampling point sent from daughter board or data, and part ginseng is carried out in FPGA and dsp chip
Several real-time processing, or store to staying post processing of doing things on the DDR memory on plate.It will be handled and tied by the cpci bus of 64
Fruit or initial data pass out to display or processing in industrial control host.Main process modules are a piece of DSP and two panels FPGA on plate,
DSP type selectings are ADSP-TS201S, and FPGA type selectings are XC5VSX95T.Between DSP and FPGA there is data/address bus to connect, while two
Also there is the high speed interconnection of definable between piece FPGA, directly data interaction between FPGA is realized；FPGA connects on bottom plate
Four external data expansion interfaces are connected to, is realized and is exchanged with external data by piggy-back daughter board mode；There is a DSP extension to connect
Mouthful, FPGA, DSP and memory are increased by piggy-back daughter board mode, so that the resource and disposal ability of expansion board clamping；FPGA with
There is self-defined bus connection between CPCI expansion slot J5, for realizing the self-defined connection between other CPCI boards.All
Self-defined bus may be defined as LVDS differential signal lines or LVCMOS single-ended signal lines, depending on FPGA design；On plate
Clock integration module can carry out parameter setting by DSP, clock frequency and phase relation required for obtaining.
Fig. 4 provides DAC subcards design module frame chart, and DAC subcards are connected at data in the form of high speed connector interface
Manage on support plate.Support plate provides power supply and clock by high speed connector for daughter board, and signal sampling point is transferred into DAC module realization
Conversion, then smothing filtering is carried out by low pass filter.D/A converter is 16 DAC chips of DAC5687 models.Should
Chip operation clock is 500MHz, and low pass filter bandwidth is maximum up to 170MHz.
As shown in figure 5, in order to be simulated while the signal produces and realizes above-mentioned AIS channel signals on card, adopting
With the multistage zero padding interpolation of signal waveform as depicted and the technology of three-level Vector Modulation.I.e. described contextual data generation software
It is while producing two channels of 161.975MHz/162.025MHz and person 156.775MHz/ producing data file
Then the base band data of two channels of 156.825MHz, its sample rate f s=122.07kHz is loaded on signal generation card card
FPGA is obtained after 4096 times of interpolations altogether according to multistage zero padding interpolation and LPF, then carries out up-conversion, realizes and AIS is believed
Number direct radio frequency analog.In order to reduce the demand of FPGA resource, the data of 87B/88B channels are modulated at same sample rate
Sampled in fs=122.07kHz I/Q data, to amount to be formed after 64 times of filtering interpolations by 8 times of interpolation of two-stage, LPF
Rate 7.8125MHz digital i/q signal, by a complex carrier signal complex modulation to 2.6MHz digital intermediate frequency；75/76 channel
Data are modulated in another fs=122.07kHz I/Q data, amount to 64 times by 8 times of interpolation of two-stage, LPF
By in a complex carrier signal complex modulation to -2.6MHz digital intermediate frequency after filtering interpolation.Two complex carriers are added and synthesized again
One sample rate is 8 times in 7.8125MHz complex carrier, then by three-level, 4 times, 2 times of interpolation, LPF formation 8 ×
4 × 2=64 times of filtering interpolation, the zero intermediate frequency i/q signal of generation 500MHz samplings, then pass through the load of Vector Modulation to 159.4MHz
On ripple, the data signal that formation two-way, the 500MHz on 4 channel/roads sample passes through the D/A on FPGA control two-way DAC subcards
Chip carries out digital to analog conversion, produces the dynamic spaceborne ais signal of binary channels examined needed for satellite AIS reconnaissance equipments.
As shown in fig. 6, certain once emulates setting AIS satellite altitude 500km, once start shooting duration 300s, and it is eastern peaceful
In the scope of reconnaissance of foreign area operation, it is assumed that 4500 ships of setting are uniformly distributed in the ocean in the covering visual field of satellite reconaissance
On, all these ships are produced in spaceborne scouting for emulating, and what satellite reception was arrived has Doppler caused by relative motion
The AIS radiofrequency signals of the dynamic change multiple target aliasing of frequency displacement, amplitude scintillation etc..
As shown in fig. 7, being dynamically superimposed Doppler's frequency for 4500 ships that Fig. 6 scenes contextual data generates Software Create
The AIS baseband signals of the dynamic change multiple target aliasing of shifting, amplitude scintillation etc., the sample rate of the baseband signal is 122.07MHz.
Claims (3)
1. the spaceborne AIS reconnaissance signals simulation system of multi-channel multi-target, including CPCI industrial computer, signal produce card, two pieces
DAC subcards,
Operation has contextual data to produce software, data loading software on the CPCI industrial computers, it is characterised in that
The contextual data generation software simulates the multiple target multichannel AIS base band for generating satellite reception according to user's set scene
Signal data file, launches ais signal scene for many ships on sea, simulates in accordance with the following steps：The first step, reads
User selection includes satellite position, the satellite ephemeris file of speed, and N ship is randomly generated on the visual visual field sea of satellite
Position, speed, operational configuration；Second step, first ship according to ITU-R M.1371-5 technical specification according to its headway and
Determines its transmission interval Tr, one time slot moment of random selection is launched in its transmission interval Tr, then according to hair
Interval and first transmission time slot are penetrated, the slot transmission moment all behind this ship are determined successively；3rd step, to follow-up n-th
Ship, it is first determined whether having ship in its visual field cell range；4th step, kth ship randomly chooses transmission time slot first and root
Following slots are determined according to transmission interval；5th step, all transmission time slots and ship in its visual field cell for judging kth ship have been distributed
Whether time slot conflicts；If not conflicting, next ship generation is carried out, if conflict, comes back to the 4th step, until not conflicting；
6th step, judges whether the time slot distribution of all ships finishes, i.e. whether k is more than N, if it is not, the 3rd step is then repeated, if so, then remembering
The emission time of every ship and the timeslot number taken, according to timeslot number and emission time, generate the AIS data of the time slot under record；
The data loading software is to select scene to generate the data file that software is produced with baseband signal, downloads to signal generation
On motherboard, for playing generated data file waveform.
2. the spaceborne AIS reconnaissance signals simulation system of multi-channel multi-target according to claim 1, it is characterised in that described
It is the CPCI plug-in cards based on 66MHz, 64 pci bus that signal, which produces card, for connecting DAC subcards and completing Some digital signal
Processing function.
3. the spaceborne AIS reconnaissance signals simulation system of multi-channel multi-target according to claim 1, it is characterised in that according to
Timeslot number and emission time, generate the AIS data of the time slot, and detailed process is as follows：
According to k-th of radiation source, i-th of timeslot number corresponding time slot moment tk,iAnd ITU-R M.1371-5 technical specifications, it is raw
Into sample rate fs1The AIS baseband signals of=GMSK modulation for being 38.4kHzWhereinExpression is rounded
Computing, the signal data of multiple time slots is multiplied by after boat-carrying transmitter antenna gain (dBi), spaceborne receiving antenna gain and Doppler frequency shift
Time domain superposition is carried out according to the corresponding moment of each time slot, you can obtain spaceborne AIS baseband signal datas：
Above formula C is P in the amplitude normalization factor, formulark,iFor the received signal power of each time slot, fdk,iFor each time slot
Doppler frequency shift, τk,iFor the time delay of each time slot, φk,iFor the phase difference of the interchannel of each time slot, Ik
For slot index, n is superposed signal quantity, and j is superposed signal sequence number；
The calculation respectively measured is as follows：
A) calculating of power (amplitude)
The P in (1) formulark,iFor received signal power space attenuation amount, according to reconnaissance equation, its calculating power is as follows
<mrow>
<msub>
<mi>P</mi>
<mrow>
<mi>r</mi>
<mi>k</mi>
<mo>,</mo>
<mi>i</mi>
</mrow>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<msub>
<mi>P</mi>
<mi>T</mi>
</msub>
<msub>
<mi>G</mi>
<mrow>
<mi>t</mi>
<mi>k</mi>
<mo>,</mo>
<mi>i</mi>
</mrow>
</msub>
<msub>
<mi>G</mi>
<mrow>
<mi>r</mi>
<mi>k</mi>
<mo>,</mo>
<mi>i</mi>
</mrow>
</msub>
<msup>
<mi>&lambda;</mi>
<mn>2</mn>
</msup>
</mrow>
<msup>
<mrow>
<mo>(</mo>
<mn>4</mn>
<msub>
<mi>&pi;r</mi>
<mrow>
<mi>k</mi>
<mo>,</mo>
<mi>i</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>2</mn>
<mo>)</mo>
</mrow>
</mrow>
(2) P in formulaTFor known onboard AIS transmission power, λ is ais signal wavelength, Gtk,iIt is onboard AIS transmitting antenna in tk,i
Gain on moment satellite pointing direction, the empirical formula of typical onboard AIS antenna gain is as follows
<mrow>
<msub>
<mi>G</mi>
<mrow>
<mi>r</mi>
<mi>t</mi>
<mo>,</mo>
<mi>i</mi>
</mrow>
</msub>
<mo>=</mo>
<mn>0.964</mn>
<mfrac>
<mn>16</mn>
<mrow>
<mn>3</mn>
<mi>&pi;</mi>
</mrow>
</mfrac>
<msup>
<mi>sin</mi>
<mn>2</mn>
</msup>
<msub>
<mi>&theta;</mi>
<mrow>
<mi>t</mi>
<mo>,</mo>
<mi>i</mi>
</mrow>
</msub>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>3</mn>
<mo>)</mo>
</mrow>
</mrow>
Wherein θt,iFor the local elevation angle of satellite place ship, its calculation formula is
<mrow>
<msub>
<mi>&theta;</mi>
<mrow>
<mi>t</mi>
<mo>,</mo>
<mi>i</mi>
</mrow>
</msub>
<mo>=</mo>
<mfrac>
<mi>&pi;</mi>
<mn>2</mn>
</mfrac>
<mo>-</mo>
<mi>arccos</mi>
<mfrac>
<mrow>
<msup>
<mrow>
<mo>(</mo>
<msub>
<mi>x</mi>
<mrow>
<mi>s</mi>
<mi>k</mi>
<mo>,</mo>
<mi>i</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>x</mi>
<mrow>
<mi>t</mi>
<mi>k</mi>
<mo>,</mo>
<mi>i</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mi>T</mi>
</msup>
<msub>
<mi>x</mi>
<mrow>
<mi>t</mi>
<mi>k</mi>
<mo>,</mo>
<mi>i</mi>
</mrow>
</msub>
</mrow>
<mrow>
<msub>
<mi>r</mi>
<mrow>
<mi>k</mi>
<mo>,</mo>
<mi>i</mi>
</mrow>
</msub>
<mo>|</mo>
<mo>|</mo>
<msub>
<mi>x</mi>
<mrow>
<mi>t</mi>
<mi>k</mi>
<mo>,</mo>
<mi>i</mi>
</mrow>
</msub>
<mo>|</mo>
<mo>|</mo>
</mrow>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>4</mn>
<mo>)</mo>
</mrow>
</mrow>
Wherein xsk,i=[xsk,i,ysk,i,zsk,i]TAnd xtk,i=[xtk,i,ytk,i,ztk,i]TRespectively satellite and kth ship is in tk,i
Three-dimensional position of the moment ground under solid (ECF) coordinate system of heart, | | | | represent 2 norms, rk,iFor the distance of satellite and ship：
rk,i=| | xsk,i-xtk,i|| (5)
(6) G in formulark,iIt is spaceborne reception antenna in tk,iGain on moment kth ship pointing direction, its typical experience meter
Calculating formula is
<mrow>
<msub>
<mi>G</mi>
<mrow>
<mi>r</mi>
<mi>k</mi>
<mo>,</mo>
<mi>i</mi>
</mrow>
</msub>
<mo>=</mo>
<msup>
<mrow>
<mo>(</mo>
<mfrac>
<mrow>
<msup>
<mrow>
<mo>(</mo>
<msub>
<mi>x</mi>
<mrow>
<mi>s</mi>
<mi>k</mi>
<mo>,</mo>
<mi>i</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>x</mi>
<mrow>
<mi>t</mi>
<mi>k</mi>
<mo>,</mo>
<mi>i</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mi>T</mi>
</msup>
<msub>
<mi>u</mi>
<mrow>
<mi>a</mi>
<mi>k</mi>
<mo>,</mo>
<mi>i</mi>
</mrow>
</msub>
</mrow>
<msub>
<mi>r</mi>
<mrow>
<mi>k</mi>
<mo>,</mo>
<mi>i</mi>
</mrow>
</msub>
</mfrac>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>6</mn>
<mo>)</mo>
</mrow>
</mrow>
Wherein uak,iThe sensing unit vector for being satellite antenna in ECF coordinate systems, typically according to the attitude of satellite and install can be with thing
First known and hypothesis is fixed in celestial body coordinate system,
B) calculating of Doppler frequency shift
(7) f in formuladk,iFor Doppler frequency shift, its calculation formula is as follows
<mrow>
<msub>
<mi>f</mi>
<mrow>
<mi>d</mi>
<mi>k</mi>
<mo>,</mo>
<mi>i</mi>
</mrow>
</msub>
<mo>=</mo>
<mo>-</mo>
<mfrac>
<mrow>
<msubsup>
<mi>v</mi>
<mrow>
<mi>k</mi>
<mo>,</mo>
<mi>i</mi>
</mrow>
<mi>T</mi>
</msubsup>
<mrow>
<mo>(</mo>
<msub>
<mi>x</mi>
<mrow>
<mi>s</mi>
<mi>k</mi>
<mo>,</mo>
<mi>i</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>x</mi>
<mrow>
<mi>t</mi>
<mi>k</mi>
<mo>,</mo>
<mi>i</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msub>
<mi>&lambda;</mi>
<mrow>
<mi>k</mi>
<mo>,</mo>
<mi>i</mi>
</mrow>
</msub>
<msub>
<mi>r</mi>
<mrow>
<mi>k</mi>
<mo>,</mo>
<mi>i</mi>
</mrow>
</msub>
</mrow>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>7</mn>
<mo>)</mo>
</mrow>
</mrow>
For T moment ship and satellite relative velocity component, λk,iFor signal wavelength；
C) calculating of binary channels phase difference
Because the line direction of dual-channel antenna is consistent with satellite motion direction, the phase difference that can be calculated between two antennas is
<mrow>
<msub>
<mi>&phi;</mi>
<mrow>
<mi>k</mi>
<mo>,</mo>
<mi>i</mi>
</mrow>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<mn>2</mn>
<msubsup>
<mi>&pi;du</mi>
<mrow>
<mi>v</mi>
<mi>k</mi>
<mo>,</mo>
<mi>i</mi>
</mrow>
<mi>T</mi>
</msubsup>
<mrow>
<mo>(</mo>
<msub>
<mi>x</mi>
<mrow>
<mi>s</mi>
<mi>k</mi>
<mo>,</mo>
<mi>i</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>x</mi>
<mrow>
<mi>t</mi>
<mi>k</mi>
<mo>,</mo>
<mi>i</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msub>
<mi>&lambda;</mi>
<mrow>
<mi>k</mi>
<mo>,</mo>
<mi>i</mi>
</mrow>
</msub>
<msub>
<mi>r</mi>
<mrow>
<mi>k</mi>
<mo>,</mo>
<mi>i</mi>
</mrow>
</msub>
</mrow>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>8</mn>
<mo>)</mo>
</mrow>
</mrow>
WhereinD is the distance between two passages, vsk,iFor horizontal velocity component, vsi,iFor vertical velocity component；
D) calculating of signal time delay
According to each target reach satellite distance, can calculate obtain each time slot signal signal time delay be
<mrow>
<msub>
<mi>&tau;</mi>
<mrow>
<mi>k</mi>
<mo>,</mo>
<mi>i</mi>
</mrow>
</msub>
<mo>=</mo>
<mfrac>
<msub>
<mi>r</mi>
<mrow>
<mi>k</mi>
<mo>,</mo>
<mi>i</mi>
</mrow>
</msub>
<msup>
<mi>c</mi>
<mo>,</mo>
</msup>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>9</mn>
<mo>)</mo>
</mrow>
</mrow>
Wherein c ' is the light velocity；
The parameter that each time slot is calculated in (1) formula can be generated according to (2)~(9) formula, many mesh can be calculated by (1) formula
The signal waveform data of superposition is marked, signal waveform data binary file is written as to disk file, you can obtain baseband signal
Data file.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710148410.4A CN107066693B (en) | 2017-03-14 | 2017-03-14 | Multichannel multi-target satellite-borne AIS reconnaissance signal simulation system |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710148410.4A CN107066693B (en) | 2017-03-14 | 2017-03-14 | Multichannel multi-target satellite-borne AIS reconnaissance signal simulation system |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN107066693A true CN107066693A (en) | 2017-08-18 |
| CN107066693B CN107066693B (en) | 2020-06-02 |
Family
ID=59622971
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201710148410.4A Expired - Fee Related CN107066693B (en) | 2017-03-14 | 2017-03-14 | Multichannel multi-target satellite-borne AIS reconnaissance signal simulation system |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN107066693B (en) |
Cited By (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108055068A (en) * | 2017-12-08 | 2018-05-18 | 上海埃威航空电子有限公司 | A kind of airborne simulation system and test method for receiving marine ais signal |
| CN109557516A (en) * | 2018-11-28 | 2019-04-02 | 四川九洲电器集团有限责任公司 | The rapid generation of multiple target spacing wave |
| CN110031814A (en) * | 2017-08-31 | 2019-07-19 | 成都玖锦科技有限公司 | A kind of continuous Multiple Target Signals synthetic method of frequency spectrum |
| CN110138416A (en) * | 2019-04-11 | 2019-08-16 | 上海卫星工程研究所 | The spaceborne wired multi-channel detection probabilistic testing method of AIS ship oceangoing ship |
| CN110224742A (en) * | 2019-06-24 | 2019-09-10 | 电子科技大学 | A kind of spaceborne ADS-B signal forming method for testing Signal separator performance |
| CN111444609A (en) * | 2020-03-24 | 2020-07-24 | 北京润科通用技术有限公司 | Data processing method and simulation system |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2012050438A1 (en) * | 2010-10-15 | 2012-04-19 | Isis - Innovative Solutions In Space B.V. | Satellite communication system for retrieving a plurality of modulated signals having a constant modulus |
| CN103117965A (en) * | 2013-01-21 | 2013-05-22 | 天津理工大学 | Joint estimation method of timing frequency offset of satellite borne automatic identification system (AIS) signals and implementation system thereof |
| CN105610761A (en) * | 2015-12-16 | 2016-05-25 | 西安空间无线电技术研究所 | Application layer system-level constraint-based satellite borne GMSK bit error rate improvement system |
| CN105634581A (en) * | 2014-10-31 | 2016-06-01 | 航天恒星科技有限公司 | Method and system for reducing AIS signal aliasing |
| CN106453178A (en) * | 2016-11-06 | 2017-02-22 | 中国电子科技集团公司第十研究所 | Satellite-based AIS signal intercept and demodulation method |
-
2017
- 2017-03-14 CN CN201710148410.4A patent/CN107066693B/en not_active Expired - Fee Related
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2012050438A1 (en) * | 2010-10-15 | 2012-04-19 | Isis - Innovative Solutions In Space B.V. | Satellite communication system for retrieving a plurality of modulated signals having a constant modulus |
| CN103117965A (en) * | 2013-01-21 | 2013-05-22 | 天津理工大学 | Joint estimation method of timing frequency offset of satellite borne automatic identification system (AIS) signals and implementation system thereof |
| CN105634581A (en) * | 2014-10-31 | 2016-06-01 | 航天恒星科技有限公司 | Method and system for reducing AIS signal aliasing |
| CN105610761A (en) * | 2015-12-16 | 2016-05-25 | 西安空间无线电技术研究所 | Application layer system-level constraint-based satellite borne GMSK bit error rate improvement system |
| CN106453178A (en) * | 2016-11-06 | 2017-02-22 | 中国电子科技集团公司第十研究所 | Satellite-based AIS signal intercept and demodulation method |
Non-Patent Citations (1)
| Title |
|---|
| 栗霖雲等: "星载AIS系统信号模拟源的仿真实现", 《无线电工程》 * |
Cited By (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110031814A (en) * | 2017-08-31 | 2019-07-19 | 成都玖锦科技有限公司 | A kind of continuous Multiple Target Signals synthetic method of frequency spectrum |
| CN110031814B (en) * | 2017-08-31 | 2023-03-24 | 成都玖锦科技有限公司 | Frequency spectrum continuous multi-target signal synthesis method |
| CN108055068A (en) * | 2017-12-08 | 2018-05-18 | 上海埃威航空电子有限公司 | A kind of airborne simulation system and test method for receiving marine ais signal |
| CN108055068B (en) * | 2017-12-08 | 2020-08-25 | 上海埃威航空电子有限公司 | Test method of airborne analog simulation system for receiving offshore AIS signals |
| CN109557516A (en) * | 2018-11-28 | 2019-04-02 | 四川九洲电器集团有限责任公司 | The rapid generation of multiple target spacing wave |
| CN110138416A (en) * | 2019-04-11 | 2019-08-16 | 上海卫星工程研究所 | The spaceborne wired multi-channel detection probabilistic testing method of AIS ship oceangoing ship |
| CN110224742A (en) * | 2019-06-24 | 2019-09-10 | 电子科技大学 | A kind of spaceborne ADS-B signal forming method for testing Signal separator performance |
| CN111444609A (en) * | 2020-03-24 | 2020-07-24 | 北京润科通用技术有限公司 | Data processing method and simulation system |
| CN111444609B (en) * | 2020-03-24 | 2023-08-08 | 北京润科通用技术有限公司 | Data processing method and simulation system |
Also Published As
| Publication number | Publication date |
|---|---|
| CN107066693B (en) | 2020-06-02 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN107066693A (en) | The spaceborne AIS reconnaissance signals simulation system of multi-channel multi-target | |
| CN106483512B (en) | A kind of general multichannel distributed object analogue echoes method and if system | |
| US7145504B1 (en) | Arbitrary radar target synthesizer (arts) | |
| CN106501782B (en) | A kind of multichannel satellite navigation and interference signal simulation source and signal imitation method | |
| EP2005210B1 (en) | Methods and systems for generating virtual radar targets | |
| CN102435984B (en) | Echo simulation device for Doppler radar and implementation method thereof | |
| CN109782241B (en) | Meteorological radar echo simulation method and system | |
| CN105333887B (en) | Radio altimeter signal imitation system | |
| CN102176029B (en) | Global positioning system (GPS) direct and multipath signal simulator and simulation method | |
| CN110515050B (en) | Satellite-borne SAR real-time echo simulator based on GPU | |
| CN107367717B (en) | Radar multi-target intermediate frequency simulation device | |
| JPH07311257A (en) | Radar target wave simulator | |
| CN113697126A (en) | Unmanned aerial vehicle anti-interference performance evaluation system and method under complex electromagnetic environment | |
| CN202794515U (en) | Navigation radar for FMCW (frequency-modulated continuous wave) ship | |
| KR101534167B1 (en) | Apparatus for analysing real time jamming effectiveness of Satellite Navigation | |
| CN115586499A (en) | Radar signal environment simulation method and system | |
| KR102172379B1 (en) | Method and apparatus for generating simulated signals in a multistatic passive coherent location(PCL) system | |
| JP2010230549A (en) | Generator for simulated clutter | |
| CN103727960B (en) | A kind of radio altimeter interference signal production method based on DRFM | |
| CN108983240A (en) | Anticollision millimetre-wave radar echo signal simulation system and method based on orthogonal modulation system | |
| US20220099798A1 (en) | Real-time closed-loop digital radar simulator | |
| CN109471110B (en) | Multi-target signal simulation method and device for millimeter wave radar for vehicle | |
| CN202093173U (en) | Multi-system navigation satellite signal simulator | |
| CN115308703B (en) | High-frequency MIMO radar auxiliary calibration source and calibration method | |
| US6128589A (en) | Method and apparatus for modelling a system which includes the transmission and reception of signals |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant | ||
| TR01 | Transfer of patent right |
Effective date of registration: 20220225 Address after: 410073 National University of Defense Technology, Hunan, Changsha Patentee after: National University of Defense Technology Address before: 410000 room 1802, financial block, Hunan wealth center, Furong Middle Road, Changsha, Hunan Patentee before: CHANGSHA DINGDUAN ELECTRONIC TECHNOLOGY CO.,LTD. |
|
| TR01 | Transfer of patent right | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20200602 |
|
| CF01 | Termination of patent right due to non-payment of annual fee |