Compensation method and device for energy leakage in signal processing and storage medium

Abstract

The invention discloses a compensation method and device for energy leakage in signal processing and a storage medium, and the method comprises the steps: calculating a noise reduction power time delay spectrum, the maximum diameter of the noise reduction power time delay spectrum and the related energy and phase information of an adjacent diameter, calculating the energy ratio and phase change based on the related energy and phase information, judging the channel type, and carrying out the compensation of the energy leakage in the signal processing. Judging whether the channel is a low-delay channel, and if not, directly performing delay spread estimation on the channel; and if yes, selecting a position point with serious energy leakage to carry out leakage energy ratio estimation, compensating the leakage energy of each position point to an original position, and then carrying out time delay spread estimation. The obtained power delay spectrum is more matched with a real channel environment, the problem of energy leakage caused by decimal sampling deviation and IFFT introduction is solved, the energy ratio and phase change are used for channel judgment, the probability that an effective path is misjudged as leaked energy can be reduced, the estimation accuracy of delay extension in an AWGN or low delay extension scene is further improved, and the estimation efficiency is improved. And the demodulation capability in the scene is improved.

Description

Technical Field

[0001] The present invention relates to the field of wireless communications, and more specifically, to a method, device and storage medium for compensating energy leakage in signal processing. Background Art

[0002] With the high integration of wireless communication technology, artificial intelligence, Internet of Things, satellite communication and other technologies, more possibilities are brought to the development and application of wireless communication, while also bringing huge challenges to the reception and processing of wireless signals. In the field of signal processing, energy leakage is a common problem, especially when performing spectrum analysis. This usually occurs when a data window of finite length (such as the window function in the fast Fourier transform FFT) is used to analyze a signal of infinite length or unclear periodicity. Due to the truncation of the data window, the energy originally concentrated on a certain frequency component will be dispersed to adjacent frequencies, resulting in distortion of the spectrum analysis results, which is the so-called "energy leakage". Energy leakage not only affects the accuracy of spectrum analysis, but may also have an adverse effect on subsequent signal processing steps such as filtering, demodulation, and recognition, thereby reducing the performance of the entire signal processing system. Therefore, finding an effective energy leakage compensation method is of great significance to improving the accuracy and reliability of signal processing.

[0003] According to the patent disclosure (announcement) number: CN 103188008 A, a wireless channel type identification method and device are disclosed, which calculates the time domain channel impulse response of the received signal and selects the effective path based on it; then calculates the channel mean square delay spread and the channel maximum delay spread value according to the position and power of the effective path; and identifies the corresponding channel type according to the channel mean square delay spread value and the channel maximum delay spread value; and estimates the channel type of each antenna separately, and performs weighted combination of the channel types identified by each antenna based on the signal quality to obtain the final identified channel type.

[0004] In the prior art of the above patent, the delay spread correlation value of the signal is estimated first, and the channel type is estimated based on the delay spread value. However, the above operation and the residual time offset will lead to fractional sampling deviations, and IFFT zero padding in the frequency domain is equivalent to time domain interpolation, both of which will lead to energy leakage. The estimation of small delay spread values ​​based on the power delay spectrum of energy leakage will result in large deviations, which will lead to errors in channel type determination. Summary of the invention

[0005] The present invention aims to overcome at least one defect in the prior art and provides a method, device and storage medium for compensating energy leakage in signal processing, which are used to solve the problem of identifying low delay spread channels and compensating for leaked energy and correcting the delay spread estimation value in this case.

