本发明涉及一种hpm产生器件优化方法。
背景技术：
随着hpm(highpowermicrowave，高功率微波)产生器件的发展，其结构越来越复杂。为进一步提高器件功率、加快起振、缩小频率变化范围，需要对器件进行优化。
采用优化算法，能提高器件优化速度与优化质量。文献1[1]公开了采用遗传算法对hpm产生器件进行优化的方法，但是文献1的方法仅能优化hpm产生器件的输出功率。而hpm产生器件的其它指标，如产生信号频率，起振时间，提取腔反射器内最大场强都是影响器件性能的重要参数，任何一项参数不符合设计要求，都有可能导致优化结果无效。因此，出现了对hpm产生器件进行多目标优化的方法。
传统的对hpm产生器件进行多目标优化方法为加权平均的多目标优化方法，是通过下面的公式(1)进行多目标优化。这种方法将多目标参数通过一定的权重，转化为单目标优化。这种多目标优化方法的缺点在于加权平均系数难以确定，不同的权重系数，对优化结果影响很大，优化结果具有很大的不确定性。
其中ωi为第i个适应性值的权重，i(x)为第i个优化目标的适应性函数，x为需要优化的向量。
[1]文献1，zaigaochen,jianguowang,andyuewang,optimizationofrelativisticbackwardwaveoscillatorwithnon‐uniformslowwavestructureandaresonantreflector,physicsofplasma,22,014502(2015)。
技术实现要素：
为了解决传统的对hpm产生器件进行多目标优化方法的优化结果具有很大不确定性的问题，本发明提供了一种基于多目标优化算法的hpm产生器件优化方法。
本发明的基本构思：
对于hpm产生器件，不同优化目标地位是不相同的。输出功率是hpm产生器件正常工作的最重要的指标，。如果hpm产生器件工作正常，其输出功率会高于一个预设值。因此，本发明将hpm产生器件的输出功率设为主优化目标，在hpm产生器件输出功率较高的情况下，再考虑其它优化目标。
本发明的技术解决方案是：
一种基于多目标优化算法的hpm产生器件优化方法，其特殊之处在于，包括以下步骤：
1)确定hpm产生器件的优化目标、优化参数及其变化范围；优化目标包括hpm产生器件的输出功率；所述优化目标还包括输出信号的频率、起振时间、输出效率和引导磁场强度中的一个或多个，优化参数包括hpm产生器件的几何参数和电参数；
2)建立总的适应性函数fit，将多目标优化转化为单目标优化；
所述总的适应性函数fit根据公式
3)利用脚本功能，建立hpm产生器件的数值模拟模型；
4)利用遗传算法对hpm产生器件进行优化，使得各个优化目标都能达到hpm产生器件所要求达到的目标值。
进一步地，所述步骤2)中其他优化目标对应的适应性函数fi(x)包括输出信号的频率适应性函数f(f)：
上式中，f为器件产生信号的频率，由粒子模拟获得的结果得到；
f1为hpm器件工作频率范围下边界；
f2为0db带宽频率下边界；
f3为0db带宽频率上边界；
f4为hpm器件工作频率范围上边界。
进一步地，所述步骤2)中其他优化目标对应的适应性函数fi(x)还包括起振时间适应性函数t(t)：
上述公式中，t为hpm器件起振时间，定义为hpm器件达到峰值功率的80％所需要的时间，由粒子模拟获得的结果得到；
t0为hpm器件起振时间优化目标，即目标起振时间；
t1为hpm器件优化目标可接受的最大起振时间。
进一步地，所述步骤2)中其他优化目标对应的适应性函数fi(x)还包括输出效率的适应性函数eta(p)：
eta(p)＝p/pin/pe，
其中，p为输出功率；
pin为输入功率；
pe为效率优化目标；
eta(p)为期望效率。
进一步地，所述步骤2)中其他优化目标对应的适应性函数fi(x)还包括引导磁场强度的适应性函数b(b)：
其中：
b为hpm器件引导磁场强度，
b0为hpm器件引导磁场强度优化目标，即目标引导磁场强度；
b1为hpm器件优化目标可接受的最高磁场强度。
进一步地，步骤1)中的hpm产生器件为相对论返波管、速调管、虚阴极振荡器、磁绝缘线振荡器、磁控管或回旋管。
与传统多目标优化方法相比，本发明的优点是：
1、本发明能实现hpm产生器件多目标同时优化，在一个优化过程中，除了可以优化hpm产生器件的输出功率外，还能优化频率、起振时间和/或场强；本发明与传统多目标转换为单目标的方法不同，本发明考虑到对于hpm产生器件而言，输出功率是最重要的优化目标，而其它例如起振时间、频率等也是需要考虑的较重要的优化目标，因此，本发明以hpm器件的输出功率作为主优化参数，将其它优化目标作为归一化修正目标，在优化输出功率时，不会影响和破坏对其它优化目标比如起振时间、频率等，既能保证输出功率这个最重要的优化目标达到指标，又能保证其它优化目标达到指标。
2、本发明适应性函数形式是根据优化目标确定的，因此各个适应性函数的物理意义明确，可以方便地写出目标优化函数，进行适应性函数参数的设置和更改，以适应不同优化要求。
附图说明
图1是本发明实施例一待优化相对论返波管(rbwo)的结构示意图，图中标号：1‐阴极，2-管头，3-反射器，4-慢波结构，5-提取腔，6-输出波导；
图2是本发明实施例一待优化rbwo输出信号的频率适应性函数曲线；
图3是本发明实施例一待优化rbwo起振时间适应性函数曲线；
图4是本发明实施例一待优化rbwo遗传算法优化曲线；
图5是本发明实施例一待优化rbow的输出功率曲线；
图6是本发明实施例一待优化rbwo的频率曲线；
图7是本发明实施例二待优化虚阴极振荡器(vco)的结构示意图，图中标号：7-阴极，8-阳极网，9-调制腔，10-收集极，11-输出波导；
图8是本发明实施例二待优化vco遗传算法优化曲线；
图9是本发明实施例二待优化vco的输出功率曲线；图10是本发明实施例二待优化vco的频率曲线。
具体实施方式
实施例一：
待优化的hpm产生器件为相对论返波管(rbwo)，由阴极1、管头2、反射器3、慢波结构4、提取腔5、输出波导6组成，如附图1所示。采用本发明的方法对其优化的过程如下：
步骤1：根据应用需求，确定相对论返波管(rbwo)的优化目标、优化参数及其变化范围。
优化目标为：(1)输出功率大于3gw；(2)输出信号的频率范围为9.9～10.1ghz；(3)起振时间小于15ns。
优化参数为：(1)反射器3的宽度w1、半径r1以及距离l1；(2)提取腔5的宽度w2、半径r2以及距离l2。
步骤2：建立相对论返波管(rbwo)总的适应性函数fit，将多目标优化转化为单目标优化。
步骤2.1：建立相对论返波管(rbwo)的频率适应性函数；
对于工程实用化的高功率微波器件，必须对微波器件产生信号的频率有要求，因此在器件优化时，需要将器件输出信号的频率限定在特定的频率范围(即优化目标)之内。目前应用于高功率微波器件的遗传算法，只能实现单目标优化。优化结果如果频率不在要求范围之内，即使功率优化结果再好，微波器件也是无效的。考虑到hpm器件输出信号的频率会随着待优化参数向量x的变化而在一定范围内变化，本发明通过hpm器件输出信号的频率适应性函数f(f)对相对论返波管(rbwo)输出信号的频率进行限定，将其输出信号的频率限定在工程实际要求范围内，这也是本发明的核心和创新点之一。
hpm器件输出信号的频率适应性函数f(f)为归一化梯形窗口函数，频率适应性函数f(f)的表达式为：
上式中，f为器件产生信号的频率，由粒子模拟获得的结果得到；
f1为hpm器件工作频率范围下边界；
f2为0db带宽频率下边界；
f3为0db带宽频率上边界；
f4为hpm器件工作频率范围上边界。
通过上述频率适应性函数f(f)可对器件输出信号的频率进行限制，如果频率在目标范围之内，f(f)值为1，如果频率不在目标范围内，f(f)值<1，从图2函数图就能看出。如果频率适应性函数f(f)值小于1，就会使得总的适应性函数fit减小，从而使得这个算例竞争优势下降。
本实施例中，根据步骤1中的优化目标(2)可以确定rbwo的频率适应性函数f(f)的中心频率为10ghz，0db带宽为0.2ghz，底宽为2ghz，从而得到f1，f2，f3，f4分别为9.0ghz，9.9ghz，10.1ghz，11ghz，如附图2所示。
步骤2.2：建立相对论返波管(rbwo)的起振时间适应性函数；
hpm器件起振时间随待优化参数向量x的变化规律目前没有理论进行精确预测，因而hpm器件起振时间主要是根据粒子模拟结果进行判断。为保证hpm器件的能量效率，hpm器件设计时，要求hpm器件起振时间小于一定时间(即优化目标)。
hpm器件起振时间适应性函数表达式为：
上述公式中，t为hpm器件起振时间，定义为hpm器件达到峰值功率的80％所需要的时间，由粒子模拟获得的结果得到；
t0为hpm器件起振时间优化目标(也就是目标起振时间)；
t1为hpm器件优化目标可接受的最大起振时间；
在本实施例中，t0，t1分别为15ns，25ns，根据上述公式建立rbwo的起振时间适应性函数，如附图3所示，图3中目标起振时间t0＝15ns，优化目标可接受的最大起振时间t1＝25ns；从图3中可以看出，当起振时间小于t0，则rbwo的起振时间适应性函数的函数值保持不变；当起振时间大于t0，则rbwo的起振时间适应性函数的函数值线性下降；当起振时间大于t1时，rbwo的起振时间适应性函数的函数值为0。
步骤2.3建立相对论返波管(rbwo)总的适应性函数fit；
hpm器件总的适应性函数fit根据公式
将步骤2.1和2.3得到的相对论返波管(rbwo)频率适应性函数和起振时间适应性函数代入
步骤3：利用脚本功能，建立相对论返波管(rbwo)的数值模拟模型，建立数值模拟模型的方法为现有公知方法。并对数值模拟模型进行数值计算。
步骤4：利用遗传算法进行输出功率、输出信号的频率、起振时间优化，使得它们都能达到步骤1所确定的目标值。
对hpm产生器件进行优化，就是通过不同的参数组合，得到不同的hpm产生器件结构。传统的遗传算法是公知技术，如文献1采用单目标遗传算法优化，但是不满足我们对hpm产生器件多目标同时优化的需求。
本发明的多目标优化过程是利用传统遗传算法进行不断迭代修正的过程，当满足优化目标之后，遗传算法停止迭代。
上述步骤2确定遗传算法的优化目标以及优化参数设置，在步骤2确定的优化目标函数的判断下，步骤3、步骤4需要循环迭代。遗传算法产生一组参数，步骤3根据这一组参数建立数值模拟模型，步骤4对步骤3建立的数值模拟模型进行计算，分析计算结果，生成新的一组参数，步骤3再根据新的一组参数建立数值模拟模型，以此类推，循环往复。
上述一组参数具体是指x的一组确定数值点(x为需要优化的向量)，对于实施例一来说，这一组参数就是反射器3的宽度w1、半径r1以及距离l1，提取腔5的宽度w2、半径r2以及距离l2这六个参数具体的一组数值。
遗传算法优化曲线如附图4所示，随着遗传算法的进行，相对论返波管(rbwo)的适应性值不断提高，表明相对论返波管(rbwo)越来越接近优化目标。附图5为经遗传算法优化后得到的相对论返波管(rbwo)输出功率图，从图5中可以看出，rbwo输出功率为3.2gw，大于目标值；起振时间为12ns，优于步骤1所确定的优化目标。
附图6给出了经遗传算法优化后相对论返波管(rbwo)的输出信号频率谱，信号频率为9.94ghz，在优化目标范围内。
从上述优化结果看出，相对论返波管(rbwo)的三个优化目标：输出功率、输出信号的频率、器件的起振时间都达到优化预期要求，实现了相对论返波管(rbwo)的多目标优化。
实施例二：
待优化的hpm产生器件为虚阴极振荡器(vco)，由阴极7、阳极网8、调制腔9、收集极10、输出波导11组成，如附图7所示。采用本发明的方法对其优化的过程如下：
步骤1：根据应用需求，确定虚阴极振荡器(vco)的优化目标、优化参数及其变化范围。
优化目标为：(1)输出功率大于6.5gw；(2)输出信号的频率范围为2.8～3.2ghz。
优化参数为：虚阴极调制腔的长度z1、z2，以及径向位置r1、r2。
步骤2：建立虚阴极振荡器(vco)总的适应性函数fit，将多目标优化转化为单目标优化。
步骤2.1：建立虚阴极振荡器(vco)的频率适应性函数；
此步骤与实施例一步骤2.1方法相同，在本实施例中，根据优化目标可以确定vco的频率适应性函数f(f)的中心频率为3ghz，0db带宽为0.2ghz，底宽为1ghz，从而得到f1，f2，f3，f4分别为2.5ghz，2.9ghz，3.1ghz，3.5ghz。
步骤2.2建立虚阴极振荡器(vco)总的适应性函数fit；
此步骤与实施例一步骤2.3方法相同，本实施例中fi(x)为虚阴极振荡器(vco)的频率适应性函数。
步骤3：利用脚本功能，建立虚阴极振荡器(vco)的数值模拟模型，建立数值模拟模型的方法为现有公知方法。
步骤4：利用遗传算法进行输出功率和输出信号的频率优化，使得它们都能达到步骤1所确定的目标值。
遗传算法优化曲线如附图8所示，随着遗传算法的进行，虚阴极振荡器(vco)的适应性值不断提高，表明虚阴极振荡器(vco)越来越接近优化目标。附图9为经遗传算法优化后得到的虚阴极振荡器(vco)输出功率图，从图9可以看出，虚阴极振荡器(vco)输出功率为7gw，大于目标值。附图10给出了经遗传算法优化后虚阴极振荡器(vco)输出信号频率谱，信号频率为3.05ghz，在优化目标范围内。
从上述优化结果看出，虚阴极振荡器(vco)的两个优化目标功率：输出功率和输出信号的频率都达到了优化预期要求，实现了虚阴极振荡器(vco)的多目标优化。