CN114759959B - Phased array beam forming method for inhibiting interference between beams - Google Patents
Phased array beam forming method for inhibiting interference between beams Download PDFInfo
- Publication number
- CN114759959B CN114759959B CN202210381030.6A CN202210381030A CN114759959B CN 114759959 B CN114759959 B CN 114759959B CN 202210381030 A CN202210381030 A CN 202210381030A CN 114759959 B CN114759959 B CN 114759959B
- Authority
- CN
- China
- Prior art keywords
- representing
- user
- interference
- main lobe
- power
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 33
- 230000002401 inhibitory effect Effects 0.000 title abstract 2
- 238000005457 optimization Methods 0.000 claims abstract description 55
- 238000004891 communication Methods 0.000 claims abstract description 18
- 239000013598 vector Substances 0.000 claims description 31
- 238000005070 sampling Methods 0.000 claims description 24
- 230000005540 biological transmission Effects 0.000 claims description 13
- 230000008054 signal transmission Effects 0.000 claims description 9
- 239000000654 additive Substances 0.000 claims 1
- 230000000996 additive effect Effects 0.000 claims 1
- 238000013461 design Methods 0.000 abstract description 20
- 230000001629 suppression Effects 0.000 abstract description 7
- 230000000694 effects Effects 0.000 abstract 1
- 230000005764 inhibitory process Effects 0.000 abstract 1
- 239000011159 matrix material Substances 0.000 description 18
- 238000004088 simulation Methods 0.000 description 5
- 238000010586 diagram Methods 0.000 description 4
- 238000000342 Monte Carlo simulation Methods 0.000 description 3
- 238000000354 decomposition reaction Methods 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 2
- 230000017105 transposition Effects 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 1
- 238000003491 array Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 239000002360 explosive Substances 0.000 description 1
- 238000010295 mobile communication Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000010606 normalization Methods 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
Classifications
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04B—TRANSMISSION
- H04B7/00—Radio transmission systems, i.e. using radiation field
- H04B7/02—Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas
- H04B7/04—Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas using two or more spaced independent antennas
- H04B7/0413—MIMO systems
- H04B7/0426—Power distribution
- H04B7/043—Power distribution using best eigenmode, e.g. beam forming or beam steering
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04B—TRANSMISSION
- H04B17/00—Monitoring; Testing
- H04B17/30—Monitoring; Testing of propagation channels
- H04B17/391—Modelling the propagation channel
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04W—WIRELESS COMMUNICATION NETWORKS
- H04W24/00—Supervisory, monitoring or testing arrangements
- H04W24/02—Arrangements for optimising operational condition
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02D—CLIMATE CHANGE MITIGATION TECHNOLOGIES IN INFORMATION AND COMMUNICATION TECHNOLOGIES [ICT], I.E. INFORMATION AND COMMUNICATION TECHNOLOGIES AIMING AT THE REDUCTION OF THEIR OWN ENERGY USE
- Y02D30/00—Reducing energy consumption in communication networks
- Y02D30/70—Reducing energy consumption in communication networks in wireless communication networks
Landscapes
- Engineering & Computer Science (AREA)
- Computer Networks & Wireless Communication (AREA)
- Signal Processing (AREA)
- Power Engineering (AREA)
- Physics & Mathematics (AREA)
- Electromagnetism (AREA)
- Variable-Direction Aerials And Aerial Arrays (AREA)
Abstract
Description
技术领域Technical field
本发明涉及抑制波束间干扰的相控阵波束成形技术领域,特别是涉及一种抑制波束间干扰的相控阵波束成形方法。The present invention relates to the technical field of phased array beamforming for suppressing inter-beam interference, and in particular, to a phased array beamforming method for suppressing inter-beam interference.
背景技术Background technique
近年来随着无线通信的发展,无线电话业务一直在不断升级,移动通信系统面临着爆炸性的数据流量增长和海量设备连接,刺激了人们对将毫米波应用于无线通信的相关研究。毫米波的波长短,仅为数毫米,这使得毫米波通信能够在更小的天线尺寸条件下进行,收发端可以在有限的区域内配置大规模的天线阵列,毫米波系统主要利用该大规模天线阵列形成波束来实现高速率传输。但由于大量天线的存在,为每根天线分配专用的射频链路不仅在硬件上难以实现,还会造成高功耗以及高昂的射频链路成本,因此使用少量的射频链路通过一定数目的移相器驱动所有天线阵元的混合架构在现行的毫米波系统中广泛应用。In recent years, with the development of wireless communications, wireless phone services have been continuously upgraded. Mobile communication systems are facing explosive growth in data traffic and massive device connections, which has stimulated research on the application of millimeter waves in wireless communications. The wavelength of millimeter waves is short, only a few millimeters, which enables millimeter wave communications to be carried out with smaller antenna sizes. The transceiver end can configure a large-scale antenna array in a limited area. Millimeter wave systems mainly use this large-scale antenna Arrays form beams to enable high-rate transmission. However, due to the existence of a large number of antennas, assigning a dedicated RF link to each antenna is not only difficult to implement in hardware, but also causes high power consumption and high RF link costs. Therefore, using a small number of RF links to pass a certain number of mobile Hybrid architectures in which phase drivers drive all antenna elements are widely used in current millimeter-wave systems.
在毫米波多用户无线通信系统中,不同的数据流可以占用相同的时频资源,通过空分复用形成波束传输不同的数据。在实际工程中,不同的数据流之间会出现相互干扰。受到射频链数目的约束,需要通过移相器组成的相控阵网络对传输数据的波束进行设计,以达到抑制数据流间的干扰的目的。但由于移相器只有相位可调,相控阵中的每个元素均具有恒定模值的约束,因此如何在相控阵上进行抑制波束间干扰的波束成形设计是毫米波多用户无线通信系统的难点所在。In a millimeter-wave multi-user wireless communication system, different data streams can occupy the same time-frequency resources and form beams to transmit different data through space division multiplexing. In actual projects, there will be mutual interference between different data streams. Constrained by the number of radio frequency chains, it is necessary to design the beam for transmitting data through a phased array network composed of phase shifters to suppress interference between data streams. However, since the phase shifter only has an adjustable phase, each element in the phased array has a constant modulus constraint. Therefore, how to perform beamforming design on the phased array to suppress inter-beam interference is a key issue in millimeter-wave multi-user wireless communication systems. The hard part.
此外,在雷达系统中,通常希望通过最小化传输到噪声源和从噪声源接收到的信号的功率来抑制杂波和其他信号的干扰,因此,波束的干扰抑制设计对于提高雷达系统的性能具有重要意义。在波束方向图的特定方向上形成零陷是雷达系统中一种有效的抗干扰技术,纯相控阵网络因其馈电网络的经济性和简捷性在大型相控阵系统中备受关注,如何设计波束使其在特定方向形成零陷以达到抑制干扰的目的也是具有恒定模值约束的雷达系统的难点所在。In addition, in radar systems, it is usually desirable to suppress the interference of clutter and other signals by minimizing the power of signals transmitted to and received from the noise source. Therefore, the interference suppression design of the beam is useful for improving the performance of the radar system. Significance. Forming a null in a specific direction of the beam pattern is an effective anti-interference technology in radar systems. Pure phased array networks have attracted much attention in large phased array systems due to their economy and simplicity of the feed network. How to design the beam to form a null in a specific direction to suppress interference is also a difficulty in radar systems with constant modulus constraints.
发明内容Contents of the invention
有鉴于此,本发明的目的在于提供一种抑制波束间干扰的相控阵波束成形方法,该方法在满足相控阵的恒模约束下,通过引入权重因子,以最小化波束在干扰方向的功率和波束在主瓣中心区域功率的线性组合为优化目标,借助流形优化工具箱,优化设计了单个或多个波束。In view of this, the object of the present invention is to provide a phased array beamforming method that suppresses inter-beam interference. This method minimizes the interference of the beam in the interference direction by introducing a weighting factor while satisfying the constant mode constraint of the phased array. The linear combination of power and beam power in the main lobe center area is the optimization goal. With the help of the manifold optimization toolbox, single or multiple beams are optimally designed.
为了实现上述目的,本发明采用如下技术方案:In order to achieve the above objects, the present invention adopts the following technical solutions:
一种抑制波束间干扰的相控阵波束成形方法,所述方法包括如下步骤:A phased array beamforming method for suppressing inter-beam interference, the method includes the following steps:
步骤S1、针对一无线系统,构建以最小化波束在干扰方向的功率和波束在主瓣中心区域功率的线性组合为优化目标的优化问题;其中,所述的无线系统包括毫米波多用户无线通信系统和雷达系统,所述的优化问题满足相控阵中所有移相器的恒模约束;Step S1. For a wireless system, construct an optimization problem with the optimization goal of minimizing the linear combination of the power of the beam in the interference direction and the power of the beam in the main lobe center area; wherein the wireless system includes a millimeter-wave multi-user wireless communication system. and radar systems, the described optimization problem satisfies the constant modulus constraints of all phase shifters in the phased array;
步骤S2、通过求解步骤S1中构建的优化问题,为所述的无线系统设计单个或多个波束。Step S2: Design single or multiple beams for the wireless system by solving the optimization problem constructed in step S1.
进一步的,在所述步骤S1中,当所述无线系统为毫米波多用户无线通信系统时,首先构建其下行信号传输模型,以最小化当前用户波束在其主瓣中心区域功率和该波束在其余用户波束主瓣中心区域的功率的线性组合为优化目标构建优化问题。Further, in step S1, when the wireless system is a millimeter wave multi-user wireless communication system, its downlink signal transmission model is first constructed to minimize the power of the current user beam in its main lobe center area and the power of this beam in the rest of the The linear combination of the power in the center area of the main lobe of the user beam constructs an optimization problem for the optimization objective.
进一步的,在所述步骤S1中,若基站能够获得每个用户精确的出发角,则所述的主瓣中心区域简化为一个中心点。Further, in step S1, if the base station can obtain the accurate departure angle of each user, the main lobe center area is simplified to a center point.
进一步的,所述的下行信号传输模型具体表示为:Further, the downlink signal transmission model is specifically expressed as:
其中,yq表示用户接收到的信号;表示基站与第q个用户间的下行信道矢量,/>表示复数域,(·)H表示共轭转置运算;/>的K个列向量/>表示基站发送给K个用户的波束成形矢量;/>表示基站发送给K个用户的数据流,其满足的功率约束,Ps表示基站的发射功率,s的第q个元素[s]q表示基站发送给第q个用户的数据;ηq表示第q个用户接收到的加性白噪声,服从均值为0、方差为σ2的复高斯分布。Among them, y q represents the signal received by the user; Represents the downlink channel vector between the base station and the q-th user,/> Represents the complex number field, (·) H represents the conjugate transpose operation;/> K column vectors/> Represents the beamforming vector sent by the base station to K users;/> Represents the data stream sent by the base station to K users, which satisfies The power constraint of A complex Gaussian distribution with 0 and variance σ 2 .
进一步的,在所述的下行信号传输模型中,信道建模为:Further, in the downlink signal transmission model, the channel is modeled as:
其中,Nt表示基站天线数目,Lq表示基站与第q个用户间信道的多径总数,αq,l和φq,l分别表示第q个用户的第l条路径的复增益和出发角;a(φq,l)表示第q个用户的第l条路径的阵列导向矢量,其具体表达式为:Among them, N t represents the number of base station antennas, L q represents the total number of multipaths in the channel between the base station and the q-th user, α q,l and φ q,l represent the complex gain and departure value of the l-th path of the q-th user respectively. angle; a(φ q,l ) represents the array guidance vector of the l-th path of the q-th user, and its specific expression is:
其中,λc为载波波长,d表示天线阵元间距,(·)T表示转置运算。Among them, λ c is the carrier wavelength, d represents the antenna array element spacing, and (·) T represents the transposition operation.
进一步的,所述优化问题具体表示为:Further, the optimization problem is specifically expressed as:
其中,λ为一正实数权重因子,θk,j表示在第k个波束的主瓣中心区域的第j个采样点,k=1,2,…,K,j=1,2,…,J;wq表示FRF的第q列;表示第q个波束在其他波束的主瓣中心区域采样点的功率总和,定义为第q个波束对其余波束的干扰功率总和;M表示在第q个波束的主瓣中心区域的采样点数,/>表示第q个波束在其主瓣中心区域的M个采样点的功率总和;|·|表示取模值,[wq]n表示波束成形矢量wq的第n个元素,根据相控阵所有移相器的恒模约束,wq中的每一个元素均满足Among them, λ is a positive real weight factor, θ k,j represents the j-th sampling point in the main lobe center area of the k-th beam, k=1,2,…,K, j=1,2,…, J; w q represents the qth column of F RF ; Represents the total power of the sampling points of the qth beam in the main lobe center area of other beams, which is defined as the sum of the interference power of the qth beam to the other beams; M represents the number of sampling points in the main lobe center area of the qth beam, / > Represents the sum of power of the M sampling points of the q-th beam in the center area of its main lobe; |·| represents the modulus value, [w q ] n represents the n-th element of the beamforming vector w q , according to all phases of the phased array The constant modulus constraint of the phase shifter, each element in w q satisfies
进一步的,当基站能够获得每个用户精确的出发角,所述的主瓣中心区域简化为一个中心点时,其对应的优化问题具体表示为:Furthermore, when the base station can obtain the precise departure angle of each user and the main lobe center area is simplified to a center point, the corresponding optimization problem is specifically expressed as:
其中,θk表示第k个波束的主瓣中心点。Among them, θ k represents the main lobe center point of the k-th beam.
进一步的,当所述无线系统为雷达系统时,以最小化波束在其主瓣中心区域的功率和该波束在干扰区域的功率的线性组合为优化目标构建优化问题。Further, when the wireless system is a radar system, an optimization problem is constructed with the optimization goal of minimizing the linear combination of the power of the beam in its main lobe center area and the power of the beam in the interference area.
进一步的,当所述无线系统为雷达系统时,所述优化问题具体表示为:Further, when the wireless system is a radar system, the optimization problem is specifically expressed as:
其中,λ为一正实数权重因子,w表示波束成形矢量,θi,j表示在第i个干扰区域的第j个采样点,表示波束在干扰区域采样点的功率总和;M表示在波束的主瓣中心区域的采样点数,/>表示波束在其主瓣中心区域的M个采样点的功率总和;[w]n表示波束成形矢量w的第n个元素,根据相控阵所有移相器的恒模约束,w中的每一个元素均满足/> Among them, λ is a positive real weight factor, w represents the beamforming vector, θ i,j represents the j-th sampling point in the i-th interference area, Represents the total power of the beam's sampling points in the interference area; M represents the number of sampling points in the main lobe center area of the beam,/> Represents the sum of power of the M sampling points of the beam in the center area of its main lobe; [w] n represents the nth element of the beamforming vector w. According to the constant modulus constraints of all phase shifters of the phased array, each of w All elements are satisfied/>
进一步的,在所述步骤S2中,通过流形优化工具箱对优化问题进行求解。Further, in step S2, the optimization problem is solved through the manifold optimization toolbox.
本发明的有益效果是:The beneficial effects of the present invention are:
1、对于全连接架构下的毫米波多用户波束成形问题,本发明通过考虑基站发射向不同用户的波束间的干扰抑制建立优化问题模型,基于此提出了一种波束成形方法,该方法能够在只使用模拟波束成形而不使用数字波束成形的条件下,达到近似混合波束成形的平均和速率性能;1. Regarding the millimeter wave multi-user beamforming problem under the fully connected architecture, the present invention establishes an optimization problem model by considering the interference suppression between the beams transmitted by the base station to different users. Based on this, a beamforming method is proposed, which can only Achieve average and rate performance that approximates hybrid beamforming using analog beamforming instead of digital beamforming;
2、对于雷达系统中抑制干扰的唯相位波束设计问题,本发明给出了一种有效的波束设计方案,该方法相对于现有的唯相位波束零陷设计方案具有更好的性能。2. Regarding the problem of phase-only beam design to suppress interference in radar systems, the present invention provides an effective beam design scheme, which has better performance than the existing phase-only beam nulling design scheme.
附图说明Description of the drawings
图1和图2是本发明实施例1中使用的毫米波多用户无线通信系统模型的示意图;Figures 1 and 2 are schematic diagrams of a millimeter wave multi-user wireless communication system model used in Embodiment 1 of the present invention;
图3是基站端配备64阵元,阵元间隔为半波长,射频链路数目为4,基站获取到的为精确的出发角、用户与基站之间的传输路径总数等于1时,利用本发明实施例1提供的波束成形方法的用户平均和速率与混合波束成形、全数字波束成形的用户平均和速率比较图;Figure 3 shows that the base station is equipped with 64 array elements, the array element spacing is half a wavelength, the number of radio frequency links is 4, the base station obtains an accurate departure angle, and the total number of transmission paths between the user and the base station is equal to 1, using the present invention Comparison diagram of the average user sum rate of the beamforming method provided in Embodiment 1 and the average user sum rate of hybrid beamforming and fully digital beamforming;
图4是基站端配备64阵元,阵元间隔为半波长,射频链路数目为4,用户与基站之间的传输路径总数等于1时,利用本发明实施例1在不精确的出发角条件下设计的波束成形方法的用户平均和速率与在精确的出发角条件下设计的混合波束成形、全数字波束成形的用户平均和速率比较图;Figure 4 shows that when the base station is equipped with 64 array elements, the array element spacing is half a wavelength, the number of radio frequency links is 4, and the total number of transmission paths between the user and the base station is equal to 1, Embodiment 1 of the present invention is used under imprecise departure angle conditions. Comparison diagram of the user average sum rate of the beamforming method designed under the conditions of hybrid beamforming and full digital beamforming designed under precise departure angle conditions;
图5是基站端配备64阵元,阵元间隔为半波长,射频链路数目为4,用户与基站之间的传输路径总数等于3时,利用本发明实施例1在不精确的出发角条件下设计的波束成形方法的用户平均和速率与在精确的出发角条件下设计的混合波束成形、全数字波束成形的用户平均和速率比较图;Figure 5 shows that when the base station is equipped with 64 array elements, the array element spacing is half a wavelength, the number of radio frequency links is 4, and the total number of transmission paths between users and the base station is equal to 3, Embodiment 1 of the present invention is used under imprecise departure angle conditions. Comparison diagram of the user average sum rate of the beamforming method designed under the conditions of hybrid beamforming and full digital beamforming designed under precise departure angle conditions;
图6是雷达天线阵列配备32阵元,阵元间隔为半波长时,利用本发明实施例2设计的波束和利用文献[1]中的半定松弛算法、文献[2]中的克罗内克分解算法设计的波束的对比。Figure 6 shows a radar antenna array equipped with 32 array elements, and the array element spacing is half a wavelength, using the beam designed in Embodiment 2 of the present invention and using the semidefinite relaxation algorithm in the literature [1] and the Cronet in the literature [2] Comparison of beams designed by g-decomposition algorithm.
具体实施方式Detailed ways
为使本发明实施例的目的、技术方案和优点更加清楚,下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。In order to make the purpose, technical solutions and advantages of the embodiments of the present invention clearer, the technical solutions in the embodiments of the present invention will be clearly and completely described below in conjunction with the drawings in the embodiments of the present invention. Obviously, the described embodiments These are some embodiments of the present invention, rather than all embodiments. Based on the embodiments of the present invention, all other embodiments obtained by those of ordinary skill in the art without creative efforts fall within the scope of protection of the present invention.
实施例1Example 1
参见图1-图5,本实施例提供一种针对毫米波多用户无线通信系统的抑制波束间干扰的相控阵波束方法,该方法具体包括如下的步骤:Referring to Figures 1 to 5, this embodiment provides a phased array beam method for suppressing inter-beam interference for a millimeter wave multi-user wireless communication system. The method specifically includes the following steps:
步骤S1、构建基于毫米波多用户无线通信系统模型,其中该模型的结构如图1所示,具体的说:Step S1: Construct a millimeter-wave multi-user wireless communication system model. The structure of the model is shown in Figure 1. Specifically:
针对一个基站服务K个用户的下行通信场景,其用户均为单天线用户,基站采用全连接的混合波束成形架构,射频链路数目与用户数目相等即NRF=K,天线阵列为有Nt根天线、天线间隔为半波长的均匀线性阵列。For a downlink communication scenario where a base station serves K users, all of which are single-antenna users, the base station adopts a fully connected hybrid beamforming architecture, the number of radio frequency links is equal to the number of users, that is, N RF =K, and the antenna array has N t A uniform linear array with antennas and antenna spacing of half a wavelength.
具体的说,以第q个用户为例,其下行信号传输模型可建立为:Specifically, taking the q-th user as an example, its downlink signal transmission model can be established as:
其中,yq表示用户接收到的信号;表示基站与第q个用户间的下行信道矢量,/>表示复数域;/>表示基站发送给K个用户的数据流,其满足/>的功率约束,Ps表示基站的发射功率,s的第q个元素[s]q表示基站发送给第q个用户的数据;ηq表示第q个用户接收到的加性白噪声,服从均值为0,方差为σ2的复高斯分布,即/>表示基站的模拟波束成形矩阵,模拟波束成形在由移相器组成的相控阵网络上进行,受到移相器的恒定模值约束,即FRF中的每一个元素均满足/> 表示基站的数字波束成形矩阵。Among them, y q represents the signal received by the user; Represents the downlink channel vector between the base station and the q-th user,/> Represents a complex number field;/> Represents the data stream sent by the base station to K users, which satisfies/> The power constraint of is 0 and is a complex Gaussian distribution with variance σ 2 , that is,/> Represents the simulated beamforming matrix of the base station. The simulated beamforming is performed on a phased array network composed of phase shifters, and is constrained by the constant modulus value of the phase shifter, that is, each element in F RF satisfies/> Represents the digital beamforming matrix of the base station.
具体的说,在本实施例中,其采用的波束成形方法为相控阵波束成形而不使用数字波束成形,因此,数字波束成形矩阵设定为单位矩阵,即FBB=IK,则图1中的混合波束成形架构简化为图2中的模拟波束成形架构,上述的下行信号传输模型简化为:Specifically, in this embodiment, the beamforming method used is phased array beamforming instead of digital beamforming. Therefore, the digital beamforming matrix is set to the unit matrix, that is, F BB =I K , then figure The hybrid beamforming architecture in Figure 1 is simplified to the analog beamforming architecture in Figure 2. The above downlink signal transmission model is simplified to:
其中,FRF的K个列向量分别表示基站发送给K个用户的波束成形矢量。Among them, the K column vectors of F RF Represents the beamforming vectors sent by the base station to K users respectively.
具体的说,在本实施例中,上述的下行信号传输模型中,信道建模为:Specifically, in this embodiment, in the above downlink signal transmission model, the channel is modeled as:
其中,Lq表示基站与第q个用户间信道的多径总数,αq,l和φq,l分别表示第q个用户的第l条路径的复增益和出发角(Angle of Departure,AOD)。a(φq,l)表示第q个用户的第l条路径的阵列导向矢量,其具体表达式为:Among them, L q represents the total number of multipaths in the channel between the base station and the q-th user, α q,l and φ q,l respectively represent the complex gain and angle of departure (AOD) of the l-th path of the q-th user. ). a(φ q,l ) represents the array guidance vector of the l-th path of the q-th user, and its specific expression is:
其中λc为载波波长,d表示天线阵元间距,(·)T表示转置运算。where λ c is the carrier wavelength, d represents the antenna array element spacing, and (·) T represents the transposition operation.
步骤S2、针对步骤S1中构建的基于毫米波多用户无线通信系统模型,建立其具有恒模约束的抑制波束间干扰的波束成形模型;Step S2: For the millimeter-wave multi-user wireless communication system model constructed in step S1, establish a beamforming model with constant mode constraints that suppresses inter-beam interference;
具体的说,由于用户均为单天线单射频链用户,为最大化用户接收到的信号功率,模拟波束成形矩阵FRF形成的K个波束的主瓣中心需要分别对准K个用户复增益最大的路径(一般为视距路径)的AOD,记为但由于旁瓣的存在,波束之间不可避免地存在干扰,影响基站和用户之间的通信。Specifically, since the users are all single-antenna single-RF chain users, in order to maximize the signal power received by the users, the main lobe centers of the K beams formed by the simulated beamforming matrix F RF need to be aligned with the maximum complex gain of the K users respectively. The AOD of the path (usually a line-of-sight path) is recorded as However, due to the existence of side lobes, interference inevitably exists between beams, which affects the communication between the base station and the user.
步骤S201、根据基站否能获得每个用户精确的来构建不同的优化问题,若基站无法获得每个用户精确的/>则以最小化当前用户波束在其主瓣中心区域功率和该波束在其余用户波束主瓣中心区域的功率的线性组合为优化目标构建优化问题。Step S201: Depending on whether the base station can obtain accurate information for each user To construct different optimization problems, if the base station cannot obtain accurate information for each user/> Then, an optimization problem is constructed with the optimization goal of minimizing the linear combination of the power of the current user beam in the center area of its main lobe and the power of the beam in the center area of the main lobe of other user beams.
具体的说,在本实施例中,该步骤S201包括:Specifically, in this embodiment, step S201 includes:
为减小波束间干扰对通信质量的影响,需要对波束间的干扰进行抑制。此外,考虑到基站在波束扫描阶段可能受波束分辨率的影响获取不到每个用户精确的而只能获取到/>所在的一个宽度为波束分辨率的角度域范围Ωk,存在对/>定位的误差即波束扫描误差,所以基站在实际发射波束时波束中心对准的是角度域范围Ωk的中心,而非精确的/>精确的/>可能分布在Ωk的任何位置。因此为抑制发射向不同用户的波束之间的干扰,同时保证每个用户的通信质量,需要使单个波束在主瓣偏移尽可能小的条件下在其余波束的Ωk区域形成零陷。多个用户的波束的设计问题为相互独立的K个子问题。In order to reduce the impact of inter-beam interference on communication quality, it is necessary to suppress the inter-beam interference. In addition, considering that the base station may be affected by the beam resolution during the beam scanning stage, it cannot obtain accurate information for each user. And can only get/> Located in an angular domain range Ω k with a width of beam resolution, there is a pair/> The positioning error is the beam scanning error, so when the base station actually transmits the beam, the beam center is aligned with the center of the angular domain range Ω k , rather than the precise/> precise/> may be distributed anywhere in Ω k . Therefore, in order to suppress interference between beams transmitted to different users while ensuring the communication quality of each user, it is necessary to make a single beam form a null in the Ω k region of the remaining beams under the condition that the main lobe offset is as small as possible. The design problem of beams for multiple users is K sub-problems that are independent of each other.
更具体的说,在本实施例中,以发射向第q个用户的波束的设计为例,波束成形矢量wq=FRF(:,q)的设计问题可建立为波束区域零陷的优化问题:More specifically, in this embodiment, taking the design of a beam transmitted to the q-th user as an example, the design problem of the beamforming vector w q =F RF (:,q) can be established as the optimization of the null notch in the beam area question:
其中,λ为一正实数权重因子,θk,j表示在第k个波束的主瓣中心区域的第j个采样点,k=1,2,…,K,j=1,2,…,J;wq表示FRF的第q列;表示第q个用户的波束分配在其他用户波束主瓣中心区域Ωk的采样点的功率总和,即第q个用户的波束对其余用户波束造成的干扰功率总和;M表示在待设计的当前波束即第q个波束的主瓣中心区域的采样点数,/>表示第q个波束在其主瓣中心区域的M个采样点的功率总和;|·|表示取模值,[wq]n表示波束成形矢量wq的第n个元素,根据相控阵所有移相器的恒模约束,wq中的每一个元素均满足/> Among them, λ is a positive real weight factor, θ k,j represents the j-th sampling point in the main lobe center area of the k-th beam, k=1,2,…,K, j=1,2,…, J; w q represents the qth column of F RF ; Indicates the total power of the sampling points where the beam of the qth user is allocated in the main lobe center area Ω k of other users' beams, that is, the total interference power caused by the beam of the qth user to the beams of other users; M represents the current beam to be designed That is, the number of sampling points in the main lobe center area of the q-th beam,/> Represents the sum of power of the M sampling points of the q-th beam in the center area of its main lobe; |·| represents the modulus value, [w q ] n represents the n-th element of the beamforming vector w q , according to all phases of the phased array The constant modulus constraint of the phase shifter, each element in w q satisfies/>
步骤S202、若基站能获得每个用户精确的则无需对每个波束的主瓣中心区域区域Ωk进行采样,上述优化问题中的主瓣中心区域简化为一个中心点,因此,步骤S201中构建的波束区域零陷的优化问题可以简化为多点零陷的优化问题,具体表示为:Step S202: If the base station can obtain accurate information of each user Then there is no need to sample the main lobe center area Ω k of each beam. The main lobe center area in the above optimization problem is simplified to a center point. Therefore, the optimization problem of the beam area null constructed in step S201 can be simplified to multiple The optimization problem of point zero depression is specifically expressed as:
其中,θk表示第k个波束的主瓣中心点。Among them, θ k represents the main lobe center point of the k-th beam.
步骤S203、针对上述构建的波束区域零陷的优化问题或者多点零陷的优化问题进行求解,其包括:设定权重因子λ,通过流形优化工具箱,对优化问题进行求解可得到波束成形矢量wq。以此类推,可依次解出wk(1≤k≤K),并得到模拟波束成形矩阵FRF,其中FRF(:,k)=wk。Step S203: Solve the optimization problem of the beam area nulling or the optimization problem of multi-point nulling constructed above, which includes: setting the weight factor λ, and solving the optimization problem through the manifold optimization toolbox to obtain the beam forming Vector w q . By analogy, w k (1≤k≤K) can be solved in sequence, and the simulated beamforming matrix F RF can be obtained, where F RF (:,k)=w k .
实施例2Example 2
参见图6,本实施例提供一种针对雷达系统的抑制波束间干扰的相控阵波束设计方法,该方法基于一个阵元间隔为半波长的雷达天线阵列,其阵元数目为Nt,阵列排布为均匀线性阵列。该方法具体包括如下步骤:Referring to Figure 6, this embodiment provides a phased array beam design method for suppressing inter-beam interference in a radar system. This method is based on a radar antenna array with an element spacing of half a wavelength. The number of array elements is N t . Arranged in a uniform linear array. The method specifically includes the following steps:
步骤S1、确定该雷达系统,构建其波束方向图表达式,具体包括:Step S1: Determine the radar system and construct its beam pattern expression, specifically including:
步骤S101、针对该均匀线性阵列,其在θ方向的阵列导向矢量表示为:Step S101. For the uniform linear array, its array steering vector in the θ direction is expressed as:
步骤S102、设波束成形矢量为w,则其对应的波束方向图可以表示为:Step S102. Assume the beamforming vector is w, then its corresponding beam pattern can be expressed as:
f(θ)=|wHa(θ)|2 f(θ)=|w H a(θ)| 2
其中,|·|表示取模值。Among them, |·| represents the modulus value.
步骤S2、针对步骤S1中构建的波束方向图,对其进行设计,该设计包括:以最小化波束在其主瓣中心区域的功率和该波束在干扰区域的功率的线性组合为优化目标构建优化问题,再求解该优化问题,得到波束成形矢量。Step S2: Design the beam pattern constructed in step S1. The design includes: minimizing the linear combination of the power of the beam in its main lobe center area and the power of the beam in the interference area as the optimization goal. problem, and then solve the optimization problem to obtain the beamforming vector.
具体的说,在本实施例中,该步骤S2具体包括:Specifically, in this embodiment, step S2 specifically includes:
步骤S201、为了抑制特定方向的干扰,需要对波束进行设计,使波束方向图在特定方向形成零陷,该问题可以转化为最小化波束分配在特定方向的功率。同时,为了尽可能减小主瓣的偏移,还需要保证波束在主瓣中心区域的功率最大。波束成形矢量w的设计问题可建立为如下的优化问题:Step S201. In order to suppress interference in a specific direction, the beam needs to be designed so that the beam pattern forms a null in a specific direction. This problem can be transformed into minimizing the power allocated to the beam in a specific direction. At the same time, in order to minimize the offset of the main lobe, it is also necessary to ensure that the power of the beam is maximized in the center area of the main lobe. The design problem of the beamforming vector w can be established as the following optimization problem:
其中,λ为一正实数权重因子,θi,j表示在第i个干扰区域的第j个采样点,i=1,2,…,I,表示波束在干扰区域采样点的功率总和。M表示在波束的主瓣中心区域的采样点数,/>表示波束在其主瓣中心区域的M个采样点的功率总和。|·|表示取模值,[w]n表示波束成形矢量w的第n个元素,根据相控阵所有移相器的恒模约束,w中的每一个元素均满足/> Among them, λ is a positive real weight factor, θ i,j represents the j-th sampling point in the i-th interference area, i=1,2,...,I, Indicates the total power of the beam at the sampling points in the interference area. M represents the number of sampling points in the center area of the main lobe of the beam,/> Represents the total power of M sample points of the beam in the center area of its main lobe. |·| represents the modulus value, [w] n represents the nth element of the beamforming vector w. According to the constant modulus constraints of all phase shifters in the phased array, each element in w satisfies/>
步骤S202、设定权重因子λ,通过流形优化工具箱对优化问题进行求解可得到波束成形矢量w。Step S202: Set the weight factor λ, and solve the optimization problem through the manifold optimization toolbox to obtain the beamforming vector w.
为了验证上述的实施例1和实施例2中方法的正确性和先进性,因此进行了仿真实验,具体包括:In order to verify the correctness and advancement of the methods in Example 1 and Example 2 above, simulation experiments were conducted, specifically including:
图3仿真参数为:基站天线数Nt为64,射频链路数NRF为4,基站一共服务K=4个用户。用户与基站之间的传输路径总数Lk等于1,仅包含1条主径,主径路径增益服从复高斯分布,即基站端的信道状态信息为精确的路径AOD,权重因子λ=1000。图3中,首先利用实施例1中的步骤S202建立优化问题模型,再通过流形优化工具箱得到每个用户的波束成形矢量,形成模拟波束成形矩阵。结合实际信道,改变传输信噪比,进行2000次蒙特卡罗仿真,绘出用户平均和速率与SNR的关系曲线,如图3中圆形实线所示。同时,画出混合波束成形条件下的用户平均和速率与SNR的关系曲线,如图3中星号实线所示,以及全数字波束成形条件下的用户平均和速率与SNR的关系曲线,如图3中虚线所示。用户平均和速率的表达式为The simulation parameters in Figure 3 are: the number of base station antennas N t is 64, the number of radio frequency links N RF is 4, and the base station serves a total of K = 4 users. The total number of transmission paths L k between the user and the base station is equal to 1, including only 1 main path, and the main path gain obeys the complex Gaussian distribution, that is The channel state information at the base station is the precise path AOD, and the weight factor λ=1000. In Figure 3, step S202 in Embodiment 1 is first used to establish an optimization problem model, and then the beamforming vector of each user is obtained through the manifold optimization toolbox to form a simulated beamforming matrix. Combined with the actual channel, the transmission signal-to-noise ratio was changed, 2000 Monte Carlo simulations were performed, and the relationship curve between the user's average sum rate and SNR was drawn, as shown by the circular solid line in Figure 3. At the same time, draw the relationship curve between the user average sum rate and SNR under hybrid beamforming conditions, as shown by the asterisk solid line in Figure 3, and the relationship curve between the user average sum rate and SNR under fully digital beamforming conditions, as shown Shown by the dotted line in Figure 3. The expression for the average sum rate of users is
其中Rk表示第k个用户的可达速率,其具体表达式为Among them, R k represents the reachable rate of the k-th user, and its specific expression is
其中,F=FRFFBB,记混合波束成形、全数字波束成形和本发明提出的相控阵波束成形的F分别为F2、F3和F1。相控阵波束成形的F1=FRF;混合波束成形即星号实线的由对等效信道矩阵Heq=HFRF进行迫零预编码得到,即FBB=(Heq)-1=(HFRF)-1,其中H=[h1,h2,…,hK]H,初步的混合波束成形矩阵FHB=FRFFBB,由于混合波束成形并不具有功率增益,因此还需对FHB每一列做能量归一化的处理,得到最终的F2,即/> 全数字波束成形的/>为迫零预编码器,F3=HH(HHH)-1。对比上述三条曲线可以发现,本发明提出的波束成形方法可以达到与混合波束成形相同的性能,并且与全数字波束成形相比,二者的用户平均和速率在SNR为15dB时只相差了0.0457bps/Hz。这是因为本发明提出的相控阵波束成形方法在对用户的波束进行设计时,使得当前用户的波束分配在其他用户波束中心位置的能量尽可能接近零,等效信道矩阵Heq=HFRF非对角元素的模值接近0,完成了用户间干扰抑制的任务,因此能够达到与混合波束成形相同的性能。Among them, F=F RF F BB , and F for hybrid beamforming, full digital beamforming and phased array beamforming proposed by the present invention are respectively F 2 , F 3 and F 1 . F 1 =F RF for phased array beamforming; hybrid beamforming is the asterisk solid line It is obtained by performing zero-forcing precoding on the equivalent channel matrix Heq = HF RF , that is, F BB = ( Heq ) -1 = (HF RF ) -1 , where H = [h 1 ,h 2 ,...,h K ] H , preliminary hybrid beamforming matrix F HB =F RF F BB . Since hybrid beamforming does not have power gain, it is necessary to perform energy normalization processing on each column of F HB to obtain the final F 2 , that is /> Fully digital beamforming/> For a zero-forcing precoder, F 3 = HH (HH H ) -1 . Comparing the above three curves, it can be found that the beamforming method proposed by the present invention can achieve the same performance as hybrid beamforming, and compared with all-digital beamforming, the average user sum rate of the two only differs by 0.0457bps when the SNR is 15dB. /Hz. This is because the phased array beamforming method proposed by the present invention designs the user's beam so that the energy allocated to the current user's beam at the center position of other users' beams is as close to zero as possible, and the equivalent channel matrix Heq = HF RF The modulus value of the non-diagonal elements is close to 0, which completes the task of inter-user interference suppression, so it can achieve the same performance as hybrid beamforming.
图4仿真参数为:基站天线数Nt为64,射频链路数NRF为4,基站共服务K=4个用户。用户与基站之间的传输路径总数Lk等于1,主径的路径增益服从复高斯分布,即基站端的信道状态信息为不精确的路径出发角,权重因子λ=1000。图4中,首先利用实施例1中的步骤S201建立优化问题模型,再通过流形优化工具箱得到每个用户的波束成形矢量,形成相控阵波束成形矩阵。结合实际信道,改变传输信噪比,进行2000次蒙特卡罗仿真,绘出用户平均和速率与SNR的关系曲线,如图4中圆形实线所示。同时,画出混合波束成形条件下的用户平均和速率与SNR的关系曲线,如图4中星号实线所示,以及全数字波束成形条件下的用户平均和速率与SNR的关系曲线,如图4中虚线所示。其中,为表明性能上界,在计算混合波束成形矩阵和全数字波束成形矩阵时使用的H为精确的信道状态信息即精确的路径出发角。对比上述三条曲线可以发现,本发明提出的波束成形方法可以在SNR<20dB时达到与混合波束成形相同的性能,并且与全数字波束成形相比,二者的用户平均和速率在SNR为15dB时只相差了1.13bps/Hz。这是因为本发明提出的相控阵波束成形方法在对用户的波束进行设计时,使得当前用户的波束分配在其他用户波束中心区域的能量尽可能小,等效信道矩阵Heq=HFRF接近对角矩阵,起到了干扰抑制的作用,因此能够在SNR较小即噪声为影响和速率的主要因素时达到与混合波束成形相同的性能。The simulation parameters in Figure 4 are: the number of base station antennas N t is 64, the number of radio frequency links N RF is 4, and the base station serves a total of K = 4 users. The total number of transmission paths L k between the user and the base station is equal to 1, and the path gain of the main path obeys the complex Gaussian distribution, that is The channel state information at the base station is an inaccurate path departure angle, and the weight factor λ=1000. In Figure 4, step S201 in Embodiment 1 is first used to establish an optimization problem model, and then the beamforming vector of each user is obtained through the manifold optimization toolbox to form a phased array beamforming matrix. Combined with the actual channel, the transmission signal-to-noise ratio is changed, and 2000 Monte Carlo simulations are performed to draw the relationship curve between the user's average sum rate and SNR, as shown by the circular solid line in Figure 4. At the same time, draw the relationship curve between the user average sum rate and SNR under hybrid beamforming conditions, as shown by the asterisk solid line in Figure 4, and the relationship curve between the user average sum rate and SNR under fully digital beamforming conditions, as shown Shown by the dotted line in Figure 4. Among them, in order to indicate the performance upper bound, H used when calculating the hybrid beamforming matrix and the fully digital beamforming matrix is the precise channel state information, that is, the precise path departure angle. Comparing the above three curves, it can be found that the beamforming method proposed by the present invention can achieve the same performance as hybrid beamforming when SNR is less than 20dB, and compared with full digital beamforming, the average user rate of the two is when SNR is 15dB. The difference is only 1.13bps/Hz. This is because the phased array beamforming method proposed by the present invention designs the user's beam so that the energy allocated to the current user's beam in the center area of other users' beams is as small as possible, and the equivalent channel matrix Heq = HF RF is close to The diagonal matrix plays the role of interference suppression, so it can achieve the same performance as hybrid beamforming when the SNR is small, that is, when noise is the main factor affecting the sum rate.
图5仿真参数为:基站天线数Nt为64,射频链路数NRF为4,基站共服务K=4个用户。用户与基站之间的传输路径总数Lk等于3,包含1条主径、2条从径,其中主径的路径增益服从复高斯分布,即从径的路径增益也服从复高斯分布,且能量为主径的1/100,即/>基站端的信道状态信息为不精确的路径出发角,权重因子λ=1000。图5中,首先利用实施例1中的步骤S201建立优化问题模型,通过流形优化工具箱得到每个用户的波束成形矢量,形成模拟波束成形矩阵。结合实际信道,改变传输信噪比,进行2000次蒙特卡罗仿真,绘出用户平均和速率与SNR的关系曲线,如图5中圆形实线所示。同时,画出混合波束成形条件下的用户平均和速率与SNR的关系曲线,如图5中星号实线所示,以及全数字波束成形条件下的用户平均和速率与SNR的关系曲线,如图5中虚线所示。其中,为表明性能上界,在计算混合波束成形矩阵和全数字波束成形矩阵时使用的H为精确的信道状态信息即精确的路径出发角。对比上述三条曲线可以发现,由于受到2条从径的干扰,本发明提出的相控阵波束成形的性能低于结合了迫零预编码的混合波束成形的性能,但仍可以在SNR<10dB时达到与混合波束成形相同的性能。这是因为本发明提出的相控阵波束成形方法在对用户的波束进行设计时,使得当前用户的波束分配在其他用户波束中心区域的能量尽可能小,起到了干扰抑制的作用。The simulation parameters in Figure 5 are: the number of base station antennas N t is 64, the number of radio frequency links N RF is 4, and the base station serves a total of K = 4 users. The total number of transmission paths L k between the user and the base station is equal to 3, including 1 main path and 2 slave paths. The path gain of the main path obeys the complex Gaussian distribution, that is The path gain of the slave path also obeys the complex Gaussian distribution, and the energy is 1/100 of the main path, that is/> The channel state information at the base station is an inaccurate path departure angle, and the weight factor λ=1000. In Figure 5, step S201 in Embodiment 1 is first used to establish an optimization problem model, and the beamforming vector of each user is obtained through the manifold optimization toolbox to form a simulated beamforming matrix. Combined with the actual channel, the transmission signal-to-noise ratio was changed, 2000 Monte Carlo simulations were performed, and the relationship curve between the user's average sum rate and SNR was drawn, as shown by the circular solid line in Figure 5. At the same time, draw the relationship curve between the user average sum rate and SNR under hybrid beamforming conditions, as shown by the asterisk solid line in Figure 5, and the relationship curve between the user average sum rate and SNR under fully digital beamforming conditions, as shown Shown by the dotted line in Figure 5. Among them, in order to indicate the performance upper bound, H used when calculating the hybrid beamforming matrix and the fully digital beamforming matrix is the precise channel state information, that is, the precise path departure angle. Comparing the above three curves, it can be found that due to the interference from the two slave paths, the performance of the phased array beamforming proposed by the present invention is lower than the performance of the hybrid beamforming combined with zero-forcing precoding, but it can still be achieved when SNR<10dB Achieve the same performance as hybrid beamforming. This is because the phased array beamforming method proposed by the present invention, when designing the user's beam, makes the energy allocated by the current user's beam in the center area of other users' beams as small as possible, which plays the role of interference suppression.
图6仿真参数为:雷达均匀线性阵列的天线数目Nt=32,阵元间隔为半波长。波束中心角度θ0=0°,主瓣中心区域对应的区间为[-1.7°,1.7°],采样点数M=5。干扰区域即零陷区间为[-10.8°,-14.45°]∪[14.9°,18.2°],即I=2,每个干扰区域的采样点数为J=10,权重因子λ=5000。图6中,首先利用实施例2中的步骤S201建立优化问题模型,通过流形优化工具箱解出波束成形矢量w,绘出w对应的波束,如图6中实线所示。同时,画出文献[1]的半定松弛算法设计的波束、文献[2]中的克罗内克分解算法设计的波束以及没有进行波束零陷设计的准静态波束,如图6中点划线、点线和虚线所示。对比上述4条曲线可以发现,在零陷区间,本发明设计的波束的阵列增益均达到了-60dB以下,低于文献[1]和文献[2]的阵列增益,且与准静态波束在该区间的阵列增益峰值相比下降了40dB左右;在波束中心位置,本发明设计的波束与准静态波束相比,阵列增益只下降了1.16dB,与文献[1]设计的波束相比,阵列增益只相差了0.26dB,而与文献[2]设计的波束相比,阵列增益要高出10.06dB。由此可以看出,本发明设计的波束能够在保证主瓣中心区域阵列增益的前提下,取得最优的干扰抑制的性能。The simulation parameters in Figure 6 are: the number of antennas in the radar uniform linear array N t =32, and the array element spacing is half a wavelength. The beam center angle θ 0 =0°, the interval corresponding to the main lobe center area is [-1.7°, 1.7°], and the number of sampling points M=5. The interference area, that is, the null interval is [-10.8°, -14.45°] ∪ [14.9°, 18.2°], that is, I = 2. The number of sampling points in each interference area is J = 10, and the weight factor λ = 5000. In Figure 6, step S201 in Embodiment 2 is first used to establish an optimization problem model, the beamforming vector w is solved through the manifold optimization toolbox, and the beam corresponding to w is drawn, as shown by the solid line in Figure 6. At the same time, draw the beam designed by the semidefinite relaxation algorithm in literature [1], the beam designed by the Kronecker decomposition algorithm in literature [2], and the quasi-static beam without beam null design, as shown in the dotted line in Figure 6 Shown as lines, dotted lines and dashed lines. Comparing the above four curves, it can be found that in the null interval, the array gain of the beam designed in the present invention reaches below -60dB, which is lower than the array gain of literature [1] and literature [2], and is different from the quasi-static beam in this range. The array gain peak value in the interval has dropped by about 40dB; at the beam center position, compared with the quasi-static beam, the array gain of the beam designed in the present invention has only dropped by 1.16dB. Compared with the beam designed in the literature [1], the array gain The difference is only 0.26dB, and compared with the beam designed in literature [2], the array gain is 10.06dB higher. It can be seen from this that the beam designed in the present invention can achieve optimal interference suppression performance on the premise of ensuring the array gain in the main lobe center region.
上述的文献[1]为:P.J.Kajenski.Phase Only Antenna Pattern Notching Viaa Semidefinite Programming Relaxation[J].IEEE Transactions on Antennas andPropagation,2012,60(5):2562-2565.The above-mentioned documents [1] are: P.J.Kajenski.Phase Only Antenna Pattern Notching Viaa Semidefinite Programming Relaxation[J]. IEEE Transactions on Antennas andPropagation, 2012, 60(5):2562-2565.
上述的文献[2]为:Gu T,Zhang X,He Z,et al.Phase-Only Nulling forUniform Linear Array via Kronecker Decomposition[A].In:2021XXXIVth GeneralAssembly and Scientific Symposium of the International Union of Radio Science(URSI GASS)[C].2021,pp.1-4.The above-mentioned documents [2] are: Gu T, Zhang X, He Z, et al. Phase-Only Nulling for Uniform Linear Array via Kronecker Decomposition [A]. In: 2021 GASS)[C].2021,pp.1-4.
综上所述,本发明给出的抑制波束间干扰的相控阵波束成形方法,能够在只使用模拟波束成形而不使用数字波束成形的条件下,达到近似混合波束成形的平均和速率性能,对用户波束间的干扰起到了有效的抑制作用。同时,在雷达系统中对单个波束进行设计时,本发明给出的波束零陷设计方案比起现有方法具有更好的性能。In summary, the phased array beamforming method for suppressing inter-beam interference provided by the present invention can achieve the average sum rate performance of approximate hybrid beamforming under the condition that only analog beamforming is used instead of digital beamforming. It effectively suppresses interference between user beams. At the same time, when designing a single beam in a radar system, the beam null design scheme provided by the present invention has better performance than the existing method.
本发明未详述之处,均为本领域技术人员的公知技术。Everything that is not described in detail in the present invention is a well-known technology for those skilled in the art.
以上详细描述了本发明的较佳具体实施例。应当理解,本领域的普通技术人员无需创造性劳动就可以根据本发明的构思作出诸多修改和变化。因此,凡本技术领域中技术人员依本发明的构思在现有技术的基础上通过逻辑分析、推理或者有限的实验可以得到的技术方案,皆应在由权利要求书所确定的保护范围内。The preferred embodiments of the present invention are described in detail above. It should be understood that those skilled in the art can make many modifications and changes based on the concept of the present invention without creative efforts. Therefore, any technical solutions that can be obtained by those skilled in the art through logical analysis, reasoning or limited experiments based on the concept of the present invention and on the basis of the prior art should be within the scope of protection determined by the claims.
Claims (4)
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210381030.6A CN114759959B (en) | 2022-04-12 | 2022-04-12 | Phased array beam forming method for inhibiting interference between beams |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210381030.6A CN114759959B (en) | 2022-04-12 | 2022-04-12 | Phased array beam forming method for inhibiting interference between beams |
Publications (2)
Publication Number | Publication Date |
---|---|
CN114759959A CN114759959A (en) | 2022-07-15 |
CN114759959B true CN114759959B (en) | 2024-01-30 |
Family
ID=82329088
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202210381030.6A Active CN114759959B (en) | 2022-04-12 | 2022-04-12 | Phased array beam forming method for inhibiting interference between beams |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN114759959B (en) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN116015375A (en) * | 2022-12-12 | 2023-04-25 | 东南大学 | Millimeter wave array communication multi-user beam tracking method based on auxiliary beam |
CN118054825B (en) * | 2024-02-21 | 2024-10-18 | 浙江大学 | A multi-arm beam design method based on full antenna method |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2018094565A1 (en) * | 2016-11-22 | 2018-05-31 | 深圳大学 | Method and device for beamforming under pulse noise |
CN110365388A (en) * | 2019-07-31 | 2019-10-22 | 东南大学 | A low-complexity millimeter-wave multicast beamforming method |
CN110988816A (en) * | 2019-11-12 | 2020-04-10 | 西安电子科技大学 | Coherent MIMO radar transmitting waveform design method based on space-time phase coding |
CN113721198A (en) * | 2021-09-09 | 2021-11-30 | 哈尔滨工程大学 | Physical layer security combined beam forming method for dual-function MIMO radar communication system |
-
2022
- 2022-04-12 CN CN202210381030.6A patent/CN114759959B/en active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2018094565A1 (en) * | 2016-11-22 | 2018-05-31 | 深圳大学 | Method and device for beamforming under pulse noise |
CN110365388A (en) * | 2019-07-31 | 2019-10-22 | 东南大学 | A low-complexity millimeter-wave multicast beamforming method |
CN110988816A (en) * | 2019-11-12 | 2020-04-10 | 西安电子科技大学 | Coherent MIMO radar transmitting waveform design method based on space-time phase coding |
CN113721198A (en) * | 2021-09-09 | 2021-11-30 | 哈尔滨工程大学 | Physical layer security combined beam forming method for dual-function MIMO radar communication system |
Also Published As
Publication number | Publication date |
---|---|
CN114759959A (en) | 2022-07-15 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Li et al. | Beyond diagonal reconfigurable intelligent surfaces: A multi-sector mode enabling highly directional full-space wireless coverage | |
CN111294096A (en) | A Channel Capacity Optimization Method of Smart Reflector MISO Wireless Communication System | |
CN114759959B (en) | Phased array beam forming method for inhibiting interference between beams | |
Liu et al. | Intelligent reflecting surface aided NOMA for millimeter-wave massive MIMO with lens antenna array | |
CN104617994B (en) | A kind of 3D beam-forming methods based on horizontal and vertical combined optimization | |
CN114095318A (en) | A smart metasurface-assisted channel estimation method for a millimeter-wave communication system with a hybrid configuration | |
Wang et al. | Electromagnetic-compliant channel modeling and performance evaluation for holographic MIMO | |
CN117240330A (en) | Multi-ISAC user terminal transmission precoding method based on MIMO radar and communication | |
US11956058B1 (en) | Method and system for mobile device signal to interference plus noise ratio (SINR) improvement via polarization adjusting/optimization | |
Ding et al. | Secure full-duplex communication via movable antennas | |
Kang et al. | Smart antenna for emerging 5G and application | |
Enahoro et al. | Massive multiple-input multiple-output antenna architecture for multiband 5G adaptive beamforming applications | |
CN116260496A (en) | Beamforming method for communication-aware integrated system assisted by intelligent reflective surface | |
CN110233649B (en) | Dynamic subarray design method based on millimeter wave MIMO system | |
CN114337902B (en) | IRS-assisted millimeter wave multi-cell interference suppression method | |
CN113595604B (en) | Multi-user millimeter wave communication beam forming method under partial connection architecture | |
CN116054894B (en) | A hybrid precoding method to suppress beam squint effects | |
Reichardt et al. | Capacity maximizing MIMO antenna design for Car-to-car communication | |
CN112636800B (en) | Hybrid precoding method based on mmWave massive MIMO multi-user scenario | |
Chopra et al. | Capacity analysis of hybrid MIMO using sparse signal processing in mmW 5G heterogeneous wireless networks | |
CN115694582A (en) | Robust Transmit Beamforming Method Based on NOMA-DFRC System | |
CN114584994A (en) | Beam forming and power distribution method and device for satellite-ground integrated communication network | |
Agarwal | Fairness driven TDMA with joint phase and beamforming optimization for ND-IRS assisted MU-MISO communication systems | |
Xiang et al. | Cooperative double-IRS assisted integrated sensing and communication under NLoS conditions | |
CHEN et al. | Joint Beamforming Design for Dual-Functional Radar-Communication Systems Under Beampattern Gain Constraints |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |