CN106533300A - 一种基于速度环模糊控制和高频注入法的无传感器控制系统 - Google Patents

Abstract

本发明公开了一种基于速度环模糊控制和高频注入法的无传感器控制系统，PMSM模块、Clark变换模块、Park变换模块、转子参数估算模块、高频信号注入模块、第一比较器模块、模糊控制器模块、MTPA模块、第二比较器模块、第一PI调节模块、第三比较器模块、第二PI调节模块、Park反变换模块、SVPWM模块和逆变器模块，在转子角度和转速的估算中，高频信号注入法尤其适合零速和低速，同时采用模糊控制器替代传统PI速度调节器。在永磁同步电机矢量控制系统中，模糊控制的鲁棒性强，干扰和参数变化对控制效果的影响被大大减弱，尤其适合于非线性、时变及纯滞后系统的控制，有利于模拟人工控制的过程和方法，增强控制系统的适应能力。

Description

一种基于速度环模糊控制和高频注入法的无传感器控制系统

技术领域

本发明涉及无速度传感器测速技术领域，具体涉及一种基于速度环模糊控制和高频注入法的无传感器控制系统。

背景技术

永磁同步电机(Permanent Magnet Synchronous Motor，简称PMSM)具有功率密度高、能量转换效率高、调速范围广、体积小、重量轻等优点，在工业、民用、军事等领域得到广泛的应用。

永磁同步电机的控制需要获得电机转子的位置和速度信息，目前应用比较普遍的位置传感器包括光电编码器、旋转变压器等装置，而这些装置的使用不但增加了系统的体积和成本，降低了系统的可靠性，也限制了永磁同步电机在特殊环境下的应用，为解决机械传感器带来的诸多缺陷，无传感器控制技术的研究已成为国内外的研究热点，并取得了一定成果，但还存在许多问题。最重要的是目前还没有一种单一的无传感器技术能够适用于在各种运行条件下有效地控制电机。现有技术中，或适用于低速运行，或适用于高速运行，或受电机参数影响较大，或计算量很大、结构复杂，或稳定性不是很好。

在电机速度检测过程中，机械传感器存在很多难以解决的缺点。如：在一些特殊的工作环境下(高温，高压)，其提供的信息精度不值得信赖；同时使用机械传感器使电机控制系统成本的增加、维护困难等。此外，因为常规PI控制器一般都会存在一个问题——积分饱和。所谓积分饱和，是指系统存在一个方向的偏差时，PI控制器的积分环节不断累加，最终到达控制器的限幅值，即使继续积分作用，控制器输出不变，所以出现了积分饱和。一旦系统出现反向偏差，控制器反向积分，控制器输出逐渐从饱和区退出，退出的时间与之间积分饱和的深度有关。但是，在退饱和的时间内，控制器输出还是在限幅值，此时容易出现调节滞后，导致系统性能变差。

发明内容

为了克服现有的基于无速度传感器的永磁同步电机的转子角度、转速估计方法存在的原理复杂、计算量大以及积分饱和的问题，现在特别提出一种具有较高动态性能且易于工程实现的一种基于速度环模糊控制和高频注入法的无传感器控制系统，通过模糊控制器调整PI调节器的比例积分系数，以使PI调节器能在电机很宽的速度范围内都具有良好的动稳态性能。

为了达到上述发明目的，解决其技术问题所采用的技术方案如下：

一种基于速度环模糊控制和高频注入法的无传感器控制系统，包括PMSM模块、Clark变换模块、Park变换模块、转子参数估算模块、高频信号注入模块、第一比较器模块、模糊控制器模块、MTPA模块、第二比较器模块、第一PI调节模块、第三比较器模块、第二PI调节模块、Park反变换模块、SVPWM模块和逆变器模块，其中：

所述PMSM模块，用于检测输出三相电流i_a、i_b和i_c；

所述Clark变换模块，用于将所述PMSM模块输出的三相电流i_a、i_b和i_c通过Clark变换后输出两相静止直角坐标系α-β下的两相定子电流i_α和i_β；

所述Park变换模块，用于将所述Clark变换模块输出的两相定子电流i_α和i_β通过Park变换后输出两相同步旋转坐标系d-q下的两相电流i_d和i_q；

所述转子参数估算模块，用于将所述Clark变换模块输出的两相定子电流i_α和i_β、所述高频信号注入模块注入的旋转两相高频电压信号u_asi和u_βsi与所述PMSM模块输出的转矩T_e一并输入转子参数估算模块内的全维观测器中进行估算处理，估算出转子转速的估计值和转子位置的估计值估算出转子转速的估计值乘以一常数得到估算的转子转速n；

所述第一比较器模块，用于将估算的转子转速n与实际的转子转速n*进行作差运算；

所述模糊控制器模块，用于将所述第一比较器模块比较的差值通过PI调节后输出参考转矩

所述MTPA模块，用于将所述模糊控制器模块输出的参考转矩通过最大转矩电流比控制后得到q轴参考电流和d轴参考电流

所述第二比较器模块，用于将所述MTPA模块输出的q轴参考电流与所述Park变换模块中输出的电流i_q进行作差运算；

所述第一PI调节模块，用于将所述第二比较器模块比较的差值通过PI调节后输出q轴参考电压u_q；

所述第三比较器模块，用于将所述MTPA模块输出的d轴参考电流与所述Park变换模块中输出的电流i_d进行作差运算；

所述第二PI调节模块，用于将所述第三比较器模块比较的差值通过PI调节后输出d轴参考电压u_d；

所述Park反变换模块，用于将所述第一PI调节模块输出的q轴参考电压u_q和所述第二PI调节模块输出的d轴参考电压u_d通过Park反变换后输出两相静止直角坐标系α-β下的两相控制电压u_α和u_β；

