永磁同步电机MTPA控制策略及matlab仿真

本文主要是介绍永磁同步电机MTPA控制策略及matlab仿真，希望对大家解决编程问题提供一定的参考价值，需要的开发者们随着小编来一起学习吧！

一、目的

二、主体框架

三、转PI速PI

四、电流PI

五、MPTA环节

六、仿真结果

七、附加问题

一、目的

1、给自己写一份学习笔记

2、自己对MTPA控制中感到困惑的点的说明

二、主体框架

内置式永磁同步电机，那么就有Ld != Lq，

对于表贴式Ld=Lq,id等于0时，就是MPTA控制了（从转矩方程分析）。

所以表贴式一般用id等于0控制，内置式更适合用mpta控制。

我是在双PI控制的基础上加入MPTA控制的，如下图一

图一

图二加入MPTA后

三、转PI速PI

图三转速PI

PI参数自己试凑了一些

图四 PI参数

四、电流PI

图五

电流PI参数

Id和Iq的PI环节是一样的

后面环节SVPWM，坐标变换等不具体展开

PMSM电机参数

五、MPTA环节

原理：参考文献《电动汽车用 IPMSM 矢量控制策略研究》郑飞

以上是文献的原文内容，但是我感觉其中有几个错误之处。第一个是图中红色底纹的部分，id不应该全部提出来；第二个是id*的公式，在其他文献上找到的是下面这个：

另附网上找的一篇公式的详细推导：https://www.cnblogs.com/iable/archive/2012/03/09/4206899.html

有了以上的公式，就可以搭建MPTA了，可以用模块，这里为了更快采用了M函数。

外部：

内部：用M函数写的代码如下

function [iq,id] = fcn(Te,Pn,flux,Ld,Lq)% id、iq是电流环的给定，也就是MTPA的输出，无法用星号表示% 第一步 计算iq
num1 = 8*Te*flux/(3*Pn);
num2 = flux^2-4*(Ld-Lq)^2;
num4 = num1^2-4*num2*((0.5*num1/flux)^2-flux^2);
iq = (num1+sqrt(num4))/(2*num2);%第二步 计算id
num3 = flux*0.5/(Lq-Ld);
num5 = 0.25*flux^2/(Ld-Lq)^2+iq^2;
id = num3-sqrt(num5);% num = flux^2+4*(Ld-Lq)^2*iq^2;
% id = (-flux+sqrt(num))/(2*(Ld-Lq));
end