CN107222226B

CN107222226B - The method and apparatus of harmonic carcellation interference

Info

Publication number: CN107222226B
Application number: CN201610160793.2A
Authority: CN
Inventors: 张志伟
Original assignee: Huawei Technologies Co Ltd
Current assignee: Huawei Technologies Co Ltd
Priority date: 2016-03-21
Filing date: 2016-03-21
Publication date: 2019-08-27
Anticipated expiration: 2036-03-21
Also published as: CN107222226A

Abstract

Present applicant proposes a kind of methods of harmonic carcellation interference, comprising: determines that the harmonic wave amount of moving, the harmonic wave amount of moving are the difference between the harmonic wave interference frequency point and the center frequency point for receiving signal；According to the reception signal and the harmonic wave amount of moving, the harmonic amplitude and harmonic phase of the harmonic interference signals for including in the reception signal are determined；According to the harmonic amplitude and the harmonic phase, the harmonic interference signals are determined；The harmonic interference signals are subtracted from the reception signal, so that the harmonic interference signals received in signal are eliminated.Therefore, method described herein, the harmonic signal efficiently solved is interfered signal bring is received, and improves the sensitivity behaviour of receiver, influences the performance of receiver by harmonic signal intensity, robustness with higher.

Description

The method and apparatus of harmonic carcellation interference

Technical field

The present invention relates to the method and apparatus of the harmonic carcellation interference in the communications field more particularly to the communications field.

Background technique

Narrowband interference is a major issue for needing to solve when wireless reception of signals processing.Clock harmonic interference is narrowband One kind of interference, it is a kind of interference signal of single-frequency point.The producing cause of harmonic wave interference have radio frequency (Radio Frequency, Referred to as " RF ") sources such as crystal oscillator harmonic wave and Analog Baseband (Analog Base Band, referred to as " ABB ") sampling clock harmonic wave, when The harmonic component of crystal oscillator or sampling clock is coupled to RF high frequency input port, then by after mixer (Mixer), reforming into reception letter Number with interior interference.

As RF chip area becomes, the wiring of small and smart terminal product is closer and closer, and harmonic wave interference problem can be cured Send out intractable.In general, the energy of harmonic wave interference is fixed, thus when receiver farther out from base station when, signal energy can be remote Lower than interfering energy.Therefore harmonic wave interference will lead to receiving signal sensitivity severe exacerbation, and net etc. is caused to seriously affect user The problem of experience.Meanwhile handling capacity when middle high s/n ratio also has greater loss, harmonic wave interference is stronger, caused performance loss It is bigger.

Summary of the invention

In view of this, the embodiment of the invention provides a kind of method and apparatus of harmonic carcellation, to solve harmonic wave to reception The interference problem of signal.So as to avoid receiving signal sensitivity caused by harmonic wave interference deteriorate, fall net etc. influence user's body The problem of testing, and reduce performance loss.

In a first aspect, providing a kind of method of harmonic carcellation interference, which is characterized in that harmonic wave interference frequency point, which is located at, to be received In the reception bandwidth of signal, which comprises

Determine that the harmonic wave amount of moving, the harmonic wave amount of moving are the harmonic wave interference frequency point and the center frequency for receiving signal Difference between point；

According to the reception signal and the harmonic wave amount of moving, the harmonic interference signals for including in the reception signal are determined Harmonic amplitude and harmonic phase；

According to the harmonic amplitude and the harmonic phase, the harmonic interference signals are determined；

The harmonic interference signals are subtracted from the reception signal, so that the harmonic wave received in signal is dry Signal is disturbed to be eliminated.

Therefore, method described herein, the harmonic signal efficiently solved is interfered signal bring is received, and is mentioned The high sensitivity behaviour of receiver, influences the performance of receiver by harmonic signal intensity, robustness with higher.

It is described according to the reception signal and the harmonic wave amount of moving as another embodiment, determine the reception letter The harmonic amplitude and harmonic phase for the harmonic interference signals for including in number, comprising:

According to the harmonic wave amount of moving, the harmonic wave interference frequency point is moved to the position of the center frequency point for receiving signal It sets；

According to the signal value of each sampled point of the reception signal in unit length, the reception signal is calculated Average value；

According to the average value for receiving signal, the institute for the harmonic interference signals for including in the reception signal is determined State harmonic amplitude and the harmonic phase.

Optionally, it is assumed that crystal oscillator fundamental frequency is F, and the center frequency point for receiving signal is f_c, harmonic wave interference frequency point is f_i, then f= f_i-f_cIt can be defined as the amount of moving of Harmonics elimination.Wherein, center frequency point f_cWith harmonic wave interference frequency point f_iIt can respectively indicate Are as follows:

f_c=F × (N+1/M) (1)

f_i=F × N (2)

Wherein, N indicates that the integer frequency of crystal oscillator, 1/M indicate decimal frequency multiplication, frequency unit Hz.

Optionally, before determining the harmonic wave amount of moving, the method also includes: determine that the harmonic wave interference frequency point is located at institute It states in the reception bandwidth for receiving signal.

Assuming that reception bandwidth is BW, then when interference frequency point and center frequency point interval | f_i-f_cWhen |≤BW/2, it is believed that interference is fallen In reception bandwidth.At this point, the baseband receiving signals with harmonic wave interference can indicate are as follows:

Y (t)=x (t)+Ae^j(2πf+θ)+n(t) (3)

Wherein, f=f_i-f_cFor the amount of moving of Harmonics elimination, A is harmonic amplitude, and θ is harmonic phase.In addition, n (t) is indicated The noise section of signal is received, x (t) indicates noiseless and noise signal section, Ae^j(2πf+θ)Indicate that the harmonic wave for receiving signal is dry Disturb signal.

By the way that signal y (t) and e will be received^-j2πftComplex multiplication is carried out, realizes and harmonic wave frequency point is moved into center frequency point position. Reception signal y'(t after being moved).

Y'(t)=y (t) × e^-j2πft=[x (t)+n (t)] × e^-j2πft+Ae^jθ (7)

In formula (7), [x (t)+n (t)] × e-^j2πftFor single-tone modulation signal, if [x (t)+n (t)] is relative to list Sound amplitude very little (sensitivity scene meets the hypothesis), then after the single-tone number of cycles that adds up is secondary, the result of this is intended to 0, i.e., Obtain result shown in formula (8):

According to formula (8), it is estimated that the harmonic amplitude and harmonic phase of harmonic interference signals.

Optionally, for demodulation modes and alien frequencies mode, it is assumed that frequency deviation is Δ f, then receives the center frequency point f of signal_cWith Harmonic wave interference frequency point f_iIt is respectively as follows:

f_c=(F+ Δ f) × (N+1/M) (4)

f_i=(F+ Δ f) × N (5)

