Broadband non-cooperative interference sampling array manifold design method

Abstract

The invention belongs to the technical field of broadband non-cooperative interference resistance, and particularly discloses a broadband non-cooperative interference sampling array manifold design method, which comprises the following steps: establishing a broadband non-cooperative interference sampling array model, and determining the number of array elements; according to the broadband non-cooperative interference cancellation sampling array model, calculating a transient solution in a cancellation weight convergence process and a steady-state solution after cancellation weight convergence under broadband interference; determining a performance parameter according to the transient solution and the steady-state solution, and determining a corresponding relation between the performance parameter and an array manifold; determining an initial sampling array manifold by combining the corresponding relationship between the performance parameters and the array manifold with a performance index given by a cancellation system; and performing test verification on the initial sampling array manifold, and comparing test data with a performance index to iteratively adjust the array manifold until the performance index is met. According to the invention, the coping capability of the sampling array to broadband non-cooperative interference and the performance of a cancellation system can be improved.

Description

Technical Field

[0001] The present application belongs to the technical field of anti-wideband non-cooperative interference, and more specifically, to a method for designing a wideband non-cooperative interference sampling array manifold. Background Art

[0002] At present, the electromagnetic environment is becoming increasingly complex. For many communication systems, anti-interference has become one of the issues that must be considered. Interference signals may come from cooperative systems on the same platform, or they may be high-power interference imposed intentionally or unintentionally by non-cooperative systems. And with the continuous evolution of interference technology, interference signals tend to develop in a broadband manner. Broadband interference has a wider bandwidth than narrowband interference and can interfere with more communication frequency bands, making anti-interference technology have to consider broadband interference. Interference cancellation technology is an important anti-interference technology. Sampling interference signals is an important part of interference cancellation technology. Interference cancellation requires synthesizing a cancellation signal through the sampled signal, which cancels out the interference signal with equal amplitude and anti-phase.

[0003] For cooperative interference, the interference signal is known and can be sampled by coupling, but for non-cooperative interference, the interference signal is unknown, so it is necessary to sample the interference signal through a sampling array. The manifold of the sampling array directly affects the amplitude and phase of the sampled signal, and thus affects the effect of interference cancellation. The manifold of the sampling array needs to be designed to meet the cancellation performance of the interference cancellation device.

[0004] In the process of interference cancellation, the most important role of the cancellation device is the update of the weight vector. According to the convergence state of the weight vector, it can be divided into the transient state during the convergence process of the weight vector and the steady state after the weight vector converges. At different stages, the weight vector exhibits different cancellation performance. In the transient process during the convergence of the weight vector, the convergence speed and convergence stability are crucial; in the steady state process after the convergence of the weight vector, the system's interference cancellation ratio, the array's resolution and the cancellation blind area become the main manifestations of the cancellation performance. These performances are closely related to the sampling array manifold, and the quantitative relationship between the cancellation performance and the sampling array manifold needs to be calculated through modeling and solving. This quantitative relationship and system indicators are used to guide the design of the sampling array manifold.

[0005] In order to achieve the design of sampling array manifold, the objective function is often set and then the constraints of the array manifold are solved. There are mainly the following methods: analytical and semi-analytical methods, convex optimization methods, and global optimization algorithms. However, the current methods still have the defects of difficulty in solving complex optimization problems under multiple constraints or inaccurate solution results, which leads to poor accuracy and precision of the array manifold and cannot meet the requirements of broadband non-cooperative interference sampling. Summary of the invention

[0006] In view of the defects of the prior art, the purpose of this application is to provide a broadband non-cooperative interference sampling array manifold design method, which aims to solve the problem that complex optimization problems under multiple constraints in current methods are difficult to solve or the solution results are inaccurate, which in turn leads to poor accuracy and precision of the array manifold.

[0007] To achieve the above objectives, the present application provides a broadband non-cooperative interference sampling array manifold design method, comprising:

[0008] Establish a broadband non-cooperative interference sampling array model and determine the number of array elements;

[0009] Calculating the transient solution of the cancellation weight convergence process under broadband interference and the steady-state solution after the cancellation weight convergence according to the broadband non-cooperative interference cancellation sampling array model;

[0010] Determine performance parameters according to the transient solution and the steady-state solution, and determine the corresponding relationship between the performance parameters and the array manifold; the performance parameters include array resolution, cancellation blind area, interference cancellation ratio, system convergence time and system convergence stability;

[0011] The corresponding relationship between the performance parameter and the array manifold is combined with the performance index given by the cancellation system to determine the initial sampling array manifold;

[0012] The initial sampling array manifold is tested and verified, and the test data is compared with the performance index to iteratively adjust the array manifold until the performance index is met.

[0013] This application can effectively handle broadband non-cooperative interference by establishing a suitable array model. The transient solution describes the change of the weights of the cancellation system during the convergence process, while the steady-state solution reflects the stable performance of the system after long-term operation. By obtaining the correspondence between performance parameters and array manifolds, the array design and algorithm selection can be optimized, and the initial array manifold can be iteratively adjusted through test verification and performance index comparison. This application solves the complex optimization problem of array manifold solution under multiple constraints through array design optimization, cancellation weight optimization, and performance parameter optimization to improve the accuracy and precision of the array manifold, thereby improving the sampling array's ability to respond to broadband non-cooperative interference and the performance of the cancellation system.

[0014] In some embodiments, determining the number of array elements includes:

[0015] Obtaining a relationship between the number of interference sources and the number of array elements in a sampling array, and determining the number of array elements according to the relationship; the number of interference sources is the number of unrelated linear constraints imposed on the weight vector of the array;

[0016] The number of interference sources is determined according to the maximum number of broadband non-cooperative interference sources that need to be countered in the use scenario of the cancellation device; and the number of array elements is greater than the number of linear constraints.

[0017] In some embodiments, the process of obtaining the transient solution and the steady-state solution includes:

[0018] Determine an initial weight vector, obtain instantaneous data of the initial weight vector changing with sampling time during an iterative convergence process, and use the instantaneous data as a transient solution;

[0019] Determine a converged weight vector corresponding to the initial weight vector after convergence, and use the converged weight vector as the steady-state solution;

[0020] Obtaining an equivalent amplitude vector and an equivalent phase vector of the convergence weight vector;

[0021] The relationship between the convergence weight vector and the sampling array manifold parameter is determined according to the equivalent amplitude vector and the equivalent phase vector.

[0022] This application can effectively handle broadband non-cooperative interference, improve the system's signal processing capability and interference suppression capability, and optimize the cancellation weight convergence process through transient and steady-state solution analysis, thereby improving the system's response speed and stability.

[0023] In some embodiments, determining the initial sampling array manifold includes:

[0024] The initial array aperture, the farthest distance of the boundary array elements and the array element position range are obtained, and the initial sampling array manifold is determined according to the array aperture, the farthest distance of the boundary array elements and the array element position range.

[0025] In some embodiments, the method for determining the array aperture includes:

[0026] Utilizing the steady-state solution to obtain a first relationship between array resolution and interference sampling array manifold;

[0027] Determine all interference incident angles, input the interference incident angles into the first relationship, and obtain candidate array apertures corresponding to the interference incident angles;

[0028] Determine whether the maximum array aperture corresponding to each interference incident angle meets the resolution index, and screen the candidate array apertures according to the determination result until the final array aperture is obtained;

[0029] The resolution index is the array resolution required by the cancellation device at all interference incident angles.

[0030] In some embodiments, the method for determining the farthest distance of the boundary element includes:

[0031] Using the steady-state solution to obtain a second relationship between a broadband interference cancellation ratio and an interference sampling array manifold;

[0032] Determine the farthest distance of the boundary array element according to the second relationship and the cancellation ratio index;

[0033] The wideband interference cancellation ratio is related to the signal bandwidth, and the wideband interference cancellation ratio is negatively correlated to the array element spacing.

[0034] In some embodiments, the method for determining the array element position range includes:

[0035] Using the steady-state solution, a third relationship between the cancellation blind area and the interference sampling array manifold is obtained;

[0036] According to the third relationship and the cancellation blind area index, the array element position range is determined when the cancellation blind area accounts for the smallest proportion.

[0037] In some embodiments, the method for reducing the range of array element positions includes:

[0038] Using the transient solution, obtaining a fourth relationship between the interference cancellation convergence time and the interference sampling array manifold, and a fifth relationship between the convergence stability and the interference sampling array manifold;

[0039] According to the fourth relationship and the convergence time index, and the fifth relationship and the convergence stability index, selecting an array element position combination with the shortest convergence time and the best convergence stability from the array element position range, and combining the array element positions to form a new array element position range to achieve range reduction;

[0040] The interference cancellation convergence time is the time required for the initial weight vector to converge to a steady state; and the convergence stability is the convergence stability of the weight vector taking into account the convergence step size and loop gain.

[0041] In some embodiments, the iteratively adjusting the array manifold until the performance indicator is met includes:

[0042] The initial sampling array manifold is installed on a cancellation device, and the measurement data of the initial sampling array manifold is obtained by using the cancellation device; the measurement data includes the array resolution and cancellation blind area actually measured, and the cancellation ratio, convergence time and convergence stability of the cancellation system;

[0043] Determine whether the measurement data can meet the index requirements. If all measurement data meet the performance index, determine the final sampling array manifold; if any measurement data does not meet the performance index, readjust the array manifold until all performance indicators are met.

[0044] This application reduces the impact of blind spot elimination and improves the accuracy of the system in target detection and positioning through appropriate array manifold design and optimization. Optimize array manifold and algorithm selection to accelerate the system's ability to adapt to interference in complex environments and improve real-time performance and operability.

[0045] In a second aspect, the present application provides an electronic device comprising: at least one memory for storing programs; and at least one processor for executing the programs stored in the memory. When the program stored in the memory is executed, the processor is used to execute the method described in the first aspect or any possible implementation of the first aspect.

[0046] In a third aspect, the present application provides a computer-readable storage medium, which stores a computer program. When the computer program runs on a processor, the processor executes the method described in the first aspect or any possible implementation of the first aspect.

[0047] In a fourth aspect, the present application provides a computer program product. When the computer program product runs on a processor, the processor executes the method described in the first aspect or any possible implementation of the first aspect.

[0048] It can be understood that the beneficial effects of the second to fifth aspects mentioned above can be found in the relevant description of the first aspect mentioned above, and will not be repeated here.

[0049] In general, the above technical solutions conceived by this application have the following beneficial effects compared with the prior art:

[0050] (1) This application can effectively handle broadband non-cooperative interference by establishing a suitable array model. The transient solution describes the change of the weights of the cancellation system during the convergence process, while the steady-state solution reflects the stable performance of the system after long-term operation. By obtaining the corresponding relationship between performance parameters and array manifolds, the array design and algorithm selection can be optimized, and the initial array manifold can be iteratively adjusted through test verification and performance index comparison. This application solves the solution problems existing in the complex optimization of array manifolds under multiple constraints of array design optimization, cancellation weight optimization, and performance parameter optimization to improve the accuracy and precision of the array manifold, thereby improving the sampling array's ability to respond to broadband non-cooperative interference and the performance of the cancellation system.

[0051] (2) This application is modeled from the perspective of array manifold, and the weight vector is solved according to the two states of the system. Then, a set of quantitative mapping relationships between cancellation performance and sampling array manifold is established based on the weight vector. It is applicable to different broadband non-cooperative interference cancellation systems, and different broadband non-cooperative interference sampling arrays can be designed according to different index requirements. In addition, this application is aimed at the design of broadband non-cooperative interference sampling array manifold. The interference signal can be a broadband signal. When the interference signal is a narrowband signal, the model is also applicable by simplifying the signal form. Therefore, this application improves the applicability of the non-cooperative interference cancellation system.

[0052] (3) This application provides a theoretical basis for the design of sampling arrays for non-cooperative interference cancellation, solves the difficulties of non-cooperative interference sampling, and takes into account multiple performance indicators such as interference cancellation ratio, array resolution and cancellation blind area, convergence time and convergence stability. The designed array can better improve the overall performance of the interference cancellation system. BRIEF DESCRIPTION OF THE DRAWINGS

[0053] Figure 1 This is one of the flow charts of the broadband non-cooperative interference sampling array manifold design method provided in the embodiment of the present application;

[0054] Figure 2 This is the second flow chart of the broadband non-cooperative interference sampling array manifold design method provided in the embodiment of the present application;

[0055] Figure 3 It is a schematic diagram of a broadband non-cooperative interference sampling array and cancellation principle;

[0056] Figure 4 is a graph showing the relationship between array resolution and array aperture under different interference incident angles in this embodiment;

[0057] Figure 5 is a graph showing the relationship between array resolution and interference incident angle under different numbers of array elements in this embodiment;

[0058] Figure 6 is a relationship diagram of the number of selected array elements and the interference incident angle for verifying the array resolution in this embodiment;

[0059] Figure 7 is a simulation diagram of the relationship between the interference cancellation ratio and the distance between the array element and the main antenna in this embodiment;

[0060] Figure 8 Schematic diagram of the cancellation blind zone of this embodiment, wherein (a) is a three-dimensional diagram of the array output signal-to-interference-noise ratio, and (b) is a diagram of the array 3dB cancellation blind zone;

[0061] Fig. 9Schematic diagrams of weight vector convergence speeds at different array element positions of this embodiment, wherein (a) is a weight vector phase convergence time diagram of a “zero-sum” array, (b) is a weight vector amplitude convergence time diagram of a “zero-sum” array, (c) is a weight vector phase convergence time diagram of a non-“zero-sum” array, and (d) is a weight vector amplitude convergence time diagram of a non-“zero-sum” array;

[0062] Fig.10 It is a schematic diagram of the structure of an electronic device provided in an embodiment of the present application. DETAILED DESCRIPTION

[0063] In order to make the purpose, technical solution and advantages of the present application more clearly understood, the present application is further described in detail below in conjunction with the accompanying drawings and embodiments. It should be understood that the specific embodiments described herein are only used to explain the present application and are not used to limit the present application.

[0064] The term "and / or" in this article is a description of the association relationship of associated objects, indicating that there can be three relationships. For example, A and / or B can represent: A exists alone, A and B exist at the same time, and B exists alone. The symbol " / " in this article indicates that the associated objects are in an or relationship, for example, A / B means A or B.

[0065] The terms "first" and "second" in the specification and claims herein are used to distinguish different objects rather than to describe a specific order of objects. For example, a first relationship and a second relationship are used to distinguish different relationships rather than to describe a specific order of relationships.

[0066] In the embodiments of the present application, words such as "exemplary" or "for example" are used to indicate examples, illustrations or descriptions. Any embodiment or design described as "exemplary" or "for example" in the embodiments of the present application should not be interpreted as being more preferred or more advantageous than other embodiments or designs. Specifically, the use of words such as "exemplary" or "for example" is intended to present related concepts in a specific way.

[0067] In the description of the embodiments of the present application, unless otherwise specified, "multiple" means two or more than two. For example, multiple processing units refer to two or more processing units, etc.; multiple elements refer to two or more elements, etc.

[0068] Next, the technical solutions provided in the embodiments of the present application are introduced.

[0069] Reference Figure 1 The present application provides a broadband non-cooperative interference sampling array manifold design method, comprising:

[0070] S1. Establish a broadband non-cooperative interference sampling array model and determine the number of array elements;

[0071] S2. Calculating the transient solution of the cancellation weight convergence process under broadband interference and the steady-state solution after the cancellation weight convergence according to the broadband non-cooperative interference cancellation sampling array model;

[0072] S3. Determine performance parameters according to the transient solution and the steady-state solution, and determine the corresponding relationship between the performance parameters and the array manifold; the performance parameters include array resolution, cancellation blind area, interference cancellation ratio, system convergence time and system convergence stability;

[0073] S4. Combining the corresponding relationship between the performance parameters and the array manifold with the performance index given by the cancellation system to determine the initial sampling array manifold;

[0074] S5. Testing and verifying the initial sampling array manifold, comparing the test data with the performance index, and iteratively adjusting the array manifold until the performance index is met.

[0075] Reference Figure 2 The complete process of the broadband non-cooperative interference sampling array manifold design method of the embodiment of the present application includes:

[0076] Array Population Modeling;

[0077] Determine the number of array elements;

[0078] Solve the steady-state solution and transient solution of the weight vector;

[0079] Calculate the relationship between resolution and array size; input resolution index; determine array aperture;

[0080] Calculate the relationship between the cancellation ratio and array current; input the cancellation ratio index; determine the farthest distance of the boundary array element;

