CN113804985B - 一种基于混合屏蔽室的抗干扰天线方向图测量方法 - Google Patents
一种基于混合屏蔽室的抗干扰天线方向图测量方法 Download PDFInfo
- Publication number
- CN113804985B CN113804985B CN202111001299.9A CN202111001299A CN113804985B CN 113804985 B CN113804985 B CN 113804985B CN 202111001299 A CN202111001299 A CN 202111001299A CN 113804985 B CN113804985 B CN 113804985B
- Authority
- CN
- China
- Prior art keywords
- antenna
- calibration
- coefficient
- ref
- measured
- 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 62
- 238000012360 testing method Methods 0.000 claims abstract description 74
- 238000010586 diagram Methods 0.000 claims abstract description 56
- 239000011159 matrix material Substances 0.000 claims abstract description 42
- 230000005540 biological transmission Effects 0.000 claims abstract description 36
- 230000005855 radiation Effects 0.000 claims abstract description 21
- 238000005259 measurement Methods 0.000 claims description 11
- 230000005684 electric field Effects 0.000 claims description 9
- 239000000523 sample Substances 0.000 claims description 6
- 238000005070 sampling Methods 0.000 claims description 4
- 238000004364 calculation method Methods 0.000 claims description 3
- 230000002452 interceptive effect Effects 0.000 claims 1
- 230000000694 effects Effects 0.000 description 7
- 239000006096 absorbing agent Substances 0.000 description 2
- 239000011358 absorbing material Substances 0.000 description 2
- 238000004891 communication Methods 0.000 description 2
- 238000007796 conventional method Methods 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 239000000284 extract Substances 0.000 description 2
- 238000000691 measurement method Methods 0.000 description 2
- 238000012545 processing Methods 0.000 description 2
- 238000003756 stirring Methods 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 230000002950 deficient Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 230000001771 impaired effect Effects 0.000 description 1
- 238000013507 mapping Methods 0.000 description 1
- 238000011056 performance test Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 230000035945 sensitivity Effects 0.000 description 1
- 238000010561 standard procedure Methods 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 238000010998 test method Methods 0.000 description 1
- 238000012795 verification Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R29/00—Arrangements for measuring or indicating electric quantities not covered by groups G01R19/00 - G01R27/00
- G01R29/08—Measuring electromagnetic field characteristics
- G01R29/10—Radiation diagrams of antennas
- G01R29/105—Radiation diagrams of antennas using anechoic chambers; Chambers or open field sites used therefor
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R29/00—Arrangements for measuring or indicating electric quantities not covered by groups G01R19/00 - G01R27/00
- G01R29/08—Measuring electromagnetic field characteristics
- G01R29/0864—Measuring electromagnetic field characteristics characterised by constructional or functional features
- G01R29/0892—Details related to signal analysis or treatment; presenting results, e.g. displays; measuring specific signal features other than field strength, e.g. polarisation, field modes, phase, envelope, maximum value
Landscapes
- Physics & Mathematics (AREA)
- Electromagnetism (AREA)
- General Physics & Mathematics (AREA)
- Monitoring And Testing Of Transmission In General (AREA)
Abstract
一种基于混合屏蔽室的抗干扰天线方向图测量方法,选取两个天线,一个待测天线,另一个校准天线,在标准暗室环境下测量校准天线的辐射方向图并求得校准天线的传输系数;在混合屏蔽室环境下,先将待测天线放置在室内的转台上,控制待测天线进行旋转,得到待测天线旋转矩阵;再将校准天线放置在转台上进行旋转,得到校准天线的旋转矩阵;计算待测天线和校准天线的相关性系数矩阵,以及待测天线和校准天线在不同转台位置时的旋转方向图和球面波系数矩阵;利用球面波系数矩阵和相关性系数矩阵重构出待测天线的方向图,同时利用校准天线传输系数来补偿待测天线方向图,完成待测天线重构和补偿工作。本发明改善混合测试环境下的待测天线辐射方向图。
Description
技术领域
本发明属于微波测量技术领域,具体涉及一种基于混合屏蔽室的抗干扰天线方向图测量方法。
背景技术
混合屏蔽室通过在室内放置负载或吸收体从而形成静态、各向异性以及多路径的环境,并能模拟出特定的空间信道模型。而在无线通信测试中,每个新提出的应用模型均需要在定向信道的环境下进行测试,对此混合屏蔽室提供了必要的测试基础和标准化方法,并且混合屏蔽室因其自身的构建成本、可操作性以及***的灵活性使其在无线通信产业中具有广泛的应用前景和潜力。相关微波测量研究已经证明混合屏蔽室在创建样本空间信道环境方面的可行性。
在考虑混合屏蔽室测量环境及特性的前提下,通过借鉴暗室及混响室对天线效率、方向图以及辐射功率等参数测量的方法,可以看出混合屏蔽室具有在多径环境下完成天线相关参数及性能测试的应用前景。其中,天线方向图是判断并衡量待测天线的方向性、半功率波束宽度以及增益等参数的重要指标和参考,因此在混合屏蔽室中实现对待测天线方向图的测量是一项具有重要意义的工作。
在理想暗室环境下可以通过多探头球面近场测试***得到待测天线的球面近场数据,并利用球面波理论得到天线的辐射方向图;而在理想混响室环境下由于存在视线传输距离及多反射径干扰的影响,因此很难直接对待测天线方向图进行测量。目前,南京航空航天大学的徐千团队提出了利用混响室测量待测天线自相关系数反向重构出天线辐射方向图的方法,该方法从视线传输距离以及搅拌响应中提取待测天线的三维方向信息,其中传输天线不需要指向待测天线,利用混响室的丰富多径环境,通过旋转搅拌器或移动待测天线来提取混响室内未搅拌部分的信息。
但是传统微波暗室或混响室无法满足在特定信道环境下天线方向图的测试工作,并且在实际测试环境中,非理想暗室或混响室会因吸波材料、腔室尺寸以及测量器件等因素而无法满足理想测试环境的要求,从而造成测试结果上的误差。对此,在混合屏蔽室的测试环境下,可以通过调整室内吸波负载位置及数量的方式完成在特定空间信道模型下的测试任务需求,此外由于混和屏蔽室引入了多径效应的问题,因此需要通过借鉴重构及补偿的方法对混和屏蔽室的特定信道环境下的方向图进行恢复。
现有的补偿或重构的方法主要包括测试区域场补偿方法、反卷积方法、天线方向图校准方法以及相关系数重构方法,但上述方法都只针对非理想暗室或非理想混响室的其中一种测试场景具有较好的方向图补偿和重构效果。例如在非理想暗室测试环境中,当室内吸波材料性能较差或在某些角落未安装吸波结构且球面近场测试还存在球面反射径干扰的问题时,就会导致在天线辐射方向图的测量结果上出现由腔室内多反射径所造成的干扰,并损害方向图的主瓣、旁瓣以及方向性等特性参数,此时用传统的测试区域场补偿方法就无法较好地补偿待测天线的辐射特性;此外,在非理想混响室中,如果转台上的待测天线未能按照既定角度进行旋转或室内搅拌不均匀,则会引起球面散射信号不完整或缺损的问题,进而通过矢量网络分析仪所得到的S参数就无法准确计算出待测天线的相关性系数,因此用相关系数重构的方法只能得到含有反射径干扰的待测天线方向图,而无法反向重构出完整的天线远场辐射特性。
发明内容
为了克服以上技术问题,本发明提出一种基于混合屏蔽室的抗干扰天线方向图测量方法,能够降低测试场景中多反射径干扰并能重构和补偿待测天线球面波系数以及辐射特性,本方法将测试区域场补偿技术应用于混合屏蔽室当中,利用区域补偿方法和相关性系数反向重构矩阵解决多反射径干扰的影响,从而改善混合测试环境下的待测天线辐射方向图。
为了实现上述目的,本发明采用的技术方案是:
一种基于混合屏蔽室的抗干扰天线方向图测量方法,包括以下步骤;
步骤1.选取两个天线,将其中一个天线作为待测天线105,另一个天线作为校准天线106,在标准暗室环境下测量校准天线106的辐射方向图并求得校准天线106的传输系数;
步骤2.在混合屏蔽室环境下,先将待测天线105放置在室内的转台103上,控制待测天线105进行旋转,得到待测天线105的S参数矩阵;再将校准天线106放置在转台上进行旋转,得到校准天线106的S参数矩阵;
步骤3.计算待测天线105和校准天线106的相关性系数矩阵和校准天线106的自相关系数矩阵;
步骤4.利用所述球面波系数矩阵和相关性系数矩阵重构出待测天线的方向图,同时利用步骤1所述的校准天线106传输系数来补偿待测天线方向图,完成待测天线105重构和补偿工作。
所述步骤2中,转台103能够控制待测天线105和校准天线106沿X、Y和Z方向进行旋转,三个方向的旋转范围是0°-360°,角度间隔是1度;
根据测试区域场补偿理论,利用式(2)将混合屏蔽室中含反射径干扰的待测天线105远场分量表示为传输系数与信源相乘的形式:
其中v表示待测天线上的发射信源,Tmn表示待测天线的传输系数,与式(1)中的和等价,均表示球面波生成函数。同理,通过式(1)(2)的求解同样可以得到所述校准天线106的远场电场分量以及球面波系数amn和bmn;
将待测天线105在含多径反射影响下计算得到的球面波系数记为a'AUT,mn和b'AUT,mn,而待求的理论无干扰的待测天线球面波系数记为aAUT,mn和bAUT,mn,待求的理论无干扰的待测天线传输系数记为TAUT,mn;将所述已知的校准天线106传输系数记为TREF,mn,含多径反射影响的校准天线106球面波系数记为a'REF,mn和b'REF,mn,已知无干扰的信源信号记为v,含多径反射影响的信源信号记为v';
通过式(1)(2)可以看出,含多径反射影响的信源信号v'可以表示为式(3)的形式:
(a'REF,mn,b'REF,mn)=TREF,mnv' (3)
而待测天线105和校准天线106是在相同的混合屏蔽室环境下进行的测试,故两种天线测试结果下的信号v'一致,因此通过式(4)可以根据所述已知的校准天线106传输系数TREF,mn和计算得到的校准天线106球面波系数a'REF,mn和b'REF,mn推得待测天线的理论传输系数TAUT,mn:
TAUT,mn=(a'AUT,mn,b'AUT,mn)((TREF,mn)+·(a'REF,mn,b'REF,mn))+ (4)
其中+表示矩阵的广义逆矩阵。
所述步骤3中根据式5计算旋转待测天线105与校准天线106之间的相关性系数以及校准天线106自身的自相关系数:
其中amn1和bmn1表示含反射径影响的待测天线105的球面波系数,amn2和bmn2表示含反射径影响的校准天线106的球面波系数,T(α,β,γ)表示待测天线105的旋转矩阵,通过T(α,β,γ)与amn2和bmn2的矩阵相乘,得到校准天线106的旋转方向图,利用式(5)得到旋转校准天线106与待测天线105之间的相关系数,记为ρ1;同理,将式(5)中的待测天线105换成校准天线106,即amn1和bmn1用amn2和bmn2替换,得到校准天线106的旋转方向图以及校准天线106与其自身旋转方向图之间的自相关系数,记为ρ2;
为了在混响室中得到所述ρ1和ρ2,通过测量不同旋转位置处的待测天线105和校准天线106的S参数,并计算S参数之间的相关性,从而得到对应的相关系数ρ1和自相关系数ρ2,如式(6)所示:
其中,S21,k,REF和S21,k,AUT是发射天线到校准天线106与待测天线105之间的S参数,<S21,REF>和<S21,AUT>则对应上述两种参数的所有采样值的平均值,利用(6)式即可得到相关系数ρ1;若将(6)式中的S21,k,AUT和<S21,AUT>用S'21,k,REF和<S'21,REF>进行替换,其中S'21,k,REF和<S'21,REF>表示校准天线106在未旋转状态下的S参数和对应的平均值,则可得到自相关系数ρ2。
利用上述已知的ρ1和ρ2相关性系数以及旋转矩阵T(α,β,γ),通过式(7)所示的反向重构法,得到校准天线106和待测天线105的含反射径干扰的球面波系数:
其中Tα,Tβ和Tγ分别对应所述旋转矩阵T(α,β,γ)在X,Y和Z上的投影矩阵,ρs,α,ρs,β和ρs,γ分别表示相关性系数ρ1(下标s取1)或自相关系数ρ2(下标s取2)在X,Y和Z上的投影矩阵,利用广义逆矩阵的方法反向推导出待测天线105在含球面反射径干扰时的方向图,同理,将式(7)中的待测天线105和校准天线106互换,便可得到含反射径干扰的校准天线方向图。
所述步骤4中,对已求得的待测天线方向图进行补偿,在混合屏蔽室的条件下,待测天线105和校准天线106的测试环境保持一致,应用测试区域场补偿,即将式(3)(4)中的信号v'看作混合屏蔽室中因球面反射径干扰所产生的信源信号,进而通过式(8)的逆矩阵补偿方式得到去除球面反射径干扰的待测天线传输系数:
TAUT,mn=(amn1,bmn1)((TREF,mn)+·(amn2,bmn2))+ (8)
式(8)中利用了上述计算所得的含反射径干扰的待测天线球面波系数amn1和bmn1,以及含反射径干扰的校准天线球面波系数amn2和bmn2,并利用所述已知的校准天线传输系数TREF,mn补偿得到改善后的待测天线方向图。
本发明的有益效果:
本发明所述的球面波系数重构和补偿方法将测试区域场补偿技术引入到了混合屏蔽室测试当中,通过球面近场模式敏感度和校准天线传输系数测量,补偿待测天线在混合屏蔽室中因近场球面反射径干扰所造成的方向图测量损害。
本发明所述的球面波系数重构和补偿方法将混合测试环境下的多反射径干扰进行了划分,其中通过传统混响室测试方案得到的重构球面波系数解决了反射径干扰的影响,而测试区域场补偿方法则改善了球面近场反射对方向图的影响,相比于传统单一的测试区域补偿和重构算法,本发明所述方法更加有效地划分了混合环境下的干扰特性,并用针对性的方法有效地重构并补偿了待测天线辐射方向图。
本发明所述的球面波系数重构和补偿方法适用于混合屏蔽室的测试环境,相比于传统的测试区域场补偿方法所应用的非理想暗室及其他复杂环境,所述测试方法的应用测试环境在设计结构、成本以及操作性上都具有显著优势,整体测试过程只需通过旋转转台来控制待测天线和校准天线沿不同轴向进行旋转,从而得到待补偿的待测天线重构球面波系数,最后利用测试区域补偿算法即可得到改善后的待测天线远场辐射方向图。
附图说明
图1是本发明实施例采用的整体混合屏蔽室测试环境结构示意图。
图2是本发明实施例将混合屏蔽室测试环境下利用传统测试区域场补偿方法得到的待测天线方向图与理想待测天线方向图在三维及phi等于0°和90°的二维切面上进行对比的示意图。
图3是本发明实施例采用的球面波系数重构和补偿方法的算法流程图。
图4是本发明实施例将混合屏蔽室测试环境下利用本发明所提出的球面波系数重构补偿方法得到的待测天线方向图与理想待测天线方向图在三维及phi等于0°和90°的二维切面上进行对比的示意图。
图5是本发明实施例将理想待测天线辐射方向图与利用传统测试区域场补偿方法和利用所提出的球面波系数重构补偿方法得到的待测天线辐射方向图在phi等于0°和90°的二维切面上进行对比的示意图。
图6是本发明实施例分别将传统测试区域补偿法和本发明所述方法所得待测天线方向图与理论待测天线辐射方向图在phi=(0°-360°)和theta在phi=(0°-180°)的整个二维平面上进行误差对比的示意图。
具体实施方式
下面结合附图对本发明作进一步详细说明。
本发明实施例的整体测试环境如图1所示,混合屏蔽室配置有两个机械搅拌器101、102和一个转台103,转台可以进行旋转,从而使安装在转台上的天线沿X、Y和Z进行旋转。以天线自身建立三维直角坐标系,当天线的Z轴垂直转台放置并转动转台时,天线可以完成绕Z轴旋转的效果,同理当天线的Y轴和X轴垂直转台放置并转动转台时,天线可以实现绕Y和X轴的旋转效果。同时在测试环境中,安装有负载吸收体104,待测天线105和校准天线106,并且校准天线106的传输系数已知。
为验证所述球面波系数重构和补偿方法的有效性,本发明实施例对混合屏蔽室环境下的传统测试区域场补偿方法以及所述球面波系数重构和补偿方法分别进行验证:
1.对测试区域边界为电壁的多探头球面近场***所得的数据进行传统测试区域场补偿算法处理,得到补偿后的待测天线方向图;
2.对混合屏蔽室内的待测天线进行旋转,得到待测天线与校准天线间的相关系数,通过重构球面波系数和测试区域场补偿处理,得到改善后的待测天线方向图。
所述验证内容所需的数据包括混合屏蔽室内待测天线近场采样数据、待测天线与校准天线间的相关性系数矩阵以及待测天线和校准天线经过旋转后的球面波系数矩阵。
根据测试区域场补偿理论,可以利用式(2)将混合屏蔽室中含反射径干扰的待测天线105远场分量表示为传输系数与信源相乘的形式:
其中v表示待测天线上的发射信源,Tmn表示待测天线的传输系数,与式(1)中的和等价,均表示球面波生成函数。同理,通过式(1)(2)的求解同样可以得到所述校准天线106的远场电场分量以及球面波系数amn和bmn。
现将待测天线105在含多径反射影响下计算得到的球面波系数记为a'AUT,mn和b'AUT,mn,而待求的理论无干扰的待测天线球面波系数记为aAUT,mn和bAUT,mn,待求的理论无干扰的待测天线传输系数记为TAUT,mn;将所述已知的校准天线106传输系数记为TREF,mn,含多径反射影响的校准天线106球面波系数记为a'REF,mn和b'REF,mn,已知无干扰的信源信号记为v,含多径反射影响的信源信号记为v'。
通过式(1)(2)可以看出,含多径反射影响的信源信号v'可以表示为式(3)的形式:
(a'REF,mn,b'REF,mn)=TREF,mnv' (3)
而待测天线105和校准天线106是在相同的混合屏蔽室环境下进行的测试,故两种天线测试结果下的信号v'一致,因此通过式(4)可以根据所述已知的校准天线106传输系数TREF,mn和计算得到的校准天线106球面波系数a'REF,mn和b'REF,mn推得待测天线的理论传输系数TAUT,mn:
TAUT,mn=(a'AUT,mn,b'AUT,mn)((TREF,mn)+·(a'REF,mn,b'REF,mn))+ (4)
其中+表示矩阵的广义逆矩阵。
上述方法完成了传统测试区域场补偿算法,图2是传统补偿方法所得待测天线辐射方向图与理想待测天线辐射方向图进行对比的结果。
在上述补偿算法的基础上,本发明提供了如图3所示的球面波系数重构和补偿方法的流程示意图。在混合屏蔽室的环境下首先将所述待测天线放置在转台上,并在待测天线本身建立三维直角坐标系,分别在天线的X、Y和Z方向与转台水平面垂直时让转台转动360度,并且每一次的角度间隔为1度,同理将待测天线换成校准天线重复相同操作,就可得到如图3所述的待测天线和校准天线的旋转方向图。此时,需要根据式(5)计算旋转待测天线与校准天线之间的相关性系数以及校准天线自身的自相关系数:
其中amn1和bmn1表示含反射径影响的待测天线的球面波系数,amn2和bmn2表示含反射径影响的校准天线的球面波系数,T(α,β,γ)表示待测天线的旋转矩阵,通过T(α,β,γ)与amn2和bmn2的矩阵相乘,就可以得到图3所述校准天线的旋转方向图,进而利用式(5)得到旋转校准天线与待测天线之间的相关系数,记为ρ1。同理,将式(5)中的待测天线换成校准天线,即amn1和bmn1用amn2和bmn2替换,则可以得到图3所述校准天线的旋转方向图以及校准天线与其自身旋转方向图之间的自相关系数,记为ρ2。
为了在混响室中得到所述ρ1和ρ2,可以通过测量不同旋转位置处的待测天线和校准天线的S参数,并计算S参数之间的相关性,从而得到对应的相关系数ρ1和自相关系数ρ2,如式(6)所示:
其中,S21,k,REF和S21,k,AUT是发射天线到校准天线与待测天线之间的S参数,<S21,REF>和<S21,AUT>则对应上述两种参数的所有采样值的平均值,利用(6)式即可得到相关系数ρ1;若将(6)式中的S21,k,AUT和<S21,AUT>用S'21,k,REF和<S'21,REF>进行替换,其中S'21,k,REF和<S'21,REF>表示校准天线在未旋转状态下的S参数和对应的平均值,则可得到自相关系数ρ2。
下面利用上述已知的ρ1和ρ2相关性系数以及旋转矩阵T(α,β,γ),便可通过式(7)所示的反向重构法,得到对应图3所述的校准天线和待测天线的含反射径干扰的球面波系数:
其中Tα,Tβ和Tγ分别对应所述旋转矩阵T(α,β,γ)在X,Y和Z上的投影矩阵,ρs,α,ρs,β和ρs,γ分别表示相关性系数ρ1(下标s取1)或自相关系数ρ2(下标s取2)在X,Y和Z上的投影矩阵。利用广义逆矩阵的方法就可以反向推导出待测天线在含球面反射径干扰时的方向图。同理,将式(7)中的待测天线和校准天线互换,便可得到含反射径干扰的校准天线方向图。
现在通过测试区域场补偿的方法,对已求得的待测天线方向图进行补偿。在混合屏蔽室的条件下,待测天线和校准天线的测试环境保持一致,因此可以应用测试区域场补偿的思想,即将式(3)(4)中的信号v'看作混合屏蔽室中因球面反射径干扰所产生的信源信号,进而通过式(8)的逆矩阵补偿方式得到如图3所述的去除球面反射径干扰的待测天线传输系数:
TAUT,mn=(amn1,bmn1)((TREF,mn)+·(amn2,bmn2))+ (8)
式(8)中利用了上述计算所得的含反射径干扰的待测天线球面波系数amn1和bmn1,以及含反射径干扰的校准天线球面波系数amn2和bmn2,并利用所述已知的校准天线传输系数TREF,mn便可补偿得到改善后的待测天线方向图。
图4是将利用本发明方法重构补偿所得的待测天线方向图与理想待测天线方向图进行对比的示意图,图5则是将利用传统补偿方法和本发明方法所得的待测天线方向图跟理想待测天线方向图进行对比的示意图,其中实线表示理想待测方向图,点线表示传统方法的效果,虚线表示本发明方法的效果,通过效果曲线对比,可以看出,本发明所述方法对待测天线的主瓣恢复效果明显,并且能够较好地补偿因球面反射径干扰所造成的方向图损害。最后将所述两种方法得到的待测天线方向图与理论方向图在整个球面上作差,并映射到俯仰角θ和方位角的取值范围为(0°-180°)和(0°-360°)的二维平面上,从而得到如图6所示的误差分析图。由图6的误差对比可以看出,本发明所述的球面波系数重构和补偿方法在整个球面上的误差要小于传统测试区域场补偿技术,同时通过计算两种重构补偿方法在整个球面上的误差总和,可知传统方法的整体误差为389.9402dB,而本文所述方法的重构补偿误差为301.8276dB,因此通过本发明所述方法所得的待测天线方向图与理论方向图之间的误差更为准确。
本发明在混合屏蔽室中通过反向重构的方法降低反射径的影响,同时将测试区域补偿方法应用到混和屏蔽室的待测天线方向图补偿问题中,相关研究内容将会对完善混和屏蔽室下的天线测量技术具有十分积极的意义。
以上所述,仅为本发明较佳的具体实施方式,但本发明的保护范围并不局限于此,任何熟悉本技术领域的技术人员在本发明披露的技术范围内,可轻易想到的变化或替换,都应涵盖在本发明的保护范围之内。因此,本发明的保护范围应该以权利要求书的保护范围为准。
Claims (3)
1.一种基于混合屏蔽室的抗干扰天线方向图测量方法,其特征在于,包括以下步骤;
步骤1.选取两个天线,将其中一个天线作为待测天线(105),另一个天线作为校准天线(106),在标准暗室环境下测量校准天线(106)的辐射方向图并求得校准天线(106)的传输系数;
步骤2.在混合屏蔽室环境下,先将待测天线(105)放置在室内的转台(103)上,控制待测天线(105)进行旋转,得到待测天线(105)S参数矩阵;再将校准天线(106)放置在转台上进行旋转,得到校准天线(106)的S参数矩阵;
步骤3.计算待测天线(105)和校准天线(106)的相关性系数矩阵和校准天线(106)的自相关系数矩阵;
步骤4.利用球面波系数矩阵和相关性系数矩阵重构出待测天线的方向图,同时利用步骤1所述的校准天线(106)传输系数来补偿待测天线方向图,完成待测天线(105)重构和补偿工作;
所述步骤2中,转台(103)能够控制待测天线(105)和校准天线(106)沿X、Y和Z方向进行旋转,三个方向的旋转范围是0°-360°,角度间隔是1度;
根据测试区域场补偿理论,利用式(2)将混合屏蔽室中含反射径干扰的待测天线105远场分量表示为传输系数与信源相乘的形式:
其中v表示待测天线上的发射信源,Tmn表示待测天线的传输系数,与式(1)中的和等价,均表示球面波生成函数,同理,通过式(1)(2)的求解同样可以得到所述校准天线(106)的远场电场分量以及球面波系数amn和bmn;
将待测天线(105)在含多径反射影响下计算得到的球面波系数记为a'AUT,mn和b'AUT,mn,而待求的理论无干扰的待测天线球面波系数记为aAUT,mn和bAUT,mn,待求的理论无干扰的待测天线传输系数记为TAUT,mn;将所述已知的校准天线(106)传输系数记为TREF,mn,含多径反射影响的校准天线(106)球面波系数记为a'REF,mn和b'REF,mn,已知无干扰的信源信号记为v,含多径反射影响的信源信号记为v';
通过式(1)(2)可以看出,含多径反射影响的信源信号v'可以表示为式(3)的形式:
(a'REF,mn,b'REF,mn)=TREF,mnv' (3)
而待测天线(105)和校准天线(106)是在相同的混合屏蔽室环境下进行的测试,故两种天线测试结果下的信号v'一致,因此通过式(4)可以根据所述已知的校准天线(106)传输系数TREF,mn和计算得到的校准天线(106)球面波系数a'REF,mn和b'REF,mn推得待测天线的理论传输系数TAUT,mn:
TAUT,mn=(a'AUT,mn,b'AUT,mn)((TREF,mn)+·(a'REF,mn,b'REF,mn))+ (4)
其中+表示矩阵的广义逆矩阵。
2.根据权利要求1所述的一种基于混合屏蔽室的抗干扰天线方向图测量方法,其特征在于,所述步骤3中根据式(5)计算旋转待测天线(105)与校准天线(106)之间的相关性系数以及校准天线(106)自身的自相关系数:
其中amn1和bmn1表示含反射径影响的待测天线(105)的球面波系数,amn2和bmn2表示含反射径影响的校准天线(106)的球面波系数,T(α,β,γ)表示待测天线(105)的旋转矩阵,通过T(α,β,γ)与amn2和bmn2的矩阵相乘,得到校准天线(106)的旋转方向图,利用式(5)得到校准天线(106)与待测天线(105)之间的相关系数,记为ρ1;同理,将式(5)中的待测天线(105)换成校准天线(106),即amn1和bmn1用amn2和bmn2替换,得到校准天线(106)的旋转方向图以及校准天线(106)与其自身旋转方向图之间的自相关系数,记为ρ2;
为了在混响室中得到所述ρ1和ρ2,通过测量不同旋转位置处的待测天线(105)和校准天线(106)的S参数,并计算S参数之间的相关性,从而得到对应的相关系数ρ1和自相关系数ρ2,如式(6)所示:
其中,S21,k,REF和S21,k,AUT是发射天线到校准天线(106)与待测天线(105)之间的S参数,<S21,REF>和<S21,AUT>则对应S21,k,REF和S21,k,AUT两种参数的所有采样值的平均值,利用(6)式即可得到相关系数ρ1;若将(6)式中的S21,k,AUT和<S21,AUT>用S'21,k,REF和<S'21,REF>进行替换,其中S'21,k,REF和<S'21,REF>表示校准天线(106)在未旋转状态下的S参数和对应的平均值,则可得到自相关系数ρ2;
利用上述已知的ρ1和ρ2相关性系数以及旋转矩阵T(α,β,γ),通过式(7)所示的反向重构法,得到校准天线(106)和待测天线(105)的含反射径干扰的球面波系数:
其中Tα,Tβ和Tγ分别对应所述旋转矩阵T(α,β,γ)在X,Y和Z上的投影矩阵,ρs,α,ρs,β和ρs,γ分别表示相关性系数ρ1或自相关系数ρ2在X,Y和Z上的投影矩阵,利用广义逆矩阵的方法反向推导出待测天线(105)在含球面反射径干扰时的方向图,同理,将式(7)中的待测天线(105)和校准天线(106)互换,便可得到含反射径干扰的校准天线方向图。
3.根据权利要求1所述的一种基于混合屏蔽室的抗干扰天线方向图测量方法,其特征在于,所述步骤4中,对已求得的待测天线方向图进行补偿,在混合屏蔽室的条件下,待测天线(105)和校准天线(106)的测试环境保持一致,应用测试区域场补偿,即将式(3)(4)中的信号v'看作混合屏蔽室中因球面反射径干扰所产生的信源信号,进而通过式(8)的逆矩阵补偿方式得到去除球面反射径干扰的待测天线传输系数:
TAUT,mn=(amn1,bmn1)((TREF,mn)+·(amn2,bmn2))+ (8)
式(8)中利用了上述计算所得的含反射径干扰的待测天线球面波系数amn1和bmn1,以及含反射径干扰的校准天线球面波系数amn2和bmn2,并利用所述已知的校准天线传输系数TREF,mn补偿得到改善后的待测天线方向图。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111001299.9A CN113804985B (zh) | 2021-08-30 | 2021-08-30 | 一种基于混合屏蔽室的抗干扰天线方向图测量方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111001299.9A CN113804985B (zh) | 2021-08-30 | 2021-08-30 | 一种基于混合屏蔽室的抗干扰天线方向图测量方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113804985A CN113804985A (zh) | 2021-12-17 |
CN113804985B true CN113804985B (zh) | 2022-07-26 |
Family
ID=78942282
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202111001299.9A Active CN113804985B (zh) | 2021-08-30 | 2021-08-30 | 一种基于混合屏蔽室的抗干扰天线方向图测量方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113804985B (zh) |
Families Citing this family (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114487618B (zh) * | 2022-01-27 | 2022-08-23 | 北京航空航天大学 | 一种复合材料低频电磁参数等效提取装置及方法 |
CN115508626A (zh) * | 2022-10-25 | 2022-12-23 | 西安交通大学 | 一种在混响室中唯幅值测量的天线方向图重构方法及*** |
CN115524541B (zh) * | 2022-11-28 | 2023-03-10 | 陕西拾贝通讯技术有限公司 | 基于准远场模式滤波的任意极化线天线阵方向图测量方法 |
CN117590090B (zh) * | 2024-01-18 | 2024-04-02 | 中国计量科学研究院 | 一种电磁混响室场均匀性快速确认装置、方法及设备 |
CN117890689B (zh) * | 2024-03-13 | 2024-05-17 | 成都华兴大地科技有限公司 | 一种适用于无控幅相控阵的校准方法 |
Citations (17)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2006029906A (ja) * | 2004-07-14 | 2006-02-02 | Maspro Denkoh Corp | 妨害排除能力試験装置 |
CN102818942A (zh) * | 2012-08-24 | 2012-12-12 | 湖北航天技术研究院计量测试技术研究所 | 天线远场参数校准装置及校准方法 |
CN103364646A (zh) * | 2013-07-25 | 2013-10-23 | 合肥师范学院 | 微波暗室天线远场快速测量方法 |
CN104515909A (zh) * | 2015-01-20 | 2015-04-15 | 中国人民解放军国防科学技术大学 | 一种基于相关法的大天线方向图测量方法 |
CN105548729A (zh) * | 2016-02-22 | 2016-05-04 | 石家庄世联达科技有限公司 | 一种阵列天线辐射特性的快速测量方法 |
CN106291133A (zh) * | 2016-10-12 | 2017-01-04 | 中国科学院电子学研究所 | 一种uhf波段宽频带dbf阵列天线测试方法 |
CN106291129A (zh) * | 2015-06-01 | 2017-01-04 | 北京空间飞行器总体设计部 | 相控阵天线远场方向图测量方法 |
CN106546828A (zh) * | 2016-10-31 | 2017-03-29 | 西北工业大学 | 一种基于微波网络分析的球面多探头天线测试数据处理方法 |
CN107796995A (zh) * | 2017-11-30 | 2018-03-13 | 上海英恒电子有限公司 | 微波暗室及相应的天线测试*** |
CN109541330A (zh) * | 2018-11-12 | 2019-03-29 | 北京航空航天大学 | 一种平面波模拟器的阵列天线通道校准*** |
CN109581078A (zh) * | 2018-11-30 | 2019-04-05 | 电子科技大学 | 一种适用于半空间环境中天线的方向图测量***及方法 |
CN110260797A (zh) * | 2019-06-18 | 2019-09-20 | 西安交通大学 | 一种应用于恒/变速光栅信号的自适应滤波方法 |
CN111337758A (zh) * | 2020-03-11 | 2020-06-26 | 南京航空航天大学 | 一种基于混响室的天线辐射效率测量方法 |
CN111707877A (zh) * | 2020-07-20 | 2020-09-25 | 广东圣大电子有限公司 | 一种射频发射机杂散辐射测试***及测试方法 |
CN112130007A (zh) * | 2020-09-25 | 2020-12-25 | 四川九洲电器集团有限责任公司 | 一种天线远场方向图测试***及测试方法 |
KR20210009612A (ko) * | 2019-07-17 | 2021-01-27 | 재단법인 자동차융합기술원 | 알터네이터 및 전압조정기의 전자기적합성 성능평가 시스템 |
CN112462168A (zh) * | 2020-11-05 | 2021-03-09 | 西安交通大学 | 一种基于可电调吸波超表面的快速空口测试方法 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US11671144B2 (en) * | 2019-03-28 | 2023-06-06 | Intel Corporation | Near-field test apparatus for far-field antenna properties |
US11175337B2 (en) * | 2019-08-29 | 2021-11-16 | Rohde & Schwarz Gmbh & Co. Kg | Over-the-air measurement system |
-
2021
- 2021-08-30 CN CN202111001299.9A patent/CN113804985B/zh active Active
Patent Citations (17)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2006029906A (ja) * | 2004-07-14 | 2006-02-02 | Maspro Denkoh Corp | 妨害排除能力試験装置 |
CN102818942A (zh) * | 2012-08-24 | 2012-12-12 | 湖北航天技术研究院计量测试技术研究所 | 天线远场参数校准装置及校准方法 |
CN103364646A (zh) * | 2013-07-25 | 2013-10-23 | 合肥师范学院 | 微波暗室天线远场快速测量方法 |
CN104515909A (zh) * | 2015-01-20 | 2015-04-15 | 中国人民解放军国防科学技术大学 | 一种基于相关法的大天线方向图测量方法 |
CN106291129A (zh) * | 2015-06-01 | 2017-01-04 | 北京空间飞行器总体设计部 | 相控阵天线远场方向图测量方法 |
CN105548729A (zh) * | 2016-02-22 | 2016-05-04 | 石家庄世联达科技有限公司 | 一种阵列天线辐射特性的快速测量方法 |
CN106291133A (zh) * | 2016-10-12 | 2017-01-04 | 中国科学院电子学研究所 | 一种uhf波段宽频带dbf阵列天线测试方法 |
CN106546828A (zh) * | 2016-10-31 | 2017-03-29 | 西北工业大学 | 一种基于微波网络分析的球面多探头天线测试数据处理方法 |
CN107796995A (zh) * | 2017-11-30 | 2018-03-13 | 上海英恒电子有限公司 | 微波暗室及相应的天线测试*** |
CN109541330A (zh) * | 2018-11-12 | 2019-03-29 | 北京航空航天大学 | 一种平面波模拟器的阵列天线通道校准*** |
CN109581078A (zh) * | 2018-11-30 | 2019-04-05 | 电子科技大学 | 一种适用于半空间环境中天线的方向图测量***及方法 |
CN110260797A (zh) * | 2019-06-18 | 2019-09-20 | 西安交通大学 | 一种应用于恒/变速光栅信号的自适应滤波方法 |
KR20210009612A (ko) * | 2019-07-17 | 2021-01-27 | 재단법인 자동차융합기술원 | 알터네이터 및 전압조정기의 전자기적합성 성능평가 시스템 |
CN111337758A (zh) * | 2020-03-11 | 2020-06-26 | 南京航空航天大学 | 一种基于混响室的天线辐射效率测量方法 |
CN111707877A (zh) * | 2020-07-20 | 2020-09-25 | 广东圣大电子有限公司 | 一种射频发射机杂散辐射测试***及测试方法 |
CN112130007A (zh) * | 2020-09-25 | 2020-12-25 | 四川九洲电器集团有限责任公司 | 一种天线远场方向图测试***及测试方法 |
CN112462168A (zh) * | 2020-11-05 | 2021-03-09 | 西安交通大学 | 一种基于可电调吸波超表面的快速空口测试方法 |
Non-Patent Citations (4)
Title |
---|
A Method of Antenna Gain Testing without Standard Gain Antenna;Kuo Gao等;《 2020 International Conference on Sensing, Measurement & Data Analytics in the era of Artificial Intelligence (ICSMD)》;20201124;全文 * |
Evaluation of the Purity of OAM Modes Using the Reverberation Chamber Technique;Wei Xue等;《2020 14th European Conference on Antennas and Propagation (EuCAP)》;20200708;全文 * |
低增益GPS天线方向图的微波暗室测试方法;贾耀成等;《西安工业大学学报》;20150831;第35卷(第8期);全文 * |
混响室内空间特性测量;王玥等;《中国传媒大学学报自然科学版》;20090330;第16卷(第1期);全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN113804985A (zh) | 2021-12-17 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN113804985B (zh) | 一种基于混合屏蔽室的抗干扰天线方向图测量方法 | |
Kildal et al. | Characterization of reverberation chambers for OTA measurements of wireless devices: Physical formulations of channel matrix and new uncertainty formula | |
US8502546B2 (en) | Multichannel absorberless near field measurement system | |
Newell | Error analysis techniques for planar near-field measurements | |
CN110618408B (zh) | 一种精密测距***天线相位中心的***标定方法 | |
US9615274B2 (en) | Plane wave generation within a small volume of space for evaluation of wireless devices | |
CN113225147A (zh) | 一种阵列天线总辐射功率的测量方法、装置和*** | |
CN111366793B (zh) | 一种减小截断误差的平面天线近场测量方法及*** | |
CN107783086B (zh) | 用于诊断天线阵口径幅相场的畸变位置的方法 | |
WO2009046516A1 (en) | Multichannel absorberless near field measurement system | |
CN112083234A (zh) | 阵列天线总辐射功率测量方法、装置以及计算机存储介质 | |
CN111766455B (zh) | 基于口径电流法的相控阵天线方向图预测方法及*** | |
US7336230B2 (en) | Method for measuring electromagnetic radiation pattern and gain of radiator using term waveguide | |
CN113381187B (zh) | 球面相控阵天线坐标远近场比较修正方法 | |
Expósito et al. | Uncertainty assessment of a small rectangular anechoic chamber: From design to operation | |
Toivanen et al. | Modified test zone field compensation for small-antenna measurements | |
Joy | Near-field qualification methodology | |
Leatherwood et al. | Plane wave, pattern subtraction, range compensation | |
Dobychina et al. | Unmanned aerial vehicle antenna measurement using anechoic chamber | |
CN111385036B (zh) | 无线设备的射频性能测试方法、装置及测试仪 | |
KR100284134B1 (ko) | 지.티.이.엠셀출력단을이용한복사전계측정방법 | |
CN111487474B (zh) | 数值孪生电磁测量*** | |
CN104852775A (zh) | 一种大型相控阵天线波束性能的快速评估方法 | |
CN115508626A (zh) | 一种在混响室中唯幅值测量的天线方向图重构方法及*** | |
Fedorov et al. | Comparison of the Measurement Accuracy of Material Sample Specular Reflection Coefficient for Two Types of Measuring Facilities |
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 | ||
TR01 | Transfer of patent right |
Effective date of registration: 20240702 Address after: D503, Research Building 1, No. 6, Kejiao Road, Shishan Town, Nanhai District, Foshan City, Guangdong Province, 528225 Patentee after: Foshan lanpuda Technology Co.,Ltd. Country or region after: China Address before: Beilin District Xianning West Road 710049, Shaanxi city of Xi'an province No. 28 Patentee before: XI'AN JIAOTONG University Country or region before: China |