Huapu-CP: From Knowledge Graphs to a Data Central-Platform
-
摘要: 针对碎片化的各姓氏家谱数据, 华谱系统通过构建家谱知识图谱的数据中台, 能够解决数据孤岛、烟囱式开发等问题. “数据中台”是一个源自国内的新近技术概念, 在华谱系统建设中, 我们通过家谱知识图谱的构建和应用, 对这个概念进行了正式定义. 基于这个定义和对应的7项核心功能, 本文提出一种用于家谱数据分析的数据中台建设架构Huapu-CP (华谱系统), 并通过该架构详细介绍面向家谱领域的数据中台核心技术, 分析数据中台构建的关键问题.Abstract: With fragmented family tree data, Huapu aims to solve the problems of data islands and chimney development by building a data central-platform based on genealogical knowledge graphs. Data central-platform is a new technology that has recently emerged in China. During the construction of the Huapu system, we started with knowledge graph construction and analytics, and then formulated a formal definition for this technology concept. Based on this definition and its corresponding 7 key components, this paper proposes a Huapu-CP (or Huapu for short) framework for genealogical data analytics and presents the core techniques and key challenges in building such a genealogical data platform.
-
Key words:
- Genealogy /
- data central-platform /
- data governance /
- knowledge graph
-
感应电机交流传动系统广泛应用于各工业领域, 其中新能源汽车等场合普遍采用转子磁链定向的控制方式, 获得精确的转子磁链是实现交流传动系统高性能控制的关键[1].矢量控制策略是当前应用最为广泛的感应电机控制方法, 虽然采用矢量控制的感应电机交流传动系统具有很好的控制性能, 但在设计控制器过程中需要精确的电机参数来实现定、转子控制的解耦.无论是感应电机转子磁链定向矢量控制, 还是其他的非线性控制策略, 都需要转子磁链矢量的幅值和相位[2].转子磁链矢量的检测和获取方法分为直接法和间接法.直接法是在感应电机定子内表面装贴霍尔元件或者在电机槽内埋设探测线圈等直接检测转子磁链, 但由于工艺和技术难度较大, 实际的矢量控制系统中不适用直接法[3].间接法是检测感应电机的定子电压、电流及转速等容易获得的物理量, 利用转子磁链观测模型, 实时计算转子磁链的幅值和相位[4].而由于观测模型不够精确, 控制系统中的延迟问题以及电机参数变化的影响等, 使提高转子磁链观测精度成为提高交流传动系统控制性能的关键问题之一.
为了提高转子磁链观测的精度, 很多专家学者不懈努力, 进行了深入的研究.提出了电压模型法和电流模型法的转子磁链观测器、U-I法磁链观测器、全阶磁链观测器、扩展卡尔曼滤波器、自适应观测器和滑模观测器等方法[5-13].电压模型法模型结构简单, 计算过程是纯积分, 其估算结果受积分初值和输入信号的直流偏移影响很大, 导致结果存在误差, 且电压模型法依赖电机的定子电阻参数, 其受温度等因素影响较大, 也会产生估算误差[5].电流模型法依赖电机的定子电流和转速参数, 同时磁链估算过程需要转子参数, 鲁棒性差[6].电压模型法和电流模型法都是基于开环算法的磁链估算, 观测精度受限.文献[7]提出了一种基于改进U-I法的磁链观测方法, U-I法不需要转速和转子参数, 具有较好的鲁棒性, 但其对定子电阻的摄动较敏感且存在积分漂移, 造成观测误差.文献[8]提出了一种全阶磁链观测器设计方法, 通过观测得到的电机转速等相关参数, 设计合理的估算模型, 得到磁链观测值, 设计过程的极点配置以及受电机参数影响较大, 限制了其观测的精度.文献[9]提出了一种基于扩展Kalman滤波器的转子磁链观测方法, 能够有效减少噪声对磁链观测精度的影响, 但其对参数变化的敏感性及估算过程需要大量的数学计算, 限制了其在实际工程中应用.文献[10]提出了一种基于模型参考自适应系统的自适应磁链观测方法, 改善磁链观测的精度, 但其受系统参数影响较大.文献[12]提出了一种非线性滑模磁链观测方法, 具有较强的抗干扰性, 但其存在的抖振问题无法消除, 极大地限制了其使用范围.
转子磁链幅值和相位的准确估计是构建感应电机交流传动矢量控制系统的关键环节.磁链幅值估计实现系统的磁链控制, 转子位置观测实现矢量控制系统的坐标变换, 从而完成感应电机转矩和励磁控制的解耦[14].中立型系统理论是基于中立型延迟系统的一种理论, 而中立型延迟系统是一种能够精确描述延迟系统的模型, 模型中既包括状态延迟, 也包括状态微分延迟, 使得对延迟系统的描述更加精确[15].本文将中立型系统理论应用到感应电机控制系统中, 解决由于观测模型不够精确、系统控制中的延迟问题以及电机参数变化的影响, 导致转子磁链观测精度不高的问题, 实现系统的高性能控制.
1. 感应电机数学模型及中立型系统理论原理
1.1 感应电机数学模型
感应电机的数学模型具有高阶、非线性、强耦合等特征.矢量控制系统建立在感应电机的动态模型上, 在同步旋转$ M$-$T$坐标系中, 当电机转子磁链矢量与$T$轴重合时, 即为按转子磁链定向.通过同步旋转坐标系按转子磁链方向定向, 实现感应电机转矩和磁通的解耦控制.
本文基于感应电机在$ M$-$T$坐标系下的状态方程, 提出一种电机转子磁链观测方法.选择电机的定子电流和转子磁链为状态变量, 以电压变量为输入, 转子磁链为输出, 建立感应电机数学模型[16].
$$ \begin{eqnarray} \left\{ \begin{aligned} &\dot{\Psi}_r = -\lambda_1\Psi_r + \lambda_2i_{sm} \\ &\dot{i}_{sm} = \lambda_3\Psi_r - \lambda_5i_{sm} + \lambda_6i_{st} + \lambda_7u_{sm} \\ &\dot{i}_{st} = -\lambda_4\Psi_r - \lambda_5i_{st} - \lambda_6i_{sm} + \lambda_7u_{st} \end{aligned} \right. \end{eqnarray} $$ (1) $$ \lambda_1 = \frac{1}{T_r}, \lambda_2 = \frac{L_m}{T_r}, \lambda_6=\omega_s, \lambda_7 = \frac{1}{\sigma L_s} $$ $$ \lambda_3 = \frac{L_m}{\sigma L_s L_r T_r}, \lambda_4 = \frac{L_m}{\sigma L_s T_r}, \lambda_5 = \frac{R_s L^{2}_r + R_r L^{2}_m}{\sigma L_s L^{2}_r} $$ 式中, $u_{sm}$和$u_{st}$为定子电压; $i_{sm}$和$i_{st}$为定子电流; $R_s$和$R_r$分别为定子和转子绕组电阻; $L_m$为定子和转子同轴等效绕组互感; $L_s$和$L_r$分别为定子和转子等效两相绕组自感; $\Psi_r$为转子磁链矢量; $n_p$为极对数; $T_L$为负载转矩; $T_e$为电磁转矩; $J$为电机的转动惯量; $\sigma$为电机的漏磁系数, $\sigma=1-L^{2}_m/(L_sL_r)$; $T_r$为转子电磁时间常数, $T_r=L_r/R_r$; $\omega$为转子转动角速度; $\omega_s$为同步角速度, 感应电机矢量控制辅助方程为
$$ \begin{eqnarray} \left\{ \begin{aligned} &\dot{\omega} = \frac{n^{2}_p L_m}{JL_r}i_{st}\Psi_r - \frac{n_p}{J}T_L \\ &T_e = \frac{n_p L_m}{L_r}i_{st}\Psi_r \\ &\omega_s = \omega + \frac{L_m}{T_r\Psi_r}i_{st}(t) \end{aligned} \right. \end{eqnarray} $$ (2) 1.2 中立型系统理论
中立型系统理论是一种针对解决工程实践中延迟问题的理论, 中立型系统模型对系统描述的精确度要优于当前其他的建模理论.标称延迟系统方程为:
$$ \begin{eqnarray*}\dot{{\pmb x}}(t) = f({\pmb x}(t), {\pmb X}(t-d), {\pmb U}(t))\end{eqnarray*} $$ 而中立型延迟系统方程为:
$$ \begin{eqnarray*}\dot{{\pmb x}}(t) - G\dot{{\pmb x}}(t-d) = f({\pmb x}(t), {\pmb X}(t-d), {\pmb U}(t))\end{eqnarray*} $$ 因为系统含有微分差分算子$\dot{D}(t, {\pmb x}_t) = \dot{{\pmb x}}(t) - G\dot{{\pmb x}}(t-d)$, 使得标称延迟系统的许多理论成果无法简单地推广到中立型延迟系统中去.
中立型延迟系统模型方程[15]:
$$ \begin{eqnarray} \left\{ \begin{aligned} &\dot{{\pmb x}}(t) - G\dot{{\pmb x}}(t-d) = A{\pmb x}(t) + C {\pmb x}(t-d) + \\&f(t, {\pmb x}(t), {\pmb X}(t-d)) + B {\pmb u}(t) \\ &{\pmb y}(t) = D {\pmb x}(t) \\ &{\pmb x}(t) = {\pmb \phi}(t), t\in[-d, 0] \end{aligned} \right. \end{eqnarray} $$ (3) 其中, ${\pmb x}(t)\in {\bf R}^{n}$为状态向量; $d>0$为延时时间; $A, C, G$和$B$为维数适当的常数矩阵, 且满足$\parallel G\parallel < 1$; ${\pmb \phi}(t)\in([-d, 0], {\bf R}^{n})$为向量初值函数; $f\in C([0, +\infty], {\bf R}^{n})$为不可观测的非线性不确定扰动, 且满足:
$$ \begin{eqnarray*}\left\{ \begin{aligned} &f(t, 0, 0) = 0 \\ &\parallel f(t, {\pmb x}(t-d))\parallel \leq \alpha\parallel {\pmb x} \parallel + \beta\parallel {\pmb x}(t-d)\parallel \end{aligned} \right.\end{eqnarray*} $$ 其中, $\alpha, \beta$为已知常数.
1.3 中立型系统理论在感应电机中的应用
电机控制系统中, 由于电力电子器件的存在, 以及控制器响应时间等, 控制信号会产生延迟现象, 出现开关动作、信号响应等和系统模态不同步, 这种现象称为异步切换[17].当开关频率较高时, 异步切换现象对系统性能的影响可以忽略.而在大功率低开关频率条件下, 异步现象会造成电机的电流畸变、发热等问题.中立型系统在建模时不需要考虑异步切换现象, 故本文将中立型系统理论应用到感应电机系统建模中, 建立感应电机中立型延迟系统数学模型.
$$ \begin{equation} \dot{{\pmb x}}(t) = \begin{bmatrix} A &B \end{bmatrix}\begin{bmatrix} {\pmb x}(t) \\{\pmb u}(t) \end{bmatrix} \end{equation} $$ (4) 令式(3)中系数矩阵$A = C$.其中
$$ \begin{align*} &\dot{{\pmb x}}(t) = \begin{bmatrix} \dot{\Psi}_r(t) \\\dot{i}_{sm}(t) \\\dot{i}_{st}(t) \end{bmatrix}, \quad {\pmb x}(t) = \begin{bmatrix} \Psi_r(t) \\ i_{sm}(t) \\i_{st}(t) \end{bmatrix}\\&{\pmb u}(t) = \begin{bmatrix} u_{sm}(t) \\u_{st}(t) \end{bmatrix}, \quad A = \begin{bmatrix} -\lambda_1 & \lambda_2 & 0 \\ \lambda_3 & -\lambda_5 & \lambda_6 \\ -\lambda_4 & -\lambda_6 & -\lambda_5 \end{bmatrix} \\&B = \begin{bmatrix} 0 & 0 \\\lambda_7 & 0 \\0 & \lambda_7 \end{bmatrix}\end{align*} $$ 在电机控制系统中, 由于异步切换现象的存在, 状态滞后对系统的影响具有不确定性.但在一定范围内, 影响程度取决于对应的系数矩阵元素的大小, 本文定义$\mu(\chi)$为影响因子, 表示系统中状态滞后对系统影响程度的度量的物理量, 其中$\chi$为对应的系数矩阵.这里:
$$ \begin{equation*}\mu(\chi) = \frac{1}{\parallel \chi\parallel_\infty}\end{equation*} $$ 考虑状态滞后对电机控制系统的影响, 得到感应电机标称延迟系统方程(5), 其中$\tau$为总延时时间.
$$ \begin{equation} \dot{{\pmb x}}(t) = \begin{bmatrix} (1-\mu)A & \mu C & B \end{bmatrix}\begin{bmatrix} {\pmb x}(t) \\{\pmb x}(t-\tau) \\{\pmb u}(t) \end{bmatrix} \end{equation} $$ (5) 式(5)由莱布尼兹公式得:
$$ \begin{align} &\left[ {\begin{array}{*{20}{c}} 1&-\mu C\int \end{array}} \right]\begin{bmatrix} \dot{{\pmb x}}(t) \\\dot{{\pmb x}}(t-d) \end{bmatrix} =\nonumber\\& \qquad \left[ {\begin{array}{*{20}{c}} (1-\mu)A & \mu C & B \end{array}} \right]\left[ \begin{array}{c} {\pmb x}(t) \\{\pmb x}(t-d) \\{\pmb u}(t) \end{array} \right] \end{align} $$ (6) 式中, $d$为系统状态延时时间, $d=\tau/2$.
表达式(6)中含有积分项, 这里用$N$个小直角梯形的面积之和极限逼近积分项的值, $N\in Z$, 又由数值均值定理, 得到表达式(7).
$$ \begin{align} &\int^{-\frac{\tau}{2}}_0\dot{{\pmb x}}(t-d+\theta){\rm d}\theta\approx\nonumber\\& \qquad\frac{d}{N}\sum\limits_{i=1}^N\dot{{\pmb x}}\left(t-d-\frac{(2i-1)\cdot\tau}{4N}\right) \end{align} $$ (7) 对于任意的$d$, 存在$N\geq M, M\in Z$, 使不等式$\mu\frac{d}{N}\leq\mu(\chi)$成立.
由上式得到感应电机中立型转子磁链模型方程(8).
$$ \begin{align} &\left[ {\begin{array}{*{20}{c}} 1 & -\mu C \cdot \frac{d}{N}\end{array}} \right]\begin{bmatrix} \dot{{\pmb x}}(t) \\\sum\limits_{i=1}^N\dot{{\pmb x}}(t-d-\varepsilon) \end{bmatrix} =\nonumber\\& \qquad \left[ {\begin{array}{*{20}{c}} (1-\mu)A & \mu C & B \end{array}} \right]\left[ \begin{array}{c} {\pmb x}(t) \\{\pmb x}(t-d) \\{\pmb u}(t) \end{array} \right] \end{align} $$ (8) 其中, $\varepsilon=\frac{(2i-1)d}{2N}$.
2. 中立型磁链观测器设计
考虑标准型中立型系统[18]:
$$ \begin{equation} \left\{\begin{aligned} &\dot{{\pmb x}}(t)-G\dot{{\pmb x}}(t-d)=A{\pmb x}(t)+C{\pmb x}(t-d)\\ &{\pmb y}(t)=D{\pmb x}(t) \end{aligned} \right. \end{equation} $$ (9) 式中
$$ \begin{align*} A =\, & (1-\mu)\times\\ &\left[{\begin{array}{*{20}{c}} -\frac{1}{T_r}&\frac{L_m}{T_r}&0 \\ \frac{L_m}{\sigma L_sL_rT_r}&-\frac{R_sL^{2}_r+R_rL^{2}_m}{\sigma L_sL^{2}_r}&\omega_s\\ -\frac{L_m}{\sigma L_sL_r}\omega&-\omega_s&-\frac{R_sL^{2}_r+R_rL^{2}_m}{\sigma L_sL^{2}_r}\end{array}} \right]\\ C =\, & \mu\left[ {\begin{array}{*{20}{c}} -\frac{1}{T_r}&\frac{L_m}{T_r}&0 \\ \frac{L_m}{\sigma L_sL_rT_r}&-\frac{R_sL^{2}_r+R_rL^{2}_m}{\sigma L_sL^{2}_r}&\omega_s\\ -\frac{L_m}{\sigma L_sL_r}\omega&-\omega_s&-\frac{R_sL^{2}_r+R_rL^{2}_m}{\sigma L_sL^{2}_r}\end{array}} \right]\end{align*} $$ $$ \begin{align*} G =\, & \mu\frac{d}{N}\times\\ &\left[ {\begin{array}{*{20}{c}} -\frac{1}{T_r}&\frac{L_m}{T_r}&0 \\ \frac{L_m}{\sigma L_sL_rT_r}&-\frac{R_sL^{2}_r+R_rL^{2}_m}{\sigma L_sL^{2}_r}&\omega_s\\ -\frac{L_m}{\sigma L_sL_r}\omega&-\omega_s&-\frac{R_sL^{2}_r+R_rL^{2}_m}{\sigma L_sL^{2}_r}\end{array}} \right]\\ D =\, & \begin{bmatrix} 1&0&0 \end{bmatrix}\end{align*} $$ 设计如下状态观测器:
$$ \begin{align} &\hat{\dot{{\pmb x}}}(t)-G\hat{\dot{{\pmb x}}}(t-d) = A\hat{{\pmb x}}(t)+C\hat{{\pmb x}}(t-d)+\nonumber\\& \qquad L({\pmb y}(t)-D\hat{{\pmb x}}(t)) \end{align} $$ (10) 使得误差动态系统方程(11)渐近稳定.
$$ \begin{equation} \dot{{\pmb e}}(t)-G\dot{{\pmb e}}(t-d) = H {\pmb e}(t) + C {\pmb e}(t-d) \end{equation} $$ (11) 式中, $H=A-LD, {\pmb e}(t)={\pmb x}(t)-\hat{{\pmb x}}(t)$为误差向量, $\hat{{\pmb x}}(t)\in R^{n}$是观测状态, $L$为$n\times q$阶观测器增益矩阵.
设误差动态系统方程的初值条件为
$$ \left\{ \begin{aligned} &{\pmb e}(t_r)={\pmb \varphi}(t_r), -d\leq t_r\leq0 \\ &{\pmb \varphi}(t_r)\in C([-d, 0], R^{n}) \end{aligned} \right. $$ 若存在正定阵$P, Q, R$和矩阵$K$, 使得线性矩阵不等式(12)成立, 则误差动态系统方程(11)渐近稳定.
$$ \begin{equation}\begin{aligned} \begin{bmatrix} 2PA-2PLD+Q+R & D^{\rm T}L^{\rm T}PG-A^{\rm T}PG & PC\\ G^{\rm T}PLD-G^{\rm T}PA & -Q & G^{\rm T}PC \\ C^{\rm T}P & C^{\rm T}PG & -R\end{bmatrix}<0 \end{aligned}\end{equation} $$ (12) 令$K=PL$, 则线性矩阵不等式写为
$$ \begin{equation} \begin{aligned} \begin{bmatrix} 2PA-2KD+Q+R & D^{\rm T}K^{\rm T}G-A^{\rm T}PG & PC\\ G^{\rm T}KD-G^{\rm T}PA & -Q & G^{\rm T}PC \\ C^{\rm T}P & C^{\rm T}PG & -R\end{bmatrix} <0 \end{aligned} \end{equation} $$ (13) 即若存在正定阵$P, Q, R$, 满足线性矩阵不等式(13), 则存在状态观测器(10), 解得观测器增益矩阵$L=P^{-1}K$.
这里以感应电机转子磁链为观测目标, 由感应电机中立型系统方程(8), 设计中立型转子磁链观测器模型, 如图 1.其中系数$A_{11}, $ $A_{12}, $ $C_{11}, $ $C_{12}, $ $G_{11}, $ $G_{12}$为系统方程系数矩阵$A, C, G$的对应元素.
图 1 M-T坐标系下中立型转子磁链观测模型Fig. 1 Observation model of neutral rotor flux linkage in M-T coordinate system3. 稳定性分析
已知状态观测器方程(10), 假设观测器增益系数$L$已知, 证明误差动态系统方程(11)是渐近稳定的.
为此本文提出了一种基于线性矩阵不等式(Linear matrix inequality, LMI)的Lyapunov泛函[19-22], 用来证明方程(11)的稳定性.
$$ \begin{align} V(t, {\pmb e})=\, &[{\pmb e}(t)-G{\pmb e}(t-d)]^{\rm T}P[{\pmb e}(t)-G{\pmb e}(t-d)]+\nonumber\\&\int^{t}_{t-d}{\pmb e}^{\rm T}(\rho)Q{\pmb e}(\rho){\rm d}\rho+\int^{t}_{t-d}{\pmb e}^{\rm T}(\rho)R{\pmb e}(\rho){\rm d}\rho \end{align} $$ (14) 其中, $P, Q, R$为正定对称矩阵.
令${\pmb\eta}_1={\pmb e}(t), {\pmb\eta}_2={\pmb\eta}_3={\pmb e}(t-d)$, 则表达式为
$$ \begin{align} V(t, {\pmb e})=\, &({\pmb\eta}_1-G{\pmb\eta}_2)^{\rm T}P({\pmb\eta}_1-G{\pmb\eta}_2)+\nonumber\\&\int^{t}_{t-d}{\pmb e}^{\rm T}(\rho)Q{\pmb e}(\rho){\rm d}\rho+\int^{t}_{t-d}{\pmb e}^{\rm T}(\rho)R{\pmb e}(\rho){\rm d}\rho \end{align} $$ (15) 对式(15)求导, 得
$$ \begin{align} \dot{V}(t, {\pmb e})=\, &2({\pmb\eta}_1-G{\pmb\eta}_2)^{\rm T}P(\dot{{\pmb\eta}_1}-G\dot{{\pmb\eta}_2})+\nonumber\\&{\pmb\eta}_1^{\rm T}Q{\pmb\eta}_1-{\pmb\eta}_2^{\rm T}Q{\pmb\eta}_2+{\pmb\eta}_1^{\rm T}R{\pmb\eta}_1-{\pmb\eta}_3^{\rm T}R{\pmb\eta}_3 \end{align} $$ (16) 由中立型系统方程式(9), 式(16)可以表达为
$$ \begin{align} \dot{V}(t, {\pmb e})=\, &2({\pmb\eta}_1-G{\pmb\eta}_2)^{\rm T}P(H{\pmb\eta}_1-C{\pmb\eta}_2)+\nonumber\\&{\pmb\eta}_1^{\rm T}(Q+R){\pmb\eta}_1-{\pmb\eta}_2^{\rm T}Q{\pmb\eta}_2-{\pmb\eta}_3^{\rm T}R{\pmb\eta}_3 \end{align} $$ (17) $$ \begin{align} \dot{V}(t, {\pmb e})=\, &2{\pmb\eta}_1^{\rm T}PH{\pmb\eta}_1+2{\pmb\eta}_1^{\rm T}PC{\pmb\eta}_3-2{\pmb\eta}_2^{\rm T}G^{\rm T}PH{\pmb\eta}_1-\nonumber\\&2{\pmb\eta}_2^{\rm T}G^{\rm T}PC{\pmb\eta}_3+{\pmb\eta}_1^{\rm T}(Q+R){\pmb\eta}_1-\nonumber\\&{\pmb\eta}_2^{\rm T}Q{\pmb\eta}_2-{\pmb\eta}_3^{\rm T}R{\pmb\eta}_3 =\nonumber\\&{\pmb\eta}_1^{\rm T}(PH+H^{\rm T}P+Q+R){\pmb\eta}_1+\nonumber\\&2{\pmb\eta}_1^{\rm T}PC{\pmb\eta}_3-2{\pmb\eta}_2^{\rm T}G^{\rm T}PH{\pmb\eta}_1-\nonumber\\&2{\pmb\eta}_2^{\rm T}G^{\rm T}PC{\pmb\eta}_3-{\pmb\eta}_2^{\rm T}Q{\pmb\eta}_2-{\pmb\eta}_3^{\rm T}R{\pmb\eta}_3 \end{align} $$ (18) 将式(18)等价为线性矩阵不等式的形式为式(19).其中${\pmb\eta}(t)=\begin{bmatrix}{\pmb e}(t) & {\pmb e}(t-d) & {\pmb e}(t-d)\end{bmatrix}$, $H=A-LD$.令$K=PL$则式(19)表达为式(20).
若线性矩阵不等式(13)有解, 即存在正定阵$P, Q, R$和矩阵$K$, 得$\dot{V}(t, {\pmb e}) < 0$由Razumikhin型定理[20]得误差动态系统方程(11)渐近稳定.
$$ \begin{align} \dot{V}(t, {\pmb e})={\pmb\eta}^{\rm T}(t)\times \begin{bmatrix} PH+H^{\rm T}P+Q+R & -H^{\rm T}PG & PC\\ -G^{\rm T}PH & -Q & -G^{\rm T}PC \\C^{\rm T}P & -C^{\rm T}PG & -R \end{bmatrix}{\pmb\eta}(t) \end{align} $$ (19) $$ \begin{align} \dot{V}(t, {\pmb e})={\pmb\eta}^{\rm T}(t)\times \begin{bmatrix} 2PA-2KD+Q+R & D^{\rm T}K^{\rm T}G-A^{\rm T}PG & PC\\ G^{\rm T}KD-G^{\rm T}PA & -Q & G^{\rm T}PC \\C^{\rm T}P & C^{\rm T}PG & -R \end{bmatrix}{\pmb\eta}(t) \end{align} $$ (20) 另外, 这里假设式(1)中的感应电机定子和转子绕组互感和绕组电阻等电机参数在一定范围内发生变化时, 记为$L^{\ast}_s, L^{\ast}_r, R^{\ast}_s, R^{\ast}_r$等, 分别对应变化后的方程(8)系数矩阵$A^{\ast}, B^{\ast}, C^{\ast}, D^{\ast}$, 把变化后的系数矩阵带入线性矩阵不等式(13), 求解不等式有解, 即存在正定阵$P^{\ast}, Q^{\ast}, R^{\ast}$和矩阵$K^{\ast}$, 观测器增益矩阵$L^{\ast}=P^{\ast-1}K^{\ast}$, 得误差动态系统方程(11)依然渐近稳定.所以通过分析说明了在一定范围内, 所设计的中立型磁链观测器具有较好的鲁棒性.
4. 仿真及实验研究
针对感应电机中立型转子磁链观测模型, 利用Matlab/Simulink搭建仿真模型, 借助Simulink/S-函数对中立型转子磁链观测器增益矩阵进行求解, 并从实际应用的角度出发, 设计实验方案, 借助DSP电力电子与电气传动实验平台, 对中立型转子磁链观测器进行实验验证.仿真和实验使用的电机参数一样, 电机参数见表 1.
表 1 电机参数Table 1 Motor parameters参数 数值 额定功率$P_N$ /kW 4 额定电压$U_N$/V 380 额定频率$f_N$/Hz 50 额定电流$I_N$/A 8.8 定子电阻$R_s/{\rm \Omega}$ 1.405 定子电感$L_s$/mH 178 转子电阻$R_r/{\rm \Omega}$ 1.395 转子电感$L_r/$mH 178 定转子互感$L_m/$mH 172.2 极对数$n_p$ 2 额定转速$n/{\rm r}\cdot {\rm min}_{-1}$ 1 440 将表 1电机参数带入感应电机中立型延迟系统方程, 得到方程各项系数矩阵, 利用MATLAB中的LMI工具箱, 求解线性矩阵不等式(13), 解得正定矩阵$P, Q, R$以及矩阵$K$, 即可求出观测器增益矩阵$L$, 如下:
$$ \begin{align*}&P=\begin{bmatrix}19.0204&0.1979&0.0079\\ 0.1979&0.1770&0.0275\\ 0.0079&0.02575&0.0073\end{bmatrix}\\& Q=\begin{bmatrix}39.2067&0.7058&0.9499\\ 0.7058&19.4849&1.1232\\ 0.9499&1.1232&8.2866\end{bmatrix}\\& R=\begin{bmatrix}18.8541&-0.0330&-0.0008\\ -0.0330&20.6600&0.8838\\ -0.0008&0.8838&8.2839\end{bmatrix}\\& K=\begin{bmatrix}-6.1763\\ -8.8302\\ 13.7074\end{bmatrix}, L=\begin{bmatrix}0.0063\\ -0.8371\\ 5.0218\end{bmatrix}\times10^{3}\end{align*} $$ 因为存在正定矩阵$P, Q, R$, 已知$C=\begin{bmatrix}1 & 0 & 0\end{bmatrix}$, 因此, 矩阵$L$中的元素$L_1$即是图 1中的观测器增益.
上述所求正定矩阵及$L$是实时变化的, 上面是电机达到稳定后的一组解, 证明了感应电机中立型延迟系统模型的合理性.
本文利用S-函数, 模拟微处理器运行, 建立基于中立型转子磁链观测器的三相感应电机矢量控制系统仿真模型, 如图 2.仿真时间设为2.0 s.仿真和实验中电流控制器均采用的是PI控制器.
图 2 基于中立型转子磁链观测器的感应电机矢量控制系统仿真模型Fig. 2 Simulation model of vector control system of induction motor based on the observer of neutral rotor flux linkage仿真中转子磁链的真实值由电机模型给出.给定电机转速为500 r/m, 给定转子磁链为0.96 Wb, 电流采样控制频率为4 kHz, 开关频率为500 Hz, 即开关周期为0.002 s, 系统总延时时间[2]为$\tau=1.5T_s=3$ ms.图 3是基于中立型转子磁链观测器的电机转速响应图, 从图 3中可以看出, 电机转速在0.15 s即达到给定值, 响应速度快, 且曲线平滑, 稳定性好; 图 4是转子磁链幅值响应图, 从图 4中可知, 转子磁链在小于0.03 s的时间内即达到稳定值, 响应迅速, 且稳定性好.
图 3 中立型磁链观测方法转速响应图Fig. 3 Speed response diagram of flux linkage observation of the neutral method
图 4 中立型磁链观测方法转子磁链幅值响应图Fig. 4 Magnitude response diagram of rotor flux linkage of flux linkage observation of the neutral method分别采用中立型转子磁链观测器、电压模型法观测器和二阶滑模方法观测器观测转子磁链, 并进行稳态观测对比.电机稳态运行时的转子磁链观测及观测误差波形如图 5和图 6所示.
图 5 感应电机稳态运行时磁链观测波形Fig. 5 Observation waveform of flux linkage of induction motor in steady state operation
图 6 中立型观测方法和二阶滑模观测方法观测磁链对比图Fig. 6 Observational flux linkage contrast diagram of neutral observation method and second order sliding-mode method从图 5中可知, 电机稳态运行时采用中立型转子磁链观测器的磁链观测曲线平滑, 波动小, 观测误差峰-峰值为0.02 Wb, 采用电压模型法的磁链观测曲线波动较大, 观测误差峰-峰值为0.1 Wb, 采用二阶滑模方法的磁链观测误差峰-峰值为0.06 Wb.从图 6中可以看出, 中立型观测方法和二阶滑模观测方法观测磁链的观测值都始终收敛于单位圆内, 满足预先给定值, 而采用中立型观测方法观测磁链的曲线密集程度明显优于采用二阶滑模方法的磁链观测, 从对应波形的放大图中看更加显著.所以本文提出的基于中立型系统理论的转子磁链观测方法的观测精度优于电压模型法和二阶滑模方法的观测精度.
考虑电机参数变化对磁链观测精度的影响, 对设计的中立型转子磁链观测器进行鲁棒性研究.图 7是感应电机转子电阻突变为真实值的1.5倍的磁链观测波形.图 8是感应电机转子电阻突变为真实值的0.5倍的磁链观测波形.
图 7 $R_r$值变为真实值的1.5倍磁链观测图Fig. 7 Observation diagram of flux linkage with $R_r$ value changed to 1.5-fold real value
图 8 $R_r$值变为真实值的0.5倍磁链观测图Fig. 8 Observation diagram of flux linkage with $R_r$ value changed to 0.5-fold real value由图 7和图 8可知, 当感应电机转子电阻发生突变后, 电机稳态运行时的磁链观测误差在小范围内波动, 有微小增大, 变化率均在5 %以内; 由图 9可知, 当感应电机转子电感发生突变后, 电机稳态运行时的磁链观测误差变化率低于8 %; 此外经验证把上述求得的正定矩阵$P, Q, R$以及矩阵$K$和$L$带入线性矩阵不等式(13), 把转子电阻变化后记为$R_r{'}$, 其作为变量, $R_r{'}\in[0.5R_r, 1.5R_r]$, 带入式(13), 线性矩阵不等式仍然成立, 说明了转子电阻变化对中立型转子磁链观测方法的影响很小, 证明了中立型磁链观测器具有较强的鲁棒性.
图 9 $L_r$值变为真实值的0.5倍磁链观测图Fig. 9 Observation diagram of flux linkage with $L_r$ value changed to 0.5-fold real value综合图 7~图 9以及对应的分析, 可以得到当电机参数发生突变时, 电机稳态运行时对应的磁链观测误差的变化范围见表 2.
表 2 磁链观测误差Table 2 Observation error of flux linkage参数变化量 $-0.5R_r$ $+0.5R_r$ $-0.5L_r$ $+0.5L_r$ $\Delta\Psi_{r\alpha}$变化量$10^{-3}$/kW + 0.9 + 1.0 +1.5 +1.6 $\Delta\Psi_{r\alpha}$变化率 $\leq 5 %$ $\leq 8 %$ 为了研究系统的动态性能, 初始转速给定值为500 r/m, 空载情况下, 在0.5 s时, 调节转速从500 r/m升高到1 000 r/m, 图 10是基于中立型观测方法的电机调速波形.电机以给定转速稳态运行时, 在0.5 s加35 N$\cdot $ m的负载, 图 11是负载阶跃的动态响应波形.
图 10 基于中立型观测方法的电机调速波形Fig. 10 Motor speed regulation waveform based on neutral observation method
图 11 基于中立型观测方法的负载阶跃的响应波形Fig. 11 Response waveform of load step based on neutral observation method从图 10 (a)中可以看出, 电机转速由500 r/m升高到1 000 r/m所用时间约为0.2 s, 说明系统具有较好的转速动态响应性能; 由图 10 (b)和(c)可知, 磁链观测响应速度快, 观测幅值不变, 观测误差很快收敛到稳态误差值0.02 Wb; 电机电流响应迅速, 说明系统具有很好的动态响应性能.
从图 11中可以看出, 当系统外部负载发生变化时, 转子磁链观测幅值在负载突变时有短时微小波动, 很快稳定到给定值; 电机转矩电流响应迅速, 曲线波动较小; 电机三相电流曲线平滑, 响应迅速; 电机实际转矩迅速稳定在给定值35 N$\cdot $ m, 表明系统对外部负载变化具有良好的抗干扰能力, 鲁棒性强, 即采用本文提出的中立型转子磁链观测器的系统具有很好的动态特性.
为了研究应用中立型磁链观测器的感应电机在低速条件下的性能, 这里给出电源频率为5 Hz时的仿真.给定电机转速为150 r/m, 空载条件下, 在0.5 s时让电机停止, 即使得转速降为0 r/m, 图 12是基于中立型观测方法的电机低速条件下的仿真波形.
图 12 基于中立型观测方法的感应电机低速运行波形Fig. 12 Low speed operation waveform of induction motor based on neutral observation method从图 12 (a)可以看出, 在低频状态下电机转速很快稳定到给定值150 r/m, 0.5 s电机制动, 转速很快降为0, 曲线平滑, 说明了低频下系统具有较好的转速动态响应性能; 由图 12 (b)可知, 低频条件下, 电机转速达到稳定值时, 磁链观测误差峰值约为0.02 Wb, 当电机低速制动时, 磁链观测误差开始减小, 且磁链观测误差收敛至0 Wb, 证明了中立型磁链观测方法在低速条件下的有效性.
为了研究延迟问题对磁链观测精度的影响, 这里考虑在仿真中加入延迟时间.因为上述仿真中采用电压模型法观测磁链时并未考虑延迟时间, 而中立型观测方法中已考虑延迟时间.当采用电压模型法的电机仿真中加入延时时间$d$时, 得到如图 13的磁链观测波形.由图 13可知, 当采用电压模型法的仿真中考虑延时时间$d$时, 转子磁链观测波动较大, 观测幅值波动高于$\pm0.24$ Wb, 而磁链观测误差由0.1 Wb增大到大于0.15 Wb, 误差波动率大于50 %, 观测精度显著降低, 电机的稳态及动态性能将受到影响.
为验证中立型转子磁链观测器的可行性, 本文使用"电力电子与电气传动综合实验台"进行实验验证.实验平台及其结构原理图如图 14所示.实验平台由电机、负载、主回路、PC机、TMS320F2812 (DSP)控制板和保护电路等部分组成, 可以选择按转子磁链定向、空间矢量以及脉宽调制等模式实验, 本文选择按转子磁链定向模式, 具有信息采集功能, 具有绘制电机三相电流、转速和磁链的实时波形的能力, DSP控制器具有编程下载执行能力, 电机使用的是三相鼠笼式感应电机, 电机参数、观测器参数及实验给定参数与仿真时相同.实验结果如图 15和图 16所示.
由图 15 (a)可知, 低速条件下, 电机从启动加速到给定值500 r/min, 上升时间约为0.1 s, 这表明系统具有较好的速度动态响应性能; 图 15 (b)说明在转速阶跃过程中电流变化稳定, 说明了中立型转子磁链观测方法的有效性; 由图 15 (c)可知电机转速从500 r/min调节至650 r/min稳定转速所用时间约为0.04 s, 表明电机具有较好的低速动态响应性能.图 16 (a)是高速条件下的实验调速波形图, 电机转速从750 r/min调节至1500 r/min, 所用时间约为0.14 s, 表明电机具有较好的高速动态响应性能; 由图 16 (b)中可以看出, 高速条件下, 转子磁链观测波形近似单位圆, 观测值收敛于圆内, 表明中立型磁链观测器对转子磁链观测的准确性.实验结果和仿真结果的曲线趋势一致, 证明了中立型转子磁链观测方法的切实可行性.
5. 结论
本文提出了一种基于中立型系统理论的感应电机磁链观测方法, 将中立型延迟系统引入到感应电机磁链观测模型中, 运用线性矩阵不等式理论证明了中立型转子磁链观测器的稳定性.通过仿真分析和实验验证, 得出所提方法有效提高了磁链观测精度, 削弱了电机参数变化对磁链观测精度的影响, 解决系统控制延时对磁链观测的影响问题, 增强了系统观测的鲁棒性, 且该观测方法具有参数自整定, 时效性好, 使用范围广的优点, 证明了所提方法和设计转子磁链观测器的可行性.
-
表 1 角色表
Table 1 Role table
角色类型 角色诞生及身份转变方式 普通用户 注册华谱系统的普通用户, 可进行创建家谱、查看公开数据等. 普通家谱成员 普通用户向某一共建家谱申请成为家谱成员, 只拥有针对该家谱的最基本权限, 例如查看该家谱中的基本信息. 家谱共建者 家谱成员向家谱创建者或核心修谱成员申请成为家谱共建者, 拥有上传数据、对本人上传的数据以及其他用户分享的数据拥有基本修改权限. 若涉及家谱主树结构变化, 需经过审核. 数据合作拥有者 包含部分、全部合作拥有者, 家谱成员向数据拥有者申请数据增加、修改等权限, 成为数据合作拥有者. 数据拥有者 家谱成员保存个人家谱数据, 成为数据拥有者, 拥有数据的全部权限, 且在共建家谱中其权限可转让或共享. 核心修谱成员 家谱成员向家谱创建者申请成为核心修谱成员. 拥有对该家谱数据的大部分权限, 包括查看、编辑、审核本家谱所有家谱人物. 不具备指定核心修谱成员、修改家谱名称等少数信息的权限. 家谱创建者 拥有该家谱的所有权限, 可指定每位家谱成员的权限级别. 
下载: 导出CSV
表 2 数据权限表
Table 2 Data authority table
数据类型公开级 相关描述 完全公开10 如公开家谱、百度百科人物;
所有人可进行查看、编辑完全公开9 如公开家谱、百度百科人物;
所有人可进行查看.8 针 对 共
建 家 谱
而 言部分用户可查看 7 部分用户可修改人物信息 6 部分用户可修改人物关系 5 部分用户可删除人物 4 仅数据录入者可修改 3 私有家谱 2 普通管理员可修改查看 1 超级管理员可修改、查看
下载: 导出CSV
-
[1] 陈宁宁. 家谱研究历史现状. 图书馆杂志, 1998, 17(2): 12−1318Chen Ning-Ning. History and current status of genealogical research. Library Journal, 1998, 17(2): 12−1318 [2] Carolina N, Nils G, Hilary C, and Alexander L. Lineage: Visualizing Multivariate Clinical Data in Genealogy Graphs. IEEE Transactions on Visualization and Computer Graphics, 2019, 25(1): 544−554 [3] Hayden E C. Colossal family tree reveals environment's influence on lifespan. Nature, 2018. [4] 湛庐. 家谱中的文献问题. 北京大学学报(哲学社会科学版), 2007, (1): 150−151Zhan Lu. Literature questions in genealogy. Journal of Peking University (Philosophy and Social Sciences), 2007, (1): 150−151 [5] 欧阳康. 大数据与人文社会科学研究的变革与创新. 光明日报, 2016-11-10(016).Ou Yang-Kang. The reform and innovation of big data and humanities and social science research. Guangming Daily, 2016-11-10(016). [6] 孙建军. 大数据时代人文社会科学如何发展. 光明日报, 2014-07-07(011).Sun Jian-Jun. How to develop humanities and social sciences in the age of big data. Guangming Daily, 2014-07-07(011). [7] Wu X D, Chen H H, Wu G Q, et al. Knowledge engineering with big data. IEEE Intelligent Systems, 2015, 30(5): 46−55 [8] Wu X D, Zhu X Q, Wu G Q, Ding W. Data mining with big data. IEEE Transactions on Knowledge and Data Engineering. 2014, 36(1): 97−107. [9] 吴信东, 何进, 陆汝钤, 郑南宁. 从大数据到大知识: HACE+BigKE. 自动化学报, 2016, 42(7): 965−982Wu Xin-Dong, He Jin, Lu Ru-Qian, Zheng Nan-Ning. From big data to big knowledge: HACE + BigKE. Acta Automatica Sinica, 2016, 42(7): 965−982 [10] Wu M H, Wu X D. On big wisdom. Knowledge and Information Systems, 2018, 58(1): 1−8 [11] 钟华. 企业IT架构转型之道: 阿里巴巴中台战略思想与架构实战. 北京: 机械工业出版社, 2017.Zhong Hua. The Transformation of IT Framework in Enterprises: the Strategic Thinking and Framework of Alibaba. Beijing: China Machine Press, 2017. [12] 付登坡, 江敏, 任寅姿等. 数据中台: 让数据用起来. 北京: 机械工业出版社, 2020.Fu Deng-Po, Jiang Min, Ren Yin-Zi, etc. Data Middle Office: Make Data Valuable. Beijing: China Machine Press, 2020. [13] 陈新宇, 罗家鹰, 邓通, 江威. 中台战略: 中台建设与数字商业. 北京: 机械工业出版社, 2019.Chen Xin-Yu, Luo Jia-Ying, Deng Tong, Jiang Wei. Middle-Platform Strategy: Middle-Platform Construction and Digital Commerce. Beijing: China Machine Press, 2019. [14] 吴信东, 董丙冰, 堵新政, 杨威. 数据治理技术. 软件学报, 2019, 30(9): 2830−2856Wu Xin-Dong, Dong Bing-Bing, Du Xin-Zheng, Yang Wei. Data governance technology. Ruan Jian Xue Bao/Journal of Software, 2019, 30(9): 2830−2856 [15] 吴信东, 嵇圣硙. MapReduce与Spark用于大数据分析之比较. 软件学报, 2018, 29(6): 260−281Wu Xin-Dong, Ji Sheng-Wei. Comparative study on MapReduce and Spark for big data analytics. Ruan Jian Xue Bao/Journal of Software, 2018, 29(6): 260−281 [16] Ji S W, Bu C Y, Li L, W XD. Local graph edge partitioning with a two-Stage heuristic method. In: Proceedings of the 39th IEEE International Conference on Distributed Computing Systems. Dallas, Texas, USA. IEEE, 2019. [17] 吴共庆, 胡骏, 李莉, 徐喆昊, 刘鹏程, 胡学钢, 吴信东. 基于标签路径特征融合的在线Web新闻内容抽取. 软件学报, 2016, 27(3): 714−735Wu Gong-Qing, Hu Jun, Li Li, Xu Zhe-Hao, Liu Peng-Cheng, Hu Xue-Gang, Wu Xin-Dong. Online web news extraction via tag path feature fusion. Ruan Jian Xue Bao/Journal of Software, 2016, 27(3): 714−735 [18] 吴信东, 李娇, 周鹏, 卜晨阳. 碎片化家谱数据的融合技术. 软件学报. http://www.jos.org.cn/1000-9825/6010.htm.Wu Xin-Dong, Li Jiao, Zhou Peng, Bu Chen-Yang. A fusion technique for fragmented genealogy data. Ruan Jian Xue Bao/Journal of Software, 2020 http://www.jos.org.cn/1000-9825/6010.htm. [19] Liu X J, Zhu Y, Ji S W. Web log analysis in genealogy system. In: Proceedings of the 11th IEEE International Conference on Knowledge Graph. Nanjing, China. IEEE, 2020. 期刊类型引用(7)
1. 陈健，万佳泽，林丽，李佐勇. 自适应语义感知网络的盲图像质量评价. 中国图象图形学报. 2023(11): 3400-3414 . 
百度学术2. 程茹秋，余烨，石岱宗，蔡文. 图像与视频质量评价综述. 中国图象图形学报. 2022(05): 1410-1429 . 百度学术3. 刘兴奥，周日贵，郭文宇. 量子线性卷积及其在图像处理中的应用. 自动化学报. 2022(06): 1504-1519 . 本站查看4. 汪杰，陈曼龙，李奎，杨帆，燕立志. 机器视觉螺纹图像评价方法. 应用光学. 2022(05): 904-912 . 百度学术5. 温静，白鑫. 自适应融合局部和全局特征的图像质量评价. 计算机技术与发展. 2022(11): 50-57 . 百度学术6. 邓杰航，袁仲鸣，林好润，顾国生. 协同超像素和视觉显著性的图像质量评价. 广东工业大学学报. 2021(05): 33-39 . 百度学术7. 邓杰航，袁仲鸣，刘栋濠，顾国生. 基于视觉显著性的自适应图像质量评价方法. 软件导刊. 2021(11): 191-196 . 百度学术其他类型引用(15)
-
下载:
计量
- 文章访问数: 4123
- HTML全文浏览量: 647
- PDF下载量: 643
- 被引次数: 22
