留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

姓名
邮箱
手机号码
标题
留言内容
验证码

基于微波网络参数的非均匀传输线终端响应研究

鲁合德 朱峰 刘庆想 李鑫

鲁合德, 朱峰, 刘庆想, 李鑫. 基于微波网络参数的非均匀传输线终端响应研究[J]. 电子科技大学学报, 2018, 47(5): 698-704. doi: 10.3969/j.issn.1001-0548.2018.05.010
引用本文: 鲁合德, 朱峰, 刘庆想, 李鑫. 基于微波网络参数的非均匀传输线终端响应研究[J]. 电子科技大学学报, 2018, 47(5): 698-704. doi: 10.3969/j.issn.1001-0548.2018.05.010
LU He-de, ZHU Feng, LIU Qing-xiang, LI Xin. Research on Terminal Response of Abruptly Non-Uniform Transmission Line Based on Microwave Network Parameters[J]. Journal of University of Electronic Science and Technology of China, 2018, 47(5): 698-704. doi: 10.3969/j.issn.1001-0548.2018.05.010
Citation: LU He-de, ZHU Feng, LIU Qing-xiang, LI Xin. Research on Terminal Response of Abruptly Non-Uniform Transmission Line Based on Microwave Network Parameters[J]. Journal of University of Electronic Science and Technology of China, 2018, 47(5): 698-704. doi: 10.3969/j.issn.1001-0548.2018.05.010

基于微波网络参数的非均匀传输线终端响应研究

doi: 10.3969/j.issn.1001-0548.2018.05.010
基金项目: 

国家自然科学基金 U1434203

详细信息
    作者简介:

    鲁合德(1988-), 男, 博士生, 主要从事电磁兼容场线耦合相关方面的研究

  • 中图分类号: TN97

Research on Terminal Response of Abruptly Non-Uniform Transmission Line Based on Microwave Network Parameters

图(7) / 表(2)
计量
  • 文章访问数:  7138
  • HTML全文浏览量:  2179
  • PDF下载量:  79
  • 被引次数: 0
出版历程
  • 收稿日期:  2017-06-21
  • 修回日期:  2018-04-08
  • 刊出日期:  2018-09-30

基于微波网络参数的非均匀传输线终端响应研究

doi: 10.3969/j.issn.1001-0548.2018.05.010
    基金项目:

    国家自然科学基金 U1434203

    作者简介:

    鲁合德(1988-), 男, 博士生, 主要从事电磁兼容场线耦合相关方面的研究

  • 中图分类号: TN97

摘要: 为了研究电磁场对有损大地非均匀传输线终端响应的影响,以传输线模型为基础,应用微波网络参数转移矩阵推导了平面波入射条件下有损大地上突变型非均匀传输线终端响应公式。利用该公式可精确计算突变型非均匀传输线的终端响应大小。为了证明该理论的正确性,对有损大地上接有转换器的突变型非均匀传输线结构进行了理论分析计算。分别对不同电导率有损大地的情况进行了理论分析计算:不同电导率对终端感应电流的影响相差小于3 dB。最后结合实际选取大地电导率的典型值0.01S/m对该非均匀传输线结构终端电流响应进行了实测。实测结果与理论结果对比,两者吻合较好,因此该微波网络转移矩阵方法特别适合对突变型非均匀传输线场线耦合终端响应的求解计算。

English Abstract