[0006] The technical method adopted by the present invention is a method for compensating energy leakage in signal processing, comprising the following steps: Step S100: Calculate the noise reduction power delay spectrum; Step S200: Calculating the relevant energy and phase information of the maximum path and adjacent paths of the noise reduction power delay spectrum; Step S300: Calculate energy proportion and phase change based on the correlation energy and the phase information to determine the channel type; Step S400: Selecting a location point with serious energy leakage to estimate the leakage energy ratio; Step S500: compensating the leakage energy of each position point back to the original position; Step S600: performing delay spread estimation; Before step S400, it is necessary to determine whether the channel is a low-delay channel. If so, continue to execute the remaining steps. If not, directly execute step S600. By estimating the channel type first, and then performing power delay profile (power delay profile) energy correction to obtain a more accurate delay spread value.

[0007] Step S100 specifically includes: estimating the estimated matrix according to the pilot position on the configured frequency domain resource block Perform fast Fourier transform to obtain the time domain impulse response , and then calculate its power delay spectrum :

[0008] represents the conjugate transposed time domain impulse response, for IFFT The number of points, Based on LS The frequency domain impulse response obtained by the channel estimation algorithm; Get the power delay spectrum Then, the noise reduction power delay spectrum is calculated :

[0009] is the noise power.

[0010] Step S200 specifically includes searching for the energy of the maximum path in the noise reduction power delay spectrum. and its index , and the energy of the two adjacent sample points on the left and right of the maximum diameter ( , ) and index ( , ), and obtain the time domain impulse response values ​​corresponding to the above three points ( , , ).

[0011] In step S300, the energy proportion and phase change are obtained by calculating the phase change value according to the phase information of the maximum diameter and the two sample points adjacent to the maximum diameter, and calculating the energy proportion according to the energy information of the maximum diameter and the two sample points adjacent to the maximum diameter. First, calculate the phase change:

[0012]

[0013] is the phase value, is the phase change value, is the time domain impulse response matrix, is the four-quadrant inverse tangent function; Secondly, calculate the energy ratio :

[0014] In step S300, the channel type is determined, specifically, based on the phase change value and energy ratio Determine the channel type: There are two bases for judging AWGN channels: The first judgment basis is: the energy of the maximum path and its adjacent sample points accounts for more than 98% of the total signal path energy, that is, ; The second type of judgment requires two points to be met at the same time: ① The maximum path and its adjacent paths on both sides account for more than 90% and less than 98% of the total signal path energy, that is, ; ② Affected by the fractional sampling deviation, the sample points on both sides of the maximum diameter are regarded as the energy leakage of the maximum diameter, and their phases will not change. The phase difference between the first two points and the last two points among the three points should be less than 0.03, that is, , when both ① and ② are satisfied, it is determined to be an AWGN channel; In addition to AWGN channels, for low delay spread channels where the sampling rate is too low to distinguish multiple adjacent signal paths, or the signal path energy is concentrated near the maximum path, it is also necessary to correct the accuracy of delay spread estimation through maximum path leakage energy compensation. The judgment condition for such channels is that the maximum path and the adjacent paths on both sides occupy more than 90% of the total signal path energy and less than 98%, that is, And the phase difference between the first two points and the last two points among the three points is greater than 0.03, that is, ; After the two channels are successfully determined, leakage energy compensation is required, otherwise no compensation is performed. 98%, 90%, and 0.03 are recommended values, and the specific data can be adaptively modified according to the specific situation.

[0015] S400 specifically includes: first, performing linear interpolation on the energy of the maximum diameter and its adjacent sample points, and calculating the fractional time deviation :

[0016] Secondly, calculate the leakage energy power of the maximum diameter :

[0017] represents the imaginary part, and the energy leakage of the maximum diameter obeys function; Therefore, the leakage energy for

[0018] in, is the number of pilot subcarriers, for The number of points, , is the distance between the sample point and the maximum diameter. The leakage energy is normalized to obtain the ratio of leakage energy to the maximum diameter energy:

[0019] when is the maximum diameter position.

[0020] Step S500 specifically includes, first, selecting some locations with larger leakage energy and calculating their leakage energy size:

[0021] Secondly, the energy of the maximum diameter is corrected and location :

[0022]

[0023]