[0081] Calculate the relationship between the cancellation blind area and the array current; input the cancellation blind area index; determine the internal array element position range;

[0082] Calculate the relationship between convergence speed and array popularity; compare convergence speeds; narrow the range of internal array element positions;

[0083] Calculate the relationship between convergence stability and array popularity; compare convergence stability; determine the combination of array element positions.

[0084] The above steps are described in detail below.

[0085] The process of determining the number of array elements includes:

[0086] Obtaining a relationship between the number of interference sources and the number of array elements in a sampling array, and determining the number of array elements according to the relationship; the number of interference sources is the number of unrelated linear constraints imposed on the weight vector of the array;

[0087] The number of interference sources is determined according to the maximum number of broadband non-cooperative interference sources that need to be countered in the use scenario of the cancellation device; and the number of array elements is greater than the number of linear constraints.

[0088] To determine the number of sampling array elements, it is necessary to determine the maximum number N of broadband non-cooperative interference sources that may occur. The number of elements must be greater than N. However, considering the complexity of the backend channel, the number of sampling array elements is generally N+1.

[0089] Reference Figure 3 , Figure 3 The broadband non-cooperative interference sampling array and cancellation principle diagram of the embodiment of the present application are shown in FIG. The sampling array is modeled, and a three-dimensional rectangular coordinate system is established with the main antenna as the coordinate origin, where the plane where the main antenna is located is the xoy plane, the coordinates of the main antenna are [0,0,0], and the position coordinates of the mth array element are P M =[x m ,y m ,z m ], the delay of the mth array element relative to the main antenna is τ m ,θ j represents the interference signal incident angle, θ s Represents the incident angle of the useful signal.

[0090] Specifically, the model building process is as follows: a broadband interference signal model is established at a frequency of f l The signal is Broadband interference signals are characterized by multi-tone signals The signal received by the mth sampling array is The signal received by the sampling array is x=[x 1 ,x 2 ,…,x m ,…,x M ] H .

[0091] It includes the following sub-steps:

[0092] Sub-step S11: Establish a three-dimensional coordinate system with the main antenna position as the origin, and the positions of the M array elements are P = [P1, P2, ... P M ] T , where P M =[x m ,y m ,z m] is the position coordinate of the mth array element, and the delay of the M array elements relative to the main antenna is τ=[τ1,τ1,…,τ M ], where τ m is the delay of the mth array element relative to the main antenna, as shown in the following formula:

[0093]

[0094] Where x, y, z are the x-axis coordinates, y-axis coordinates, and z-axis coordinates of the m-th array element position, and θ is the pitch angle of the signal incident direction. is the azimuth of the signal incident direction, and c is the propagation speed of the electromagnetic wave.

[0095] Sub-step S12: Characterizing the broadband interference signal with a multi-tone signal x l The determination method is shown in the following formula:

[0096]

[0097] Among them, l represents polyphonic f l The frequency number, γ1 represents the input interference-to-noise ratio, a j represents the steering vector of the array to the interference signal. s(t) represents the useful signal, γ2 represents the input signal-to-noise ratio, a s represents the steering vector of the array for the useful signal; a l represents the amplitude of the lth tone signal, represents the phase of the lth single tone signal.

[0098] Sub-step S13: Obtain the interference signal received by the m-th array element, which can be expressed as:

[0099]

[0100] Among them, τ m is the delay of the mth array element relative to the main antenna;

[0101] The weight vector of m array elements can be expressed as:

[0102] w=[w 1 ,w 2 ,…,w m ,…,w M ] H

[0103] Wherein, the superscript H represents the conjugate transpose;

[0104] The remaining signal after cancellation can be expressed as:

[0105]

[0106] Where K is a constant, including the loop gain of each device in the channel, x H It is the sampling signal received by the sampling array.

[0107] Furthermore, the transient solution and steady-state solution solving process includes:

[0108] Determine an initial weight vector, obtain instantaneous data of the initial weight vector changing with sampling time during an iterative convergence process, and use the instantaneous data as a transient solution;

[0109] Determine a converged weight vector corresponding to the initial weight vector after convergence, and use the converged weight vector as the steady-state solution;

[0110] Obtaining an equivalent amplitude vector and an equivalent phase vector of the convergence weight vector;

[0111] The relationship between the convergence weight vector and the sampling array manifold parameter is determined according to the equivalent amplitude vector and the equivalent phase vector.

[0112] Specifically, the solution process of the weight vector steady-state solution and transient solution is as follows:

[0113] Sub-step S21: According to the above modeling, the transient difference equation during the weight convergence process is:

[0114]

[0115] Where μ is the convergence step size, is the sampling signal received by the sampling antenna, x e (n) is the residual signal.

[0116] Solve the difference equation and substitute w(0)=0 to obtain the weighted transient solution:

[0117]

[0118] in:

[0119] R=UAU H

[0120] Λ=diag{λ i}

[0121] i=1,2,3,…,M

[0122] q is the cross-correlation vector between the interference signal received by the sampling array and the interference signal received by the main antenna, expressed as:

[0123]

[0124] R is the autocorrelation matrix of the received signal of the sampling array, which is expressed as:

[0125]

[0126] Sub-step S22: After the weights converge, they can be expressed as a complex constant w m , which can be expressed as:

[0127] w m =k m exp(-jω m τ m )

[0128] The weight vector is expressed as

[0129] The residual interference signal can be expressed as

[0130] By using the residual interference signal to derivate the weight equivalent amplitude and equivalent angular frequency respectively, the optimal weight equivalent angular frequency can be obtained to satisfy

[0131] The optimal weight equivalent amplitude satisfies the following formula:

[0132]