所述SVPWM模块，用于将所述Park反变换模块输出的两相控制电压u_α和u_β与所述高频信号注入模块注入的旋转两相高频电压信号u_asi和u_βsi进行叠加后进行空间矢量调制，输出PWM波形至所述逆变器模块，所述逆变器模块向所述PMSM模块输入三相电压u_a、u_b和u_c，从而控制所述PMSM模块。

进一步的，还包括A/D转换器模块和D/A转换器模块，其中：

所述A/D转换器模块，用于将所述第一比较器模块作差运算得到精确值e经过A/D转换后把模拟量转换成数字量并送入所述模糊控制器模块；

所述D/A转换器模块，用于将所述A/D转换器模块中得到的数字量经过所述模糊控制器模块模糊处理后输出的精确值u通过D/A转换后把数字量转换为模拟量，并输出参考转矩

进一步的，所述模糊控制器模块包括模糊量化处理子模块、推理机子模块、规则库子模块和去模糊化处理子模块，其中：

所述模糊量化处理子模块，用于将所述A/D转换器模块中得到的数字量经过模糊量化处理，得到一模糊值e；

所述推理机子模块，用于将上述模糊值e结合所述规则库子模块中的模糊控制规则R根据推理合成规则进行模糊决策，得到模糊控制量u，模糊值u＝e*R；

所述去模糊化处理子模块，用于将所述推理机子模块中得出的模糊值u进行去模糊化处理，得到精确值u。

进一步的，所述转子参数估算模块包括同步旋转高通滤波器子模块、外差计算子模块和全维观测器子模块，其中：

所述同步旋转高通滤波器子模块，用于将所述Clark变换模块输出的两相定子电流i_α和i_β通过同步旋转滤波后，剩下的电流分量只包含高频电流负序成分i_αi-in和i_βi-in；

所述外差计算器子模块，用于将所述同步旋转高通滤波器子模块滤波后得到的高频电流负序成分i_αi-in和i_βi-in与所述高频信号注入模块注入的旋转两相高频电压信号u_asi和u_βsi进行外差法运算，得出转子位置的误差角度θ_e；

所述全维观测器子模块，用于将所述外差计算器子模块得到的误差角度θ_e与所述PMSM模块输出的转矩T_e一并输入进行估算处理，得到估计角度和估计速度

进一步的，所述同步旋转高通滤波器子模块具体包括以下步骤：

首先，建立交流永磁同步电机在两相静止直角坐标系α-β中的数学模型：

u_βs＝R_Si_βs+Pψ_βs (1)

u_αs＝R_Si_αs+Pψ_αs (2)

式中，u_αs和u_βs为两相静止直角坐标系α-β中电压，R_s为定子电阻，i_αs和i_βs为两相静止直角坐标系α-β中电流，P为微分算子，ψ_αs和ψ_βs代表定子磁链；

其中，磁链方程为：

其中：

式中，为平均电感，为调制电感，θ_r为d轴领先A相相轴的空间电角度，L_md、L_mq为阻尼绕组归算到定子侧的d、q分量，i_Q、i_D分别为归算后的转子交、直轴阻尼绕组电流，ψ_f代表转子永磁磁链。

进一步的，所述同步旋转高通滤波器子模块中，通过同步旋转滤波后，剩下的电流分量只含高频电流负序成分，其矢量表达式为：

式中，θ_r为d轴领先A相相轴的空间电角度，θ_i＝ω_it，ω_i代表该注入电压信号的角频率，θ_i代表该注入电压信号的角度，i_in代表电流负序的幅值。

进一步的，所述外差计算器子模块中注入的电压信号：

式中，U_si代表在静止坐标系上注入高频旋转电压的幅值，ω_i代表注入电压信号u_αsi的角频率；

载波信号注入后，电机坐标下的电压方程为：

式中，U_se代表正序电流幅值，ω_r代表转子角频率；

在此高频电压注入下，产生的电流将由三部分组成：第一部分是与注入的电压旋转方向相同的正序电流，第二部分是与旋转电压方向相反的负序电流，第三部分是由三相绕组不对称产生的零序电流，电流响应可以表示为：

其中，

式中，θ_r为d轴领先A相相轴的空间电角度，θ_i代表该注入电压信号的角频率为ω_i，i_in代表电流负序的幅值，U_si代表在静止坐标系上注入高频旋转电压的幅值，ω_i代表注入电压信号的角频率，L代表平均电感，ΔL代表空间调制电感；

从公式(8)中得出，只有高频响应电流的负序成分中含有转子位置信息，通过滤波器将电源产生的频率成分和正序电流分量滤除，然后用外差法得出转子位置的误差角度θ_e，再利用全维观测器提取出转子位置信息。

进一步的，所述外差计算器子模块中外差法运算包括将公式(9)中i_αi、i_βi分别乘以和然后作差：

式中，θ_r为d轴领先A相相轴的空间电角度，代表高频电压注入法获得的转子初判角，ω_i代表注入电压信号的角频率；

其中，第一项为含电流的高频分量，第二项为仅含转子位置的信息，通过低通滤波可得转子位置的误差信号，从而：

在角度误差很小的情况下，

进一步的，所述全维观测器中转子转速的估计值通过以下公式求得：

交流永磁同步电机的运动方程可以表示为：

式中,J为转动惯量，T_L代表负载转矩；

电机转子在一个采样周期T_s上的角位移公式是:

式中，t₀代表转子开始时间，T代表转子经过时间；

采样周期很短，上式表示为:

式中，ω_r代表转子角速度；

由式(13)和(15)可以得到：

电机系统中负载变换缓慢，所以可认为：

将式(13)、(16)及(17)改写成矩阵形式:

式中，l₁、l₂和l₃三个表示的是在观测器中的增益值；

通过极点配置的方式来设置合理的全维观测器，离散化后的全维观测器方程为：

