发明名称 一种基于NSGA‑Ⅱ的微波窗快速优化设计方法
摘要 本发明属于微波电真空领域微波窗的模拟设计与优化技术领域,具体涉及一种基于NSGA‑Ⅱ的微波窗快速优化设计方法。本发明结合模式匹配法与NSGA‑Ⅱ,并根据微波窗的性能要求,提出了两个特有的互不矛盾,相互协调的目标优化函数,充分发挥了GA算法优胜劣汰的机制,保证算法逼近全局最优解,以实现各类微波窗的快速优化设计。本发明的有益效果是:实现各类微波窗的快速优化设计,解决了微波窗设计过程中采用纯数值电磁仿真软件存在的耗时长、计算机资源消耗大的问题,以及采用模式匹配法结合传统单目标优化GA算法不能满足微波窗宽频带的要求的局限。同时,本发明也适用于其他宽带微波器件的优化设计。
申请公布号 CN106503359A 申请公布日期 2017.03.15
申请号 CN201610948622.6 申请日期 2016.10.26
申请人 电子科技大学 发明人 巩进杰;朱小芳;胡权;胡玉禄;李斌
分类号 G06F17/50(2006.01)I 主分类号 G06F17/50(2006.01)I
代理机构 电子科技大学专利中心 51203 代理人 闫树平
主权项 一种基于NSGA‑Ⅱ的微波窗快速优化设计方法,包括以下步骤:S1、初始化微波窗结构参数,确定优化参数及其范围,设定模式匹配法计算相关参数:在对确定微波窗结构进行优化计算时,首先初始化微波窗结构参数,同时确定优化参数以及每个优化参数的范围,优化参数构成式1中决策向量X;给出模式匹配法计算所需的微波窗工作频率范围[f<sub>min</sub>,f<sub>max</sub>],频率采样步长Δf和模式数m;<maths num="0001"><math><![CDATA[<mrow><mi>g</mi><mo>=</mo><munder><mrow><mi>m</mi><mi>i</mi><mi>n</mi></mrow><mi>x</mi></munder><mo>{</mo><munder><mrow><mi>m</mi><mi>a</mi><mi>x</mi></mrow><mrow><msub><mi>f</mi><mi>min</mi></msub><mo>&lt;</mo><msub><mi>f</mi><mi>i</mi></msub><mo>&lt;</mo><msub><mi>f</mi><mrow><mi>m</mi><mi>a</mi><mi>x</mi></mrow></msub></mrow></munder><mo>&lsqb;</mo><mo>|</mo><mi>&Gamma;</mi><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>f</mi><mo>)</mo></mrow><mo>|</mo><mo>&rsqb;</mo><mo>}</mo><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001141579980000011.GIF" wi="1310" he="119" /></maths>S2、设定NSGA‑Ⅱ算法参数:设定NSGA‑Ⅱ算法参数包括设定种群数目N、最大进化代数M、交叉率P<sub>c</sub>和变异率P<sub>m</sub>,设置变量s记录迭代次数,即遗传算法基本参数;然后再结合S1,NSGA‑Ⅱ算法将在给定的优化参数范围内随机产生N组决策向量x的集合,将其记为P={x<sub>k</sub>|k=1,2,…,N},其中x<sub>k</sub>=[x<sub>1</sub>,x<sub>2</sub>,…,x<sub>v</sub>]<sup>T</sup>,v表示所要优化的参数数目;S3、确定目标函数:根据电压驻波比VSWR的定义<maths num="0002"><math><![CDATA[<mrow><mi>V</mi><mi>S</mi><mi>W</mi><mi>R</mi><mo>=</mo><mfrac><mrow><mn>1</mn><mo>+</mo><mrow><mo>|</mo><msub><mi>&Gamma;</mi><mi>L</mi></msub><mo>|</mo></mrow></mrow><mrow><mn>1</mn><mo>-</mo><mrow><mo>|</mo><msub><mi>&Gamma;</mi><mi>L</mi></msub><mo>|</mo></mrow></mrow></mfrac><mo>&GreaterEqual;</mo><mn>1</mn><mo>,</mo><mn>0</mn><mo>&le;</mo><mo>|</mo><msub><mi>&Gamma;</mi><mi>L</mi></msub><mo>|</mo><mo>&le;</mo><mn>1</mn><mo>,</mo><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001141579980000012.GIF" wi="1366" he="134" /></maths>其中,Γ<sub>L</sub>是器件端口反射系数,选择第一个目标函数为g<sub>1</sub><maths num="0003"><math><![CDATA[<mrow><msub><mi>g</mi><mn>1</mn></msub><mo>=</mo><mi>m</mi><mi>i</mi><mi>n</mi><mrow><mo>(</mo><munder><mo>&Sigma;</mo><mrow><msub><mi>f</mi><mrow><mi>m</mi><mi>a</mi><mi>x</mi></mrow></msub><mo>&le;</mo><msub><mi>f</mi><mi>i</mi></msub><mo>&le;</mo><msub><mi>f</mi><mi>min</mi></msub></mrow></munder><msub><mi>VSWR</mi><msub><mi>f</mi><mi>i</mi></msub></msub><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>3</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001141579980000013.GIF" wi="1253" he="159" /></maths>式中,f<sub>min</sub>,f<sub>max</sub>,f<sub>i</sub>同式1,<img file="FDA0001141579980000019.GIF" wi="154" he="63" />表示采样频点f<sub>i</sub>对应的电压驻波比;该目标函数使得在频带范围内,各个采样点对应电压驻波比之和具有最小值;为了保证带宽,在第一个目标函数的基础上,增加第二个目标函数g<sub>2</sub><maths num="0004"><math><![CDATA[<mrow><msub><mi>g</mi><mn>2</mn></msub><mo>=</mo><mi>m</mi><mi>i</mi><mi>n</mi><mrow><mo>(</mo><munder><mo>&Sigma;</mo><mrow><msub><mi>f</mi><mrow><mi>m</mi><mi>a</mi><mi>x</mi></mrow></msub><mo>&le;</mo><msub><mi>f</mi><mi>i</mi></msub><mo>&le;</mo><msub><mi>f</mi><mi>min</mi></msub></mrow></munder><msub><mi>&delta;</mi><msub><mi>f</mi><mi>i</mi></msub></msub><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>4</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001141579980000014.GIF" wi="1241" he="159" /></maths>其中,f<sub>min</sub>,f<sub>max</sub>,f<sub>i</sub>同式1,<img file="FDA00011415799800000110.GIF" wi="52" he="63" />定义为<maths num="0005"><math><![CDATA[<mrow><msub><mi>&delta;</mi><msub><mi>f</mi><mi>i</mi></msub></msub><mo>=</mo><mfenced open = "{" close = ""><mtable><mtr><mtd><mrow><mn>1</mn><mo>,</mo><msub><mi>VSWR</mi><msub><mi>f</mi><mi>i</mi></msub></msub><mo>&GreaterEqual;</mo><mi>C</mi></mrow></mtd></mtr><mtr><mtd><mrow><mn>0</mn><mo>,</mo><msub><mi>VSWR</mi><msub><mi>f</mi><mi>i</mi></msub></msub><mo>&lt;</mo><mi>C</mi></mrow></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>5</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001141579980000015.GIF" wi="1214" he="166" /></maths>式6中1<C≤1.1,当某一频点f<sub>i</sub>对应电压驻波比<img file="FDA0001141579980000016.GIF" wi="150" he="63" />大于给定阈值时,<img file="FDA0001141579980000017.GIF" wi="54" he="71" />取值为1,反之,<img file="FDA0001141579980000018.GIF" wi="52" he="70" />取值为0;目标函数g<sub>2</sub>使得在采样频点内VSWR超过阈值的频点数目最少;取目标函数g<sub>2</sub>的阈值为Q<sub>0</sub>,0.2Q<Q<sub>0</sub><Q,<img file="FDA0001141579980000021.GIF" wi="387" he="127" />Q即采样频点数目;S4、并行调用模式匹配法计算集合P中N组决策向量的目标函数值:NSGA‑2算法调用模式匹配法计算集合P中每个决策向量的目标函数值时,将此过程并行化处理,其中x的下标i1、i2、i3、i4、i5、i6表示并行处理时对集合P中x的分组边界序号;S5、根据S4的结果,取当前集合P中的最优解,判断其对应的目标函数值g<sub>2</sub>是否小于等于Q<sub>0</sub>;若是,输出结果,结束;若否,进一步判断NSGA‑II算法是否迭代到指定次数,若达到,结束当前迭代,调整优化参数范围,或增加迭代次数,重复S1‑S5,直到找到满足要求的解为止;若未迭代到指定次数,NSGA‑II算法对集合P执行遗传操作,形成新的优化参数集合P',令P=P',转到S4,重复执行S4、S5。
地址 611731 四川省成都市高新区(西区)西源大道2006号