CN103474737A - 支持向量机对膜片建模的毫米波e面滤波器及膜片建模方法 - Google Patents

支持向量机对膜片建模的毫米波e面滤波器及膜片建模方法 Download PDF

Info

Publication number
CN103474737A
CN103474737A CN2013103649654A CN201310364965A CN103474737A CN 103474737 A CN103474737 A CN 103474737A CN 2013103649654 A CN2013103649654 A CN 2013103649654A CN 201310364965 A CN201310364965 A CN 201310364965A CN 103474737 A CN103474737 A CN 103474737A
Authority
CN
China
Prior art keywords
diaphragm
filter
rectangular waveguide
svms
face
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.)
Granted
Application number
CN2013103649654A
Other languages
English (en)
Other versions
CN103474737B (zh
Inventor
孙璐
詹劲松
文枫
王家礼
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Xidian University
Original Assignee
Xidian University
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Xidian University filed Critical Xidian University
Priority to CN201310364965.4A priority Critical patent/CN103474737B/zh
Publication of CN103474737A publication Critical patent/CN103474737A/zh
Application granted granted Critical
Publication of CN103474737B publication Critical patent/CN103474737B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Prostheses (AREA)
  • Control Of Motors That Do Not Use Commutators (AREA)

Abstract

本发明涉及毫米波E面滤波器的设计方法,特别是支持向量机对膜片建模的毫米波E面滤波器及膜片建模方法。该支持向量机对膜片建模的毫米波E面滤波器,其特征是,至少包括:矩形波导(1),两个或两个以上的金属膜片(2),矩形波导(1)中间右E面腔;E面腔与矩形波导长度方向平行,其中包含金属膜片(2);所述金属膜片(2)位于垂直矩形波导E面中央,并与E面平行,相邻金属膜片(2)不等间隔分布,相邻金属膜片在沿矩形波导长度方向一条直线,膜片面与矩形波导长度方向一条直线垂直。本发明大大的提高了滤波器的设计速度和精度,避免了重复使用电磁计算的复杂方法。

Description

支持向量机对膜片建模的毫米波E面滤波器及膜片建模方法
技术领域
本发明涉及毫米波E面滤波器的设计方法,特别是支持向量机对膜片建模的毫米波E面滤波器及膜片建模方法。
背景技术
在毫米波E面滤波器的设计中,对E面膜片的建模和分析是一个关键的环节。以往采用的数值分析方法虽然精度较高,但建模速度慢。
发明内容
本发明的目的是提供一种建模精度高,建模速度快的支持向量机对膜片建模的毫米波E面滤波器及膜片建模方法,以便大大提高滤波器的设计速度。
本发明的目的是这样实现的,支持向量机对膜片建模的毫米波E面滤波器,其特征是,至少包括:矩形波导,两个或两个以上的金属膜片,矩形波导中间右E面腔;E面腔与矩形波导长度方向平行,其中包含金属膜片;所述金属膜片位于垂直矩形波导E面中央,并与E面平行,相邻金属膜片不等间隔分布,相邻金属膜片在沿矩形波导长度方向一条直线,膜片面与矩形波导长度方向一条直线垂直。
所述金属膜片起耦合作用,相邻膜片之间构成谐振腔。
支持向量机对膜片建模的毫米波E面滤波器的膜片建模方法,其特征是:至少包括如下步骤:
步骤600:根据提出的指标综合得到滤波器的低通原型,计算出K变换器的归一化阻抗,代入滤波器设计公式,以计算得到各个并联阻抗值;
步骤601:调用已经建立好的支持向量机模型,送入滤波器设计频率和多个金属膜片的长度,可以快速仿真得到对应的并联阻抗值;
步骤602:将两组并联阻抗值经过对比,找出最接近最优的膜片长度和谐振腔长度;
步骤603:确定好金属膜片的长度和其在波导中的位置,以确定E面膜片滤波器的基本结构;
步骤604:利用高频电磁仿真软件验证滤波器的设计结果,并优化设计;
步骤605:对滤波器进行试用和调试。
所述的步骤601至少包括如下步骤:
步骤500:选定矩形波导的型号,并确定E面腔中每个金属膜片的厚度t和中心频率f0以及膜片在矩形波导中所处的位置;
步骤501:以金属膜片的长度w为变量,以模式匹配法为基础求得相应的S参数,其中S包括S11,S12,S21,S22
步骤502:使用场参数和电路参数之间的变换公式,求得对应T-型等效电路的参数Xs和Xp;所用公式如下所示:
jX s = 1 - s 12 + s 11 1 - s 11 + s 12 jX p = 2 s 12 ( 1 - s 11 ) 2 - s 12 2 - - - ( 1 )
则相应的K阻抗变换器为:
φ=-tan-1(2Xp+Xs)-tan-1XsK=|tan(φ/2+tan-1Xs)|  (2)
其中,Xp和Xs表示T型等效电路电抗值;φ表示等效电路中的虚拟电长度;K为K变换器的变换比;S表示金属膜片的散射参数;j表示电抗元件值的虚部单位。
步骤503:采用步骤501和502的方法,计算出几组频率、金属膜片的长度w和等效电路参数之间的对应关系数据,将这些数据分两部分,以作为支持向量机的样本参数和实验数据;
步骤504:对步骤503所获取的样本数据进行归一化处理,选择支持向量机结构的核函数;
步骤505:将样本数据送到支持向量机,调整参数,建立精确的支持向量机模型;
步骤506:对支持向量机模型送入验证数据,验证支持向量机模型的精确度;
步骤507:结束支持向量机对E波膜片的建模过程。
本发明的优点是:通过支持向量机建立起来的膜片模型,在波导工作频带内的任意频率,都可以快速得到膜片长度w和T型等效电路参数Xs和Xp之间的对应关系。在设计滤波器时,首先根据提出的指标综合出滤波器的低通原型,计算出K变换器的归一化阻抗,带入公式计算出各个并联阻抗值,通过支持向量机模型结果经过搜索就可以得到与各个并联阻抗值对应的膜片长度w。确定了耦合段的大小后,谐振腔长度也就确定了,从而可以得到滤波器的原始尺寸,大大的提高了滤波器的设计速度和精度,避免了重复使用电磁计算的复杂方法。
附图说明
下面结合具体实施例附图对本发明作进一步说明:
图1是矩形波导E面膜片滤波器的结构示意图;
图2是图1的参数示意图;
图3是金属膜片用所示的T型等效电路;
图4支持向量机建模过程。
图5支持向量机对E面膜片建模流程图;
图6支持向量机对膜片建模的毫米波E面滤波器设计流程图。
图中,1、矩形波导;2、金属膜片。
具体实施方式
如图1、图2所示,支持向量机对膜片建模的毫米波E面滤波器,其特征是,至少包括:矩形波导1,两个或两个以上的金属膜片2,矩形波导1中间右E面腔;E面腔与矩形波导长度方向平行,其中包含金属膜片2;所述金属膜片2位于垂直矩形波导E面中央,并与E面平行,相邻金属膜片2不等间隔分布,相邻金属膜片在沿矩形波导长度方向一条直线,膜片面与矩形波导长度方向一条直线垂直。
所述金属膜片2起耦合作用,相邻膜片之间构成谐振腔。
支持向量机对膜片建模的毫米波E面滤波器设计流程如图6所示,其特征是,至少包括如下步骤:
步骤600:根据提出的指标综合得到滤波器的低通原型,计算出K变换器的归一化阻抗,代入滤波器设计公式,以计算得到各个并联阻抗值;
步骤601:调用已经建立好的支持向量机模型,送入滤波器设计频率和多个金属膜片2的长度,可以快速仿真得到对应的并联阻抗值;
步骤602:将两组并联阻抗值经过对比,找出最接近最优的膜片长度和谐振腔长度;
步骤603:确定好金属膜片2的长度和其在波导中的位置,以确定E面膜片滤波器的基本结构;
步骤604:利用高频电磁仿真软件验证滤波器的设计结果,并优化设计;
步骤605:对滤波器进行试用和调试。
如图3所示,在本发明中,支持向量机对E面膜片建模是一个非常关键的步骤,其基本过程如图5所示,所述的步骤601包括至少包括如下步骤:
步骤500:选定矩形波导1的型号,并确定E面腔中每个金属膜片2的厚度t和中心频率f0以及膜片在矩形波导中所处的位置;
步骤501:以金属膜片2的长度w为变量,以模式匹配法为基础求得相应的S(包括S11,S12,S21,S22)散射参数;
步骤502:使用场参数和电路参数之间的变换公式,求得对应T-型等效电路的参数Xs和Xp;所用公式如下所示:
jX s = 1 - s 12 + s 11 1 - s 11 + s 12 jX p = 2 s 12 ( 1 - s 11 ) 2 - s 12 2 - - - ( 1 )
则相应的K阻抗变换器为:
φ=-tan-1(2Xp+Xs)-tan-1XsK=|tan(φ/2+tan-1Xs)|   2)
其中,Xp和Xs表示T型等效电路电抗值;φ表示等效电路中的虚拟电长度;K为K变换器的变换比;S(包括S11,S12,S21,S22)表示金属膜片2的散射参数;j表示电抗元件值的虚部单位。
步骤503:采用步骤501和502的方法,计算出几组频率、金属膜片的长度w和等效电路参数之间的对应关系数据,将这些数据分两部分,以作为支持向量机的样本参数和实验数据;
步骤504:对步骤503所获取的样本数据进行归一化处理,选择支持向量机结构的核函数;
步骤505:将样本数据送到支持向量机,调整参数,建立精确的支持向量机模型;
步骤506:对支持向量机模型送入验证数据,验证支持向量机模型的精确度;
步骤507:结束支持向量机对E波膜片的建模过程。
支持向量机的基本思想是通过用内积函数定义的非线性变换将输入空间变换到一个高维空间,在这个高维空间中寻找输入变量和输出变量之间的一种非线性关系,其基本结构如图4所示。支持向量机算法是一个凸二次优化问题,保证找到的解是全局最优解,能较好地解决小样本、非线性和高维数、局部极小点等实际问题。
将已经采用模式匹配法计算出的几组频率值和膜片长度w作为输入样本,输入支持向量机结构。首先对样本数据做归一化,将所有建模数据归一化到[-1,1]范围内,最后将建模预测结果再进行反归一处理即可。归一化的主要优点是有效的限制数据变化范围,提高建模效率,同时还可以避免大量计算的困难。由于E面膜片模型的建模关系是一组非线性映射关系,因此选取了ε-SVR型支持向量机完成函数拟合。
核函数的选择是支持向量机理论研究的一个核心问题,但是目前还没有一种针对具体问题构造出合适的核函数的有效方法。本发明选择的核函数为径向基核函数,原因是由于①RBF核把样本非线性映射到更高维的空间上,因而可以处理类标签和属性之间的非线性关系。线性核是RBF核的特例,含有惩罚参数C的线性核同有参数(C,γ)(其中C为惩罚因子,γ为核参数)的RBF核具有相同的性能。另外,对于一定的参数,sigmoid核的性能与RBF核差不多。②超参数的个数影响模型选择的复杂性。多项式核的超参数比RBF核的要多,选择起来比较复杂。③RBF核的计算难度比较小。如果多项式的次数比较高的话,在运算时可能会出现溢出的现象。
在使用径向基核函数时,有两个参数需要调整,分别是C和γ,事先并不能知道对某一问题哪一个C和γ是合适的,因此需要对这两个参数进行合适的选择,目标是能利用支持向量机最小误差的预测未知的非线性散射函数数据。
通常的方式是将训练数据分成两个部分,一部分作为训练数据来训练支持向量机,调整C和γ,获得最好的训练精度,使得模拟结果和训练数据吻合良好;另一部分作为验证数据,预测出未知频率和功率点的模型输出,并与实际测量值做比较,得出预测的精确度高低指标。
如对于BJ320波导,对三个频率样点即f=33.0GHz、35.0GHz、37.0GHz,采用模式匹配方法计算出膜片长度与T-型等效电路参数值作为输出样本集。确定其中35.0GHz和37.0GHz数据作为样本点来建立支持向量机模型,33.0GHz数据作为网络仿真验证误差使用。可以看出建立起的支持向量机模型很好的预测了33.0GHz频率点的等效电路参数。
本实施例没有详细叙述的部分、常用结构及英文缩写属本行业的公知常识,这里不一一叙述。

