一种基于粒子群算法的无线感应链路参数实时估算方法技术

本发明专利技术提出了一种基于粒子群算法的无线感应链路参数实时估算方法。本发明专利技术输入粒子总数，迭代次数以及设置参数的搜索边界；输入作为粒子的各系统参数集中的系统参数，根据输入电流估算算法得到各粒子的计算输入电流和实际输入电流；计算各粒子计算输入电流与实际输入电流的均方差，比较各均方差，其中均方差最小的最优粒子即为估算得出的无线感应链路系统参数集。本发明专利技术可以对无线感应链路的系统参数进行实施估算，以帮助发射端控制电力传输，这比仅根据离线测量的参数数据控制电力传输更精确也更方便。本发明专利技术实施测量系统参数的方式，还能有效避免了电气元件的制造公差对系统性能的影响，以保证系统高效运行。

A real-time estimation method of wireless inductive link parameters based on particle swarm optimization

【技术实现步骤摘要】
一种基于粒子群算法的无线感应链路参数实时估算方法
本专利技术属于无线电能传输
，尤其涉及一种基于粒子群算法的无线感应链路参数实时估算方法。
技术介绍
无线感应链路已广泛应用于不同功率水平的场景，如生物医学植入，照明，运输和感应加热等。装置在设计时通常强调在最佳电能传输效率的情况下运行系统。为了调节负载功率，需要进行负载侧电参数的即时测量，并建立发射端和接收端之间的通信链路，以向发送端提供接收端的工作情况。为了最大限度地降低电路复杂性，目前提出了一些利用发送器侧电压和电流信息来估计两个线圈之间的互感，获得负载侧条件以及调节负载功率的计算算法。然而，这些方法通常假设元件值是先验已知的，或者在估计过程中采用标称元件值。由于制造公差，温度效应和老化，电子元件的参数通常都会受到影响，这种不可避免的问题会降低链路的系统性能。众所周知，无线电感链路的功率传输效率和功率传输特性由若干内在和外在因素决定，如耦合系数，品质因数，发射和接收线圈的匹配情况以及工作频率等。这些参数的变化将改变发射端和接收端侧的谐振频率。而发射和接收端谐振频率的不匹配会导致负载功率大幅度下降。
技术实现思路
为了解决上述技术问题，本专利技术提出了一种基于粒子群算法的无线感应链路参数实时估算方法，以实时监控无线感应链路系统参数的变化，以发射端控制电力传输。为了实现上述目的，本专利技术采用的技术方案为一种基于粒子群算法的无线感应链路参数实时估算方法。本专利技术通过改进的粒子群优化算法用于最小化实际输入电流和计算输入电流之间的误差来估算系统参数；通过输入电流估算算法计算输入电流的时间序列，其中包括实际输入电压的时间序列和改进的粒子群优化算法提供的系统参数，具体包括以下步骤：步骤1：输入粒子总数，迭代次数，设置参数的搜索边界；步骤2：输入各粒子作为系统参数，根据输入电流估算算法得到各粒子的计算输入电流和实际输入电流；步骤3：计算各粒子的均方差，比较各粒子的均方差，得出当前最优粒子；步骤4：将当前最优粒子引入扰动得到扰动后新粒子；步骤5：分别计算扰动后新粒子和当前最优粒子的均方差，计算最优粒子集被扰动集；步骤6：根据粒子的速度更新公式和位置更新公式更新粒子，得到估算得出的无线感应链路系统参数；作为优选，步骤1中所述输入粒子总数为N，迭代次数为δ；步骤1中所述的参数为：L1为发射线圈自感，ω1为发射线圈的谐振频率，ω2为接收线圈的谐振频率，k为发射和接收线圈的耦合系数；步骤1中所述设置参数的搜索边界为：L1的搜索上界L1max和搜索下界L1min，ω1的搜索上界ω1max和搜索下界ω1min，ω2的搜索上界ω2max和搜索下界ω2min，k的搜索上界kmax和搜索下界kmin；初始化全部N个粒子，得粒子集SP＝{P1,P2,...,PN}，i＝1,2,3,...,N，Pi为第i个粒子，Pi中的参数都在其相应的搜索边界内随机初始化，Pi的具体定义为：Pi＝{L1,i,ω1,i,ω2,i,ki}其中，L1,i为第i个粒子发射线圈自感，ω1,i为第i个粒子发射线圈的谐振频率，ω2,i为第i个粒子接收线圈的谐振频率，ki为第i个粒子发射和接收线圈的耦合系数；作为优选，步骤2输入电流估算算法步骤如下：步骤2.1：使用电流表和电压表测量得输入电流，输入电压以及输入电压的幅值；步骤2.1中所述输入电流为i1，所述输入电压为vin，所述输入电压的幅值为Vin；步骤2.2：根据输入电流和输入电压，通过Goertzel算法计算输入电压的基波和输入电流的基波，并计算输入电压与输出电入之间的相位角；所述输入电压的基波为vin.fund，所述输入电流的基波为i1.fund；所述输入电压与输出电入之间的相位角为：其中，∠i1,fund为实际输入电流的基波相位角，∠vin,fund为实际输入电压的基波相位角；步骤2.3：根据无线感应链路时域模型，输入改进的粒子群优化算法中的粒子；步骤2.3中所述改进的粒子群优化算法的粒子为：Pi＝{L1i,ω1i,ω2i,ki}其中，L1i为第i个粒子发射线圈自感，ω1i为第i个粒子发射线圈的谐振频率，ω2i为第i个粒子接收线圈的谐振频率，ki为第i个粒子发射和接收线圈的耦合系数。输入设为定值的其他系统参数：发射电路等效电容C1，接收电路等效电容C2，发射电路等效电阻r1，接收电路等效电阻r2，发射电路品质因数Q1，接收电路品质因数Q2，系统开关频率ωs，以及无线感应链路的互感系数M；以及输入给定的初始条件，求解无线感应链路时域模型的电流系数K′1～K′8，接收电路等效电容电压非齐次微分方程的一般解λ1～λ4，以及链路感应系数GL1～GL4；上述初始条件：其中vL1(ωt)为发射电路等效电感电压的时间函数，vL2(ωt)为接收电路等效电感电压的时间函数，vC1(ωt)为发射电路等效电容电压的时间函数，vC2(ωt)为接收电路等效电容电压的时间函数；步骤2.4：根据i1i,cal(ωst)的计算公式：其中，K′1～K′8为无线感应链路时域模型的电流系数，λ1～λ4为接收电路等效电容电压非齐次微分方程的一般解，GL1～GL4为链路感应系数，L1i为第i个粒子发射线圈自感，ki为第i个粒子发射和接收线圈的耦合系数，ωs系统开关频率。输入已知的系统参数和步骤4求得的无线感应链路时域模型参数，计算相位为ωst＝0时的输入电流值i1i,cal(0)。步骤2.5：为ωst增加一个根据i1i,cal(ωst)的计算公式计算该相位的输入电流值i1i,cal(ωst)，其中m为任意有理数；。步骤2.6：重复步骤2.5直到相位为ωst＝π，得到ωst∈[0,π]时计算输入电流的时间序列i1i,cal(ωst)；步骤2.7：ωst∈[π,2π]时计算输入时间序列i1i,cal(ωst)，ωst∈[0,π]时计算输入电流时间序列i1i,cal(ωst)的负值即-i1i,cal(ωst)，于是得到ωst∈[0,2π]时计算输入电流的时间序列i1i,cal(ωst)。作为优选，步骤3中所述得出的当前最优粒子为：输入步骤2得到实际输入电流i1(ωst)和各粒子的计算输入电流i1i,cal(ωst)，求各粒子的均方差为：根据各粒子的均方差建立粒子均方差值集合：SMSE＝{MSE(P1),MSE(P2),...,MSE(Pi),...,MSE(PN)}设当前迭代次数为g，g＝1,2,...,δ，当前最优粒子为根据目标函数选择均方差集合SMSE中均方差最小的粒子作为当前第g代最优粒子作为优选，步骤4中所述将当前最优粒子引入扰动得到扰动后新粒子具体为：根据式为当前最优粒子引入一个扰动ΔPB(g)，产生扰动后新粒子作为优选，步骤5中所述计算最优粒子集被扰动集为：步骤5中所述扰动后新粒子的均方差为步骤5中所述当前最优粒子的均方差为MSE(PG(g))；若则当前最优粒子集PG(g)被扰动集：SB＝{PG(g)+ΔPB(g),PG(g)+2ΔPB(g),PG(g)+4ΔPB(g),...,PG(g)+2jΔPB(g)}中均方差最小的一个粒子取代，即其中，p＝PG(g)+2rΔPB(g)，r＝0,1,2,...,j，PG(g)+2jΔPB(g)中元素搜索边界之内而PG(g)+2j+1ΔPB(g)中元素本文档来自技高网...

