一种基于近场定位的磁信道估计方法与流程

本发明属于无线通信领域，具体涉及一种基于近场定位的磁信道估计方法。

背景技术：

随着社会进步与科技发展,我们的日常生活愈发依赖于众多个人移动设备，例如手机，平板电脑和可穿戴设备。尽管每一种设备都可以使生活变得更加方便，但是我们必须记住每天都要给这些设备充电，这是一个经常性且日益重要的负担。如果可以对这些设备进行无线充电，这样就能够减轻这种日常焦虑，也会大大减少对各种充电器的的数量。无线充电技术的主要优势在于它的便捷性和通用性，通过采用无线充电技术，公共移动设备充电站将有可能成为现实。缺点在于无线充电的效率相对于有线充电来说有些偏低，但是，移动设备对低功耗的追求为无线充电技术提供了广阔前景。

现有的无线能量传输技术根据其物理机理可以分为三类：电感耦合，磁共振耦合和电磁辐射。前两类主要利用了电磁的近场非辐射特性。由于磁感效应随距离增加快速减弱，基于电感耦合的无线能量传输通常只能达到厘米级的充电距离。磁共振耦合可以实现米级的能量传输，但对线圈间的距离和对齐角度有严格要求。在采用磁波束成形的无线能量传输电路中，发送端需要先估计磁信道(由线圈间的互感系数决定)，然后再根据信道调整发送线圈上的电流值。发送端的信道估计精度将直接影响接收端的能量接收效率。在多发送线圈的能量传输电路中，现有的两种磁信道估计方法包括1)逐一闭合发送线圈(同时断开其他发送线圈)，根据发送线圈上的观测电压和电流值估计互感系数；2)接收端将感应电流值通过通信链路反馈给发送端，发送端根据kvl方程计算互感系数。

现有的估计方法中，方法(1)需要对所有发送线圈的电压和电流进行观测，计算量较大；方法(2)需要建立通信反馈电路，受用面较窄。

技术实现要素：

本发明的目的在于提出一种基于近场定位的磁信道估计方法，先利用在部分发送线圈上的电压和感应电流估计出接收线圈的位置，再根据接收线圈的位置以及其他发送线圈与接收线圈的相对位移，估计出该接收线圈和其他发送线圈之间的互感系数，计算量较小且不需要建立反馈链路。

本发明采用的技术方案为一种基于近场定位的磁信道估计方法，包括如下步骤：

(1)初始时从发送线圈阵列中选择3个发送线圈，加载带有固定频率的交流信号，并假设唯一的接收线圈摆放在任意位置；

(2)依次观测3个发送线圈上的电流，然后观测实际加载到发送线圈上的电压，每个发送线圈的观测实施多次；

(3)根据发送线圈与接收线圈之间的互感与磁场强度的理论关系，利用最小二乘法，估计每个发送线圈与接收线圈间的磁场强度，继而估计接收线圈的位置；

(4)利用所述步骤(3)估计的位置以及其他发送线圈与接收线圈的相对位移，算出其他发送线圈与接收线圈间的互感，根据计算出的互感调整所有发送线圈上电压或电流的幅度与相位，计算接收功率与发送功率，进而计算能量传输效率，使调整后的能量传输效率优于调整前。

进一步的，所述步骤(2)中，包括如下步骤：

不同时刻的发送线圈上的电流向量是正交的：每个时刻只有一个发送线圈闭合，观测的是发送线圈上的电流，发送线圈上的总电压，以及实际加载到发送线圈上的电压，所述实际加载到发送线圈上的电压指的是总电压减去电路中电阻上的分压，每个发送线圈的观测实施多次。

进一步的，所述步骤(3)中，发送线圈与接收线圈间的kvl(基尔霍夫定律)方程表述为：

式中：ir为接收线圈上的电流，rr为接收线圈上的负载阻抗，j为复数的虚部，ω为施加在发送线圈上的交流信号的角频率，为第n个发送线圈与接收线圈间的互感，为第n个发送线圈上的电流，为第n个发送线圈上的阻抗，vn为加载到第n个发送线圈上的总电压；

上面两个式子简化为：

令其中，l表示第l次观测，yn(l)为第l次观测时实际加载到第n个发送线圈上的真实电压，vn(l)为第l次观测时加载到第n个发送线圈电路上的总电压，表示第l次观测时第n个发送线圈上的电流，将式子表示为：