Claims (4)

1.支持向量机对膜片建模的毫米波E面滤波器,其特征是,至少包括:矩形波导(1),两个或两个以上的金属膜片(2),矩形波导(1)中间右E面腔;E面腔与矩形波导长度方向平行,其中包含金属膜片(2);所述金属膜片(2)位于垂直矩形波导E面中央,并与E面平行,相邻金属膜片(2)不等间隔分布,相邻金属膜片在沿矩形波导长度方向一条直线,膜片面与矩形波导长度方向一条直线垂直。
2.依据权利要求书1所述支持向量机对膜片建模的毫米波E面滤波器,其特征是:所述金属膜片(2)起耦合作用,相邻膜片之间构成谐振腔。
3.支持向量机对膜片建模的毫米波E面滤波器的膜片建模方法,其特征是:至少包括如下步骤:
步骤600:根据提出的指标综合得到滤波器的低通原型,计算出K变换器的归一化阻抗,代入滤波器设计公式,以计算得到各个并联阻抗值;
步骤601:调用已经建立好的支持向量机模型,送入滤波器设计频率和多个金属膜片(2)的长度,可以快速仿真得到对应的并联阻抗值;
步骤602:将两组并联阻抗值经过对比,找出最接近最优的膜片长度和谐振腔长度;
步骤603:确定好金属膜片(2)的长度和其在波导中的位置,以确定E面膜片滤波器的基本结构;
步骤604:利用高频电磁仿真软件验证滤波器的设计结果,并优化设计;
步骤605:对滤波器进行试用和调试。
4.根据权利要求3所述的支持向量机对膜片建模的毫米波E面滤波器的膜片建模方法,其特征是,所述的步骤601至少包括如下步骤:
步骤500:选定矩形波导(1)的型号,并确定E面腔中每个金属膜片(2)的厚度t和中心频率f0以及膜片在矩形波导中所处的位置;
步骤501:以金属膜片(2)的长度w为变量,以模式匹配法为基础求得相应的S参数,其中S包括S11,S12,S21,S22
步骤502:使用场参数和电路参数之间的变换公式,求得对应T-型等效电路的参数Xs和Xp;所用公式如下所示:
jX s = 1 - s 12 + s 11 1 - s 11 + s 12 jX p = 2 s 12 ( 1 - s 11 ) 2 - s 12 2 - - - ( 1 )
则相应的K阻抗变换器为:
φ=-tan-1(2Xp+Xs)-tan-1XsK=|tan(φ/2+tan-1Xs)|  (2)
其中,Xp和Xs表示T型等效电路电抗值;φ表示等效电路中的虚拟电长度;K为K变换器的变换比;S表示金属膜片(2)的散射参数;j表示电抗元件值的虚部单位。
步骤503:采用步骤501和502的方法,计算出几组频率、金属膜片的长度w和等效电路参数之间的对应关系数据,将这些数据分两部分,以作为支持向量机的样本参数和实验数据;
步骤504:对步骤503所获取的样本数据进行归一化处理,选择支持向量机结构的核函数;
步骤505:将样本数据送到支持向量机,调整参数,建立精确的支持向量机模型;
步骤506:对支持向量机模型送入验证数据,验证支持向量机模型的精确度;
步骤507:结束支持向量机对E波膜片的建模过程。
CN201310364965.4A 2013-08-20 2013-08-20 支持向量机对膜片建模的毫米波e面滤波器及膜片建模方法 Expired - Fee Related CN103474737B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201310364965.4A CN103474737B (zh) 2013-08-20 2013-08-20 支持向量机对膜片建模的毫米波e面滤波器及膜片建模方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201310364965.4A CN103474737B (zh) 2013-08-20 2013-08-20 支持向量机对膜片建模的毫米波e面滤波器及膜片建模方法