[0133] Further, performance parameters are determined according to the transient solution and the steady-state solution, and the corresponding relationship between the performance parameters and the array manifold is determined; the performance parameters include array resolution, cancellation blind area, interference cancellation ratio, system convergence time and system convergence stability, specifically including the following steps:

[0134] Sub-step S31: Calculate the relationship between the sampling array resolution and the array manifold.

[0135] The directional pattern formed by the main antenna and the sampling array is expressed as |1-WA′|, where WA′ represents the product of the weight vector of the sampling array and the steering vector. The directional pattern is in the interference incident direction θ j The main null is formed at the point where the incident angle of the useful signal is θ s When approaching the main zero, the useful signal will be attenuated. j The minimum value that satisfies the useful signal attenuation constraint of αdB is the αdB attenuation resolution of the sampling array. Substituting the steady-state solution of the weight vector and the array steering vector into the sample array, the attenuation resolution and the maximum aperture of the array satisfy , M represents the number of sampling antennas, and r represents the spacing between the sampling antennas.

[0136] Sub-step S32: Calculate the relationship between the interference cancellation ratio and the sampling array manifold. The interference signal at each frequency point f l, which can be regarded as a vector with an amplitude of 1 and a phase of 0°. The cancellation of multiple sampling antennas can be regarded as the fitting process of multiple vectors to the vector 1∠0° at multiple frequency points. By observing the equivalent angular frequency of the steady-state weights, it can be found that the equivalent frequency ω corresponding to the optimal weight of each antenna m Similarly, the sampling array elements form a positive and negative delay distribution relative to the main antenna, and when the cumulative delay of each sampling array element is zero, the phase of the fitted cancellation signal and the interference signal can be exactly opposite. Therefore, the sampling array elements should be arranged with the main antenna as the center to form a delay sum of zero, which will improve the cancellation performance. And the interference cancellation ratio and the delay of the sampling array element from the main antenna meet relationship;

[0137] Sub-step S33: Calculate the relationship between the cancellation blind area and the sampling array manifold. During the cancellation process, the directional pattern will form a grid null except for the interference direction. When the useful signal is incident from these angles, it will cause the useful signal to be attenuated and reduce the output signal-to-noise ratio. When the incident angle of the interference signal is fixed and the incident angle of the useful signal changes, the relationship between the output signal-to-noise ratio and the array manifold satisfies the following relationship:

[0138]

[0139] in,

[0140]

[0141] Sub-step S44: Calculate the relationship between the convergence time and the sampling array manifold. Observing the transient solution of the weight vector, it can be seen that during the weight convergence process, the weight expectation is a linear combination of M exponential decay functions. The convergence time is determined by the weight with the slowest convergence speed, and the exponential decay function decays to e -1 The time required is:

[0142]

[0143] Among them, μ represents the convergence step size, λ i is the i-th eigenvalue of the autocorrelation matrix of the sampling array, K is a constant, and contains the loop gain of each device in the channel.

[0144] Obviously, the minimum eigenvalue corresponds to the exponential function with the slowest convergence speed, which is the time required for the system to converge. The eigenvalue is determined by the autocorrelation matrix of the sampling array, which is:

[0145]

[0146] Sub-step S45: Calculate the relationship between convergence stability and sampling array manifold. Since the weight update adopts a closed-loop algorithm, the solution is:

[0147]

[0148] in, For w m (n) is irrelevant, and as n increases, G(n) tends to be stable. As the iteration time n increases, if you want the system to be stable, you need to satisfy F(n) n+1 Approximately zero.

[0149]

[0150] The solution is It can be seen that the stability condition of the system is determined by the sampling array manifold. m related.

[0151] Further, the corresponding relationship between the performance parameter and the array manifold is combined with the performance index given by the cancellation system to determine the initial sampling array manifold;

[0152] In this embodiment, the initial array aperture, the farthest distance of the boundary elements and the array element position range are obtained, and the initial sampling array manifold is determined according to the array aperture, the farthest distance of the boundary elements and the array element position range.

[0153] First, the method for determining the array aperture includes:

[0154] Utilizing the steady-state solution to obtain a first relationship between array resolution and interference sampling array manifold;

[0155] Determine all interference incident angles, input the interference incident angles into the first relationship, and obtain candidate array apertures corresponding to the interference incident angles;

[0156] Determine whether the maximum array aperture corresponding to each interference incident angle meets the resolution index, and screen the candidate array apertures according to the determination result until the final array aperture is obtained;

[0157] The resolution index is the array resolution required by the cancellation device at all interference incident angles.

[0158] The maximum aperture of the array should be determined based on the broadband interference sampling array resolution |θ-θ j | min Relationship with array aperture to determine, where λ is the wavelength corresponding to the center frequency of the broadband interference signal, and Mr is the maximum aperture of the array.

[0159] Reference Figure 4 , Figure 4 θ is the relationship between the array resolution and the array aperture under different interference incident angles in this embodiment. The black line in the figure represents θ j =90°; the red line indicates θ j=60°; the blue line indicates θ j =45°; the green line indicates θ j =30°. As can be seen from the figure, the array aperture and resolution are negatively correlated.

[0160] Reference Figure 5 , Figure 5 : is a graph showing the relationship between array resolution and interference incident angle under different numbers of array elements in this embodiment. In the graph, the black line indicates M=4; the red line indicates M=6; the blue line indicates M=8; and the green line indicates M=10. It can be seen from the graph that the resolution is negatively correlated with the interference incident angle.

[0161] Reference Figure 6 , Figure 6 This is a relationship diagram of the number of selected array elements and the interference incident angle to verify the array resolution in this embodiment; the black line in the figure indicates that M=4, θ j =90°; the red line indicates M=6, θ j = 60°; the blue line indicates M = 8, θ j =45°; the green line indicates M=10, θ j =30°.

[0162] Secondly, the method for determining the farthest distance of the boundary element includes:

[0163] Using the steady-state solution to obtain a second relationship between a broadband interference cancellation ratio and an interference sampling array manifold;

[0164] Determine the farthest distance of the boundary array element according to the second relationship and the cancellation ratio index;

[0165] Reference Figure 7 , Figure 7 This is a simulation diagram of the relationship between the interference cancellation ratio and the distance between the array element and the main antenna in this embodiment. In the figure, the black line indicates bandwidth = 10MHz, the red line indicates bandwidth = 20MHz, the blue line indicates bandwidth = 40MHz, the green line indicates bandwidth = 60MHz, and the purple line indicates bandwidth = 80MHz. It can be seen from the figure that the interference cancellation ratio is negatively correlated with the distance between the array element and the main antenna.

[0166] The maximum distance between array elements is determined according to the cancellation ratio required by the system. According to the steady-state weight, the relationship between the interference cancellation ratio and the sampling array manifold can be obtained as follows:

[0167]

[0168] Among them, τ m is the delay of the mth array element relative to the main antenna, ω l is the angular frequency of the lth tone signal, ω mis the equivalent angular frequency of the mth sampling antenna weight, k m The equivalent amplitude of the mth sampling antenna weight. When the array elements are arranged with the main antenna as the center, forming a positive and negative delay distribution relative to the main antenna, and the cumulative sum of the delays of each sampling antenna is zero, the cancellation ratio performance is optimal. The cancellation ratio is negatively correlated with the delay of the array element relative to the main antenna. Given the minimum required cancellation ratio of the system, determine the maximum spacing of the array elements.

[0169] Thirdly, the method for determining the array element position range includes:

[0170] Using the steady-state solution, a third relationship between the cancellation blind area and the interference sampling array manifold is obtained;

[0171] According to the third relationship and the cancellation blind area index, the array element position range is determined when the cancellation blind area accounts for the smallest proportion.

[0172] Reference Figure 8 , Figure 8 Schematic diagram of the cancellation blind zone of this embodiment, (a) is a three-dimensional diagram of the array output signal to interference noise ratio, and (b) is a diagram of the array 3dB cancellation blind zone.

[0173] The range of array element positions should be determined so that the cancellation blind area ratio meets the system's required indicators. According to the above requirements, the array boundary element positions have been determined, and the array internal element positions need to be adjusted to reduce the cancellation blind area. The cancellation blind area is determined based on the output signal-to-interference-to-noise ratio. The relationship between the output signal-to-interference-to-noise ratio and the array manifold satisfies:

[0174]

[0175] Among them, γ1 is the input signal-to-noise ratio, γ2 is the input interference-to-noise ratio, a s is the steering vector of the useful signal, a j is the steering vector of the interference signal.

[0176] Fourthly, the method for reducing the position range of the array element includes:

[0177] Using the transient solution, obtaining a fourth relationship between the interference cancellation convergence time and the interference sampling array manifold, and a fifth relationship between the convergence stability and the interference sampling array manifold;

[0178] According to the fourth relationship and the convergence time index, and the fifth relationship and the convergence stability index, selecting an array element position combination with the shortest convergence time and the best convergence stability from the array element position range, and combining the array element positions to form a new array element position range to achieve range reduction;

[0179] The interference cancellation convergence time is the time required for the initial weight vector to converge to a steady state; and the convergence stability is the convergence stability of the weight vector taking into account the convergence step size and loop gain.

[0180] Reference Fig. 9 , Fig. 9 Schematic diagrams of the weight vector convergence speed at different array element positions of this embodiment, (a) is a weight vector phase convergence time diagram of a "zero-sum" array, (b) is a weight vector amplitude convergence time diagram of a "zero-sum" array, (c) is a weight vector phase convergence time diagram of a non-"zero-sum" array, and (d) is a weight vector amplitude convergence time diagram of a non-"zero-sum" array.

[0181] In the figure, the blue line represents auxiliary antenna 1, the red line represents auxiliary antenna 2, and the yellow line represents auxiliary antenna 3. It can be seen that in the "zero-sum" array, the relationship between the weight amplitude, weight phase and time of the three auxiliary antennas remains stable. In the "non-zero-sum" array, the weight amplitude of auxiliary antenna 2 remains stable, the weight amplitude of auxiliary antenna 1 gradually increases and then remains stable, and the weight amplitude of auxiliary antenna 3 first decreases and then gradually increases, and then remains stable. The weight phases of auxiliary antenna 1 and auxiliary antenna 2 remain stable, and the weight of auxiliary antenna 3 gradually increases and then remains stable.

[0182] Specifically, the range of array element positions is narrowed, and array element positions with better cancellation convergence time and cancellation convergence stability are selected within the determined array element range, and several sets of better array element position combinations are determined. The relationship between the cancellation convergence time and the sampling array manifold satisfies the following relationship:

[0183]

[0184] Among them, μ is the convergence step size, K is the loop gain, λ min is the minimum eigenvalue of the autocorrelation matrix of the sampling array. The relationship between convergence stability and the sampling array manifold satisfies:

[0185]

[0186] Among them, a l is the amplitude of the lth tone signal, is the phase of the lth tone signal.

[0187] Specifically, obtaining the initial array aperture, the farthest distance of the boundary array element, and the array element position range includes the following steps:

[0188] Sub-step S41: Determine system index parameters, for example, the array αdB attenuation resolution must be less than a°, the cancellation ratio must reach above bdB, the cancellation blind area range must be less than c%, and the convergence time must be less than dμs;