鲁合德, 朱峰, 刘庆想, 李鑫. 基于微波网络参数的非均匀传输线终端响应研究[J]. 电子科技大学学报, 2018, 47(5): 698-704. doi: 10.3969/j.issn.1001-0548.2018.05.010
引用本文: 鲁合德, 朱峰, 刘庆想, 李鑫. 基于微波网络参数的非均匀传输线终端响应研究[J]. 电子科技大学学报, 2018, 47(5): 698-704. doi: 10.3969/j.issn.1001-0548.2018.05.010
LU He-de, ZHU Feng, LIU Qing-xiang, LI Xin. Research on Terminal Response of Abruptly Non-Uniform Transmission Line Based on Microwave Network Parameters[J]. Journal of University of Electronic Science and Technology of China, 2018, 47(5): 698-704. doi: 10.3969/j.issn.1001-0548.2018.05.010
Citation: LU He-de, ZHU Feng, LIU Qing-xiang, LI Xin. Research on Terminal Response of Abruptly Non-Uniform Transmission Line Based on Microwave Network Parameters[J]. Journal of University of Electronic Science and Technology of China, 2018, 47(5): 698-704. doi: 10.3969/j.issn.1001-0548.2018.05.010
  • 无论是电力系统还是电子信息系统,都离不开各类传输线来进行能量或信号的传输,而传输线系统极易受到外部电磁场的影响。外部电磁场在传输线上产生感应电流,从而将能量施加到传输线上,影响传输线正常的响应,因此电磁场与传输线之间的耦合作用是电磁兼容的关键问题之一[1-3]。传输线模型在电磁兼容计算模型中扮演着重要角色[4],如多导体传输线间的串扰[5-7]、外电磁场照射下的传输线响应[8-10]等一直是电磁兼容研究的热点。对于场线耦合的研究,文献[6, 11]分别给出了不同的计算模型;文献[12-13]研究了大地对传输线的影响效果,得到地电阻的不同求解形式。文献[14]对传输线和散射模型在计算高空核电磁脉冲对地面上长导线的响应进行了比较,并对两种模型进行了数值计算。文献[15-16]应用细线天线理论推导了电磁场对有损地面上有限长传输线场线耦合电报方程的一般形式及其时域方程,并同传输线模型进行了对比计算。文献[17]使用时域BLT方程对非线性终端传输线进行了分析研究。文献[18]对不同传输线模型进行了比较分析。文献[19]利用电磁拓扑学使用扩展的BLT方程得到了电偶极子激励下架空传输线的瞬态响应。文献[20-21]对有损大地架空屏蔽线缆耦合特性进行分析研究,并对传输线共模和异模分量进行了分析计算。

    不同传输线之间经常会使用转换器进行连接,而转换器的加入会造成传输线特性的突变,使得传输线整体不均匀。文献[122-25]对非均匀传输线的研究作了大量工作。文献[1]基于Sturm-Liouville特征函数解及有限差分法得到有场激励情况下的非均匀传输线响应。文献[26]对非均匀传输线的场路特性进行了详细的研究。文献[27]利用等效电路法解决电磁脉冲注入下有损非均匀传输线及其网络端接非线性负载的时域响应问题。

    微波网络是用参数化矩阵表示任意微波端口的特性,不用考虑网络内部详细的场结构,只需要知道端口处的参数即可,对于只关心终端响应的情况非常适用。转移矩阵(或链参数矩阵)是微波网络参数矩阵的一种,它给出的是传输线输入端电压电流与输出端电压电流间的关系,对于由几部分串联而成的传输线,其总转移矩阵是每个部分转移矩阵的乘积[5, 28]。因此对于突变型非均匀传输线可以等效为几段均匀传输线的级联,通过转移矩阵相乘得到总转移矩阵,进而得到最终的终端响应。总结文献发现,对于非均匀传输线的计算,文献[5]首先提出使用微波网络的转移矩阵法计算突变非均匀传输线,只是针对无源情况下的突变非均匀传输线的计算。有源传输线的转移矩阵法的理论计算及实验研究却很少涉及。本文以传输线模型为基础,结合微波网络的转移矩阵方法,推导了平面波照射下有损大地上突变型非均匀传输线的终端响应。然后结合推导公式

    在不同电导率情况下,对接有转换器的传输线结构进行了理论计算分析。最后选取大地电导率典型值0.01 S/m,通过实验实测了接有转换器的传输线终端电流响应结果,与通过转移矩阵级联得到的结果对比,两者吻合较好。该计算方法特别适合突变型非均匀传输线终端响应的计算,为计算接有转接器的非均匀传输线提供了理论方法。

    • 微波网络的转移矩阵是用来描述二端口网络输入端口电压和电流与输出端口电压和电流关系的。二端口网络结构如图 1所示。转移矩阵元素无明确的物理意义,但特别适合分析二端口网络,尤其对于级联网络的计算,可以直接通过矩阵相乘得到,其转移矩阵为[5]

      图  1  二端口网络结构

      $$\left[ {\begin{array}{*{20}{c}} {{V_2}}\\ {{I_2}} \end{array}} \right] = {\mathit{\boldsymbol{ \boldsymbol{\varPhi} }}_1}\left[ {\begin{array}{*{20}{c}} {{V_1}}\\ {{I_1}} \end{array}} \right]$$ (1)

      式中,Φ1为二端口网络的转移矩阵,是一个二维方阵。对于图 2所示的级联形式,可以用转移矩阵表示为:

      图  2  N个二端口网络的级联

      $$\left[ {\begin{array}{*{20}{c}} {{V_N}}\\ {{I_N}} \end{array}} \right] = \prod\limits_{i = 1}^N {{\mathit{\boldsymbol{ \boldsymbol{\varPhi} }}_i}} \left[ {\begin{array}{*{20}{c}} {{V_1}}\\ {{I_1}} \end{array}} \right] = {\mathit{\boldsymbol{ \boldsymbol{\varPhi} }}_N}\left[ {\begin{array}{*{20}{c}} {{V_1}}\\ {{I_1}} \end{array}} \right]$$ (2)

      对于突变型非均匀传输线,可以进行分段处理,每一段都是均匀传输线,可以得到每段均匀传输线的转移矩阵形式,然后各均匀段的转移矩阵相乘,进而得到突变型非均匀传输线终端响应结果。

    • 计算模型如图 3所示,传输线长度为L,导线半径为a,导线位于有损地面上方h处,其中大地的电导率为σg,相对介电常数为εr。入射波为任意极化平面波以任意角度入射,入射电磁场在导线上感应出沿线的分布电流源和电压源。

      图  3  平面波照射有损大地传输线模型图

      图 3所示的传输线模型的电报方程为[4]

      $$\frac{{{\rm{d}}V(x)}}{{{\rm{d}}x}} + Z'_\omega I(x) = {V_s}(x)$$ (3)
      $$\frac{{{\rm{d}}I(x)}}{{{\rm{d}}x}} + Y'_\omega V(x) = {I_s}(x)$$ (4)

      式中,${Z'_\omega }$和${Y'_\omega }$分别为有损大地情况下传输线的单位长度特性阻抗和单位长度特性导纳;Vs(x)和Is(x)是电磁场在导线表面激励起的分布源,通过求解麦克斯韦方程组可以得到。根据电路理论可以得到传输线的边界条件为:

      $$V(0) = - {R_1}I(0)$$ (5)
      $$V(L) = {R_2}I(L)$$ (6)

      由电报方程式(3)和式(4)及它们的边界条件可以使用格林函数方法得到电报方程的解:通过传输线电流和电压的表达式求出各点的电压和电流,以此作为格林函数的激励源在源分布范围内积分,就可以得到传输线上任意一点的电压和电流值。

      当只考虑传输线终端电压电流的响应时,将方程式(3)和式(4)的解写成转移矩阵的形式如下[5, 28]

      $$\left[ {\begin{array}{*{20}{c}} {V(L)}\\ {I(L)} \end{array}} \right] = \mathit{\boldsymbol{ \boldsymbol{\varPhi} }}(L)\left[ {\begin{array}{*{20}{c}} {V(0)}\\ {I(0)} \end{array}} \right] + \left[ {\begin{array}{*{20}{c}} {{{\hat V}_s}(L)}\\ {{{\hat I}_s}(L)} \end{array}} \right]$$ (7)

      式中,${\hat V_s}(L)$、${\hat I_s}(L)$分别为平面波在均匀传输线上各点激励起的分布源项在终端L处产生的电压和电流。转移矩阵中各元素的表达式分别为[5]

      $${\phi _{11}} = \cosh (\gamma L) = \frac{{{{\rm{e}}^{\gamma L}} + {{\rm{e}}^{ - \gamma L}}}}{2}$$ (8)
      $${\phi _{12}} = - {Z'_\omega }\sinh (\gamma L) = - {Z'_\omega }\frac{{{{\rm{e}}^{\gamma L}} - {{\rm{e}}^{ - \gamma L}}}}{2}$$ (9)
      $${\phi _{21}} = - \frac{1}{{{{Z'}_\omega }}}\sinh (\gamma L) = - \frac{1}{{{{Z'}_\omega }}}\frac{{{{\rm{e}}^{\gamma L}} - {{\rm{e}}^{ - \gamma L}}}}{2}$$ (10)
      $${\phi _{22}} = \cosh (\gamma L) = \frac{{{{\rm{e}}^{\gamma L}} + {{\rm{e}}^ - }^{\gamma L}}}{2}$$ (11)

      式中,γ为有损大地情况下平面波的传播常数。

    • 用菲涅尔反射系数表示电磁波入射到有损大地上产生的反射激励场,如图 4所示。入射电磁波为平面波,振幅为E0。任意极化方向入射的平面波都可以分解为垂直极化分量(Einc=E0cosα)和水平极化分量(Einc=E0sinα),其中α为极化角,入射电场及反射电场的表达式不再详细推导,具体推导详见文献[4]。

      图  4  有损大地表面平面波入射和反射示意图

      为求得分布电压源和电流源需要先求得电磁场的电磁场分量HyEZ[4]。参照图 4所示,入射场的入射角为φθ,地上部分电场及磁场各分量为入射场和反射场的叠加:

      $$\begin{gathered} E_z^{{\rm{ex}}}(x, z) = E_z^{{\rm{inc}}} + E_z^{{\rm{ref}}}{\rm{ = }} \\ {E_0}\cos \alpha \cos \varphi {{\rm{e}}^{ - {\rm{j}}kx\cos \varphi \sin \theta }}[{{\rm{e}}^{{\rm{j}}kz\sin \varphi }} + {R_v}{{\rm{e}}^{ - {\rm{j}}kz\sin \varphi }}] \\ \end{gathered} $$ (12)
      $$\begin{gathered} E_x^{{\rm{ex}}}(x, 0) = E_x^{{\rm{inc}}} + E_x^{{\rm{ref}}}{\rm{ = }} \\ {E_0}[\cos \alpha \sin \varphi \cos \theta (1 - {R_v}) + \\ {\rm{ }}\sin \alpha \sin \varphi (1 + {R_h})]{{\rm{e}}^{ - {\rm{j}}kx\cos \varphi \cos \theta }} \\ \end{gathered} $$ (13)
      $$\begin{gathered} H_y^{{\rm{ex}}}(x, z) = H_y^{{\rm{inc}}}(x, z) + H_y^{{\rm{ref}}}(x, z) = \\ - \frac{{{E_0}}}{{{{Z'}_\omega }}}[\cos \alpha \cos \theta ({{\rm{e}}^{{\rm{j}}kz\sin \varphi }} + {R_v}{{\rm{e}}^{ - {\rm{j}}kz\sin \varphi }}) + \\ \sin \alpha \sin \varphi \sin \theta ({{\rm{e}}^{{\rm{j}}kz\sin \varphi }} + {R_h}{{\rm{e}}^{ - {\rm{j}}kz\sin \varphi }})]{{\rm{e}}^{ - {\rm{j}}kx\cos \varphi \cos \theta }} \\ \end{gathered} $$ (14)

      利用式(12)~(14)式可以通过求解麦克斯韦方程得到传输线上分布源,进而通过式(7)得到传输线的终端响应值。

      对于有损大地情况,可以将地面的影响考虑成附加了单位地电阻Zg,地电阻的大小用文献[13]的近似表达式得到:

      $${Z_g} = \frac{{{\rm{j}}\omega \mu {}_0}}{{2{\rm{ \mathsf{ π} }}}}\ln \left( {\frac{{1 + {\gamma _g}h}}{{{\gamma _g}h}}} \right)$$ (15)

      式中,γg为电磁波在大地中的传播常数。有损大地的影响是通过影响电磁波的传播常数影响地电阻的大小,从而对传输线的响应产生影响。不同的大地电导率的电阻大小不同。利用该地阻抗可以求得有损大地情况下传输线的特性阻抗和特性导纳。

    • 为了证明上述理论的正确性,本节对接有转换器的非均匀传输线结构进行理论计算,讨论不同大地电导率对理论计算结果的影响。

      非均匀传输线结构如图 5所示,两段特性不同的均匀传输线中间通过转换器进行连接。该结构在连接处传输线特性突变导致传输线整体不均匀,该结构可以等效为三段均匀传输线的级联,从而利用上面推导的公式进行理论计算。其中传输线1和传输线2的长度均为4 m(L1=L2=4 m),半径分别为a1=4 mm,a2=2 mm。传输线距离地面的高度h=0.3 m,平面波采用垂直计划照射(即α=0),入射角度φ=45°,θ=45°,传输线两端电阻R1=0 Ω,R2=50 Ω。

      图  5  传输线与转换器连接示意图

      通过式(8)~式(11)可以得到传输线1和传输线2的网络参数转移矩阵Φ1Φ3,选用的转换器的特性阻抗为Zc2=50 Ω,其转移矩阵为Φ2。则传输线2终端的响应结果由式(7)可以得到:

      $$\left[ {\begin{array}{*{20}{c}} {{V_N}}\\ {{I_N}} \end{array}} \right] = {\mathit{\boldsymbol{ \boldsymbol{\varPhi} }}_1}*{\mathit{\boldsymbol{ \boldsymbol{\varPhi} }}_2}*{\mathit{\boldsymbol{ \boldsymbol{\varPhi} }}_3}\left[ {\begin{array}{*{20}{c}} {V(0)}\\ {I(0)} \end{array}} \right] + \left[ {\begin{array}{*{20}{c}} {{{\hat V}_s}}\\ {{{\hat I}_s}} \end{array}} \right]$$ (16)

      根据上述推导结果,利用MATLAB编程理论计算不同大地电导率情况下传输线2终端响应的频谱结果,如图 6所示。通过图 6可以看出,大地电导率对结果的影响很小,通过图 6b可以看出只在低频端几个特殊频点影响较大,最大相差2.978 dB,因此在实验验证时选取电导率的典型值σg=0.01 S/m;传输线终端感应电流呈现谐振现象,这是由于传输线两端阻抗不匹配反射的结果,图中第一个谐振点对应传输线非均匀处半波长谐振;在高频端,传输线终端感应电流整体大于低频端感应电流。

      图  6  不同大地电导率情况下传输线2终端电流响应结果

    • 为了验证上述数值计算结果的正确性,对接有转接头的传输线的电流响应进行实测验证。实验中用到的实验设备及参数如表 1所示。实验测试平台布置示意图如图 7所示,天线距离线缆的距离为25 m,因此线缆位置属于天线测试各频点的远场区。信号发生器产生特定频率的电磁波通过天线进行发射以此来近似平面波照射。

      表 1  实验设备及参数

      实验设备 型号 技术指标
      信号源 Agilent N9310A 9 kHz~3 GHz
      发射天线 R & S:HK116 20 MHz~300 MHz
      电流探头 BK-CP-02 10 Hz~100 MHz
      接收机 R & S ESCI 9 kHz~3 GHz
      线缆、转接头 HUBER+SUHNER RG_213、32-N-50
      场强仪 TES-593 10 MHz~8 GHz

      图  7  实验测试平台示意图

      首先在没有传输线存在时使用场强仪测量传输线1始端位置处的场强大小,从而可以以此场强值计算得到传输线2终端电流响应的理论值。然后使用电流探头连接到接收机上,对传输线2的终端电流响应进行实测。接收机上得到的读数为电压值,需要加修正因子转换为电流值,同时电缆存在一定的损耗,还要加上电缆的损耗才是实测的电流值。

      选取感应电流极小值和极大值点附近的频率,同时选取几个典型频率值进行实测,根据场强大小,利用微波网络参数转移矩阵理论方法计算得到的电流响应值与实测电流响应值对比见表 2。其中所用线缆为RG_213型号线缆,单位长度的衰减系数为:

      表 2  感应电流理论值与实测值对比结果

      频率/MHz 场强/V*m-1 理论值/dBμA 接收机读数/dBμV 修正因子/dB 线缆损耗/dB 实测值/dBμA 误差/dB
      27.5 0.56 27.46 43.73 -17.0 0.15 26.88 0.58
      30 0.84 44.25 60.28 -17.0 0.15 43.43 0.82
      36.3 0.77 48.17 63.68 -17.0 0.17 46.83 1.34
      43 0.68 36.91 53.07 -17.0 0.19 36.24 0.67
      48 0.74 29.84 47.48 -17.0 0.20 30.63 0.79
      60 0.61 34.88 50.42 -17.0 0.22 33.64 1.24
      75 0.68 44.54 62.01 -16.9 0.25 45.26 0.72
      90 0.66 39.43 53.87 -16.8 0.28 37.45 1.98
      115 0.73 48.87 52.97 -16.5 0.32 46.62 2.25
      $${\rm{Att = }}a\sqrt f + bf$$ (17)

      式中,a=0.167 9;b=0.058 5;f 的单位为GHz;Att单位为dB/m。其中所用线缆的长度为5 m,由此可以得到不同频率下线缆的损耗值。其中,电流实测值=接收机电压读数+修正因子+线缆损耗。表 2给出了感应电流实测值与理论计算值的对比结果。

      通过上表可以看出,实测结果与理论值的误差范围为0.58~2.25 dB。其中造成误差的原因有以下3点:1)大地电导率的不均匀性会造成地电阻的不均匀,从而影响感应电场的大小;2)线缆接头处接触不好也会对测试结果造成一定影响;3)天线产生的电磁场在远场处可以近似看作平面波,此处的近似也会对测试结果产生一定的影响。同时接近100 MHz的两个频点90 MHz和115 MHz实测值和理论值相差较大,这是由于电流探头工作频带范围为10 Hz~100 MHz,超出电流探头工作频段上限会造成较大误差。

      综上所述,通过表 2可以看出通过微波网络参数转移矩阵理论计算得到的理论值与实测值之间误差在3 dB以内,证明了使用微波网络参数方法计算突变型非均匀传输线的正确性。

    • 本文将微波网络中的转移矩阵应用到传输线模型中,对有损大地上突变型非均匀传输线进行数值计算,并通过实验进行了验证,实验结果与理论计算结果吻合较好。该方法特别适合外场激励下有明显突变的非均匀传输线终端响应的计算,省去了传输线中间非均匀段的繁琐计算,只要知道传输线一端的响应参数,另一端的响应参数可以很容易地得到。为计算接有转接器等类型非均匀传输线提供了理论方法。

参考文献 (28)

目录

    /

    返回文章
    返回