A time delay correction method and system based on passive intermodulation correction algorithm
By using basis functions and dichotomy in passive intermodulation correction, the problem of high complexity of delay estimation calculation method is solved, and fast and accurate delay estimation is achieved, reducing the computational complexity and implementation difficulty.
Patent Information
- Application Number
- CN202211333821.8
- Authority / Receiving Office
- CN · China
- Patent Type
- Patents(China)
- Current Assignee / Owner
- Filing Date
- 2022-10-28
- Publication Date
- 2025-01-10
- Estimated Expiration
- 2042-10-28
AI Technical Summary
In the prior art, the delay estimation calculation method is very complex and it is difficult to effectively solve the delay problem in passive intermodulation correction, especially in the decimal delay estimation process, which requires traversing multiple delay values to find the largest. The cross-correlation value of the value.
The baseband signal generator sends two send signals to the PIM signal model, obtains multiple basis functions, and uses the PIM signal model and the basis function to linearly combine it, selects any basis function as the reference signal and the received signal for cross-correlation operation, and estimates Integer delay, dichotomy is used to quickly estimate the decimal delay, reducing the complexity of the algorithm.
This method quickly approximates the maximum value of the decimal delay through dichotomy, reducing the complexity and implementation difficulty of the algorithm, and improving the efficiency and performance of delay estimation.
Smart Images

