CN114915225A - 一种基于优化bp神经网络的永磁同步电机参数辨识方法 - Google Patents

Abstract

本发明涉及一种基于优化BP神经网络的永磁同步电机参数辨识方法，具体包括以下步骤：(1)搭建FOC控制系统控制已知参数永磁同步电机运行；(2)采集电机运行信号；(3)坐标变换；(4)PI控制；(5)改进鲸鱼优化算法优化BP神经网络训练；(6)实现电机参数识别。本发明提出的永磁同步电机参数辨识方法引入了一种智能寻优算法与神经网络的结合算法，相比传统的参数辨识方法，其优势是收敛速度快，无需进行公式推导和计算，模型搭建简便快捷，鲁棒性强，辨识精度高，辨识结果波动小。

Description

一种基于优化BP神经网络的永磁同步电机参数辨识方法

技术领域

本发明属于电机参数辨识技术领域，尤其涉及一种基于优化BP神经网络的永磁同步电机参数辨识方法。

背景技术

永磁同步电机因具有高效率、高转矩电流比、高功率密度、可靠性等优点，已在电动汽车、航空工业等领域得到了广泛应用。在永磁同步电机的控制过程中，闭环控制、减小运行损耗方法以及无传感器控制等方面，都比较依赖永磁同步电机自身参数，输入不准确参数的控制会导致永磁同步电机闭环控制效果差，系统波动大；引起减小运行损耗方法效果不佳导致电机损耗增大；无传感器控制中如果电机参数输入不准确，会引起电机运行不畅甚至无法运行。因此高精度的永磁同步电机参数辨识方法是很有必要的。

目前多采用参数辨识方法包括最小二乘法、模型参考自适应法、神经网络法等。其中，最小二乘法的参数辨识无法实现多参数同时辨识，在辨识某一参数时，其他参数是假定为已知值，因此最小二乘法参数辨识的辨识效果较差，在对电机参数有较高要求的领域难以胜任；模型参考自适应法一定程度上解决了最小二乘法无法同时实现多参数辨识的问题，且参数辨识精度更高，但是需要进行大量的公式推导和计算，同时还需要人工调节模型参考自适应的自适应律以及搭建精准的电机控制参考模型和可调模型，使用不够便捷；神经网络法无需进行复杂的公式推导及计算，也无需搭建精准的电机模型，使用方便，且电机参数辨识结果远好于其他方法，但是存在辨识初始阶段波动较大，鲁棒性较差，辨识速度慢等问题。

发明内容

本发明所要解决的技术问题是提供一种基于改进鲸鱼优化算法优化BP神经网络的永磁同步电机参数辨识方法，以解决现有技术中存在的无法进行多参数辨识、需要进行复杂的公式推导及运算、人工调节自适应律、需要搭建精确模型等问题。

为解决以上技术问题，本发明的技术方案为：提供一种基于优化BP神经网络的永磁同步电机参数辨识方法，其创新点在于：具体包括以下步骤：

(1)搭建FOC控制系统控制已知参数永磁同步电机运行：根据永磁同步电机FOC矢量控制方法搭建FOC控制系统，并控制电机在小于该电机最大转矩、最大转速的工况下运行；

(2)采集电机运行信号：通过电流传感器检测永磁同步电机的三相电流i_a、i_b、i_c；通过转速传感器检测电机的机械角速度ω_m，通过计算得到电机电角速度ω_e；通过位置传感器检测电机机械角度θ_m，通过计算得到电机电角度θ_e；

(3)坐标变换：将步骤(2)得到的电机的三相电流i_a、i_b、i_c通过Clarke变换转化为静止坐标系下的电流i_α、i_β，再通过Park变换转化为同步旋转坐标系下的电流i_d、i_q；

(4)PI控制：将步骤(2)得到的电角速度ω_e和参考转速ω_{e_ref}输入速度环PI模块，速度环PI模块输出q轴电流参考值

至q轴电流环PI模块；将参考电流i_{d_ref}输入d轴电流环PI模块，同时将步骤(3)得到的同步旋转坐标系下的电流i_d、i_q分别输入d轴电流环PI模块以及q轴电流环PI模块进行计算，d轴电流环PI模块以及q轴电流环PI模块分别输出d轴电压u_d和q轴电压u_q；

(5)改进鲸鱼优化算法优化BP神经网络训练：改进鲸鱼优化算法优化BP神经网络TCWOA-BPNN，BP神经网络TCWOA-BPNN由输入层、隐含层及输出层构成，使用改进的鲸鱼优化算法获取BP神经网络的初始权值和阈值，将步骤(3)所得到的d、q轴电流i_d、i_q、步骤(4)得到的d、q轴电压u_d、u_q以及步骤(2)得到的电机电角速度ω_e作为TCWOA-BPNN的输入，将步骤(1)中所述的已知参数的永磁同步电机参数作为输出，进行训练，获得用于永磁同步电机参数辨识的BP神经网络TCWOA-BPNN；

