CN106289239A - A kind of method eliminating the interference of wideband time domain in the pulsar data time of advent - Google Patents

Abstract

The present invention relates to astrophysics technical field, it is specifically related to a kind of method eliminating the interference of wideband time domain in the pulsar data time of advent, including: the pulse signal of observation pulsar, obtain the recurrent pulses profile of pulsar, carry out eliminating dispersion delay disposal to the recurrent pulses profile of pulsar, wideband radio interference in observation pulsar pulse data signal, eliminates the wideband time domain interference in the pulsar data time of advent.The present invention utilizes the dependency that wideband time domain is disturbed, subtracted each other by passage and eliminate the wideband time domain interference that the time domain interference of wideband eliminates in observations of pulsar data, obtain real pulsar profile, improving the precision of the pulsar time of advent, the observational study time of advent for pulsar provides more reliable guarantee.

Description

A kind of method eliminating the interference of wideband time domain in the pulsar data time of advent

Technical field

The present invention relates to astrophysics technical field, eliminate the pulsar data time of advent more particularly, it relates to a kind of The method of middle wideband time domain interference.

Background technology

Pulsar is one of 20th century big astronomical discovery sixties four, and its discovery development astrophysical to the present age produces Tremendous influence, correlational study work once obtained Nobel Prize in physics, this observation showing pulsar and theory twice respectively Research has great scientific meaning, and the observational study of pulsar is always the focus of field in astrophysics, and it is that research is extreme Physical property and " natural laboratory " of physics law under physical condition.

The most significant feature of pulsar is exactly to have the extremely stable rotation period, and long term monitoring result shows, its The long-time stability in cycle can compare favourably with atomic clock, even better than atomic clock.Radio astronomy observation is to utilize radio to look in the distance Mirror receives the pulse signal of pulsar from cosmic space, and the cycle that carries out the data of pulsar folds the spoke obtaining pulsar The signal penetrating profile and pulsar arrives the time of the earth, i.e. pulsar time of advent.The research of pulsar is concentrated mainly on arteries and veins Rushing radiation profile and the research time of advent of star, pulsar profile is studied, can be studied radiation characteristic and the spoke of pulsar itself It is emitted through journey；The pulsar research time of advent can study the stability of pulsar and the applied research of pulsar, such as pulsar meter Time and pulsar navigation research.

Owing to pulsar signal comprises substantial amounts of system noise and radio interference noise, the research to pulsar, first have to Consider interference is eliminated.

In digital processing field, conventional interference elimination method mainly has filtering and smooths, and pulsar signal Inherently a kind of pulse, relatively with data acquisition and data analysis for, itself be exactly one noise, conventional filtering and Smooth simultaneously can causing the signal of pulsar filters or smoothing effect, and the essence of pulsar research is through studying pulse The profile of star studies every character of pulsar, it is desirable to can not affect the wheel of pulsar itself while eliminating radio interference Profile shape and feature.So conventional filtering eliminates requirement with the smooth interference that can not meet observations of pulsar.

Have at observations of pulsar research time domain wideband radio interference elimination method the most in the world: self-adaptive routing；Zero Dispersion passage integration, eliminate capable and experienced interference method.

Self-adaptive routing is maximally effective interference elimination method in the observation of current PRF star, but self-adaptive routing needs Complicated hardware environment, needs two telescopes, receiver and Real-time Collection terminal, two telescopes be respectively primary telescope with With reference to telescope, this method is effective, but implements extremely difficult, it is necessary first to set up complicated hardware environment, simultaneously Adaptive filter algorithm needs debugging repeatedly, is the most only applied in Parkes observatory of Australia.

Zero dispersion passage integration, eliminate capable and experienced interference method be current pulsar Data processing eliminate wideband radio interference Main method, the method is to utilize passage zero dispersion amount integration, finds out noisy time point, then utilizes surrounding time point Data carry out interpolation and substitute the data of this time point, but, the method is the most a certain degree of shadow while eliminating radio interference Ring the contour shape of pulsar.

Therefore, those skilled in the art need urgently to solve the technical problem that exactly: the proposition how to innovate A kind of measure, eliminates the radio interference of time domain wideband while not affecting pulsar contour shape.

