发明名称 频率分叉时电动汽车无线充电桩效率极值点跟踪方法
摘要 本发明公开了一种频率分叉时电动汽车无线充电桩效率极值点跟踪方法,将一般粒子群算法中的粒子群规模分开设定,分别为最大粒子群规模<i>Nmax</i>=30和最小粒子群规模<i>Nmin</i>=2,算法在运行中随着迭代次数的增加逐渐减小粒子群规模,精简了算法,加快了算法后期收敛速度。本发明粒子群算法使得粒子规模选取有据可依,在搜索后期,加快收敛速度,搜索时间减小。
申请公布号 CN105141016A 申请公布日期 2015.12.09
申请号 CN201510559055.0 申请日期 2015.09.06
申请人 河南师范大学 发明人 王萌;孙长兴;施艳艳;梁洁
分类号 H02J7/02(2006.01)I;H02J5/00(2006.01)I;B60L11/18(2006.01)I;G06Q50/06(2012.01)I 主分类号 H02J7/02(2006.01)I
代理机构 新乡市平原专利有限责任公司 41107 代理人 路宽
主权项 频率分叉时电动汽车无线充电桩效率极值点跟踪方法,其特征在于:将一般粒子群算法中的粒子群规模分开设定,分别为最大粒子群规模Nmax=30和最小粒子群规模Nmin=2,粒子群规模随着迭代次数增加而逐渐减小,其具体实施步骤为:(1)、初始化算法,包括设定粒子种群维数D,最大迭代次数MaxNum,同时限定粒子最大速度v<sub>max</sub>,初始化惯性权重w,随机初始化粒子的速度v和粒子的位置;(2)、初始化粒子群规模Nmax为30,并设定最小粒子群规模Nmin=2,初始化迭代次数t=1;(3)、计算当前种群每个粒子的适应度函数值f<sub>i</sub>,f<sub>i</sub>表示第i个粒子的适应度函数值,其中适应度函数值f<sub>i</sub>根据适应度函数<img file="FDA0000796438740000011.GIF" wi="471" he="164" />计算得到,式中<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><msub><mi>Z</mi><mn>1</mn></msub><mo>=</mo><msub><mi>R</mi><mn>1</mn></msub><mo>+</mo><msub><mi>R</mi><mi>s</mi></msub><mo>+</mo><msub><mi>j&omega;L</mi><mn>1</mn></msub><mo>+</mo><mfrac><mn>1</mn><mrow><msub><mi>j&omega;C</mi><mn>1</mn></msub></mrow></mfrac><mo>,</mo><msub><mi>Z</mi><mn>2</mn></msub><mo>=</mo><msub><mi>R</mi><mn>2</mn></msub><mo>+</mo><msub><mi>R</mi><mi>L</mi></msub><mo>+</mo><mi>j</mi><mrow><mo>(</mo><msub><mi>&omega;L</mi><mn>2</mn></msub><mo>-</mo><mfrac><mn>1</mn><mrow><msub><mi>&omega;C</mi><mn>2</mn></msub></mrow></mfrac><mo>)</mo></mrow><mo>,</mo><mi>&omega;</mi><mo>=</mo><mn>2</mn><msub><mi>&pi;f</mi><mi>r</mi></msub><mo>,</mo></mrow>]]></math><img file="FDA0000796438740000012.GIF" wi="1461" he="165" /></maths>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>=0和全局极值f<sub>i‑gbest</sub>=0,个体极值用f<sub>i‑best</sub>表示,所代表的含义是第i个粒子截止到第t次迭代时搜寻到的最优适应度函数值,全局极值用f<sub>i‑gbest</sub>表示,所代表的含义是截止到第t次迭代时,全部粒子搜索到的最优适应度函数值,将步骤(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)、用公式N<sub>present</sub>=Nmax‑(Nmax‑Nmin)*t/(MaxNum)生成下一代粒子种群,其中Npresent为粒子群当前规模,Nmax为最大粒子群规模,Nmin为最小粒子群规模,MaxNum为最大迭代次数,t为当前迭代次数,按公式<maths num="0002" id="cmaths0002"><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="FDA0000796438740000013.GIF" wi="1202" he="109" /></maths>和公式<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><msubsup><mi>x</mi><mi>i</mi><mrow><mi>t</mi><mo>+</mo><mn>1</mn></mrow></msubsup><mo>=</mo><msubsup><mi>x</mi><mi>i</mi><mi>t</mi></msubsup><mo>+</mo><msubsup><mi>v</mi><mi>i</mi><mrow><mi>t</mi><mo>+</mo><mn>1</mn></mrow></msubsup></mrow>]]></math><img file="FDA0000796438740000014.GIF" wi="288" he="79" /></maths>更新下一代粒子种群每个粒子的速度和位置,然后令迭代次数t=t+1,转向步骤(6),其中<img file="FDA0000796438740000015.GIF" wi="94" he="77" />代表t+1次迭代第i个粒子的速度,<img file="FDA0000796438740000016.GIF" wi="67" he="71" />代表当前第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="FDA0000796438740000017.GIF" wi="99" he="79" />代表t+1次迭代第i个粒子位置,<img file="FDA0000796438740000018.GIF" wi="60" he="75" />代表第t次迭代第i个粒子当前位置,w代表惯性权重;(6)、根据公式<img file="FDA0000796438740000021.GIF" wi="680" he="316" />计算粒子适应度函数值的方差之和,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,判断方差是否小于某一精度值ε或者粒子群算法是否达到最大迭代次数,ε根据寻优对象确定,如果否,则转向步骤(3),如果是则转向步骤(7);(7)、输出搜索到最优频率值,即粒子适应度值为全局最优值f<sub>i‑gbest</sub>的粒子位置p<sub>g</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>2mαx</sub>(k)|&gt;Δ是否成立,如果判断结果为是,则转向步骤(1),算法重启,如果判断结果为否,算法转向步骤(7)。
地址 453007 河南省新乡市牧野区建设东路46号
您可能感兴趣的专利