进一步的，所述高频信号注入模块向两相静止直角坐标系α-β中注入高频旋转电压信号u_asi和u_βsi为：

u_asi＝v_si sinω_it (20)

u_βsi＝v_si cosω_it (21)

其中，v_si是注入的高频电压信号的幅值,ω_i为注入的高频电压信号的角频率。

本发明由于采用以上技术方案，使之与现有技术相比，具有以下的优点和积极效果：

1、本发明一种基于速度环模糊控制和高频注入法的无传感器控制系统，对系统扰动、参数摄动等不确定性因素具有鲁棒性，因此可以更好的实现永磁同步电机的无传感器控制；

2、本发明所设计旋转高频注入法和模糊控制相结合下，能及时并准确的跟踪电动机的转速和转角变化，具有快速性好，控制准确性高，动态性能好，鲁棒性强的特点，而且所设计的观测器无论在硬件和软件上实施起来都较为方便，具有一定的实用性；

3、本发明通过采用全维观测器实现状态估计，显著提高了转子位置与速度的估计精确度；

4、本发明应用模糊控制器调整PI调节器的比例积分系数，使PI自适应调节器在电机很宽的速度范围内都有良好的动稳态性能，从而使观测器在低速时可以抑制检测的转子位置角度的小幅振荡，高速时减小其角度的相位延迟，提高了转子位置的检测精度；

5、本发明模糊控制的鲁棒性强，干扰和参数变化对控制效果的影响被大大减弱，尤其适合于非线性、时变及纯滞后系统的控制，模糊控制是基于启发性的知识及语言决策规则设计的，这有利于模拟人工控制的过程和方法，增强控制系统的适应能力，具有一定的智能水平，对那些数学模型难以获取，动态特性不易掌握或变化非常显著的对象非常适用；

6、本发明具有低成本、控制算法简单、转速及位置的估算速度及精度高等优点。

附图说明

为了更清楚地说明本发明实施例的技术方案，下面将对实施例描述中所需要使用的附图作简单的介绍。显而易见，下面描述中的附图仅仅是本发明的一些实施例，对于本领域技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。附图中：

图1是本发明一种基于速度环模糊控制和高频注入法的无传感器控制系统的整体系统架构图；

图2是本发明一种基于速度环模糊控制和高频注入法的无传感器控制系统中的模糊控制器结构图；

图3是本发明一种基于速度环模糊控制和高频注入法的无传感器控制系统中e的隶属函数图；

图4是本发明一种基于速度环模糊控制和高频注入法的无传感器控制系统中de的隶属函数图；

图5是本发明一种基于速度环模糊控制和高频注入法的无传感器控制系统中du的隶属函数图；

图6是本发明一种基于速度环模糊控制和高频注入法的无传感器控制系统的实际角度和估计角度仿真图；

图7是本发明一种基于速度环模糊控制和高频注入法的无传感器控制系统的转角误差图；

【主要符号标记】

1-PMSM模块；

2-Clark变换模块；

3-Park变换模块；

4-转子参数估算模块；

5-高频信号注入模块；

6-第一比较器模块；

7-模糊控制器模块；

8-MTPA模块；

9-第二比较器模块；

10-第一PI调节模块；

11-第三比较器模块；

12-第二PI调节模块；

13-Park反变换模块；

14-SVPWM模块；

15-逆变器模块；

16-A/D转换器模块；

17-D/A转换器模块；

71-模糊量化处理子模块；

72-推理机子模块；

73-规则库子模块；

74-去模糊化处理子模块。

具体实施方式

以下将结合本发明的附图，对本发明实施例中的技术方案进行清楚、完整的描述和讨论，显然，这里所描述的仅仅是本发明的一部分实例，并不是全部的实例，基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动的前提下所获得的所有其他实施例，都属于本发明的保护范围。

如图1所示，本发明公开了一种基于速度环模糊控制和高频注入法的无传感器控制系统，包括PMSM(Permanent Magnet Synchronous Motor，永磁同步电机)模块1、Clark变换模块2、Park变换模块3、转子参数估算模块4、高频信号注入模块5、第一比较器模块6、模糊控制器模块7、MTPA(Maximum Torque Per Ampere，最大转矩电流比)模块8、第二比较器模块9、第一PI调节模块10、第三比较器模块11、第二PI调节模块12、Park反变换模块13、SVPWM(Space Vector Pulse Width Modulation，空间矢量脉宽调制)模块14和逆变器模块15，其中：

所述PMSM模块1，用于检测输出三相电流i_a、i_b和i_c；

所述Clark变换模块2，用于将所述PMSM模块1输出的三相电流i_a、i_b和i_c通过Clark变换后输出两相静止直角坐标系α-β下的两相定子电流i_α和i_β；

所述Park变换模块3，用于将所述Clark变换模块2输出的两相定子电流i_α和i_β通过Park变换后输出两相同步旋转坐标系d-q下的两相电流i_d和i_q；

所述转子参数估算模块4，用于将所述Clark变换模块2输出的两相定子电流i_α和i_β、所述高频信号注入模块5注入的旋转两相高频电压信号u_asi和u_βsi与所述PMSM模块1输出的转矩T_e一并输入转子参数估算模块4内的全维观测器中进行估算处理，估算出转子转速的估计值和转子位置的估计值估算出转子转速的估计值乘以一常数得到估算的转子转速n；

所述第一比较器模块6，用于将估算的转子转速n与实际的转子转速n*进行作差运算；

所述模糊控制器模块7，用于将所述第一比较器模块6比较的差值通过PI调节后输出参考转矩

所述MTPA模块8，用于将所述模糊控制器模块7输出的参考转矩通过最大转矩电流比控制后得到q轴参考电流和d轴参考电流

所述第二比较器模块9，用于将所述MTPA模块8输出的q轴参考电流与所述Park变换模块3中输出的电流i_q进行作差运算；

