主权项 |
一种基于微波网络分析的球面多探头天线测试数据的快速处理方法,其特征在于步骤如下:步骤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>μ</mi></msub><mrow><mo>(</mo><mi>A</mi><mo>,</mo><mi>θ</mi><mo>,</mo><mi>φ</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mn>1</mn><mrow><mn>2</mn><mi>π</mi></mrow></mfrac><munderover><mo>∫</mo><mrow><mi>χ</mi><mo>=</mo><mn>0</mn></mrow><mrow><mn>2</mn><mi>π</mi></mrow></munderover><mi>w</mi><mrow><mo>(</mo><mi>A</mi><mo>,</mo><mi>χ</mi><mo>,</mo><mi>θ</mi><mo>,</mo><mi>φ</mi><mo>)</mo></mrow><msup><mi>e</mi><mrow><mo>-</mo><mi>i</mi><mi>μ</mi><mi>χ</mi></mrow></msup><mi>d</mi><mi>χ</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>μ</mi><mi>m</mi></mrow></msub><mrow><mo>(</mo><mi>A</mi><mo>,</mo><mi>θ</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mn>1</mn><mrow><mn>2</mn><mi>π</mi></mrow></mfrac><munderover><mo>∫</mo><mrow><mi>φ</mi><mo>=</mo><mn>0</mn></mrow><mrow><mn>2</mn><mi>π</mi></mrow></munderover><msub><mi>w</mi><mi>μ</mi></msub><mrow><mo>(</mo><mi>A</mi><mo>,</mo><mi>θ</mi><mo>,</mo><mi>φ</mi><mo>)</mo></mrow><msup><mi>e</mi><mrow><mo>-</mo><mi>i</mi><mi>m</mi><mi>φ</mi></mrow></msup><mi>d</mi><mi>φ</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>μ</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>∫</mo><mrow><mi>θ</mi><mo>=</mo><mn>0</mn></mrow><mi>π</mi></msubsup><msub><mi>w</mi><mrow><mi>μ</mi><mi>m</mi></mrow></msub><mrow><mo>(</mo><mi>A</mi><mo>,</mo><mi>θ</mi><mo>)</mo></mrow><msubsup><mi>d</mi><mrow><mi>μ</mi><mi>m</mi></mrow><mi>n</mi></msubsup><mrow><mo>(</mo><mi>θ</mi><mo>)</mo></mrow><mi>sin</mi><mi>θ</mi><mi>d</mi><mi>θ</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>μ</mi><mi>m</mi></mrow><mi>n</mi></msubsup><mrow><mo>(</mo><mi>θ</mi><mo>)</mo></mrow><mo>=</mo><msup><mi>i</mi><mrow><mi>μ</mi><mo>-</mo><mi>m</mi></mrow></msup><munderover><mo>Σ</mo><mrow><mi>m</mi><mo>=</mo><mo>-</mo><mi>n</mi></mrow><mi>n</mi></munderover><msubsup><mi>Δ</mi><mrow><msup><mi>m</mi><mo>′</mo></msup><mi>μ</mi></mrow><mi>n</mi></msubsup><msubsup><mi>Δ</mi><mrow><msup><mi>m</mi><mo>′</mo></msup><mi>m</mi></mrow><mi>n</mi></msubsup><msup><mi>e</mi><mrow><mo>-</mo><msup><mi>im</mi><mo>′</mo></msup><mi>θ</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>χ</mi><mo>,</mo><mi>θ</mi><mo>,</mo><mi>φ</mi><mo>)</mo></mrow><mo>=</mo><mi>v</mi><munder><munder><mo>Σ</mo><mrow><mi>s</mi><mi>m</mi><mi>n</mi></mrow></munder><mi>μ</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>φ</mi></mrow></msup><msubsup><mi>d</mi><mrow><mi>μ</mi><mi>m</mi></mrow><mi>n</mi></msubsup><mrow><mo>(</mo><mi>θ</mi><mo>)</mo></mrow><msup><mi>e</mi><mrow><mi>i</mi><mi>μ</mi><mi>χ</mi></mrow></msup><msubsup><mi>P</mi><mrow><mi>s</mi><mi>μ</mi><mi>n</mi></mrow><mi>∞</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>→</mo></mover><mrow><mo>(</mo><mi>θ</mi><mo>,</mo><mi>φ</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>χ</mi><mn>1</mn></msub><mo>,</mo><mi>θ</mi><mo>,</mo><mi>φ</mi><mo>)</mo></mrow><mover><mi>θ</mi><mo>^</mo></mover><mo>+</mo><mi>W</mi><mrow><mo>(</mo><mi>A</mi><mo>,</mo><msub><mi>χ</mi><mn>2</mn></msub><mo>,</mo><mi>θ</mi><mo>,</mo><mi>φ</mi><mo>)</mo></mrow><mover><mi>φ</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" />为球坐标系下φ方向的单位向量。 |