matlab atem(),一种直接计算ATEM传输室谐振频率的简单方法

横电磁波传输室(TEM)是近20年发展起来的一项新技术,已应用于高频微波电磁近区场可计算的计量标准.它的使用范围受最高使用上限频率(第一谐振频率)和室内场均匀性的限制,尺寸越大,最高使用频率越低,为此,应用时一般采用加填吸波材料的方法提高上限使用频率.同时,为了扩大有效测试空间以便进行设备的电磁干扰(EMI)和电磁抗扰度(EMS)测量,一般采用非对称的横电磁波传输室(ATEM).ATEM的传输主模是TEM模,然而由于角锥段两端渐近缩小,形成近似谐振腔,高阶模在其中发生谐振,这样,它的最高使用频率就受到第一谐振频率的限制.而目前的研究工作[114]无法准确计算谐振频率,只能通过事后的测量手段或经验公式计算获得,这给设计制造带来诸多不便.目前国内已有人采用时域有限差分(FDTD)方法对此进行了研究,但过于繁杂,且计算结果不理想.因此,有必要进行深入研究,寻找可实用的计算谐振频率的方法.1等效传输线模型典型的非对称横电磁波室如图1、2所示.其中L为ATEM室的长度,单位m;w为内导体板的半宽度,单位m.它主要包括2部分:主段和图1ATEM室的俯视图图2ATEM室的横截面过渡段.谐振发生时,中垂面MN对于横电场(TE2nmp)或横磁场(TM(2n+1)mp)模式相当于电墙,等效于传输线的短路状态.而对于TE(2n+1)mp或TM2nmp模式则相当于磁墙,等效于传输线的开路状态.这样,沿平行于MN且距离为d(m)的任一平面AA与ATEM室相截所得截面如图3所示,沿x方向可等效为一非均匀传输线,如图4所示.图3沿x方向截面图4沿x方向等效非均匀传输线通过几何计算,对于1.0m0.6m2.0m且w/a=0.36的ATEM室,h=0.6L,B=0.583(1-2d),Lc=2.0-1.08(L/w),单位m.这样,就可采用矩量法(点选配)计算图3所表示的结构电容,而通过Z=00C=0C(1)计算图4所标示的等效非均匀传输线的沿线阻抗分布.式(1)中,Z为传输线特性阻抗,单位;C为分布电容,单位F;0和0为自由空间的介电常数和磁导率;为自由空间的波阻抗,大小约为377.将图4非均匀传输线等分多段,通过传输线转移阻抗计算公式(2)计算终端输入阻抗,进而,联立式(3)和(2)计算其谐振频率.Zin=ZcZL+jZctg(kd)Zc+jZLtg(kd)(2)其中,k为波数;Zc为终端输入阻抗,单位;d为AA面距中垂面MN的距离,单位为m.Zg=-j2fCg(3)其中,f为频率,单位MHz;Cg为终端气隙电容(F/m).气隙电容的确定采用文献[15]中图4的奇模电容Cf0,在内导体板厚度t=0,且w/a=0.36的条件下查得单位长气隙电容为1.60(F/m),那么沿线气隙电容Cg为1.470(F/m).2数值计算结果采用矩量法(点选配)计算图3所表示的结构电容并计算沿线阻抗,计算结果如图5所示.利用式(2)、(3)及图5结果,采用Matlab编制程序,计算的谐振频率与文献[13]的对比如表1所示.表1谐振频率计算结果fR/MHzTE(101)TM(110)TE(101)TE(111)TE014本文结果151.5205.5186.5292.5346.5文献[13]1563-1833--文献[15]156.3-184.3272.7351.3本文的计算方法所得的结果与文献吻合较好,实现了谐振频率的直接计算.克服了以往理论计算谐振频率需依靠实验测定的经验系数的方法.图5X方向截面阻抗分布3ATEM室的最高使用上限频率的结论由于ATEM是两端渐进缩小,形成近似的谐振腔,电磁波在倾斜