所述第一PI调节模块10，用于将所述第二比较器模块9比较的差值通过PI调节后输出q轴参考电压u_q；

所述第三比较器模块11，用于将所述MTPA模块8输出的d轴参考电流与所述Park变换模块3中输出的电流i_d进行作差运算；

所述第二PI调节模块12，用于将所述第三比较器模块11比较的差值通过PI调节后输出d轴参考电压u_d；

所述Park反变换模块13，用于将所述第一PI调节模块10输出的q轴参考电压u_q和所述第二PI调节模块12输出的d轴参考电压u_d通过Park反变换后输出两相静止直角坐标系α-β下的两相控制电压u_α和u_β；

所述SVPWM模块14，用于将所述Park反变换模块13输出的两相控制电压u_α和u_β与所述高频信号注入模块5注入的旋转两相高频电压信号u_asi和u_βsi进行叠加后进行空间矢量调制，输出PWM波形至所述逆变器模块14，所述逆变器模块14向所述PMSM模块1输入三相电压u_a、u_b和u_c，从而控制所述PMSM模块1。

具体的，在所述Clark变换模块2中将所述PMSM模块1输出的三相电流i_a、i_b和i_c通过Clark变换后输出两相静止直角坐标系α-β下的两相定子电流i_α和i_β，具体涉及的换算公式如下：

具体的，在所述Park变换模块3中将所述Clark变换模块2输出的两相定子电流i_α和i_β通过Park变换后输出两相同步旋转坐标系d-q下的两相电流i_d和i_q，具体涉及的换算公式如下：

其中，为估算的转子角。

具体的，在所述转子参数估算模块4中，估算出转子转速的估计值与估算的转子转速n之间的关系为：

即，所述常数为9.55。

图2是本发明中模糊控制系统框图，给定值为实际的给定速度，与全维观测器反馈的速度作差，得到了速度的差值即精确值e，精确值e经过A/D转换器把模拟量转换成数字量，送入模糊控制器，经过模糊控制器处理后输出精确值u，精确值u经过D/A转化器吧数字量转换为模拟量。

其中，模糊控制器的控制规律由计算机的程序实现，实现一步模糊控制算法的过程是：微机采样获取被控制对象的精确值，然后将此量与给定值比较得到误差信号e；一般选误差信号e作为模糊控制器的一个输入量，把e的精确量进行模糊量化变成模糊量，误差e的模糊量可用相应的模糊语言表示；从而得到误差e的模糊语言集合的一个子集e(实际上是一个模糊向量)；再由模糊向量e和模糊控制规则R(模糊关系)根据推理的合成规则进行模糊决策，得到模糊控制量u为u＝e·R。

式中u为一个模糊量；为了对被控对象(PMSM)施加精确的控制，还需要将模糊量u进行去模糊化处理转换为精确量：得到精确数字量后，经数模转换变为精确的模拟量送给执行机构(包括PI调节器、Park反变换和空间矢量调制SVPWM)，对被控对象进行一步控制；然后，进行第二次采样，完成第二步控制，这样循环下去，就实现了被控对象的模糊控制。

本实施例中，结合图2进一步还包括A/D转换器模块16和D/A转换器模块17，其中：

所述A/D转换器模块16，用于将所述第一比较器模块6作差运算得到精确值e经过A/D转换后把模拟量转换成数字量并送入所述模糊控制器模块7；

所述D/A转换器模块17，用于将所述A/D转换器模块中得到的数字量经过所述模糊控制器模块7模糊处理后输出的精确值u通过D/A转换后把数字量转换为模拟量，并输出参考转矩

一实施例中，所述模糊控制器模块7包括模糊量化处理子模块71、推理机子模块72、规则库子模块73和去模糊化处理子模块74，其中：

所述模糊量化处理子模块71，用于将所述A/D转换器模块16中得到的数字量经过模糊量化处理，得到一模糊值e；

所述推理机子模块72，用于将上述模糊值e结合所述规则库子模块73中的模糊控制规则R根据推理合成规则进行模糊决策，得到模糊控制量u，模糊值u＝e*R；

所述去模糊化处理子模块74，用于将所述推理机子模块72中得出的模糊值u进行去模糊化处理，得到精确值u。

具体的，在所述Park反变换模块13中，将所述第一PI调节模块10输出的q轴参考电压u_q和所述第二PI调节模块12输出的d轴参考电压u_d通过Park反变换后输出两相静止直角坐标系α-β下的两相控制电压u_α和u_β，具体涉及以下换算公式：

其中，为估算的转子角。

本实施例中，所述转子参数估算模块4包括同步旋转高通滤波器子模块、外差计算子模块和全维观测器子模块，其中：

所述同步旋转高通滤波器子模块，用于将所述Clark变换模块2输出的两相定子电流i_α和i_β通过同步旋转滤波后，剩下的电流分量只包含高频电流负序成分i_αi-in和i_βi-in；

所述外差计算器子模块，用于将所述同步旋转高通滤波器子模块滤波后得到的高频电流负序成分i_αi-in和i_βi-in与所述高频信号注入模块5注入的旋转两相高频电压信号u_asi和u_βsi进行外差法运算，得出转子位置的误差角度θ_e；

所述全维观测器子模块，用于将所述外差计算器子模块得到的误差角度θ_e与所述PMSM模块1输出的转矩T_e一并输入进行估算处理，得到估计角度和估计速度

进一步的，所述同步旋转高通滤波器子模块具体包括以下步骤：

首先，建立交流永磁同步电机在两相静止直角坐标系α-β中的数学模型：

u_βs＝R_Si_βs+Pψ_βs (1)

u_αs＝R_Si_αs+Pψ_αs (2)