用表示第l次观测时实际加载到第n个发送线圈上的电压的观测值，将该问题描述为最小二乘法：

该式表示观测值与实际值误差的2-范数平方和，式中，l为观测次数，上述式子就是估算磁场强度所利用的最小二乘法表达式。

进一步的，所述步骤(3)中，互感与磁场强度之间的关系为：

式中：vind为发送线圈上的电流在接收线圈上产生的感应电压，μ0为空气的磁导率，ntx为发送线圈的匝数，nrx为接收线圈的匝数，arx为接收线圈的面积，arx＝πb²，b是接收线圈的半径，it为发送线圈上的电流，hint为磁场强度；根据上式和得到磁场强度与互感之间的关系为：

进一步的，所述步骤(3)中发送线圈与接收线圈间的磁场强度与接收线圈的位置的关系为：

式中：a为发送线圈的半径，δ、d分别为接收线圈相对于发送线圈的横向位移与纵向位移，m为模量，且0≤m≤1，k、e分别为第一、第二类完全椭圆积分且均与m有关。

进一步的，所述步骤(4)中根据发送线圈与接收线圈间的互感调整发送线圈上的电流，调整参数为：

其中：mti,r为第i个发送线圈与接收线圈间的互感，rl为接收线圈上的负载电阻，zl为接收线圈上的负载阻抗，mi为磁信道参数，βi为波束成形向量，*表示共轭；在多发单收无线能量传输系统中，通过磁波束成形计算出波束成形向量(β1,β2,...βi...βn)，继而调整发送线圈上的电流。

进一步的，所述固定频率为1mhz。

本发明的有益效果为：

本发明提出了基于近场定位的磁信道估计方法。先利用三个发送线圈上的互感系数的估计值估计接收线圈的位置，再利用估计得到的位置和其他发送线圈互感系数之间的关系获得接收线圈和其他发送线圈之间的互感系数。相比现有方法，该方法不需要接收线圈反馈电流值给发送线圈，也不需要逐一观测每个发送线圈上的电压和感应电流，简化了现有的磁信道估计方法。

附图说明

图1是本发明的系统框架示意图；

图2是本发明的线圈阵列示意图；

图3是5个发送线圈估计接收线圈位置的示意图；

图4是不同估计精度下，发送功率与接收功率的关系图。

具体实施方式

下面结合附图对本发明内容作进一步详细说明。

我们所考虑的无线能量传输系统如图1所示，是在近场下工作在谐振状态中的多发单收系统。如图2所示，左侧为发送线圈阵列，右侧为接收线圈。建立电路模型图时，为求简化分析，忽略线圈间的交叉耦合，只考虑两线圈间的直接耦合作用。从图上可以看出，发送线圈回路参数为l1、l2···ln，c1、c2···cn，r1、r2···rn。rl为系统等效负载，v1、v2···vn为激励电源。因为系统是工作在谐振状态下的，所以在实验中我们需要将发送与接收线圈调谐且工作频率相同为1mhz。

如图3所示，我们假设共有5个发送线圈，1个接收线圈。发送线圈的半径均为0.035m，线圈间的距离为0.085m。5个发送线圈的坐标分别为(0,0,0)，(8.5cm，0，0)，(0,8.5cm,0)，(8.5cm,8.5cm,0)，(0,17cm,0)。发送线圈上的电阻rt＝2.2ω，接收线圈上的电阻rr＝10ω。给每个发送线圈加载幅度5v、频率1mhz的激励信号。

1)闭合第n(1≤n≤3)个发送线圈，同时断开其他发送线圈。令ir表示接收线圈上的电流，rr表示接收线圈上的负载阻抗，表示第n个发送线圈与接收线圈间的互感，表示第n个发送线圈上的电流，rt表示发送线圈上的阻抗，vn为加载到第n个发送线圈上的总电压，ω表示施加在发送线圈上的信号的角频率。根据发送线圈和接收线圈上的kvl方程，可得：

我们对第n个发送线圈上的电压vn和电流itn进行l次观测。令其中，l表示第l次观测，yn(l)为第l次观测时实际加载到第n个发送线圈上的真实电压，vn(l)为第l次观测时加载到第n个发送线圈电路上的总电压，表示第l次观测时第n个发送线圈上的电流。则第l次观测时，式(1)可以表示为令表示第l次观测时实际加载到第n个发送线圈上的电压的观测值，即其中，en(l)表示观测误差。对观测值进行最小二乘估计，获得的最小二乘估计值即从获得存在相位模糊的问题。为了消除相位模糊，我们需要进一步观测其他断开的发送线圈上的感应电压。根据kvl方程，在第n'个发送线圈上的感应电压为

