一种行波管注波互作用分布结构的优化方法技术

该发明专利技术公开了一种行波管注波互作用分布结构的优化方法，目的是为了提高行波管注波互作用优化算法的优化效率。该发明专利技术将优化算法与注波互作用内在机理结合，涉及嵌套算法，利用注波互作用过程顺序计算特性以及相对相位角Φ是否满足或的条件提前结束计算。对于普通的高效率螺距分布结构，该发明专利技术相对于普通优化算法可以提高效率3-4倍。

【技术实现步骤摘要】
一种行波管注波互作用分布结构的优化方法
本专利技术属于行波管
，具体涉及一种行波管注波互作用分布结构优化方法。
技术介绍
行波管是使用最广泛的真空电子器件之一，广泛应用于卫星通讯、雷达、电子对抗等领域。目前，采用计算机辅助设计(ComputerAidedDesign,CAD)技术是节约成本，改进设计和提高行波管整体性能的主要手段之一。在行波管CAD技术中，对行波管内电子注与高频电磁场的相互作用（即注波互作用）进行分析具有重要意义。行波管中的注波互作用是一个信号从左往右能量放大的过程：信号在行波管高频结构内的左端激励高频电磁场，然后在一定长度的互作用区域内与被约束的电子注相互作用（注波互作用过程），期间高频电磁场的能量被逐步放大，最终在右端通过输出窗输出高频电磁场的能量，得到放大的信号。并且通过设计合理的注波互作用分布结构，控制注波互作用过程，能够优化最终输出的功率、电子效率等行波管性能参数。在行波管设计过程中目前主要采用半解析半数值的参数理论模型（参见“一种行波管注波互作用的模拟方法”专利CN201110236508.8，2011，胡玉禄，彭维峰，李斌等）进行注波互作用过程计算，单次计算时间大约几十秒钟乃至几分钟左右，但是复杂的行波管注波互作用分布结构设计（见图1）往往需要采用优化算法大量繁琐的迭代进行优化计算，普通的优化算法视一次注波互作用过程计算为整体（见图2），多达数百上千次的迭代次数乘以单次注波互作用计算时间，也将带来数十分钟至数小时的计算时间。然而，行波管的注波互作用计算本身具有顺序计算特性，互作用过程往往从左往右顺序计算（见图1），右边的结构变化带来的性能改变将不影响左边的计算结果。同时可以利用注波互作用过程中相对相位角Φ是否满足或的条件判断是否性能开始下降，从而提前结束互作用计算。
技术实现思路
本专利技术的目的是为了提高行波管注波互作用优化算法的优化效率。基于此目的，本专利技术提出了一种行波管注波互作用分布结构的优化方法。本专利技术的核心是将优化算法与注波互作用内在机理结合，利用注波互作用过程顺序计算特性以及相对相位角Φ是否满足或的条件提前结束计算，实现行波管注波互作用高效率优化的目的。为了实现上述目的，本专利技术的技术方案是：一种行波管注波互作用分布结构高效率优化方法。涉及嵌套算法，包括以下步骤：A、计算注波互作用至第一个备份点，备份电子注所有粒子的运动状态以及高频电磁场状态数据；B、从备份点恢复电子注所有粒子的运动状态以及高频电磁场状态数据，并计算后续注波互作用；C、计算期间判断相对相位角Φ是否满足或的条件；D、当前计算停止，分为三种情况：①相对相位角Φ满足或②达到注波互作用输出端、③计算至下一备份点进行嵌套优化；E、在①和②的情况下，计算停止并返回计算结果；在③的情况下，备份数据，重复B、C、D，进行嵌套优化并返回优化结果；F、反复调用步骤B、C、D、E，根据计算输出结果，利用优化算法策略，反馈优化调整输入参数，直至达到优化目标要求，优化计算结束。步骤A中，备份点可以设置两个或两个以上，其中第一个备份点之前的结构是不需要优化。本专利技术的有益效果：利用本专利技术提出的波管注波互作用分布结构优化方法优化注波互作用分布，优化效率视备份点位置而定，对于普通的高效率螺距分布结构（如图1），本专利技术相对于普通优化算法可以提高效率3-4倍。附图说明图1是行波管注波互作用区域结构示意图；图2是行波管注波互作用螺距分布结构示意图；图3是普通优化算法实施过程示意图；图4是本专利技术实施例优化算法实施过程示意图；图5是注波互作用从左往右求解过程示意图。附图标号说明：1为电子注，2为输入信号，3为集中衰减器，4为输出信号，5为夹持杆，6为管壳。具体实施方式下面以图1互作用分布结构优化为例，结合附图和具体实施例子对本专利技术作进一步说明。具体包含6个变量（即四个分布长度变量z1、z2、z3、z4和两个螺距变量p1、p2，其中p1为分布长度z2部分的螺距大小，p2为分布长度z4部分的螺距大小）的优化，p0（分布长度z1部分的螺距大小）变量通过经验或者适当的扫描确定，可以节约优化时间。将这6个互作用分布结构的变量分为区域1和区域2两部分进行优化，其中区域2嵌套调用区域1内变量的优化算法，将其优化的最佳结果返回给区域1内的优化算法。1、计算注波互作用至备份点A，备份电子注所有粒子的运动状态以及高频电磁场状态数据。行波管中的注波互作用计算采用如图4所示的流程实现从左往右顺序计算，在信号进入的输入端初始化电子注内所有粒子的位置、速度等运动状态以及输入信号决定的高频电磁场的场幅值以及相位状态，利用粒子运动的微分方程以及高频电磁场微分方程（见公式1.1以及公式1.2，数值计算上采用差分方程），求解下一个位置的电子注运动状态以及高频电磁场状态，重复上述过程逐步计算求解，最终可以获得输出口的电子注运动状态以及高频电磁场状态。根据输出口的电磁场状态获得信号放大以后的功率。典型的电子运动微分方程：公式（1.1）m0表示电子质量，e表示电子的电荷，c表示光速，γν表示电子v的相对论因子，为本征功率。a(z)为高频场幅值，是轴向位置变量z的复函数,矢量x表示空间位置矢量，下标“⊥”表示横向分量，变量z和t分别表示轴向位置和时间。e(x⊥,z)为无源的高频电场的分布函数，满足无源麦克斯韦方程组。<·>beam表示所有电子的平均值，ψv表示电子v的相位，i为虚数因子，ω＝2πf为角频率，f表示本征频率，Ib表示电流大小，rbi和rbo分别表示电子注内外半径，通常rbi为零表示实心电子注。R表示等离子体降低因子。下标m表示空间谐波次数。典型的高频电磁场微分方程：公式（1.2）其中，是关于轴向位置z的微分运算。σ(z)为单位长度高频场的衰减量，即衰减系数；l为高频结构的轴向周期长度，T＝2π/ω为与信号频率相对应的时间周期，表示沿高频结构横截面积分。j(x⊥,z)为电流密度分布，上标“*”表示取共轭。2、优化算法确定注波互作用分布结构变量z1，p1和z2初始值，从备份点A恢复电子注所有粒子的运动状态以及高频电磁场状态，为后续互作用计算提供起始计算数据。3、与步骤1类似，利用粒子运动的微分方程以及高频电磁场微分方程,计算后续注波互作用至备份点B（典型位置为互作用长度的75%-80%处）。计算过程中如果相对相位角Φ满足或的条件，计算停止，将当前计算结果作为优化目标返回给优化算法程序进行优化调整。如果相对相位角没有满足上述条件，计算至备份点B，备份当前电子注运动状态以及高频电磁场状态数据，嵌套优化变量注波互作用分布结构变量z3，p2和z4。将优化后的计算结果作为优化目标返回给优化算法程序进行优化调整。计算结果可以是单个频率的输出功率、输出电子效率，从而优化单频率的输出功率、输出电子效率最大，也可以是带宽内多个频率的输出功率、输出电子效率的最小值，从而优化整个带宽内输出功率、电子效率最大。4、在计算过程中，加入相对相位角Φ是否满足或的条件判断是否提前结束计算。场相对相位角Φ定义为场F与单位距离场增量之间的相位差，即：公式（1.3）公式（1.3）中，场相对相位角Φ表征了场与电子注相互作用的同步状态，反映了场幅值的增减变化。当Φ近似等于零时，场增加最快；当或本文档来自技高网...