式中，u_αs和u_βs为两相静止直角坐标系α-β中电压，R_s为定子电阻，i_αs和i_βs为两相静止直角坐标系α-β中电流，P为微分算子，ψ_αs和ψ_βs代表定子磁链；

其中，磁链方程为：

其中：

式中，为平均电感，为调制电感，θ_r为d轴领先A相相轴的空间电角度，L_md、L_mq为阻尼绕组归算到定子侧的d、q分量，i_Q、i_D分别为归算后的转子交、直轴阻尼绕组电流，ψ_f代表转子永磁磁链。

进一步的，所述同步旋转高通滤波器子模块中，通过同步旋转滤波后，剩下的电流分量只含高频电流负序成分，其矢量表达式为：

式中，θ_r为d轴领先A相相轴的空间电角度，θ_i＝ω_it，ω_i代表该注入电压信号的角频率，θ_i代表该注入电压信号的角度，i_in代表电流负序的幅值。

进一步的，所述外差计算器子模块中注入的电压信号：

式中，U_si代表在静止坐标系上注入高频旋转电压的幅值，ω_i代表注入电压信号u_αsi的角频率；

载波信号注入后，电机坐标下的电压方程为：

式中，U_se代表正序电流幅值，ω_r代表转子角频率；

在此高频电压注入下，产生的电流将由三部分组成：第一部分是与注入的电压旋转方向相同的正序电流，第二部分是与旋转电压方向相反的负序电流，第三部分是由三相绕组不对称产生的零序电流，电流响应可以表示为：

其中，

式中，θ_r为d轴领先A相相轴的空间电角度，θ_i代表该注入电压信号的角频率为ω_i，i_in代表电流负序的幅值，U_si代表在静止坐标系上注入高频旋转电压的幅值，ω_i代表注入电压信号的角频率，L代表平均电感，ΔL代表空间调制电感；

从公式(8)中得出，只有高频响应电流的负序成分中含有转子位置信息，通过滤波器将电源产生的频率成分和正序电流分量滤除，然后用外差法得出转子位置的误差角度θ_e，再利用全维观测器提取出转子位置信息。

进一步的，所述外差计算器子模块中外差法运算包括将公式(9)中i_αi、i_βi分别乘以和然后作差：

式中，θ_r为d轴领先A相相轴的空间电角度，代表高频电压注入法获得的转子初判角，ω_i代表注入电压信号的角频率；

其中，第一项为含电流的高频分量，第二项为仅含转子位置的信息，通过低通滤波可得转子位置的误差信号，从而：

在角度误差很小的情况下，

进一步的，所述全维观测器中转子转速的估计值通过以下公式求得：

交流永磁同步电机的运动方程可以表示为：

式中,J为转动惯量，T_L代表负载转矩；

电机转子在一个采样周期T_s上的角位移公式是:

式中，t₀代表转子开始时间，T代表转子经过时间；

采样周期很短，上式表示为:

式中，ω_r代表转子角速度；

由式(13)和(15)可以得到：

电机系统中负载变换缓慢，所以可认为：

将式(13)、(16)及(17)改写成矩阵形式:

式中，l₁、l₂和l₃三个表示的是在观测器中的增益值；

根据控制原理知识可知：系统稳定的条件是系统闭环传递函数的所有零极点必须在s平面的左半平面。但又考虑到系统的动态性能的要求，通常取零极点远离虚轴。所以要综合以上因素，可以通过极点配置的方式来设置合理的全维观测器，离散化后的全维观测器方程为：

进一步的，所述高频信号注入模块5向两相静止直角坐标系α-β中注入高频旋转电压信号u_asi和u_βsi为：

u_asi＝v_si sinω_it (20)

u_βsi＝v_si cosω_it (21)

其中，v_si是注入的高频电压信号的幅值,ω_i为注入的高频电压信号的角频率。

附图3、图4和图5的所有模糊集合的论域选为[-1,1]。权衡控制精度和计算复杂度，模糊集合子元素选为7个，分别为NL、NM、NS、ZO、PS、PM、PL。量化因子K_e、K_i的选择，实际中应考虑性能需求以及e和de的变化情况，选取合理的调节范围。假设e和de的论域范围分别为[-m,m]和[-n,n]，其中满足隶属函数的选择三角形和梯形隶属函数，因为相对而言选择三角形和梯形隶属函数控制器有较好的性能。推理和解模糊方法选择MAMDANI模糊推理和重心解模糊法。

模糊规则库通常是基于专家经验或过程知识生成的控制规则集合。对于永磁同步电机调速系统，设计的模糊控制器是针对速度控制，所以控制规则也是基于速度响应过程。

如果e>0、de<0，此时速度趋向给定值，应该给以较小的控制器输出；

如果e<0、de<0，此时出现速度超调，应该尽快通过控制器抑制超调；

如果e<0、de>0，此时抑制起到作用，速度回归给定值，控制器输出应该较小；

如果e>0、de>0，此时速度跟踪不上给定，控制器应该给以较大输出。

图6是本发明一种基于旋转高频注入法和模糊PI控制的无速度传感器控制方法的实际角度和估计角度仿真图，虚线代表的是实际角度，实线代表的是估计的角度。从图中可以看出本发明的转子位置跟踪效果很好，快速性好，角度得波形在1s时有所波动，原因是在1s时负载转矩从3N.m增加到5N.m，并且很快就稳定。从整体看，实际角度和估计角度的波动都比较小。

图7是本发明一种基于旋转高频注入法和模糊PI控制的无速度传感器控制方法的转角误差图，该图表示实际转角与估计转角之差，从图中可以看出转角误差几乎稳定在-0.1到0.1之间，表明了转角跟踪效果好。

以上所述，仅为本发明较佳的具体实施方式，但本发明的保护范围并不局限于此，任何熟悉本技术领域的技术人员在本发明揭露的技术范围内，可轻易想到的变化或替换，都应涵盖在本发明的保护范围之内。因此，本发明的保护范围应该以权利要求的保护范围为准。

