发明名称 基于适应性教学优化的变压器优化设计方法
摘要 本发明公开了一种基于适应性教学优化的变压器优化设计方法,利用适应性教学优化算法来优化设计变压器。在适应性教学优化算法中设计了基于适应性惯性权重和搜索步长的教学搜索算子,并将适应值较优的个体融合到教学搜索算子中。另一方面,在自学搜索算子中同时将最优个体和随机个体的信息引入到高斯变异策略中,以此协调收敛速度与种群多样性之间的平衡。本发明能够提高变压器优化设计的精度,提高变压器优化设计的效率。
申请公布号 CN106528996A 申请公布日期 2017.03.22
申请号 CN201610958380.9 申请日期 2016.11.04
申请人 江西理工大学 发明人 郭肇禄;李大海;杨火根;周才英;鄢化彪;余法红;李康顺
分类号 G06F17/50(2006.01)I;G06N3/12(2006.01)I 主分类号 G06F17/50(2006.01)I
代理机构 赣州凌云专利事务所 36116 代理人 曾上;卢和炳
主权项 一种基于适应性教学优化的变压器优化设计方法,其特征在于:包括以下步骤:步骤1,根据需要优化设计变压器的物理和电气特性构建以下形式的最小化优化目标的数学模型:最小化变压器的优化设计目标函数f(X),并满足优化设计约束条件:h<sub>k</sub>(X)≤0,k=1,2,...,G,其中X=[x<sub>1</sub>x<sub>2</sub>...x<sub>D</sub>]为变压器的优化设计变量组成的向量;D为变压器需要优化设计的变量个数;h<sub>k</sub>(X)≤0为第k个优化设计约束条件,G为优化设计约束条件个数;步骤2,用户初始化参数,所述初始化参数包括变压器所需优化设计的变量个数D,种群大小Popsize,最大评价次数MAX_FEs;步骤3,令当前演化代数t=0;步骤4,令当前评价次数FEs=0;步骤5,随机产生初始种群<img file="FDA0001143651320000011.GIF" wi="666" he="86" />其中:个体下标i=1,2,...,Popsize,并且<img file="FDA0001143651320000012.GIF" wi="578" he="71" />为种群P<sub>t</sub>中的第i个个体,其随机初始化公式为:<maths num="0001"><math><![CDATA[<mrow><msubsup><mi>B</mi><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow><mi>t</mi></msubsup><mo>=</mo><msub><mi>LB</mi><mi>j</mi></msub><mo>+</mo><mi>r</mi><mi>a</mi><mi>n</mi><mi>d</mi><mrow><mo>(</mo><mn>0</mn><mo>,</mo><mn>1</mn><mo>)</mo></mrow><mo>&times;</mo><mrow><mo>(</mo><msub><mi>UB</mi><mi>j</mi></msub><mo>-</mo><msub><mi>LB</mi><mi>j</mi></msub><mo>)</mo></mrow><mo>,</mo></mrow>]]></math><img file="FDA0001143651320000013.GIF" wi="830" he="71" /></maths>其中维数下标j=1,2,...,D,并且D表示变压器有多少个需要优化设计的变量;<img file="FDA0001143651320000014.GIF" wi="62" he="63" />为种群P<sub>t</sub>中的第i个个体,存储了D个优化设计变量的值,rand(0,1)为在[0,1]之间服从均匀分布的随机实数产生函数,LB<sub>j</sub>和UB<sub>j</sub>分别为第j个优化设计变量的取值范围的下界和上界;步骤6,采用智能优化算法中的惩罚函数法计算种群P<sub>t</sub>中每个个体<img file="FDA0001143651320000015.GIF" wi="58" he="62" />的适应值<img file="FDA0001143651320000016.GIF" wi="195" he="62" />其中个体下标i=1,2,...,Popsize;步骤7,当前评价次数FEs=FEs+Popsize;步骤8,保存种群P<sub>t</sub>中的最优个体Best<sup>t</sup>,然后令搜索步长均值因子FSU=0.5,并令惯性权重均值因子WSU=0.5;步骤9,执行基于适应性惯性权重和搜索步长的教学搜索算子,具体步骤如下:步骤9.1,根据适应值从好到差对种群P<sub>t</sub>中所有个体进行排序;步骤9.2,记录每个个体<img file="FDA0001143651320000021.GIF" wi="53" he="62" />在排序后的种群中的序号<img file="FDA0001143651320000022.GIF" wi="147" he="63" />其中个体下标i=1,2,...,Popsize;步骤9.3,按以下公式计算每个个体<img file="FDA0001143651320000023.GIF" wi="53" he="62" />的选择概率SelP<sub>i</sub><sup>t</sup>:<maths num="0002"><math><![CDATA[<mrow><msup><msub><mi>SelP</mi><mi>i</mi></msub><mi>t</mi></msup><mo>=</mo><mfrac><mrow><mi>P</mi><mi>o</mi><mi>p</mi><mi>s</mi><mi>i</mi><mi>z</mi><mi>e</mi><mo>-</mo><msubsup><mi>SRN</mi><mi>i</mi><mi>t</mi></msubsup><mo>+</mo><mn>1</mn></mrow><mrow><mi>P</mi><mi>o</mi><mi>p</mi><mi>s</mi><mi>i</mi><mi>z</mi><mi>e</mi></mrow></mfrac><mo>,</mo></mrow>]]></math><img file="FDA0001143651320000024.GIF" wi="623" he="135" /></maths>其中个体下标i=1,2,...,Popsize;步骤9.4,令搜索步长列表FUList为空,并令惯性权重列表WUList为空;步骤9.5,计算种群P<sub>t</sub>中所有个体的平均值得到均值个体MeanB<sup>t</sup>;步骤9.6,令计数器ki=1;步骤9.7,令搜索步长FS=NormRand(FSU,0.2),其中NormRand(FSU,0.2)表示以FSU为均值,0.2为标准差的高斯随机数产生函数;步骤9.8,令惯性权重WS=NormRand(WSU,0.2),其中NormRand(WSU,0.2)表示以WSU为均值,0.2为标准差的高斯随机数产生函数;步骤9.9,根据种群P<sub>t</sub>中每个个体的选择概率采用轮盘赌策略选择出个体<img file="FDA0001143651320000025.GIF" wi="116" he="63" />步骤9.10,令均值个体因子TF=round(1+rand(0,1)),其中round为取四舍五入的函数;步骤9.11,令试验个体<img file="FDA0001143651320000031.GIF" wi="1035" he="63" />步骤9.12,计算试验个体<img file="FDA0001143651320000032.GIF" wi="130" he="63" />的适应值<img file="FDA0001143651320000033.GIF" wi="266" he="63" />步骤9.13,如果试验个体<img file="FDA0001143651320000034.GIF" wi="130" he="62" />比个体<img file="FDA0001143651320000035.GIF" wi="67" he="62" />更优,则转到步骤9.14,否则转到步骤9.15;步骤9.14,将FS添加到搜索步长列表FUList中,并将WS添加到惯性权重列表WUList中;步骤9.15,在试验个体<img file="FDA0001143651320000036.GIF" wi="131" he="61" />和个体<img file="FDA0001143651320000037.GIF" wi="60" he="62" />之间执行选择操作算子;步骤9.16,令计数器ki=ki+1;步骤9.17,如果计数器ki小于或等于Popsize,则转到步骤9.7,否则转到步骤9.18;步骤9.18,计算搜索步长列表FUList中所有数据的平均值NFU,然后在[0.5,1.0]之间产生一个随机实数KR1;步骤9.19,令搜索步长均值因子FSU=FSU×KR1+NFU×(1‑KR1);步骤9.20,计算惯性权重列表WUList中所有数据的平均值NWU,然后在[0.5,1.0]之间产生一个随机实数KR2;步骤9.21,令惯性权重均值因子WSU=WSU×KR2+NWU×(1‑KR2);步骤9.22,转到步骤10;步骤10,令当前评价次数FEs=FEs+Popsize;步骤11,执行基于高斯变异策略的自学搜索算子,具体步骤如下:步骤11.1,令计数器km=1;步骤11.2,在[1,Popsize]之间随机产生一个正整数RI;步骤11.3,令计数器j=1;步骤11.4,令均值<img file="FDA0001143651320000041.GIF" wi="571" he="71" />步骤11.5,令方差<img file="FDA0001143651320000042.GIF" wi="617" he="70" />其中abs表示取绝对值的函数;步骤11.6,<img file="FDA0001143651320000043.GIF" wi="773" he="71" />其中NormRand(MBV,SDFV)表示以MBV为均值,SDFV为方差的高斯随机数产生函数;步骤11.7,令计数器j=j+1;步骤11.8,如果计数器j小于或等于D,则转到步骤11.4,否则转到步骤11.9;步骤11.9,计算试验个体<img file="FDA0001143651320000044.GIF" wi="139" he="62" />的适应值<img file="FDA0001143651320000045.GIF" wi="282" he="63" />步骤11.10,在试验个体<img file="FDA0001143651320000046.GIF" wi="139" he="63" />和个体<img file="FDA0001143651320000047.GIF" wi="74" he="63" />之间执行选择操作算子;步骤11.11,令计数器km=km+1;步骤11.12,如果计数器km小于或等于Popsize,则转到步骤11.2,否则转到步骤11.13;步骤11.13,转到步骤12;步骤12,令当前评价次数FEs=FEs+Popsize;步骤13,保存种群P<sub>t</sub>中的最优个体Best<sup>t</sup>;步骤14,当前演化代数t=t+1;步骤15,重复步骤9至步骤14直至当前评价次数FEs达到MAX_FEs后结束,执行过程中得到的最优个体Best<sup>t</sup>即为变压器优化设计的结果。
地址 341000 江西省赣州市红旗大道86号