A channel angle estimation method, device, electronic equipment and readable storage medium
Patent Information
- Authority / Receiving Office
- CN · China
- Patent Type
- Applications(China)
- Current Assignee / Owner
- BEIJING UNIV OF POSTS & TELECOMM
- Filing Date
- 2026-03-19
- Publication Date
- 2026-06-26
Smart Images

Figure CN122293469A_ABST
Abstract
Description
Technical Field
[0001] This application relates to the field of wireless communication and channel measurement technology, and in particular to a channel angle estimation method, apparatus, electronic device and readable storage medium. Background Technology
[0002] In millimeter wave (MMW) and terahertz (THz) communication systems, the high operating frequency bands, short wavelengths, and available bandwidth enable ultra-high data rates and low latency services. However, these frequency bands exhibit characteristics such as sparse multipath propagation, high propagation loss, and sensitivity to obstruction, requiring high-gain narrow beams to maintain sufficient link budget. In this context, the angular information of the multipath components (Angle of Departure (AoD) and Angle of Arrival (AoA)) becomes a critical parameter, directly affecting beamforming and beam alignment performance, and playing a vital role in channel modeling, user localization, and environmental awareness.
[0003] In Multiple Input Multiple Output (MIMO) communication systems, acquiring accurate angular information using the spatial resolution capability of the antenna array is a core element for achieving efficient high-frequency communication. However, high-frequency system hardware is susceptible to phase noise, frequency offset, and array calibration errors. Furthermore, limited by finite radio frequency links, the system often only acquires low-dimensional power or projected signals. This makes achieving high-precision and robust AoD / AoA estimation with limited training and measurement overhead a technical challenge and one of the key directions for current millimeter-wave / terahertz MIMO channel research and engineering implementation. Summary of the Invention
[0004] The technical objective of this application is to provide a channel angle estimation method, apparatus, electronic device, and readable storage medium to solve the problem that MIMO systems in the prior art are difficult to achieve high-precision and high-robust channel angle estimation.
[0005] To address the aforementioned technical problems, embodiments of this application provide a channel angle estimation method, including:
[0006] Obtain the discrete beam set of the preset angular domain and the pointing angle, radiation pattern and received power of each beam in the discrete beam set;
[0007] Determine the center beam and obtain the power difference ratio between the center beam and the target neighboring beam, wherein the center beam is the beam corresponding to the larger value of the received power, and the target neighboring beam is one of the adjacent beams that surround the angular domain of the center beam;
[0008] The power difference ratio is mapped to an angular offset based on the radiation pattern corresponding to the center beam;
[0009] The pointing angle corresponding to the center beam is corrected based on the angle offset to obtain the estimated angle.
[0010] Specifically, in the channel angle estimation method described above, obtaining the power difference ratio between the center beam and the target adjacent beam includes:
[0011] The beam corresponding to the larger value of the received power among the adjacent beams is identified as the target adjacent beam.
[0012] The power difference ratio is obtained by calculating the ratio of the power difference between the center beam and the adjacent beam to the sum of their powers.
[0013] Furthermore, in the channel angle estimation method described above, obtaining the power difference ratio between the center beam and the target adjacent beam includes:
[0014] The ratio of the power difference between the center beam and each adjacent beam to the sum of their powers is calculated to obtain the reference power difference ratio.
[0015] Based on the radiation pattern corresponding to the center beam, the power difference ratio is determined from the reference power difference ratio, and the adjacent beam corresponding to the power difference ratio is determined as the target adjacent beam.
[0016] Specifically, in the channel angle estimation method described above, mapping the power difference ratio to an angle offset based on the radiation pattern corresponding to the center beam includes:
[0017] The power difference ratio is inverted into the angular offset in a closed loop based on the radiation pattern;
[0018] Alternatively, the power difference ratio can be looked up and interpolated on a pre-constructed monotonic table to obtain the angle offset, wherein the monotonic table is calculated based on the radiation pattern.
[0019] Preferably, the channel angle estimation method described above further includes:
[0020] The estimated angles corresponding to multiple preset angle domains are integrated to form an angle set corresponding to multiple paths.
[0021] Preferably, the channel angle estimation method described above further includes:
[0022] The estimated angle is smoothly updated and / or Kalman updated according to a time sliding window of a preset size.
[0023] Another embodiment of this application also provides a control device, including:
[0024] The first processing module is used to obtain a set of discrete beams in a preset angular domain, as well as the pointing angle, radiation pattern and received power of each beam in the set of discrete beams;
[0025] The second processing module is used to determine the center beam and obtain the power difference ratio between the center beam and the target neighboring beam. The center beam is the beam corresponding to the larger value of the received power, and the target neighboring beam is one of the adjacent beams that surround the angular domain of the center beam.
[0026] The third processing module is used to map the power difference ratio into an angular offset according to the radiation pattern corresponding to the center beam;
[0027] The fourth processing module is used to correct the pointing angle corresponding to the center beam based on the angle offset to obtain an estimated angle.
[0028] Another embodiment of this application provides an electronic device, including a processor, a memory, and a computer program stored in the memory and executable on the processor, wherein the computer program, when executed by the processor, implements the steps of the channel angle estimation method as described above.
[0029] Another embodiment of this application provides a computer-readable storage medium storing a computer program that, when executed by a processor, implements the steps of the channel angle estimation method as described above.
[0030] Another embodiment of this application provides a computer program product including computer instructions that, when executed by a processor, implement the steps of the channel angle estimation method as described above.
[0031] Compared with existing technologies, the channel angle estimation method, apparatus, electronic device, and readable storage medium provided in this application have at least the following advantages: This application utilizes the correspondence between the power difference between the center beam and adjacent beams and the radiation pattern to perform main lobe interpolation correction on the angle; without encrypted scanning / codebook amplification, the angle error can be reduced from the "step upper limit" to a smaller range within the main lobe, improving the angle resolution while maintaining the existing step size. Simultaneously, it only requires amplitude / power data to operate, does not impose strict requirements on array phase consistency, covariance estimation, and precise amplitude-phase calibration, and has low dependence on phase and strict calibration, making it suitable for horn directional scanning to detect channels. Furthermore, each candidate direction mainly involves three power readouts and one one-dimensional inversion, requiring no feature decomposition or large-scale iteration, facilitating real-time implementation on existing host computers or embedded processors. Attached Figure Description
[0032] Figure 1 This is one of the flowcharts illustrating the channel angle estimation method of this application;
[0033] Figure 2 This is the second flowchart illustrating the channel angle estimation method of this application;
[0034] Figure 3 This is the third flowchart illustrating the channel angle estimation method of this application;
[0035] Figure 4 This is a schematic diagram of the monotonic table in this application;
[0036] Figure 5 This is a schematic diagram of the control device of this application. Detailed Implementation
[0037] To make the technical problems, technical solutions, and advantages of this application clearer, a detailed description will be provided below in conjunction with the accompanying drawings and specific embodiments. In the following description, specific details such as particular configurations and components are provided merely to aid in a comprehensive understanding of the embodiments of this application. Therefore, those skilled in the art should understand that various changes and modifications can be made to the embodiments described herein without departing from the scope and spirit of this application. Furthermore, for clarity and brevity, descriptions of known functions and structures have been omitted.
[0038] It should be understood that the phrase "one embodiment" or "an embodiment" throughout the specification means that a specific feature, structure, or characteristic related to the embodiment is included in at least one embodiment of this application. Therefore, "in one embodiment" or "in an embodiment" appearing throughout the specification does not necessarily refer to the same embodiment. Furthermore, these specific features, structures, or characteristics can be combined in any suitable manner in one or more embodiments.
[0039] In the various embodiments of this application, it should be understood that the sequence number of each process described below does not imply the order of execution. The execution order of each process should be determined by its function and internal logic, and should not constitute any limitation on the implementation process of the embodiments of this application.
[0040] It should be understood that the term "and / or" in this article is merely a description of the relationship between related objects, indicating that three relationships can exist. For example, A and / or B can represent: A existing alone, A and B existing simultaneously, and B existing alone. Additionally, the character " / " in this article generally indicates that the preceding and following related objects have an "or" relationship.
[0041] In the embodiments provided in this application, it should be understood that "B corresponding to A" means that B is associated with A, and B can be determined based on A. However, it should also be understood that determining B based on A does not mean that B is determined solely based on A, but can also be determined based on A and / or other information.
[0042] See Figure 1 One embodiment of this application provides a channel angle estimation method, including:
[0043] Step S101: Obtain the discrete beam set of the preset angular domain and the pointing angle, radiation pattern and received power of each beam in the discrete beam set;
[0044] Step S102: Determine the center beam and obtain the power difference ratio between the center beam and the target neighboring beam. The center beam is the beam corresponding to the larger value of the received power, and the target neighboring beam is one of the adjacent beams that surround the angular domain of the center beam.
[0045] Step S103: Map the power difference ratio to an angle offset according to the radiation pattern corresponding to the center beam;
[0046] Step S104: Correct the pointing angle corresponding to the center beam according to the angle offset to obtain the estimated angle.
[0047] In this embodiment, when performing channel angle estimation, the discrete beam set of the preset angle domain, i.e., the target scanning angle domain, is first obtained. This discrete beam set can be represented as... ,in, Let be the number of beams. Then, based on the beams in this discrete beam set, obtain the pointing angle, radiation pattern, and received power corresponding to each beam. The set of pointing angles can be represented as... The pointing interval between adjacent beams can be expressed as At this point, the beams can be evenly spaced or non-uniformly distributed. To ensure the monotonically reversible nature of subsequent mapping, in a specific embodiment, it is preferable that the pointing interval is not less than or equal to the main lobe half-width, i.e., the half-power beamwidth (HPBW), which can be expressed as: .
[0048] Radiation patterns corresponding to each beam (with the beam center as the zero point,) The relative offset angle can be determined through one of the following: a simulation model (e.g., Gaussian / von Mises approximation), a measured discrete radiation pattern of the anechoic chamber, or a codebook directional response derived from array geometry and phase shift accuracy. This radiation pattern is only used to map the power difference into an angular offset. In one embodiment, to ensure the smooth progress of subsequent mapping, when acquiring the radiation pattern, it is preferable to have a single peak within the main lobe, approximately symmetrical, and monotonically attenuating. The sidelobe information can be retained for anomaly detection but does not participate in subsequent inversion.
[0049] Received power corresponding to each beam The beam pointing angle can be measured sequentially by directional scanning sounding (DSS), i.e., by an electric turntable driving a directional antenna to measure each beam pointing angle; or by MIMO hybrid beamforming, i.e., by measuring the beam pointing angle corresponding to the array activation codebook.
[0050] After obtaining the above parameters, the center beam is determined based on the received power of each beam, and the power difference ratio between the center beam and the target adjacent beams in the two adjacent beams on the left and right is obtained. Specifically, when determining the center beam, the received power of all beams is compared, and the beam corresponding to the maximum value is selected. The center beam is determined, and its corresponding pointing angle is... That is, the initial angle. Thus, based on the angular domain surround, the two adjacent beams to the left and right of the central beam are determined. and This method establishes beam index relationships using azimuth / elevation coordinates, rather than relying solely on beam numbers or physical locations. This ensures the availability, accuracy, and applicability of the obtained adjacent beams. Then, based on a pre-defined algorithm, one of the two adjacent beams is selected as the target neighbor beam, and the normalized power difference ratio between the center beam and the target neighbor beam is calculated. Specifically, it can be expressed as: ,in, That is, the received power corresponding to the center beam. That is, the received power corresponding to the target adjacent beam, and the power difference ratio. It is robust to noise, and amplitude normalization can effectively suppress overall gain drift.
[0051] Furthermore, based on the radiation pattern corresponding to the center beam, this power difference ratio can be mapped into an angular offset relative to the center beam within the main lobe. Then, by correcting the pointing angle corresponding to the center beam based on this angular offset, the required estimated angle can be obtained, which can be specifically expressed as: ,in, That is, the center beam The corresponding pointing angle. In another embodiment, ,in This indicates mapping the angle to the target interval (e.g.) or ).
[0052] In summary, this application utilizes the correspondence between the power difference between the central beam and adjacent beams and the radiation pattern to perform main lobe interpolation correction on the angle. Without encrypted scanning or codebook amplification, the angle error can be reduced from the "step upper limit" to a smaller range within the main lobe, improving angular resolution while maintaining the existing step size. Furthermore, it only requires amplitude / power data to operate, without strict requirements on array phase consistency, covariance estimation, and precise amplitude and phase calibration, exhibiting low dependence on phase and rigorous calibration, making it suitable for horn-based directional scanning detection channels. Moreover, each candidate direction mainly involves three power readouts and one one-dimensional inversion, requiring no eigenvalue decomposition or large-scale iterations, facilitating real-time implementation on existing host computers or embedded processors.
[0053] It should be noted that, in one specific embodiment, after obtaining the radiation patterns corresponding to each beam, amplitude normalization processing is performed on each radiation pattern, wherein the main lobe peak value is denoted as... and store it On a high-resolution grid, it serves as the basis for subsequent table lookup / interpolation.
[0054] In one embodiment, when acquiring the received power, the received power corresponding to each beam can be acquired within one time slot through a single detection. If the system supports multiple snapshots / multiple pilots, then without changing the original process, multiple measurements of the same beam are performed using a small number of snapshots, and energy averaging is used to improve the signal-to-noise ratio (SNR) and stability. Preferably, for broadband systems, the power of each subcarrier can be incoherently combined before proceeding to subsequent steps (maintaining phase independence).
[0055] In one embodiment, when determining the center beam, if there are multiple similar peaks, multiple candidate peaks can be retained, and the final center beam can be determined by refining them one by one (e.g., by interpolation algorithm).
[0056] See Figure 2 Specifically, in the channel angle estimation method described above, obtaining the power difference ratio between the center beam and the target adjacent beam includes:
[0057] Step S201: Determine the beam corresponding to the larger value of the received power among the adjacent beams as the target adjacent beam;
[0058] Step S202: Calculate the ratio of the power difference between the center beam and the adjacent beam to the sum of their powers to obtain the power difference ratio.
[0059] In this embodiment, when obtaining the power difference ratio between the center beam and the target adjacent beam, it can be based on comparing the received power of two adjacent beams to determine the beam corresponding to the larger value as the target adjacent beam. The received power of the two adjacent beams are respectively... and Specifically, this can be expressed as: , ,in, Used to record "which side of the center beam the corresponding beam position is located on". Indicates the right side. This indicates the left side. After determining the target adjacent beam, the final power difference ratio can be calculated based on the above formula for obtaining the power difference.
[0060] See Figure 3 Furthermore, in the channel angle estimation method described above, obtaining the power difference ratio between the center beam and the target adjacent beam includes:
[0061] Step S301: Calculate the ratio of the power difference between the center beam and each adjacent beam to the sum of their powers to obtain the reference power difference ratio;
[0062] Step S302: Based on the radiation pattern corresponding to the center beam, determine the power difference ratio from the reference power difference ratio, and determine the adjacent beam corresponding to the power difference ratio as the target adjacent beam.
[0063] In another embodiment, when obtaining the power difference ratio, the power difference ratio can be determined by a two-sided ratio method. Specifically, this can be done by first calculating the power difference ratio between the center beam and each adjacent beam, and recording it as the reference power difference ratio, as follows: and Then, based on the power difference ratio and the corresponding radiation pattern of the center beam, the consistency between the two adjacent beams and their radiation patterns is obtained, and the power difference ratio required for the reference power difference ratio corresponding to the side with better consistency is determined. The corresponding adjacent beam is the required target neighbor beam. This method not only ensures the accuracy of the obtained target neighbor beam and power difference ratio, but also helps to avoid symbol sensitivity in symmetrical cases and improves response speed.
[0064] Specifically, in the channel angle estimation method described above, mapping the power difference ratio to an angle offset based on the radiation pattern corresponding to the center beam includes:
[0065] The power difference ratio is inverted into the angular offset in a closed loop based on the radiation pattern;
[0066] Alternatively, the power difference ratio can be looked up and interpolated on a pre-constructed monotonic table to obtain the angle offset, wherein the monotonic table is calculated based on the radiation pattern.
[0067] In this embodiment, we first demonstrate the forward model within the main lobe, assuming a relative angular offset. , The amplitude pattern of the main lobe is given by the adjacent beam angle distance. The ideal power difference ratio between the center beam and the target adjacent beam is: Therefore, it can be seen that within the main lobe Follow Monotony (strict condition by) (The monotonicity of the single peak is guaranteed), which can be used for inversion.
[0068] Therefore In the case of resolvability, the power difference ratio can be inverted into angular offset using a closed-loop inversion method. Specifically, this can be done when... (When using the von Mises approximation) ,in The amplitude response function of the beam pattern describes the angular shift of the signal relative to the beam center. Changes in receiver gain over time; This indicates the maximum amplitude gain at the beam center direction; These are the lumped parameters of the von Mises pattern, used to characterize the convergence of the main lobe of the beam. Their values are usually determined by the half-power beamwidth. The larger the beam size, the narrower the beam and the stronger the directionality. This represents the true angular offset to be determined, which is the fine deviation of the target signal relative to the direction of the center beam; This is the direction sign factor, which takes a value of +1 or -1, and is used to distinguish whether the target is located on the left or right side of the center beam, adjacent to the beam direction. This indicates the angular spacing between the center beam and adjacent beams; The normalized power difference ratio, constructed based on the received power of the center beam and the adjacent beams, is used to reflect the relative positional relationship of the target between the two beams. Represents the natural logarithm operation. This represents the arcsine operation. The closed-form inversion described above can be performed with constant computational complexity, making it suitable for online implementation.
[0069] When using measured discrete radiation patterns, the angular offset can be determined by interpolation through a lookup table. That is, it is pre-defined in... Calculate on the grid The resulting monotonic table can be like this: Figure 4 As shown, the online power difference ratio The angle offset can be obtained by performing a binary search followed by spline / linear interpolation. .
[0070] In another embodiment, the angular offset can also be obtained through numerical solution, specifically by first constructing an objective function. Then, in Perform a one-dimensional search with a single peak, for example, by using the golden ratio to ensure convergence.
[0071] Preferably, the channel angle estimation method described above further includes:
[0072] The estimated angles corresponding to multiple preset angle domains are integrated to form an angle set corresponding to multiple paths.
[0073] In this embodiment, when it is necessary to obtain the angle estimate of the multipath, the above steps need to be performed on the preset angle domain corresponding to each path to obtain the estimated angle corresponding to each preset angle. Then, the estimated angles can be integrated to obtain the angle set corresponding to the multipath.
[0074] In another embodiment, if AoA / AoD need to be estimated separately, the same process can be performed at both the receiving and sending ends.
[0075] Preferably, the channel angle estimation method described above further includes:
[0076] The estimated angle is smoothly updated and / or Kalman updated according to a time sliding window of a preset size.
[0077] In this implementation, channel angle estimation in dynamic scenarios is exemplified. Specifically, the estimated angle of consecutive time slots is smoothly updated and / or Kalman updated according to a preset time sliding window. This can significantly reduce jitter without introducing other biases and ensure the timing tracking of angle estimation in dynamic scenarios.
[0078] See Figure 5 Another embodiment of this application also provides a control device, including:
[0079] The first processing module 501 is used to obtain a set of discrete beams in a preset angular domain and the pointing angle, radiation pattern and received power of each beam in the set of discrete beams;
[0080] The second processing module 502 is used to determine the center beam and obtain the power difference ratio between the center beam and the target neighboring beam. The center beam is the beam corresponding to the larger value of the received power, and the target neighboring beam is one of the adjacent beams that surround the angular domain of the center beam.
[0081] The third processing module 503 is used to map the power difference ratio into an angular offset according to the radiation pattern corresponding to the center beam;
[0082] The fourth processing module 504 is used to correct the pointing angle corresponding to the center beam according to the angle offset to obtain an estimated angle.
[0083] Specifically, in the control device described above, the second processing module includes:
[0084] The first processing unit is configured to determine the beam corresponding to the larger value of the received power among the adjacent beams as the target adjacent beam;
[0085] The second processing unit is used to calculate the ratio of the power difference between the center beam and the adjacent beam to the sum of their powers, and obtain the power difference ratio.
[0086] Furthermore, in the control device described above, the second processing module includes:
[0087] The third processing unit is used to calculate the ratio of the power difference between the center beam and each adjacent beam to the sum of their powers, and to obtain the reference power difference ratio.
[0088] The fourth processing unit is configured to determine the power difference ratio from the reference power difference ratio based on the radiation pattern corresponding to the center beam, and to determine the adjacent beam corresponding to the power difference ratio as the target adjacent beam.
[0089] Specifically, in the control device described above, the three processing modules include:
[0090] The fifth processing unit is used to convert the power difference ratio into the angle offset in a closed loop based on the radiation pattern.
[0091] Alternatively, the sixth processing unit is used to look up and interpolate the power difference ratio on a pre-constructed monotonic table to obtain the angle offset, wherein the monotonic table is calculated based on the radiation pattern.
[0092] Preferably, the control device described above further includes:
[0093] The fifth processing module is used to integrate the estimated angles corresponding to multiple preset angle domains to form an angle set corresponding to multiple paths.
[0094] Preferably, the control device described above further includes:
[0095] The sixth processing module is used to perform smooth updates and / or Kalman updates on the estimated angle according to a time sliding window of a preset size.
[0096] The apparatus embodiments of this application are apparatuses corresponding to the embodiments of the methods described above. All implementation means in the method embodiments described above are applicable to the apparatus embodiments and can achieve the same technical effects. The apparatus provided in this application embodiments can implement all the method steps implemented in the method embodiments described above and can achieve the same technical effects. Therefore, the parts and beneficial effects that are the same as those in the method embodiments in this embodiment will not be described in detail here.
[0097] Another embodiment of this application provides an electronic device, including a processor, a memory, and a computer program stored in the memory and executable on the processor. When the computer program is executed by the processor, it implements the steps of the channel angle estimation method described above and achieves the same technical effect. To avoid repetition, it will not be described again here.
[0098] Another embodiment of this application provides a computer-readable storage medium storing a computer program. When executed by a processor, the computer program implements the steps of the channel angle estimation method described above and achieves the same technical effect. To avoid repetition, it will not be described again here. The computer-readable storage medium may be a read-only memory (ROM), a random access memory (RAM), a magnetic disk, or an optical disk.
[0099] Another embodiment of this application provides a computer program product, including computer instructions that, when executed by a processor, implement the steps of the channel angle estimation method described above and achieve the same technical effect. To avoid repetition, it will not be described again here.
[0100] Furthermore, reference numerals and / or letters may be repeated in different examples within this application. Such repetition is for the purpose of simplification and clarity and does not in itself indicate a relationship between the various embodiments and / or settings discussed.
[0101] It should also be noted that, in this document, relational terms such as "first" and "second" are used only to distinguish one entity or operation from another, and do not necessarily require or imply any such actual relationship or order between these entities or operations. Furthermore, the terms "comprising," "including," or any other variations thereof are intended to cover non-exclusive inclusion.
[0102] The above description is the preferred embodiment of this application. It should be noted that for those skilled in the art, several improvements and modifications can be made without departing from the principles described in this application, and these improvements and modifications should also be considered within the scope of protection of this application.
Claims
1. A method of channel angle estimation, the method comprising: The method comprises: acquiring a set of discrete beams in a preset angular domain and a pointing angle, a directional diagram and a received power corresponding to each beam in the set of discrete beams; determining a center beam and acquiring a power difference ratio of the center beam and a target adjacent beam, the center beam being a beam corresponding to a larger value in the received power, and the target adjacent beam being one of adjacent beams surrounding the center beam in the angular domain; mapping the power difference ratio to an angular offset according to the directional diagram corresponding to the center beam; correcting the pointing angle corresponding to the center beam according to the angular offset to obtain an estimated angle.
2. The channel angle estimation method of claim 1, wherein, The acquiring of the power difference ratio of the center beam and the target adjacent beam comprises: determining the target adjacent beam as a beam corresponding to a larger value in the received power among the adjacent beams; calculating a ratio of a power difference between the center beam and the adjacent beam to a power sum to obtain the power difference ratio.
3. The channel angle estimation method of claim 2, wherein, The acquiring of the power difference ratio of the center beam and the target adjacent beam comprises: calculating a ratio of a power difference between the center beam and each adjacent beam to a power sum to obtain a reference power difference ratio; determining the power difference ratio from the reference power difference ratio according to the directional diagram corresponding to the center beam, and determining the adjacent beam corresponding to the power difference ratio as the target adjacent beam.
4. The channel angle estimation method of claim 1, wherein The mapping of the power difference ratio to the angular offset according to the directional diagram corresponding to the center beam comprises: inverting the power difference ratio to the angular offset according to the directional diagram; or, looking up and interpolating the power difference ratio on a pre-constructed monotonic table to obtain the angular offset, the monotonic table being calculated according to the directional diagram.
5. The channel angle estimation method of claim 1, wherein, The method further comprises: integrating the estimated angles corresponding to a plurality of preset angular domains to form an angle set corresponding to multipath.
6. The channel angle estimation method of claim 1, wherein, The method further comprises: performing smooth updating and / or Kalman updating on the estimated angles according to a preset size of a time sliding window.
7. A control device characterized by comprising: The method comprises: a first processing module, configured to acquire a set of discrete beams in a preset angular domain and a pointing angle, a directional diagram and a received power corresponding to each beam in the set of discrete beams; a second processing module, configured to determine a center beam and acquire a power difference ratio of the center beam and a target adjacent beam, the center beam being a beam corresponding to a larger value in the received power, and the target adjacent beam being one of adjacent beams surrounding the center beam in the angular domain; a third processing module, configured to map the power difference ratio to an angular offset according to the directional diagram corresponding to the center beam; a fourth processing module, configured to correct the pointing angle corresponding to the center beam according to the angular offset to obtain an estimated angle.
8. An electronic device, comprising: The method comprises a processor, a memory and a computer program stored on the memory and executable on the processor, and the computer program is executed by the processor to implement the steps of the channel angle estimation method according to any one of claims 1 to 6.
9. A computer-readable storage medium, characterized in that, A computer program is stored on the computer readable storage medium, and the computer program is executed by the processor to implement the steps of the channel angle estimation method according to any one of claims 1 to 6.
10. A computer program product, characterised in that, It includes computer instructions that, when executed by a processor, implement the steps of the channel angle estimation method as described in any one of claims 1 to 6.