CN107272033B - A kind of random delay modeling method of the airborne reflection signal of Beidou MEO satellite - Google Patents
A kind of random delay modeling method of the airborne reflection signal of Beidou MEO satellite Download PDFInfo
- Publication number
- CN107272033B CN107272033B CN201710392310.6A CN201710392310A CN107272033B CN 107272033 B CN107272033 B CN 107272033B CN 201710392310 A CN201710392310 A CN 201710392310A CN 107272033 B CN107272033 B CN 107272033B
- Authority
- CN
- China
- Prior art keywords
- signal
- random
- reflection
- reflection signal
- delay
- 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
Links
- 238000000034 method Methods 0.000 title claims abstract description 26
- 239000011159 matrix material Substances 0.000 claims abstract description 19
- 230000008569 process Effects 0.000 claims abstract description 7
- 238000009826 distribution Methods 0.000 claims description 15
- 238000011161 development Methods 0.000 abstract description 4
- 230000007246 mechanism Effects 0.000 abstract description 3
- 238000011160 research Methods 0.000 abstract description 3
- 238000004088 simulation Methods 0.000 abstract description 3
- 238000010586 diagram Methods 0.000 description 4
- 230000008901 benefit Effects 0.000 description 3
- 238000012935 Averaging Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000013178 mathematical model Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000000149 penetrating effect Effects 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 230000002035 prolonged effect Effects 0.000 description 1
Classifications
-
- 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/01—Satellite radio beacon positioning systems transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
- G01S19/13—Receivers
- G01S19/35—Constructional details or hardware or software details of the signal processing chain
- G01S19/37—Hardware or software details of the signal processing chain
Landscapes
- Engineering & Computer Science (AREA)
- Signal Processing (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Computer Networks & Wireless Communication (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Radar Systems Or Details Thereof (AREA)
- Radio Relay Systems (AREA)
Abstract
The present invention provides a kind of random delay modeling methods of the airborne reflection signal of Beidou MEO satellite, belong to modeling operation field, selection including random delay, the reflection signal of multiple groups random delay switches, the selection of random amplitude, and reflection signal combining output, original fixed delay model is replaced by using random delay model, range is established by what the statistical property for reflecting signal and dynamic characteristic were included in model, reflection signal matrix is established in modeling process simultaneously, and introduce the handover mechanism of matrix, so that the signal generated has more randomness, to improve the validity of output signal, it is provided safeguard for the research of the airborne reflection analog simulation of Beidou MEO satellite and the development of simulator.
Description
Technical field
The invention belongs to model operation field, the in particular to random delay of a kind of airborne reflection signal of Beidou MEO satellite
Modeling method.
Background technique
With the development of technology, satellite navigation has incorporated the every aspect of work and life, is the trip and work of the people
Provide guarantee.The especially Beidou Navigation System of China's voluntarily family, shows more excellent in terms of navigation accuracy
Characteristic.
But at present in terms of the signal processing of navigation satellite, it would be highly desirable to a kind of better modeling pattern is established, thus raising pair
The reliability of satellite-signal.
Summary of the invention
In order to solve shortcoming and defect existing in the prior art, the present invention provides defend for improving a kind of Beidou MEO
The random delay modeling method of the airborne reflection signal of star.
In order to reach above-mentioned technical purpose, the present invention provides a kind of Beidou MEO satellite it is airborne reflection signal it is random when
Prolong modeling method, the modeling method, comprising:
Step 1, the selection of random delay, comprising:
Determine specular reflection point, that chooses sea echo relative specular reflection point time delay is distributed as [0, K] chip, with
Every reflection signal randomly selects in the section in machine Time Delay Model；
It is assumed that reflection signal strips number is N, the time delay central point of N item reflection signal can be expressed as in random delay model
In formula (1): K is the reflected signal delay distribution upper limit；N is reflection signal strips number；I is reflection signal serial number；
TiIt is the time delay central point of i-th reflection signal；
In view of the particularity of specular reflection point, the time delay value of signal is not involved at random at specular reflection point, except mirror surface is anti-
Other (N-1) items reflection signal at exit point is in respective time delay central pointA bit, then N item reflects for random selection in chip
Signal random delay changing value Δ tiAnd random delay value tiIt is respectively as follows:
ti=Ti+Δti, i=0,1 ..., N-1 (3)
In formula (2), (3): rand () function is the equally distributed random number of 0-1；ΔtiWhen random for i-th reflection signal
Prolong changing value；tiFor i-th reflection signal random delay value；
Utilize tiRandom reflected signal time-domain position is fitted when as fitting, obtains the anti-of N random delay distribution
Penetrate signal；
Step 2, the reflection signal switching of multiple groups random delay, comprising:
It repeats above method M times, obtains M group, every group of N random delay reflection signal, define M'=M+O, O is 0 signal
Group group number, defines random delay model and reflects signal matrix SM'NAre as follows:
In formula (4), M is that random delay reflects signal group number；N is the item number of the reflection signal in every group；M' is that joined O
Reflection signal group number after 0 signal of group；SM'NOne group of random delay of each behavior reflects signal in matrix, each to be classified as same time delay
Center point reflection signal；sijSignal is reflected for i-th group of reflection signal group j-th strip；Wherein 0 signal then represents corresponding time delay central point
Channel interruption of the neighbouring signal to reflection signal receiver；
Channel interruption characteristic caused by 0 signal imitation sea dynamic of O group is introduced, the signal at every group of same delay center is (same
One column signal) introduce switching frequency fi, i=1,2 ..., N are randomly selected from M bars when this group of switching signal arrives
Z item output, for realizing the output of the reflection signal of random delay；
Step 3, the selection of random amplitude, comprising:
The Z*N as shown in formula (5) mutually independent random value A are introduced, A obeys exponential distribution (mean value 1):
A~E (λ), λ=1 (5)
The reflection signal of each output will be multiplied by the value, and the signal amplitude is with switching frequency (fi, i=1,
2 ..., N) it is randomly selected；
Step 4, reflection signal combining output, comprising:
The time delay of the reflection signal of any time output is random, amplitude is random, then letter is reflected in the combining of any time output
Number indicate are as follows:
In formula (6): R_Signals is the output signal of random delay model；Randint () is random whole to be uniformly distributed
Number function；ixThe random integers for arriving M' for 1；For SM'NI in matrixxThe reflection signal of row j column position；Ax,jTo obey mean value
For 1 mutually independent exponential distribution；Z is the reflection signal strips number that each time delay center position randomly selects output；
The time delay for reflecting signal in output signal R_Signals hasKind random combine, random magnitude are bigger
The randomness of signal is increased greatly；SM'NActually random delay and random power module reflect signal library, row (M) and column (N)
Unlimited extension is carried out to cover all possible reflection signal, is then taken out at random from the reflection signal matrix being randomly derived again
The number of winning the confidence is exported multiplied by combining after random magnitude；
Wherein, the value range of K, N are natural number.
Technical solution provided by the invention has the benefit that
Original fixed delay model is replaced by using random delay model, the statistical property and dynamic of signal will be reflected
Characteristic is included in the range of establishing of model, while reflection signal matrix is established in modeling process, and introduce the switching machine of matrix
System, so that the signal generated has more randomness, it is that Beidou MEO satellite is airborne anti-to improve the validity of output signal
The development of the research and simulator of penetrating analog simulation provides safeguard.
Detailed description of the invention
It, below will be to attached drawing needed in embodiment description in order to illustrate more clearly of technical solution of the present invention
It is briefly described, it should be apparent that, drawings in the following description are only some embodiments of the invention, general for this field
For logical technical staff, without creative efforts, it is also possible to obtain other drawings based on these drawings.
Fig. 1 is a kind of process of the random delay modeling method of the airborne reflection signal of Beidou MEO satellite provided by the invention
Schematic diagram；
Fig. 2 is the random range schematic diagram of reflected signal delay provided by the invention；
Fig. 3 is random delay model reflection signal distributions schematic diagram provided by the invention.
Specific embodiment
To keep structure and advantage of the invention clearer, structure of the invention is made further below in conjunction with attached drawing
Description.
Embodiment one
The present invention provides a kind of random delay modeling method of the airborne reflection signal of Beidou MEO satellite, the modeling sides
Method, as shown in Figure 1, comprising:
Step 1, the selection of random delay, comprising:
Determine specular reflection point, that chooses sea echo relative specular reflection point time delay is distributed as [0, K] chip, with
Every reflection signal randomly selects in the section in machine Time Delay Model；
It is assumed that reflection signal strips number is N, the time delay central point of N item reflection signal can be expressed as in random delay model
In formula (1): K is the reflected signal delay distribution upper limit；N is reflection signal strips number；I is reflection signal serial number；
TiIt is the time delay central point of i-th reflection signal；The random range of time delay for delimiting every reflection signal is as shown in Figure 2.
In view of the particularity of specular reflection point, the time delay value of signal is not involved at random at specular reflection point, except mirror surface is anti-
Other (N-1) items reflection signal at exit point is in respective time delay central pointA bit, then N item reflects for random selection in chip
Signal random delay changing value Δ tiAnd random delay value tiIt is respectively as follows:
ti=Ti+Δti, i=0,1 ..., N-1 (3)
In formula (2), (3): rand () function is the equally distributed random number of 0-1；ΔtiWhen random for i-th reflection signal
Prolong changing value；tiFor i-th reflection signal random delay value；
Utilize tiRandom reflected signal time-domain position is fitted when as fitting, obtains the anti-of N random delay distribution
Signal is penetrated, time delay distribution schematic diagram is as shown in Figure 3.
Step 2, the reflection signal switching of multiple groups random delay, comprising:
It repeats above method M times, obtains M group, every group of N random delay reflection signal, define M'=M+O, O is 0 signal
Group group number, defines random delay model and reflects signal matrix SM'NAre as follows:
In formula (4), M is that random delay reflects signal group number；N is the item number of the reflection signal in every group；M' is that joined O
Reflection signal group number after 0 signal of group；SM'NOne group of random delay of each behavior reflects signal in matrix, each to be classified as same time delay
Center point reflection signal；sijSignal is reflected for i-th group of reflection signal group j-th strip；Wherein 0 signal then represents corresponding time delay central point
Channel interruption of the neighbouring signal to reflection signal receiver；
Channel interruption characteristic caused by 0 signal imitation sea dynamic of O group is introduced, the signal at every group of same delay center is (same
One column signal) introduce switching frequency fi, i=1,2 ..., N are selected from M' bars at random when this group of switching signal arrives
Z item is taken to export, for realizing the output of the reflection signal of random delay；
Step 3, the selection of random amplitude, comprising:
The Z*N as shown in formula (5) mutually independent random value A are introduced, A obeys exponential distribution (mean value 1):
A~E (λ), λ=1 (5)
The reflection signal of each output will be multiplied by the value, and the signal amplitude is with switching frequency (fi, i=1,
2 ..., N) it is randomly selected；
Step 4, reflection signal combining output, comprising:
The time delay of the reflection signal of any time output is random, amplitude is random, then letter is reflected in the combining of any time output
Number indicate are as follows:
In formula (6): R_Signals is the output signal of random delay model；Randint () is random whole to be uniformly distributed
Number function；ixThe random integers for arriving M' for 1；For SM'NI in matrixxThe reflection signal of row j column position；Ax,jTo obey mean value
For 1 mutually independent exponential distribution；Z is the reflection signal strips number that each time delay center position randomly selects output；
The time delay for reflecting signal in output signal R_Signals hasKind random combine, random magnitude is more significantly
Increase the randomness of signal；SM'NActually random delay and random power module reflect signal library, row (M) and arrange (N) can be into
Then row unlimited extension is randomly selected from the reflection signal matrix being randomly derived again with covering all possible reflection signal
Signal is exported multiplied by combining after random magnitude；
Wherein, the value range of K, N are natural number.
The technical scheme is that
(1) it according to reflection signal inverse analog basic principle, randomly selects time delay and carries out reflection Signal averaging.Delimit every
The random range of time delay for reflecting signal, obtains dynamic random delay.
(2) above procedure is repeated, obtains multiple groups random delay reflection signal, and the reflection signal group that 0 signal is added carrys out mould
Channel interruption characteristic caused by quasi- sea dynamic introduces the output that switching frequency carries out random delay reflection signal.
(3) random value for obeying exponential distribution is introduced, to control the amplitude of output reflection signal.
(4) by the above process, the time delay and amplitude of the reflection signal of any time output are random, therefore are combined
Output reflection signal.
The invention has the advantages that
(1) this method is improved to random delay model on the basis of fixed delay model, to not only allow for reflecting
The statistical property of signal, it is also contemplated that the dynamic characteristic of reflection signal.
(2) reflection signal matrix is established, and introduces handover mechanism, so that the composition of output signal is more complicated, more at random
Property.
(3) this method is the proposition based on mathematical model, is applicable to various programming languages and completes software and hardware realization.
The present invention provides a kind of Beidou MEO satellite it is airborne reflection signal random delay modeling method, including it is random when
The switching of reflection signal, the selection of random amplitude and the reflection signal combining output of the selection, multiple groups random delay prolonged, pass through
Original fixed delay model is replaced using random delay model, the statistical property for reflecting signal and dynamic characteristic are included in model
Establish range, while reflection signal matrix is established in modeling process, and introduce the handover mechanism of matrix, so that generating
Signal has more randomness, is the airborne reflection analog simulation of Beidou MEO satellite to improve the validity of output signal
Research and the development of simulator provide safeguard.
Each serial number in above-described embodiment is for illustration only, the assembling for not representing each component or the elder generation in use process
Sequence afterwards.
The above description is only an embodiment of the present invention, is not intended to limit the invention, all in the spirit and principles in the present invention
Within, any modification, equivalent replacement, improvement and so on should all be included in the protection scope of the present invention.
Claims (1)
1. a kind of random delay modeling method of the airborne reflection signal of Beidou MEO satellite, which is characterized in that the modeling method,
Include:
Step 1, the selection of random delay, comprising:
Determine specular reflection point, that chooses sea echo relative specular reflection point time delay is distributed as [0, K] chip, when random
Prolong every reflection signal in model to randomly select in the section；
It is assumed that reflection signal strips number is N, the time delay central point of N item reflection signal can be expressed as in random delay model
In formula (1): K is the reflected signal delay distribution upper limit；N is reflection signal strips number；I is reflection signal serial number；TiIt is
The time delay central point of i item reflection signal；
In view of the particularity of specular reflection point, the time delay value of signal is not involved in the process of randomly selecting at specular reflection point, removes mirror
Other (N-1) items reflection signal at the reflection point of face is in respective time delay central pointIt is randomly choosed in chip a bit, then N item
Reflect signal random delay changing value △ tiAnd random delay value tiIt is respectively as follows:
ti=Ti+△ti, i=0,1 ..., N-1 (3)
In formula (2), (3): rand () function is the equally distributed random number of 0-1；△tiBecome for i-th reflection signal random delay
Change value；tiFor i-th reflection signal random delay value；
Utilize tiRandom reflected signal time-domain position is fitted when as fitting, obtains the reflection letter of N random delay distribution
Number；
Step 2, the reflection signal switching of multiple groups random delay, comprising:
It repeating step 1 M times, obtains M group, every group of N random delay reflection signal, define M'=M+O, O is 0 signal group group number,
It defines random delay model and reflects signal matrix SM'NAre as follows:
In formula (4), M is that random delay reflects signal group number；N is the item number of the reflection signal in every group；M' is that joined O group 0
Reflection signal group number after signal；SM'NOne group of random delay of each behavior reflects signal in matrix, each to be classified as while the Yanzhong heart
Point reflection signal；sijSignal is reflected for i-th group of reflection signal group j-th strip；Wherein 0 signal then represents near corresponding time delay central point
Channel interruption of the signal to reflection signal receiver；
Introduce channel interruption characteristic caused by 0 signal imitation sea dynamic of O group, the signal, that is, same row at every group of same delay center
Signal introduces switching frequency fi, i=1,2 ..., N randomly select Z item when this group of switching signal arrives from M' bars
Output, for realizing the output of the reflection signal of random delay；
Step 3, the selection of random amplitude, comprising:
The Z*N as shown in formula (5) mutually independent random value A are introduced, A obeys exponential distribution, mean value 1:
A~E (λ), λ=1 (5)
The reflection signal of each output will be multiplied by the value, and the signal amplitude is with switching frequency fi, i=1,2 ..., N is carried out
It randomly selects；
Step 4, reflection signal combining output, comprising:
The time delay of the reflection signal of any time output is random, amplitude is random, then signal table is reflected in the combining of any time output
It is shown as:
In formula (6): R_Signals is the output signal of random delay model；Randint () is to be uniformly distributed random integers letter
Number；ixThe random integers for arriving M' for 1；For SM'NI in matrixxThe reflection signal of row j column position；Ax,jIt is 1 to obey mean value
Mutually independent exponential distribution；Z is the reflection signal strips number that each time delay center position randomly selects output；
The time delay for reflecting signal in output signal R_Signals hasKind random combine, random magnitude increase more greatly
The randomness of signal is added；SM'NActually random delay and random power module reflect signal library, and it is unlimited that M row and N column can carry out
Then system extension is randomly selected signal from the reflection signal matrix being randomly derived again and is multiplied to cover all possible reflection signal
To be combined output after random magnitude；
Wherein, the value range of K, N are natural number.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710392310.6A CN107272033B (en) | 2017-05-27 | 2017-05-27 | A kind of random delay modeling method of the airborne reflection signal of Beidou MEO satellite |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710392310.6A CN107272033B (en) | 2017-05-27 | 2017-05-27 | A kind of random delay modeling method of the airborne reflection signal of Beidou MEO satellite |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN107272033A CN107272033A (en) | 2017-10-20 |
| CN107272033B true CN107272033B (en) | 2019-07-02 |
Family
ID=60065770
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201710392310.6A Active CN107272033B (en) | 2017-05-27 | 2017-05-27 | A kind of random delay modeling method of the airborne reflection signal of Beidou MEO satellite |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN107272033B (en) |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109856646B (en) * | 2018-12-24 | 2021-03-26 | 西安空间无线电技术研究所 | Method for controlling relative time delay between high-precision satellite navigation signals |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5546087A (en) * | 1993-11-05 | 1996-08-13 | Agence Spatiale Europeene | Altimetry method |
| CN2700895Y (en) * | 2004-04-29 | 2005-05-18 | 杨东凯 | Delay mapping receiver based on global satellite positioning system |
| CN101349749A (en) * | 2008-09-10 | 2009-01-21 | 北京航空航天大学 | Generating device for simulating navigation satellite reflection signal |
| CN105182308A (en) * | 2015-08-27 | 2015-12-23 | 北京航空航天大学 | On-board GNSS marine reflection signal generation method |
-
2017
- 2017-05-27 CN CN201710392310.6A patent/CN107272033B/en active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5546087A (en) * | 1993-11-05 | 1996-08-13 | Agence Spatiale Europeene | Altimetry method |
| CN2700895Y (en) * | 2004-04-29 | 2005-05-18 | 杨东凯 | Delay mapping receiver based on global satellite positioning system |
| CN101349749A (en) * | 2008-09-10 | 2009-01-21 | 北京航空航天大学 | Generating device for simulating navigation satellite reflection signal |
| CN105182308A (en) * | 2015-08-27 | 2015-12-23 | 北京航空航天大学 | On-board GNSS marine reflection signal generation method |
Non-Patent Citations (2)
| Title |
|---|
| A simulator for GNSS-R polarimetric observations over the ocean;Domenico Schiavulli 等;《IEEE》;20141231;3802-3805 |
| 机载GNSS海洋反射信号的建模与仿真;祁永强 等;《北京航空航天大学学报》;20170331;第43卷(第3期);567-572 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN107272033A (en) | 2017-10-20 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Athens et al. | A paleoenvironmental and archaeological model-based age estimate for the colonization of Hawai’l | |
| Shashkov | Closure models for multimaterial cells in arbitrary Lagrangian–Eulerian hydrocodes | |
| CN107589410B (en) | A method for multi-target signal synthesis without breakpoints | |
| GB2384558A (en) | Neural net prediction of seismic streamer shape | |
| CN105553447B (en) | Clock switch circuit | |
| CN107272033B (en) | A kind of random delay modeling method of the airborne reflection signal of Beidou MEO satellite | |
| CN106291011B (en) | A kind of random pulse generator based on FPGA | |
| Trindade et al. | Observed and dark diversity dynamics over millennial time scales: fast life-history traits linked to expansion lags of plants in northern Europe | |
| JPS5386136A (en) | Control circuit | |
| Harper | The Ordovician brachiopod radiation: roles of alpha, beta, and gamma diversity | |
| US6028988A (en) | System for logic synthesis-for-testability capable of improving testability for an FSM having an asynchronous reset state | |
| CN212255684U (en) | Extensible DBF combined system | |
| CN109344328A (en) | Obtain the method and device of recommender system best parameter group | |
| Ryabov et al. | Depressed emission between magnetic arcades near a sunspot | |
| Thiery et al. | Universality in the mean spatial shape of avalanches | |
| CN114706067B (en) | Target detection method, device, equipment and computer readable storage medium | |
| Darbo et al. | A three-dimensional ripple trigger for the DELPHI Time Projection Chamber | |
| Cho et al. | Boosting invisible Higgs boson searches by tagging a gluon jet for the gluon fusion process | |
| Ghalwash | Extended Petri-Net theory for modeling and analysis of decision-making organizations. | |
| Brau et al. | Higgs branching ratio measurements at a future Linear Collider | |
| Sjöstrand | Monte Carlo event generation for LHC | |
| Altaber et al. | Readout system test benches | |
| SU1262518A1 (en) | Device for investigating graphs | |
| Wang et al. | Inversion of the latitudinal diversity gradient at high taxonomic level in liverworts revealed by a phylogenetically deconstructive approach | |
| SU1228112A1 (en) | Device for studying paths in graphs |
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 |