其中，表示第n个发送线圈和第n'个发送线圈之间的互感系数。通过(2)式可以看出对应的实部，因此可以通过感应电压vn'和电流的观测值确定其实部的符号。假设为正数，我们可通过上述方法确定其余互感系数的符号。

根据每个发送线圈的观测电流与观测电压值算出的如表1-3所示。xn(l)的单位是安培，观测电压的单位是伏特。l＝10。

表1根据发送线圈1上的电压电流值计算出的x1(l)，

表2根据发送线圈2上的电压电流值计算出的x2(l)，

表3根据发送线圈3上的电压电流值计算出的x3(l)，

根据表中数据，可以得到根据上面提到的去模糊方法，我们可以进一步得到至此，我们估计出了发送线圈与接收线圈之间的互感值。

2)接下来，我们将根据互感系数和接收线圈的位置之间的关系，估计接收线圈的位置参数。互感系数与磁场强度的关系为：

其中，μ0为空气的磁导率，为第n个发送线圈的匝数，nrx为接收线圈的匝数，arx为接收线圈的面积，an为第n个发送线圈的半径，δn和dn分别为接收线圈相对于第n个发送线圈的横向位移与纵向位移，k(mn)、e(mn)分别为与参数mn相关的第一、第二完全椭圆积分。

令(xr,yr,zr)表示接收线圈的位置，表示第n个发送线圈的位置。当接收线圈和发送线圈平行放置时，将dn和δn的表达式带入(3)式，在发送线圈的位置已知的条件下，互感系数只和接收线圈的位置相关。在1)中我们已经获得了互感系数的最小二乘估计值令

现有3个发送线圈，分别对应3个非线性方程。通过求解非线性方程组可以得到(xr,yr,zr)的估计值估计出接收线圈的位置坐标为接收线圈的实际位置坐标为接收线圈位置的估计值和真实值之间的误差为{ex＝0.0003,ey＝0.0004,ez＝0}。

由于和相关，因此，通过只能获得的估计值。因此，对于接收线圈的位置仍然存在模糊性，但这不影响对于磁信道的估计。

3)根据和其他发送线圈(比如第n'个发送线圈，n'＞3)的位置可以计算出发送线圈和接收线圈之间的相对位移(相对纵向位移dn'和相对横向位移δn')，再将dn'和δn'带入(3)式，便可得到通过上述方法，我们可以获得所有发送线圈和接收线圈之间的互感系数的估计值。

4)在得到了互感系数的估计值后，我们根据互感值调整发送线圈中的电流的幅度和相位，使接收线圈上的功率最大。令it表示所有发送线圈上的电流构成的电流向量，rti表示对角线元素为发送线圈电阻的对角阵(i为单位阵)。在发送功率一定的条件下，最大化接收功率的问题可以表示为

其中，s.t.表示约束条件，h表示矩阵的共轭转置，tr()表示求矩阵的迹，表示发送功率，令(4)转化为

其中，rank(si)表示矩阵si的秩，先去掉约束条件rank(si)＝1，并令则(5)进一步转化为

由于m为秩为1的对称阵，m的特征值分解可以表示为m＝λuu^h，其中，(||m||表示对m求2-范数)，表示矩阵m的特征值，表示特征值对应的特征向量。(6)的最优解为因此，(5)的最优解为si＝(rti+m)^-1由于si为秩为1的矩阵，(4)的最优解为

根据it中的值调整各个发送线圈中的电流的幅度和相位。

本发明的估计方法中，观测次数l的值和位置的估计精度相关，l越大，位置的估计精度越高，发送线圈的能量调整精度越高，能量传输效率也越高。图4为l取不同值时，发送功率和接收功率之间的关系图。

从图4中可以看出，在相同的估计精度下，发送功率与接收功率几乎呈线性关系；l＝1时，位置估计误差较大，能量传输效率较低，随着观测次数l的增大大，接收功率也增大，当l＝10时，位置估计精度较高，得到的接收功率与实际值十分接近，即位置的估计精度越高，能量传输效率也越高。