(6)实现电机参数识别：将未知参数的永磁同步电机使用FOC控制系统控制运行，重复步骤(2)、步骤(3)和步骤(4)，将得到的数据i_d、i_q、u_d、u_q、ω_e输入至步骤(5)得出的BP神经网络TCWOA-BPNN中，即可输出永磁同步电机PMSM的各项参数，则完成对永磁同步电机的参数识别。

进一步的，所述步骤(1)中根据永磁同步电机FOC矢量控制方法搭建FOC控制系统包括速度环PI模块、q轴电流环PI模块、d轴电流环PI模块、Clarke变换模块、Park变换模块，SVPWM模块、逆Park变换模块，逆变器、永磁同步电机PMSM模块、电流传感器、速度/位置传感器、TCWOA-BPNN参数辨识模块；还包括需要给定速度环PI模块的参考转速ω_{e_ref}以及d轴电流环PI模块的参考电流i_{d_ref}；

所述逆变器和永磁同步电机PMSM模块连接来驱动永磁同步电机PMSM模块运行，所述电流传感器模块设于逆变器和永磁同步电机PMSM模块之间来采集逆变器输出到永磁同步电机PMSM模块的三相电流i_a、i_d、i_c并输出到Clark变换模块进行运算；

所述Clark变换模块分别和Park变换模块、速度/位置传感器连接，用于将三相电流i_a、i_d、i_c通过Clark变换模块得到的静止坐标系下的电流i_α、i_β分别输入Park变换模块、速度/位置传感器中；

所述速度/位置传感器分别与速度环PI模块、Park变换模块、逆Park变换模块、永磁同步电机PMSM模块连接，从永磁同步电机PMSM模块采集永磁同步电机机械角速度ω_m并计算得到电角速度ω_e，输入速度环PI模块对永磁同步电机转速进行控制；从永磁同步电机PMSM模块采集永磁同步电机机械角度θ_m并计算得到电角度θ_e，分别输入Park变换模块、逆Park变换模块参与计算；

所述速度环PI模块将永磁同步电机电角速度ω_e以及参考转速ω_{e_ref}进行运算得到q轴电流参考值i_q*，输入到q轴电流环PI模块；

所述Park变换模块分别与q轴电流环PI模块、d轴电流环PI模块连接；所述Park变换模块将静止坐标系下的电流i_α、i_β以及速度/位置传感器得到的电角度θ_e输入进行运算转换为同步旋转坐标系下的电流i_d、i_q；将i_d、i_q分别输入到d轴电流环PI模块以及q轴电流环PI模块，完成永磁同步电机d、q轴电流控制；

所述q轴电流环PI模块将q轴电流参考值i_q*以及i_q进行运算得到u_q，所述d轴电流环PI模块将输入的参考电流i_{d_ref}以及i_d进行运算得到u_d；

所述逆Park变换模块分别与q轴电流环PI模块、d轴电流环PI模块、SVPWM模块连接，用于将从q轴电流环PI模块和d轴电流环PI模块输出的u_q、u_d以及速度/位置传感器得到的电角度θ_e输入进行运算，得到静止坐标系下的电压u_α、u_β并输入SVPWM模块；

所述SVPWM模块和逆变器模块连接，用于将静止坐标系下的电压u_α、u_β输入到SVPWM模块运算，得到三相电压u_a、u_b、u_c，并输入到逆变器模块；

所述TCWOA-BPNN参数辨识模块与Park变换模块、d轴电流环PI模块、q轴电流环PI模块以及速度/位置传感器模块连接，将i_d、i_q、u_d、u_q、ω_e输入TCWOA-BPNN参数辨识模块中，TCWOA-BPNN参数辨识模块计算输出永磁同步电机的线电阻R_s、d轴电感L_d、q轴电感L_q以及磁链ψ_f。

进一步的，所述逆Park变换描述如下：

通过逆Park变换将同步旋转坐标系下的电压u_d、u_q转化为静止坐标系下的电压u_α、u_β，具体方法为：

u_α＝u_dcosθ_e-u_qsinθ_e

u_β＝u_dsinθ_e+u_qcosθ_e。

进一步的，所述步骤(1)中已知参数永磁同步电机所需要已知的参数包括同步旋转坐标系下的d轴电感L_d、q轴电感L_q、永磁同步电机极对数p、永磁同步电机线电阻R_s、永磁同步电机磁链ψ_f、永磁同步电机的转动惯量J、阻尼系数B。

进一步的，所述步骤(2)中转速传感器检测电机的机械角速度ω_m，通过计算得到的电机电角速度ω_e的具体方法为：

ω_e＝p*ω_m

通过位置传感器检测电机机械角度θ_m，并通过计算得到的电机电角度θ_e，其具体方法为：

θ_e＝p*θ_m

其中，p表示永磁同步电机PMSM极对数。