Claims (10)

1.一种基于速度环模糊控制和高频注入法的无传感器控制系统，其特征在于，包括PMSM模块、Clark变换模块、Park变换模块、转子参数估算模块、高频信号注入模块、第一比较器模块、模糊控制器模块、MTPA模块、第二比较器模块、第一PI调节模块、第三比较器模块、第二PI调节模块、Park反变换模块、SVPWM模块和逆变器模块，其中：

所述PMSM模块，用于检测输出三相电流i_a、i_b和i_c；

所述Clark变换模块，用于将所述PMSM模块输出的三相电流i_a、i_b和i_c通过Clark变换后输出两相静止直角坐标系α-β下的两相定子电流i_α和i_β；

所述Park变换模块，用于将所述Clark变换模块输出的两相定子电流i_α和i_β通过Park变换后输出两相同步旋转坐标系d-q下的两相电流i_d和i_q；

所述转子参数估算模块，用于将所述Clark变换模块输出的两相定子电流i_α和i_β、所述高频信号注入模块注入的旋转两相高频电压信号u_asi和u_βsi与所述PMSM模块输出的转矩T_e一并输入转子参数估算模块内的全维观测器中进行估算处理，估算出转子转速的估计值和转子位置的估计值估算出转子转速的估计值乘以一常数得到估算的转子转速n；

所述第一比较器模块，用于将估算的转子转速n与实际的转子转速n*进行作差运算；

所述模糊控制器模块，用于将所述第一比较器模块比较的差值通过PI调节后输出参考转矩

所述MTPA模块，用于将所述模糊控制器模块输出的参考转矩通过最大转矩电流比控制后得到q轴参考电流和d轴参考电流

所述第二比较器模块，用于将所述MTPA模块输出的q轴参考电流与所述Park变换模块中输出的电流i_q进行作差运算；

所述第一PI调节模块，用于将所述第二比较器模块比较的差值通过PI调节后输出q轴参考电压u_q；

所述第三比较器模块，用于将所述MTPA模块输出的d轴参考电流与所述Park变换模块中输出的电流i_d进行作差运算；

所述第二PI调节模块，用于将所述第三比较器模块比较的差值通过PI调节后输出d轴参考电压u_d；

所述Park反变换模块，用于将所述第一PI调节模块输出的q轴参考电压u_q和所述第二PI调节模块输出的d轴参考电压u_d通过Park反变换后输出两相静止直角坐标系α-β下的两相控制电压u_α和u_β；

所述SVPWM模块，用于将所述Park反变换模块输出的两相控制电压u_α和u_β与所述高频信号注入模块注入的旋转两相高频电压信号u_asi和u_βsi进行叠加后进行空间矢量调制，输出PWM波形至所述逆变器模块，所述逆变器模块向所述PMSM模块输入三相电压u_a、u_b和u_c，从而控制所述PMSM模块。

2.根据权利要求1所述的一种基于速度环模糊控制和高频注入法的无传感器控制系统，其特征在于，还包括A/D转换器模块和D/A转换器模块，其中：

所述A/D转换器模块，用于将所述第一比较器模块作差运算得到精确值e经过A/D转换后把模拟量转换成数字量并送入所述模糊控制器模块；

所述D/A转换器模块，用于将所述A/D转换器模块中得到的数字量经过所述模糊控制器模块模糊处理后输出的精确值u通过D/A转换后把数字量转换为模拟量，并输出参考转矩

3.根据权利要求2所述的一种基于速度环模糊控制和高频注入法的无传感器控制系统，其特征在于，所述模糊控制器模块包括模糊量化处理子模块、推理机子模块、规则库子模块和去模糊化处理子模块，其中：

所述模糊量化处理子模块，用于将所述A/D转换器模块中得到的数字量经过模糊量化处理，得到一模糊值e；

所述推理机子模块，用于将上述模糊值e结合所述规则库子模块中的模糊控制规则R根据推理合成规则进行模糊决策，得到模糊控制量u，模糊值u＝e*R；

所述去模糊化处理子模块，用于将所述推理机子模块中得出的模糊值u进行去模糊化处理，得到精确值u。

4.根据权利要求1所述的一种基于速度环模糊控制和高频注入法的无传感器控制系统，其特征在于，所述转子参数估算模块包括同步旋转高通滤波器子模块、外差计算子模块和全维观测器子模块，其中：

所述同步旋转高通滤波器子模块，用于将所述Clark变换模块输出的两相定子电流i_α和i_β通过同步旋转滤波后，剩下的电流分量只包含高频电流负序成分i_αi-in和i_βi-in；

所述外差计算器子模块，用于将所述同步旋转高通滤波器子模块滤波后得到的高频电流负序成分i_αi-in和i_βi-in与所述高频信号注入模块注入的旋转两相高频电压信号u_asi和u_βsi进行外差法运算，得出转子位置的误差角度θ_e；

所述全维观测器子模块，用于将所述外差计算器子模块得到的误差角度θ_e与所述PMSM模块输出的转矩T_e一并输入进行估算处理，得到估计角度和估计速度

5.根据权利要求4所述的一种基于速度环模糊控制和高频注入法的无传感器控制系统，其特征在于，所述同步旋转高通滤波器子模块具体包括以下步骤：

首先，建立交流永磁同步电机在两相静止直角坐标系α-β中的数学模型：

u_βs＝R_Si_βs+Pψ_βs (1)

u_αs＝R_Si_αs+Pψ_αs (2)

式中，u_αs和u_βs为两相静止直角坐标系α-β中电压，R_s为定子电阻，i_αs和i_βs为两相静止直角坐标系α-β中电流，P为微分算子，ψ_αs和ψ_βs代表定子磁链；

