4.2 传统观测器与抗扰技术
在永磁同步电机（PMSM）高性能控制系统中，为实现对转矩、转速及位置的精确闭环控制，必须获取准确的状态反馈信息。物理传感器（如电流传感器、转速编码器）虽能直接测量，但存在成本、可靠性及安装限制。此外，系统运行中不可避免的负载扰动、模型参数变化等不确定性因素会直接影响控制精度与动态性能。传统观测器与抗扰技术为此提供了有效的解决方案，它们基于系统的数学模型，利用可直接测量的电气量（如端电压、相电流）来实时估算无法直接测量或变化缓慢的状态量（如磁链、扰动转矩），并通过前馈补偿等方式提升系统的鲁棒性。本节将系统阐述以龙伯格观测器、滑模观测器为代表的传统状态观测器设计方法，并深入分析基于扰动观测的负载转矩补偿等关键抗扰技术。
4.2.1 状态观测的基本原理与龙伯格观测器
状态观测器是一种基于系统动态模型的软测量技术。对于一个线性系统，其状态空间方程可描述为：
x˙=Ax+Buy=Cx egin{aligned} dot{mathbf{x}} &= mathbf{A}mathbf{x} + mathbf{B}mathbf{u} \ mathbf{y} &= mathbf{C}mathbf{x} end{aligned}x˙y=Ax+Bu=Cx
其中，xmathbf{x}x为状态变量，umathbf{u}u为输入，ymathbf{y}y为输出。若系统(A,C)(mathbf{A}, mathbf{C})(A,C)可观测，则可构造一个与原系统并行运行的动态系统——观测器，其输入为原系统的实际输入umathbf{u}u和输出ymathbf{y}y，输出为状态估计值x^hat{mathbf{x}}x^。
龙伯格观测器是最经典的全维状态观测器，其结构如下：
x^˙=Ax^+Bu+L(y−y^) dot{hat{mathbf{x}}} = mathbf{A}hat{mathbf{x}} + mathbf{B}mathbf{u} + mathbf{L}(mathbf{y} - hat{mathbf{y}})x^˙=Ax^+Bu+L(y−y^)
其中，y^=Cx^hat{mathbf{y}} = mathbf{C}hat{mathbf{x}}y^=Cx^，Lmathbf{L}L为观测器增益矩阵。观测误差x~=x−x^ ilde{mathbf{x}} = mathbf{x} - hat{mathbf{x}}x~=x−x^的动态方程为x~˙=(A−LC)x~dot{ ilde{mathbf{x}}} = (mathbf{A} - mathbf{L}mathbf{C}) ilde{mathbf{x}}x~˙=(A−LC)x~。通过合理设计Lmathbf{L}L，可将矩阵(A−LC)(mathbf{A} - mathbf{L}mathbf{C})(A−LC)的特征值配置在复平面左半平面的期望位置，从而确保观测误差以指数速度收敛至零。收敛速度通常快于原系统动态，这是观测器设计的一般原则[1]。
在PMSM控制中，龙伯格观测器常被用于转子磁链或反电动势的观测，进而估算转速与位置。以在静止αβalphaetaαβ坐标系下为例，将反电动势e=[eα,eβ]Tmathbf{e} = [e_alpha, e_eta]^Te=[eα,eβ]T扩展为状态变量，构建包含电流与反电动势的状态方程。通过设计观测器增益，可实现对反电动势e^hat{mathbf{e}}e^的渐近跟踪。由于在稳态下，反电动势与转子位置角θe heta_eθe直接相关（eα∝−ωesinθe,eβ∝ωecosθee_alpha propto -omega_e sin heta_e, e_eta propto omega_e cos heta_eeα∝−ωesinθe,eβ∝ωecosθe），因此可通过θ^e=arctan(−e^α/e^β)hat{ heta}_e = arctan(-hat{e}_alpha / hat{e}_eta)θ^e=arctan(−e^α/e^β)估算位置，并通过微分或计算得到转速ω^ehat{omega}_eω^e。标准龙伯格观测器对模型参数准确性较为敏感，参数失配会导致稳态观测误差。
4.2.2 滑模观测器及其改进
滑模观测器 是一种基于变结构控制理论的非线性观测器，因其对参数摄动和外部扰动具有强鲁棒性而被广泛应用于PMSM无位置传感器控制。其基本思想是设计一个不连续的反馈控制律，迫使系统状态轨迹在有限时间内到达并保持在预定的滑模面上运动，此时等效控制量包含了待观测信息。
1. 传统滑模观测器设计
以PMSM在两相静止坐标系下的电流方程为例：
ddt[iαiβ]=1Ls([vαvβ]−Rs[iαiβ]−[eαeβ]) frac{d}{dt}egin{bmatrix} i_alpha \ i_eta end{bmatrix} = frac{1}{L_s} left( egin{bmatrix} v_alpha \ v_eta end{bmatrix} - R_s egin{bmatrix} i_alpha \ i_eta end{bmatrix} - egin{bmatrix} e_alpha \ e_eta end{bmatrix}
ight)dtd[iαiβ]=Ls1([vαvβ]−Rs[iαiβ]−[eαeβ])
构造滑模电流观测器：
ddt[i^αi^β]=1Ls([vαvβ]−Rs[i^αi^β]−[zαzβ]) frac{d}{dt}egin{bmatrix} hat{i}_alpha \ hat{i}_eta end{bmatrix} = frac{1}{L_s} left( egin{bmatrix} v_alpha \ v_eta end{bmatrix} -