[0189] Sub-step S42: According to the array resolution requirement, the maximum aperture of the sampling array may be determined;

[0190] Sub-step S43: According to the requirement of broadband interference cancellation ratio, it can be determined that the sampling array elements should be placed around the main antenna, and the delay sum of all sampling array elements relative to the main antenna should be as zero as possible. The position of the boundary array element can be determined, and multiple combinations of the array element positions within the boundary can be obtained;

[0191] Sub-step S44: calculating the cancellation blind area of ​​the array under each combination, and adjusting the position of the internal array element on this basis to find the internal array element position range with a smaller cancellation blind area;

[0192] Sub-step S45: within the internal array element position range, calculate the convergence time of each array element position within the range, and determine whether the system stability is satisfied. Select the array element position combination with the shorter convergence time as the preliminary sampling array manifold.

[0193] Finally, the initial sampling array manifold is tested and verified, and the test data is compared with the performance index to iteratively adjust the array manifold until the performance index is met.

[0194] Specifically, the iterative process includes the following steps:

[0195] Sub-step S51: Arrange the sampling array of the cancellation device according to the array manifold preliminarily determined. Connect the transmitting antenna with a signal source as the useful signal and the interference signal, and incident from different angles;

[0196] Sub-step S52: measuring whether the cancellation system can converge stably under different interference signal and useful signal incident angles, and measuring the resolution, cancellation blind area, interference cancellation ratio and convergence time of the array;

[0197] Sub-step S53: Determine whether the index of the sampling array manifold meets the standard. If the index requirements are met, the array manifold can be used as the sampling array manifold of the system; if any index cannot meet the standard, the above steps need to be repeated to adjust the array element position until the conditions are met.

[0198] By inputting useful signals and interference signals at different angles, the embodiments of the present application can comprehensively evaluate the performance of the cancellation system in actual operation, and confirm through testing whether the cancellation system can stably converge in a complex signal environment and whether it can effectively suppress interference signals, thereby improving the overall system performance. And by judging whether the indicators of the sampling array manifold meet the standards and adjusting the array element positions as needed, the system optimization cycle can be effectively reduced, and possible problems can be quickly identified and resolved, thereby speeding up the deployment and application process of the system. By repeatedly testing and adjusting to ensure that the sampling array manifold meets the performance indicator requirements, the stability and reliability of the system can be improved.

[0199] Reference Fig.10Based on the method in the above embodiment, the embodiment of the present application provides an electronic device, which may include: a processor (processor) 1010, a communication interface (Communications Interface) 1020, a memory (memory) 1030 and a communication bus 1040, wherein the processor 1010, the communication interface 1020, and the memory 1030 communicate with each other through the communication bus 1040. The processor 1010 can call the logic instructions in the memory 1030 to execute the method in the above embodiment.

[0200] In addition, the logic instructions in the above-mentioned memory 1030 can be implemented in the form of software functional units and can be stored in a computer-readable storage medium when sold or used as an independent product. Based on this understanding, the technical solution of the present application, or the part that contributes to the prior art, or the part of the technical solution, can be embodied in the form of a software product, which is stored in a storage medium and includes a number of instructions for 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 method described in each embodiment of the present application.

[0201] Based on the method in the above embodiment, an embodiment of the present application provides a computer-readable storage medium, which stores a computer program. When the computer program runs on a processor, the processor executes the method in the above embodiment.

[0202] Based on the method in the above embodiment, an embodiment of the present application provides a computer program product. When the computer program product runs on a processor, the processor executes the method in the above embodiment.

[0203] It is understandable that the processor in the embodiment of the present application can be a central processing unit (CPU), or other general-purpose processors, digital signal processors (DSP), application specific integrated circuits (ASIC), field programmable gate arrays (FPGA) or other programmable logic devices, transistor logic devices, hardware components or any combination thereof. The general-purpose processor can be a microprocessor or any conventional processor.

[0204] The method steps in the embodiments of the present application can be implemented by hardware or by a processor executing software instructions. The software instructions can be composed of corresponding software modules, and the software modules can be stored in random access memory (RAM), flash memory, read-only memory (ROM), programmable read-only memory (PROM), erasable programmable read-only memory (EPROM), electrically erasable programmable read-only memory (EEPROM), registers, hard disks, mobile hard disks, CD-ROMs, or any other form of storage medium known in the art. An exemplary storage medium is coupled to a processor so that the processor can read information from the storage medium and write information to the storage medium. Of course, the storage medium can also be a component of the processor. The processor and the storage medium can be located in an ASIC.

[0205] In the above embodiments, it can be implemented in whole or in part by software, hardware, firmware or any combination thereof. When implemented using software, it can be implemented in whole or in part in the form of a computer program product. The computer program product includes one or more computer instructions. When the computer program instructions are loaded and executed on a computer, the process or function described in the embodiment of the present application is generated in whole or in part. The computer may be a general-purpose computer, a special-purpose computer, a computer network, or other programmable device. The computer instructions may be stored in a computer-readable storage medium or transmitted through the computer-readable storage medium. The computer instructions may be transmitted from a website site, a computer, a server or a data center to another website site, a computer, a server or a data center by wired (e.g., coaxial cable, optical fiber, digital subscriber line (DSL)) or wireless (e.g., infrared, wireless, microwave, etc.) means. The computer-readable storage medium may be any available medium that a computer can access or a data storage device such as a server or a data center that includes one or more available media integrated. The available medium may be a magnetic medium (e.g., a floppy disk, a hard disk, a tape), an optical medium (e.g., a DVD), or a semiconductor medium (e.g., a solid state disk (SSD)), etc.

