CN103199889A - Field programmable gata array (FPGA) implementation method of iteration frequency domain anti-interference algorithm - Google Patents
Field programmable gata array (FPGA) implementation method of iteration frequency domain anti-interference algorithm Download PDFInfo
- Publication number
- CN103199889A CN103199889A CN2013100770581A CN201310077058A CN103199889A CN 103199889 A CN103199889 A CN 103199889A CN 2013100770581 A CN2013100770581 A CN 2013100770581A CN 201310077058 A CN201310077058 A CN 201310077058A CN 103199889 A CN103199889 A CN 103199889A
- Authority
- CN
- China
- Prior art keywords
- data
- sum
- ram1
- read
- module
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Images
Landscapes
- Complex Calculations (AREA)
- Noise Elimination (AREA)
Abstract
The invention provides a field programmable gata array (FPGA) implementation method of an iteration frequency domain anti-interference algorithm. The FPGA implementation method of the iteration frequency domain anti-interference algorithm includes the steps: conducting windowing operation after AD sample is conducted to satellite signals which are converted down into intermediate frequency; selecting appropriate threshold optimal coefficients after a fast Fourier transform (FFT) algorithm is conducted; deducing a self-adaptive iteration threshold and reducing the value of frequency points with frequency domain amplitudes larger than the threshold to threshold value; and finally conducting inverse fast Fourier transform (IFFT) and merging and outputting the results after the IFFT. By means of the FPGA implementation method of the iteration frequency domain anti-interference algorithm, an interference detection threshold is set self-adaptively, dynamic interference rejection capability is good, the phenomenon that hardware resources are occupied due to the fact that middle process data need to be stored in repeated judgments, the FPGA, namely the technology of exchanging time for resources, is fully used, and hardware resource consumption is reduced greatly.
Description
Technical field
The hardware that the present invention relates to a kind of anti-interference algorithm realizes, particularly a kind of FPGA of the anti-interference algorithm of frequency domain based on the iteration thresholding realizes technology, is used for judgement and inhibition that spread spectrum radio communications system disturbs.
Background technology
Wireless communication technology plays an increasingly important role in the modern communication field, but because wireless signal is extremely faint, be subjected to various natural or artificial interference easily and can't use, therefore, need to introduce the antijamming capability that interference mitigation technology improves wireless communication system.Because artificial disturbance mostly is the arrowband and disturbs, effective Suppression of narrow band interference technology can greatly be improved the performance of communication system.
In anti-interference algorithm, the judgement of interference and inhibition are the cores of whole algorithm.Wherein the trap method is one of algorithm of using always.So-called trap method, exactly signal is transformed into frequency domain from time domain after, amplitude and the thresholding of frequency domain are made comparisons, sagging the signal frequency point amplitude greater than thresholding, finish interference inhibition work.
The document 2 Suppression of narrow band interference research [electronic information countermeasure techniques; pp; 51-53; 2009] of adaptive threshold " under the low signal-to-noise ratio based on " proposes be approximately the narrowband Gaussian signal at spread-spectrum signal under the noiseless condition through behind the DFT, its envelope Rayleigh distributed, square obeys index distribution of envelope.By the derivation of exponential distribution probability density formula, propose to optimize COEFFICIENT K 〉=8 o'clock at single order thresholding in the threshold algorithm, under the condition of noiseless existence, can not cause damage to useful signal.Simultaneously, disclose a kind of segmentation interference suppression algorithm, chosen different thresholdings by the power that calculates current demand signal.But this algorithm complex height, expense is big.
In sum, existing document has much been used the narrowband Gaussian model at the anti-interference algorithm of frequency domain, and choose decision threshold according to this characteristic of square obeys index distribution of its envelope, but when specific implementation, find, because of square numerical value of envelope often very big, can take very many storage resources, this is to be reluctant to see in the specific implementation, and existing document fails to provide the jamproof FPGA realization of the iterative frequency-domain of disturbing at the arrowband.
Summary of the invention
In order to overcome the deficiencies in the prior art, the present invention propose a kind of based on single order apart from the jamproof FPGA realization of the frequency domain of iteration self-adapting thresholding technology.From this angle of narrowband Gaussian signal envelope Rayleigh distributed, choose suitable thresholding and optimize coefficient, derive the adaptive iteration thresholding, the frequency of frequency domain amplitude greater than thresholding sunk to threshold value, can guarantee that useful signal can all pass through, can reduce simultaneously the storage resources of hardware greatly, finish the inhibition work of interference.
The technical solution adopted for the present invention to solve the technical problems may further comprise the steps:
(1) will carry out the AD sampling to the satellite-signal that changes to intermediate frequency down, the data that obtain are output as two-way, and one the tunnel carries out 1/2 window length delay, and two paths of data is input to the windowing module respectively;
(2) windowing module is carried out the windowing operation for the data of input, and window function adopts broad sense hamming window, and window length is L, adopts the quantification of 8 bit data bit wides; After the windowing operation is finished, data input FFT module;
(3) use a degree of depth to store as the data that the RAM of L checks input in the FFT module, read data afterwards and be transported to FFT nuclear, carry out the FFT computing.FFT is set to Pipelined, and Streaming I/O pattern, operating frequency are f
In, f
InIdentical with the AD sampling clock, the output of its computing has real part Re and imaginary part Im two parts, merges into one road signal with these two parts with by the mould value R that the real part imaginary part calculates jointly according to senior middle school's low level among the present invention, i.e. { Re, Im, R} is transported to and disturbs identification and suppress module;
(4) each disturbs identification as follows with the course of work that suppresses module:
A. two-port RAM nuclear RAM1 and RAM2 that to set up two degree of depth be L=512, data input selected cell is selected RAM1 work, and the data of FFT output are with clock f
InStore among the RAM1, the L point that calculates mould value R in institute's deposit data simultaneously add up and
And threshold value
K is that thresholding is optimized coefficient;
B. after RAM1 was filled with, data input selected cell was selected RAM2 work, and the data of FFT output are deposited among the RAM2, and the L point that calculates the deposit data mould value R of institute simultaneously among the repeating step a adds up and SUM and threshold value TH
1Work;
C. in the process of RAM2 data storages, to the data among the RAM1 with clock f
sRead f
s=N*f
In, N is iterations; Can get, can be to the data read among the RAM1 3 times in this following period of time of RAM2 data storages; A flag register that bit wide is identical with the RAM1 degree of depth is set, in to RAM1 first pass process of reading, will represents a part of R of mould value in the data that read
jWith TH
1Compare, j=1, if 2,3...L is greater than TH
1Then the corresponding position of flag register puts 1, i.e. flag[j-1]=1 (because the bit wide of flag is since 0, so need to use j-1 to represent the response position of flag), SUM=SUM-(R
j-TH
1); Otherwise flag[j-1]=0, SUM remains unchanged; Read the moment that finishes at first pass, calculate thresholding TH by up-to-date SUM
2, begin RAM1 is read for second time; With the data R that reads
jWith TH
2Relatively, if R
jGreater than thresholding TH
2, check whether the relevant position of flag is 1 this moment, if 1, then the flag relevant position remains unchanged, SUM=SUM-(TH
1-TH
2); Otherwise flag puts the relevant position 1, SUM=SUM-(R
j-TH
2); If R
jLess than thresholding TH
2, then flag corresponding positions and SUM remain unchanged; Read the moment that finishes at second time, calculate thresholding TH according to up-to-date SUM value
3, begin the 3rd time read; The 3rd time to N-1; All over reading the operation done and to read the operation of doing for second time identical, to the last read data one time, namely N can obtain threshold T H all over reading data by N-1 time read operation
N, if this moment is the data R that reads
jGreater than TH
N, then with the real part Re that reads
jBe set to TH
N, imaginary part Im
jBe set to 0 output, i.e. Re
j=TH
N, Im
j=; Otherwise with the real part Re that reads
jWith imaginary part Im
jData are in statu quo exported;
D. when the data decision of RAM1 finished, the data among the RAM2 were also stored end, and begin to store new data to the work of the data repeating step c among the RAM2 to RAM1 this moment;
(5) will disturb the data of identification and the output of inhibition module with clock f
sBe input to the IFFT module, in this module by two-port RAM storage and with clock f
InRead, finish the conversion of data rate.The data input IFFT nuclear that reads carries out inverse transformation and output;
(6) the direct addition of the output of two-way IFFT module is merged, flow to external interface, anti-interference end-of-job.
The invention has the beneficial effects as follows: the FPGA based on the anti-interference algorithm of iteration thresholding that the present invention proposes realizes, can adaptive setting interference detection threshold on algorithm, and have good dynamic disturbance and suppress ability.Realizing technical employing high-frequency clock realization reading data, judgement and inhibition, avoided to take hardware resource to the storage of middle process data because repeatedly adjudicating, taken full advantage of the technology that the FPGA time exchanges resource for, significantly reduced the consumption to hardware resource.
Description of drawings
Fig. 1 is based on the FPGA realization flow figure of the anti-interference algorithm of iteration thresholding;
Fig. 2 disturbs to suppress and the judging module flow chart;
Fig. 3 is frequency spectrum and the correlation schematic diagram that does not add the spread-spectrum signal when disturbing;
Fig. 4 is frequency spectrum and the correlation schematic diagram that adds the spread-spectrum signal after the interference;
Fig. 5 is thresholding iterative process Modelsim analogous diagram;
Fig. 6 is through the frequency spectrum of the spread-spectrum signal after anti-interference and correlation schematic diagram;
Fig. 7 is the RTL figure that FPGA realizes.
Embodiment
The present invention is further described below in conjunction with drawings and Examples.
The present invention proposes a kind of FPGA realization technology based on the anti-interference algorithm of iterative frequency-domain.In the algorithm of adaptive threshold is derived, make improvement, can reduce the consumption of storage resources greatly.In its nucleus module (disturb identification with suppress module), adopt ping-pong operation to receive data, prevent that data from overflowing processing, and use high-frequency clock to finish RAM data read and follow-up thresholding iterative computation and interference decision operation.
(I) first moment iteration threshold algorithm is theoretical derives:
Signal can be expressed as S (k)+N (k) under glitch-free condition, and wherein S (k) represents useful signal, and N (k) represents noise signal.Through becoming S (n)+N (n) after the N point DFT conversion, it is a narrowband Gaussian signal.The envelope of S (n)+N (n) | S (n)+N (n) | Rayleigh distributed, and envelope square | S (n)+N (n) |
2Obeys index distribution.A lot of documents all are that the probability distribution from envelope square begins to derive, and provide the iteration thresholding, and the present invention is from its envelope | S (n)+N (n) | distribution situation derive.
Suppose signal | S (n)+N (n) | obeying parameter is the rayleigh distributed of σ, and its desired value is
Probability density function is
Choose the first moment threshold value
Wherein thresholding is optimized COEFFICIENT K=1,2,3,4,5.Then | S (n)+N (n) | the probability of<TH
Can obtain the probability of different K values, as shown in the table:
The probability distribution situation of table 1 different K values
K=1 | K=2 | K=3 | K=4 | K=5 |
0.5441 | 0.9568 | 0.9991 | 1.0000 | 1.0000 |
The present invention finally chooses K=4, can guarantee that useful signal can all pass through, and the parameter of Xuan Zeing is fit to the hardware realization simultaneously, and multiply operation only can realize with displacement.
When DFT counts more for a long time (N〉256) desired value
Can use the average statistical of signal to substitute, namely
Thresholding is chosen work and is finished.
Adaptive threshold is by repeatedly calculating the average statistical of signal, with the threshold value of adaptively changing correspondence, makes signal pass through the judgement of iteration thresholding, reaches and disturbs the effect that suppresses.Concrete grammar can repeat no more here referring to following performing step.
(II) a kind of FPGA of the anti-interference algorithm based on the iteration thresholding realizes, principle as shown in Figure 1, characterization step is as follows:
(1) data that the AD sampling is obtained are divided into two-way, and one the tunnel carries out 1/2 window length delay (the present invention postpones 256 points), and two paths of data is input to the windowing module respectively.
(2) windowing module is carried out the windowing operation for the data of input, and window function adopts broad sense hamming window, and window length is 512 points, adopts the quantification of 8 bit data bit wides.After the windowing operation is finished, data input FFT module.
(3) use a degree of depth to be L(L=512 in the FFT module) the RAM data of checking input store, read data afterwards and be transported to FFT nuclear, carry out the FFT computing.FFT is set to flowing water type base 2 mode of operations, and operating frequency is f
In, the output of its computing has real part Re and imaginary part Im two parts, among the present invention this two paths of signals and Qi Mo value R is merged into one road signal according to senior middle school's low level, and namely { R} is transported to and disturbs identification and suppress module for Re, Im.
(4) disturbing identification is Key Implementation part of the present invention with suppressing module, its flow process as shown in Figure 2, detailed process is as follows:
A. setting up two degree of depth is for example L=512 of L() two-port RAM nuclear RAM1 and RAM2, data input selected cell is selected RAM1 work, the data of FFT output are with clock f
InStore among the RAM1, the L point that calculates mould value R in institute's deposit data simultaneously adds up and SUM and threshold value TH
1
B. after RAM1 was filled with, data input cell was selected RAM2 work, and the data of FFT output are deposited among the RAM2, and the L point that calculates the deposit data mould value R of institute simultaneously in the repeating step 1 adds up and SUM and threshold value TH
1Work.
C. in the process of RAM2 data storages, to the data among the RAM1 with clock f
sRead.f
sCan determine that for example iterations is 3 times, then f according to required iterations
s=3*f
InCan get, can be to the data read among the RAM1 3 times in this following period of time of RAM2 data storages.A bit wide and RAM1(RAM2 are set) flag register (for example [511:0] flag) that the degree of depth is identical, in to RAM1 first pass process of reading, will represent a part of R of mould value in the data that read
j(j=0,1 ... 511) and TH
1Compare, if greater than TH
1Then the corresponding position of flag register puts 1, i.e. flag[j]=1, SUM=SUM-(R
j-TH
1).Otherwise flag[j]=0, SUM remains unchanged.Read the moment that finishes at first pass, calculate thresholding TH by up-to-date SUM
2, begin RAM1 is read for second time.With the data R that reads
jWith TH
2Relatively, if R
jGreater than thresholding TH
2, check whether the relevant position of flag is 1 this moment, if 1, then the flag relevant position remains unchanged, SUM=SUM-(TH
1-TH
2); Otherwise flag puts the relevant position 1, SUM=SUM-(R
j-TH
2).If R
jLess than thresholding TH
2, then flag corresponding positions and SUM remain unchanged.Read the moment that finishes at second time, calculate thresholding TH according to up-to-date SUM value
3, begin the 3rd time read.If the data R that reads
jGreater than TH
3, then with the real part Re that reads
jBe set to TH
3, imaginary part Im
jBe set to 0 output, i.e. Re
j=TH
3, Im
j=0; Otherwise with the real part Re that reads
jWith imaginary part Im
jData are in statu quo exported.
D. when the data decision of RAM1 finished, the data among the RAM2 were also stored end, and begin to store new data to the work of the data repeating step 3 among the RAM2 to RAM1 this moment.
(5) will disturb the data of identification and the output of inhibition module with clock f
sBe input to the IFFT module, in this module by two-port RAM storage and with clock f
InRead, finish the conversion of data rate.The data input IFFT nuclear that reads carries out inverse transformation and output.
(6) output of two-way IFFT module is merged flow to external interface, anti-interference end-of-job.
Now in conjunction with certain actual wireless band spread receiver the present invention is described further:
The hardware platform that algorithm is realized: the fpga chip of Xilinx x5v1x155 model.
Exploitation and simulated environment: ISE10.1 and Modelsim SE6.5c.
The spread-spectrum signal bandwidth of using is 2MHz, and interference signal is 0.2MHz, and jamming-to-signal ratio is 45dB.
The global clock of receiver system is 32MHz and 96MHz.Specific implementation adopts three iteration.The data bit width of AD input is 4.The IP kernel that carries out FFT and IFFT computing all adopts fixedly bit wide (scaled) pattern, and namely importing bit wide is 8, and the real part of output and imaginary part bit wide are 8.RAM nuclear in the realization all adopts 512 storage depth, and wherein RAM nuclear bit wide is 4 in the FFT module; Disturb to suppress and judging module in the RAM bit wide be 25, the real part data of most-significant byte storage FFT output, in the imaginary data of 8 storage FFT output, store the mould values for low 9; RAM nuclear bit wide in the IFFT module is 20, high 12 storage real part data, least-significant byte storage imaginary data.
The concrete steps that realize are as follows:
(1) data that will obtain the satellite-signal AD sampling that changes to intermediate frequency down are output as two-way, and one the tunnel carries out 1/2 window length delay (the present invention postpones 256 points), and two paths of data is input to the windowing module respectively.
(2) windowing module is carried out the windowing operation for the data of input, and window function adopts broad sense hamming window, and window length is 512 points, adopts the quantification of 8 bit data bit wides.After the windowing operation is finished, data input FFT module.
(3) using a degree of depth in the FFT module is that the data that 512 RAM checks input are stored, and reads data afterwards and is transported to FFT nuclear, carries out the FFT computing.FFT is set to Pipelined, and Streaming I/O pattern, its operating frequency are f
In=32MHz(f
InIdentical with the AD sampling clock), the output of its computing has real part Re and imaginary part Im two parts, among the present invention with these two parts and the mould value that calculated jointly by the real part imaginary part
Merge into one road signal according to senior middle school's low level, namely { R} is transported to and disturbs identification and suppress module for Re, Im.
(4) each disturbs identification as follows with the course of work that suppresses module:
A. set up two degree of depth and be 512 two-port RAM nuclear RAM1 and RAM2, data input selected cell is selected RAM1 work, and the data of FFT output are with clock f
In=32MHz stores among the RAM1, calculate simultaneously that 512 of mould value R in institute's deposit data add up and
And threshold value
B. after RAM1 was filled with, data input selected cell was selected RAM2 work, and the data of FFT output are deposited among the RAM2, and 512 that calculate the deposit data mould value R of institute simultaneously among the repeating step a are added up and SUM and threshold value TH
1Work.
C. in the process of RAM2 data storages, to the data among the RAM1 with clock f
s=3*32=96MHz reads.A bit wide is set is 512 flag register, in to RAM1 first pass process of reading, will represent a part of R of mould value in the data that read
j, j=1,2,3...512 is with TH
1Compare, if greater than TH
1Then the corresponding position of flag register puts 1, i.e. flag[j-1]=1 (because the bit wide of flag is since 0, so need to use j-1 to represent the relevant position of flag bit wide), SUM=SUM-(R
j-TH
1); Otherwise flag[j-1]=0, SUM remains unchanged.Read the moment that finishes at first pass, calculate thresholding TH by up-to-date SUM
2(computing formula is shown in step a) begins RAM1 is read for second time.With the data R that reads
jWith TH
2Relatively, if R
jGreater than thresholding TH
2, check whether the relevant position of flag is 1 this moment, if 1, then the flag relevant position remains unchanged, SUM=SUM-(TH
1-TH
2); Otherwise flag puts the relevant position 1, SUM=SUM-(R
j-TH
2); If R
jLess than thresholding TH
2, then flag corresponding positions and SUM remain unchanged.Read the moment that finishes at second time, calculate thresholding TH according to up-to-date SUM value
3, begin the 3rd time read.If this moment is the data R that reads
jGreater than TH
3, then with the real part Re that reads
jBe set to TH
3, imaginary part Im
jBe set to 0 output, i.e. Re
j=TH
3, Im
j=0; Otherwise with the real part Re that reads
jWith imaginary part Im
jData are in statu quo exported
D. when the data decision of RAM1 finished, the data among the RAM2 were also stored end, and begin to store new data to the work of the data repeating step c among the RAM2 to RAM1 this moment.
(5) data of the inhibition that will disturb and judging module output are with clock f
s=96MHz is input to the IFFT module, in this module by two-port RAM storage and with clock f
In=32MHz reduction of speed rate reads, and finishes the conversion of data rate.The data input IFFT nuclear that reads carries out inverse transformation and output.
(6) the direct addition of the output of two-way IFFT module is merged, flow to external interface, anti-interference end-of-job.
According to the characteristic of spread-spectrum signal, can judge the degree that interference signal is suppressed by the correlation peak of spread-spectrum signal and local code signal.
The frequency spectrum of the spread-spectrum signal when not adding interference and correlation are as shown in Figure 3.
Add the frequency spectrum of spread-spectrum signal after the interference signal that bandwidth 2MHz jamming-to-signal ratio is 40dB and correlation as shown in Figure 4.
The Modelsim analogous diagram of thresholding iterative process as shown in Figure 5 when FPGA realized.
Threshold1 in the signal shown in Figure 5, threshold2, the respectively corresponding iteration of the threshold3 first time, the second time and threshold value for the third time.
From overall procedure, after every judgement through 512 point data, each threshold value (as threshold1) can obtain upgrading according to the statistics of the data of new input, and namely threshold value can be according to according to the difference of input data and dynamic change.Once complete judging process from 512 point data, through declaring for the first time iteration, threshold value drops to threshold2 by threshold1, through the iteration second time, threshold value drops to threshold3 by threshold2, in the process of iteration for the third time, with data decision output, finish the interference judgement work of 512 point data.Specific implementation conforms to theory.
After the interference that signal is realized by FPGA suppressed module, the frequency spectrum of spread-spectrum signal and correlation were as shown in Figure 6.
According to Fig. 4, Fig. 5 and result shown in Figure 6, can not get tangible correlation peak after adding interference signal as can be seen, receiver can't be worked.Through after disturb suppressing, correlation peak is apparent in view, and the position is with not add interference signal identical, and interference suppresses successfully.
After the comprehensive and placement-and-routing, the RTL that FPGA realizes schemes as shown in Figure 7 through ISE.
Shown in the concrete resource consumption information slip 2.
Table 2 resource consumption figure
Claims (1)
1. the FPGA implementation method of the anti-interference algorithm of iterative frequency-domain is characterized in that comprising the steps:
(1) will carry out the AD sampling to the satellite-signal that changes to intermediate frequency down, the data that obtain are output as two-way, and one the tunnel carries out 1/2 window length delay, and two paths of data is input to the windowing module respectively;
(2) windowing module is carried out the windowing operation for the data of input, and window function adopts broad sense hamming window, and window length is L, adopts the quantification of 8 bit data bit wides; After the windowing operation is finished, data input FFT module;
(3) use a degree of depth to store as the data that the RAM of L checks input in the FFT module, read data afterwards and be transported to FFT nuclear, carry out the FFT computing; FFT is set to Pipelined, and Streaming I/O pattern, operating frequency are f
In, f
InIdentical with the AD sampling clock, the output of its computing has real part Re and imaginary part Im two parts, merges into one road signal with these two parts with by the mould value R that the real part imaginary part calculates jointly according to senior middle school's low level among the present invention, i.e. { Re, Im, R} is transported to and disturbs identification and suppress module;
(4) each disturbs identification as follows with the course of work that suppresses module:
A. two-port RAM nuclear RAM1 and RAM2 that to set up two degree of depth be L=512, data input selected cell is selected RAM1 work, and the data of FFT output are with clock f
InStore among the RAM1, the L point that calculates mould value R in institute's deposit data simultaneously add up and
And threshold value
K is that thresholding is optimized coefficient;
B. after RAM1 was filled with, data input selected cell was selected RAM2 work, and the data of FFT output are deposited among the RAM2, and the L point that calculates the deposit data mould value R of institute simultaneously among the repeating step a adds up and SUM and threshold value TH
1Work;
C. in the process of RAM2 data storages, to the data among the RAM1 with clock f
sRead f
s=N*f
In, N is iterations; Can get, can be to the data read among the RAM1 3 times in this following period of time of RAM2 data storages; A flag register that bit wide is identical with the RAM1 degree of depth is set, in to RAM1 first pass process of reading, will represents a part of R of mould value in the data that read
jWith TH
1Compare, j=1, if 2,3...L is greater than TH
1Then the corresponding position of flag register puts 1, i.e. flag[j-1]=1, SUM=SUM-(R
j-TH
1); Otherwise flag[j-1]=0, SUM remains unchanged; Read the moment that finishes at first pass, calculate thresholding TH by up-to-date SUM
2, begin RAM1 is read for second time; With the data R that reads
jWith TH
2Relatively, if R
jGreater than thresholding TH
2, check whether the relevant position of flag is 1 this moment, if 1, then the flag relevant position remains unchanged, SUM=SUM-(TH
1-TH
2); Otherwise flag puts the relevant position 1, SUM=SUM-(R
j-TH
2); If R
jLess than thresholding TH
2, then flag corresponding positions and SUM remain unchanged; Read the moment that finishes at second time, calculate thresholding TH according to up-to-date SUM value
3, begin the 3rd time read; The 3rd time to N-1; All over reading the operation done and to read the operation of doing for second time identical, to the last read data one time, namely N can obtain threshold T H all over reading data by N-1 time read operation
N, if this moment is the data R that reads
jGreater than TH
N, then with the real part Re that reads
jBe set to TH
N, imaginary part Im
jBe set to 0 output, i.e. Re
j=TH
N, Im
j=0; Otherwise with the real part Re that reads
jWith imaginary part Im
jData are in statu quo exported;
D. when the data decision of RAM1 finished, the data among the RAM2 were also stored end, and begin to store new data to the work of the data repeating step c among the RAM2 to RAM1 this moment;
(5) will disturb the data of identification and the output of inhibition module with clock f
sBe input to the IFFT module, in this module by two-port RAM storage and with clock f
InRead, finish the conversion of data rate; The data input IFFT nuclear that reads carries out inverse transformation and output;
(6) the direct addition of the output of two-way IFFT module is merged, flow to external interface, anti-interference end-of-job.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310077058.1A CN103199889B (en) | 2013-03-11 | 2013-03-11 | Field programmable gata array (FPGA) implementation method of iteration frequency domain anti-interference algorithm |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310077058.1A CN103199889B (en) | 2013-03-11 | 2013-03-11 | Field programmable gata array (FPGA) implementation method of iteration frequency domain anti-interference algorithm |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103199889A true CN103199889A (en) | 2013-07-10 |
CN103199889B CN103199889B (en) | 2014-12-10 |
Family
ID=48722271
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310077058.1A Active CN103199889B (en) | 2013-03-11 | 2013-03-11 | Field programmable gata array (FPGA) implementation method of iteration frequency domain anti-interference algorithm |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103199889B (en) |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103926595A (en) * | 2014-04-22 | 2014-07-16 | 中国电子科技集团公司第二十研究所 | Interference signal detection method based on small satellite navigation terminal |
CN104076369A (en) * | 2014-07-10 | 2014-10-01 | 西北工业大学 | Frequency domain anti-interference method and device based on adaptive threshold judgment |
CN104716982A (en) * | 2015-03-25 | 2015-06-17 | 北京北斗星通导航技术股份有限公司 | Robust anti-interference processing method and device of spread frequency system |
CN104777491A (en) * | 2015-04-22 | 2015-07-15 | 北京北斗星通导航技术股份有限公司 | Blind beam broadband interference suppression method and device |
CN109815877A (en) * | 2019-01-17 | 2019-05-28 | 北京邮电大学 | A kind of method and device for noise reduction of satellite-signal |
CN111585593A (en) * | 2020-03-25 | 2020-08-25 | 北京理工大学 | Ultra-wideband signal interference suppression method and system |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101154983A (en) * | 2006-09-27 | 2008-04-02 | 上海微小卫星工程中心 | Implementing method for anti-interference technology of satellite communication based on single-time sampling |
CN101212235A (en) * | 2007-12-24 | 2008-07-02 | 北京邮电大学 | Low-complexity frequency domain implementation method for combined frequency domain equalization interference elimination detection algorithm |
US7720134B2 (en) * | 2004-05-10 | 2010-05-18 | Stmicroelectronics S.R.L. | Frequency-domain multi-user access interference cancellation and nonlinear equalization in CDMA receivers |
-
2013
- 2013-03-11 CN CN201310077058.1A patent/CN103199889B/en active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7720134B2 (en) * | 2004-05-10 | 2010-05-18 | Stmicroelectronics S.R.L. | Frequency-domain multi-user access interference cancellation and nonlinear equalization in CDMA receivers |
CN101154983A (en) * | 2006-09-27 | 2008-04-02 | 上海微小卫星工程中心 | Implementing method for anti-interference technology of satellite communication based on single-time sampling |
CN101212235A (en) * | 2007-12-24 | 2008-07-02 | 北京邮电大学 | Low-complexity frequency domain implementation method for combined frequency domain equalization interference elimination detection algorithm |
Non-Patent Citations (5)
Title |
---|
CATHARINA CARLEMALM, H. VINCENT POOR,ANDREW LOGOTHETIS: "《Suppression of Multiple Narrowband Interferers in a Spread-Spectrum Communication System》", 《IEEE JOURNAL ON SELECTED AREAS IN COMMUNICATIONS》 * |
张爱民,胡艳龙,韩方景: "《低信噪比下基于自适应门限的窄带干扰抑制研究》", 《电子信息对抗技术》 * |
沈志伟: "《导航接收机迭代处理抗干扰技术》", 《电子科技大学硕士学位论文》 * |
王耿: "《无线通信中迭代干扰抑制技术研究》", 《电子科技大学硕士学位论文》 * |
邹宁,徐松涛,刘明园,刘凯: "《一种重叠加窗频域抑制窄带干扰算法及研究》", 《现代防御技术》 * |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103926595A (en) * | 2014-04-22 | 2014-07-16 | 中国电子科技集团公司第二十研究所 | Interference signal detection method based on small satellite navigation terminal |
CN104076369A (en) * | 2014-07-10 | 2014-10-01 | 西北工业大学 | Frequency domain anti-interference method and device based on adaptive threshold judgment |
CN104716982A (en) * | 2015-03-25 | 2015-06-17 | 北京北斗星通导航技术股份有限公司 | Robust anti-interference processing method and device of spread frequency system |
CN104777491A (en) * | 2015-04-22 | 2015-07-15 | 北京北斗星通导航技术股份有限公司 | Blind beam broadband interference suppression method and device |
CN109815877A (en) * | 2019-01-17 | 2019-05-28 | 北京邮电大学 | A kind of method and device for noise reduction of satellite-signal |
CN111585593A (en) * | 2020-03-25 | 2020-08-25 | 北京理工大学 | Ultra-wideband signal interference suppression method and system |
Also Published As
Publication number | Publication date |
---|---|
CN103199889B (en) | 2014-12-10 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103199889B (en) | Field programmable gata array (FPGA) implementation method of iteration frequency domain anti-interference algorithm | |
CN102664657B (en) | Method for self-adaption mid-value threshold frequency domain anti-jamming | |
CN102420671B (en) | Self-adaptive Turbo decoding iterations selection method | |
CN104007453A (en) | Frequency domain and space domain joint anti-interference method aided by probability search | |
CN110932807B (en) | Spectrum sensing method and system for MIMO (multiple input multiple output) system under non-Gaussian noise | |
CN107425929B (en) | Non-auxiliary data equalization method for fading channel under Alpha stable distributed noise | |
CN103926595A (en) | Interference signal detection method based on small satellite navigation terminal | |
CN105635009A (en) | Self-adaptive MIMO pre-distortion method for hybrid compensation of multi-branch crosstalk and IQ imbalance | |
Sui et al. | Jointly optimized extreme learning machine for short-term prediction of fading channel | |
CN112014801A (en) | Composite interference identification method based on SPWVD and improved AlexNet | |
Alam et al. | Performance comparison of STFT, WT, LMS and RLS adaptive algorithms in denoising of speech signal | |
CN101711049A (en) | Routing method and device based on interference elimination | |
Li et al. | Lightweight channel estimation networks for OFDM systems | |
Liu et al. | Complex adaptive LMS algorithm employing the conjugate gradient principle for channel estimation and equalization | |
CN102045290B (en) | Gray modeling-based OFDM narrow-band slow-fading slowly time-varying channel estimation method | |
CN110248325B (en) | Bluetooth indoor positioning system based on signal multiple noise elimination | |
Dai et al. | NLMS adaptive algorithm implement based on FPGA | |
CN101651900B (en) | Channel estimation method and device for eliminating interference within frequency domain | |
Liao et al. | Determining neighborhoods of image pixels automatically for adaptive image denoising using nonlinear time series analysis | |
CN110516566B (en) | Filtering method and device based on convolutional layer | |
Zhu et al. | A novel wavelet denoising pre-processing algorithm for TDOA localization | |
Cai et al. | Application of three-threshold FCME and extended interpolation algorithm in narrowband interference suppression | |
CN106559364A (en) | A kind of iterative channel estimation method and device | |
CN102185808B (en) | Rapidly-convergent immune-clone-based orthogonal wavelet transform constant modulus blind equalization algorithm | |
Yao et al. | Iterative threshold algorithm for narrow-band interference suppression |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |