一种电动汽车直流充电畸变信号自适应识别方法与流程

本发明属于电动汽车充电技术领域，具体涉及一种充电过程中畸变信号的识别方法。

背景技术：

近年来，随着全球能源危机问题日益严重，可替代能源与能源可持续发展问题迫在眉睫，电能以其清洁的可再生能源优势成为众多国家的首选。目前，纯电动汽车已经开始产业化进入市场，配套的电动汽车充电机也大规模生产建设。

电动汽车直流充电机和蓄电池属于大功率非线性设备，导致电动汽车充电负荷含有大量的纹波，体现明显的非线性。尤其是充电站广泛应用的快速充电模式，采用大功率直流充电，密集的大规模集中充电导致瞬时负荷过大，时常出现骤变、非稳态信号，相较于常规交流充电模式，非线性特征更加突出，对于电动汽车充电电能的精确计量、充电故障的准确诊断以及畸变信号的防治产生不利影响。

针对电动汽车直流充电的特殊环境，急需针对充电负荷中存在的纹波、非稳态和骤变信号等畸变信号开展识别方法研究。目前本领域采用的傅里叶算法存在“吉布斯效应”，无法在时域上精确检测骤变和非稳态信号，且在骤变和非稳态信号的影响下，纹波的分解精度也会大幅下降。

技术实现要素：

发明目的：针对现有技术中存在的问题，本发明公开了一种电动汽车直流充电畸变信号自适应识别方法，该方法能够在电动汽车复杂的充电环境下同时精确检测各类畸变信号。

技术方案：本发明采用如下技术方案：

一种电动汽车直流充电畸变信号自适应识别方法，包括如下步骤：

(1)获取电动汽车充电的电信号f(t)，并进行等间隔采样，得到离散的电信号f(n)，采样频率为fs；

(2)对离散电信号f(n)进行小波包分解，获取小波包系数，具体包括如下步骤：

(2.1)小波包分解递推公式为：

其中，为小波包分解的第j层上第γ个小波包分解系数，j∈[1,m0],γ∈[1,2^j-1]，m0为小波包分解层数；h0(2n-z)、g0(2n-z)分别为进行小波包分解时的低通、高通滤波器；z取整数，是对连续小波包变换离散化的参数；

(2.2)对小波包分解系数进行修正，修正函数为：

其中，和分别为修正前和修正后的小波包系数，j∈[1,m0]；ε为修正阈值；n∈r为小波包比例系数；

(3)根据分解后的小波包系数判断骤变信号，并对骤变信号和非稳态信号进行重构；

判断直流充电信号中第j层小波包是否出现骤变信号和非稳态信号：

其中σ1为预设的信号突变阈值，△t＝1/fs为信号采样间隔，t0为当前离散序列时刻，t为信号突变判定所需的离散序列间隔。

如果上式成立，第j层小波包有骤变信号和非稳态信号，对骤变信号和非稳态信号进行重构：

其中，为重构后的小波包系数，h0(n-2z)和g0(n-2z)小波包重构的低通和高通滤波器组，与h0(2n-z)、g0(2n-z)无关；f骤(n)和f非稳态(n)为重构的骤变信号和非稳态信号；j骤，γ骤为骤变信号小波包的层次和所在层次的序号，位于小波包系数的频域最低层，即0hz层次；j非稳态，γ非稳态非稳态信号小波包的层次和所在层次的序号，位于小波包系数的频域最高层。

还包括根据分解后的小波包系数获取电动汽车直流充电信号的基波和纹波，具体包括如下步骤：

(a1)获取k时刻电动汽车充电的稳态电信号s(k)：

s(k)＝f(k)-f非稳态(k)-f骤(k)

(a2)建立基于变步长自适应卡尔曼滤波算法的矩阵模型识别基波和纹波信号，包括如下步骤：

(a2.1)构建的k时刻系统的状态向量:

x(k)＝[x01,x02,xi1,xi2…,xi1,xi2,…,xm1,xm2]^t

