发明名称 高传输效率无线电能发射系统的频率快速搜索方法
摘要 本发明公开了一种高传输效率无线电能发射系统频率快速搜索方法,将一般算法中的粒子群规模分开设定,分别为最大粒子群规模<i>Nmax</i>=30和最小粒子群规模<i>Nmin</i>=2,粒子群规模随着迭代次数增加而逐渐减小,其中粒子群规模的总体减小方式类似于指数型变化的曲线。本算法在搜索前期,粒子群规模变化较慢,有利于全局搜索,在搜索后期,粒子群规模变为最小,减去了冗余粒子,精简了算法,加快了算法后期收敛速度。本发明算法不但使得粒子规模选取有据可依,且算法在搜索前期具有较大自我学能力和社会学能力,在搜索后期,加快收敛速度,算法搜索时间减小。
申请公布号 CN105140972A 申请公布日期 2015.12.09
申请号 CN201510559086.6 申请日期 2015.09.06
申请人 河南师范大学 发明人 王萌;孙长兴;施艳艳;梁洁
分类号 H02J5/00(2006.01)I;G06N3/00(2006.01)I 主分类号 H02J5/00(2006.01)I
代理机构 新乡市平原专利有限责任公司 41107 代理人 路宽
主权项 高传输效率无线电能发射系统的频率快速搜索方法,其特征在于:将一般粒子群算法中的粒子群规模分开设定,分别为最大粒子群规模Nmax=30和最小粒子群规模Nmin=2,粒子群规模随着迭代次数增加而沿着指数曲线的方式逐渐减小,其具体实施步骤为:(1)、初始化算法,包括设定粒子种群维数D,最大迭代次数MaxNum,同时限定粒子最大速度v<sub>max</sub>,初始化惯性权重w;(2)、直接设定粒子群最大规模Nmax为30和粒子群最小规模Nmin为2,随机初始化粒子的速度v和粒子的位置,设定初始粒子群规模为最大规模Nmax=30,初始化迭代次数t=1;(3)、采用适应度函数<img file="FDA0000796517730000011.GIF" wi="461" he="161" />计算当前种群每个粒子的适应度函数值f<sub>i</sub>,f<sub>i</sub>表示第i个粒子的适应度函数值,其中<img file="FDA0000796517730000012.GIF" wi="606" he="147" /><img file="FDA0000796517730000013.GIF" wi="631" he="166" />ω=2πf<sub>r</sub>,f<sub>r</sub>为当前激励频率,ω为激励电源的角频率,M为发射和接收线圈之间的互感,L<sub>1</sub>,L<sub>2</sub>为发射线圈和接收线圈电感,C<sub>1</sub>,C<sub>2</sub>为电容,R<sub>s</sub>为电源内阻,R<sub>L</sub>为负载电阻,R<sub>1</sub>,R<sub>2</sub>为回路中电阻;(4)、用f<sub>i‑best</sub>表示第i个粒子截止到第t次迭代时搜寻到的最优适应度函数值,用f<sub>i‑gbest</sub>表示截止到第t次迭代时,全部粒子搜索到的最优适应度函数值,在粒子群算法开始迭代之前,设定f<sub>i‑best</sub>=0,f<sub>i‑gbest</sub>=0,将步骤(3)中得到的粒子适应度函数值f<sub>i</sub>和个体极值f<sub>i‑best</sub>及全局极值f<sub>i‑gbest</sub>相比较,如果f<sub>i</sub>≤f<sub>i‑best</sub>,那么f<sub>i‑best</sub>=f<sub>i</sub>,p<sub>i</sub>=x<sub>i</sub>,p<sub>i</sub>表示适应度函数值为f<sub>i‑best</sub>的粒子位置,x<sub>i</sub>是所对适应度函数值为f<sub>i</sub>粒子的位置,如果f<sub>i</sub>≤f<sub>i‑gbest</sub>,那么f<sub>i‑gbest</sub>=f<sub>i</sub>,p<sub>g</sub>=x<sub>i</sub>,p<sub>g</sub>是粒子种群中全局最优值为f<sub>i‑gbest</sub>的粒子位置;(5)、按公式<img file="FDA0000796517730000014.GIF" wi="661" he="123" />更新粒子群规模,其中Npresent为粒子群当前规模,Nmax为最大粒子群规模,Nmin为最小粒子群规模,MaxNum为最大迭代次数,t为当前迭代次数,n为控制粒子群规模变化规律的幂指数,通过参数n调节粒子群规模变化的快慢程度,按公式<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><msubsup><mi>v</mi><mi>i</mi><mrow><mi>t</mi><mo>+</mo><mn>1</mn></mrow></msubsup><mo>=</mo><mi>w</mi><mo>*</mo><msubsup><mi>v</mi><mi>i</mi><mi>t</mi></msubsup><mo>+</mo><msub><mi>c</mi><mn>1</mn></msub><mo>*</mo><mi>r</mi><mi>a</mi><mi>n</mi><mi>d</mi><mo>*</mo><mrow><mo>(</mo><msub><mi>p</mi><mi>i</mi></msub><mo>-</mo><msubsup><mi>x</mi><mi>i</mi><mi>t</mi></msubsup><mo>)</mo></mrow><mo>+</mo><msub><mi>c</mi><mn>2</mn></msub><mo>*</mo><mi>r</mi><mi>a</mi><mi>n</mi><mi>d</mi><mo>*</mo><mrow><mo>(</mo><msub><mi>p</mi><mi>g</mi></msub><mo>-</mo><msubsup><mi>x</mi><mi>i</mi><mi>t</mi></msubsup><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000796517730000015.GIF" wi="1208" he="104" /></maths>和公式<img file="FDA0000796517730000016.GIF" wi="292" he="82" />更新各个粒子的速度和位置,然后令迭代次数t=t+1,转向步骤(6),其中<img file="FDA0000796517730000017.GIF" wi="99" he="79" />代表t+1次迭代第i个粒子的速度,<img file="FDA0000796517730000018.GIF" wi="63" he="74" />代表当前第t次迭代第i个粒子的速度,c<sub>1</sub>和c<sub>2</sub>代表学习因子,rand代表[01]之间的随机数,p<sub>i</sub>表示适应度函数值为f<sub>i‑best</sub>的粒子位置,p<sub>g</sub>是粒子种群中全局最优值为f<sub>i‑gbest</sub>的粒子位置,<img file="FDA0000796517730000019.GIF" wi="102" he="77" />代表t+1次迭代第i个粒子位置,<img file="FDA00007965177300000110.GIF" wi="68" he="75" />代表第t次迭代第i个粒子当前位置,w代表惯性权重;(6)、根据公式<img file="FDA0000796517730000021.GIF" wi="677" he="318" />计算粒子适应度函数值的方差之和,f<sub>avg</sub>为全部粒子适应度函数值的平均值,其中如果有(f<sub>i</sub>‑f<sub>avg</sub>)&gt;1,则a=max(f<sub>i</sub>‑f<sub>avg</sub>),否则,a=1,判断方差是否等于0或者算法是否达到最大迭代次数,如果否,则转向步骤(3),如果是则转向步骤(7);(7)、输出搜索到的全局最优值p<sub>g</sub>,p<sub>g</sub>是粒子种群中全局最优值为f<sub>i‑gbest</sub>的粒子位置,即搜索到的最优值对应的频率值;(8)、用电流传感器检测负载电流i<sub>2</sub>的峰值,设Δ为设定的最大电流峰值波动范围,i<sub>2max</sub>为所检测的负载电流峰值,i<sub>2max</sub>(k)为负载的第k个电流周期电流峰值,i<sub>2max</sub>(k+1)为负载的第k+1个电流周期的电流峰值,判断|i<sub>2max</sub>(k+1)|‑|i<sub>2max</sub>(k)|&gt;Δ是否成立,如果判断结果为是,则转向步骤(1),算法重启;如果判断结果为否,算法转向步骤(7)。
地址 453007 河南省新乡市牧野区建设东路46号
您可能感兴趣的专利