进一步的，所述步骤(3)中，Clarke变换、Park变换描述如下：

将电机的三相电流i_a、i_b、i_c通过Clarke变换转化为静止坐标系下的电流i_α、i_β的具体方法为：

将静止坐标系下的电流i_α、i_β通过Park变换转化为同步旋转坐标系下的电流i_d、i_q的具体方法为：

i_d＝i_αcosθ_e+i_βsinθ_e

i_q＝-i_αsinθ_e+i_βcosθ_e。

进一步的，所述步骤(5)中，改进鲸鱼优化算法优化BP神经网络训练TCWOA-BPNN，使用改进的鲸鱼优化算法获取BP神经网络的初始权值和阈值具体步骤如下：

A、鲸鱼优化算法的改进包括Tent混沌映射以及精英学习机制，Tent混沌映射表示为：

式中x₀为鲸鱼优化算法种群中的初始个体，x_n和x_n+1为种群中生成的第n个个体和第n+1个个体，rand[0,1]表示从0到1之间取随机数，精英学习机制表示为：

X_C＝j·k·[X_i-X^*]

式中X_C为鲸鱼优化算法生成的精英个体，j和k为随机参数值，X_i为当前鲸鱼优化算法种群中的个体，X^*为种群中当前位置最优个；

B、由于改进鲸鱼优化算法用于找到BP神经网络训练结果偏差最小的值，自变量数量为1，则设置改进鲸鱼优化算法的自变量维度为1，种群规模为N，最大迭代次数设置为T次；

C、设置概率p₁＝0.5为种群会按当前最优值继续寻优，设置概率p₂＝0.5为种群会从头开始寻优；

D、设置改进鲸鱼优化算法再得到的BP神经网络训练误差小于等于2％或迭代达到最大迭代次数时停止，改进鲸鱼优化算法停止工作后，得到的BP神经网络训练结果偏差最小的值对应的初始权值j₀和阈值k₀，将初始权值j₀和阈值k₀赋予BP神经网络作为初始权值和阈值进行训练。

进一步的，所述步骤(5)中，改进鲸鱼优化算法优化BP神经网络训练，BP神经网络的参数设置如下：

A、BP神经网络的输入数据集为{i_d、i_q、u_d、u_q、ω_e}，则设置BP神经网络的输入数为5；输出数据集为{R_s、L_d、L_q、ψ_f}，则设置输出数为4；

B、BP神经网络的隐含层数量选择m层，隐含层神经元数量根据经验及实践确定为n个。

进一步的，所述步骤(6)中，输出未知参数的永磁同步电机各项参数，具体可输出的参数包括：永磁同步电机线电阻R_s、d轴电感L_d、q轴电感L_q以及磁链ψ_f。

与现有技术相比，本发明具有以下优点：

1.本发明提出的永磁同步电机参数辨识方法通过引入的一种改进的鲸鱼优化算法与BP神经网络的相结合，相比传统的参数辨识方法，其优势是收敛速度快，无需进行公式推导和计算，模型搭建简便快捷，鲁棒性强，辨识精度高，辨识结果波动小。

2.本发明提出的一种基于改进鲸鱼优化算法优化BP神经网络的永磁同步电机参数辨识方法，相比最小二乘参数辨识法，能够直接进行多个参数的辨识，且辨识精度比最小二乘法高。

3.本发明提出的一种基于改进鲸鱼优化算法优化BP神经网络的永磁同步电机参数辨识方法，相比模型参考自适应参数辨识法，辨识精度更高，且无需进行大量的公式推导和计算，同时省略了人工调节模型参考自适应的自适应律以及搭建精准的电机控制参考模型和可调模型的步骤，使用起来更加方便。

4.本发明提出的一种基于改进鲸鱼优化算法优化BP神经网络的永磁同步电机参数辨识方法，相比传统神经网络参数辨识法，辨识初始阶段波动小，鲁棒性较好，辨识速度更快。

附图说明

以下附图是本发明利用FOC控制系统所得出的部分仿真结果图，显而易见的，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图，以下实施例为了更直观表现本发明对永磁同步电机参数辨识的效果，加入了未作改进的鲸鱼优化算法优化BP神经网络(WOA-BPNN)以及未做优化的BP神经网络(BPNN)两种辨识法的对比，在正常运用时，无需进行对比。

图1是一种高精度的永磁同步电机的参数辨识方法流程图。

图2是本发明的基于改进鲸鱼优化算法优化BP神经网络的永磁同步电机参数辨识方法的结构框图。

图3为永磁同步电机电阻、d轴、q轴电感以及磁链变化情况。

图4是给定工况下，永磁同步电机线电阻R_s基于TCWOA-BPNN、WOA-BPNN与BPNN的辨识结果对比。

图5是给定工况下，永磁同步电机d轴电感L_d基于TCWOA-BPNN、WOA-BPNN与BPNN的辨识结果对比。