It is f=f that the nominal harmonic wave amount of moving, which can be calculated, according to formula (1) and (2)_i-f_c=-F/M, but there are frequency deviation Δ f When, it is f '=- (F+ Δ f)/M, and the frequency deviation value of AFC estimation is f that the interval is practical_AFC=Δ f (N+1/M), frequency deviation scale factor For Fratio=Δ f/F.Therefore, it can use Given information and calculate the harmonic wave amount of moving that actual needs uses are as follows:

F '=- (F+ Δ f)/M=- (F+F*Fratio) */M=f* (1+Fratio) (6)

When there are frequency deviation, it is believed that harmonic wave interference frequency point is f+ Δ f relative to the distance of center frequency point, then:

Y'(t)=y (t) × e^-j2πft=[x (t)+n (t)] × e^-j2πft+Ae^j2πΔfte^jθ (9)

Since Δ f is relatively small, e^j2πΔftPeriod it is very long, it is believed that in one section of sampling point, the tone signal e^j2πΔftAmplitude and phase be constant, therefore the segment data can be added up and is averaging, to estimate that harmonic wave interference is believed Number amplitude and phase, and for eliminating the harmonic wave interference of lower one piece of data.

Optionally, before determining the harmonic wave amount of moving, the method also includes: determine that the unit for receiving signal is long Degree, and the reception signal is segmented according to the unit length, which indicates any in the reception signal The sampling number for including in one segment signal, such as the n-th segment signal, the n is positive integer and n can less than or equal to reception signal The total number of segment being divided into, or receive total number of segment in signal for harmonic wave estimation.

It should be understood that can be estimated according to using the n-th segment signal after determining the harmonic interference signals in n-th segment signal Harmonic interference signals, to eliminate the harmonic interference signals in the (n+1)th segment signal, the harmonic wave estimated according to the (n+1)th segment signal Interference signal, to eliminate the harmonic interference signals in the n-th+2 segment signal, and so on.It can also be estimated according to the n-th segment signal Harmonic interference signals, to eliminate subsequent all signals or any one segment signal or multi-segment signal.

Therefore, because we by the way of segmentation, will be used for lower one piece of data when the result of leading portion estimation, in realization It is a kind of mode of stream treatment, such processing mode is very low to hard-wired timing requirements, simplifies module and realizes again Miscellaneous degree.

On the other hand, it in addition to the first segment data does not have estimated result, can not carry out outside Harmonics elimination, follow-up data estimated value It all comes into force in real time, therefore solves the problems, such as to need that parameter is waited to restrain in the prior art.

As another embodiment, before the determining harmonic wave amount of moving, the method also includes:

According to the size of the frequency shift (FS) for receiving signal, Harmonics elimination mode, the Harmonics elimination mode packet are determined Include grade mode, demodulation modes or alien frequencies mode.

According to the Harmonics elimination mode, the number of the sampled point in the unit length is determined, in the unit length The number of sampled point and the size of the frequency shift (FS) be inversely proportional.

It is divided in more detail it should be understood that the Harmonics elimination mode can according to need when product is realized, it is described Determine that Harmonics elimination mode institute's basis principle is the size of residual frequency offset to be treated.

Wherein, when the Harmonics elimination mode is demodulation modes, unit length can be 1024；When the Harmonics elimination When mode is alien frequencies mode, unit length can be 256；When the Harmonics elimination mode is grade mode, unit length can Think 128.

As another embodiment, when determining that the Harmonics elimination mode is demodulation modes or alien frequencies mode or described When the frequency shift (FS) is not present in harmonic wave interference frequency point, the determining harmonic wave amount of moving, comprising:

According to the size of the frequency shift (FS), the harmonic wave interference frequency point and the center frequency point for receiving signal are determined；

According to the harmonic wave interference frequency point and the center frequency point for receiving signal, the above-mentioned harmonic wave amount of moving is determined.

As another embodiment, when determining the Harmonics elimination mode is grade mode, the determining harmonic wave is moved Amount, comprising:

According to the size of the frequency shift (FS) of the harmonic wave interference frequency point, the harmonic wave interference frequency point and center frequency are determined Point；

According to the harmonic wave interference frequency point and the center frequency point, the first harmonic amount of moving, the second harmonic amount of moving are determined With the third harmonic amount of moving, the second harmonic amount of moving is less than the first harmonic amount of moving, the third harmonic amount of moving Greater than the first harmonic amount of moving, and the difference of the first harmonic amount of moving and the second harmonic amount of moving is equal to described The difference of the third harmonic amount of moving and the first harmonic amount of moving；

From the first harmonic amount of moving, the second harmonic amount of moving and the third harmonic amount of moving, institute is selected State the harmonic wave amount of moving.

As another embodiment, it is described from the first harmonic amount of moving, the second harmonic amount of moving and described In the three harmonic wave amounts of moving, the harmonic wave amount of moving is selected, comprising:

According to the reception signal and the first harmonic amount of moving, the second harmonic amount of moving and the third The harmonic wave amount of moving, determining and the first harmonic amount of moving, the second harmonic amount of moving and the third harmonic amount of moving one One corresponding first harmonic interference signal, second harmonic interference signal and third harmonic interference signal；

According to the first harmonic interference signal, the second harmonic interference signal and the third harmonic interference signal, It is determining to be corresponded with the first harmonic interference signal, the second harmonic interference signal and the third harmonic interference signal The first signal power, second signal power and third signal power；

Determine the peak signal function in first signal power, the second signal power and the third signal power The corresponding harmonic wave amount of moving of rate is the target harmonic wave amount of moving；

In multiple unit lengths that the reception signal includes, the first harmonic amount of moving, described the are counted respectively The two harmonic wave amounts of moving and the third harmonic amount of moving are targeted the number of the harmonic wave amount of moving；

Determine that the corresponding harmonic wave amount of moving of the number maximum value is the harmonic wave amount of moving.

It optionally, is f, f- Δ f' and f+ Δ f' tri- according to the harmonic wave amount of moving according to formula (6) for grade mode Frequency point determines the harmonic wave amount of moving.Such as Δ f'=10KHz, it is assumed that sampling rate 30.72Msps, 1ms include 30720 A sampling point, unit length 128, then (30720-4096)/128=208 sections of calculating can be carried out in total by receiving data.For Every one piece of data is utilized respectively [f-10KHz, f, f+10KHz] these three amounts of moving and calculates corresponding harmonic interference signals Harmonic amplitude and harmonic phase, and corresponding harmonic interference signals of these three available amounts of moving.Utilize formula P= I²+Q²The corresponding power of three harmonic interference signals is calculated separately, wherein I indicates the imaginary part of plural harmonic signal, and Q indicates plural number The real part of harmonic signal.Alpha (Alpha) filtering is carried out later, compares maximum value therein, and the maximum value is corresponding The harmonic wave amount of moving be determined as the target harmonic wave amount of moving.It is simultaneously that the harmonic wave amount of moving for being targeted the harmonic wave amount of moving is corresponding Counter add 1.In multiple unit lengths that the reception signal includes, the first harmonic amount of moving, described the are counted The two harmonic wave amounts of moving and the third harmonic amount of moving are targeted the number of the harmonic wave amount of moving, and by its corresponding harmonic wave Harmonics elimination of the amount of moving as the gear estimated, for subsequent sampling point.Where determining the harmonic wave amount of moving after gear, then The harmonic wave estimation of receipt of subsequent is carried out in the way of the fixed harmonic wave amount of moving and is eliminated.