其中，基波平行分量x01＝f0(k△t)，正交分量x12＝0；xi1＝fi(k△t)cosθi，为第i次纹波电信号的平行分量；xi2＝-fi(k△t)sinθi，为第i次纹波电信号的正交分量；fi(k△t)第i次纹波随时间变化的幅值，△t＝1/fs为信号采样间隔，θi为充电信号第i次纹波初相位。

(a2.2)建立状态方程和观测方程：

其中w(k)和v(k)分别表示状态方程和观测方程的噪声；a(k)为状态转移矩阵，h(k)为观测矩阵；

(a2.3)进行自适应步长状态更新：

k时刻的步长μ(k)按下式计算，自适应修正(a2.4)的每一次迭代更新的k值：

其中s为遗忘因子，0<s<1；ρ为常数；是代价函数的瞬时估计值；dk为k时刻数据的弱化系数；tr(·)为矩阵求迹运算，ceil(·)为取最小近似整数函数。

(a2.4)更新系统噪声和量测噪声：

其中q(k)和q(k)分别为k时刻系统噪声的均值和方差；r(k)和r(k)分别为k时刻量测噪声的均值和方差；e(k)为观察误差；kg(k)为卡尔曼增益；p(k)为x(k)的方差；x(k|k-1)为通过k-1时刻的数据预测的k时刻系统的状态向量；p(k/k-1)为x(k/k-1)的方差。

(a2.5)计算基波和各次纹波信号的幅值和相位：

其中fi(k)为k时刻第i次纹波电信号的幅值；θi(k)为k时刻第i次纹波电信号的初相位。

所述步骤(2.1)中小波包分解中采用db40作为小波基函数。

所述步骤(2.1)中小波包分解层数m0的取值为：

其中，f0为各次纹波频率的最大公约数；ceil(·)为取最小近似整数函数。

本发明中，电动汽车充电的电信号f(t)为电压信号或电流信号。

状态转移矩阵a(k)为2(m+1)×2(m+1)的单位矩阵；

观测矩阵h(k)为2(m+1)×2(m+1)的方阵，取值为：

其中，f0为各次纹波频率的最大公约数，i＝0,1,3,5...,19，m＝10。

有益效果：与现有技术相比，本发明公开的电动汽车直流充电畸变信号自适应识别方法可以有效提取出非稳态信号以及骤变信号，相对于传统小波变换精确度更高；在消除前两类畸变信号的影响之后，采用变步长自适应卡尔曼滤波算法对基波和纹波建立相关矩阵模型，通过时变噪声的递推和预测估计器，随输入信号实时修正协方差参数，优化标准卡尔曼滤波算法的自适应性，遗忘因子、可变步长可以增加非稳态环境下信号识别的收敛速度；两种算法的组合应用能够解决现有单一算法难以同时精确检测各类畸变信号的问题，在电动汽车复杂的充电环境下，具有较高的正确性与适用性。

附图说明

图1为本发明公开的电动汽车直流充电畸变信号自适应识别方法的流程图。

具体实施方式

为使本发明的目的、技术方案和优点更加清楚，下面结合附图对本发明的具体实施案例做说明。

本发明公开了一种电动汽车直流充电畸变信号自适应识别方法，包括如下步骤：

(1)获取电动汽车充电的电信号f(t)，并进行等间隔采样，得到离散的电信号f(n)，采样频率为fs；

(2)对离散电信号f(n)进行小波包分解，获取小波包系数，具体包括如下步骤：

(2.1)小波包分解递推公式为：

其中，为小波包分解的第j层上第γ个小波包分解系数，j∈[1,m0],γ∈[1,2^j-1]，m0为小波包分解层数；h0(2n-z)、g0(2n-z)分别为进行小波包分解时的低通、高通滤波器；z取整数，是对连续小波包变换离散化的参数；

(2.2)对小波包分解系数进行修正，修正函数为：