图6是给定工况下，永磁同步电机q轴电感L_q基于TCWOA-BPNN、WOA-BPNN与BPNN的辨识结果对比。

图7是给定工况下，永磁同步电机磁链ψ_f基于TCWOA-BPNN、WOA-BPNN与BPNN的辨识结果对比。

具体实施方式

为了说明本发明所述的技术方案，结合附图进一步说明本发明。

本发明提供基于改进鲸鱼优化算法优化BP神经网络的永磁同步电机参数辨识方法，如图1和2所示，包括以下步骤：

(1)搭建FOC控制系统控制已知参数永磁同步电机运行：根据永磁同步电机FOC矢量控制方法搭建FOC控制系统，并控制电机在小于本电机最大转矩、最大转速的工况下运行；在本发明的实施方式中，本发明使用Matlab/Simulink软件进行FOC控制系统搭建，通过所使用的永磁同步电机的型号获取永磁同步电机同步旋转坐标系下的d轴电感L_d、q轴电感L_q、永磁同步电机极对数p、永磁同步电机线电阻R_s、永磁同步电机磁链ψ_f、永磁同步电机的转动惯量J、阻尼系数B。使用Matlab/Simulink软件搭建速度环PI模块、q轴电流环PI模块、d轴电流环PI模块、Clarke变换模块、Park变换模块，SVPWM模块、逆Park变换模块，逆变器模块、永磁同步电机PMSM模块、电流传感器、速度/位置传感器、WOA-BPNN参数辨识模块，并给定速度环PI模块的参考转速ω_{e_ref}以及d轴电流环PI模块的参考电流i_{d_ref}，其中Clarke变换模块、Park变换模块，SVPWM模块、逆Park变换模块、永磁同步电机PMSM模块以及WOA-BPNN参数辨识模块使用Matlab/Simulink中的S-function进行程序编写，其余模块则是使用Matlab/Simulink模块库自带的模块；使用Matlab/Simulink软件搭建的FOC控制系统的具体如下：

逆变器模块和永磁同步电机PMSM模块连接来驱动永磁同步电机PMSM模块运行，电流传感器模块设于逆变器和永磁同步电机PMSM模块之间采集逆变器输出到永磁同步电机PMSM模块的三相电流i_a、i_b、i_c并输出到Clark变换模块进行运算；

Clark变换模块分别和Park变换模块、速度/位置传感器连接，将三相电流i_a、i_b、i_c通过Clark变换模块得到的静止坐标系下的电流i_α、i_β分别输入Park变换模块、速度/位置传感器中；

速度/位置传感器分别与速度环PI模块、Park变换模块、逆Park变换模块、永磁同步电机PMSM模块连接，从永磁同步电机PMSM模块采集永磁同步电机机械角速度ω_m并计算得到电角速度ω_e，输入速度环PI模块对永磁同步电机转速进行控制；从永磁同步电机PMSM模块采集永磁同步电机机械角度θ_m并计算得到电角度θ_e，分别输入Park变换模块、逆Park变换模块参与计算；

速度环PI模块将永磁同步电机电角速度ω_e以及参考转速ω_{e_ref}进行运算得到q轴电流参考值i_q*，输入到q轴电流环PI模块；

Park变换模块分别与q轴电流环PI模块、d轴电流环PI模块连接；Park变换模块将静止坐标系下的电流i_α、i_β以及速度/位置传感器得到的电角度θ_e输入进行运算转换为同步旋转坐标系下的电流i_d、i_q；将i_d、i_q分别输入到d轴电流环PI模块以及q轴电流环PI模块，完成永磁同步电机d、q轴电流控制；

q轴电流环PI模块将q轴电流参考值i_q*和i_q进行运算得到u_q，d轴电流环PI模块将输入的参考电流i_{d_ref}和i_d进行运算得到u_d；

逆Park变换模块分别与q轴电流环PI模块、d轴电流环PI模块、SVPWM模块连接，将分别从q轴电流环PI模块和d轴电流环PI模块得到的u_q、u_d以及速度/位置传感器得到的电角度θ_e输入进行运算，得到静止坐标系下的电压u_α、u_β并输入SVPWM模块；

SVPWM模块和逆变器模块连接，用于将静止坐标系下的电压u_α、u_β输入到SVPWM模块运算，得到三相电压u_a、u_b、u_c，并输入到逆变器模块；

WOA-BPNN参数辨识模块与Park变换模块、d轴电流环PI模块、q轴电流环PI模块以及速度/位置传感器模块连接，将i_d、i_q、u_d、u_q、ω_e输入WOA-BPNN参数辨识模块中，计算输出永磁同步电机的线电阻R_s、d轴电感L_d、q轴电感L_q以及磁链ψ_f。

通过逆Park变换将同步旋转坐标系下的电压u_d、u_q转化为静止坐标系下的电压u_α、u_β，具体方法为：

u_α＝u_dcosθ_e-u_qsinθ_e

u_β＝u_dsinθ_e+u_qcosθ_e。