其中，磁链方程为：

$\left[\begin{array}{c}{\mathrm{\&psi;}}_{\mathrm{\&beta;}s}\\ {\mathrm{\&psi;}}_{\mathrm{\&alpha;}s}\end{array}\right]=L_s\left(2{\mathrm{\&theta;}}_r\right)\left[\begin{array}{c}{i}_{\mathrm{\&beta;}s}\\ i_{\mathrm{\&alpha;}s}\end{array}\right]+L_{sr}\left({\mathrm{\&theta;}}_r\right)\left[\begin{array}{c}{i}_Q\\ i_D\end{array}\right]+{\mathrm{\&psi;}}_f\left[\begin{array}{c}sin{\mathrm{\&theta;}}_r\\ {\mathrm{cos\&theta;}}_r\end{array}\right]---\left(3\right)$

其中：

$L_s\left(2{\mathrm{\&theta;}}_r\right)=\left[\begin{array}{cc}L+{\mathrm{\&Delta;Lcos\&theta;}}_r& -{\mathrm{Lsin\&theta;}}_r\\ -{\mathrm{Lsin\&theta;}}_r& L-{\mathrm{\&Delta;Lcos\&theta;}}_r\end{array}\right]---\left(4\right)$

$L_{sr}\left({\mathrm{\&theta;}}_r\right)=\left[\begin{array}{cc}{L}_{mq}{\mathrm{cos\&theta;}}_r& L_{md}{\mathrm{sin\&theta;}}_r\\ -L_{mq}{\mathrm{sin\&theta;}}_r& L_{md}{\mathrm{cos\&theta;}}_r\end{array}\right]---\left(5\right)$

式中，为平均电感，为调制电感，θ_r为d轴领先A相相轴的空间电角度，L_md、L_mq为阻尼绕组归算到定子侧的d、q分量，i_Q、i_D分别为归算后的转子交、直轴阻尼绕组电流，ψ_f代表转子永磁磁链。

6.根据权利要求4所述的一种基于速度环模糊控制和高频注入法的无传感器控制系统，其特征在于，所述同步旋转高通滤波器子模块中，通过同步旋转滤波后，剩下的电流分量只含高频电流负序成分，其矢量表达式为：

$i_{\mathrm{\&alpha;}\mathrm{\&beta;}i-in}=i_{in}{e}^{j(2{\mathrm{\&theta;}}_r-{\mathrm{\&theta;}}_i+\mathrm{\&pi;}/2)}---\left(6\right)$

式中，θ_r为d轴领先A相相轴的空间电角度，θ_i＝ω_it，ω_i代表该注入电压信号的角频率，θ_i代表该注入电压信号的角度，i_in代表电流负序的幅值。

7.根据权利要求4所述的一种基于速度环模糊控制和高频注入法的无传感器控制系统，其特征在于，所述外差计算器子模块中注入的电压信号：

$\left[\begin{array}{c}{u}_{\mathrm{\&alpha;}si}\\ u_{\mathrm{\&beta;}si}\end{array}\right]=U_{si}\left[\begin{array}{c}-sin\left({\mathrm{\&omega;}}_{i}t\right)\\ \cos \left({\mathrm{\&omega;}}_{i}t\right)\end{array}\right]=U_{si}{e}^{{\mathrm{j\&omega;}}_{i}t}---\left(7\right)$

式中，U_si代表在静止坐标系上注入高频旋转电压的幅值，ω_i代表注入电压信号u_αsi的角频率；

载波信号注入后，电机坐标下的电压方程为：

$\left[\begin{array}{c}{u}_{\mathrm{\&alpha;}s}\\ u_{\mathrm{\&beta;}s}\end{array}\right]=U_{se}\left[\begin{array}{c}\cos \left({\mathrm{\&omega;}}_{r}t\right)\\ -\sin \left({\mathrm{\&omega;}}_{r}t\right)\end{array}\right]+U_{si}\left[\begin{array}{c}\cos \left({\mathrm{\&omega;}}_{i}t\right)\\ -\sin \left({\mathrm{\&omega;}}_{i}t\right)\end{array}\right]=U_{se}{e}^{{\mathrm{j\&omega;}}_{r}t}+U_{si}{e}^{{\mathrm{j\&omega;}}_{i}t}---\left(8\right)$

式中，U_se代表正序电流幅值，ω_r代表转子角频率；

在此高频电压注入下，产生的电流将由三部分组成：第一部分是与注入的电压旋转方向相同的正序电流，第二部分是与旋转电压方向相反的负序电流，第三部分是由三相绕组不对称产生的零序电流，电流响应可以表示为：

$i_{\mathrm{\&alpha;}\mathrm{\&beta;}si}=i_{ip}{e}^{j({\mathrm{\&theta;}}_i-\mathrm{\&pi;}/2)}+i_{in}{e}^{j(2{\mathrm{\&theta;}}_r-{\mathrm{\&theta;}}_i+\mathrm{\&pi;}/2)}---\left(9\right)$

其中，

式中，θ_r为d轴领先A相相轴的空间电角度，θ_i代表该注入电压信号的角频率为ω_i，i_in代表电流负序的幅值，U_si代表在静止坐标系上注入高频旋转电压的幅值，ω_i代表注入电压信号的角频率，L代表平均电感，ΔL代表空间调制电感；

从公式(8)中得出，只有高频响应电流的负序成分中含有转子位置信息，通过滤波器将电源产生的频率成分和正序电流分量滤除，然后用外差法得出转子位置的误差角度θ_e，再利用全维观测器提取出转子位置信息。

8.根据权利要求7所述的一种基于速度环模糊控制和高频注入法的无传感器控制系统，其特征在于，所述外差计算器子模块中外差法运算包括将公式(9)中i_αi、i_βi分别乘以和然后作差：