As another embodiment, the method also includes:

According to the signal value of each sampled point of the direct current signal in the reception signal in unit length, institute is calculated State the average value of direct current signal；

The average value of the direct current signal is subtracted from the reception signal, so that described straight in the reception signal Stream signal is eliminated.

Second aspect provides a kind of device of harmonic carcellation interference, for executing above-mentioned first aspect or first aspect Any possible implementation in method.Specifically, the device of harmonic carcellation interference includes:

First determining module, for determining that the harmonic wave amount of moving, the harmonic wave amount of moving are the harmonic wave interference frequency point and institute State the difference between the center frequency point for receiving signal；

First determining module is also used to:

Processing module, the harmonic wave interference letter determined for subtracting first determining module from the reception signal Number, so that the harmonic interference signals received in signal are eliminated.

Therefore, the device of the interference of harmonic carcellation described in the embodiment of the present invention can efficiently solve harmonic signal docking Collection of letters bring interference problem, and the sensitivity behaviour of receiver is improved, make the performance of receiver not by harmonic signal The influence of intensity, robustness with higher.

The third aspect provides a kind of device of harmonic carcellation interference, including processor, memory, bus system, reception Device and transmitter, the transmitter and the receiver are respectively used to send and receive information during communication, described to deposit For storing instruction, the processor is used to execute the instruction of memory storage to reservoir, and to storing in the memory Instruction execution so that the processor execute either first aspect or first aspect face possibility implementation in side Method.Specifically, the processor is specifically used for:

Fourth aspect provides a kind of computer-readable medium, and for storing computer program, which includes The instruction of the method in any possible implementation for executing first aspect or first aspect.

Based on the above-mentioned technical proposal, the method for the harmonic carcellation interference of the embodiment of the present invention, the harmonic wave letter efficiently solved Number to the interference of signal bring is received, and the sensitivity behaviour of receiver is improved, believe that the performance of receiver by harmonic wave The influence of number intensity, robustness with higher.

Detailed description of the invention

In order to illustrate the technical solution of the embodiments of the present invention more clearly, will make below to required in the embodiment of the present invention Attached drawing is briefly described, it should be apparent that, drawings described below is only some embodiments of the present invention, for For those of ordinary skill in the art, without creative efforts, it can also be obtained according to these attached drawings other Attached drawing.

Fig. 1 is the scheme that the frequency domain harmonic of the prior art is eliminated.

Fig. 2 is the scheme that the time harmonic of the prior art is eliminated.

Fig. 3 is the structural schematic diagram of the wireless receiver of the embodiment of the present invention.

Fig. 4 is the schematic flow chart of the method for the harmonic carcellation interference of the embodiment of the present invention.

Fig. 5 is the schematic flow chart of the method for the harmonic carcellation interference of the embodiment of the present invention.

Fig. 6 is the schematic block diagram of the gear estimation of the grade mode of the embodiment of the present invention.

Fig. 7 is the schematic flow chart of the amplitude and phase of the determination harmonic interference signals of the embodiment of the present invention.

Fig. 8 is the schematic block diagram of the method for the harmonic carcellation interference of the embodiment of the present invention.

Fig. 9 is the schematic diagram of the device of the harmonic carcellation interference of the embodiment of the present invention.

Figure 10 is the schematic diagram of the device of the harmonic carcellation interference of the embodiment of the present invention.

Specific embodiment

Following will be combined with the drawings in the embodiments of the present invention, and technical solution in the embodiment of the present invention carries out clear, complete Site preparation description, it is clear that described embodiment is a part of the embodiments of the present invention, rather than whole embodiments.Based on this hair Embodiment in bright, those of ordinary skill in the art's every other reality obtained without making creative work Example is applied, all should belong to the scope of protection of the invention.

Typical Harmonics elimination scheme at present, is divided into frequency domain and two kinds of time domain.Fig. 1 is that the frequency domain harmonic of the prior art is eliminated Scheme.As shown in Figure 1, frequency domain harmonic eliminate first will receive signal by Fourier transformation (Fourier Transform, Referred to as " FFT ") frequency domain is transformed to, the power of each frequency point signal is counted later and is judged with preset threshold value, is more than Threshold value is then considered to interfere.The frequency point that finally will be deemed as interference is set to zero, or is set to the mean power of adjacent frequency.But The frequency domain harmonic cancellation scheme has the disadvantages that the processing of 1, frequency domain harmonic only inhibits interference, and there is no completely It eliminates, therefore the modules such as channel estimation and demodulation is also needed to do corresponding processing 2, due to being the processing on frequency domain, receiver The interference effect that signal procesing in time domain part such as AGC, cell searching etc. is subject to cannot be eliminated.3, frequency domain harmonic elimination is only used for Harmonic wave is much higher than the scene of signal, because the fluctuation of signal itself just has 8 decibels of (Decible, abbreviations under fading channel " dB ") more than, the probability of miscarriage of justice of Interference Detection is very big.4, when product is realized, cost of implementation and complexity, the point of FFT are limited to Number can not be too big, and when harmonic wave is not the integral multiple of FFT subcarrier, the influence of harmonic wave interference can be diffused into multiple subcarriers On.

Fig. 2 is the scheme that the time harmonic of the prior art is eliminated.Typical method is by adaptive iteration algorithm, training The corresponding notch filter (Notch Filter) of interference frequency point, time-domain signal reach the mesh for eliminating interference after filtering out , the notch filter that should be obtained by adaptive iteration algorithm can be referred to as adaptive notch filter.But this method has Following disadvantage: 1, the depth and width needs of Notch Filter are adjusted according to interference strength.Because of different interference intensity The sub-carrier number of influence is different, and setting is not allowed to will lead to interference frequency point and adjacent subcarrier is eliminated not enough or crossed and eliminates, algorithm Robustness it is difficult to ensure that.2, the convergence time of adaptive algorithm and interference strength are related, when interference strength is higher than signal on frequency domain When less than 10dB, convergence time needs access to a subframe, and interference cannot be completely eliminated at this time.

In consideration of it, the problems such as embodiment of the present invention is inadequate for the too many robustness of existing scheme constraint condition, proposes one kind Realize that simple, performance robust Time Domain Piecewise is averaged Harmonics elimination scheme.

Since clock harmonic frequency point is that crystal oscillator fundamental frequency or Analog Baseband (Analog Base Band, referred to as " ABB ") are adopted The integral multiple of sample clock, therefore the frequency point position that the harmonic wave frequency point may be fallen in is fixed.Therefore, on actual product plate only There are fixed several frequency points that can there are problems that harmonic wave interference.

