电机仿真系列-基于最小二乘法的永磁同步电机参数辨识
基于最小二乘法的永磁同步电机参数辨识
永磁同步电机具有体积小、转动惯量低、结构简单等优点,被广泛应用于控制系统中。然而在实际应用过程中,控制系统会受到高温、负载等外界因素的影响,永磁同步电机的电感、转子磁链等参数会发生变化,导致系统振荡,影响实际控制效果。因此精确的参数辨识是为了达到更好地控制效果的必要条件。目前常用的电机参数辨识方法有:卡尔曼滤波算法、最小二乘法、遗传算法、粒子群算法等。考虑到最小二乘法结构简单、易于实现的优点,本期采用基于最小二乘法的永磁同步电机参数辨识方法。基于对永磁同步电机的推导,采用定子电压等可测变量作为输入输出,辨识出电机的定子电阻及d、q轴电枢电感。
1、最小二乘法
最小二乘参数辨识方法可以解决线性时变系统、线性定常系统、含有噪声的线性系统等问题。常见的最小二乘法有:递推最小二乘法、遗忘因子递推最小二乘法、修正补偿最小二乘法等,这些方法均可以用于系统参数辨识中。递推最小二乘法的计算公式可以参考往期的文章。
2、辨识模型的建立
永磁同步电机是一个非线性的多变量系统。在d、q旋转坐标下的电压方程为:
式中,iq和id分别为q、d的电流;uq和ud分别为q、d的电压;Rs、Lq、Ld分别为定子绕组的电阻和q、d轴的电感;Ψq、Ψd分别为q、d磁链的分量;ω为转子的电角速度。
将上述公式修改为最小二乘法的表达形式:
3、仿真实验
将S函数编写的参数辨识模块加入到永磁同步电动机控制系统模型中进行仿真实验。
function [sys,x0,str,ts] = Synchronous_demarcate2(t,x,u,flag)
%此程序为辨识电机参数的转动惯量
switch flagcase 0 %初始化[sys,x0,str,ts]=mdlInitializeSizes;case 2 %离散状态计算,下一步仿真时刻,终止仿真设定sys=[];%mdlUpdates(t,x,u);case 3 %输出信号计算sys=mdlOutputs(t,x,u);case {1,4,9} %输出信号计算sys=[];otherwiseDAStudio.error('Simulink:blocks:unhandledFlag', num2str(flag));
endfunction [sys,x0,str,ts]=mdlInitializeSizes %系统的初始化
sizes = simsizes;
sizes.NumContStates = 0; %设置系统连续状态的变量
sizes.NumDiscStates = 0; %设置系统离散状态的变量
sizes.NumOutputs = 1; %设置系统输出的变量
sizes.NumInputs = 3; %设置系统输入的变量
sizes.DirFeedthrough = 1; %如果在输出方程中显含输入变量u,则应该将本参数设置为1,输入不直接传到输出口
sizes.NumSampleTimes = 1; % 模块采样周期的个数% 需要的样本时间,一般为1.% 猜测为如果为n,则下一时刻的状态需要知道前n个状态的系统状态
sys = simsizes(sizes);
x0 = []; % 系统初始状态变量
str = []; % 保留变量,保持为空
ts = [-1 0]; % 采样时间[t1 t2] t1为采样周期,如果取t1=-1则将继承输入信号的采样周期;参数t2为偏移量,一般取为0global P_past2 theta_past2
P_past2 = 1e4 * eye(2,2); %一般取1e4 - 1e10
theta_past2 = [0.0001; 0.0001]; %一般取一个极小的正实向量function sys=mdlOutputs(t,x,u) %产生(传递)系统输出
%初值的确定
lambda = 0.99; %遗忘因子0-1
global P_past2 theta_past2
xt = [u(1) u(2)]; %1*2 fait
y = u(3);
I = [1 0;0 1];K = P_past2*xt'/(lambda + xt * P_past2*xt'); %2*1
P_new = 1/lambda*(I - K*xt) * P_past2; %2*2
theta_new = theta_past2 + K*(y-xt*theta_past2); %2*1P_past2 = P_new ;
theta_past2 = theta_new;sys(1) = theta_new(1);
function [sys,x0,str,ts] = Synchronous_demarcate(t,x,u,flag)
%此程序为辨识电机参数的定子电阻Rs, 永磁体磁通faif, dq轴电感L
switch flagcase 0 %初始化[sys,x0,str,ts]=mdlInitializeSizes;case 2 %离散状态计算,下一步仿真时刻,终止仿真设定sys=[];%mdlUpdates(t,x,u);case 3 %输出信号计算sys=mdlOutputs(t,x,u);case {1,4,9} %输出信号计算sys=[];otherwiseDAStudio.error('Simulink:blocks:unhandledFlag', num2str(flag));
endfunction [sys,x0,str,ts]=mdlInitializeSizes %系统的初始化
sizes = simsizes;
sizes.NumContStates = 0; %设置系统连续状态的变量
sizes.NumDiscStates = 0; %设置系统离散状态的变量
sizes.NumOutputs = 3; %设置系统输出的变量
sizes.NumInputs = 5; %设置系统输入的变量
sizes.DirFeedthrough = 1; %如果在输出方程中显含输入变量u,则应该将本参数设置为1,输入不直接传到输出口
sizes.NumSampleTimes = 1; % 模块采样周期的个数% 需要的样本时间,一般为1.% 猜测为如果为n,则下一时刻的状态需要知道前n个状态的系统状态
sys = simsizes(sizes);
x0 = []; % 系统初始状态变量
str = []; % 保留变量,保持为空
ts = [0 0]; % 采样时间[t1 t2] t1为采样周期,如果取t1=-1则将继承输入信号的采样周期;参数t2为偏移量,一般取为0global P_past theta_past
P_past = 1e4 * eye(3,3);
theta_past = [0;0;10000];function sys=mdlOutputs(t,x,u) %产生(传递)系统输出
%初值的确定
lambda = 1; %遗忘因子, 选择不遗忘 0.8~1
global P_past theta_past
id = u(1);
iq = u(2);
w = u(3); %电角速度 = 转子速度 * 极对数
diq = u(4);
uq = u(5);
xt = [-iq -w uq];
y = diq+w*id;
P_new = 1/lambda*(P_past-(P_past*(xt'*xt)*P_past)/(lambda+xt*P_past*xt'));
L = P_past*xt'/(lambda+xt*P_past*xt');
theta_new = theta_past + L*(y-xt*theta_past);
a = theta_new(1);
b = theta_new(2);
c = theta_new(3);
R = a/c;
Ke = b/c;
Ld = 1/c;
P_past = P_new ;
theta_past = theta_new;
sys(1) = R;
sys(2) = Ke;
sys(3) = Ld;
转动惯量参数辨识结果:
电磁电感参数辨识结果:
电角速度观测值:
负载观测结果:
三相电流观测结果:
仿真结果表明,采用递推最小二乘法的参数辨识算法能够准确辨识出电机的电磁参数,能够很好地跟随电机参数的变化,具有快速的收敛性和良好的辨识结果。如需了解更多仿真专题,可以关注我的微信公众号:xinnengyuanqiche666。
电机仿真系列-基于最小二乘法的永磁同步电机参数辨识相关推荐
- matlab最小二乘法参数辨识,基于最小二乘法的永磁同步电机在线参数辨识的仿真...
摘 要:较高性能的永磁同步电机矢量控制系统需要实时更新电机参数,文章中采用一种在线辨识永磁同步电机参数的方法.这种基于最小二乘法参数辨识方法是在转子同步旋转坐标系下进行的,通过MATLAB/SIMU ...
- 基于改进MRAS算法的永磁同步电机参数辨识
基于改进MRAS算法的永磁同步电机参数辨识 摘 要 永磁同步点电机参数辨识算法介绍 永磁同步电机数学模型 改进MRAS参数辨识算法 递推最小二乘法辨识原理 递推最小二乘法结合MRAS算法原理 仿真结果 ...
- 基于模型参考自适应的永磁同步电机参数辨识模型
基于模型参考自适应的永磁同步电机参数辨识模型 编号:7540655961097765三人班
- 永磁同步电机参数辨识模型,在线辨识,离线辨识
永磁同步电机参数辨识模型,在线辨识,离线辨识,电参数机械参数均可辨识,基于最小二乘法,滑模观测,电压注入,模型参考自适应等 机械参数在线 离线 ,电气参数在线 离线 指明要哪个 q ID:695066 ...
- 电机仿真系列-基于LabVIEW的电机测试系统研究
基于LabVIEW的电机测试系统研究 本设计基于LabVIEW设计一个直流有刷电机PID调试及控制助手,实现PID参数调节及电机转速及位置的控制.系统下位机采用STM32F103作为直流有刷电机控 ...
- 同步电机调速matlab,基于matlab的永磁同步电机调速系统的仿真word格式
基于matlab的永磁同步电机调速系统的仿真word格式 摘 要本文首先介绍了永磁同步电机的国内外发展状况,然后介绍了永磁同步电机的结构及原理,接着建立了永磁同步电机的数学模型,并在此基础上用MATL ...
- 基于ESO的永磁同步电机无感FOC 采用线性扩张状态观测器(LESO)估计电机反电势,利用锁相环从反电势中提取位置和转速信息
基于ESO的永磁同步电机无感FOC 1.采用线性扩张状态观测器(LESO)估计电机反电势,利用锁相环从反电势中提取位置和转速信息,从而实现无位置传感器控制: 2.提供算法对应的参考文献和仿真模型. 购 ...
- 同步电机调速matlab,基于matlab的永磁同步电机调速系统的仿真.doc
摘 要 本文首先介绍了永磁同步电机的国内外发展状况,然后介绍了永磁同步电机的结构及原理,接着建立了永磁同步电机的数学模型,并在此基础上用MATLAB进行了仿真,最后进行了仿真及仿真结果的分析. 永磁同 ...
- 基于 ESO-PLL 的永磁同步电机无位置传感器控制
基于 ESO-PLL 的永磁同步电机无位置传感器控制 1.PMSM 的无位置传感器控制方法分为两类,一类是适用于零.低速范围的高频注入法,另一类是适用于中.高速范围的观测器法.在中.高速范围,最常见的 ...
最新文章
- Tomcat的Classpath-常见问题以及如何解决
- 转载 - 背景透明,文字不透明【学习】
- PHP实现简易版区块链
- c语言宏定义比较三个数大小,C语言中两个宏进行大小对比,其中一个没有定义,这种行为如何定义。...
- leetcode —— 15. 三数之和
- mysql面试题sql语句_数据库MySQL经典面试题之SQL语句
- matlab中输出汉字的函数,matlab中的输出显示函数
- 软件测试用例优先级,软件测试用例的优先级划分方法
- 这就是艺术「GitHub 热点速览 v.22.25」
- 一周文章导读:fork() 成为负担;网络协议;(CPRI)带宽计算;互联网架构“高并发”;从MCU到FPGA;
- HTTP headers 详解 headers 常用属性
- (附源码)springboot教材订购系统的开发毕业设计081419
- 软件开发中的SD、SE、QA和RD是什么意思?
- 软考高级信息系统项目管理师如何备考?
- 弹性盒子 -- flex
- Maven的三种packaging方式(pom、jar、war)
- 第9章第19节:继续商业模式幻灯片的制作 [PowerPoint精美幻灯片实战教程]
- FM调制的FPGA实现
- Oracle常用命令(一)
- 我的springBoot框架学习总结