CN112600592A - 巴特勒矩阵相位加权优化方法及巴特勒矩阵 - Google Patents

巴特勒矩阵相位加权优化方法及巴特勒矩阵 Download PDF

Info

Publication number
CN112600592A
CN112600592A CN202011352213.2A CN202011352213A CN112600592A CN 112600592 A CN112600592 A CN 112600592A CN 202011352213 A CN202011352213 A CN 202011352213A CN 112600592 A CN112600592 A CN 112600592A
Authority
CN
China
Prior art keywords
input signal
butler matrix
phase
phase difference
vector synthesis
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.)
Granted
Application number
CN202011352213.2A
Other languages
English (en)
Other versions
CN112600592B (zh
Inventor
包晓军
刘远曦
李琳
王永刚
刘会涛
王育才
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Guangdong Narui Radar Technology Co ltd
Original Assignee
Guangdong Narui Radar Technology Co ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Guangdong Narui Radar Technology Co ltd filed Critical Guangdong Narui Radar Technology Co ltd
Priority to CN202011352213.2A priority Critical patent/CN112600592B/zh
Publication of CN112600592A publication Critical patent/CN112600592A/zh
Priority to PCT/CN2021/095624 priority patent/WO2022110699A1/zh
Application granted granted Critical
Publication of CN112600592B publication Critical patent/CN112600592B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

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/0408—Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas using two or more spaced independent antennas using two or more beams, i.e. beam diversity
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10—Complex mathematical operations
    • G06F17/16—Matrix or vector computation, e.g. matrix-matrix or matrix-vector multiplication, matrix factorization
    • H—ELECTRICITY
    • H01—ELECTRIC ELEMENTS
    • H01Q—ANTENNAS, i.e. RADIO AERIALS
    • H01Q3/00—Arrangements for changing or varying the orientation or the shape of the directional pattern of the waves radiated from an antenna or antenna system
    • H01Q3/26—Arrangements for changing or varying the orientation or the shape of the directional pattern of the waves radiated from an antenna or antenna system varying the relative phase or relative amplitude of energisation between two or more active radiating elements; varying the distribution of energy across a radiating aperture
    • H01Q3/30—Arrangements for changing or varying the orientation or the shape of the directional pattern of the waves radiated from an antenna or antenna system varying the relative phase or relative amplitude of energisation between two or more active radiating elements; varying the distribution of energy across a radiating aperture varying the relative phase between the radiating elements of an array
    • H01Q3/34—Arrangements for changing or varying the orientation or the shape of the directional pattern of the waves radiated from an antenna or antenna system varying the relative phase or relative amplitude of energisation between two or more active radiating elements; varying the distribution of energy across a radiating aperture varying the relative phase between the radiating elements of an array by electrical means
    • H01Q3/36—Arrangements for changing or varying the orientation or the shape of the directional pattern of the waves radiated from an antenna or antenna system varying the relative phase or relative amplitude of energisation between two or more active radiating elements; varying the distribution of energy across a radiating aperture varying the relative phase between the radiating elements of an array by electrical means with variable phase-shifters
    • 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

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Mathematical Physics (AREA)
  • Pure & Applied Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Mathematical Optimization (AREA)
  • Signal Processing (AREA)
  • Computational Mathematics (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Data Mining & Analysis (AREA)
  • Mathematical Analysis (AREA)
  • Software Systems (AREA)
  • General Engineering & Computer Science (AREA)
  • Databases & Information Systems (AREA)
  • Algebra (AREA)
  • Computing Systems (AREA)
  • Power Engineering (AREA)
  • Variable-Direction Aerials And Aerial Arrays (AREA)

Abstract

本发明公开了一种巴特勒矩阵相位加权优化方法及可以实现该方法的巴特勒矩阵，该方法基于n×n阶巴特勒矩阵，其中n＝2N，N为正整数，包括：在每个输入端口处分别设置一个可调移相器；n个输入信号分别经过对应的输入端口和对应的可调移相器后，进行N级矢量合成，并从n个输出端口处输出n个输出波束；获取对于任意一个输出波束，在每一级矢量合成中进行矢量合成的两个输入信号之间的相位差理论值和相位差实际值，进而获得两个输入信号之间的相位差误差值；根据相位差误差值，分别计算每个可调移相器所需的调节量并进行调节。根据本发明的巴特勒矩阵相位加权优化方法，能够降低巴特勒矩阵的相位偏差，提高巴特勒矩阵波束合成性能。

Description

巴特勒矩阵相位加权优化方法及巴特勒矩阵
技术领域
本发明涉及微波通信技术领域，尤其是涉及一种巴特勒矩阵相位加权优化方法及可以实现该方法的巴特勒矩阵。
背景技术
巴特勒矩阵是一种等幅相位加权网络，用于多波束合成，通过3dB正交分支线耦合器(或者180°环形耦合器)和若干个固定可调移相器实现不同的相位加权，形成多个不同指向的波束。采用巴特勒矩阵方式实现多波束，每一个波束均利用了整个阵面的所有辐射单元，因此对于相控阵波束合成，是一种无损的多波束合成方法。理想的巴特勒矩阵形成的多个波束之间具有正交性，每个波束的峰值和其他波束零值重合，波束之间相互干扰最小。
但是现有的巴特勒矩阵，在加工、制造过程中的各种误差很容易带来幅度和相位的偏差，进而影响到相控阵关键技术指标，如波束指向、最大旁瓣、积分旁瓣、波束隔离度、抗干扰能力等。这也就使得巴特勒矩阵目前主要还是应用在规模较小，而且对波束合成精确度要求不高的领域，如5G通信领域，而巴特勒矩阵在微波频段、大规模有源相控阵雷达领域的应用则还是非常少见的，因此需要提升巴特勒矩阵的波束合成精度。
发明内容
本发明旨在至少解决现有技术中存在的技术问题之一。为此，本发明提出了一种巴特勒矩阵相位加权优化方法，能够降低巴特勒矩阵的相位偏差，提高巴特勒矩阵波束合成性能。
本发明还提出了一种能够实现上述巴特勒矩阵相位加权优化方法的巴特勒矩阵。
根据本发明第一方面实施例的巴特勒矩阵相位加权优化方法，基于n×n阶巴特勒矩阵，所述n×n阶巴特勒矩阵具有n个输入端口及n个输出端口，其中n＝2N，N为正整数，N表示n个输入信号进行矢量合成的级数，其特征在于，包括：在每个所述输入端口处分别设置一个可调移相器，每个所述可调移相器的输入端分别与对应的所述输入端口电性连接；n个输入信号分别经过对应的所述输入端口和对应的所述可调移相器后，进行N级矢量合成，并从所述n个输出端口处输出n个输出波束；获取对于任意一个所述输出波束，在每一级矢量合成中进行矢量合成的两个输入信号之间的相位差理论值和相位差实际值，进而获得所述两个输入信号之间的相位差误差值；根据所述相位差误差值，计算每个所述可调移相器所需的调节量并进行调节。
根据本发明实施例的巴特勒矩阵相位加权优化方法，至少具有如下有益效果：通过算法计算对于每一个输出波束，在每一级矢量合成中，进行矢量合成的两个输入信号的相位差理论值，再通过实际测试获取这两个输入信号的相位差实际值，进而获得两个输入信号之间的相位差误差值，再计算每个可调移相器所需的调节量并进行调节，从而修正巴特勒矩阵的相位误差。根据本发明实施例的巴特勒矩阵相位加权优化方法，能够获得成品巴特勒矩阵各矢量合成点的理论值与实际值之间的偏差，再通过可调移相器统一进行修正，能够有效修正巴特勒矩阵在制造、加工环节产生的各种相位偏差，提高了巴特勒矩阵波束合成性能，进而增强巴特勒矩阵在大规模阵列波束合成领域的适用性。
根据本发明的一些实施例，n个输入信号进行N级矢量合成的方法为：在第m级矢量合成中，第i个输入信号与第i+n/2m个输入信号进行矢量合成，且在每一级矢量合成中，每个输入信号只进行一次矢量合成，其中m为不大于N的正整数，i为小于n的正整数。
根据本发明的一些实施例，在每一级矢量合成中，分别设置有n/2个耦合器，在第m级矢量合成中，第i个输入信号与第i+n/2m个输入信号通过对应的所述耦合器进行矢量合成。
根据本发明的一些实施例，所述耦合器为3dB耦合器或者环形耦合器。
根据本发明的一些实施例，所述n个输出波束呈中心对称，两侧各有n/2个输出波束，对于第一侧的第k个输出波束，在第m级矢量合成中，第i个输入信号与第i+n/2m个输入信号之间的相位差理论值为δ＝(n/2m)×(2k-1)π/n，对于第二侧的第k个输出波束，在第m级矢量合成中，第i个输入信号与第i+n/2m个输入信号之间的相位差理论值为δ＝(-1)×(n/2m)×(2k-1)π/n，其中k为不大于n/2的正整数。
根据本发明的一些实施例，将第j个输出端口的输出波束记为波束j，计算波束j对应于所述第一侧的第几个输出波束或者所述第二侧的第几个波束，其中j为不大于n的正整数。
根据本发明的一些实施例，对于波束j，在第m级矢量合成中，将第i个输入信号与第i+n/2m个输入信号之间的相位差理论值记为Amji，将第i个输入信号与第i+n/2m个输入信号之间的相位差实际值记为Bmji，则第i个输入信号与第i+n/2m个输入信号之间的相位差误差值为Cmji＝Amji-Bmji。
根据本发明的一些实施例，计算每个所述可调移相器所需的调节量的方法为：
Figure BDA0002801637090000031
Figure BDA0002801637090000032
则第i个可调移相器所需的调节量为Ci。
根据本发明第二方面实施例的巴特勒矩阵，可以实现如本发明上述第一方面实施例所述的巴特勒矩阵相位加权优化方法。
根据本发明实施例的巴特勒矩阵，至少具有如下有益效果：通过采用本发明上述第一方面实施例所述的巴特勒矩阵相位加权优化方法，能够有效修正巴特勒矩阵在制造、加工环节产生的各种相位偏差，提高了巴特勒矩阵波束合成性能，进而增强巴特勒矩阵在大规模阵列波束合成领域的适用性。
本发明的附加方面和优点将在下面的描述中部分给出，部分将从下面的描述中变得明显，或通过本发明的实践了解到。
附图说明
本发明的上述和/或附加的方面和优点从结合下面附图对实施例的描述中将变得明显和容易理解，其中：
图1为本发明实施例的巴特勒矩阵相位加权优化方法的流程示意图；
图2为现有技术中的8×8阶巴特勒矩阵的结构示意图；
图3为本发明实施例的8×8阶巴特勒矩阵的结构示意图。
具体实施方式
本部分将详细描述本发明的具体实施例，本发明之较佳实施例在附图中示出，附图的作用在于用图形补充说明书文字部分的描述，使人能够直观地、形象地理解本发明的每个技术特征和整体技术方案，但其不能理解为对本发明保护范围的限制。
在本发明的描述中，需要理解的是，涉及到方位描述，例如上、下、前、后、左、右等指示的方位或位置关系为基于附图所示的方位或位置关系，仅是为了便于描述本发明和简化描述，而不是指示或暗示所指的装置或元件必须具有特定的方位、以特定的方位构造和操作，因此不能理解为对本发明的限制。
在本发明的描述中，若干的含义是一个或者多个，多个的含义是两个以上，大于、小于、超过等理解为不包括本数，以上、以下、以内等理解为包括本数。如果有描述到第一、第二只是用于区分技术特征为目的，而不能理解为指示或暗示相对重要性或者隐含指明所指示的技术特征的数量或者隐含指明所指示的技术特征的先后关系。
本发明的描述中，除非另有明确的限定，设置、安装、连接等词语应做广义理解，所属技术领域技术人员可以结合技术方案的具体内容合理确定上述词语在本发明中的具体含义。
如图1所示，根据本发明第一方面实施例的巴特勒矩阵相位加权优化方法，基于n×n阶巴特勒矩阵，n×n阶巴特勒矩阵具有n个输入端口及n个输出端口，其中n＝2N，N为正整数，N表示n个输入信号进行矢量合成的级数，包括如下四个步骤：
S100：在每个输入端口处分别设置一个可调移相器，每个可调移相器的输入端分别与对应的输入端口电性连接；
S200：n个输入信号分别经过对应的输入端口和对应的可调移相器后，进行N级矢量合成，并从n个输出端口处输出n个输出波束；
S300：获取对于任意一个输出波束，在每一级矢量合成中进行矢量合成的两个输入信号之间的相位差理论值和相位差实际值，进而获得两个输入信号之间的相位差误差值；
S400：根据相位差误差值，分别计算每个可调移相器所需的调节量并进行调节。
根据本发明实施例的巴特勒矩阵相位加权优化方法，能够获取成品巴特勒矩阵的整体相位差误差值，再通过可调移相器进行修正，有效修正巴特勒矩阵在制造、加工环节中产生的各种相位偏差，提高巴特勒矩阵波束合成性能，进而增强巴特勒矩阵在大规模阵列波束合成领域的适用性。
下面以8×8阶巴特勒矩阵为例具体说明本发明实施例的巴特勒矩阵相位加权优化方法，可以理解的是，下述描述仅是示例性说明，而不是对发明的具体限制。
如图2和图3所示，对于8×8阶巴特勒矩阵，n＝8，N＝3，8个输入端口分别为P1至P8，8个输出端口分别为Q1至Q8，从8个输出端口输出的输出波束分别为波束1至波束8。在8个输入端口处设置8个可调移相器，分别为PH1至PH8，其中，图2示出的PH9至PH16为传统8×8阶巴特勒矩阵的固定式移相器，H1至H12为传统8×8阶巴特勒矩阵的矢量合成节点。8个输入信号经过相应的输入端口和相应的可调移相器后，进行3级矢量合成，并从8个输出端口处输出8个输出波束，随后获取对于任意一个输出波束，在每一级矢量合成中进行矢量合成的两个输入信号之间的相位差理论值和相位差实际值，进而获得两个输入信号之间的相位差误差值，进而计算每个可调移相器所需的调节量并进行调节。
n个输入信号进行N级矢量合成的方法为：在第m级矢量合成中，第i个输入信号与第i+n/2m个输入信号进行矢量合成，且在每一级矢量合成中，每个输入信号只进行一次矢量合成，其中m为不大于N的正整数，i为小于n的正整数。如图3所示，仍以8×8阶巴特勒矩阵为例，在第一级矢量合成中，m＝1，第1个输入信号与第5个输入信号进行矢量合成、第2个输入信号与第6个输入信号进行矢量合成，第3个输入信号与第7个输入信号进行矢量合成，第4个输入信号与第8个输入信号进行矢量合成；在第二级矢量合成中，m＝2，第1个输入信号与第3个输入信号进行矢量合成，第2个输入信号与第4个输入信号进行矢量合成，第5个输入信号与第7个输入信号进行矢量合成，第6个输入信号与第8个输入信号进行矢量合成；在第三级矢量合成中，m＝3，第1个输入信号与第2个输入信号进行矢量合成，第3个输入信号与第4个输入信号进行矢量合成，第5个输入信号与第6个输入信号进行矢量合成，第7个输入信号与第8个输入信号进行矢量合成。上述描述的方法同样适用于4×4阶巴特勒矩阵、16×16阶巴特勒矩阵等，而不限于此。
如图2和图3所示，在每一级矢量合成中，分别设置有n/2个耦合器，在第m级矢量合成中，第i个输入信号与第i+n/2m个输入信号通过对应的耦合器进行矢量合成。仍以8×8阶巴特勒矩阵为例，每一级矢量合成中，分别设置有4个耦合器，在第一级矢量合成中，4个耦合器为H1至H4，在第二级矢量合成中，4个耦合器为H5至H8，在第三级矢量合成中，4个耦合器为H9至H12。可以理解的是，对于巴特勒矩阵，两个输入信号通过耦合器进行矢量合成的原理和过程，为本领域技术人员所应该知道或者可以知道的公知常识，因而在此并未对矢量合成的原理和过程进行详细地阐述。耦合器可以采用3dB耦合器或者环形耦合器。
对于n×n阶巴特勒矩阵，n个输出波束会呈中心对称，两侧各有n/2个输出波束，对于第一侧的第k个输出波束，在第m级矢量合成中，第i个输入信号与第i+n/2m个输入信号之间的相位差理论值为δ＝(n/2m)×(2k-1)π/n，对于第二侧的第k个输出波束，在第m级矢量合成中，第i个输入信号与第i+n/2m个输入信号之间的相位差理论值为δ＝(-1)×(n/2m)×(2k-1)π/n，其中k为不大于n/2的正整数。以8×8阶巴特勒矩阵为例，8个输入信号经过巴特勒矩阵后会生成8个呈中心对称的输出波束，两侧各有4个输出波束，对于第一侧的第1个波束，在第一级矢量合成中，第1个输入信号与第5个输入信号的相位差理论值为π/2，第2个输入信号与第6个输入信号的相位差理论值为π/2，第3个输入信号与第7个输入信号的相位差理论值为π/2，第4个输入信号与第8个输入信号的相位差理论值为π/2；在第二级矢量合成中，第1个输入信号与第3个输入信号的相位差理论值为π/4，第2个输入信号与第4个输入信号的相位差理论值为π/4，第5个输入信号与第7个输入信号的相位差理论值为π/4，第6个输入信号与第8个输入信号的相位差理论值为π/4；在第三级矢量合成中，第1个输入信号与第2个输入信号的相位差理论值为π/8，第3个输入信号与第4个输入信号的相位差理论值为π/8，第5个输入信号与第6个输入信号的相位差理论值为π/8，第7个输入信号与第8个输入信号的相位差理论值为π/8。依次类推，可以获取对于任意输出波束，在每一级矢量合成中进行矢量合成的两个输入信号之间的相位差理论值，再通过对巴特勒矩阵进行实际测试，便可以获得两个输入信号之间的相位差实际值，进而获得两个输入信号之间的相位差误差值。
将第j个输出端口的输出波束记为波束j，计算波束j分别是对应于第一侧的第几个波束或者第二侧的第几个波束，其中j为不大于n的正整数。具体的计算方法如下：
(1)获取波束j对应的序号L，L＝j-1；
(2)将序号L转换为对应的二进制数字，并将该二进制数字逐位颠倒，再将颠倒后得到的数字转换为对应的十进制数字G；
(3)若G小于n/2，则波束j为第一侧的第(G+1)个波束，若G大于等于n/2，则波束j为第二侧的第(n-G)个波束。
参照下面的表1，以8×8阶巴特勒矩阵为例，波束1对应的序号L为0，0对应的二进制数字为000，逐位颠倒后的数字为000，对应的十进制数字G为0，则波束1对应于第一侧的第1个波束；波束2对应的序号L为1，1对应的二进制数字为001，逐位颠倒后的数字为100，对应的十进制数字G为4，则波束1对应于第二侧的第4个波束；依次类推，可以得到波束1至波束8分别对应于第一侧或者第二侧的第几个波束。可以理解的是，上述方法同样适用于4×4阶巴特勒矩阵、16×16阶巴特勒矩阵等，而不限于此。
Figure BDA0002801637090000081
表1
对于波束j，在第m级矢量合成中，将第i个输入信号与第i+n/2m个输入信号之间的相位差理论值记为Amji，将第i个输入信号与第i+n/2m个输入信号之间的相位差实际值记为Bmji，则第i个输入信号与第i+n/2m个输入信号之间的相位差误差值为Cmji＝Amji-Bmji。仍以8×8阶巴特勒矩阵为例，对于波束1，在第一级矢量合成中，第1个输入信号与第5个输入信号之间的相位差理论值为A111，相位差实际值为B111，相位差误差值C111＝A111-B111；第2个输入信号与第6个输入信号之间的相位差理论值为A112，相位差实际值为B112，相位差误差值C112＝A112-B112；第3个输入信号与第7个输入信号之间的相位差理论值为A113，相位差实际值为B113，相位差误差值C113＝A113-B113；第4个输入信号与第8个输入信号之间的相位差理论值为A114，相位差实际值为B114，相位差误差值C114＝A114-B114；在第二级矢量合成中，第1个输入信号与第3个输入信号之间的相位差理论值为A211，相位差实际值为B211，相位差误差值C211＝A211-B211；第2个输入信号与第4个输入信号之间的相位差理论值为A212，相位差实际值为B212，相位差误差值C212＝A212-B212；第5个输入信号与第7个输入信号之间的相位差理论值为A215，相位差实际值为B215，相位差误差值C215＝A215-B215；第6个输入信号与第8个输入信号之间的相位差理论值为A216，相位差实际值为B216，相位差误差值C216＝A216-B216；在第三级矢量合成中，第1个输入信号与第2个输入信号之间的相位差理论值为A311，相位差实际值为B311，相位差误差值C311＝A311-B211；第3个输入信号与第4个输入信号之间的相位差理论值为A313，相位差实际值为B313，相位差误差值C313＝A313-B213；第5个输入信号与第6个输入信号之间的相位差理论值为A315，相位差实际值为B315，相位差误差值C315＝A315-B315；第7个输入信号与第8个输入信号之间的相位差理论值为A317，相位差实际值为B317，相位差误差值C317＝A317-B317。同理，可以计算出C121、C122、C123、C124、C221、C222、C225、C226、C321、C323、C325和C327等
计算每个可调移相器所需的调节量的方法为：
Figure BDA0002801637090000091
Figure BDA0002801637090000101
则第i个可调移相器的所需的调节量为Ci。以8×8阶巴特勒矩阵为例：
C11＝1/8×(C111+C121+C131+C141+C151+C161+C171+C181)，
C12＝1/8×(C112+C122+C132+C142+C152+C162+C172+C182)，
C13＝1/8×(C113+C123+C133+C143+C153+C163+C173+C183)，
C14＝1/8×(C114+C124+C134+C144+C154+C164+C174+C184)，
C21＝1/8×(C211+C221+C231+C241+C251+C261+C271+C281)，
C22＝1/8×(C212+C222+C232+C242+C252+C262+C272+C282)，
C25＝1/8×(C215+C225+C235+C245+C255+C265+C275+C285)，
C26＝1/8×(C216+C226+C236+C246+C256+C266+C276+C286)，
C31＝1/8×(C311+C321+C331+C341+C351+C361+C371+C381)，
C33＝1/8×(C313+C323+C333+C343+C353+C363+C373+C383)，
C35＝1/8×(C315+C325+C335+C345+C355+C365+C375+C385)，
C37＝1/8×(C317+C327+C337+C347+C357+C367+C377+C387)，
根据上述式子，可以计算得到：
C1＝C11+C21+C31，
C2＝C12+C22，
C3＝C13+C33，
C4＝C14，
C5＝C25+C35，
C6＝C26，
C7＝C37，
C8＝0，
由上述式子可知，可调移相器PH1的调节量为C1，可调移相器PH1的调节量为C1，可调移相器PH2的调节量为C2，可调移相器PH3的调节量为C3，可调移相器PH4的调节量为C4，可调移相器PH5的调节量为C5，可调移相器PH6的调节量为C6，可调移相器PH7的调节量为C7，可调移相器PH8的调节量为C8。根据每个可调移相器的调节量进行调节，便能修正巴特勒矩阵的整体相位差误差值，最大限度降低巴特勒矩阵在制造、加工环节产生的各种相位偏差，提高了巴特勒矩阵波束合成性能，进而增强巴特勒矩阵在大规模阵列波束合成领域的适用性。
参考图1至图3，下面以8×8阶巴特勒矩阵为例，以一个具体的实施例详细描述根据本发明实施例的巴特勒矩阵相位加权优化方法。值得理解的是，下述描述仅是示例性说明，而不是对发明的具体限制。
对于8×8阶巴特勒矩阵，n＝8，N＝3，在8个输出端口P1至P8处设置8个可调移相器PH1至PH8，可调移相器可以采用市面上常见的模拟可调移相器或数字可调移相器。8个输入信号经过对应的输入端口和对应的可调移相器后，进行3级矢量合成，在第一级矢量合成中，m＝1，第1个输入信号与第5个输入信号进行矢量合成，对于8个输出波束，相位差误差值分别为C111至C181；第2个输入信号与第6个输入信号进行矢量合成，对于8个输出波束，相位差误差值分别为C112至C182；第3个输入信号与第7个输入信号进行矢量合成，对于8个输出波束，相位差误差值分别为C113至C183；第4个输入信号与第8个输入信号进行矢量合成，对于8个输出波束，相位差误差值分别为C114至C184；在第二级矢量合成中，m＝2，第1个输入信号与第3个输入信号进行矢量合成，对于8个输出波束，相位差误差值分别为C211至C281；第2个输入信号与第4个输入信号进行矢量合成，对于8个输出波束，相位差误差值分别为C212至C282；第5个输入信号与第7个输入信号进行矢量合成，对于8个输出波束，相位差误差值分别为C215至C285；第6个输入信号与第8个输入信号进行矢量合成，对于8个输出波束，相位差误差值分别为C216至C286；在第三级矢量合成中，m＝3，第1个输入信号与第2个输入信号进行矢量合成，对于8个输出波束，相位差误差值分别为C311至C381；第3个输入信号与第4个输入信号进行矢量合成，对于8个输出波束，相位差误差值分别为C313至C383；第5个输入信号与第6个输入信号进行矢量合成，对于8个输出波束，相位差误差值分别为C315至C385；第7个输入信号与第8个输入信号进行矢量合成，对于8个输出波束，相位差误差值分别为C317至C387，则由
Figure BDA0002801637090000121
可以得出：
C11＝1/8×(C111+C121+C131+C141+C151+C161+C171+C181)，
C12＝1/8×(C112+C122+C132+C142+C152+C162+C172+C182)，
C13＝1/8×(C113+C123+C133+C143+C153+C163+C173+C183)，
C14＝1/8×(C114+C124+C134+C144+C154+C164+C174+C184)，
C21＝1/8×(C211+C221+C231+C241+C251+C261+C271+C281)，
C22＝1/8×(C212+C222+C232+C242+C252+C262+C272+C282)，
C25＝1/8×(C215+C225+C235+C245+C255+C265+C275+C285)，
C26＝1/8×(C216+C226+C236+C246+C256+C266+C276+C286)，
C31＝1/8×(C311+C321+C331+C341+C351+C361+C371+C381)，
C33＝1/8×(C313+C323+C333+C343+C353+C363+C373+C383)，
C35＝1/8×(C315+C325+C335+C345+C355+C365+C375+C385)，
C37＝1/8×(C317+C327+C337+C347+C357+C367+C377+C387)，
再根据
Figure BDA0002801637090000122
可以得到：
C1＝C11+C21+C31，
C2＝C12+C22，
C3＝C13+C33，
C4＝C14，
C5＝C25+C35，
C6＝C26，
C7＝C37，
C8＝0，
综上可以得到，可调移相器PH1的调节量为C1，可调移相器PH1的调节量为C1，可调移相器PH2的调节量为C2，可调移相器PH3的调节量为C3，可调移相器PH4的调节量为C4，可调移相器PH5的调节量为C5，可调移相器PH6的调节量为C6，可调移相器PH7的调节量为C7，可调移相器PH8的调节量为C8。根据每个可调移相器的调节量进行调节，便能修正巴特勒矩阵的整体相位差误差值，修正巴特勒矩阵在制造、加工环节产生的各种相位偏差，提高了巴特勒矩阵波束合成性能，进而增强巴特勒矩阵在大规模阵列波束合成领域的适用性。
根据本发明第二方面实施例的巴特勒矩阵，通过采用上述巴特勒矩阵相位加权优化方法，提高了波束合成性能，能够适用于大规模阵列波束合成领域。
在本说明书的描述中，参考术语“一个实施例”、“进一步实施例”、“一些具体实施例”或“一些示例”等的描述意指结合该实施例或示例描述的具体特征、结构或者特点包含于本发明的至少一个实施例或示例中。在本说明书中，对上述术语的示意性表述不一定指的是相同的实施例或示例。而且，描述的具体特征、结构、材料或者特点可以在任何的一个或多个实施例或示例中以合适的方式结合。
尽管已经示出和描述了本发明的实施例，本领域的普通技术人员可以理解：在不脱离本发明的原理和宗旨的情况下可以对这些实施例进行多种变化、修改、替换和变型，本发明的范围由权利要求及其等同物限定。

Claims (9)

1.一种巴特勒矩阵相位加权优化方法，基于n×n阶巴特勒矩阵，所述n×n阶巴特勒矩阵具有n个输入端口及n个输出端口，其中n＝2N，N为正整数，N表示n个输入信号进行矢量合成的级数，其特征在于，包括：
在每个所述输入端口处分别设置一个可调移相器，每个所述可调移相器的输入端分别与对应的所述输入端口电性连接；
n个输入信号分别经过对应的所述输入端口和对应的所述可调移相器后，进行N级矢量合成，并从所述n个输出端口处输出n个输出波束；
获取对于任意一个所述输出波束，在每一级矢量合成中进行矢量合成的两个输入信号之间的相位差理论值和相位差实际值，进而获得所述两个输入信号之间的相位差误差值；
根据所述相位差误差值，分别计算每个所述可调移相器所需的调节量并进行调节。
2.根据权利要求1所述的巴特勒矩阵相位加权优化方法，其特征在于，n个输入信号进行N级矢量合成的方法为：在第m级矢量合成中，第i个输入信号与第i+n/2m个输入信号进行矢量合成，且在每一级矢量合成中，每个输入信号只进行一次矢量合成，其中m为不大于N的正整数，i为小于n的正整数。
3.根据权利要求2所述的巴特勒矩阵相位加权优化方法，其特征在于，在每一级矢量合成中，分别设置有n/2个耦合器，在第m级矢量合成中，第i个输入信号与第i+n/2m个输入信号通过对应的所述耦合器进行矢量合成。
4.根据权利要求3所述的巴特勒矩阵相位加权优化方法，其特征在于，所述耦合器为3dB耦合器或者环形耦合器。
5.根据权利要求2所述的巴特勒矩阵相位加权优化方法，其特征在于，所述n个输出波束呈中心对称，两侧各有n/2个输出波束，对于第一侧的第k个输出波束，在第m级矢量合成中，第i个输入信号与第i+n/2m个输入信号之间的相位差理论值为δ＝(n/2m)×(2k-1)π/n，对于第二侧的第k个输出波束，在第m级矢量合成中，第i个输入信号与第i+n/2m个输入信号之间的相位差理论值为δ＝(-1)×(n/2m)×(2k-1)π/n，其中k为不大于n/2的正整数。
6.根据权利要求5所述的巴特勒矩阵相位加权优化方法，其特征在于，将第j个输出端口的输出波束记为波束j，计算波束j对应于所述第一侧的第几个输出波束或者所述第二侧的第几个波束，其中j为不大于n的正整数。
7.根据权利要求6所述的巴特勒矩阵相位加权优化方法，其特征在于，对于波束j，在第m级矢量合成中，将第i个输入信号与第i+n/2m个输入信号之间的相位差理论值记为Amji，将第i个输入信号与第i+n/2m个输入信号之间的相位差实际值记为Bmji，则第i个输入信号与第i+n/2m个输入信号之间的相位差误差值为Cmji＝Amji-Bmji。
8.根据权利要求7所述的巴特勒矩阵相位加权优化方法，其特征在于，计算每个所述可调移相器所需的调节量的方法为：
Figure FDA0002801637080000021
Figure FDA0002801637080000022
则第i个可调移相器所需的调节量为Ci。
9.一种巴特勒矩阵，其特征在于，可以实现如权利要求1至8任一项所述的巴特勒矩阵相位加权优化方法。
CN202011352213.2A 2020-11-27 2020-11-27 巴特勒矩阵相位加权优化方法及巴特勒矩阵 Active CN112600592B (zh)

Priority Applications (2)

Application Number Priority Date Filing Date Title
CN202011352213.2A CN112600592B (zh) 2020-11-27 2020-11-27 巴特勒矩阵相位加权优化方法及巴特勒矩阵
PCT/CN2021/095624 WO2022110699A1 (zh) 2020-11-27 2021-05-25 巴特勒矩阵相位加权优化方法及巴特勒矩阵

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202011352213.2A CN112600592B (zh) 2020-11-27 2020-11-27 巴特勒矩阵相位加权优化方法及巴特勒矩阵

Publications (2)

Publication Number Publication Date
CN112600592A true CN112600592A (zh) 2021-04-02
CN112600592B CN112600592B (zh) 2021-10-08

Family

ID=75184303

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202011352213.2A Active CN112600592B (zh) 2020-11-27 2020-11-27 巴特勒矩阵相位加权优化方法及巴特勒矩阵

Country Status (2)

Country Link
CN (1) CN112600592B (zh)
WO (1) WO2022110699A1 (zh)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113659353A (zh) * 2021-08-02 2021-11-16 电子科技大学 一种小型化输出相差360°连续可调的巴特勒矩阵
WO2022110699A1 (zh) * 2020-11-27 2022-06-02 广东纳睿雷达科技股份有限公司 巴特勒矩阵相位加权优化方法及巴特勒矩阵

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20230188187A1 (en) * 2021-12-14 2023-06-15 Qualcomm Incorporated Beamforming circuit using configurable phase shifters for a butler matrix architecture

Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4973971A (en) * 1989-12-18 1990-11-27 Allied-Signal Inc. Broadband circular phased array antenna
CN203071234U (zh) * 2012-11-28 2013-07-17 寰波科技股份有限公司 巴特勒矩阵
CN103594801A (zh) * 2013-11-21 2014-02-19 天津中兴智联科技有限公司 一种巴特勒矩阵结构
CN104934675A (zh) * 2015-06-23 2015-09-23 西安空间无线电技术研究所 一种用于星载轻型sar的大功率合成器和合成方法
CN106025475A (zh) * 2016-06-13 2016-10-12 北京航空航天大学 一种基于Butler阵的新型6x6多波束网络
CN109613512A (zh) * 2018-12-06 2019-04-12 上海交通大学 基于诺伦矩阵的n×m集成多波束激光雷达发射系统
CN109888507A (zh) * 2018-12-22 2019-06-14 中国电波传播研究所（中国电子科技集团公司第二十二研究所） 一种紧凑型16×16 Butler矩阵多波束馈电网络
CN110034415A (zh) * 2019-03-07 2019-07-19 中山大学 一种具有宽带特性的诺兰矩阵及其制造方法
CN111812627A (zh) * 2020-07-20 2020-10-23 电子科技大学 一种基于fpga的可控多波束形成方法

Family Cites Families (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
FR2750258B1 (fr) * 1996-06-24 1998-08-21 Europ Agence Spatiale Systeme de conformation de faisceau zonal reconfigurable pour une antenne embarquee sur un satellite en orbite et procede d'optimisation de la reconfiguration
KR101809383B1 (ko) * 2016-12-19 2017-12-14 전자부품연구원 하이브리드 버틀러 매트릭스 및 이를 이용한 통신장치
IT201800006163A1 (it) * 2018-06-08 2019-12-08 Rete di ripartizione a radiofrequenza riconfigurabile
CN112600592B (zh) * 2020-11-27 2021-10-08 广东纳睿雷达科技股份有限公司 巴特勒矩阵相位加权优化方法及巴特勒矩阵

Patent Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4973971A (en) * 1989-12-18 1990-11-27 Allied-Signal Inc. Broadband circular phased array antenna
CN203071234U (zh) * 2012-11-28 2013-07-17 寰波科技股份有限公司 巴特勒矩阵
CN103594801A (zh) * 2013-11-21 2014-02-19 天津中兴智联科技有限公司 一种巴特勒矩阵结构
CN104934675A (zh) * 2015-06-23 2015-09-23 西安空间无线电技术研究所 一种用于星载轻型sar的大功率合成器和合成方法
CN106025475A (zh) * 2016-06-13 2016-10-12 北京航空航天大学 一种基于Butler阵的新型6x6多波束网络
CN109613512A (zh) * 2018-12-06 2019-04-12 上海交通大学 基于诺伦矩阵的n×m集成多波束激光雷达发射系统
CN109888507A (zh) * 2018-12-22 2019-06-14 中国电波传播研究所（中国电子科技集团公司第二十二研究所） 一种紧凑型16×16 Butler矩阵多波束馈电网络
CN110034415A (zh) * 2019-03-07 2019-07-19 中山大学 一种具有宽带特性的诺兰矩阵及其制造方法
CN111812627A (zh) * 2020-07-20 2020-10-23 电子科技大学 一种基于fpga的可控多波束形成方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
ERIC T.DER等: "Miniaturized 4 × 4 Butler Matrix and Tunable Phase Shifter Using Ridged Half-Mode Substrate Integrated Waveguide", 《IEEE TRANSACTIONS ON MICROWAVE THEORY AND TECHNIQUES》 *
PEDRAM YAZDANBAKHSH等: "Phase error corrections for Microstrip 8×8 Butler Matrix used in 7 Tesla MRI", 《2011 GERMAN MICROWAVE CONFERENCE》 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2022110699A1 (zh) * 2020-11-27 2022-06-02 广东纳睿雷达科技股份有限公司 巴特勒矩阵相位加权优化方法及巴特勒矩阵
CN113659353A (zh) * 2021-08-02 2021-11-16 电子科技大学 一种小型化输出相差360°连续可调的巴特勒矩阵

Also Published As

Publication number Publication date
WO2022110699A1 (zh) 2022-06-02
CN112600592B (zh) 2021-10-08

Similar Documents

Publication Publication Date Title
CN112600592B (zh) 巴特勒矩阵相位加权优化方法及巴特勒矩阵
CN106656099B (zh) 数字移相器
JP2014220812A (ja) マルチポート増幅器を較正する方法、このような方法の実装を可能にするマルチポート増幅器、およびこのような増幅器を備える衛星
WO2019214257A1 (zh) 波束指向可调天线的方向图的计算方法及装置
CN113296049B (zh) 互质阵列脉冲环境下非圆信号的共轭增广doa估计方法
CN113343588A (zh) 一种多约束毫米波车载mimo雷达天线阵列设计方法与装置
CN109752687B (zh) 一种降低阵元互耦效应的改进型互质阵列设置方法
KR20030009216A (ko) 시스토릭 어레이 장치
Yang et al. High-order cumulants based sparse array design via fractal geometries—Part I: Structures and DOFs
US5777579A (en) Low cost butler matrix modeformer circuit
CN107889542A (zh) 信号分配网络
US5101211A (en) Closed loop RF power amplifier output correction circuit
CN113783550B (zh) 一种用于k波段的高精度数控移相器及其移相方法
CN110516286A (zh) 基于稀疏重构算法的阵列幅相误差估计方法
CN110361697A (zh) 一种基于协方差矩阵混合重构的稳健波束形成方法
CN109871580A (zh) 一种基于简化实频方法的放大器匹配网络设计方法
CN112347681B (zh) 一种基于宏基函数阵列互耦特性预测的稳健波束形成方法
WO2015140533A1 (en) Isolation in a multi-port amplifier
CN110333504B (zh) 一种空时二维滤波的快速宽带波束形成方法
JP3113837B2 (ja) バトラーマトリクス形式のｎポートモードフォーマのバイアスエラー減少方法及び装置
CN112100837A (zh) 一种用于非圆信号波达方向估计的二级扩展变换嵌套阵设计方法
TWI775296B (zh) 使用環張量之深度神經網路加速方法及其系統
CN115378446B (zh) 基于神经网络的宽带数字预失真系统和方法
CN117491751A (zh) 一种数字有源相控阵的阵列方向图校准算法
CN117706467B (zh) 一种干涉仪布阵方法

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