CN111934651A - Space-time frequency cascade adaptive filtering processing method, device and equipment - Google Patents

Space-time frequency cascade adaptive filtering processing method, device and equipment Download PDF

Info

Publication number
CN111934651A
CN111934651A CN202010552672.9A CN202010552672A CN111934651A CN 111934651 A CN111934651 A CN 111934651A CN 202010552672 A CN202010552672 A CN 202010552672A CN 111934651 A CN111934651 A CN 111934651A
Authority
CN
China
Prior art keywords
reference channel
domain data
channel
frequency
frequency domain
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Pending
Application number
CN202010552672.9A
Other languages
Chinese (zh)
Inventor
任超
屠恩源
吴红甲
郑永翔
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
BEIJING BDSTAR NAVIGATION Co.,Ltd.
Original Assignee
Unicorecomm Shanghai Technology Co ltd
Hexin Xingtong Technology Beijing Co ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Unicorecomm Shanghai Technology Co ltd, Hexin Xingtong Technology Beijing Co ltd filed Critical Unicorecomm Shanghai Technology Co ltd
Priority to CN202010552672.9A priority Critical patent/CN111934651A/en
Publication of CN111934651A publication Critical patent/CN111934651A/en
Pending legal-status Critical Current

Links

Images

Classifications

    • H—ELECTRICITY
    • H03—ELECTRONIC CIRCUITRY
    • H03H—IMPEDANCE NETWORKS, e.g. RESONANT CIRCUITS; RESONATORS
    • H03H21/00—Adaptive networks
    • H03H21/0012—Digital adaptive filters
    • H03H21/0043—Adaptive algorithms
    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04B—TRANSMISSION
    • H04B1/00—Details of transmission systems, not covered by a single one of groups H04B3/00 - H04B13/00; Details of transmission systems not characterised by the medium used for transmission
    • H04B1/06—Receivers
    • H04B1/10—Means associated with receiver for limiting or suppressing noise or interference
    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04B—TRANSMISSION
    • H04B1/00—Details of transmission systems, not covered by a single one of groups H04B3/00 - H04B13/00; Details of transmission systems not characterised by the medium used for transmission
    • H04B1/06—Receivers
    • H04B1/10—Means associated with receiver for limiting or suppressing noise or interference
    • H04B1/12—Neutralising, balancing, or compensation arrangements
    • H04B1/123—Neutralising, balancing, or compensation arrangements using adaptive balancing or compensation means
    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04B—TRANSMISSION
    • H04B17/00—Monitoring; Testing
    • H04B17/20—Monitoring; Testing of receivers
    • H04B17/21—Monitoring; Testing of receivers for calibration; for correcting measurements

Abstract

A space-time-frequency cascade adaptive filtering processing method, a device, equipment and a computer-readable storage medium are provided, wherein the method comprises the following steps: carrying out frequency conversion amplification and analog-to-digital conversion on an input signal to obtain a plurality of paths of digital intermediate frequency signals; selecting a reference channel from channels corresponding to the multi-channel digital intermediate frequency signals; converting the multi-channel digital intermediate frequency signals into frequency domain data through Fourier transformation, and determining an amplitude-phase difference correction coefficient of each non-reference channel based on the frequency domain data of the reference channel; correcting the frequency domain data of each non-reference channel according to the amplitude difference correction coefficient; and carrying out inverse Fourier transform on the corrected frequency domain data of the non-reference channel to obtain time domain data. The embodiment of the application carries out frequency domain amplitude-phase error correction based on the reference channel, can improve the broadband interference suppression capability and simultaneously enhances the robustness of self-adaptive filtering.

Description