[0024] In step S100, the time domain noise power is calculated There are many methods, including but not limited to the following three: ①In A noise window with a fixed position is selected to calculate the noise power; ②In Calculate the power in segments on the spectrum, and select the power mean of 1 to 2 intervals with smaller power values ​​as the noise power; ③ The noise power calculated by the frequency domain method is converted to the time domain noise power; Select any of the above to calculate the time domain noise power .

[0025] The present invention also provides a device for compensating energy leakage in signal processing, which is constructed based on the above-mentioned method for compensating energy leakage in signal processing and includes the following modules: A noise reduction power delay spectrum calculation module, used to calculate the noise reduction power delay spectrum; Energy and phase calculation module, used to calculate the relevant energy and phase information of the maximum diameter and adjacent diameters; Channel determination module, determines the channel type based on energy proportion and phase change; A low-latency channel determination module determines whether the channel type is a low-latency channel, and if so, continues to execute; if not, jumps to a delay spread estimation module; The leakage energy estimation module is used to select the location points with serious energy leakage to estimate the leakage energy ratio; A compensation module is used to compensate the leakage energy at each position point back to the original position; Delay spread estimation module: estimates the delay spread of the channel.

[0026] The present invention also provides a computer-readable storage medium having a computer program stored thereon, and when the computer program is executed by a processor, the computer program implements the above-mentioned method for compensating energy leakage in signal processing.

[0027] Compared with the prior art, the present invention has the following beneficial effects: 1. Using energy proportion and phase change to judge the channel scene can reduce the probability of the effective path being misjudged as leakage energy, while ensuring that leakage energy compensation is only performed for scenes where the maximum path energy proportion is particularly large (when the proportion is small, the leakage energy can be ignored); 2. The power delay spectrum better matches the actual channel environment, solving the problem of fractional sampling deviation and IFFT The energy leakage problem introduced; 3. Improve the estimation accuracy of delay spread in AWGN or low delay spread scenarios, thereby improving the demodulation capability in this scenario. BRIEF DESCRIPTION OF THE DRAWINGS

[0028] The drawings of the present invention are only used for illustrative purposes and cannot be understood as limiting the present invention. In order to more clearly illustrate the technical methods of the embodiments of the present invention, the drawings required for describing the embodiments are briefly introduced below. Obviously, the drawings described below are only some embodiments of the present invention. For ordinary technicians in this field, other drawings can be obtained based on these drawings without creative work.

[0029] Figure 1 It is a flow chart of the present invention.

[0030] Figure 2 It is a system module diagram of the present invention. DETAILED DESCRIPTION

[0031] The following will be combined with the drawings in the embodiments of the present invention to clearly and completely describe the technical methods in the embodiments of the present invention. Obviously, the described embodiments are only part of the embodiments of the present invention, not all of the embodiments. Based on the embodiments of the present invention, all other embodiments obtained by ordinary technicians in this field without creative work are within the scope of protection of the present invention.

[0032] Example 1 like Figure 1 As shown, the present embodiment provides a method for compensating energy leakage in signal processing, including: step S100, calculation of noise reduction power delay spectrum; S200, calculation of maximum path and adjacent path related energy and phase and other information; S300, determining the channel type based on energy proportion and phase change; after determining the channel type, determining whether the channel is a low delay extension channel, if not, executing step S600; if so, continuing to execute the following steps; S400, selecting a location point with severe energy leakage to estimate the leakage energy proportion; S500, compensating the leakage energy of each location point back to the original position; S600, delay extension estimation. By first estimating the channel type and then performing power delay spectrum Energy correction, and based on the corrected The delay spread is estimated based on the spectrum to obtain a more accurate delay spread value.

[0033] In this embodiment, the noise reduction power delay spectrum calculation in step S100 is specifically as follows: the estimation matrix estimated according to the pilot position on the configured frequency domain resource block Perform fast Fourier transform to obtain the time domain impulse response , and then calculate its power delay spectrum :