Figure CN115664623B_ABST
Abstract
Description
Technical Field
[0001] The present invention relates to the technical field of passive intermodulation correction, and in particular to a time delay correction method and system based on a passive intermodulation correction algorithm. Background Art
[0002] LTE-A introduces carrier aggregation (CA) technology to meet the requirements of throughput, data capacity, etc. in communication systems. Some passive components in RF transceivers, such as duplexers, will introduce passive intermodulation (PIM). The nonlinear characteristics of passive components cause intermodulation distortion of signals on different carriers. In addition, device connection problems or surface debris will also introduce PIM. Some PIM distortion signals will fall within the receiving frequency band, affecting the receiving performance. The PIM correction algorithm is mainly based on parameter estimation based on the correlation between the PIM in the transmitted signal and the received signal. There is a loop delay between the transmitted signal and the PIM distortion, so the transmit and receive signals must be synchronized before PIM parameter estimation, and then the correlation can be used to correct PIM distortion.
[0003] The delay can be divided into two parts: integer delay and fractional delay. The delay estimation is based on the correlation between the transmitted and received signals. When estimating fractional delay, the usual practice is to traverse the fractional delay value space in a certain step, calculate the cross-correlation function value corresponding to each delay value, find the maximum value of the correlation function, and the corresponding delay value is the estimated fractional delay. For example, the range of fractional delay is [-0.5t s , 0.5t s ), t s is the sampling period. For fractional delay, interpolation operation is required, and Lagrange interpolation is used to obtain the signal after the fractional delay. The specific interpolation formula is as follows: Among them, f(x) is the signal obtained after a fractional delay, x is the fractional delay, y i is the signal before the fractional delay, x i and x j is the sampling point value, the value range is [-N / 2, N / 2-1], N is the Lagrange interpolation order. The transmitted signal obtained from the transmitted signal is cross-correlated with the original received signal after a fractional delay to obtain the corresponding cross-correlation value. The usual solution is to set the fractional delay x from -0.5t s Traverse to 0.5t according to certain steps s, the corresponding cross-correlation value is calculated for each value, and then the maximum value is obtained by comparison. The fractional delay corresponding to the maximum cross-correlation value is the final fractional delay estimate. This solution has a large amount of computation, and it is necessary to calculate the correlation function for each value, and it is necessary to find the maximum value among all the final values, which greatly increases the complexity of the algorithm. Summary of the invention
[0004] The purpose of the present invention is to provide a delay correction method and system based on a passive intermodulation correction algorithm, so as to improve the problems of high algorithm complexity and high difficulty in implementing the algorithm in the prior art for delay estimation schemes.
[0005] The embodiment of the present invention is achieved as follows:
[0006] In a first aspect, an embodiment of the present application provides a delay correction method based on a passive intermodulation correction algorithm, which comprises the following steps:
[0007] Sending two transmission signals to the PIM signal model through a baseband signal generator, and obtaining a plurality of basis functions based on the two transmission signals and the PIM signal model;
[0008] Based on the PIM signal model and multiple basis functions, a corresponding PIM signal is obtained;
[0009] Select any basis function as a reference signal, perform a cross-correlation operation between the reference signal and the received signal containing the PIM signal, and obtain an integer delay between the received signal and the transmitted signal;
[0010] The integer delay is used to perform corresponding delay processing on the transmitted signal. After the integer delay compensation is performed on the transmitted signal, the remaining part is the fractional delay.
[0011] The fractional delay of the remaining part is estimated using the binary search method.
[0012] In some embodiments of the present invention, the step of estimating the fractional delay of the remaining part by using the binary search method includes:
[0013] Obtain and define two fractional delay values, left and right, according to the value range of the fractional delay, and calculate the difference between the two fractional delay values;
[0014] If the difference is less than the preset accuracy, left or right is directly output as a decimal delay estimate;
[0015] If the difference is greater than the preset accuracy, the corresponding cross-correlation value R is calculated respectively left and R right , and compare R left and R right The size of Rleft is the cross-correlation value corresponding to left, R right is the cross-correlation value corresponding to right;
[0016] If R left Greater than R right , then use Update the right parameter, if R left Less than R right , then use Update the left parameter and repeat this step to recursively update the left parameter or the right parameter until the difference between the two fractional delay values of left and right is less than the preset precision, and output the updated left parameter or the updated right parameter as the fractional delay estimate.
[0017] In some embodiments of the present invention, the step of selecting any basis function as a reference signal and performing a cross-correlation operation using the reference signal and a received signal including a PIM signal comprises:
[0018] The reference signal and the received signal are calculated using the cross-correlation function formula A cross-correlation operation is performed to obtain the position corresponding to the maximum value of the cross-correlation function to determine the integer delay between the received signal and the transmitted signal, where R is the cross-correlation function, N is the signal length, x is the reference signal, y is the received signal, n is the time domain sampling point of the transmitted and received signals, m is the independent variable of the cross-correlation function, and * represents the conjugate operation.
[0019] In some embodiments of the present invention, if the difference is greater than the preset accuracy, the corresponding cross-correlation value R is calculated respectively. left and R right The steps include:
[0020] Using the formula Calculate the cross-correlation value R corresponding to left left , where R is the cross-correlation function, N is the signal length, x is the reference signal, y is the received signal, n is the time domain sampling point of the transmitted and received signals, m is the independent variable of the cross-correlation function, and * represents the conjugate operation;
[0021] Using the formula Calculate the cross-correlation value R corresponding to right right , where R is the cross-correlation function, N is the signal length, x is the reference signal, y is the received signal, n is the time domain sampling point of the transmitted and received signals, m is the independent variable of the cross-correlation function, and * represents the conjugate operation.
[0022] In some embodiments of the present invention, the step of obtaining the corresponding PIM signal based on the PIM signal model and the multiple basis functions includes:
[0023] The PIM signal is obtained by linearly combining multiple basis functions.
[0024] In some embodiments of the present invention, the delay correction method based on the passive intermodulation correction algorithm further includes:
[0025] In response to the user's modification operation, the parameters of the preset precision are modified.
[0026] In a second aspect, an embodiment of the present application provides a delay correction system based on a passive intermodulation correction algorithm, which includes:
[0027] A basis function obtaining module is used to send two transmission signals to the PIM signal model through a baseband signal generator, and obtain multiple basis functions based on the two transmission signals and the PIM signal model;
[0028] A PIM signal obtaining module, used for obtaining a corresponding PIM signal based on a PIM signal model and a plurality of basis functions;
[0029] A cross-correlation operation module is used to select any basis function as a reference signal, perform a cross-correlation operation using the reference signal and a received signal containing a PIM signal, and obtain an integer delay between the received signal and the transmitted signal;
[0030] An integer delay compensation module is used to perform corresponding delay processing on the transmitted signal using the integer delay. After the transmitted signal is compensated for the integer delay, the remaining part is a fractional delay.
[0031] The fractional delay estimation module is used to estimate the fractional delay of the remaining part by using the binary method.
[0032] In some embodiments of the present invention, the above-mentioned fractional delay estimation module includes:
[0033] A difference calculation unit, used to obtain and define two fractional delay values of left and right according to the value range of the fractional delay, and calculate the difference between the two fractional delay values of left and right;
[0034] A difference comparison unit, used to directly output left or right as a decimal delay estimation value if the difference is less than a preset accuracy;
[0035] The cross-correlation value calculation unit is used to calculate the corresponding cross-correlation value R if the difference is greater than the preset accuracy. left and R right , and compare R left and R right The size of R left is the cross-correlation value corresponding to left, R right is the cross-correlation value corresponding to right;
[0036] Recursive update unit, used if R left Greater than R right , then use Update the right parameter, if R left Less than R right , then use Update the left parameter, repeatedly execute the recursive update unit to recursively update the left parameter or the right parameter until the difference between the two fractional delay values of left and right is less than the preset accuracy, and output the updated left parameter or the updated right parameter as the fractional delay estimate.
[0037] In a third aspect, an embodiment of the present application provides an electronic device, comprising a memory for storing one or more programs and a processor. When the one or more programs are executed by the processor, any method in the first aspect is implemented.
[0038] In a fourth aspect, an embodiment of the present application provides a computer-readable storage medium having a computer program stored thereon, which, when executed by a processor, implements a method as described in any one of the first aspects above.
[0039] Compared with the prior art, the embodiments of the present invention have at least the following advantages or beneficial effects:
[0040] The present invention provides a delay correction method and system based on a passive intermodulation correction algorithm, which comprises the following steps: sending two transmission signals to a PIM signal model through a baseband signal generator, and obtaining multiple basis functions based on the two transmission signals and the PIM signal model. Based on the PIM signal model and the multiple basis functions, a corresponding PIM signal is obtained. Any basis function is selected as a reference signal, and a cross-correlation operation is performed using the reference signal and a received signal containing the PIM signal to obtain an integer delay between the received signal and the transmission signal. The transmission signal is subjected to corresponding delay processing using the integer delay, and after the transmission signal is compensated for the integer delay, the remaining part is a fractional delay. The fractional delay of the remaining part is estimated using a dichotomy method. Aiming at the delay estimation scheme, the method and system firstly obtain multiple basis functions by using two transmission signals and a PIM signal model, obtain the corresponding PIM signal based on the PIM signal model and the multiple basis functions, and then select any basis function in the PIM signal model as a reference signal to perform a cross-correlation operation with a received signal containing the PIM signal for integer delay estimation, and use the binary search method to estimate the fractional delay, thereby reducing the algorithm complexity and the difficulty of algorithm implementation. BRIEF DESCRIPTION OF THE DRAWINGS
[0041] In order to more clearly illustrate the technical solutions of the embodiments of the present invention, the drawings required for use in the embodiments are briefly introduced below. It should be understood that the following drawings only show certain embodiments of the present invention and therefore should not be regarded as limiting the scope. For ordinary technicians in this field, other related drawings can be obtained based on these drawings without creative work.
[0042] Figure 1 A flowchart of a delay correction method based on a passive intermodulation correction algorithm provided by an embodiment of the present invention;
[0043] Figure 2 A signal processing flow chart provided by an embodiment of the present invention;
[0044] Figure 3 A flowchart of a fractional delay estimation method provided by an embodiment of the present invention;
[0045] Figure 4 A structural block diagram of a time delay correction system based on a passive intermodulation correction algorithm provided by an embodiment of the present invention;
[0046] Figure 5 A schematic structural block diagram of an electronic device provided by an embodiment of the present invention.
[0047] Icons: 110 - basis function acquisition module; 120 - PIM signal acquisition module; 130 - cross-correlation operation module; 140 - integer delay compensation module; 150 - fractional delay estimation module; 101 - memory; 102 - processor; 103 - communication interface. DETAILED DESCRIPTION
[0048] In order to make the purpose, technical solution and advantages of the embodiments of the present application clearer, the technical solution in the embodiments of the present application will be clearly and completely described below in conjunction with the drawings in the embodiments of the present application. Obviously, the described embodiments are part of the embodiments of the present application, rather than all the embodiments. The components of the embodiments of the present application described and shown in the drawings here can be arranged and designed in various different configurations.
[0049] Therefore, the following detailed description of the embodiments of the present application provided in the accompanying drawings is not intended to limit the scope of the present application for which protection is sought, but merely represents selected embodiments of the present application. Based on the embodiments in the present application, all other embodiments obtained by ordinary technicians in the field without creative work are within the scope of protection of the present application.
[0050] It should be noted that similar reference numerals and letters represent similar items in the following drawings, so once an item is defined in one drawing, it does not need to be further defined and explained in the subsequent drawings. At the same time, in the description of this application, if the terms "first", "second", etc. appear, they are only used to distinguish the description and cannot be understood as indicating or implying relative importance.
[0051] It should be noted that, in this article, relational terms such as first and second, etc. are only used to distinguish one entity or operation from another entity or operation, and do not necessarily require or imply any such actual relationship or order between these entities or operations. Moreover, if the terms "include", "comprise" or any other variant thereof appear, it is intended to cover non-exclusive inclusion, so that the process, method, article or equipment including a series of elements includes not only those elements, but also includes other elements not explicitly listed, or also includes elements inherent to such process, method, article or equipment. In the absence of further restrictions, if an element defined by the sentence "comprises a..." appears, it does not exclude the presence of other identical elements in the process, method, article or equipment including the element.
[0052] In the description of the present application, it should be noted that if the terms "upper", "lower", "inside", "outside", etc. appear to indicate an orientation or position relationship, they are based on the orientation or position relationship shown in the drawings, or are the orientation or position relationship in which the product of the application is usually placed when used. This is only for the convenience of describing the present application and simplifying the description, and does not indicate or imply that the device or element referred to must have a specific orientation, be constructed and operated in a specific orientation, and therefore cannot be understood as a limitation on the present application.
[0053] In the description of this application, it is also necessary to explain that, unless otherwise clearly specified and limited, the terms "set" and "connection" should be understood in a broad sense, for example, it can be a fixed connection, a detachable connection, or an integral connection; it can be a mechanical connection or an electrical connection; it can be a direct connection, or it can be an indirect connection through an intermediate medium, or it can be the internal connection of two components. For ordinary technicians in this field, the specific meanings of the above terms in this application can be understood according to specific circumstances.
[0054] In conjunction with the accompanying drawings, some embodiments of the present application are described in detail below. In the absence of conflict, the following embodiments and features in the embodiments can be combined with each other.
[0055] Example
[0056] Please refer to Figure 1 , Figure 1The figure is a flow chart of a delay correction method based on a passive intermodulation correction algorithm provided by an embodiment of the present invention. The embodiment of the present application provides a delay correction method based on a passive intermodulation correction algorithm, which comprises the following steps:
[0057] S110: Send two transmission signals to the PIM signal model through a baseband signal generator, and obtain multiple basis functions based on the two transmission signals and the PIM signal model;
[0058] Specifically, there are multiple models of PIM signals, each model includes its corresponding basis functions, and the multiple basis functions are obtained from two transmitted signals based on the PIM signal model.
[0059] S120: Obtain a corresponding PIM signal based on the PIM signal model and multiple basis functions;
[0060] Specifically, based on the PIM signal model, multiple basis functions are linearly combined to obtain a PIM signal, and the PIM signal will fall within the receiving frequency band, so that the received signal contains the PIM signal.
[0061] S130: selecting any basis function as a reference signal, performing a cross-correlation operation using the reference signal and a received signal including a PIM signal, and obtaining an integer delay between the received signal and the transmitted signal;
[0062] Specifically, since the PIM signal model contains the high-order components of the transmitted signal, the transmitted signal and the received signal cannot be directly cross-correlated. Instead, the corresponding PIM signal must be obtained based on the PIM signal model and the PIM distortion order. The PIM signal will fall within the receiving frequency band, so that the received signal contains the PIM signal. Then, one of the multiple basis functions is selected as a reference signal. After the reference signal is obtained, the reference signal is cross-correlated with the received signal to find the position corresponding to the maximum value of the cross-correlation function, and the integer delay between the transmitted and received signals can be estimated.
[0063] It should be noted that the above PIM distortion order is related to signal modeling. The corresponding PIM distortion order is the same as the order of PIM distortion modeling.
[0064] S140: Performing corresponding delay processing on the transmitted signal by using the integer delay, and after the transmitted signal is compensated for the integer delay, the remaining part is the fractional delay;
[0065] S150: Use the binary search method to estimate the fractional delay of the remaining part.
[0066] In the above implementation process, the method, for the delay estimation scheme, first uses two transmitted signals and a PIM signal model to obtain multiple basis functions, obtains the corresponding PIM signal based on the PIM signal model and the multiple basis functions, and then selects any basis function in the PIM signal model as a reference signal to perform a cross-correlation operation with a received signal containing the PIM signal for integer delay estimation, and uses a binary search method to estimate fractional delay, thereby reducing the algorithm complexity and the difficulty of algorithm implementation.
[0067] Please refer to Figure 2 , Figure 2 The figure shows a signal processing flow chart provided by an embodiment of the present invention. First, two transmission signals TX1 and TX2 are sent through a baseband signal generator. The two transmission signals will be transmitted to the PIM signal model, and a PIM signal will be generated through the PIM signal model. The PIM signal will fall within the receiving frequency band, so that the receiving signal RX contains the PIM signal. Any basis function in the PIM signal model is used to perform a cross-correlation operation with the receiving signal to perform delay estimation. That is, the delay estimation is completed using the cross-correlation of the PIM distortion signal part in the transmission signal and the receiving signal.
[0068] Please refer to Figure 3 , Figure 3 FIG. 1 is a flow chart of a fractional delay estimation method provided by an embodiment of the present invention. In some implementations of this embodiment, the step of using the binary method to estimate the fractional delay of the remaining part includes:
[0069] Obtain and define two fractional delay values, left and right, according to the value range of the fractional delay, and calculate the difference between the two fractional delay values;
[0070] If the difference is less than the preset accuracy, left or right is directly output as a decimal delay estimate;
[0071] If the difference is greater than the preset accuracy, the corresponding cross-correlation value R is calculated respectively left and R right , and compare R left and R right The size of R left is the cross-correlation value corresponding to left, R right is the cross-correlation value corresponding to right;
[0072] If R left Greater than R right , then use Update the right parameter, if R left Less than R right , then use Update the left parameter and repeat this step to recursively update the left parameter or the right parameter until the difference between the two fractional delay values of left and right is less than the preset precision, and output the updated left parameter or the updated right parameter as the fractional delay estimate.
[0073] Specifically, since the relationship between the fractional delay and the cross-correlation function value satisfies a parabolic function, the binary search method can be used to find the fractional delay corresponding to the maximum value. First, define the two fractional delay values of left and right, and calculate the difference between the two. If the difference is less than the precision, then left or right can be directly output as the final fractional delay estimate. If the difference between the two is greater than the precision, then the corresponding cross-correlation values need to be calculated separately, and then the correlation values are compared. If the correlation value R corresponding to left is left Greater than the corresponding correlation value R right , then use the average of the two to update the right parameter, and vice versa, use the average of the two to update the left parameter. Repeat the above steps and recursively update the left or right parameter until the difference between the two is less than the precision, then output one of them as the estimated value of the fractional delay. This solution can more quickly approximate the fractional delay corresponding to the maximum value, and there is no need to find the maximum value of the related function, which greatly improves the algorithm performance and reduces the algorithm complexity. If you want to improve the search accuracy, you can directly modify the precision parameter. It only takes a few additional iterations to converge to below the precision. For the traversal method, the improvement of the computational complexity by modifying the precision is much greater than the binary search solution, so using the binary search method for fractional delay estimation can greatly improve the algorithm performance.
[0074] In some implementations of this embodiment, the step of selecting any basis function as a reference signal and performing a cross-correlation operation using the reference signal and a received signal including a PIM signal includes:
[0075] The reference signal and the received signal are calculated using the cross-correlation function formula A cross-correlation operation is performed to obtain the position corresponding to the maximum value of the cross-correlation function to determine the integer delay between the received signal and the transmitted signal, where R is the cross-correlation function, N is the signal length, x is the reference signal, y is the received signal, n is the time domain sampling point of the transmitted and received signals, m is the independent variable of the cross-correlation function, and * represents the conjugate operation.
[0076] In some implementations of this embodiment, if the difference is greater than the preset accuracy, the corresponding cross-correlation value R is calculated respectively. left and R right The steps include:
[0077] Using the formula Calculate the cross-correlation value R corresponding to left left , where R is the cross-correlation function, N is the signal length, x is the reference signal, y is the received signal, n is the time domain sampling point of the transmitted and received signals, m is the independent variable of the cross-correlation function, and * represents the conjugate operation;
[0078] Using the formula Calculate the cross-correlation value R corresponding to right right , where R is the cross-correlation function, N is the signal length, x is the reference signal, y is the received signal, n is the time domain sampling point of the transmitted and received signals, m is the independent variable of the cross-correlation function, and * represents the conjugate operation.
[0079] In some implementations of this embodiment, the step of obtaining the corresponding PIM signal based on the PIM signal model and the multiple basis functions includes:
[0080] The PIM signal is obtained by linearly combining multiple basis functions.
[0081] In some implementations of this embodiment, the delay correction method based on the passive intermodulation correction algorithm further includes:
[0082] In response to the user's modification operation, the parameters of the preset precision are modified.
[0083] Please refer to Figure 4 , Figure 4 The figure shows a structural block diagram of a time delay correction system based on a passive intermodulation correction algorithm provided by an embodiment of the present invention. The embodiment of the present application provides a time delay correction system based on a passive intermodulation correction algorithm, which includes:
[0084] A basis function obtaining module 110 is used to send two transmission signals to the PIM signal model through a baseband signal generator, and obtain a plurality of basis functions based on the two transmission signals and the PIM signal model;
[0085] A PIM signal obtaining module 120 is used to obtain a corresponding PIM signal based on a PIM signal model and a plurality of basis functions;
[0086] A cross-correlation operation module 130 is used to select any basis function as a reference signal, perform a cross-correlation operation using the reference signal and a received signal including a PIM signal, and obtain an integer delay between the received signal and the transmitted signal;
[0087] An integer delay compensation module 140 is used to perform corresponding delay processing on the transmission signal by using the integer delay, and the remaining part of the transmission signal after the integer delay compensation is the fractional delay;
[0088] The fractional delay estimation module 150 is used to estimate the fractional delay of the remaining part by using the binary method.
[0089] In the above implementation process, the system first uses two transmitted signals and a PIM signal model to obtain multiple basis functions for the delay estimation scheme, obtains the corresponding PIM signal based on the PIM signal model and the multiple basis functions, and then selects any basis function in the PIM signal model as a reference signal to perform a cross-correlation operation with a received signal containing the PIM signal for integer delay estimation, and uses a binary search method to estimate fractional delays, thereby reducing the algorithm complexity and the difficulty of algorithm implementation.
[0090] In some implementations of this embodiment, the fractional delay estimation module 150 includes:
[0091] A difference calculation unit, used to obtain and define two fractional delay values of left and right according to the value range of the fractional delay, and calculate the difference between the two fractional delay values of left and right;
[0092] A difference comparison unit, used to directly output left or right as a decimal delay estimation value if the difference is less than a preset accuracy;
[0093] The cross-correlation value calculation unit is used to calculate the corresponding cross-correlation value R if the difference is greater than the preset accuracy. left and R right , and compare R left and R right The size of R left is the cross-correlation value corresponding to left, R right is the cross-correlation value corresponding to right;
[0094] Recursive update unit, used if R left Greater than R right , then use Update the right parameter, if R left Less than R right , then use Update the left parameter, repeatedly execute the recursive update unit to recursively update the left parameter or the right parameter until the difference between the two fractional delay values of left and right is less than the preset accuracy, and output the updated left parameter or the updated right parameter as the fractional delay estimate.
[0095] Specifically, since the relationship between the fractional delay and the cross-correlation function value satisfies a parabolic function, the binary search method can be used to find the fractional delay corresponding to the maximum value. First, define the two fractional delay values of left and right, and calculate the difference between the two. If the difference is less than the precision, then left or right can be directly output as the final fractional delay estimate. If the difference between the two is greater than the precision, then the corresponding cross-correlation values need to be calculated separately, and then the correlation values are compared. If the correlation value R corresponding to left is left Greater than the corresponding correlation value R right , then use the average of the two to update the right parameter, and vice versa, use the average of the two to update the left parameter. Repeat the above steps and recursively update the left or right parameter until the difference between the two is less than the precision, then output one of them as the estimated value of the fractional delay. This solution can more quickly approximate the fractional delay corresponding to the maximum value, and there is no need to find the maximum value of the related function, which greatly improves the algorithm performance and reduces the algorithm complexity. If you want to improve the search accuracy, you can directly modify the precision parameter. It only takes a few additional iterations to converge to below the precision. For the traversal method, the improvement of the computational complexity by modifying the precision is much greater than the binary search solution, so using the binary search method for fractional delay estimation can greatly improve the algorithm performance.
[0096] Please refer to Figure 5 , Figure 5 A schematic structural block diagram of an electronic device provided in an embodiment of the present application. The electronic device includes a memory 101, a processor 102 and a communication interface 103, and the memory 101, the processor 102 and the communication interface 103 are electrically connected to each other directly or indirectly to realize data transmission or interaction. For example, these elements can be electrically connected to each other through one or more communication buses or signal lines. The memory 101 can be used to store software programs and modules, such as program instructions / modules corresponding to a delay correction system based on a passive intermodulation correction algorithm provided in an embodiment of the present application. The processor 102 executes various functional applications and data processing by executing the software programs and modules stored in the memory 101. The communication interface 103 can be used to communicate signaling or data with other node devices.
[0097] Among them, the memory 101 can be, but is not limited to, a random access memory (RAM), a read only memory (ROM), a programmable read-only memory (PROM), an erasable programmable read-only memory (EPROM), an electrically erasable read-only memory (EEPROM), etc.
[0098] The processor 102 may be an integrated circuit chip with signal processing capability. The processor 102 may be a general-purpose processor, including a central processing unit (CPU), a network processor (NP), etc.; it may also be a digital signal processor (DSP), an application specific integrated circuit (ASIC), a field programmable gate array (FPGA), or other programmable logic devices, discrete gate or transistor logic devices, or discrete hardware components.
[0099] Understandably, Figure 5 The structure shown is for illustration only. The electronic device may also include Figure 5 More or fewer components as shown, or with Figure 5 Different configurations shown. Figure 5 Each component shown in the figure can be implemented by hardware, software or a combination thereof.
[0100] In the embodiments provided in the present application, it should be understood that the disclosed devices and methods can also be implemented in other ways. The device embodiments described above are merely schematic. For example, the flowcharts and block diagrams in the accompanying drawings show the possible architecture, functions and operations of the devices, methods and computer program products according to the multiple embodiments of the present application. In this regard, each box in the flowchart or block diagram can represent a module, a program segment or a part of a code, and the module, a program segment or a part of a code contains one or more executable instructions for implementing the specified logical function. It should also be noted that in some alternative implementations, the functions marked in the box can also occur in a different order from the order marked in the accompanying drawings. For example, two consecutive boxes can actually be executed substantially in parallel, and they can sometimes be executed in the opposite order, depending on the functions involved. It should also be noted that each box in the block diagram and / or the flowchart, and the combination of boxes in the block diagram and / or the flowchart can be implemented with a dedicated hardware-based system that performs a specified function or action, or can be implemented with a combination of dedicated hardware and computer instructions.
[0101] In addition, the functional modules in the various embodiments of the present application may be integrated together to form an independent part, or each module may exist separately, or two or more modules may be integrated to form an independent part.
[0102] If the functions are implemented in the form of software function modules and sold or used as independent products, they can be stored in a computer-readable storage medium. Based on this understanding, the technical solution of the present application can be essentially or partly embodied in the form of a software product that contributes to the prior art. The computer software product is stored in a storage medium, including several instructions to enable a computer device (which can be a personal computer, a server, or a network device, etc.) to perform all or part of the steps of the methods described in each embodiment of the present application. The aforementioned storage medium includes: various media that can store program codes, such as a USB flash drive, a mobile hard disk, a read-only memory (ROM), a random access memory (RAM), a magnetic disk or an optical disk.
[0103] The above description is only the preferred embodiment of the present application and is not intended to limit the present application. For those skilled in the art, the present application may have various modifications and variations. Any modification, equivalent replacement, improvement, etc. made within the spirit and principle of the present application shall be included in the protection scope of the present application.
[0104] It will be apparent to those skilled in the art that the present application is not limited to the details of the exemplary embodiments described above, and that the present application can be implemented in other specific forms without departing from the spirit or essential features of the present application. Therefore, the embodiments should be considered exemplary and non-limiting in all respects, and the scope of the present application is defined by the appended claims rather than the above description, and it is intended that all changes falling within the meaning and scope of the equivalent elements of the claims be included in the present application. Any reference numeral in a claim should not be considered as limiting the claim to which it relates.
Claims
1. A delay correction method based on a passive intermodulation correction algorithm, characterized in that: The steps include: Sending two transmission signals to the PIM signal model through a baseband signal generator, and obtaining a plurality of basis functions based on the two transmission signals and the PIM signal model; Based on the PIM signal model and the multiple basis functions, obtain a corresponding PIM signal; Selecting any basis function as a reference signal, performing a cross-correlation operation using the reference signal and a received signal including the PIM signal, and obtaining an integer delay between the received signal and the transmitted signal; Performing corresponding delay processing on the transmission signal by using the integer delay, wherein after the transmission signal is compensated for the integer delay, the remaining part is a fractional delay; Using a binary search method to estimate the fractional delay of the remaining part; The step of estimating the fractional delay of the remaining part by using the dichotomy method comprises: Obtain and define two fractional delay values, left and right, according to the value range of the fractional delay, and calculate the difference between the two fractional delay values; If the difference is less than the preset accuracy, left or right is directly output as a decimal delay estimate; If the difference is greater than the preset accuracy, the corresponding cross-correlation value R is calculated respectively. left and R right , and compare R left and R right The size of R left is the cross-correlation value corresponding to left, R right is the cross-correlation value corresponding to right; If R left Greater than R right , then use Update the right parameter, if R left Less than R right , then use Update the left parameter, repeat this step to recursively update the left parameter or the right parameter until the difference between the two fractional delay values of left and right is less than the preset accuracy, and output the updated left parameter or the updated right parameter as the fractional delay estimate.
2. The delay correction method based on the passive intermodulation correction algorithm according to claim 1, characterized in that: The step of selecting any basis function as a reference signal and performing a cross-correlation operation using the reference signal and a received signal including the PIM signal comprises: The reference signal and the received signal are calculated using the cross-correlation function formula A cross-correlation operation is performed to obtain the position corresponding to the maximum value of the cross-correlation function to determine the integer delay between the received signal and the transmitted signal, where R is the cross-correlation function, N is the signal length, x is the reference signal, y is the received signal, n is the time domain sampling point of the transmitted and received signals, m is the independent variable of the cross-correlation function, and * represents the conjugate operation.
3. The delay correction method based on the passive intermodulation correction algorithm according to claim 1, characterized in that: If the difference is greater than the preset accuracy, the corresponding cross-correlation value R is calculated respectively. left and R right The steps include: Using the formula Calculate the cross-correlation value R corresponding to left left , where R is the cross-correlation function, N is the signal length, x is the reference signal, y is the received signal, n is the time domain sampling point of the transmitted and received signals, m is the independent variable of the cross-correlation function, and * represents the conjugate operation; Using the formula Calculate the cross-correlation value R corresponding to right right , where R is the cross-correlation function, N is the signal length, x is the reference signal, y is the received signal, n is the time domain sampling point of the transmitted and received signals, m is the independent variable of the cross-correlation function, and * represents the conjugate operation.
4. The delay correction method based on the passive intermodulation correction algorithm according to claim 1, characterized in that: Based on the PIM signal model and the multiple basis functions, the step of obtaining a corresponding PIM signal includes: The PIM signal is obtained by linearly combining multiple basis functions.
5. The delay correction method based on the passive intermodulation correction algorithm according to claim 1, characterized in that: Also includes: In response to the user's modification operation, the parameters of the preset precision are modified.
6. A delay correction system based on a passive intermodulation correction algorithm, characterized in that: include: A basis function obtaining module, used for sending two transmission signals to the PIM signal model through a baseband signal generator, and obtaining a plurality of basis functions based on the two transmission signals and the PIM signal model; A PIM signal obtaining module, configured to obtain a corresponding PIM signal based on the PIM signal model and the multiple basis functions; A cross-correlation operation module, used to select any basis function as a reference signal, perform a cross-correlation operation using the reference signal and a received signal including the PIM signal, and obtain an integer delay between the received signal and the transmitted signal; An integer delay compensation module, used to perform corresponding delay processing on the transmission signal by using the integer delay, wherein the remaining part of the transmission signal after the integer delay compensation is a fractional delay; A fractional delay estimation module, used for performing fractional delay estimation on the remaining part by using a binary method; The step of estimating the fractional delay of the remaining part by using the dichotomy method comprises: Obtain and define two fractional delay values, left and right, according to the value range of the fractional delay, and calculate the difference between the two fractional delay values; If the difference is less than the preset accuracy, left or right is directly output as a decimal delay estimate; If the difference is greater than the preset accuracy, the corresponding cross-correlation value R is calculated respectively. left and R right , and compare R left and R right The size of R left is the cross-correlation value corresponding to left, R right is the cross-correlation value corresponding to right; If R left Greater than R right , then use Update the right parameter, if R left Less than R right , then use Update the left parameter, repeat this step to recursively update the left parameter or the right parameter until the difference between the two fractional delay values of left and right is less than the preset accuracy, and output the updated left parameter or the updated right parameter as the fractional delay estimate.
7. The delay correction system based on the passive intermodulation correction algorithm according to claim 6, characterized in that: The fractional delay estimation module comprises: A difference calculation unit, used to obtain and define two fractional delay values of left and right according to the value range of the fractional delay, and calculate the difference between the two fractional delay values of left and right; A difference comparison unit, configured to directly output left or right as a decimal delay estimation value if the difference is less than a preset accuracy; A cross-correlation value calculation unit is used to calculate the corresponding cross-correlation value R respectively if the difference is greater than a preset accuracy. left and R right , and compare R left and R right The size of R left is the cross-correlation value corresponding to left, R right is the cross-correlation value corresponding to right; Recursive update unit, used if R left Greater than R right , then use Update the right parameter, if R left Less than R right , then use Update the left parameter, repeatedly execute the recursive updating unit to recursively update the left parameter or the right parameter until the difference between the two fractional delay values of left and right is less than the preset accuracy, and output the updated left parameter or the updated right parameter as the fractional delay estimate.
8. An electronic device, characterized in that: include: A memory for storing one or more programs; processor; When the one or more programs are executed by the processor, the method according to any one of claims 1 to 5 is implemented.
9. 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 according to any one of claims 1 to 5 is implemented.
Citation Information
Patent Citations
Ultrasonic signal time delay estimation method and system
CN104605888A
Channel estimation method and evolved Node B
CN106161296A