Space-time frequency cascade adaptive filtering processing method, device and equipment
Technical Field
The present disclosure relates to the field of communications, and in particular, to a method, an apparatus, a device, and a computer-readable storage medium for space-time-frequency cascaded adaptive filtering.
Background
In order to suppress forced pressure type interference in a complex electromagnetic environment, the array antenna adaptive filtering technology is widely applied to a spread spectrum communication system.
However, in practical engineering application, the anti-interference performance of the spatial filter is seriously affected by the amplitude-phase error of the array channel, which may cause the failure of the filter.
Disclosure of Invention
The application provides a space-time-frequency cascading adaptive filtering processing method, a device, equipment and a computer readable storage medium, so as to increase the robustness of filtering.
The embodiment of the application provides a space-time frequency cascade adaptive filtering processing method, which comprises the following steps:
carrying out frequency conversion amplification and analog-to-digital conversion on an input signal to obtain a plurality of paths of digital intermediate frequency signals;
selecting a reference channel from channels corresponding to the multi-channel digital intermediate frequency signals;
converting the multi-channel digital intermediate frequency signals into frequency domain data through Fourier transformation, and determining an amplitude-phase difference correction coefficient of each non-reference channel based on the frequency domain data of the reference channel;
correcting the frequency domain data of each non-reference channel according to the amplitude difference correction coefficient;
and carrying out inverse Fourier transform on the corrected frequency domain data of the non-reference channel to obtain time domain data.
The embodiment of the present application further provides a space-time frequency cascade adaptive filtering processing apparatus, including:
the first signal conversion module is used for carrying out frequency conversion amplification and analog-to-digital conversion on the input signal to obtain a plurality of paths of digital intermediate frequency signals;
the reference channel selection module is used for selecting a reference channel from the channels corresponding to the multi-channel digital intermediate frequency signals;
the second signal conversion module is used for converting the multi-channel digital intermediate frequency signals into frequency domain data through Fourier transform, and determining the amplitude-phase difference correction coefficient of each non-reference channel based on the frequency domain data of the reference channel;
the correction module is used for correcting the frequency domain data of each non-reference channel according to the amplitude-phase difference correction coefficient;
and the third signal conversion module is used for performing inverse Fourier transform on the corrected frequency domain data of the non-reference channel to obtain time domain data.
The embodiment of the present application further provides a space-time frequency cascade adaptive filtering processing apparatus, including: the space-time-frequency cascade adaptive filtering processing method comprises a memory, a processor and a computer program which is stored on the memory and can run on the processor, wherein the processor realizes the space-time-frequency cascade adaptive filtering processing method when executing the program.
The embodiment of the application also provides a computer-readable storage medium, which stores computer-executable instructions, and the computer-executable instructions are used for executing the space-time-frequency cascading adaptive filtering processing method.
Compared with the related art, the embodiment of the application comprises the following steps: carrying out frequency conversion amplification and analog-to-digital conversion on an input signal to obtain a plurality of paths of digital intermediate frequency signals; selecting a reference channel from channels corresponding to the multi-channel digital intermediate frequency signals; converting the multi-channel digital intermediate frequency signals into frequency domain data through Fourier transformation, and determining an amplitude-phase difference correction coefficient of each non-reference channel based on the frequency domain data of the reference channel; correcting the frequency domain data of each non-reference channel according to the amplitude difference correction coefficient; and carrying out inverse Fourier transform on the corrected frequency domain data of the non-reference channel to obtain time domain data. The embodiment of the application carries out frequency domain amplitude-phase error correction based on the reference channel, can improve the broadband interference suppression capability and simultaneously enhances the robustness of self-adaptive filtering.
Additional features and advantages of the application will be set forth in the description which follows, and in part will be obvious from the description, or may be learned by the practice of the application. Other advantages of the present application may be realized and attained by the instrumentalities and combinations particularly pointed out in the specification and the drawings.
Drawings
The accompanying drawings are included to provide an understanding of the present disclosure and are incorporated in and constitute a part of this specification, illustrate embodiments of the disclosure and together with the examples serve to explain the principles of the disclosure and not to limit the disclosure.
Fig. 1 is a flowchart of a space-time frequency cascade adaptive filtering processing method according to an embodiment of the present application;
fig. 2 is a schematic diagram of a space-time-frequency cascaded adaptive filtering processing structure according to an application example of the present application;
fig. 3 is a schematic diagram of space-time frequency cascade adaptive filtering processing according to an embodiment of the present application.
Detailed Description
The present application describes embodiments, but the description is illustrative rather than limiting and it will be apparent to those of ordinary skill in the art that many more embodiments and implementations are possible within the scope of the embodiments described herein. Although many possible combinations of features are shown in the drawings and discussed in the detailed description, many other combinations of the disclosed features are possible. Any feature or element of any embodiment may be used in combination with or instead of any other feature or element in any other embodiment, unless expressly limited otherwise.
The present application includes and contemplates combinations of features and elements known to those of ordinary skill in the art. The embodiments, features and elements disclosed in this application may also be combined with any conventional features or elements to form a unique inventive concept as defined by the claims. Any feature or element of any embodiment may also be combined with features or elements from other inventive aspects to form yet another unique inventive aspect, as defined by the claims. Thus, it should be understood that any of the features shown and/or discussed in this application may be implemented alone or in any suitable combination. Accordingly, the embodiments are not limited except as by the appended claims and their equivalents. Furthermore, various modifications and changes may be made within the scope of the appended claims.
Further, in describing representative embodiments, the specification may have presented the method and/or process as a particular sequence of steps. However, to the extent that the method or process does not rely on the particular order of steps set forth herein, the method or process should not be limited to the particular sequence of steps described. Other orders of steps are possible as will be understood by those of ordinary skill in the art. Therefore, the particular order of the steps set forth in the specification should not be construed as limitations on the claims. Further, the claims directed to the method and/or process should not be limited to the performance of their steps in the order written, and one skilled in the art can readily appreciate that the sequences may be varied and still remain within the spirit and scope of the embodiments of the present application.
As shown in fig. 1, an embodiment of the present application provides a space-time frequency cascaded adaptive filtering processing method, including:
step 101, performing frequency conversion amplification and analog-to-digital conversion on an input signal to obtain a plurality of paths of digital intermediate frequency signals.
Wherein the input signal may comprise a navigation signal.
Setting a receiving antenna array to have M array elements, receiving an input signal through the M array elements, performing frequency conversion amplification processing through a radio frequency circuit, and setting the transmission function of the M channels as H due to different amplitude-frequency and phase-frequency responses of each radio frequency channelm(ω), M ═ 1, 2. The output M paths of analog intermediate frequency signals pass through AD (analog revolution)Word) sampling to obtain M digital intermediate frequency signals x1(h),...,xM(h) And h represents a discrete point number of time. The array input vector is x ═ x1,...,xM]TThe weight vector is b ═ b1,...bM]T。
And 102, selecting a reference channel from the channels corresponding to the multipath digital intermediate frequency signals.
In an embodiment, the output power of each channel when the channel is the main channel is calculated, and the channel corresponding to the minimum output power is taken as the reference channel.
In one embodiment, the constraint vector when m channel is the main channel is denoted as sm＝[0,...,1,...,0]TWherein the m-channel correspondence is position set to 1. The output power when each channel is the main channel is calculated respectively,
Figure BDA0002543073510000041
selecting a minimum output power PmThe corresponding channel m is the reference channel. Where μ is the constant coefficient, RxxA data correlation matrix is received for a channel.
And 103, converting the multi-channel digital intermediate frequency signals into frequency domain data through Fourier transform, and determining the amplitude-phase difference correction coefficient of each non-reference channel based on the frequency domain data of the reference channel.
The fourier transform may be a Fast Fourier Transform (FFT) or a Discrete Fourier Transform (DFT).
In this step, M channels of digital intermediate frequency signals are subjected to K-point fourier transform, and each channel of digital intermediate frequency signal is divided into K narrow sub-bands in the frequency domain. Obtaining a frequency domain signal xmk(l) The corresponding frequency point f of the mth array element digital intermediate frequency data after FFT is shownkThe ith fast beat number of (c).
In an embodiment, the amplitude ratio α of the frequency domain data of the reference channel to the frequency domain data of the non-reference channel m is determinedmkAnd a phase difference of betamkThe amplitude-phase difference correction coefficient w is determined bymk：
Figure BDA0002543073510000051
Where k represents a frequency point. Let the reference channel be channel a, xakIs frequency domain data of a reference channel, xmkFrequency domain data for non-reference channels.
And step 104, correcting the frequency domain data of each non-reference channel according to the amplitude difference correction coefficient.
Wherein, x is passedmkwmkThe frequency domain data for each non-reference channel is corrected separately.
And 105, performing inverse Fourier transform on the corrected frequency domain data of the non-reference channel to obtain time domain data.
Wherein, the channel corrected data xmkwmkAnd performing Inverse Fast Fourier Transform (IFFT) at K points to obtain time domain data.
In an embodiment, after step 105, the method further includes:
and performing FIR (Finite Impulse Response) filtering on the time domain data of the non-reference channel, and taking the time domain data obtained after filtering and the time domain data of the reference channel as output data.
In an embodiment, the time domain data of the non-reference channel is input into an N-order FIR filter bank, and filter coefficient vectors adopted by the N-order FIR filter bank are:
Figure BDA0002543073510000052
wherein the content of the first and second substances,
Figure BDA0002543073510000053
e represents a weighted average, yaTime domain data, Y, representing the reference channel aMRepresenting the input signals of all taps of the FIR filter bank.
Wherein, the time domain data y of the reference channel aaThe digital intermediate frequency signal of the reference channel obtained in the step 101 to 102。
The embodiment of the application carries out frequency domain amplitude-phase error correction based on the reference channel, can improve the broadband interference suppression capability and simultaneously enhances the robustness of self-adaptive filtering.
Referring to fig. 2, the input signal is taken as a navigation signal, and the reference channel is taken as a first channel (channel 1) for explanation. In the figure, F denotes a frequency domain interval and T denotes a time domain interval.
The receiving antenna array has M array elements, the navigation signal is received by M array elements and is processed by frequency conversion and amplification of the radio frequency circuit, and the transmission function of M channels is set as H due to different amplitude-frequency and phase-frequency responses of each radio frequency channelm(ω)，m＝1,2,...,M。
The M paths of output analog intermediate frequency signals are subjected to AD sampling to obtain M paths of digital intermediate frequency signals x1(h),...,xM(h) And h represents a discrete point number of time. Each path of digital intermediate frequency enters the reference channel selection processing.
When the array input vector is x ═ x1,...,xM]TThe weight vector is b ═ b1,...bM]TThe array output power is:
Pout＝E{(bHx)(bHx)*}＝bHRxxb
since the practical situation is that the spatial interference is uncorrelated with the noise, and the interference power is much larger than the noise power, the correlation matrix R of the received dataxxCan be expressed as follows:
Figure BDA0002543073510000061
where xi is 1/2 pi B, sigmajFor interference power, B for reception bandwidth,
Figure BDA0002543073510000062
under the condition that the array weight vector b is normalized, the smaller the array output power is, the stronger the anti-interference capability is, so that the channel with the minimum output power is selected as a reference channel.
Let the constraint vector s take the value of channel m as the coefficient of 1, and the other coefficients are all 0, representing the output of channel m, and recording as sm＝[0,...,1,...,0]T. The output power when the channel m is the main channel is
Pm＝bHRxxb
The spatial filter is based on the minimum output power to obtain
Figure BDA0002543073510000063
Where μ is a constant coefficient. The compound is substituted into the formula to obtain the compound,
Figure BDA0002543073510000064
selecting a minimum output power PmThe corresponding channel is the reference channel. As shown in fig. 2, let the reference channel be channel 1.
The spatial filter expects the amplitude of each channel at the AD front end to be consistent, i.e. the transmission function is HmAnd (ω), M is 1,2, and M is the same, and in actual conditions, the amplitude-frequency response and the phase-frequency response of each channel are different greatly, so that the channels are subjected to amplitude phase correction in a digital domain.
And respectively carrying out K-point Fast Fourier Transform (FFT) on the M paths of digital intermediate frequency signals, and dividing each path of digital intermediate frequency signal into K narrow sub-bands on a frequency domain. Obtaining a frequency domain signal xmk(l) The corresponding frequency point f of the mth array element digital intermediate frequency data after FFT is shownkThe ith fast beat number of (c).
The same broadband signal passes through the transfer function H simultaneouslymAfter (ω), M is 1,2, M, the frequency domain signal x of each channel can be corrected bymk(l) To achieve correction of the channel.
In order to obtain the degree of inconsistency between channels, the data of the reference channel is taken as a reference signal, and then the corresponding frequency point fkThe amplitude-phase difference correction coefficient between channel m and channel 1 at
Figure BDA0002543073510000071
Wherein m ═2,...,M，k＝1,2,...,K，αmkIs the amplitude ratio, betamkIs the phase difference.
Second to M channel intermediate frequency signal x2(h),...,xM(h) After being respectively corrected by the amplitude-phase correction coefficient, the data are processed by K-point Inverse Fast Fourier Transform (IFFT) to obtain time domain data output which is recorded as y2(h),...,yM(h)。
y1(h) Is equal to x1(h) In that respect In addition to the reference channel, y2(n),...,yM(N) each path of signal is designed with an N-order FIR filter, and then each tap input signal of FIR behind array element m is ym1(h)＝ym(h)，ym2(h)＝ym(h-1)，……，ymN(h)＝ym(h-N + 1). And performing interference cancellation processing on the reference channel through M-1 FIR filter groups. Let the coefficients of the FIR filter bank be denoted qmnN, the coefficient vector is written as:
Q＝[q21,...,q2N,...,qM1,...,qMN]T
the optimization criterion can be summarized as the minimum power after cancellation processing of the reference channel:
Figure BDA0002543073510000072
in the formula
Figure BDA0002543073510000073
YM＝[y2(h),y2(h-1),...,y2(h-N+1),...,yM(h-N+1)]T。
PoutThe optimal weight Q for the minimum value can be solved for the gradient to be zero:
Figure BDA0002543073510000074
when in use
Figure BDA0002543073510000075
When the order is full, the only solution of normal equation can be solved
Figure BDA0002543073510000076
Therefore, multi-channel adaptive filtering is realized, and the final output of the space-time frequency cascade adaptive filter is as follows: y (n) ═ y1(h)+q21y2(h)+...+q2Ny2(h-N+1)+...+qM1yM(h)...+qMNyM(h-N+1)。
As shown in fig. 3, an embodiment of the present application further provides a space-time frequency cascade adaptive filtering processing apparatus, including:
the first signal conversion module 21 is configured to perform frequency conversion amplification and analog-to-digital conversion on an input signal to obtain a plurality of paths of digital intermediate frequency signals;
a reference channel selection module 22, configured to select a reference channel from channels corresponding to the multiple paths of digital intermediate frequency signals;
the second signal conversion module 23 is configured to perform fourier transform on the multiple paths of digital intermediate frequency signals to convert the multiple paths of digital intermediate frequency signals into frequency domain data, and determine an amplitude-phase difference correction coefficient of each non-reference channel based on the frequency domain data of the reference channel;
the correcting module 24 is configured to correct the frequency domain data of each non-reference channel according to the amplitude-phase difference correction coefficient;
and a third signal conversion module 25, configured to perform inverse fourier transform on the frequency domain data of the corrected non-reference channel to obtain time domain data.
In one embodiment, the reference channel selection module 22 is configured to:
and calculating the output power when each channel is a main channel, and taking the channel corresponding to the minimum output power as a reference channel.
In an embodiment, the reference channel selecting module 22 is configured to calculate the output power P when the m channel is the main channel according to the following formulam：
Figure BDA0002543073510000081
Where μ is a constant coefficient, and m is a constraint vector in the case of a main channel, RxxReceiving a data correlation matrix, s, for a channelmIs a constraint vector.
In an embodiment, the second signal conversion module 23 is configured to:
determining an amplitude ratio alpha of the frequency domain data of the reference channel to the frequency domain data of the non-reference channel mmkAnd a phase difference of betamkThe amplitude-phase difference correction coefficient w is determined bymk：
wmk＝αmkexp(jβmk)
Where k represents a frequency point.
In one embodiment, the apparatus further comprises:
and the output module is used for performing FIR filtering on the time domain data of the non-reference channel and taking the time domain data obtained after filtering and the time domain data of the reference channel as output data.
In one embodiment, the output module is configured to:
inputting the time domain data of the non-reference channel into an N-order FIR filter group, wherein the N-order FIR filter group adopts a filter coefficient vector as follows:
Figure BDA0002543073510000091
wherein the content of the first and second substances,
Figure BDA0002543073510000092
e represents a weighted average, yaTime domain data representing said reference channel, YMRepresenting the input signals of all taps of the FIR filter bank.
The embodiment of the application provides a space-time frequency cascade adaptive filtering processing method and device. Because the transmission functions of all channels are different, a reference channel is selected through data preprocessing, a frequency domain amplitude-phase error correction method based on the reference channel is provided, and interference cancellation self-adaptive filtering design of the reference channel is carried out on the premise of improving the consistency of a space-time filter channel model. The space-time frequency cascade adaptive filtering provided by the embodiment of the application effectively solves the problem of failure of a broadband beam former, improves the broadband interference suppression capability, and simultaneously enhances the robustness of an adaptive filtering algorithm.
The embodiment of the present application further provides a space-time frequency cascade adaptive filtering processing apparatus, including: the space-time-frequency cascade adaptive filtering processing method comprises a memory, a processor and a computer program which is stored on the memory and can run on the processor, wherein the processor realizes the space-time-frequency cascade adaptive filtering processing method when executing the program.
The embodiment of the present application further provides a computer-readable storage medium, which stores computer-executable instructions, where the computer-executable instructions are used to execute the space-time frequency cascade adaptive filtering processing method.
In this embodiment, the storage medium may include, but is not limited to: a U-disk, a Read-Only Memory (ROM), a Random Access Memory (RAM), a removable hard disk, a magnetic or optical disk, and other various media capable of storing program codes.
It will be understood by those of ordinary skill in the art that all or some of the steps of the methods, systems, functional modules/units in the devices disclosed above may be implemented as software, firmware, hardware, and suitable combinations thereof. In a hardware implementation, the division between functional modules/units mentioned in the above description does not necessarily correspond to the division of physical components; for example, one physical component may have multiple functions, or one function or step may be performed by several physical components in cooperation. Some or all of the components may be implemented as software executed by a processor, such as a digital signal processor or microprocessor, or as hardware, or as an integrated circuit, such as an application specific integrated circuit. Such software may be distributed on computer readable media, which may include computer storage media (or non-transitory media) and communication media (or transitory media). The term computer storage media includes volatile and nonvolatile, removable and non-removable media implemented in any method or technology for storage of information such as computer readable instructions, data structures, program modules or other data, as is well known to those of ordinary skill in the art. Computer storage media includes, but is not limited to, RAM, ROM, EEPROM, flash memory or other memory technology, CD-ROM, Digital Versatile Disks (DVD) or other optical disk storage, magnetic cassettes, magnetic tape, magnetic disk storage or other magnetic storage devices, or any other medium which can be used to store the desired information and which can accessed by a computer. In addition, communication media typically embodies computer readable instructions, data structures, program modules or other data in a modulated data signal such as a carrier wave or other transport mechanism and includes any information delivery media as known to those skilled in the art.