[0034] represents the conjugate transposed time domain impulse response, for IFFT The number of points, is the frequency domain impulse response obtained based on the LS channel estimation algorithm; Get the power delay spectrum Then, calculate the time domain noise power , calculate the time domain noise power There are many ways to do this, for example, 1) Select a noise window with a fixed position to calculate the noise power; 2) Calculate the power in segments on the spectrum, and select the power average of 1 to 2 intervals with smaller power values ​​as the noise power; 3) Convert the noise power calculated by the frequency domain method to obtain the time domain noise power; any one of the methods can be used to calculate the time domain noise power .

[0035] Finally, the noise reduction power delay spectrum is calculated : .

[0036] In this embodiment, in step S200, the calculation of the maximum path and the adjacent path related energy and phase information is specifically performed by searching the energy of the maximum path in the noise reduction power delay spectrum. and its index , and the energy of the two adjacent sample points on the left and right of the maximum diameter ( , ) and index ( , ), and obtain the time domain impulse response values ​​corresponding to the above three points ( , , ).

[0037] In this embodiment, step S300, determining the channel type based on energy proportion and phase change is specifically to determine the channel type based on the maximum diameter and the phase information of the two adjacent sample points on the left and right of the maximum diameter. Calculate the phase change value, and calculate the energy ratio based on the energy information of the maximum diameter and the two adjacent sample points on the left and right of the maximum diameter. First, calculate the phase change value :

[0038]

[0039] is the four-quadrant inverse tangent function; Secondly, calculate the energy ratio :

[0040] According to the phase change value and energy ratio Determine the channel type: There are two bases for judging AWGN channels: The first judgment basis is: the energy of the maximum path and its adjacent sample points accounts for more than 98% of the total signal path energy, that is, ; The second type of judgment requires two points to be met at the same time: ① The maximum path and its adjacent paths on both sides account for more than 90% and less than 98% of the total signal path energy, that is, ; ② Affected by the fractional sampling deviation, the sample points on both sides of the maximum diameter are regarded as the energy leakage of the maximum diameter, and their phases will not change. The phase difference between the first two points and the last two points among the three points should be less than 0.03, that is, , when both ① and ② are satisfied, it is determined to be an AWGN channel; In addition to AWGN channels, for low delay spread channels where the sampling rate is too low to distinguish multiple adjacent signal paths, or the signal path energy is concentrated near the maximum path, it is also necessary to correct the accuracy of delay spread estimation through maximum path leakage energy compensation. The judgment conditions for such channels are: The maximum path and its adjacent paths on both sides account for more than 90% and less than 98% of the total signal path energy, that is, And the phase difference between the first two points and the last two points among the three points is greater than 0.03, that is, ; After the above two channels are successfully determined, leakage energy compensation is required. Otherwise, no compensation is performed and low delay spread estimation is performed directly. 98%, 90%, and 0.03 are recommended values, and the specific values ​​can be adaptively modified according to specific circumstances.

[0041] In this embodiment, step S400, selecting a location point with serious energy leakage to estimate the leakage energy ratio is specifically as follows: First, the energy of the maximum diameter and its adjacent sample points is linearly interpolated to calculate the fractional time deviation :

[0042] Secondly, the leakage energy of the maximum path is calculated, assuming that the number of pilot subcarriers is , The points are , then in the case of decimals The maximum diameter distance is The leakage power from the maximum diameter can be calculated as follows:

[0043]

[0044]

[0045]

[0046]

[0047]

[0048] represents the imaginary part, and the energy leakage of the maximum diameter obeys function, Therefore, the leakage energy is:

[0049] The leakage energy is normalized to obtain the ratio of leakage energy to maximum diameter energy:

[0050] when is the maximum diameter position.

[0051] In this embodiment, step S500, compensating the leakage energy of each position point back to the original position is specifically as follows: First, select some locations with large leakage energy and calculate their leakage energy:

[0052] Secondly, the energy of the maximum diameter is corrected and location :

[0053]

[0054]

[0055] In addition, for the scenario where the channel type is determined to be an AWGN channel, since the maximum diameter is the first diameter, the fractional multiple of the time offset obtained in step S300 is The accuracy of subsequent first path estimation can be further improved.

[0056] Example 2 like Figure 2 As shown, a compensation device for energy leakage in signal processing provided in this embodiment is constructed based on a compensation method for energy leakage in signal processing in embodiment 1, and includes the following modules: A noise reduction power delay spectrum calculation module, used to calculate the noise reduction power delay spectrum; Energy and phase calculation module, used to calculate the relevant energy and phase information of the maximum diameter and adjacent diameters; Channel determination module, determines the channel type based on energy proportion and phase change; A low-latency channel determination module determines whether the channel type is a low-latency channel, and if so, continues to execute; if not, jumps to a delay spread estimation module; The leakage energy estimation module is used to select the location points with serious energy leakage to estimate the leakage energy ratio; A compensation module is used to compensate the leakage energy at each position point back to the original position; Delay spread estimation module: estimates the delay spread of the channel.

[0057] Example 3 Another embodiment provided by the present invention is a computer-readable storage medium, on which a computer program is stored. When the computer program is executed by a processor, the processor executes the steps of a method for compensating for energy leakage in signal processing in Example 1.

[0058] The preferred embodiments of the present invention disclosed above are only used to help explain the present invention. The preferred embodiments do not describe all the details in detail, nor do they limit the invention to specific implementation methods. Obviously, many modifications and changes can be made according to the content of this specification. This specification selects and specifically describes these embodiments in order to better explain the principles and practical applications of the present invention, so that those skilled in the art can understand and use the present invention well. The present invention is limited only by the claims and their full scope and equivalents.

Claims

1. A method for compensating energy leakage in signal processing, characterized in that: The following steps are involved: Step S100: Calculate the noise reduction power delay spectrum; Step S200: Calculating the relevant energy and phase information of the maximum path and adjacent paths of the noise reduction power delay spectrum; Step S300: Calculate energy proportion and phase change based on the correlation energy and the phase information to determine the channel type; Step S400: Selecting a location point with serious energy leakage to estimate the leakage energy ratio; Step S5 00: Compensate the leakage energy of each position point back to the original position; Step S600: performing delay spread estimation; Before step S400, it is necessary to determine whether the channel is a low-latency channel. If so, continue to execute the remaining steps; if not, directly execute step S600.

2. A method for compensating energy leakage in signal processing according to claim 1, characterized in that: The step S100 is specifically to estimate the estimated matrix according to the pilot position on the configured frequency domain resource block Perform fast Fourier transform to obtain the time domain impulse response , and then calculate its power delay spectrum : represents the conjugate transposed time domain impulse response, for IFFT The number of points, is the frequency domain impulse response matrix obtained based on the LS channel estimation algorithm; Get the power delay spectrum Then, the noise reduction power delay spectrum is calculated : is the noise power.

3. The method for compensating energy leakage in signal processing according to claim 2, characterized in that: The step S200 specifically includes searching for the energy of the maximum path in the noise reduction power delay spectrum. and its index , and the energy of the two adjacent sample points on the left and right of the maximum diameter ( , ) and index ( , ), and obtain the time domain impulse response values ​​corresponding to the above three points ( , , ).