[0206] It should be understood that the various numerical numbers involved in the embodiments of the present application are only used for the convenience of description and are not used to limit the scope of the embodiments of the present application.

[0207] It will be easily understood by those skilled in the art that the above description is only a preferred embodiment of the present application and is not intended to limit the present application. Any modifications, equivalent substitutions and improvements made within the spirit and principles of the present application shall be included in the scope of protection of the present application.

Claims

1. A broadband non-cooperative interference sampling array manifold design method, characterized in that: include: Establish a broadband non-cooperative interference sampling array model and determine the number of array elements; Calculating the transient solution of the cancellation weight convergence process under broadband interference and the steady-state solution after the cancellation weight convergence according to the broadband non-cooperative interference cancellation sampling array model; Determining performance parameters according to the transient solution and the steady-state solution, and determining a corresponding relationship between the performance parameters and the array manifold; The performance parameters include array resolution, cancellation blind area, interference cancellation ratio, system convergence time and system convergence stability; The corresponding relationship between the performance parameter and the array manifold is combined with the performance index given by the cancellation system to determine the initial sampling array manifold; The initial sampling array manifold is tested and verified, and the test data is compared with the performance index to iteratively adjust the array manifold until the performance index is met.

2. The broadband non-cooperative interference sampling array manifold design method according to claim 1 is characterized in that: The determining the number of array elements includes: Obtaining a relationship between the number of interference sources and the number of array elements in a sampling array, and determining the number of array elements according to the relationship; the number of interference sources is the number of unrelated linear constraints imposed on the weight vector of the array; The number of interference sources is determined according to the maximum number of broadband non-cooperative interference sources that need to be countered in the use scenario of the cancellation device; and the number of array elements is greater than the number of linear constraints.

3. The broadband non-cooperative interference sampling array manifold design method according to claim 1 is characterized in that: The process of obtaining the transient solution and the steady-state solution includes: Determine an initial weight vector, obtain instantaneous data of the initial weight vector changing with sampling time during an iterative convergence process, and use the instantaneous data as a transient solution; Determine a converged weight vector corresponding to the initial weight vector after convergence, and use the converged weight vector as the steady-state solution; Also includes: Obtaining an equivalent amplitude vector and an equivalent phase vector of the convergence weight vector; The relationship between the convergence weight vector and the sampling array manifold parameter is determined according to the equivalent amplitude vector and the equivalent phase vector.

4. The broadband non-cooperative interference sampling array manifold design method according to claim 3 is characterized in that: The determining of the initial sampling array manifold comprises: The initial array aperture, the farthest distance of the boundary array elements and the array element position range are obtained, and the initial sampling array manifold is determined according to the array aperture, the farthest distance of the boundary array elements and the array element position range.

5. The broadband non-cooperative interference sampling array manifold design method according to claim 4 is characterized in that: The method for determining the array aperture comprises: Utilizing the steady-state solution to obtain a first relationship between array resolution and interference sampling array manifold; Determine all interference incident angles, input the interference incident angles into the first relationship, and obtain candidate array apertures corresponding to the interference incident angles; Determine whether the maximum array aperture corresponding to each interference incident angle meets the resolution index, and screen the candidate array apertures according to the determination result until the final array aperture is obtained; The resolution index is the array resolution required by the cancellation device at all interference incident angles.

6. The broadband non-cooperative interference sampling array manifold design method according to claim 4, characterized in that: The method for determining the farthest distance of the boundary element includes: Using the steady-state solution to obtain a second relationship between a broadband interference cancellation ratio and an interference sampling array manifold; Determine the farthest distance of the boundary array element according to the second relationship and the cancellation ratio index; The wideband interference cancellation ratio is related to the signal bandwidth, and the wideband interference cancellation ratio is negatively correlated to the array element spacing.

7. The broadband non-cooperative interference sampling array manifold design method according to claim 4, characterized in that: The method for determining the array element position range includes: Using the steady-state solution, a third relationship between the cancellation blind area and the interference sampling array manifold is obtained; According to the third relationship and the cancellation blind area index, the array element position range is determined when the cancellation blind area accounts for the smallest proportion.

8. The broadband non-cooperative interference sampling array manifold design method according to claim 4, characterized in that: The method for reducing the array element position range includes: Using the transient solution, obtaining a fourth relationship between the interference cancellation convergence time and the interference sampling array manifold, and a fifth relationship between the convergence stability and the interference sampling array manifold; According to the fourth relationship and the convergence time index, and the fifth relationship and the convergence stability index, selecting an array element position combination with the shortest convergence time and the best convergence stability from the array element position range, and combining the array element positions to form a new array element position range to achieve range reduction; The interference cancellation convergence time is the time required for the initial weight vector to converge to a steady state; and the convergence stability is the convergence stability of the weight vector taking into account the convergence step size and loop gain.

9. The broadband non-cooperative interference sampling array manifold design method according to claim 1, characterized in that: The iterative adjustment of the array manifold until the performance index is met includes: The initial sampling array manifold is installed on a cancellation device, and the measurement data of the initial sampling array manifold is obtained by using the cancellation device; the measurement data includes the array resolution and cancellation blind area actually measured, and the cancellation ratio, convergence time and convergence stability of the cancellation system; Determine whether the measurement data can meet the index requirements. If all measurement data meet the performance index, determine the final sampling array manifold; if any measurement data does not meet the performance index, readjust the array manifold until all performance indicators are met.

10. An electronic device, characterized in that: include: at least one memory for storing a computer program; At least one processor is used to execute the program stored in the memory. When the program stored in the memory is executed, the processor is used to execute the method according to any one of claims 1 to 9.

Images