Claims (10)

1. A space-time frequency cascade adaptive filtering processing method comprises the following steps:
carrying out frequency conversion amplification and analog-to-digital conversion on an input signal to obtain a plurality of paths of digital intermediate frequency signals;
selecting a reference channel from channels corresponding to the multi-channel digital intermediate frequency signals;
converting the multi-channel digital intermediate frequency signals into frequency domain data through Fourier transformation, and determining an amplitude-phase difference correction coefficient of each non-reference channel based on the frequency domain data of the reference channel;
correcting the frequency domain data of each non-reference channel according to the amplitude difference correction coefficient;
and carrying out inverse Fourier transform on the corrected frequency domain data of the non-reference channel to obtain time domain data.
2. The method according to claim 1, wherein the selecting a reference channel from the channels corresponding to the plurality of digital intermediate frequency signals comprises:
and calculating the output power when each channel is a main channel, and taking the channel corresponding to the minimum output power as a reference channel.
3. The method of claim 2, wherein the step of generating the second signal comprises generating a second signal based on the first signal and the second signalIn the step of calculating the output power when each channel is a main channel, the output power P when m channels are main channels is calculated according to the following formulam：
Figure FDA0002543073500000011
Where μ is a constant coefficient, and m is a constraint vector in the case of a main channel, RxxReceiving a data correlation matrix, s, for a channelmIs a constraint vector.
4. The method of claim 1, wherein determining magnitude-phase difference correction coefficients for each non-reference channel based on frequency domain data for the reference channel comprises:
determining an amplitude ratio alpha of the frequency domain data of the reference channel to the frequency domain data of the non-reference channel mmkAnd a phase difference of betamkThe amplitude-phase difference correction coefficient w is determined bymk：
wmk＝αmkexp(jβmk)
Where k represents a frequency point.
5. The method according to any one of claims 1 to 4, wherein after performing inverse Fourier transform on the corrected frequency domain data of the non-reference channel to obtain time domain data, the method further comprises:
and performing finite long single-bit impulse response FIR filtering on the time domain data of the non-reference channel, and taking the time domain data obtained after filtering and the time domain data of the reference channel as output data.
6. The method of claim 5, wherein the FIR filtering the time domain data of the non-reference channel comprises:
inputting the time domain data of the non-reference channel into an N-order FIR filter group, wherein the N-order FIR filter group adopts a filter coefficient vector as follows:
Figure FDA0002543073500000021
wherein the content of the first and second substances,
Figure FDA0002543073500000022
e represents a weighted average, yaTime domain data representing said reference channel, YMRepresenting the input signals of all taps of the FIR filter bank.
7. A space-time frequency cascade adaptive filtering processing device is characterized by comprising:
the first signal conversion module is used for carrying out frequency conversion amplification and analog-to-digital conversion on the input signal to obtain a plurality of paths of digital intermediate frequency signals;
the reference channel selection module is used for selecting a reference channel from the channels corresponding to the multi-channel digital intermediate frequency signals;
the second signal conversion module is used for converting the multi-channel digital intermediate frequency signals into frequency domain data through Fourier transform, and determining the amplitude-phase difference correction coefficient of each non-reference channel based on the frequency domain data of the reference channel;
the correction module is used for correcting the frequency domain data of each non-reference channel according to the amplitude-phase difference correction coefficient;
and the third signal conversion module is used for performing inverse Fourier transform on the corrected frequency domain data of the non-reference channel to obtain time domain data.
8. The apparatus of claim 7, further comprising:
and the output module is used for performing FIR filtering on the time domain data of the non-reference channel and taking the time domain data obtained after filtering and the time domain data of the reference channel as output data.
9. A space-time-frequency cascaded adaptive filtering processing device, comprising: memory, processor and computer program stored on the memory and executable on the processor, characterized in that the processor implements the method according to any of claims 1 to 6 when executing the program.
10. A computer-readable storage medium storing computer-executable instructions for performing the method of any one of claims 1-6.
CN202010552672.9A 2020-06-17 2020-06-17 Space-time frequency cascade adaptive filtering processing method, device and equipment Pending CN111934651A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010552672.9A CN111934651A (en) 2020-06-17 2020-06-17 Space-time frequency cascade adaptive filtering processing method, device and equipment

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010552672.9A CN111934651A (en) 2020-06-17 2020-06-17 Space-time frequency cascade adaptive filtering processing method, device and equipment