Summary of the invention

In order to solve the problems referred to above, the pulsar eliminating bigger dispersion disturbs in the time domain of radio wavelength band observation data, obtains Obtain more real pulsar profile and the time of advent.The technical solution adopted in the present invention is as follows: a kind of pulsar that eliminates arrives Reach the method for wideband time domain interference in time data, comprise the following steps:

S1: the pulse signal of observation pulsar；

Radio telescope is utilized to receive from pulsar signal in cosmic space, observation pulsar time of advent, described arteries and veins Rush star signal through radio telescope focus on after, after being amplified by high sensitive receiver mixing after, be digitized adopting Collection, in observations of pulsar, owing to pulsar is simultaneously emitted by multi-frequency pulse signal, all pulse signals are through interstellar space After propagation, the time that the signal of different frequency arrives the earth is different；

S2: obtain the recurrent pulses profile of pulsar；

Under the conditions of wideband, the high frequency resolution of acquisition pulse star signal, high time resolution data simultaneously, by right Passage carries out the dispersion delay laggard line period superposition of correction and obtains periodically pulsing profile；

S3: carry out eliminating dispersion delay disposal to the recurrent pulses profile of pulsar；

Due to the effect of dispersion, do not eliminate different frequency signals and postpone, the profile deformation of pulsar will be caused, or can't see Pulsar signal, observations of pulsar system is divided into multiple passage the broadband signal received, and calculates the electromagnetic wave of each passage Postpone, after deducting the signal delay of each passage, then carry out passage superposition, i.e. realize the recurrent pulses profile of pulsar is entered Row eliminates dispersion delay disposal；

S4: the wideband radio interference in observation pulsar pulse data signal；

After carrying out the recurrent pulses profile of pulsar eliminating dispersion delay disposal, have very in the place having pulse signal Many radio disturb, and radio interference is superimposed upon on pulse signal, on the profile of pulsar with cause serious impact the time of advent；

S5: eliminate the wideband time domain interference in the pulsar data time of advent；

Utilize the characteristic that the dispersion of pulsar postpones, by interference cancellation algorithm, eliminate in the pulsar data time of advent Wideband time domain interference.

Further, described step S2 includes:

S21: astronomical receiver receives the pulsar signal of the interstellar matter by cosmic space, by the frequency of pulsar signal Signal in rate bandwidth deltaf f amplifies, mixes down to intermediate frequency；

S22: using immediate frequency cable that intermediate-freuqncy signal is injected into terminal system, terminal system is the letter that frequency bandwidth is Δ f Number it is divided into v₁、v₂、v₃.....v_nN subchannel, each subchannel band a width of Δ f/n；

S23: with subchannel v_n/2On the basis of, calculate subchannel v₁、v₂、v₃.....v_nTo subchannel v_n/2Time delay Δ t₁、Δt₂、Δt₃......Δt_n；

The data of each subchannel are overlapped, it is thus achieved that pulsar by S24: correct the time delay of each subchannel Profile.

Further, the interstellar matter in described cosmic space act as the pulse signal of pulsar is caused delay.

Further, interference cancellation algorithm described in step S5 includes:

In actual observation, the numerical value of each sample record is designated as p_ij, wherein i is channel position, and scope is 0～n, and j is Sampling sequence number, the cycle of pulsar is T, and sample rate is Δ τ, and the employing in a cycle is counted as T/ Δ τ.

S51: passage corrects

Being normalized the data of τ the sampling of T/ Δ of n subchannel, the system correcting n sub-channel data is corresponding And high frequency attenuation, calculate the power response of each passage, counting cycle of passage carries out arithmetic mean, with any one passage Mean power response on the basis of carry out passage correction, recalculate the numerical value p ' of the every time sampling of each passage_ij；

S52: calculate dispersion delay and pulsar profile width

Calculate the dispersion time delay in whole bandwidth by dispersion delay time formula, obtain arteries and veins by pulsar star catalogue Rush the nominal contour of star, or utilize the nominal contour of historical data to calculate the profile width W of pulsar_psr；

S53: interference eliminates