It should be understood that these frequency points can be obtained in laboratory test, and as Given information for carrying out Harmonics elimination, It can be obtained using other estimation methods, which is not limited by the present invention.

Fig. 3 is the structural schematic diagram of the wireless receiver of the embodiment of the present invention.Fig. 3 shows the RF module of wireless receiver 310, analog-digital converter (Analog-to-Digital Converter, referred to as " ADC ") 320, harmonic management module 330, automatic Gain control (Automatic Gain Control, referred to as " AGC ") module 340, channel estimation module 350, cell searching mould Block 360, demodulation module 370, automatic frequency control (Automatic Frequency Control, referred to as " AFC ") module 380 With decoding module 390.The signal that RF module 310 receives carries out analog-to-digital conversion in analog-digital converter 320, then at harmonic wave Manage the elimination that harmonic interference signals are carried out in module 330.Harmonic management module 330 is located at after ADC module 320, is used for basis The present invention implements the harmonic Elimination Method and eliminates to the harmonic wave interference received in signal.AGC module 340 is for docking The amplification factor of the collection of letters number automatically adjusts, and channel estimation module 350 is used for when receiving signal, to the relevant parameter of channel Estimated.It then receives signal and carries out demodulation coding in demodulation module 370 and decoding module 390 respectively.In addition, if with Family equipment is in cell search phase, then executing the work of cell search phase, such as user by cell search module 360 Initial synchronisation etc. between equipment and base station, wherein AFC module 380 is used to estimate the frequency shift (FS) situation for receiving signal Meter, to obtain the frequency offset estimate f and frequency shift (FS) scale factor that receive channel.

Fig. 4 is the schematic flow chart of the method for the harmonic carcellation interference of the embodiment of the present invention.Method shown in Fig. 4 can be with It is executed by harmonic management module 330, the method for the Harmonics elimination includes:

410, determine that the harmonic wave amount of moving, the harmonic wave amount of moving are in the harmonic wave interference frequency point and the reception signal Difference between heart frequency point.

420, according to the reception signal and the harmonic wave amount of moving, determine the harmonic wave interference for including in the reception signal The harmonic amplitude and harmonic phase of signal.

430, according to the harmonic amplitude and the harmonic phase, determine the harmonic interference signals.

440, the harmonic interference signals are subtracted from the reception signal, so that described humorous in the reception signal Interference wave signal is eliminated.

Specifically, position of the harmonic management module 330 according to harmonic wave interference frequency point and the center frequency point for receiving signal, meter Calculate the harmonic wave amount of moving, difference of the harmonic wave amount of moving between the harmonic wave interference frequency point and the center frequency point.Then harmonic management Module 330 estimates the harmonic wave width for receiving the harmonic interference signals for including in signal according to signal and the harmonic wave amount of moving is received Degree and harmonic phase, and according to the harmonic amplitude estimated and harmonic phase, reconstruct harmonic interference signals.So as to from rear The harmonic interference signals are subtracted in continuous reception signal, to realize the elimination to harmonic interference signals.Wherein, the harmonic wave amount of moving It is referred to as harmonic wave and moves frequency point.

Optionally, before carrying out the method, harmonic management module 330 can first determine that the unit of the reception signal is long Degree, and the reception signal is segmented according to the unit length, which indicates any one section in the reception signal The sampling number for including in signal, such as the n-th segment signal, the n is positive integer and n is less than or equal to reception signal and can be divided At total number of segment, or receive signal in for harmonic wave estimation total number of segment.Such as unit length is indicated with L here, it indicates It include L sampled point in every segment signal.

Specifically, harmonic management module 330 determines the unit length of the reception signal, and according to the unit length to reception Signal is segmented.Then to the n-th segment signal therein, the harmonic wave width for the harmonic interference signals for including in the n-th segment signal is estimated Degree and harmonic phase, and according to the harmonic amplitude estimated and harmonic phase, reconstruct harmonic interference signals.So as to from connecing The harmonic interference signals are subtracted in the (n+1)th segment signal in the collection of letters number, to realize to the harmonic wave for receiving other segment signals in signal The elimination of interference signal.

It is a kind of stream in realization when the result of leading portion estimation is used for lower one piece of data since we are by the way of segmentation The mode of water process, such processing mode are very low to hard-wired timing requirements, simplify module implementation complexity.

On the other hand, it in addition to the first segment data does not have estimated result, can not carry out outside Harmonics elimination, follow-up data estimated value It all comes into force in real time, therefore there is no need the problem of waiting parameter to restrain in the prior art.

It should be understood that after harmonic management module 330 reconstructs the harmonic interference signals in the n-th segment signal, it can be according to utilization The harmonic interference signals that n-th segment signal estimates are believed to eliminate the harmonic interference signals in the (n+1)th segment signal according to (n+1)th section Number harmonic interference signals estimated, to eliminate the harmonic interference signals in the n-th+2 segment signal, and so on.Harmonic management mould The harmonic interference signals that block 330 can also be estimated according to the n-th segment signal, to eliminate subsequent all signals or any one section of letter Number or multi-segment signal, which is not limited by the present invention.

It should also be understood that for sake of simplicity, below just for being illustrated for the n-th segment signal in the reception signal, i.e., after Reception signal in the description of face all refers to n-th segment signal.

It is specifically described respectively to 410 to 440 below with reference to Fig. 4 to Fig. 8.

Harmonic management module 330 can determine unit length according to the business scenario being presently in, thus according to difference Unit length use different Harmonics elimination modes, thus targetedly to receive signal in harmonic wave interference disappear It removes.

Optionally, as another embodiment, before the determining harmonic wave amount of moving, the method also includes:

Although the frequency point of harmonic wave interference be it is fixed, on actual product, all there is deviation in crystal oscillator, cause to receive There are frequency shift (FS) (later referred to as frequency deviations) for machine.AFC module 380 in receiver estimates frequency deviation value by receiving signal, and Frequency deviation compensation is carried out in RF module 310 or base band.In AFC module 380 after AFC is calibrated, system residual frequency deviation generally exists Within several hundred hertz (Hertz, referred to as " Hz "), inventive algorithm performance is influenced smaller.But at the beginning of the cell in practical application Begin search, the stages such as frequency sweep, AFC be it is no convergent, frequency deviation at this time is likely to be breached 15KHz or so.For cell initial search etc. AFC does not have convergent scene, simulating, verifying, and inventive algorithm can guarantee that cell of frequency deviation when within certain frequency range is searched Without hesitation can, such as within the scope of 5KHz.Therefore frequency deviation stepping algorithm is used, is respectively f, f- Δ f' and f+ Δ according to the harmonic wave amount of moving Tri- frequency points of f' estimate harmonic wave, and wherein Δ f' is stepping step-length, then select interference estimate power most strong from three frequency points Carry out eliminate.Here the processing energy of the 5KHz to illustrate can be configured according to the concrete condition of real application systems. In addition, stepping step delta f' can also be configured according to the concrete condition of real application systems, such as here with Δ f'= It is illustrated for 10KHz, i.e., estimates harmonic wave according to the harmonic wave amount of moving is tri- frequency points of f, f-10KHz and f+10KHz respectively.

Therefore, in embodiments of the present invention, according to the size of residual frequency departure, it is (corresponding that Harmonics elimination is divided into grade mode The application scenarios such as initial cell search, frequency sweep and BUCKET-BASED SEARCHING BACKGROUNDSearch), demodulation modes (scenes such as demodulation, common frequency measurement) and alien frequencies mould The Harmonics eliminations modes such as formula (scenes such as alien-frequency district search).As shown in table 1, corresponding unit under different Harmonics elimination modes Length L is different.Table 1 is the parameter configuration table under different Harmonics elimination modes.Wherein, table is first is that with long term evolution (Long Term Evolution, referred to as " LTE ") receiver for, parameter setting corresponds to the sampling rate of 30.72Msps.

Parameter configuration table under the different Harmonics elimination modes of table 1

It should be understood that product can be divided in more detail as needed when realizing, basis principle is exactly to be treated The size of residual frequency departure, the present invention to Harmonics elimination mode division without limitation.

Optionally, as another embodiment, before executing 410, the method for harmonic carcellation interference further include:

Determine that harmonic wave interference frequency point is located in the reception bandwidth for receiving signal.

Specifically, when the determining working frequency points and reception bandwidth for receiving signal of harmonic management module 330, and determine that harmonic wave is dry The interference frequency point for disturbing signal is fallen in reception bandwidth.Suitable Harmonics elimination mould is selected further according to different frequency deviation sizes later Formula eliminates harmonic wave interference.

It to sum up describes, the detailed process of the method for the Harmonics elimination of the embodiment of the present invention is summarized as Fig. 5 by we.Fig. 5 is The schematic flow chart of the method for the harmonic carcellation interference of the embodiment of the present invention.Method shown in fig. 5 can be by harmonic management mould Block 330 executes.This method specifically includes:

510, it determines in current reception bandwidth with the presence or absence of harmonic interference signals.

520, determine the working frequency points and reception bandwidth for receiving signal.

530, judge whether harmonic wave interference frequency point falls in reception bandwidth.

Specifically, it is assumed that the center frequency point for receiving signal is f_i, harmonic wave interference frequency point is f_i, receive the bandwidth of operation of signal For BW.Harmonic wave interference module 330 judges | f_i-f_c| whether it is less than or equal to BW/2, works as satisfaction | f_i-f_cWhen |≤BW/2, it is believed that humorous Wave interference frequency point is located in the reception bandwidth for receiving signal.

540, Harmonics elimination mode is determined according to current business scene.

Specifically, it if harmonic wave interference frequency point is located in the reception bandwidth for receiving signal, needs to carry out the harmonic wave interference It eliminates, harmonic management module 330 combines table 1, needs Harmonics elimination mode to be used according to the selection of current business scene.The harmonic wave The selection of elimination mode is really to be determined by the frequency deviation size under different scenes.

550, select grade mode to eliminate harmonic wave interference.

When specifically, under the scenes such as initial cell search, frequency sweep and BUCKET-BASED SEARCHING BACKGROUNDSearch, due to not having when AFC under the scene Have convergent, frequency deviation at this time is very big, therefore harmonic management module 330 selects grade mode to carry out subsequent harmonic wave interference elimination Work.

551, stepping frequency point and unit length are determined under grade mode.

Wherein, corresponding unit length is shorter under the grade mode, and needs and judge harmonic wave gear, that is, executes 552.This In different harmonic wave gear can respectively correspond these three different harmonic wave amounts of moving of f, f- Δ f' and f+ Δ f'.F, f- Δ f' and f+ Each of Δ f' harmonic wave amount of moving all corresponds to a gear.

552, grade mode judges harmonic wave gear.

Specifically, harmonic management module 330 selects interference to estimate from these three harmonic wave amounts of moving of f, f- Δ f' and f+ Δ f' Evaluation power is strongest to be eliminated.Specifically judge the mode of harmonic wave gear later to specifically describing in the elaboration of 430 parts.

553, it harmonic wave estimation and eliminates under grade mode.

560, select demodulation modes to eliminate harmonic wave interference.

Specifically, smaller in the frequency deviation generated when the work such as common frequency measurement and demodulation, therefore harmonic management module 330 It selects demodulation modes to carry out subsequent harmonic wave interference and eliminates work.

561, stepping frequency point and unit length are determined under demodulation modes.

Wherein, corresponding unit length is longer under the grade mode.

562, it harmonic wave estimation and eliminates under demodulation modes.

570, select alien frequencies mode to eliminate harmonic wave interference.

Specifically, when carrying out alien-frequency district search, harmonic management module 330 selects alien frequencies mode to determine that harmonic wave interference is believed Number and it is eliminated.

571, stepping frequency point and unit length are determined under demodulation modes.

Wherein, corresponding unit length is longer under the grade mode.

572, it harmonic wave estimation and eliminates under demodulation modes.

Harmonic management module 330 according to its corresponding unit length, and calculates after having selected Harmonics elimination mode The harmonic wave amount of moving, in 420 carry out harmonic interference signals amplitude and phase estimation.Specifically, harmonic management module 330 According to signal and the harmonic wave amount of moving is received, the harmonic amplitude and harmonic phase for receiving the harmonic interference signals for including in signal are determined.

During the elimination of harmonic interference signals, if the frequency deviation value that the AFC module 380 in Fig. 3 estimates directly exists RF module 310 adjusts phaselocked loop (Phase Locked Loop, referred to as " PLL "), then arriving harmonic management module 330 can consider There is no frequency deviation, directly calculates the harmonic wave amount of moving according to nominal frequency point.

As an example it is assumed that crystal oscillator fundamental frequency is F, the center frequency point for receiving signal is f_c, harmonic wave interference frequency point is f_i, then f =f_i-f_cIt can be defined as the amount of moving of Harmonics elimination.Wherein, center frequency point f_cWith harmonic wave interference frequency point f_iTable can be distinguished It is shown as:

f_c=F × (N+1/M) (1)

f_i=F × N (2)

Y (t)=x (t)+Ae^j(2πf+θ)+n(t) (3)

If the offset estimation value is compensated in base band, in 330 signal of harmonic management module there are frequency deviation, this The time-harmonic wave amount of moving needs to convert frequency deviation value.Assuming that frequency deviation is Δ f, then the center frequency point f of signal is received_cWith harmonic wave interference frequency Point f_iIt is respectively as follows:

f_c=(F+ Δ f) × (N+1/M) (4)

f_i=(F+ Δ f) × N (5)

F '=- (F+ Δ f)/M=- (F+F*Fratio) */M=f* (1+Fratio) (6)

