发明名称 一种基于微波网络分析的球面多探头天线测试数据处理方法
摘要 本发明涉及一种基于微波网络分析的球面多探头天线测试数据的快速处理方法,基于微波网络分析的思想,将待测天线及探头系统看成是一个开放的二端口网络,将激励等效为网络的输入,将探头接收信号等效为输出,通过坐标旋转,严格的推导出天线的传输方程,再通过探头上的接收信号和不同场区的探头响应常数,求解出待测天线的远场方向图。大幅度提高了程序运行速度,缩短了天线测试数据的处理时间。
申请公布号 CN106546828A 申请公布日期 2017.03.29
申请号 CN201610926411.2 申请日期 2016.10.31
申请人 西北工业大学 发明人 李南京;郭丽芳;胡楚锋;陈卫军;郭淑霞;郭鹏;潘世洲
分类号 G01R29/10(2006.01)I 主分类号 G01R29/10(2006.01)I
代理机构 西北工业大学专利中心 61204 代理人 刘新琼
主权项 一种基于微波网络分析的球面多探头天线测试数据的快速处理方法,其特征在于步骤如下:步骤1:利用一个已知特性的探头分别当其为θ极化和<img file="FDA0001137023500000011.GIF" wi="38" he="46" />极化时在距离待测天线r<sub>0</sub>=3~10λ的空间球面上测量得到待测天线的近场数据w(A,χ<sub>1</sub>,θ,φ)和w(A,χ<sub>2</sub>,θ,φ),并带入式(1)计算得到w<sub>μ</sub>(A,θ,φ):<maths num="0001"><math><![CDATA[<mrow><msub><mi>w</mi><mi>&mu;</mi></msub><mrow><mo>(</mo><mi>A</mi><mo>,</mo><mi>&theta;</mi><mo>,</mo><mi>&phi;</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mn>1</mn><mrow><mn>2</mn><mi>&pi;</mi></mrow></mfrac><munderover><mo>&Integral;</mo><mrow><mi>&chi;</mi><mo>=</mo><mn>0</mn></mrow><mrow><mn>2</mn><mi>&pi;</mi></mrow></munderover><mi>w</mi><mrow><mo>(</mo><mi>A</mi><mo>,</mo><mi>&chi;</mi><mo>,</mo><mi>&theta;</mi><mo>,</mo><mi>&phi;</mi><mo>)</mo></mrow><msup><mi>e</mi><mrow><mo>-</mo><mi>i</mi><mi>&mu;</mi><mi>&chi;</mi></mrow></msup><mi>d</mi><mi>&chi;</mi><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001137023500000012.GIF" wi="1406" he="151" /></maths>其中,w<sub>μ</sub>(A,θ,φ)是w(A,χ,θ,φ)关于χ有限傅里叶序列的展开系数,χ为探头旋转角,A为探头到待测天线的距离,θ为探头相对于待测天线的俯仰角,<img file="FDA0001137023500000013.GIF" wi="38" he="46" />为探头相对于待测天线的方位角,λ为待测天线的工作波长,μ=±1;步骤2:将步骤1中求得的w<sub>μ</sub>(A,θ,φ)带入式(2),并用快速傅里叶逆变换计算此积分项得到w<sub>μm</sub>(A,θ):<maths num="0002"><math><![CDATA[<mrow><msub><mi>w</mi><mrow><mi>&mu;</mi><mi>m</mi></mrow></msub><mrow><mo>(</mo><mi>A</mi><mo>,</mo><mi>&theta;</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mn>1</mn><mrow><mn>2</mn><mi>&pi;</mi></mrow></mfrac><munderover><mo>&Integral;</mo><mrow><mi>&phi;</mi><mo>=</mo><mn>0</mn></mrow><mrow><mn>2</mn><mi>&pi;</mi></mrow></munderover><msub><mi>w</mi><mi>&mu;</mi></msub><mrow><mo>(</mo><mi>A</mi><mo>,</mo><mi>&theta;</mi><mo>,</mo><mi>&phi;</mi><mo>)</mo></mrow><msup><mi>e</mi><mrow><mo>-</mo><mi>i</mi><mi>m</mi><mi>&phi;</mi></mrow></msup><mi>d</mi><mi>&phi;</mi><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001137023500000014.GIF" wi="1355" he="151" /></maths>其中,0≤n≤N,‑n≤m≤n,截断数N≥kr<sub>0</sub>+10,传播常数<img file="FDA0001137023500000015.GIF" wi="179" he="111" />r<sub>0</sub>为能包围待测天线的最小球面的半径;步骤3:将步骤2中对<img file="FDA0001137023500000016.GIF" wi="35" he="45" />的积分计算结果w<sub>μm</sub>(A,θ)带入式(3),用快速傅里叶逆变换计算此积分项得到<img file="FDA0001137023500000017.GIF" wi="179" he="63" /><maths num="0003"><math><![CDATA[<mrow><msubsup><mi>w</mi><mrow><mi>&mu;</mi><mi>m</mi></mrow><mi>n</mi></msubsup><mrow><mo>(</mo><mi>A</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mrow><mn>2</mn><mi>n</mi><mo>+</mo><mn>1</mn></mrow><mn>2</mn></mfrac><msubsup><mo>&Integral;</mo><mrow><mi>&theta;</mi><mo>=</mo><mn>0</mn></mrow><mi>&pi;</mi></msubsup><msub><mi>w</mi><mrow><mi>&mu;</mi><mi>m</mi></mrow></msub><mrow><mo>(</mo><mi>A</mi><mo>,</mo><mi>&theta;</mi><mo>)</mo></mrow><msubsup><mi>d</mi><mrow><mi>&mu;</mi><mi>m</mi></mrow><mi>n</mi></msubsup><mrow><mo>(</mo><mi>&theta;</mi><mo>)</mo></mrow><mi>sin</mi><mi>&theta;</mi><mi>d</mi><mi>&theta;</mi><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>3</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001137023500000018.GIF" wi="1317" he="119" /></maths>其中<img file="FDA0001137023500000019.GIF" wi="139" he="70" />为旋转系数:<maths num="0004"><math><![CDATA[<mrow><msubsup><mi>d</mi><mrow><mi>&mu;</mi><mi>m</mi></mrow><mi>n</mi></msubsup><mrow><mo>(</mo><mi>&theta;</mi><mo>)</mo></mrow><mo>=</mo><msup><mi>i</mi><mrow><mi>&mu;</mi><mo>-</mo><mi>m</mi></mrow></msup><munderover><mo>&Sigma;</mo><mrow><mi>m</mi><mo>=</mo><mo>-</mo><mi>n</mi></mrow><mi>n</mi></munderover><msubsup><mi>&Delta;</mi><mrow><msup><mi>m</mi><mo>&prime;</mo></msup><mi>&mu;</mi></mrow><mi>n</mi></msubsup><msubsup><mi>&Delta;</mi><mrow><msup><mi>m</mi><mo>&prime;</mo></msup><mi>m</mi></mrow><mi>n</mi></msubsup><msup><mi>e</mi><mrow><mo>-</mo><msup><mi>im</mi><mo>&prime;</mo></msup><mi>&theta;</mi></mrow></msup><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>4</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA00011370235000000110.GIF" wi="1181" he="135" /></maths>式中<img file="FDA00011370235000000111.GIF" wi="91" he="63" />和<img file="FDA00011370235000000112.GIF" wi="91" he="62" />为三维常数矩阵,i为虚数单位,‑n≤m′≤n;步骤4:将步骤3中对θ的积分计算结果<img file="FDA00011370235000000113.GIF" wi="150" he="69" />和近场探头响应常数P<sub>sμn</sub>(kA)带入式(5)即得到传输方程的传输系数vT<sub>1mn</sub>和vT<sub>2mn</sub>:<maths num="0005"><math><![CDATA[<mrow><mfenced open = "{" close = ""><mtable><mtr><mtd><mrow><msub><mi>vT</mi><mrow><mn>1</mn><mi>m</mi><mi>n</mi></mrow></msub><msub><mi>P</mi><mrow><mn>11</mn><mi>n</mi></mrow></msub><mrow><mo>(</mo><mi>k</mi><mi>A</mi><mo>)</mo></mrow><mo>+</mo><msub><mi>vT</mi><mrow><mn>2</mn><mi>m</mi><mi>n</mi></mrow></msub><msub><mi>P</mi><mrow><mn>21</mn><mi>n</mi></mrow></msub><mrow><mo>(</mo><mi>k</mi><mi>A</mi><mo>)</mo></mrow><mo>=</mo><msubsup><mi>w</mi><mrow><mn>1</mn><mi>m</mi></mrow><mi>n</mi></msubsup><mrow><mo>(</mo><mi>A</mi><mo>)</mo></mrow></mrow></mtd></mtr><mtr><mtd><mrow><msub><mi>vT</mi><mrow><mn>1</mn><mi>m</mi><mi>n</mi></mrow></msub><msub><mi>P</mi><mrow><mn>1</mn><mo>,</mo><mo>-</mo><mn>1</mn><mo>,</mo><mi>n</mi></mrow></msub><mrow><mo>(</mo><mi>k</mi><mi>A</mi><mo>)</mo></mrow><mo>+</mo><msub><mi>vT</mi><mrow><mn>2</mn><mi>m</mi><mi>n</mi></mrow></msub><msub><mi>P</mi><mrow><mn>2</mn><mo>,</mo><mo>-</mo><mn>1</mn><mo>,</mo><mi>n</mi></mrow></msub><mrow><mo>(</mo><mi>k</mi><mi>A</mi><mo>)</mo></mrow><mo>=</mo><msubsup><mi>w</mi><mrow><mo>-</mo><mn>1</mn><mi>m</mi></mrow><mi>n</mi></msubsup><mrow><mo>(</mo><mi>A</mi><mo>)</mo></mrow></mrow></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>5</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001137023500000021.GIF" wi="1284" he="165" /></maths>以m,n,θ为变量构成三维矩阵vT<sub>1mn</sub>和vT<sub>2mn</sub>,其中s=1,2,0≤n≤N,‑n≤m≤n,截断数N≥kr<sub>0</sub>+10,传播常数<img file="FDA0001137023500000022.GIF" wi="174" he="119" />r<sub>0</sub>为能包围待测天线的最小球面的半径,μ=±1;步骤5:将三维矩阵vT<sub>1mn</sub>和vT<sub>2mn</sub>的计算结果和远场探头响应常数<img file="FDA0001137023500000023.GIF" wi="163" he="62" />代入式(6)得到归一化远场探头信号:<maths num="0006"><math><![CDATA[<mrow><mi>W</mi><mrow><mo>(</mo><mi>A</mi><mo>,</mo><mi>&chi;</mi><mo>,</mo><mi>&theta;</mi><mo>,</mo><mi>&phi;</mi><mo>)</mo></mrow><mo>=</mo><mi>v</mi><munder><munder><mo>&Sigma;</mo><mrow><mi>s</mi><mi>m</mi><mi>n</mi></mrow></munder><mi>&mu;</mi></munder><msub><mi>T</mi><mrow><mi>s</mi><mi>m</mi><mi>n</mi></mrow></msub><msup><mi>e</mi><mrow><mi>i</mi><mi>m</mi><mi>&phi;</mi></mrow></msup><msubsup><mi>d</mi><mrow><mi>&mu;</mi><mi>m</mi></mrow><mi>n</mi></msubsup><mrow><mo>(</mo><mi>&theta;</mi><mo>)</mo></mrow><msup><mi>e</mi><mrow><mi>i</mi><mi>&mu;</mi><mi>&chi;</mi></mrow></msup><msubsup><mi>P</mi><mrow><mi>s</mi><mi>&mu;</mi><mi>n</mi></mrow><mi>&infin;</mi></msubsup><mrow><mo>(</mo><mi>k</mi><mi>A</mi><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>6</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001137023500000024.GIF" wi="1334" he="143" /></maths>步骤6:将归一化远场探头信号带入式(7)即得到待测天线远场方向图:<maths num="0007"><math><![CDATA[<mrow><mover><mi>K</mi><mo>&RightArrow;</mo></mover><mrow><mo>(</mo><mi>&theta;</mi><mo>,</mo><mi>&phi;</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mrow><mn>2</mn><msqrt><mn>6</mn></msqrt></mrow><mrow><mn>3</mn><mi>v</mi></mrow></mfrac><mrow><mo>(</mo><mi>W</mi><mo>(</mo><mi>A</mi><mo>,</mo><msub><mi>&chi;</mi><mn>1</mn></msub><mo>,</mo><mi>&theta;</mi><mo>,</mo><mi>&phi;</mi><mo>)</mo></mrow><mover><mi>&theta;</mi><mo>^</mo></mover><mo>+</mo><mi>W</mi><mrow><mo>(</mo><mi>A</mi><mo>,</mo><msub><mi>&chi;</mi><mn>2</mn></msub><mo>,</mo><mi>&theta;</mi><mo>,</mo><mi>&phi;</mi><mo>)</mo></mrow><mover><mi>&phi;</mi><mo>^</mo></mover><mo>)</mo><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>7</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001137023500000025.GIF" wi="1389" he="135" /></maths>其中,ν是待测天线的输入信号,<img file="FDA0001137023500000026.GIF" wi="33" he="63" />为球坐标系下θ方向的单位向量,<img file="FDA0001137023500000027.GIF" wi="38" he="70" />为球坐标系下φ方向的单位向量。
地址 710072 陕西省西安市友谊西路127号