Calculating the value of the channel data sampled point of pulsar, the value of described sampled point includes: pulsar signal value, radio are dry Disturbing value, system noise value, described pulsar signal value is only at the place with the presence of pulsar impulse radiation, described radio interference value For random disturbances noise figure, described system noise all existed within the whole cycle；

If dispersion Δ time delay t ＞ 2W_psr, in the cycle of pulsar is T, the i-th passage has pulsar signal, i.e. The value of the channel data sampled point of the i-th passage includes: pulsar signal value, radio interference value, system noise value, then the i-th+n/2 Passage does not has pulsar signal, and the value of the channel data sampled point of the i-th+n/2 passage includes: radio interference value, system noise Value；

The value of the channel data sampled point of the i-th passage is deducted the value of the channel data sampled point of the i-th+n/2 passage and get final product To pulsar signal value.

Further, described profile width W_psrPhase width for edge, signal both sides.

Further, described dispersion delay time formula is: Vt=4.148808DM [(v₀)^-2-(v_n)^-2], DM in formula For pulsar dispersion measure, unit is parsec every cubic centimetre, and described parsec is a kind of long measure astronomically, 1 second Gap is equal to 3.2615637769744079138 light-years, v₀、v_nFor frequency, unit GHz.

Compared to prior art, the invention have the benefit that

The present invention utilizes the dependency that wideband time domain is disturbed, and is subtracted each other the time domain interference elimination pulse eliminating wideband by passage Wideband time domain interference in star observation data, it is thus achieved that real pulsar profile, improves the precision of the pulsar time of advent, for arteries and veins The observational study time of advent rushing star provides more reliable guarantee.

Accompanying drawing explanation

Fig. 1 be the present invention eliminate the method flow diagram of wideband time domain interference in the pulsar data time of advent；

Fig. 2 is the pulse signal schematic diagram of the observation pulsar of the present invention；

Fig. 3 is that the pulsar pulse arrival time of the present invention is with frequency variation diagram；

Fig. 4 is the recurrent pulses profile schematic diagram obtaining pulsar of the present invention；

Fig. 5 be the present invention do not eliminate dispersion postpone time pulsar profile and time domain interferogram；

Fig. 6 is the pulsar profile after the elimination dispersion of the present invention and time domain interferogram；

Fig. 7 is the calculating pulse profile width indication figure of the present invention；

Fig. 8 embodiments of the invention do not eliminate pulsar profile when dispersion postpones and time domain interferogram；

Fig. 9 is De-dispersion but the profile diagram of pulsar after not eliminating interference in embodiments of the invention；

Figure 10 is to eliminate real pulsar profile diagram after time domain is disturbed in embodiments of the invention；

Figure 11 is to eliminate in embodiments of the invention after time domain is disturbed with true after system white noise replacement negative value pulse profile Real pulsar profile diagram.

Detailed description of the invention

Below in conjunction with specific embodiment, the present invention is further described below.

As it is shown in figure 1, the invention provides a kind of method eliminating the interference of wideband time domain in the pulsar data time of advent, Comprise the following steps:

S1: the pulse signal of observation pulsar；

As in figure 2 it is shown, utilize radio telescope to receive from pulsar signal in cosmic space, when observation pulsar arrives Between, described pulsar signal, after radio telescope focuses on, after being amplified by high sensitive receiver after mixing, is carried out Digital collection.

As it is shown on figure 3, in observations of pulsar, owing to pulsar is simultaneously emitted by multi-frequency pulse signal, all pulses Signal is after interstellar space is propagated, and the time that the signal of different frequency arrives the earth is different, and in Fig. 3, vertical coordinate is passage, Different passages is different frequency, and the signal of the signal rather low-frequency passage of hf channel arrives first takes things philosophically survey station.

S2: obtain the recurrent pulses profile of pulsar；

As shown in Figure 4, under the conditions of wideband, the high frequency resolution of acquisition pulse star signal, high time resolution simultaneously Data, obtain periodically pulsing profile by passage carries out the dispersion delay laggard line period superposition of correction；From Fig. 4 permissible Finding, pulse signal is very strong, has a lot of radio to disturb in the place having pulse signal simultaneously, and these interference are superimposed upon pulse letter On number, by the profile of pulsar with cause serious impact the time of advent, the most high-precision pulsar is seen the time of advent Survey research and the research of pulsar radiation mechanism.

Observations of pulsar system principle as shown in Figure 4, the signal of the pulsar interstellar matter (ISM) by cosmic space After, the pulsar signal of pulsar is caused delay by the effect of ISM, is that Δ f will believe by frequency bandwidth after astronomical receiver Number amplify, mix down to intermediate frequency, be injected into terminal system through immediate frequency cable, terminal system divides the signal that frequency bandwidth is Δ f Become v₁、v₂、v₃.....v_nN subchannel, each subchannel band a width of Δ f/n.With v_n/2On the basis of subchannel, calculate v₁、v₂、 v₃.....v_nEach subchannel is to v_n/2The time delay Δ t of passage₁、Δt₂、Δt₃......Δt_n, correct each passage time Between postpone after the data of each passage are overlapped, thus obtain the profile of pulsar.

S3: carry out eliminating dispersion delay disposal to the recurrent pulses profile of pulsar；

Due to the effect of dispersion, do not eliminate different frequency signals and postpone, the profile deformation of pulsar will be caused, or can't see Pulsar signal, observations of pulsar system is divided into multiple passage the broadband signal received, and calculates the electromagnetic wave of each passage Postpone, after deducting the signal delay of each passage, then carry out passage superposition, i.e. realize the recurrent pulses profile of pulsar is entered Row eliminates dispersion delay disposal.The strongest pulsar pulse signal is obtained as it is shown in figure 5, eliminate after dispersion postpones.

S4: the wideband radio interference in observation pulsar pulse data signal；

After carrying out the recurrent pulses profile of pulsar eliminating dispersion delay disposal, have very in the place having pulse signal Many radio disturb, and radio interference is superimposed upon on pulse signal, on the profile of pulsar with cause serious impact the time of advent；

S5: eliminate the wideband time domain interference in the pulsar data time of advent；

Utilize the characteristic that the dispersion of pulsar postpones, by interference cancellation algorithm, eliminate in the pulsar data time of advent Wideband time domain interference.

The size of pulsar dispersion is relevant to the distance of pulsar, and distance its dispersion values the most remote is the biggest, and the expression of dispersion is public Formula is as follows:

$DM=\underset{0}{\overset{d}{\∫}}{n}_{e}dl---\left(1\right)$

Wherein DM be pulsar dispersion, l be the distance of pulse, n_eFor milky way galaxy average electron density.

Different frequency electromagnetic wave be can be expressed as by time delay after interstellar space:

Vt=4.148808DM [(v₀)^-2-(v₁)^-2] (2)

Wherein Vt is the arrival earth receptor time delay of different frequency signals, and unit is millisecond, and DM is pulse star color Dissipating amount, unit is parsec every cubic centimetre, and parsec is a kind of long measure astronomically, and 1 parsec is equal to 3.2615637769744079138 light-year, v₀、v₁For frequency, unit GHz.

Observations of pulsar data disappear interference purpose when being to eliminate the arrowband interference in observations of pulsar band limits and wideband Territory is disturbed.The method generally eliminating pulsar interference is, finds interference on the passage of interference or time domain by certain algorithm Point, the interference on frequency domain eliminates the data of passage, time domain is use local value to substitute in time domain plus the random noise of rms Noise spot, this method can eliminate the passage of interference and the data of interference, it is thus achieved that preferably profile, but, if interference Being superimposed upon on the profile of pulsar, now the value in this moment should be:

T_sig=T_sys+T_psr+T_rfi (3)

In formula: T_sigTotal signal strength, T_psrThe signal of pulsar, T_sysFor system noise, T_rfiDisturb for radio.

In order to eliminate the impact of interference, and the value of system noise is used to substitute the sampled value in this moment, this moment Value is set to system noise, namely T_sig～=T_sys, have lost the radiant power of pulsar, this will improve the arrival of pulsar Temporal noise, can change the radiation profile of pulsar simultaneously, thus affect the precision of the time of advent of pulsar, to pulsar Radiation mechanism research impacts.

In the present invention, the numerical value recording each observed samples is p_ij, i is channel position, and scope is