${\mathrm{\&theta;}}_e=i_{\mathrm{\&beta;}i}\cos (2\widehat{{\mathrm{\&theta;}}_r}-{\mathrm{\&omega;}}_{i}t)-i_{\mathrm{\&alpha;}i}\sin (2\widehat{{\mathrm{\&theta;}}_r}-{\mathrm{\&omega;}}_{i}t)=I_{i0}\sin \&lsqb;2({\mathrm{\&omega;}}_{i}t-\widehat{{\mathrm{\&theta;}}_r})\&rsqb;+I_{i1}\sin \&lsqb;2({\mathrm{\&theta;}}_r-\widehat{{\mathrm{\&theta;}}_r})\&rsqb;---\left(10\right)$

式中，θ_r为d轴领先A相相轴的空间电角度，代表高频电压注入法获得的转子初判角，ω_i代表注入电压信号的角频率；

其中，第一项为含电流的高频分量，第二项为仅含转子位置的信息，通过低通滤波可得转子位置的误差信号，从而：

${\mathrm{\&theta;}}_e=I_{i1}sin\&lsqb;2({\mathrm{\&theta;}}_r-\widehat{{\mathrm{\&theta;}}_r})\&rsqb;---\left(11\right)$

在角度误差很小的情况下，

$\mathrm{\&Delta;}\mathrm{\&theta;}\&ap;\frac{1}{2}sin\&lsqb;2({\mathrm{\&theta;}}_r-{\widehat{\mathrm{\&theta;}}}_r)\&rsqb;.---\left(12\right)$

9.根据权利要求4所述的一种基于速度环模糊控制和高频注入法的无传感器控制系统，其特征在于，所述全维观测器中转子转速的估计值通过以下公式求得：

交流永磁同步电机的运动方程可以表示为：

$\frac{{\mathrm{d\&omega;}}_r}{dt}=\frac{T_e-T_L}{J}=a---\left(13\right)$

式中,J为转动惯量，T_L代表负载转矩；

电机转子在一个采样周期T_s上的角位移公式是:

$\mathrm{\&theta;}(t_0+T)-\mathrm{\&theta;}\left(t_0\right)=\mathrm{\&omega;}\left(t_0\right)T+\frac{1}{2}{\mathrm{aT}}^2---\left(14\right)$

式中，t₀代表转子开始时间，T代表转子经过时间；

采样周期很短，上式表示为:

$\stackrel{\&CenterDot;}{{\mathrm{\&theta;}}_r}={\mathrm{\&omega;}}_r+\frac{1}{2}aT---\left(15\right)$

式中，ω_r代表转子角速度；

由式(13)和(15)可以得到：

$\stackrel{\&CenterDot;}{{\mathrm{\&theta;}}_r}={\mathrm{\&omega;}}_r+\frac{1}{2}\frac{(T_e-T_L)}{J}---\left(16\right)$

电机系统中负载变换缓慢，所以可认为：

$\frac{{\mathrm{dT}}_L}{dt}=0---\left(17\right)$

将式(13)、(16)及(17)改写成矩阵形式:

$\left[\begin{array}{c}\stackrel{\&CenterDot;}{\widehat{{\mathrm{\&theta;}}_r}}\\ \stackrel{\&CenterDot;}{\widehat{{\mathrm{\&omega;}}_r}}\\ \stackrel{\&CenterDot;}{\widehat{T_L}}\end{array}\right]=\left[\begin{array}{ccc}0& 1& -\frac{T}{2J}\\ 0& 0& -\frac{1}{J}\\ 0& 0& 0\end{array}\right]\left[\begin{array}{c}\widehat{{\mathrm{\&theta;}}_r}\\ \widehat{{\mathrm{\&omega;}}_r}\\ \widehat{T_L}\end{array}\right]+\left[\begin{array}{c}\frac{T}{2J}\\ \frac{1}{J}\\ 0\end{array}\right]\widehat{T_e}+\left[\begin{array}{c}{l}_1\\ l_2\\ l_3\end{array}\right]({\mathrm{\&theta;}}_r-\widehat{{\mathrm{\&theta;}}_r})---\left(18\right)$

式中，l₁、l₂和l₃三个表示的是在观测器中的增益值；

通过极点配置的方式来设置合理的全维观测器，离散化后的全维观测器方程为：

$\left[\begin{array}{c}{\widehat{\mathrm{\&theta;}}}_r(k+1)\\ {\widehat{\mathrm{\&omega;}}}_r(k+1)\\ {\hat{T}}_L(k+1)\end{array}\right]=\left[\begin{array}{ccc}0& T& -\frac{T^2}{2J}\\ 0& 0& -\frac{T}{J}\\ 0& 0& 0\end{array}\right]\left[\begin{array}{c}{\widehat{\mathrm{\&theta;}}}_r\left(k\right)\\ {\widehat{\mathrm{\&omega;}}}_r\left(k\right)\\ {\hat{T}}_L\left(k\right)\end{array}\right]+\left[\begin{array}{c}\frac{T}{2J}\\ \frac{1}{J}\\ 0\end{array}\right]\widehat{T_e}\left(k\right)+\left[\begin{array}{c}{l}_1\\ l_2\\ l_3\end{array}\right]T\left({\mathrm{\&theta;}}_r\right(k)-\widehat{{\mathrm{\&theta;}}_r}(k\left)\right)---\left(19\right)$

10.根据权利要求4所述的一种基于速度环模糊控制和高频注入法的无传感器控制系统，其特征在于，所述高频信号注入模块向两相静止直角坐标系α-β中注入高频旋转电压信号u_asi和u_βsi为：

u_asi＝v_si sinω_it (20)

u_βsi＝v_si cosω_it (21)

其中，v_si是注入的高频电压信号的幅值,ω_i为注入的高频电压信号的角频率。