The above-mentioned determination harmonic wave amount of moving, the i.e. method of f' in formula (6) are suitable for demodulation modes and alien frequencies mode. It is respectively f, f- Δ f' and f+ according to the harmonic wave amount of moving using frequency deviation stepping algorithm for grade mode according to the description of front Tri- frequency points of Δ f' estimate harmonic wave, are illustrated by taking Δ f'=10KHz as an example here, and then selection interference is estimated from three frequency points The strongest frequency point of evaluation power, as the harmonic wave amount of moving that progress Harmonics elimination process is utilized in the data of receipt of subsequent.

The harmonic wave amount of moving is determined under demodulation modes and alien frequencies mode it is described above how determining.It is retouched below with reference to Fig. 6 State the determination of the harmonic wave amount of moving under grade mode.Since grade mode is not restrained for AFC such as initial cell search and frequency sweeps Scene, the accurate frequency deviation value of receiver can not be known at this time, here for example with the data of 1ms, to estimate the gear of frequency deviation. Fig. 6 is the schematic block diagram of the gear estimation of the grade mode of the embodiment of the present invention.

For example, as shown in Figure 6, it is assumed that sampling rate 30.72Msps, unit are conversion rate (Million Samples per Second, referred to as " Msps ").1ms includes 30720 sampling points, setting unit length L=128, then in Fig. 6 (30720-4096)/128=208 sections of calculating can be carried out in total by receiving data.For every one piece of data, harmonic management module 330, which are utilized respectively [f-10KHz, f, f+10KHz] these three amounts of moving, calculates harmonic wave interference estimated value, i.e. harmonic interference signals Harmonic amplitude and harmonic phase, and the corresponding harmonic interference signals of these three available amounts of moving.According to three obtained A harmonic interference signals, harmonic management module 330 utilize formula P=I²+Q²Calculate the corresponding function of three harmonic interference signals Rate, wherein I indicates the imaginary part of plural harmonic signal, and Q indicates the real part of plural harmonic signal.Carry out Alpha (Alpha) later Filtering, compares maximum value therein, and the corresponding harmonic wave amount of moving of the maximum value is determined as the target harmonic wave amount of moving.Simultaneously The corresponding counter of the harmonic wave amount of moving for being targeted the harmonic wave amount of moving is added 1, tri- frequencies of f-10KHz, f and f+10KHz The corresponding counter of point can be indicated with Cnt0, Cnt1 and Cnt2 respectively.

For receiving the current power (Power, referred to as " P ") of data, history value is set as Plod, enables Alpha here Filter factor is Alpha=0.2, then filtered value are as follows:

Pcurrent=(1-Alpha) * Plod+Alpha*P；// calculating for the first time directly enables Pcurrent=P；

Plod=Pcurrent；

Wherein, the parameter configurations such as data length, unit length L and Alpha filter factor for estimating gear are to show Example, product can be configured according to the actual situation when realizing, which is not limited by the present invention.

Define the corresponding filtered performance number of [f-10KHz, f, f-10KHz] this third gear be respectively [Pcurrent0, Pcurrent1, Pcurrent2], then every section of maximum value selected in three, and its corresponding counter is added 1.By 208 sections Statistics after, select the maximum value in Cnt0, Cnt1 and Cnt2, and by its corresponding harmonic wave amount of moving as the shelves estimated Position, the Harmonics elimination for subsequent sampling point.

Where determining the harmonic wave amount of moving after gear, receipt of subsequent is carried out according still further to the mode of the fixed harmonic wave amount of moving Harmonic wave estimation and elimination.

Optionally, as another embodiment, before the elimination of the carry out harmonic interference signals, this method further include:

Illustrate by taking grade mode as an example, before determining the harmonic wave amount of moving, harmonic management module 330 determines direct current first The length Ndc of (Direct Current, referred to as " DC ") data, and according to each sampling in the unit length of direct current signal The signal value of point, calculates the average value of direct current signal, and the average value of the direct current signal is finally subtracted from reception signal, so that The DC component received in signal is eliminated.Such as shown in Fig. 6, harmonic management module 330 utilizes the point meter of front Ndc=4096 DC component is calculated, subsequent sampling point first subtracts the DC value and carries out operation again.

Wherein, the data length Ndc for calculating direct current signal is example, and product can carry out according to the actual situation when realizing Configuration, which is not limited by the present invention.

Therefore, it by using the estimation method of input parameter DC, i.e., is counted according to statistics when module starts using front Ndc The average value of DC signal, the support DC that first disappears carry out harmonic wave estimation and elimination again, very so as to avoid harmonic wave frequency point and center frequency point When close, influence that DC component generates harmonic wave estimated accuracy.

After the harmonic wave amount of moving has been determined in 410, can the method according to 420, using the harmonic wave amount of moving to humorous The harmonic amplitude and harmonic phase of wave signal are estimated.

Fig. 7 is the schematic flow chart of the amplitude and phase of the determination harmonic interference signals of the embodiment of the present invention.Optionally, As another embodiment, as shown in fig. 7, it is described according to the reception signal and the harmonic wave amount of moving in 420, determine institute State the harmonic amplitude and harmonic phase for receiving the harmonic interference signals for including in signal, comprising:

421, according to the harmonic wave amount of moving, the harmonic wave interference frequency point is moved into the center frequency point for receiving signal Position.

422, according to the signal value of each sampled point of the reception signal in unit length, calculate the reception letter Number average value；

423, according to the average value for receiving signal, determine the harmonic interference signals for including in the reception signal The harmonic amplitude and the harmonic phase.

Below with reference to Fig. 8 specifically describe harmonic management module 330 how to estimate the harmonic interference signals harmonic amplitude and Harmonic phase.

Harmonic wave interference frequency point is moved the position for receiving the center frequency point of signal by harmonic management module 330 first, specifically As shown in Figure 8.Fig. 8 is the schematic block diagram of the method for the harmonic carcellation interference of the embodiment of the present invention.Wherein, harmonic management module The determination of harmonic wave interference amplitude and phase is completed in 330 dotted box portion in fig. 8.

As shown in figure 8, sampling rate F of the tone signal generation module according to the amount of moving f and baseband receiving signals_s, generate Plural tone signal:Wherein n=0,1,2 ..., complex table Show and plural number is taken to signal.Sampling point indexes n and receives the sampling point one-to-one correspondence of signal.Negative sign is taken to imaginary part later, obtains e^-j2πft, By the n-th segment signal y (t) and e that will receive signal^-j2πftComplex multiplication is carried out, realizes and moves harmonic wave interference frequency point to center frequency Point position.Reception signal y'(t after being moved).

Y'(t)=y (t) × e^-j2πft=[x (t)+n (t)] × e^-j2πft+Ae^jθ (7)

In formula (7), [x (t)+n (t)] × e^-j2πftFor single-tone modulation signal, if [x (t)+n (t)] is relative to list Sound amplitude very little (sensitivity scene meets the hypothesis), then after the single-tone number of cycles that adds up is secondary, the result of this is intended to 0, i.e., Obtain result shown in formula (8):

Harmonic management module 330 can estimate the harmonic amplitude and harmonic phase of harmonic interference signals according to formula (8).