Publications (2)

Publication Number Publication Date
CN103474737A true CN103474737A (zh) 2013-12-25
CN103474737B CN103474737B (zh) 2015-09-09

Family

ID=49799489

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201310364965.4A Expired - Fee Related CN103474737B (zh) 2013-08-20 2013-08-20 支持向量机对膜片建模的毫米波e面滤波器及膜片建模方法

Country Status (1)

Country Link
CN (1) CN103474737B (zh)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104241746A (zh) * 2014-09-09 2014-12-24 江苏贝孚德通讯科技股份有限公司 一种波导高频低通滤波器
CN108963403A (zh) * 2018-06-26 2018-12-07 东南大学 一种基于波导电感性窗的Doherty功率合成器
CN110543703A (zh) * 2019-08-19 2019-12-06 华南理工大学 一种考虑不同时间尺度的准谐振变换器建模分析方法

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20030174032A1 (en) * 2002-03-14 2003-09-18 Noyan Kinayman Surface mountable microwave filter configuration and method of fabricating same

Patent Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20030174032A1 (en) * 2002-03-14 2003-09-18 Noyan Kinayman Surface mountable microwave filter configuration and method of fabricating same

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
刘斌等: "一种基于支持向量机建模的新型自适应滤波器", 《大庆石油学院学报》 *
刘渝: "波导E面金属膜片的分析及其在滤波器设计中的应用", 《电子科技》 *

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104241746A (zh) * 2014-09-09 2014-12-24 江苏贝孚德通讯科技股份有限公司 一种波导高频低通滤波器
CN108963403A (zh) * 2018-06-26 2018-12-07 东南大学 一种基于波导电感性窗的Doherty功率合成器
CN110543703A (zh) * 2019-08-19 2019-12-06 华南理工大学 一种考虑不同时间尺度的准谐振变换器建模分析方法
CN110543703B (zh) * 2019-08-19 2021-05-14 华南理工大学 一种考虑不同时间尺度的准谐振变换器建模分析方法