(2)采集电机运行信号：通过电流传感器检测永磁同步电机的三相电流i_a、i_b、i_c；通过转速传感器检测电机的机械角速度ω_m，通过计算得到电机电角速度ω_e；通过位置传感器检测电机机械角度θ_m，通过计算得到电机电角度θ_e；

转速传感器检测电机的机械角速度ω_m，并通过计算得到的电机电角速度ω_e的具体方法为：

ω_e＝p*ω_m

位置传感器检测电机机械角度θ_m，并通过计算得到的电机电角度θ_e，其具体方法为：

θ_e＝p*θ_m

其中，p表示极对数。

(3)坐标变换：将步骤(2)得到的电机的三相电流i_a、i_b、i_c通过Clarke变换转化为静止坐标系下的电流i_α、i_β，再通过Park变换转化为同步旋转坐标系下的电流i_d、i_q；

将电机的三相电流i_a、i_b、i_c通过Clarke变换转化为静止坐标系下的电流i_α、i_β的具体方法为：

将静止坐标系下的电流i_α、i_β通过Park变换转化为同步旋转坐标系下的电流i_d、i_q的具体方法为：

i_d＝i_acosθ_e+i_βsinθ_e

i_q＝-i_αsinθ_e+i_βcosθ_e。

(4)PI控制：使用Matlab/Simulink中的PID模块搭建速度环PI模块、d轴电流环PI模块以及q轴电流环PI模块；将步骤(2)得到的电角速度ω_e和参考转速ω_{e_ref}输入速度环PI模块，速度环PI模块输出q轴电流参考值

至q轴电流环PI模块；将参考电流i_{d_ref}输入d轴电流环PI模块，同时将步骤(3)得到的同步旋转坐标系下的电流i_d、i_q分别输入d轴电流环PI模块以及q轴电流环PI模块进行计算，d轴电流环PI模块以及q轴电流环PI模块分别输出d轴电压u_d和q轴电压u_q；

(5)改进鲸鱼优化算法优化BP神经网络训练：改进鲸鱼优化算法优化BP神经网络TCWOA-BPNN，由输入层、隐含层及输出层构成，将步骤(3)所得到的d、q轴电流i_d、i_q、步骤(4)得到的d、q轴电压u_d、u_q以及步骤(2)得到的电机电角速度ω_e作为TCWOA-BPNN的输入，将步骤(1)中所述的已知参数的永磁同步电机参数作为输出，进行训练，获得用于永磁同步电机参数辨识的改进鲸鱼优化算法优化BP神经网络；其中改进鲸鱼优化算法的参数设置如下：

A、鲸鱼优化算法的改进包括Tent混沌映射(Tent chaotic map)以及精英学习机制(Cream learning mechanism)。Tent混沌映射表示为：

式中x₀为鲸鱼优化算法种群中的初始个体，x_n和x_n+1为种群中生成的第n个个体和第n+1个个体，rand[0,1]表示从0到1之间取随机数。精英学习机制表示为：

X_C＝j·k·[X_i-X^*]

式中X_C为鲸鱼优化算法生成的精英个体，j和k的取值为0-1之间的随机数，X_i为当前鲸鱼优化算法种群中的个体，X^*为种群中当前位置最优个体。

B、改进鲸鱼优化算法在于找到BP神经网络训练结果偏差最小的值，自变量数量为1，则设置改的鲸鱼优化算法的自变量维度等于1，设种群规模为N，最大迭代次数设置为T次，本实施方式中的N取值50，T取值30；

C、设置概率p₁＝0.5为种群会按当前最优值继续寻优，设置概率p₂＝0.5为种群会从头开始寻优；

D、设置改进鲸鱼优化算法再得到的BP神经网络训练误差小于等于2％或迭代达到最大迭代次数时停止，改进鲸鱼优化算法停止工作后，将得到的BP神经网络训练结果偏差最小的值对应的初始权值j₀和阈值k₀赋予BP神经网络作为初始权值和阈值进行训练。

BP神经网络的参数设置如下：

A、BP神经网络的输入数据集为{i_d、i_q、u_d、u_q、ω_e}，因此设置BP神经网络的输入数为5；输出数据集为{R_s、L_d、L_q、ψ_f}，因此设置输出数为4；

B、BP神经网络的隐含层数量设置为m层，隐含层神经元数量选择n个，本实施方式中m取值为1，n取值为6。

(6)实现电机参数识别：将未知参数的永磁同步电机使用FOC控制系统控制运行，重复步骤(2)、步骤(3)和步骤(4)，将得到的数据i_d、i_q、u_d、u_q、ω_e输入至步骤(5)得出的改进鲸鱼优化算法优化BP神经网络中，输出永磁同步电机PMSM的线电阻R_s、d轴电感L_d、q轴电感L_q以及磁链ψ_f。