It should be understood that, if reception signal is larger relative to harmonic power, parameter can not be accurately estimated in above-mentioned hypothesis, It receives signal and is higher than received signal strength indicator device (Received Signal Strength Indicator, abbreviation " RSSI ") instruction signal strength when, time harmonic elimination algorithm may be invalid, lead to penalty.Therefore the present invention is real Method described in example is applied for harmonic power greater than in the case where received signal power.

Y'(t)=y (t) × e^-j2πft=[x (t)+n (t)] × e^-j2πft+Ae^j2πΔfte^jθ (9)

In 430, harmonic management module 330 determines harmonic interference signals according to the harmonic amplitude and harmonic phase.

Specifically, harmonic wave interference module 330 can reconstruct Ae after determining the amplitude and phase of harmonic wave^j2πΔfte^jθAnd This is subtracted from receiving in signal, to realize the elimination for receiving the harmonic interference signals in signal to next section.

In 440, harmonic management module 330 subtracts harmonic interference signals from the reception signal, so that receiving letter The harmonic interference signals in number are eliminated.Harmonic wave shown in Fig. 8 updates estimation point and refers to that harmonic management module 330 is directed to The every L point of the harmonic amplitude and harmonic phase of harmonic interference signals updates primary, it can in each section of progress harmonic interference signals Estimation, to keep to the real-time updates of harmonic interference signals, to improve the accuracy of harmonic wave interference elimination.

Therefore, the method for the interference of harmonic carcellation described in the embodiment of the present invention, the harmonic signal efficiently solved is to reception The interference of signal bring, and the sensitivity behaviour of receiver is improved, make the performance of receiver not by harmonic signal intensity It influences, robustness with higher.

In addition, due to by the way of stream treatment, when simplifying product realization, and not needing additional convergence Between, therefore there is higher practicability.

It should be understood that magnitude of the sequence numbers of the above procedures are not meant that the order of the execution order, the execution of each process is suitable Sequence should be determined by its function and internal logic, and the implementation process of the embodiments of the invention shall not be constituted with any limitation.

Above in conjunction with Fig. 1 to Fig. 8, the method that harmonic carcellation interference according to an embodiment of the present invention is described in detail, under The device of harmonic carcellation interference according to an embodiment of the present invention will be described in detail in conjunction with Fig. 9 and Figure 10 in face.

Fig. 9 is the schematic diagram of the device of the harmonic carcellation interference of the embodiment of the present invention.Harmonic carcellation shown in Fig. 9 The device 900 of interference can be used in executing each process of the earlier figures 4 into the embodiment of the method for Fig. 8.Harmonic carcellation interference Device 900 can be above-mentioned harmonic wave interference module 330, can execute harmonic wave interference module 330 execution each process.Such as Shown in Fig. 9, the device of harmonic carcellation interference includes determining module 901 and processing module 902.Wherein, the first determining module 901 For:

First determines that the harmonic wave amount of moving, the harmonic wave amount of moving are in the harmonic wave interference frequency point and the reception signal Difference between heart frequency point；

Processing module 902, the harmonic wave interference determined for subtracting the first determining module 901 from the reception signal Signal, so that the harmonic interference signals received in signal are eliminated.

As another embodiment, the first determining module 901 is specifically used for:

As another embodiment, described device further includes the second determining module 903, is determined in the first determining module 901 Before the harmonic wave amount of moving, the second determining module 903 is used for:

As another embodiment, when the second determining module 903 determines that the Harmonics elimination mode is demodulation modes or different When the frequency shift (FS) is not present in frequency mode or the harmonic wave interference frequency point, the first determining module 901 is specifically used for:

According to the harmonic wave interference frequency point and the center frequency point for receiving signal, the harmonic wave amount of moving is determined.

As another embodiment, when the second determining module 903 determines that the Harmonics elimination mode is grade mode, the One determining module 901 is specifically used for:

Optionally, the first determining module 901 is also used to:

Processing module is also used to, and the average value of the direct current signal is subtracted from the reception signal, so that described connect The direct current signal in the collection of letters number is eliminated.

As shown in Figure 10, the embodiment of the invention also provides a kind of device of harmonic carcellation interference, harmonic carcellation interference Device include processor 1001, memory 1002, bus system 1003, receiver 1004 and transmitter 1005.Wherein, it handles Device 1001, memory 1002 and receiver 1004 are connected by bus system 1003, and memory 1002 for storing instruction, is handled Device 1001 is used to execute the instruction of the storage of memory 1002, and controls and receives device 1004 and receive information.Wherein processor 1001, deposit Reservoir 1002, bus system 1003, receiver 1004 and transmitter 1005 can be realized by one or more chips.For example, Processor 1001, memory 1002, bus system 1003, receiver 1004 and transmitter 1005 can be fully integrated in a core In piece or processor 1001, bus system 1003, receiver 1004 and transmitter 1005 can integrate in a chip and Memory 1002 is integrated in another chip, and concrete form is herein without limitation.Wherein, processor 1001 is used for:

It should be noted that the memory in the embodiment of the present invention can be volatile memory or nonvolatile memory, or can Including both volatile and non-volatile memories.Wherein, nonvolatile memory can be read-only memory (Read-Only Memory, referred to as " ROM "), programmable read only memory (Programmable ROM, referred to as " PROM "), erasable programmable Read-only memory (Erasable PROM, referred to as " EPROM "), electrically erasable programmable read-only memory (Electrically EPROM, referred to as " EEPROM ") or flash memory.Volatile memory can be random access memory (Random Access Memory, referred to as " RAM "), it is used as External Cache.By exemplary but be not restricted explanation, the RAM of many forms It can use, such as static random access memory (Static RAM, referred to as " SRAM "), dynamic random access memory (Dynamic RAM, referred to as " DRAM "), Synchronous Dynamic Random Access Memory (Synchronous DRAM, referred to as " SDRAM "), double data Rate synchronization dynamic random access memory (Double Data Rate SDRAM, referred to as " DDR SDRAM "), enhanced synchronization Dynamic random access memory (Enhanced SDRAM, referred to as " ESDRAM "), synchronized links dynamic random access memory (Synch Link DRAM, referred to as " SLDRAM ") and direct rambus random access memory (Direct Rambus RAM, Referred to as " DR RAM ").The memory of system and method described herein is intended to include but is not limited to these to be suitble to any other The memory of type.

Those of ordinary skill in the art may be aware that list described in conjunction with the examples disclosed in the embodiments of the present disclosure Member and algorithm steps can be realized with the combination of electronic hardware or computer software and electronic hardware.These functions are actually It is implemented in hardware or software, the specific application and design constraint depending on technical solution.Professional technician Each specific application can be used different methods to achieve the described function, but this realization is it is not considered that exceed The scope of the present invention.