Also Published As

Publication number Publication date
CN103474737B (zh) 2015-09-09

Similar Documents

Publication Publication Date Title
CN108170922B (zh) 一种微波滤波器的辅助调试方法、设备及存储设备
Borries et al. Multilevel fast multipole method for higher order discretizations
CN102930071B (zh) 基于非匹配网格的周期结构的三维电磁场仿真模拟方法
CN104794541B (zh) 一种基于模拟退火和锥优化的配电网运行优化方法
CN113158527B (zh) 一种基于隐式fvfd计算频域电磁场的方法
Lamecki et al. An efficient framework for fast computer aided design of microwave circuits based on the higher-order 3D finite-element method
CN112487755B (zh) 一种fltd腔体中瞬态电磁场分布的数值计算方法
CN101478069A (zh) 基于核机器学习的微波滤波器辅助调试方法
CN103474737B (zh) 支持向量机对膜片建模的毫米波e面滤波器及膜片建模方法
Pietrenko-Dabrowska et al. Fast design closure of compact microwave components by means of feature-based metamodels
CN106294894B (zh) 快速分析非均匀目标电磁散射特性的有限元边界积分方法
Flisgen et al. Compact time-domain models of complex RF structures based on the real eigenmodes of segments
CN103474989A (zh) 一种基于灵敏度分析的网络重构方法
CN106786529A (zh) 一种分布式静态安全分析方法
CN101930495A (zh) 基于网格转换的平板裂缝阵天线机电综合分析方法
CN104156504A (zh) 一种发电机励磁系统参数可辨识性判定方法
CN104809297A (zh) 一种磁场-结构场耦合计算中的异型网格间的电磁力密度传递方法
Wu et al. Parametric model for microwave filter by using multiple hidden layer output matrix extreme learning machine
CN105205299B (zh) 电大目标电磁散射特性快速降维分析方法
CN105868452B (zh) 一种基于等几何分析的多孔板力学性能分析方法
CN107526869A (zh) 一种基于函数逼近自适应三维微波管输入输出窗模型降阶的数值方法
CN105356068B (zh) 一种宽带馈源的优化设计方法
Jolly et al. Novel integrated design framework for radio frequency quadrupoles
Guo et al. Recent Advances in Gradient-based EM Optimization using Neuro-TF Surrogate
CN109670254A (zh) 机电暂态与电磁暂态混合仿真的接口位置选择方法及系统

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20150909

Termination date: 20200820

CF01 Termination of patent right due to non-payment of annual fee