Combining Global and Local Variation for Image Quality Assessment
-
摘要: 图像所包含的信息是通过灰度值在空域的变化呈现的. 梯度是度量变化的基本工具, 这使得梯度成为了目前大多数图像质量评价算法的重要组成部分. 但是梯度只能度量局部变化, 而当人类视觉系统(Human visual system, HVS)感知一幅图像时, 既能感知到局部变化, 也能感知到全局变化. 基于HVS的这一特性, 本文提出了一种结合全局与局部变化的图像质量评价算法(Global and local variation similarity, GLV-SIM). 该算法利用Grünwald-Letnikov分数阶导数来度量图像的全局变化, 利用梯度模来度量图像的局部变化. 然后结合二者计算参考图像和退化图像之间的相似度谱(Similarity map), 进而得到图像的客观评分. 在TID2013、TID2008、CSIQ与LIVE四个数据库上的仿真实验表明, 较之单一度量局部变化的方法, 本文算法能更准确地模拟HVS对图像质量的感知过程, 给出的客观评分与主观评分具有较好的一致性.
-
关键词:
- 图像质量评价 /
- 全局与局部 /
- Grünwald-Letnikov分数阶导数 /
- 梯度模
Abstract: The information contained in an image is presented by the changes of gray value in spatial domain. Gradient is the basic tool to measure changes, which makes gradient become an important ingredient of most image quality assessment algorithms. However, gradient can only measure local changes, while when human visual system (HVS) perceives an image, it can perceive both local and global changes. Based on this characteristic of HVS, this paper proposes an image quality assessment algorithm by combining global and local variations similarity (GLV-SIM). The algorithm uses Grünwald-Letnikov fractional derivative to measure the global changes and uses gradient magnitude to measure the local changes of the image. Synthesizing the two aspect changes, similarity map between reference image and distorted image is calculated, and then objective score of the image is obtained. Simulation experiments on four databases TID2013, TID2008, CSIQ and LIVE show that, comparing with the algorithm only considering local changes, the proposed algorithm can more accurately simulate the perception process of HVS on image quality and can obtain better consistency between objective scores and subjective scores. -
感应电机交流传动系统广泛应用于各工业领域, 其中新能源汽车等场合普遍采用转子磁链定向的控制方式, 获得精确的转子磁链是实现交流传动系统高性能控制的关键[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. 结论
本文提出了一种基于中立型系统理论的感应电机磁链观测方法, 将中立型延迟系统引入到感应电机磁链观测模型中, 运用线性矩阵不等式理论证明了中立型转子磁链观测器的稳定性.通过仿真分析和实验验证, 得出所提方法有效提高了磁链观测精度, 削弱了电机参数变化对磁链观测精度的影响, 解决系统控制延时对磁链观测的影响问题, 增强了系统观测的鲁棒性, 且该观测方法具有参数自整定, 时效性好, 使用范围广的优点, 证明了所提方法和设计转子磁链观测器的可行性.
-
图 3 参考图像(a)及其不同类型退化图像(b)~(f) (右下角为矩形区域局部放大图)
Fig. 3 Reference image (a) and different types of distorted images (b)~(f)
(The lower right corner is a enlarged view of the rectangular region)
评价方法 图3(b) 图3(c) 图3(d) 图3(e) 图3(f) MOS 5.0000 3.8387 4.1875 4.7667 6.2903 PSNR 30.5304 30.5784 26.1303 27.4808 27.3498 VSNR 29.7301 21.1480 20.5342 30.7072 20.2681 IFC 4.7389 3.4351 4.9319 2.9956 11.3746 SSIM 0.9250 0.8461 0.9459 0.9475 0.9568 MS-SSIM 0.9606 0.9159 0.9738 0.9727 0.9821 IW-SSIM 0.9684 0.9075 0.9645 0.9704 0.9661 GSIM 0.9958 0.9888 0.9953 0.9966 0.9979 FSIM 0.9831 0.9462 0.9538 0.9699 0.9707 GLV-SIM 0.9959 0.9845 0.9927 0.9957 0.9961 
下载: 导出CSV
下载: 导出CSV
表 3 不同IQA算法在TID2013和TID2008数据库的实验结果比较
Table 3 Comparison the performance results of different IQA algorithms on TID2013 and TID2008 databases
数据库 性能指标 PSNR VSNR IFC SSIM MS-SSIM IW-SSIM GSIM FSIM GLV-SIM TID 2013 SROCC 0.6396 0.6812 0.5389 0.7417 0.7859 0.7779 0.7946 0.8015 0.8068 KROCC 0.4698 0.5084 0.3939 0.5588 0.6047 0.5977 0.6255 0.6289 0.6381 PLCC 0.7017 0.7402 0.5538 0.7895 0.8329 0.8319 0.8464 0.8589 0.8580 RMSE 0.8832 0.8392 1.0322 0.7608 0.6861 0.6880 0.6603 0.6349 0.6368 TID 2008 SROCC 0.5531 0.7046 0.5675 0.7749 0.8542 0.8559 0.8504 0.8805 0.8814 KROCC 0.4027 0.5340 0.4236 0.5768 0.6568 0.6636 0.6596 0.6946 0.6956 PLCC 0.5734 0.6820 0.7340 0.7732 0.8451 0.8579 0.8422 0.8738 0.8648 RMSE 1.0994 0.9815 0.9113 0.8511 0.7173 0.6895 0.7235 0.6525 0.6739
下载: 导出CSV
表 4 不同IQA算法在CSIQ和LIVE数据库的实验结果比较
Table 4 Comparison the performance results of different IQA algorithms on CSIQ and LIVE databases
数据库 性能指标 PSNR VSNR IFC SSIM MS-SSIM IW-SSIM GSIM FSIM GLV-SIM CSIQ SROCC 0.8058 0.8106 0.7671 0.8756 0.9133 0.9213 0.9108 0.9242 0.9264 KROCC 0.6084 0.6247 0.5897 0.6907 0.7393 0.7529 0.7374 0.7567 0.7605 PLCC 0.8000 0.8002 0.8384 0.8613 0.8991 0.9144 0.8964 0.9120 0.9082 RMSE 0.1575 0.1575 0.1431 0.1334 0.1149 0.1063 0.1164 0.1007 0.1099 LIVE SROCC 0.8756 0.9274 0.9259 0.9479 0.9513 0.9567 0.9561 0.9634 0.9521 KROCC 0.6865 0.7616 0.7579 0.7963 0.8045 0.8175 0.8150 0.8337 0.8179 PLCC 0.8723 0.9231 0.9268 0.9449 0.9489 0.9522 0.9512 0.9597 0.9368 RMSE 13.359 10.505 10.264 8.9445 8.6188 8.3473 8.4327 7.6780 8.0864
下载: 导出CSV
表 5 不同IQA算法在TID2008数据库单一失真评价性能(SROCC)比较
Table 5 Comparison SROCC for individual distortion of different IQA algorithms on TID2008 database
数据库 失真类型 PSNR VSNR IFC SSIM MS-SSIM IW-SSIM GSIM FSIM GLV-SIM TID 2008 AWN 0.9073 0.7728 0.5806 0.8107 0.8094 0.7869 0.8573 0.8566 0.9125 ANMC 0.8994 0.7793 0.5460 0.8029 0.8064 0.7920 0.8095 0.8527 0.8979 SCN 0.9175 0.7665 0.5958 0.8144 0.8195 0.7714 0.8902 0.8483 0.9167 MN 0.8520 0.7295 0.6732 0.7795 0.8156 0.8087 0.7403 0.8021 0.8087 HFN 0.9273 0.8811 0.7318 0.8729 0.8685 0.8662 0.8932 0.9093 0.9175 IMN 0.8725 0.6470 0.5345 0.6732 0.6868 0.6465 0.7721 0.7452 0.7864 QN 0.8702 0.8271 0.5857 0.8531 0.8537 0.8177 0.8750 0.8564 0.8865 GB 0.8704 0.9330 0.8559 0.9544 0.9607 0.9636 0.9585 0.9472 0.9587 DEN 0.9422 0.9286 0.7973 0.9530 0.9571 0.9473 0.9723 0.9603 0.9666 JPEG 0.8723 0.9174 0.8180 0.9252 0.9348 0.9184 0.9391 0.9279 0.9534 JP2K 0.8131 0.9515 0.9437 0.9625 0.9736 0.9738 0.9755 0.9773 0.9751 JGTE 0.7525 0.8056 0.7909 0.8678 0.8736 0.8588 0.8832 0.8708 0.8793 J2TE 0.8312 0.7909 0.7301 0.8577 0.8522 0.8203 0.8925 0.8544 0.9021 NEPN 0.5812 0.5716 0.8418 0.7107 0.7336 0.7724 0.7372 0.7491 0.7271 BLOCK 0.6194 0.1926 0.6770 0.8462 0.7617 0.7623 0.8865 0.8492 0.8960 MS 0.6966 0.3715 0.4250 0.7231 0.7374 0.7067 0.7174 0.6698 0.6994 CTC 0.5867 0.4239 0.2713 0.5246 0.6398 0.6301 0.6736 0.6481 0.6689
下载: 导出CSV
-
[1] Mohammadi P, Ebrahimi-Moghadam A, Shirani S. Subjective and objective quality assessment of image: a survey. Majlesi Journal of Electrical Engineering, 2015, 9(1): 419−423 [2] Yang X C, Sun Q S, Wang T S. Image quality assessment improvement via local gray-scale fluctuation measurement. Multimedia Tools and Applications, 2018, 77(18): 24185−24202 doi: 10.1007/s11042-018-5740-z [3] 高敏娟, 党宏社, 魏立力, 张选德. 基于非局部梯度的图像质量评价算法. 电子与信息学报, 2019, 41(5): 1122−1129 doi: 10.11999/JEIT180597Gao Min-Juan, Dang Hong-She, Wei Li-Li, Zhang Xuan-De. Image quality assessment algorithm based on non-local gradient. Journal of Electronics and Information Technology, 2019, 41(5): 1122−1129 doi: 10.11999/JEIT180597 [4] 许丽娜, 肖奇, 何鲁晓. 考虑人类视觉特征的融合图像评价方法. 武汉大学学报(信息科学版), 2019, 44(4): 546−554Xu Li-Na, Xiao Qi, He Lu-Xiao. Fused image quality asssessment based on human visual characterisyics. Geomatics and Information Science of Wuhan University, 2019, 44(4): 546−554 [5] 曹清洁, 史再峰, 张嘉平, 李杭原, 高静, 姚素英. 分区域多标准的全参考图像质量评价算法. 天津大学学报(自然科学与工程技术版), 2019, 52(6): 625−630Cao Qing-Jie, Shi Zai-Feng, Zhang Jia-Ping, Li Hang-Yuan, Gao Jing, Yao Su-Ying. Sub-regional and multiple criteria full-reference image quality assessment. Journal of Tianjin University (Science and Technology), 2019, 52(6): 625−630 [6] Wang Z, Bovik A C. Mean squared error: love it or leave it? A new look at signal fidelity measures. IEEE Signal Processing Magazine, 2009, 26(1): 98−117 doi: 10.1109/MSP.2008.930649 [7] Huynh-Thu Q, Ghanbari M. Scope of validity of PSNR in image/video quality assessment. Electronics Letters, 2008, 44(13): 800−801 doi: 10.1049/el:20080522 [8] Wang Z, Bovik A C, Sheikh H R, Simoncelli E P. Image quality assessment: from error visibility to structural similarity. IEEE Transactions on image processing, 2004, 13(4): 600−612 doi: 10.1109/TIP.2003.819861 [9] Wang Z, Simoncelli E P, Bovik A C. Multiscale structural similarity for image quality assessment. In: Proceedings of the 37th IEEE Asilomar Conference on Signals, Systems, and Computers, Pacific Grove, USA, 2003. 1398−1402 [10] Wang Z, Li Q. Information content weighting for perceptual image quality assessment. IEEE Transactions on Image Processing, 2011, 20(5): 1185−1198 doi: 10.1109/TIP.2010.2092435 [11] Li C F, Bovik A C. Three-component weighted structural similarity index. In: Proceedings of the 2009 SPIE Conference on Image Quality and System Performance, San Jose, California, USA, 2009. 7242: 72420Q−72420Q-9 [12] 刘大瑾, 叶建兵, 刘家骏. SSIM框架下基于SVD的灰度图像质量评价算法研究. 南京师大学报(自然科学版), 2017, 40(1): 73−78Liu Da-Jin, Ye Jian-Bing, Liu Jia-Jun. SVD-based gray-scale image quality assessmentalgorithms in the SSIM perspective. Journal of Nanjing Normal University (Natural Science Edition), 2017, 40(1): 73−78 [13] 闫乐乐, 李辉, 邱聚能, 梁平. 基于区域对比度和SSIM的图像质量评价方法. 应用光学, 2015, 36(1): 58−63 doi: 10.5768/JAO201536.0102002Yan Le-Le, Li Hui, Qiu Ju-Neng, Liang Ping. Image qualty assessment method based on regional contrast and structural similarity. Journal of Applied Optics, 2015, 36(1): 58−63 doi: 10.5768/JAO201536.0102002 [14] Jin X, Jiang G Y, Chen F, Yu M, Shao F, Peng Z J, et al. Adaptive image quality assessment method based on structural similarity. Journal of Optoelectronics. Laser, 2014, 25(2): 378−385 [15] 田浩南, 李素梅. 基于边缘的SSIM图像质量客观评价方法. 光子学报, 2013, 42(1): 110−114 doi: 10.3788/gzxb20134201.0110Tian Hao-Nan, Li Su-Mei. Objective evaluation method for image quality based on edge structural similarity. Acta Photonica Sinca, 2013, 42(1): 110−114 doi: 10.3788/gzxb20134201.0110 [16] Zhang L, Zhang L, Mou X Q, Zhang D. FSIM: a feature similarity index for image quality assessment. IEEE Transactions on Image Processing, 2011, 20(8): 2378−2386 doi: 10.1109/TIP.2011.2109730 [17] Liu A M, Lin W S, Narwaria M. Image quality assessment based on gradient similarity. IEEE Transactions on Image Processing, 2012, 21(4): 1500−1512 doi: 10.1109/TIP.2011.2175935 [18] Zhang X D, Feng X C, Wang W W, Xue W F. Edge strength similarity for image quality assessment. IEEE Signal Processing Letters, 2013, 20(4): 319−322 doi: 10.1109/LSP.2013.2244081 [19] Sun W, Liao Q M, Xue J H, Zhou F. SPSIM: a superpixel-based similarity index for full-reference image quality assessment. IEEE Transactions on Image Processing, 2018, 27(9): 4232−4244 doi: 10.1109/TIP.2018.2837341 [20] Ding L, Huang H, Zang Y. Image quality assessment using directional anisotropy structure measurement. IEEE Transactions on Image Processing, 2017, 26(4): 1799−1809 doi: 10.1109/TIP.2017.2665972 [21] Xue W F, Zhang L, Mou X Q, Bovik A C. Gradient magnitude similarity deviation: a highly efficient perceptual image quality index. IEEE Transactions on Image Processing, 2014, 23(2): 684−695 doi: 10.1109/TIP.2013.2293423 [22] Pei S C, Chen L H. Image quality assessment using human visual DOG model fused with random forest. IEEE Transactions on Image Processing, 2015, 24(11): 3282−3292 doi: 10.1109/TIP.2015.2440172 [23] Reisenhofer R, Bosse S, Kutyniok G, Wiegand T. A haar wavelet-based perceptual similarity index for image quality assessment. Signal Processing: Image Communication, 2018, 61: 33−43 doi: 10.1016/j.image.2017.11.001 [24] 薛定宇. 分数阶微积分学与分数阶控制. 北京: 科学出版社, 2018. 31−72Xue Ding-Yu. Fractional Calculus and Fractional-order Control. Beijing: Publishing House of Science, 2018. 31−72 [25] 张桂梅, 孙晓旭, 刘建新, 储珺. 基于分数阶微分的TV-L1光流模型的图像配准方法研究. 自动化学报, 2017, 43(12): 2213−2224Zhang Gui-Mei, Sun Xiao-Xu, Liu Jian-Xin, Chun Jun. Research on TV-L1 optical flow model for image registration based on fractional-order differentiation. Acta Automatica Sinica, 2017, 43(12): 2213−2224 [26] 陈云, 郭宝裕, 马祥园. 基于分数阶微积分正则化的图像处理. 计算数学, 2017, (39): 393−406Chen Yun, Guo Bao-Yu, Ma Xiang-Yuan. Image processing based on regularization with fractional clculus. Mathematica Numerica Sinica, 2017, (39): 393−406 [27] Ponomarenko N, Jin L, Ieremeiev O, Lukin V, Egiazarian K, Astola J, et al. Image database TID2013: peculiarities, results and perspectives. Image Communication, 2015, 30(C): 57−77 [28] Ponomarenko N, Lukin V, Zelensky A, Egiazarian K, Astola J, Carli M, et al. TID2008: a database for evaluation of full-reference visual quality assessment metrics [Online], available: http://www.ponomarenko.info/papers/mre, November 1, 2016 [29] Larson E C, Chandler D. Categorical subjective image quality (CSIQ) database [Online], available: http://vision.okstate.edu/csiq, November 1, 2016 [30] Sheikh H R, Wang Z, Bovik A C. Image and video quality assessment research at LIVE. [Online], available: http://live.ece.utexas.edu/rese-arch/quality/, November 1, 2016 [31] Chandler D M, Hemami S S. VSNR: a wavelet-based visual signal-to-noise ratio for natural images. IEEE Transactions on Image Processing, 2007, 16(9): 2284−2298 doi: 10.1109/TIP.2007.901820 [32] Sheikh H R, Bovik A C, Veciana G D. An information fidelity criterion for image quality assessment using natural scene statistics. IEEE Transactions on Image Processing, 2005, 14(12): 2117−2128 doi: 10.1109/TIP.2005.859389 期刊类型引用(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)
-
下载:
计量
- 文章访问数: 4649
- HTML全文浏览量: 2018
- PDF下载量: 316
- 被引次数: 22