0～n, j are sampling sequence number.The cycle of pulsar is T, and sample rate is Δ τ, and the employing in a cycle is counted as T/ Δ τ。

Described step S5 includes:

S51: passage corrects

Being normalized the data of τ the sampling of T/ Δ of n subchannel, the system correcting n sub-channel data is corresponding And high frequency attenuation, calculate the power response of each passage, counting cycle of passage carries out arithmetic mean, and computing formula is:

$\stackrel{\‾}{p_i}=\frac{\underset{1}{\overset{T/\mathrm{\Δ}\mathrm{\τ}}{\Σ}}{p}_{ij}}{\frac{T}{\mathrm{\Δ}\mathrm{\τ}}}---\left(4\right)$

On the basis of the mean power of any one passage responds, carry out passage correction, recalculate each of each passage The numerical value p ' of sampling_ij；As on the basis of kth passage, then:

$p_{ij}^{\′}=\frac{{\stackrel{\‾}{p}}_k{p}_{ij}}{{\stackrel{\‾}{p}}_i}---\left(5\right)$

S52: calculate dispersion delay and pulsar profile width

Calculate the dispersion time delay in whole bandwidth by dispersion delay time formula, formula (2) can obtain:

Vt=4.148808DM [(v₀)^-2-(v_n)^-2] (6)

Obtained the nominal contour of pulsar by pulsar star catalogue, or utilize the nominal contour of historical data to calculate pulse The profile width W of star_psr；Pulse profile width is the phase width at edge, signal both sides.

As it is shown in fig. 7, profile starts from 0.485, ending at 0.519, relative width is: (0.519-0.485) × T= 0.034T, T are the cycle of pulsar.

S53: interference eliminates

In the present invention, the characteristic that the dispersion of pulsar postpones, the time of arrival (toa) of the pulsar of different frequency are utilized Different, the signal of high frequency first arrives, the rear arrival of low frequency, and the due in of note high-frequency signal is t_psr-h, arriving of low frequency signal Reaching the moment is t_psr-l, then t_psr-h＜ t_psr-l, and disturb the signal same time to arrive, t_rfi-h=t_rfi-l, by calculating pulsar The delay of observation bandwidth of radiation profile and pulsar, if pulsar postpones to be more than the arteries and veins of twice in observation bandwidth internal dispersion Rush star bandwidth deltaf t=t_psr-l-t_psr-h＞ 2W_psr, after being calibrated by passage, passage is reported to the leadship after accomplishing a task to subtract each other and can effectively be eliminated penetrating of wideband Electrical interference, and do not affect the signal of pulsar.

The value of the channel data sampled point of pulsar can be calculated by formula (3):

$P_{ij}\∝T_{ij}^{sys}+T_{ij}^{psr}+T_{ij}^{rfi}---\left(7\right)$

WhereinFor system noise, the whole cycle all exists；Pulsar signal, only has in the place having impulse radiation Signal,Disturb for radio, for random disturbances noise.

The value of the channel data sampled point of pulsar includes following several situation:

(1) there is no pulse signal and an interference sample:

$P_{ij}\∝T_{ij}^{sys}---\left(8\right)$

(2) there is no pulse signal, a noisy sampling:

$P_{ij}\∝T_{ij}^{sys}+T_{ij}^{rfi}---\left(9\right)$

(3) sampling have pulsar signal, not interfering with:

$P_{ij}\∝T_{ij}^{sys}+T_{ij}^{psr}---\left(10\right)$

(4) there are pulsar signal, a noisy sampling:

$P_{ij}\∝T_{ij}^{sys}+T_{ij}^{psr}+T_{ij}^{rfi}---\left(11\right)$

Owing to the signal of pulsar has time delay with different frequency, as Δ t ＞ 2W_psrTime, the i-th passage has pulsar to believe Number time its sampled data be:

Or

And i-th_+n2Passage does not has pulse signal, and its sampled data is:

Or

Between different passages, have after being corrected by passage:

$T_{ij}^{sys}=T_{(i+n/2)j}^{sys},T_{ij}^{rfi}=T_{(i+n/2)j}^{rfi}$

So i passage and i+n/2 passage only remain the white noise of pulsar signal and system after subtracting each other