【技术保护点】
1.一种基于粒子群算法的无线感应链路参数实时估算方法，其特征在于，包括以下步骤：步骤1：输入粒子总数，迭代次数，设置参数的搜索边界；步骤2：输入各粒子作为系统参数，根据输入电流估算算法得到各粒子的计算输入电流和实际输入电流；步骤3：计算各粒子的均方差，比较各粒子的均方差，得出当前最优粒子；步骤4：将当前最优粒子引入扰动得到扰动后新粒子；步骤5：分别计算扰动后新粒子和当前最优粒子的均方差，计算最优粒子集被扰动集；步骤6：根据粒子的速度更新公式和位置更新公式更新粒子，得到估算得出的无线感应链路系统参数。

【技术特征摘要】
1.一种基于粒子群算法的无线感应链路参数实时估算方法，其特征在于，包括以下步骤：步骤1：输入粒子总数，迭代次数，设置参数的搜索边界；步骤2：输入各粒子作为系统参数，根据输入电流估算算法得到各粒子的计算输入电流和实际输入电流；步骤3：计算各粒子的均方差，比较各粒子的均方差，得出当前最优粒子；步骤4：将当前最优粒子引入扰动得到扰动后新粒子；步骤5：分别计算扰动后新粒子和当前最优粒子的均方差，计算最优粒子集被扰动集；步骤6：根据粒子的速度更新公式和位置更新公式更新粒子，得到估算得出的无线感应链路系统参数。2.根据权利要求1所述的基于粒子群算法的无线感应链路参数实时估算方法，其特征在于：步骤1中所述输入粒子总数为N，迭代次数为δ；步骤1中所述的参数为：L1为发射线圈自感，ω1为发射线圈的谐振频率，ω2为接收线圈的谐振频率，k为发射和接收线圈的耦合系数；步骤1中所述设置参数的搜索边界为：L1的搜索上界L1max和搜索下界L1min，ω1的搜索上界ω1max和搜索下界ω1min，ω2的搜索上界ω2max和搜索下界ω2min，k的搜索上界kmax和搜索下界kmin；初始化全部N个粒子，得粒子集SP＝{P1,P2,...,PN}，i＝1,2,3,...,N，Pi为第i个粒子，Pi中的参数都在其相应的搜索边界内随机初始化，Pi的具体定义为：Pi＝{L1,i,ω1,i,ω2,i,ki}其中，L1,i为第i个粒子发射线圈自感，ω1,i为第i个粒子发射线圈的谐振频率，ω2,i为第i个粒子接收线圈的谐振频率，ki为第i个粒子发射和接收线圈的耦合系数。3.根据权利要求1所述的基于粒子群算法的无线感应链路参数实时估算方法，其特征在于：步骤2输入电流估算算法步骤如下：步骤2.1：使用电流表和电压表测量得输入电流，输入电压以及输入电压的幅值；步骤2.1中所述输入电流为i1，所述输入电压为vin，所述输入电压的幅值为Vin；步骤2.2：根据输入电流和输入电压，通过Goertzel算法计算输入电压的基波和输入电流的基波，并计算输入电压与输出电入之间的相位角；所述输入电压的基波为vin.fund，所述输入电流的基波为i1.fund；所述输入电压与输出电入之间的相位角为：其中，∠i1,fund为实际输入电流的基波相位角，∠vin,fund为实际输入电压的基波相位角；步骤2.3：根据无线感应链路时域模型，输入改进的粒子群优化算法中的粒子；步骤2.3中所述改进的粒子群优化算法的粒子为：Pi＝{L1i,ω1i,ω2i,ki}其中，L1i为第i个粒子发射线圈自感，ω1i为第i个粒子发射线圈的谐振频率，ω2i为第i个粒子接收线圈的谐振频率，ki为第i个粒子发射和接收线圈的耦合系数；输入设为定值的其他系统参数：发射电路等效电容C1，接收电路等效电容C2，发射电路等效电阻r1，接收电路等效电阻r2，发射电路品质因数Q1，接收电路品质因数Q2，系统开关频率ωs，以及无线感应链路的互感系数M；以及输入给定的初始条件，求解无线感应链路时域模型的电流系数K′1～K′8，接收电路等效电容电压非齐次微分方程的一般解λ1～λ4，以及链路感应系数GL1～GL4；上述初始条件：其中vL1(ωt)为发射电路等效电感电压的时间函数，vL2(ωt)为接收电路等效电感电压的时间函数，vC1(ωt)为发射电路等效电容电压的时间函数，vC2(ωt)为接收电路等效电容电压的时间函数；步骤2.4：根据i1i,cal(ωst)的计算公式：其中，K′1～K′8为无线感应链路时域模型的电流系数，λ1～λ4为接收电路等效电容电压非齐次微分方程...

【专利技术属性】
技术研发人员：代中余，王军华，朱仁勋，曲皓玥，周海阔，黄虹，
申请(专利权)人：武汉大学，
类型：发明
国别省市：湖北,42