使用Matlab/Simulink仿真软件运行4s，为了验证本发明所提出的一种基于鲸鱼优化算法优化BP神经网络的永磁同步电机参数辨识方法的鲁棒性，所有实施例在运行的前两秒的负载T_L＝1N·m，后两秒的负载T_L＝3N·m；此外，实施例为验证本发明所述的参数辨识方法应对电机参数发生变化时的可靠性，设置电机参数随时间发生变化，其中线电阻R_s由初始值0.180Ω，逐渐上升至0.209Ω；d轴电感L_d由初始值0.0893H，逐渐上升至0.103H；q轴电感L_q由初始值0.205H，逐渐上升至0.236H；磁链ψ_f由初始值0.0160Wb，逐渐下降至0.0136H，参数具体变化情况如图3所示；此外，电机转动惯量J＝0.003kg·m²,阻尼系数B＝0.008，极对数p＝2；所有实施例运行时都设置参考转速为1000r/min，d轴参考电流i_{d_ref}＝0。仿真结果图如图4、5、6、7所示。

为兼顾辨识的速度和精度，本实施例中，输入数据集{i_d、i_q、u_d、u_q、ω_e}中的每个元素进行了随机采样，共采样1000个数据样本，取600个样本进行TCWOA-BPNN训练，剩下400个样本作为测试。

根据附图4对比BPNN、WOA-BPNN以及TCWOA-BPNN三个参数辨识方法对永磁同步电机定子电阻R_s的辨识效果，在辨识初期以及负载发生变化时，BPNN的辨识误差极大，WOA-BPNN与TCWOA-BPNN则有较好的辨识精度，在负载发生变化时，BPNN产生的波动约为0.035Ω，WOA-BPNN产生的波动约为0.02Ω，TCWOA-BPNN产生的波动约为0.01Ω，并且从图中可以看出在稳态时，TCWOA-BPNN的定子电阻稳态辨识误差最小。

根据附图5对比BPNN、WOA-BPNN以及TCWOA-BPNN三个参数辨识方法对永磁同步电机d轴电感L_d的辨识效果，在辨识初期以及负载发生变化时，BPNN的辨识误差极大，WOA-BPNN与TCWOA-BPNN则有较好的辨识精度，在负载发生变化时，BPNN产生的波动约为0.01H，WOA-BPNN产生的波动约为0.005Ω，而TCWOA-BPNN则几乎未产生波动，并且从图中可以看出在稳态时，TCWOA-BPNN的d轴电感稳态辨识误差最小。

根据附图6对比BPNN、WOA-BPNN以及TCWOA-BPNN三个参数辨识方法对永磁同步电机q轴电感L_q的辨识效果，在辨识初期以及负载发生变化时，BPNN的辨识误差极大，WOA-BPNN与TCWOA-BPNN则有较好的辨识精度，在负载发生变化时，BPNN产生的波动约为0.05H，WOA-BPNN产生的波动约为0.03H，而TCWOA-BPNN产生的波动约为0.02H，并且从图中可以看出在稳态时，TCWOA-BPNN的q轴电感稳态辨识误差最小。

根据附图7对比BPNN、WOA-BPNN以及TCWOA-BPNN三个参数辨识方法对永磁同步电机磁链ψ_f的辨识效果，在辨识初期以及负载发生变化时，BPNN的辨识误差极大，WOA-BPNN与TCWOA-BPNN则有较好的辨识精度，在负载发生变化时，BPNN产生的波动约为0.0035Wb，WOA-BPNN产生的波动约为0.002Wb，而TCWOA-BPNN几乎不存在波动，并且从图中可以看出在稳态时，TCWOA-BPNN的磁链稳态辨识误差最小。

综合对附图4、5、6、7的对比分析可得，基于TCWOA-BPNN永磁同步电机参数辨识方法的辨识效果相较于BP神经网络(BPNN)以及未作改进的鲸鱼优化算法优化BP神经网络(WOA-BPNN)精度更高，辨识过程中波动更小，在辨识初期就有较好的辨识效果，并且在负载发生突变时，产生的辨识波动也远小于其他两种方法。因此可以验证本发明提出的一种基于改进鲸鱼优化算法优化BP神经网络的永磁同步电机参数辨识方法，辨识精度高，辨识系统的鲁棒性好。

Claims (9)

1.一种基于优化BP神经网络的永磁同步电机参数辨识方法，其特征在于：具体包括以下步骤：

(1)搭建FOC控制系统控制已知参数永磁同步电机运行：根据永磁同步电机FOC矢量控制方法搭建FOC控制系统，并控制电机在小于该电机最大转矩、最大转速的工况下运行；

(2)采集电机运行信号：通过电流传感器检测永磁同步电机的三相电流i_a、i_b、i_c；通过转速传感器检测电机的机械角速度ω_m，通过计算得到电机电角速度ω_e；通过位置传感器检测电机机械角度θ_m，通过计算得到电机电角度θ_e；

