第 34卷 第 8期 控 制 与 决 策 Vol.34 No.8
2019年 8月 Control and Decision Aug. 2019
文章编号: 1001-0920(2019)08-1702-07 DOI: 10.13195/j.kzyjc.2018.0019
广义逆张量Padé逼近的连分式递推算法
顾传青
1†
, 黄逸铮
1
, 陈之兵
2
(1. 上海大学 理学院,上海 200444;2. 深圳大学 数学与统计学院,广东 深圳 518052)
摘 要: 张量指数函数已经广泛应用于控制论、图像处理和各个工程领域. 鉴于此, 在矩阵广义逆的基础上, 首次
在张量内积空间上定义一种有效的张量广义逆, 从而构造张量 Padé逼近的一种连分式算法. 利用张量 t-积成功计
算张量的幂, 由此递推地给出张量指数函数的幂级数展开式. 在前面两个工作的基础上, 利用设计的连分式算法
逼近张量指数函数, 其特点在于, 该算法可以编程实现递推计算, 而且在计算过程中不必计算张量的乘积, 也不必
计算张量的逆. 给出的两个张量指数函数的数值实验表明, 将连分式算法与目前通常使用的截断法进行比较, 在
不降低逼近阶的条件下, 所提出算法是有效的. 如果张量的维数较大, 基于张量广义逆的连分式算法仍然具有一
定优势.
关键词: 张量;张量t-积;张量广义逆;张量Padé逼近;张量指数函数;张量连分式算法
中图分类号: O231.1 文献标志码: A
A continued fractional recurrence algorithm for generalized inverse tensor
Padé approximation
GU Chuan-qing
1†
, HUANG Yi-zheng
1
, CHEN Zhi-bing
2
(1. College of Sciences,Shanghai University,Shanghai 200444,China;2. School of Mathematics and Statistics,
Shenzhen University
,
Shenzhen 518052
,
China)
Abstract: The tensor exponential function has been widely used in cybernetics, image processing and various engineering
fields. Based on the generalized matrix inverse, an effective tensor generalized inverse is defined for the first time on the
scalar inner product space, thus constructing a continued fractional algorithm for the tensor Pad
´
e approximation. On the
other hand, we successfully use the tensor t- product to calculate the power of the tensor, and recursively giving the power
series expansion of the tensor exponential function. Based on the previous two work, the continuous fractional algorithm
designed in this paper is used to approximate the tensor exponential function. Its characteristic is that the algorithm can
be programmed to implement recursive calculations, and in the calculation process, it is not necessary to calculate the
product of the tensor and to calculate the inverse of the tensor. The numerical experiments of the two tensor exponential
functions given in this paper show that comparing the continuous fractional algorithm with the commonly used truncation
method, the proposed algorithm is effective without reducing the approximation order. If the dimension of the tensor
is relatively large, a continuous fractional algorithm based on the generalized inverse of tensors will also have certain
advantages.
Keywords: tensor;tensor t-product;tensor generalized inverse;tensor Padé approximation;tensor exponential functions;
continued fraction expression algorithm
0 引 言
张 量 动 力 系 统 已 广 泛 应 用 于 Volterra 系统识
别
[1]
、张量乘积 (TP) 模型变换
[2]
、人体动作识别
[3]
和
塑性模型
[4]
等各个领域. 由于在张量常微分方程的
求解过程中, 必须用到张量指数函数, 其计算问题已
经成为一个重要的研究领域.
考虑如下常微分方程
[4]
:
˙
Y (t) = AY (t),
Y (t
0
) = Y
0
,
(1)
其中 A 和 Y
0
是给定的张量, 一般是非对称的, 则关于
系统(1)的张量指数函数exp(·)有如下唯一解:
Y (t) = exp[(t − t
0
)A]Y
0
. (2)
对于一般的张量 A, 其指数函数通常表示为幂级
收稿日期: 2018-01-03;修回日期: 2018-03-20.
基金项目: 国家自然科学基金项目 (11371243);上海市重点学科建设项目 (S30104).
†
通讯作者. E-mail: cqgu@shu.edu.cn.