【技术保护点】
一种行波管注波互作用分布结构高效率优化方法，其特征在于，涉及嵌套算法，包括以下步骤：A、计算注波互作用至第一个备份点，备份电子注所有粒子的运动状态以及高频电磁场状态数据；B、从备份点恢复电子注所有粒子的运动状态以及高频电磁场状态数据，并计算后续注波互作用；C、计算期间判断相对相位角Φ是否满足或的条件；D、当前计算停止，分为三种情况：①相对相位角Φ满足或②达到注波互作用输出端、③计算至下一备份点进行嵌套优化；E、在①和②的情况下，计算停止并返回计算结果；在③的情况下，备份数据，重复B、C、D，进行嵌套优化并返回优化结果；F、反复调用步骤B、C、D、E，根据计算输出结果，利用优化算法策略，反馈优化调整输入参数，直至达到优化目标要求，优化计算结束。

【技术特征摘要】
1.一种行波管注波互作用分布结构高效率优化方法，其特征在于：涉及嵌套算法，包括以下步骤：A、计算注波互作用至第一个备份点，备份电子注所有粒子的运动状态以及高频电磁场状态数据；B、从备份点恢复电子注所有粒子的运动状态以及高频电磁场状态数据，并计算后续注波互作用；C、计算期间判断相对相位角Φ是否满足或的条件；D、当前计算停止，分为三种情况：①相对相位角Φ满足或②达到注波互作用输出端、③计算至计算至第N+1备份点进行第N层嵌套优化；E、在①和②的情况下，计算停止并返回计算结果；在③的情况下，备份数据，重复B、C、D，进行嵌套优化并返回优化结果；F、反复调用步骤B、C、D、E，根据计算输出结果，利用优化算法策略，反馈优化调整输入参数，直至达到优化目标要求，优化计算结束。2.如权利要求1所述的一种行波管注波互作用分布结构高效率优化方法，其特征在于，备份点可以设置两个或两个以上，其中第一个备份点之前的结构不需要优化。3.如权利要求1所述的一种行波管注波互作用分布结构高效率优化方法，其特征在于，步骤D中当N＝1，即第二个备份点的位置为互作用长度的75％-80％处。4.如权利要求1所述的一种行波管注波互作用分布结构高效率优化方法，其特征在于，在信号的输入端初始化粒子的运动状态以及高频电磁场的场幅值以及相位状态，利用粒子运动的微分方程以及高频电磁场微分方程，求解下一个位置的电子注运动状态以及高频电磁场状态，直至获得输出口的电子注运动状态以及高频电磁场状态；所述电子运动微分方程为：

【专利技术属性】
技术研发人员：胡玉禄，李斌，李建清，杨中海，
申请(专利权)人：电子科技大学，
类型：发明
国别省市：四川;51