(3)坐标变换：将步骤(2)得到的电机的三相电流i_a、i_b、i_c通过Clarke变换转化为静止坐标系下的电流i_α、i_β，再通过Park变换转化为同步旋转坐标系下的电流i_d、i_q；

(4)PI控制：将步骤(2)得到的电角速度ω_e和参考转速ω_{e_ref}输入速度环PI模块，速度环PI模块输出q轴电流参考值

至q轴电流环PI模块；将参考电流i_{d_ref}输入d轴电流环PI模块，同时将步骤(3)得到的同步旋转坐标系下的电流i_d、i_q分别输入d轴电流环PI模块以及q轴电流环PI模块进行计算，d轴电流环PI模块以及q轴电流环PI模块分别输出d轴电压u_d和q轴电压u_q；

(5)改进鲸鱼优化算法优化BP神经网络训练：改进鲸鱼优化算法优化BP神经网络TCWOA-BPNN，BP神经网络TCWOA-BPNN由输入层、隐含层及输出层构成，使用改进的鲸鱼优化算法获取BP神经网络的初始权值和阈值，将步骤(3)所得到的d、q轴电流i_d、i_q、步骤(4)得到的d、q轴电压u_d、u_q以及步骤(2)得到的电机电角速度ω_e作为TCWOA-BPNN的输入，将步骤(1)中所述的已知参数的永磁同步电机参数作为输出，进行训练，获得用于永磁同步电机参数辨识的BP神经网络TCWOA-BPNN；

(6)实现电机参数识别：将未知参数的永磁同步电机使用FOC控制系统控制运行，重复步骤(2)、步骤(3)和步骤(4)，将得到的数据i_d、i_q、u_d、u_q、ω_e输入至步骤(5)得出的BP神经网络TCWOA-BPNN中，即可输出永磁同步电机PMSM的各项参数，则完成对永磁同步电机的参数识别。

2.根据权利要求1所述的一种基于优化BP神经网络的永磁同步电机参数辨识方法，其特征在于：所述步骤(1)中根据永磁同步电机FOC矢量控制方法搭建FOC控制系统包括速度环PI模块、q轴电流环PI模块、d轴电流环PI模块、Clarke变换模块、Park变换模块，SVPWM模块、逆Park变换模块，逆变器、永磁同步电机PMSM模块、电流传感器、速度/位置传感器、TCWOA-BPNN参数辨识模块；还包括需要给定速度环PI模块的参考转速ω_{e_ref}以及d轴电流环PI模块的参考电流i_{d_ref}；

所述逆变器和永磁同步电机PMSM模块连接来驱动永磁同步电机PMSM模块运行，所述电流传感器模块设于逆变器和永磁同步电机PMSM模块之间来采集逆变器输出到永磁同步电机PMSM模块的三相电流i_a、i_b、i_c并输出到Clark变换模块进行运算；

所述Clark变换模块分别和Park变换模块、速度/位置传感器连接，用于将三相电流i_a、i_b、i_c通过Clark变换模块得到的静止坐标系下的电流i_α、i_β分别输入Park变换模块、速度/位置传感器中；

所述速度/位置传感器分别与速度环PI模块、Park变换模块、逆Park变换模块、永磁同步电机PMSM模块连接，从永磁同步电机PMSM模块采集永磁同步电机机械角速度ω_m并计算得到电角速度ω_e，输入速度环PI模块对永磁同步电机转速进行控制；从永磁同步电机PMSM模块采集永磁同步电机机械角度θ_m并计算得到电角度θ_e，分别输入Park变换模块、逆Park变换模块参与计算；

所述速度环PI模块将永磁同步电机电角速度ω_e以及参考转速ω_{e_ref}进行运算得到q轴电流参考值i_q*，输入到q轴电流环PI模块；

所述Park变换模块分别与q轴电流环PI模块、d轴电流环PI模块连接；所述Park变换模块将静止坐标系下的电流i_α、i_β以及速度/位置传感器得到的电角度θ_e输入进行运算转换为同步旋转坐标系下的电流i_d、i_q；将i_d、i_q分别输入到d轴电流环PI模块以及q轴电流环PI模块，完成永磁同步电机d、q轴电流控制；

所述q轴电流环PI模块将q轴电流参考值i_q*以及i_q进行运算得到u_q，所述d轴电流环PI模块将输入的参考电流i_{d_ref}以及i_d进行运算得到u_d；

所述逆Park变换模块分别与q轴电流环PI模块、d轴电流环PI模块、SVPWM模块连接，用于将从q轴电流环PI模块和d轴电流环PI模块输出的u_q、u_d以及速度/位置传感器得到的电角度θ_e输入进行运算，得到静止坐标系下的电压u_α、u_β并输入SVPWM模块；

所述SVPWM模块和逆变器模块连接，用于将静止坐标系下的电压u_α、u_β输入到SVPWM模块运算，得到三相电压u_a、u_b、u_c，并输入到逆变器模块；