$p_{ij}-p_{(i+n/2)j}\∝T_{ij}^{psr}$

Finally to obtaining the profile of pulsar negative value after subtracting each other, substitute with the random number of the white noise scope of system.

So, can get the following conclusions that if dispersion Δ time delay t ＞ 2W_psr, in the cycle of pulsar is T, I-th passage has pulsar signal, and the value of the channel data sampled point of the i.e. i-th passage includes: pulsar signal value, radio disturb Value, system noise value, then the i-th+n/2 passage does not has pulsar signal, and the value bag of the channel data sampled point of the i-th+n/2 passage Include: radio interference value, system noise value, the value of the channel data sampled point of the i-th passage is deducted the port number of the i-th+n/2 passage Pulsar signal value is i.e. obtained according to the value of sampled point.

So to dispersion deferred gratification Δ t ＞ 2W_psrPulsar can to correct after channel data subtract each other, non- The signal of pulsar is retained while often effectively eliminating radio interference.

The measured data of pulsar PSR J1645-0317 in actual applications, is processed by the present invention.Fig. 8 represents Owing to dispersion postpones, do not eliminate dispersion delay and the pulse profile of time domain interference；Fig. 9 represents De-dispersion but does not eliminate interference Afterpulse profile；Figure 10 represents elimination time domain interference afterpulse profile, after Figure 11 represents elimination time domain interference, uses system white noise Real pulsar profile after substituting negative value pulse profile.As can be seen from Figure 8, the signal of pulsar by dispersion postpone and time The impact of territory interference, causes pulse profile to deform, or can't see pulsar signal；Postpone from fig. 9, it can be seen that eliminate dispersion After obtain the strongest pulse signal, but have many interference on pulse profile and background；From Figure 11, Figure 10 and Fig. 8, Fig. 9 Contrast can be seen that the interference in the elimination time domain that the present invention can be extraordinary obtains real pulse profile.

Above a kind of method eliminating the interference of wideband time domain in the pulsar data time of advent provided by the present invention is entered Having gone and be discussed in detail, principle and the embodiment of the present invention are set forth by specific case used herein, above enforcement The explanation of example is only intended to help to understand method and the core concept thereof of the present invention；General technology people simultaneously for this area Member, according to the thought of the present invention, the most all will change, in sum, and this explanation Book content should not be construed as limitation of the present invention.

Finally it is noted that the foregoing is only the preferred embodiments of the present invention, it is not limited to the present invention, Although being described in detail the present invention with reference to previous embodiment, for a person skilled in the art, it still may be used So that the technical scheme described in foregoing embodiments to be modified, or wherein portion of techniques feature is carried out equivalent, All within the spirit and principles in the present invention, any modification, equivalent substitution and improvement etc. made, should be included in the present invention's Within protection domain.

Claims (6)

1. one kind eliminates the method for wideband time domain interference in the pulsar data time of advent, it is characterised in that comprise the following steps:

S1: the pulse signal of observation pulsar；

Radio telescope is utilized to receive from pulsar signal in cosmic space, observation pulsar time of advent, described pulsar Signal is after radio telescope focuses on, and after being amplified by high sensitive receiver after mixing, is digitized gathering, In observations of pulsar, owing to pulsar is simultaneously emitted by multi-frequency pulse signal, all pulse signals are propagated through interstellar space After, the time that the signal of different frequency arrives the earth is different；

S2: obtain the recurrent pulses profile of pulsar；

Under the conditions of wideband, the high frequency resolution of acquisition pulse star signal, high time resolution data simultaneously, by passage Carry out the dispersion delay laggard line period superposition of correction and obtain periodically pulsing profile；

S3: carry out eliminating dispersion delay disposal to the recurrent pulses profile of pulsar；

Due to the effect of dispersion, do not eliminate different frequency signals and postpone, the profile deformation of pulsar will be caused, or can't see pulse Star signal, observations of pulsar system is divided into multiple passage the broadband signal received, and the electromagnetic wave calculating each passage postpones, After deducting the signal delay of each passage, then carry out passage superposition, i.e. realize the recurrent pulses profile of pulsar is disappeared Except dispersion delay disposal；

S4: the wideband radio interference in observation pulsar pulse data signal；