Publications (1)

Publication Number Publication Date
CN111934651A true CN111934651A (en) 2020-11-13

Family

ID=73317608

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010552672.9A Pending CN111934651A (en) 2020-06-17 2020-06-17 Space-time frequency cascade adaptive filtering processing method, device and equipment

Country Status (1)

Country Link
CN (1) CN111934651A (en)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113630152A (en) * 2021-10-11 2021-11-09 中国人民解放军海军工程大学 Guiding type digital anti-interception anti-interference device and method
CN116643247A (en) * 2023-07-25 2023-08-25 南京航空航天大学 Method and device for in-cabin living body detection clutter suppression

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113630152A (en) * 2021-10-11 2021-11-09 中国人民解放军海军工程大学 Guiding type digital anti-interception anti-interference device and method
CN116643247A (en) * 2023-07-25 2023-08-25 南京航空航天大学 Method and device for in-cabin living body detection clutter suppression
CN116643247B (en) * 2023-07-25 2023-10-27 南京航空航天大学 Method and device for in-cabin living body detection clutter suppression

Similar Documents

Publication Publication Date Title
US8264407B2 (en) Transmitter beamforming steering matrix processing and storage
CN108881076B (en) MIMO-FBMC/OQAM system channel estimation method based on compressed sensing
JP5484682B2 (en) WIRELESS COMMUNICATION DEVICE, EQUALIZER, EQUALIZER WEIGHT COEFFICIENT PROGRAM AND Equalizer Weight Coefficient Calculation Method
CN111934651A (en) Space-time frequency cascade adaptive filtering processing method, device and equipment
CN109104389A (en) Self-adaptive correction method for channel mismatch of GNSS antenna array
CN114189293B (en) Broadband receiving array antenna channel amplitude and phase calibration method and system
KR100606790B1 (en) channel equalizer using multi antenna
CN113055317A (en) Orthogonal matching tracking channel estimation method for underwater sound OFDM system
CN111385230B (en) Channel estimation method and system based on wiener self-adaptation
CN115250216A (en) Underwater sound OFDM combined channel estimation and signal detection method based on deep learning
CN108566347B (en) Pilot frequency design method for double-selection sparse channel of multi-user OFDM system
US20020163959A1 (en) Shortening impulse reponse fliter (SIRF) and design technique therefor
CN111934728B (en) Digital multi-beam antenna processing method, device and equipment
JP2010098734A (en) Multi-carrier signal processing method and device of the same
US8160187B1 (en) Systems and methods for canceling crosstalk from digital multi-tone (DMT) signals
KR102387182B1 (en) Apparatus and method for signal detection in a wireless communication system
CN111934650A (en) Spatial domain adaptive filtering processing method, device and equipment
CN111726308B (en) Orthogonal matching pursuit channel estimation method based on frequency response pre-interpolation
Slocumb et al. A polynomial phase parameter estimation-phase unwrapping algorithm
CN111277310B (en) Blind beam pointing airspace filtering processing method, device and equipment
CN112995078B (en) Frequency offset compensation algorithm of OFDMA uplink
EP3712626B1 (en) High-rate dft-based data manipulator and data manipulation method for high performance and robust signal processing
KR101991844B1 (en) Apparatus and method for estimating time delay
CN113824488A (en) Satellite communication non-malicious interference suppression method based on decision feedback adaptive cancellation
KR101740441B1 (en) Apparatus and method for interference cancellation in wireless communication system

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
TA01 Transfer of patent application right

Effective date of registration: 20210521

Address after: 100094 Beidou Xingtong building, 7 Fengxian East Road, Yongfeng high tech industrial base, Zhongguancun, Haidian District, Beijing

Applicant after: BEIJING BDSTAR NAVIGATION Co.,Ltd.

Address before: 100094 room A301, Beidou Xingtong building, 7 Fengxian East Road, Haidian District, Beijing

Applicant before: Hexin Xingtong Technology (Beijing) Co.,Ltd.

Applicant before: UNICORECOMM (SHANGHAI) TECHNOLOGY Co.,Ltd.

TA01 Transfer of patent application right