所述TCWOA-BPNN参数辨识模块与Park变换模块、d轴电流环PI模块、q轴电流环PI模块以及速度/位置传感器模块连接，将i_d、i_q、u_d、u_q、ω_e输入TCWOA-BPNN参数辨识模块中，TCWOA-BPNN参数辨识模块计算输出永磁同步电机的线电阻R_s、d轴电感L_d、q轴电感L_q以及磁链ψ_f。

3.根据权利要求2所述的一种基于优化BP神经网络的永磁同步电机参数辨识方法，其特征在于：所述逆Park变换描述如下：

通过逆Park变换将同步旋转坐标系下的电压u_d、u_q转化为静止坐标系下的电压u_α、u_β，具体方法为：

u_α＝u_dcosθ_e-u_qsinθ_e

u_β＝u_dsinθ_e+u_qcosθ_e。

4.根据权利要求1所述的一种基于优化BP神经网络的永磁同步电机参数辨识方法，其特征在于：所述步骤(1)中已知参数永磁同步电机所需要已知的参数包括同步旋转坐标系下的d轴电感L_d、q轴电感L_q、永磁同步电机极对数p、永磁同步电机线电阻R_s、永磁同步电机磁链ψ_f、永磁同步电机的转动惯量J、阻尼系数B。

5.根据权利要求1所述的一种基于优化BP神经网络的永磁同步电机参数辨识方法，其特征在于：所述步骤(2)中转速传感器检测电机的机械角速度ω_m，通过计算得到的电机电角速度ω_e的具体方法为：

ω_e＝p*ω_m

通过位置传感器检测电机机械角度θ_m，并通过计算得到的电机电角度θ_e，其具体方法为：

θ_e＝p*θ_m

其中，p表示永磁同步电机PMSM极对数。

6.根据权利要求1所述的一种基于优化BP神经网络的永磁同步电机参数辨识方法，其特征在于：所述步骤(3)中，Clarke变换、Park变换描述如下：

将电机的三相电流i_a、i_b、i_c通过Clarke变换转化为静止坐标系下的电流i_α、i_β的具体方法为：

将静止坐标系下的电流i_α、i_β通过Park变换转化为同步旋转坐标系下的电流i_d、i_q的具体方法为：

i_d＝i_αcosθ_e+i_βsinθ_e

i_q＝-i_αsinθ_e+i_βcosθ_e。

7.根据权利要求1所述的一种基于优化BP神经网络的永磁同步电机参数辨识方法，其特征在于：所述步骤(5)中，改进鲸鱼优化算法优化BP神经网络训练TCWOA-BPNN，使用改进的鲸鱼优化算法获取BP神经网络的初始权值和阈值具体步骤如下：

A、鲸鱼优化算法的改进包括Tent混沌映射以及精英学习机制，Tent混沌映射表示为：

式中x₀为鲸鱼优化算法种群中的初始个体，x_n和x_n+1为种群中生成的第n个个体和第n+1个个体，rand[0，1]表示从0到1之间取随机数，精英学习机制表示为：

X_C＝j·k·[X_i-X^*]

式中X_C为鲸鱼优化算法生成的精英个体，j和k为随机参数值，X_i为当前鲸鱼优化算法种群中的个体，X^*为种群中当前位置最优个；

B、由于改进鲸鱼优化算法用于找到BP神经网络训练结果偏差最小的值，则自变量数量为1，则设置改进鲸鱼优化算法的自变量维度为1，设种群规模为N，最大迭代次数设置为T次；

C、设置概率p₁＝0.5为种群会按当前最优值继续寻优，设置概率p₂＝0.5为种群会从头开始寻优；

D、设置改进鲸鱼优化算法再得到的BP神经网络训练误差小于等于2％或迭代达到最大迭代次数时停止，改进鲸鱼优化算法停止工作后，得到的BP神经网络训练结果偏差最小的值对应的初始权值j₀和阈值k₀，将初始权值j₀和阈值k₀赋予BP神经网络作为初始权值和阈值进行训练。

8.根据权利要求1所述的一种基于优化BP神经网络的永磁同步电机参数辨识方法，其特征在于：所述步骤(5)中，改进鲸鱼优化算法优化BP神经网络训练，BP神经网络的参数设置如下：

A、BP神经网络的输入数据集为{i_d、i_q、u_d、u_q、ω_e}，则设置BP神经网络的输入数为5；输出数据集为{R_s、L_d、L_q、ψ_f}，则设置输出数为4；

B、BP神经网络的隐含层数量选择m层，隐含层神经元数量根据经验及实践确定为n个。

9.根据权利要求1所述的一种基于优化BP神经网络的永磁同步电机参数辨识方法，其特征在于：所述步骤(6)中，输出未知参数的永磁同步电机各项参数，具体可输出的参数包括：永磁同步电机线电阻R_s、d轴电感L_d、q轴电感L_q以及磁链ψ_f。

Images