In several embodiments provided herein, it should be understood that disclosed systems, devices and methods, it can be with It realizes by another way.For example, the apparatus embodiments described above are merely exemplary, for example, the unit It divides, only a kind of logical function partition, there may be another division manner in actual implementation, such as multiple units or components It can be combined or can be integrated into another system, or some features can be ignored or not executed.Another point, it is shown or The mutual coupling, direct-coupling or communication connection discussed can be through some interfaces, the indirect coupling of device or unit It closes or communicates to connect, can be electrical property, mechanical or other forms.

The unit as illustrated by the separation member may or may not be physically separated, aobvious as unit The component shown may or may not be physical unit, it can and it is in one place, or may be distributed over multiple In network unit.It can select some or all of unit therein according to the actual needs to realize the mesh of this embodiment scheme 's.

It, can also be in addition, the functional units in various embodiments of the present invention may be integrated into one processing unit It is that each unit physically exists alone, can also be integrated in one unit with two or more units.

It, can be with if the function is realized in the form of SFU software functional unit and when sold or used as an independent product It is stored in a computer readable storage medium.Based on this understanding, technical solution of the present invention is substantially in other words The part of the part that contributes to existing technology or the technical solution can be embodied in the form of software products, the meter Calculation machine software product is stored in a storage medium, including some instructions are used so that a computer equipment (can be a People's computer, server or network equipment etc.) it performs all or part of the steps of the method described in the various embodiments of the present invention. And storage medium above-mentioned includes: that USB flash disk, mobile hard disk, read-only memory (Read-Only Memory, ROM), arbitrary access are deposited The various media that can store program code such as reservoir (Random Access Memory, RAM), magnetic or disk.

The above description is merely a specific embodiment, but scope of protection of the present invention is not limited thereto, any Those familiar with the art in the technical scope disclosed by the present invention, can readily occur in various equivalent modifications or replace It changes, these modifications or substitutions should be covered by the protection scope of the present invention.Therefore, protection scope of the present invention should be with right It is required that protection scope subject to.

Claims

1. a kind of method of harmonic carcellation interference, which is characterized in that harmonic wave interference frequency point is located in the reception bandwidth for receiving signal, The described method includes:

Determine the harmonic wave amount of moving, the harmonic wave amount of moving be the harmonic wave interference frequency point and the center frequency point for receiving signal it Between difference；

According to the harmonic wave amount of moving, the harmonic wave interference frequency point is moved to the position of the center frequency point for receiving signal；

According to the signal value of each sampled point of the reception signal in unit length, being averaged for the reception signal is calculated Value；

According to the average value for receiving signal, the described humorous of the harmonic interference signals for including in the reception signal is determined Wave amplitude and the harmonic phase；

The harmonic interference signals are subtracted from the reception signal, so that the harmonic wave interference letter received in signal It number is eliminated.

2. the method as described in claim 1, which is characterized in that before the determining harmonic wave amount of moving, the method is also wrapped It includes:

According to the size of the frequency shift (FS) for receiving signal, determine that Harmonics elimination mode, the Harmonics elimination mode include point Shelves mode, demodulation modes or alien frequencies mode；

According to the Harmonics elimination mode, the number of the sampled point in the unit length, adopting in the unit length are determined The number of sampling point and the size of the frequency shift (FS) are inversely proportional.

3. method according to claim 2, which is characterized in that when determine the Harmonics elimination mode be demodulation modes or alien frequencies When the frequency shift (FS) is not present in mode or the harmonic wave interference frequency point, the determining harmonic wave amount of moving, comprising:

4. method according to claim 2, which is characterized in that when determining the Harmonics elimination mode is grade mode, institute State the determining harmonic wave amount of moving, comprising:

According to the size of the frequency shift (FS) of the harmonic wave interference frequency point, the harmonic wave interference frequency point and the center frequency point are determined；

According to the harmonic wave interference frequency point and the center frequency point, the first harmonic amount of moving, the second harmonic amount of moving and are determined The three harmonic wave amounts of moving, the second harmonic amount of moving are less than the first harmonic amount of moving, and the third harmonic amount of moving is greater than The first harmonic amount of moving, and the difference of the first harmonic amount of moving and the second harmonic amount of moving is equal to the third The difference of the harmonic wave amount of moving and the first harmonic amount of moving；

From the first harmonic amount of moving, the second harmonic amount of moving and the third harmonic amount of moving, select described humorous The wave amount of moving.

5. method as claimed in claim 4, which is characterized in that described to be removed from the first harmonic amount of moving, the second harmonic In shifting amount and the third harmonic amount of moving, the harmonic wave amount of moving is selected, comprising:

According to the reception signal and the first harmonic amount of moving, the second harmonic amount of moving and the third harmonic The amount of moving, it is determining a pair of with the first harmonic amount of moving, the second harmonic amount of moving and the third harmonic amount of moving one First harmonic interference signal, second harmonic interference signal and the third harmonic interference signal answered；

According to the first harmonic interference signal, the second harmonic interference signal and the third harmonic interference signal, determine With the first harmonic interference signal, the second harmonic interference signal and the third harmonic interference signal one-to-one One signal power, second signal power and third signal power；

Determine the maximum signal power pair in first signal power, the second signal power and the third signal power The harmonic wave amount of moving answered is the target harmonic wave amount of moving；

In multiple unit lengths that the reception signal includes, the first harmonic amount of moving, described second humorous is counted respectively The wave amount of moving and the third harmonic amount of moving are targeted the number of the harmonic wave amount of moving；

6. the method as described in any one of claims 1 to 5, which is characterized in that the method also includes:

According to the signal value of each sampled point of the direct current signal in the reception signal in unit length, calculate described straight Flow the average value of signal；

The average value of the direct current signal is subtracted from the reception signal, so that the direct current letter received in signal It number is eliminated.

7. a kind of device of harmonic carcellation interference, which is characterized in that harmonic wave interference frequency point is located in the reception bandwidth for receiving signal, Described device includes:

First determining module, for determining the harmonic wave amount of moving, the harmonic wave amount of moving is the harmonic wave interference frequency point and described connects Difference between the center frequency point of the collection of letters number；

First determining module is also used to:

Processing module, the harmonic interference signals determined for subtracting first determining module from the reception signal, So that the harmonic interference signals received in signal are eliminated.

8. device as claimed in claim 7, which is characterized in that described device further includes the second determining module, described first Before determining module determines the harmonic wave amount of moving, second determining module is used for:

9. device as claimed in claim 8, which is characterized in that when second determining module determines the Harmonics elimination mode When the frequency shift (FS) is not present for demodulation modes or alien frequencies mode or the harmonic wave interference frequency point, described first determines mould Block is specifically used for:

10. device as claimed in claim 8, which is characterized in that when second determining module determines the Harmonics elimination mould When formula is grade mode, first determining module is specifically used for:

11. device as claimed in claim 10, which is characterized in that first determining module is specifically used for:

12. the device as described in any one of claim 7 to 11, which is characterized in that first determining module is also used to:

The processing module is also used to, and the average value of the direct current signal is subtracted from the reception signal, so that described connect The direct current signal in the collection of letters number is eliminated.