After carrying out the recurrent pulses profile of pulsar eliminating dispersion delay disposal, have in the place having pulse signal and much penetrate Electrical interference, radio interference is superimposed upon on pulse signal, on the profile of pulsar with cause serious impact the time of advent；

S5: eliminate the wideband time domain interference in the pulsar data time of advent；

Utilize the characteristic that the dispersion of pulsar postpones, by interference cancellation algorithm, eliminate the width in the pulsar data time of advent Frequently time domain interference.

2. the method for wideband time domain interference in the elimination pulsar data time of advent as claimed in claim 1, it is characterised in that Described step S2 includes:

S21: astronomical receiver receives the pulsar signal of the interstellar matter by cosmic space, by the frequency band of pulsar signal Signal in wide Δ f amplifies, mixes down to intermediate frequency；

S22: using immediate frequency cable that intermediate-freuqncy signal is injected into terminal system, terminal system divides the signal that frequency bandwidth is Δ f Become v₁、v₂、v₃…..v_nN subchannel, each subchannel band a width of Δ f/n；

S23: with subchannel v_n/2On the basis of, calculate subchannel v₁、v₂、v₃…..v_nTo subchannel v_n/2Time delay Δ t₁、Δ t₂、Δt₃……Δt_n；

The data of each subchannel are overlapped by S24: correct the time delay of each subchannel, it is thus achieved that the profile of pulsar.

3. the method for wideband time domain interference in the elimination pulsar data time of advent as claimed in claim 2, it is characterised in that The pulse signal of pulsar is caused delay by the acting as of the interstellar matter in described cosmic space.

4. the method for wideband time domain interference in the elimination pulsar data time of advent as claimed in claim 1, it is characterised in that Described in step S5, interference cancellation algorithm includes:

In actual observation, the numerical value of each sample record is designated as p_ij, wherein i is channel position, and scope is 0～n, and j is sampling Sequence number, the cycle of pulsar is T, and sample rate is Δ τ, and the employing in a cycle is counted as T/ Δ τ.

S51: passage corrects

Being normalized the data of τ the sampling of T/ Δ of n subchannel, the system correcting n sub-channel data is corresponding and high Frequency decay, calculates the power response of each passage, and counting cycle of passage carries out arithmetic mean, flat with any one passage All carry out passage correction on the basis of power response, recalculate the numerical value p of the sampling every time of each passage_i'_j；

S52: calculate dispersion delay and pulsar profile width

Calculate the dispersion time delay in whole bandwidth by dispersion delay time formula, obtain pulsar by pulsar star catalogue Nominal contour, or utilize the nominal contour of historical data to calculate the profile width W of pulsar_psr；

S53: interference eliminates

Calculating the value of the channel data sampled point of pulsar, the value of described sampled point includes: pulsar signal value, radio disturb Value, system noise value, described pulsar signal value is only in the place with the presence of pulsar impulse radiation, and described radio interference value is Random disturbances noise figure, described system noise all existed within the whole cycle；

If dispersion Δ time delay t ＞ 2W_psr, in the cycle of pulsar is T, the i-th passage has pulsar signal, and i.e. i-th leads to The value of the channel data sampled point in road includes: pulsar signal value, radio interference value, system noise value, then the i-th+n/2 passage does not has There is pulsar signal, and the value of the channel data sampled point of the i-th+n/2 passage includes: radio interference value, system noise value；

The value of the channel data sampled point that the value of the channel data sampled point of the i-th passage is deducted the i-th+n/2 passage i.e. obtains arteries and veins Rush star signal value.

5. the method for wideband time domain interference in the elimination pulsar data time of advent as claimed in claim 4, it is characterised in that Described profile width W_psrPhase width for edge, signal both sides.

6. the method for wideband time domain interference in the elimination pulsar data time of advent as claimed in claim 4, it is characterised in that Described dispersion delay time formula is: Vt=4.148808DM [(v₀)^-2-(v_n)^-2], in formula, DM is pulsar dispersion measure, unit For parsec every cubic centimetre, described parsec is a kind of long measure astronomically, and 1 parsec is equal to 3.2615637769744079138 light-year, v₀、v_nFor frequency, unit GHz.