CN107894231A - A kind of X-ray pulsar discrimination method based on Hilbert transform - Google Patents
A kind of X-ray pulsar discrimination method based on Hilbert transform Download PDFInfo
- Publication number
- CN107894231A CN107894231A CN201711082587.5A CN201711082587A CN107894231A CN 107894231 A CN107894231 A CN 107894231A CN 201711082587 A CN201711082587 A CN 201711082587A CN 107894231 A CN107894231 A CN 107894231A
- Authority
- CN
- China
- Prior art keywords
- signal
- ray pulsar
- ray
- pulsar
- contour
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01C—MEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
- G01C21/00—Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00
- G01C21/02—Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00 by astronomical means
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Astronomy & Astrophysics (AREA)
- Automation & Control Theory (AREA)
- General Physics & Mathematics (AREA)
- Image Analysis (AREA)
Abstract
A kind of X-ray pulsar discrimination method based on Hilbert transform, it is related to the profile discrimination method in Weak Signal Processing field.It solves the problems, such as that current identification takes a large amount of operation times and memory space, while can effectively suppress the influence of high-frequency noise.The present invention step be：First, Hilbert transform is done to nominal contour and observation contour signal.2nd, the instantaneous amplitude feature of extraction standard contour signal and contour signal to be identified, and utilize the characteristic vector construction feature database of nominal contour.3rd, respectively in the characteristic vector of calculating observation contour signal and property data base each characteristic vector space length, when distance obtains minimum value, corresponding pulsar sequence number is identification result.The present invention utilizes the rejection characteristic of Hilbert transform pairs high-frequency noise, and extraction instantaneous amplitude characteristic Design minimum distance classifier completes the high speed identification of pulsar signal, suitable for Fast Identification X-ray pulsar.
Description
Technical field
The present invention relates to the profile discrimination method in Weak Signal Processing field, and in particular to one kind is based on Hilbert transform
X-ray pulsar discrimination method.
Background technology
X-ray pulsar is that radiated photons energy concentrates on X ray frequency range, the rotation period stablizes, be remote apart from the solar system
A kind of neutron star.X-ray pulsar navigation is not by artificial disturbance, it is possible to achieve entirely autonomous navigation, while X-ray pulsar
Navigation can solve the navigation problem of large-scale space-time.
X-ray pulsar is remote apart from the solar system, and the photon for traveling to spacecraft is easily flooded by cosmic background noise, energy
The photon signal enough received is very faint, and this just needs detector to be directed at the direction of X-ray pulsar radiation beam for a long time, right
The photon received is folded according to some cycles can just obtain observing profile.Complete the profile of X-ray pulsar signal
Handled after extraction, it is necessary to observe X-ray pulsar and could carry out follow-up navigation calculating after contour signal correctly recognizes.
The identification problem conventional method for X-ray pulsar observation contour signal is to utilize the amplitude for observing profile at present
Information is had a great influence compared with nominal contour by noise, and utilize observe contour signal extraction bispectrum feature or
The method of high-order statistic is, it is necessary to substantial amounts of operation time and memory space.
The content of the invention
When taking a large amount of computings when carrying out X-ray pulsar identification using high-order statistic it is an object of the invention to overcome
Between and memory space deficiency, there is provided it is a kind of based on Hilbert transform extraction instantaneous amplitude feature method to X-ray pulse
Star is recognized.It can play the inhibitory action of Hilbert transform pairs high-frequency noise, at the same take seldom calculating time and
Memory space, so as to the identification suitable for X-ray pulsar.
The purpose of the present invention is achieved through the following technical solutions：The instantaneous amplitude for extracting contour signal is special as signal
Sign, after obtaining nominal contour property data base and observing the characteristic vector of contour signal, according to identification criterion design most narrow spacing
From grader, the space length between characteristic vector and each nominal contour feature is calculated respectively, is realized to different x-ray pulse
Star is classified.
The flow chart of the present invention is as shown in figure 1, comprise the following steps that：
Step 1：Hilbert transform.
X-ray pulsar signal is that the photon number received by one section of observation time internal X-ray detector is accumulated and formed
Observation contour signal form, for different X-ray pulsar signals, because X-ray pulsar has specific thing
Feature and space characteristics are managed, traveling to interplanetary X-ray pulsar signal has geometry uniqueness.Therefore, X is selected
Feature of the instantaneous amplitude feature that ray pulse star signal obtains after Hilbert transform needed for as identification.Assuming that X is penetrated
Line pulsar signal is x (t), and the signal y (t) after Hilbert transform can be expressed asThis is
X (t) and 1/ π τ convolution.
Step 2：Characteristics extraction.
Converted by Hilbert, the real part x (t) and imaginary part y (t) of X-ray pulsar signal constitute a conjugation
Relation it is plural right, signal resolution form is expressed as S (t)=x (t)+jy (t), the instantaneous amplitude feature of X-ray pulsar signal
Function A (t) can be expressed asInstantaneous amplitude can effectively reflect that X-ray pulsar is believed
Number-detailed information of energy.Therefore, the nominal contour feature set F in navigational route databasesF can be expressed ass=[v1(t),v2
(t),...,vm(t)]T, wherein m represent navigational route database in X-ray pulsar signal quantity.Assuming that i-th of X-ray pulse
The characteristic vector of star signal is vi(t), then there is vi(t)=[xi(t),yi(t)]T.The X ray arteries and veins in navigational route database is extracted respectively
Rush the quasi- contour signal of asterisk and observe the instantaneous amplitude feature of contour signal and establish X-ray pulsar signal standards contour feature
Collect Fs, correspondingly the instantaneous amplitude-frequency characteristic vector of X-ray pulsar observation contour signal is Vs(t), then there is Vs(t)=[xs(t),
ys(t)]T。
Step 3：Design minimum distance classifier.
, can be with according to identification criterion after having obtained X-ray pulsar nominal contour and having observed the characteristic vector of contour signal
Minimum distance classifier is designed, calculates characteristic vector V respectivelys(t) with nominal contour feature set FsBetween interior each characteristic vector
Space length.Assuming that there is m X-ray pulsar in navigational route database, then the standard of i-th X-ray pulsar (1≤i≤m)
Outline Feature Vector viAnd V (t)s(t) distance can be expressed asCause
This, can obtain identification space length collection and be combined into：D={ D1,D2,...,Dm, work as DiWhen obtaining minimum value, corresponding X ray
Pulsar sequence number is identification result, i.e.,：Wherein L is the identification sequence number of X-ray pulsar.
The present invention has the following advantages that compared with prior art：
The present invention proposes a kind of X-ray pulsar discrimination method based on instantaneous amplitude feature, and based on bispectrum feature
Discrimination method compared to taking less calculating space and memory space, while the present invention utilizes Hilbert transform extraction observation
The instantaneous amplitude of profile, high-frequency noise can be suppressed.Compared to the situation for only considering observation profile amplitude information, this programme utilizes
The real and imaginary parts information structuring characteristic vector of signal, can reflect the detailed information of X-ray pulsar signal-energy, reach
To effective identification effect.X-ray pulsar signal is recognized using this programme, can improve comprehensively identification accuracy and
Speed.
Brief description of the drawings
Fig. 1 is the flow chart of the present invention.
Fig. 2 is X-ray pulsar identification result.
Fig. 3 is to contrast the X-ray pulsar J1643-1224 signal recognition times.
Fig. 4 is to contrast the X-ray pulsar B1257+12 signal recognition times.
Embodiment
With reference to embodiment and the embodiment of the brief description of the drawings present invention.
Perform step 1：X-ray pulsar nominal contour data first by EPN 11 observations provided include navigation number
According to storehouse.Therefore, the collection for observing X-ray pulsar is combined into：B1257+12, B1953+29, J0621+1002, J1012+5307,
J1518+4904, J1640+2224, J1643-1224, J1713+0747, J1730-2304, J1744-24A, J2051-0827 },
The sample of the set is attributed to the class that numbering is { 1 ..., 11 } respectively.
Perform step 2：The instantaneous amplitude characteristic vector of extraction observation contour signal, is calculated using minimum distance classifier
Its distance with each nominal contour characteristic vector, obtains recognizing the set of space length.
Perform step 3：The X-ray pulsar sequence number corresponding to minimum value in selection set is as identification result.
The X-ray pulsar of table 1 observes contour signal and the result of calculation of nominal contour signal characteristic distance
Identification result of the X-ray pulsar of table 1 in navigational route database
As it can be seen from table 1 for 11 class X-ray pulsar contour signals to be identified, same class X-ray pulsar it
Between signal recognition result of calculation by the use of instantaneous amplitude as characteristic vector it is smaller, in other X-ray pulsars it is corresponding away from
It is significantly greater from result of calculation.
Fig. 2 gives the signal recognition result of calculation of J1643-1224 and B1257+12X ray pulse stars.Wherein, horizontal seat
Mark represents the numbering of X-ray pulsar, and ordinate represents the result of calculation of feature space distance.X-ray pulsar signal recognition
As a result minimum value is numbered to the X-ray pulsar of induction signal in itself.
The calculating time based on instantaneous amplitude characteristic signal identification algorithm based on bispectrum feature signal recognition method with being carried out
Contrast is as shown in Figure 3 and Figure 4.Used computer condition is as follows：CPU is Pentium E5300, dominant frequency 2.6GHz, is inside saved as
2GB, software Matlab2010a.
From Fig. 3 and Fig. 4 as can be seen that during two X-ray pulsar observed quantities are recognized, based on instantaneous
The identification algorithm calculating time of amplitude Characteristics is much smaller than the calculating time based on bispectrum feature.With the increase of phase resolution,
Based on the calculating time of bispectrum feature discrimination method in power exponent increase, and based on the meter of instantaneous amplitude characteristic signal identification algorithm
Evaluation time linearly increases.
Summary simulation result can be obtained to draw a conclusion：
1) the signal recognition algorithm by the use of instantaneous amplitude as signal characteristic that this patent is proposed can be to X-ray pulse
Star is effectively recognized.
2) method that signal characteristic is extracted using Hilbert transform that this patent is proposed can preferably be met in real time
The demand of property.
Claims (4)
1. a kind of X-ray pulsar discrimination method based on Hilbert transform, it is characterised in that it comprises the following steps：
Step 1：Hilbert transform.
Step 2：Characteristics extraction.
Step 3：Minimum distance classifier is designed, is recognized.
A kind of 2. X-ray pulsar discrimination method based on Hilbert transform according to claim 1, it is characterised in that
Described step one is that X-ray pulsar signal is that the photon number received by one section of observation time internal X-ray detector tires out
The observation contour signal that product is formed is formed, for different X-ray pulsar signals, because X-ray pulsar has spy
Fixed physical features and space characteristics, traveling to interplanetary X-ray pulsar signal has geometry uniqueness.Cause
This, selects feature of the instantaneous amplitude feature that X-ray pulsar signal obtains after Hilbert transform needed for as identification.
Assuming that X-ray pulsar signal is x (t), the signal y (t) after Hilbert transform can be expressed asThis is x (t) and 1/ π τ convolution.
A kind of 3. X-ray pulsar discrimination method based on Hilbert transform according to claim 1, it is characterised in that
Described step two is：Converted by Hilbert, the real part x (t) and imaginary part y (t) of X-ray pulsar signal are constituted
One conjugate relation it is plural right, signal resolution form is expressed as S (t)=x (t)+jy (t), X-ray pulsar signal it is instantaneous
Amplitude Characteristics function A (t) can be expressed asInstantaneous amplitude can effectively reflect X ray arteries and veins
Rush the detailed information of star signal-energy.Therefore, the nominal contour feature set F in navigational route databasesF can be expressed ass=[v1
(t),v2(t),...,vm(t)]T, wherein m represent navigational route database in X-ray pulsar signal quantity.Assuming that i-th of X is penetrated
The characteristic vector of line pulsar signal is vi(t), then there is vi(t)=[xi(t),yi(t)]T.The X in navigational route database is extracted respectively
The instantaneous amplitude feature of the quasi- contour signal of ray pulse asterisk and observation contour signal simultaneously establishes X-ray pulsar signal standards wheel
Wide feature set Fs, correspondingly the instantaneous amplitude-frequency characteristic vector of X-ray pulsar observation contour signal is Vs(t), then there is Vs(t)=
[xs(t),ys(t)]T。
A kind of 4. X-ray pulsar discrimination method based on Hilbert transform according to claim 1, it is characterised in that
Described step three is：It is accurate according to identification after having obtained X-ray pulsar nominal contour and having observed the characteristic vector of contour signal
Minimum distance classifier can be then designed, calculates characteristic vector V respectivelys(t) with nominal contour feature set FsInterior each characteristic vector
Between space length.Assuming that there is m X-ray pulsar in navigational route database, then i-th X-ray pulsar (1≤i≤m)
Nominal contour characteristic vector viAnd V (t)s(t) distance can be expressed as
Therefore, identification space length collection can be obtained to be combined into：D={ D1,D2,...,Dm, work as DiWhen obtaining minimum value, corresponding X is penetrated
Line pulsar sequence number is identification result, i.e.,：Wherein L is the identification sequence number of X-ray pulsar.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201711082587.5A CN107894231A (en) | 2017-11-06 | 2017-11-06 | A kind of X-ray pulsar discrimination method based on Hilbert transform |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201711082587.5A CN107894231A (en) | 2017-11-06 | 2017-11-06 | A kind of X-ray pulsar discrimination method based on Hilbert transform |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN107894231A true CN107894231A (en) | 2018-04-10 |
Family
ID=61804241
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201711082587.5A Pending CN107894231A (en) | 2017-11-06 | 2017-11-06 | A kind of X-ray pulsar discrimination method based on Hilbert transform |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN107894231A (en) |
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110207689A (en) * | 2019-05-30 | 2019-09-06 | 西安电子科技大学 | A kind of pulsar signal denoising and discrimination method based on Wavelet Entropy |
| CN110986922A (en) * | 2019-12-30 | 2020-04-10 | 西安电子科技大学 | Method for acquiring X-ray pulsar short-time observation high signal-to-noise ratio contour |
| CN111189445A (en) * | 2020-01-14 | 2020-05-22 | 哈尔滨工业大学 | Pulsar identification method based on stochastic resonance |
| CN111649735A (en) * | 2020-06-12 | 2020-09-11 | 中国空间技术研究院 | Pulsar signal noise reduction method based on photon probability |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20050192719A1 (en) * | 2003-12-08 | 2005-09-01 | Suneel Ismail Sheikh | Navigational system and method utilizing sources of pulsed celestial radiation |
| CN103674031A (en) * | 2012-09-04 | 2014-03-26 | 西安电子科技大学 | Method for measuring attitude of spacecraft by using pulsar radiation vector and linear polarization information |
| CN105136138A (en) * | 2015-08-05 | 2015-12-09 | 西安电子科技大学 | X-ray pulsar photon signal identification method based on nuclear extreme learning machine |
| CN105371836A (en) * | 2015-12-18 | 2016-03-02 | 哈尔滨工业大学 | Mixed type fiber-optic gyroscope signal filtering method based on EEMD and FIR |
| CN105841714A (en) * | 2015-11-13 | 2016-08-10 | 湖南大学 | High speed X ray pulsar pulse profile delay measurement method |
| CN106840637A (en) * | 2017-03-24 | 2017-06-13 | 国网山东省电力公司莱芜供电公司 | Based on the GIS mechanical oscillation signal Time-Frequency Analysis Methods for improving HHT algorithms |
-
2017
- 2017-11-06 CN CN201711082587.5A patent/CN107894231A/en active Pending
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20050192719A1 (en) * | 2003-12-08 | 2005-09-01 | Suneel Ismail Sheikh | Navigational system and method utilizing sources of pulsed celestial radiation |
| CN103674031A (en) * | 2012-09-04 | 2014-03-26 | 西安电子科技大学 | Method for measuring attitude of spacecraft by using pulsar radiation vector and linear polarization information |
| CN105136138A (en) * | 2015-08-05 | 2015-12-09 | 西安电子科技大学 | X-ray pulsar photon signal identification method based on nuclear extreme learning machine |
| CN105841714A (en) * | 2015-11-13 | 2016-08-10 | 湖南大学 | High speed X ray pulsar pulse profile delay measurement method |
| CN105371836A (en) * | 2015-12-18 | 2016-03-02 | 哈尔滨工业大学 | Mixed type fiber-optic gyroscope signal filtering method based on EEMD and FIR |
| CN106840637A (en) * | 2017-03-24 | 2017-06-13 | 国网山东省电力公司莱芜供电公司 | Based on the GIS mechanical oscillation signal Time-Frequency Analysis Methods for improving HHT algorithms |
Non-Patent Citations (3)
| Title |
|---|
| XIAOYU LI,ET AL.: "Noise level estimation method with application to EMD-based signal denoising", 《JOURNAL OF SYSTEMS ENGINEERING AND ELECTRONICS》 * |
| 杨世锡等: "旋转机械振动信号基于EMD的希尔伯特变换和小波变换时频分析比较", 《中国电机工程学报》 * |
| 王璐等: "脉冲星辐射脉冲信号辨识的新算法", 《宇航学报》 * |
Cited By (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110207689A (en) * | 2019-05-30 | 2019-09-06 | 西安电子科技大学 | A kind of pulsar signal denoising and discrimination method based on Wavelet Entropy |
| CN110207689B (en) * | 2019-05-30 | 2022-09-16 | 西安电子科技大学 | Pulsar signal denoising and identifying method based on wavelet entropy |
| CN110986922A (en) * | 2019-12-30 | 2020-04-10 | 西安电子科技大学 | Method for acquiring X-ray pulsar short-time observation high signal-to-noise ratio contour |
| CN110986922B (en) * | 2019-12-30 | 2022-09-06 | 西安电子科技大学 | Method for acquiring X-ray pulsar short-time observation high signal-to-noise ratio contour |
| CN111189445A (en) * | 2020-01-14 | 2020-05-22 | 哈尔滨工业大学 | Pulsar identification method based on stochastic resonance |
| CN111649735A (en) * | 2020-06-12 | 2020-09-11 | 中国空间技术研究院 | Pulsar signal noise reduction method based on photon probability |
| CN111649735B (en) * | 2020-06-12 | 2021-11-16 | 中国空间技术研究院 | Pulsar signal noise reduction method based on photon probability |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN107894231A (en) | A kind of X-ray pulsar discrimination method based on Hilbert transform | |
| Jang et al. | Indoor localization with WiFi fingerprinting using convolutional neural network | |
| CN107301381A (en) | Recognition Method of Radar Emitters based on deep learning and multi-task learning strategy | |
| CN104134216B (en) | The laser point cloud autoegistration method described based on 16 dimensional features and system | |
| CN103596267A (en) | Fingerprint map matching method based on Euclidean distances | |
| CN102789578B (en) | Infrared remote sensing image change detection method based on multi-source target characteristic support | |
| CN102680825A (en) | Interference source identification method in system-grade electromagnetic compatibility fault diagnosis | |
| CN107230197A (en) | Tropical cyclone based on satellite cloud picture and RVM is objective to determine strong method | |
| CN103913166B (en) | A kind of punctate opacity of the cornea extracting method based on Energy distribution | |
| CN105654423A (en) | Area-based remote sensing image registration method | |
| CN108959182A (en) | The small feature loss gravitational field modeling method returned based on Gaussian process | |
| Akbar et al. | A novel benchmark RGBD dataset for dormant apple trees and its application to automatic pruning | |
| CN108120875A (en) | A kind of echo signal wide band detection method based on fast frequency spectrum template matches | |
| CN108107410B (en) | Abnormal radar joint judgment target detection method | |
| CN101916439A (en) | Space debris hyperspectral sequence detection method based on Hilbert-Huang transform | |
| CN101782964A (en) | Weight loss feather extraction method based on acceleration transducer | |
| CN116755005A (en) | Model-based method and system for monitoring magnetic field during lightning generation | |
| CN103645517A (en) | Comprehensive anomaly extraction method based on blind source separation technology and apparatus thereof | |
| CN113608190B (en) | Sea surface target detection method and system based on three characteristics of singular space | |
| CN113569695B (en) | Sea surface target detection method and system based on bispectrum three characteristics | |
| CN105550667A (en) | Stereo camera based framework information action feature extraction method | |
| CN112967241B (en) | Hyperspectral image anomaly detection method based on local gradient guidance | |
| Xiang et al. | Signal recognition method of x-ray pulsar based on cnn and attention module cbam | |
| CN105719272A (en) | Image characteristic point matching method for maintaining space structure | |
| CN104155692B (en) | Method and device for obtaining focus point of microseism |
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 | ||
| WD01 | Invention patent application deemed withdrawn after publication | ||
| WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20180410 |