4. The method for compensating energy leakage in signal processing according to claim 3, characterized in that: In step S300, the energy proportion and phase change are obtained by calculating the phase change value according to the phase information of the maximum diameter and the two sample points adjacent to the maximum diameter, and calculating the energy proportion according to the energy information of the maximum diameter and the two sample points adjacent to the maximum diameter. First, calculate the phase change: is the phase value, is the phase change value, is the time domain impulse response, is the four-quadrant inverse tangent function; Secondly, calculate the energy ratio : 。 5. The method for compensating energy leakage in signal processing according to claim 4, characterized in that: In step S300, the channel type is determined by, specifically, determining the channel type according to the phase change value and the energy proportion value Determine the channel type: There are two bases for judging AWGN channels: The first judgment basis is: the energy of the maximum path and its adjacent sample points accounts for more than 98% of the total signal path energy, that is, , determined to be an AWGN channel; The second type of judgment requires two points to be met at the same time: ① The maximum path and its adjacent paths on both sides account for more than 90% and less than 98% of the total signal path energy, that is, ; ② Affected by the fractional sampling deviation, the sample points on both sides of the maximum diameter are regarded as the energy leakage of the maximum diameter, and their phases will not change. Therefore, the absolute value of the phase difference between the first two points and the last two points among the three points should be less than 0.03, that is, , when both ① and ② are satisfied, it is determined to be an AWGN channel; In addition to AWGN channels, for low delay spread channels where the sampling rate is too low to distinguish multiple adjacent signal paths, or the signal path energy is concentrated near the maximum path, it is also necessary to correct the accuracy of delay spread estimation through maximum path leakage energy compensation. The judgment condition for such channels is that the maximum path and the adjacent paths on both sides occupy more than 90% of the total signal path energy and less than 98%, that is, And the absolute value of the phase difference between the first two points and the last two points among the three points is greater than 0.03, that is, ; After the two channels are determined successfully, leakage energy compensation is required, otherwise no compensation is performed.

6. A method for compensating energy leakage in signal processing according to claim 5, characterized in that: Step S400 specifically includes: First, the energy of the maximum diameter and its adjacent sample points is linearly interpolated to calculate the fractional time deviation : Secondly, calculate the leakage energy power of the maximum diameter : represents the imaginary part, is the function obeyed by the energy leakage of the maximum diameter; Therefore, the leakage energy for in, is the number of pilot subcarriers, for The number of points, , is the distance between the sample point and the maximum diameter. The leakage energy is normalized to obtain the ratio of leakage energy to the maximum diameter energy: when is the maximum diameter position.

7. A method for compensating energy leakage in signal processing according to claim 6, characterized in that: The step S500 specifically includes, firstly, selecting some locations with larger leakage energy and calculating their leakage energy. size: Secondly, the energy of the maximum diameter is corrected and location : 。 8. The method for compensating energy leakage in signal processing according to claim 7, characterized in that: In step S100, the time domain noise power is calculated There are three ways: ①In A noise window with a fixed position is selected to calculate the noise power; ②In Calculate the power in segments on the spectrum, and select the power mean of 1 to 2 intervals with smaller power values ​​as the noise power; ③ The noise power calculated by the frequency domain method is converted to the time domain noise power; Select any of the above to calculate the time domain noise power .

9. A device for compensating energy leakage in signal processing, characterized in that: The method for compensating energy leakage in signal processing according to claims 1 to 8 is constructed, and includes the following modules: A noise reduction power delay spectrum calculation module, used to calculate the noise reduction power delay spectrum; Energy and phase calculation module, used to calculate the relevant energy and phase information of the maximum diameter and adjacent diameters; Channel determination module, determines the channel type based on energy proportion and phase change; A low-latency channel determination module determines whether the channel type is a low-latency channel, and if so, continues to execute; if not, jumps to a delay spread estimation module; The leakage energy estimation module is used to select the location points with serious energy leakage to estimate the leakage energy ratio; A compensation module is used to compensate the leakage energy at each position point back to the original position; Delay spread estimation module: estimates the delay spread of the channel.

10. A computer-readable storage medium having a computer program stored thereon, characterized in that: When the computer program is executed by a processor, the method for compensating energy leakage in signal processing as claimed in any one of claims 1 to 8 is implemented.

Images