其中，和分别为修正前和修正后的小波包系数，j∈[1,m0]；ε为修正阈值；n∈r为小波包比例系数，本发明中取1.12，可根据实际情况调整取值；

(3)根据分解后的小波包系数判断骤变信号，并对骤变信号和非稳态信号进行重构；

判断直流充电信号中第j层小波包是否出现骤变信号和非稳态信号：

其中σ1为预设的信号突变阈值，△t＝1/fs为信号采样间隔，t0为当前离散序列时刻，t为信号突变判定所需的离散序列间隔。

如果上式成立，第j层小波包有骤变信号和非稳态信号，对骤变信号和非稳态信号进行重构：

其中，为重构后的小波包系数，h0(n-2z)和g0(n-2z)小波包重构的低通和高通滤波器组，与h0(2n-z)、g0(2n-z)无关；f骤(n)和f非稳态(n)为重构的骤变信号和非稳态信号；j骤，γ骤为骤变信号小波包的层次和所在层次的序号，位于小波包系数的频域最低层，即0hz层次；j非稳态，γ非稳态非稳态信号小波包的层次和所在层次的序号，位于小波包系数的频域最高层。

还包括根据分解后的小波包系数获取电动汽车直流充电信号的基波和纹波，具体包括如下步骤：

(a1)获取k时刻电动汽车充电的稳态电信号s(k)：

s(k)＝f(k)-f非稳态(k)-f骤(k)

(a2)建立基于变步长自适应卡尔曼滤波算法的矩阵模型识别基波和纹波信号，包括如下步骤：

(a2.1)构建的k时刻系统的状态向量:

x(k)＝[x01,x02,xi1,xi2…,xi1,xi2,…,xm1,xm2]^t

其中，基波平行分量x01＝f0(k△t)，正交分量x12＝0；xi1＝fi(k△t)cosθi，为第i次纹波电信号的平行分量；xi2＝-fi(k△t)sinθi，为第i次纹波电信号的正交分量；fi(k△t)第i次纹波随时间变化的幅值，△t＝1/fs为信号采样间隔，θi为充电信号第i次纹波初相位。

(a2.2)建立状态方程和观测方程：

其中w(k)和v(k)分别表示状态方程和观测方程的噪声；a(k)为状态转移矩阵，h(k)为观测矩阵；

(a2.3)进行自适应步长状态更新：

k时刻的步长μ(k)按下式计算，自适应修正(a2.4)的每一次迭代更新的k值：

其中s为遗忘因子，0<s<1；ρ为常数；是代价函数的瞬时估计值；dk为k时刻数据的弱化系数；tr(·)为矩阵求迹运算，ceil(·)为取最小近似整数函数。

(a2.4)更新系统噪声和量测噪声：

其中q(k)和q(k)分别为k时刻系统噪声的均值和方差；r(k)和r(k)分别为k时刻量测噪声的均值和方差；e(k)为观察误差；kg(k)为卡尔曼增益；p(k)为x(k)的方差；x(k|k-1)为通过k-1时刻的数据预测的k时刻系统的状态向量；p(k/k-1)为x(k/k-1)的方差。

(a2.5)计算基波和各次纹波信号的幅值和相位：

其中fi(k)为k时刻第i次纹波电信号的幅值；θi(k)为k时刻第i次纹波电信号的初相位。

所述步骤(2.1)中小波包分解中采用db40作为小波基函数。

所述步骤(2.1)中小波包分解层数m0的取值为：

其中，f0为各次纹波频率的最大公约数；ceil(·)为取最小近似整数函数。

本发明中，电动汽车充电的电信号f(t)为电压信号或电流信号。

状态转移矩阵a(k)为2(m+1)×2(m+1)的单位矩阵；

观测矩阵h(k)为2(m+1)×2(m+1)的方阵，取值为：

其中，f0为各次纹波频率的最大公约数，i＝0,1,3,5...,19，m＝10。