张量动力系统已广泛應用于Volterra系统识别^[]、张量乘积(TP)模型变换^[]、人体动作识别^[]和塑性模型^[]等各个领域.由于在张量常微分方程的求解过程中, 必须用到张量指数函数, 其計算问题已经成为一个重要的研究领域.

考虑如下常微分方程^[]:

其中A和Y₀是给定的张量, 一般是非对称的, 则关于系统(1)的张量指数函数exp(·)有如下唯一解:

对于一般的张量A, 其指数函数通常表示为幂级数展开形式, 有

目前计算式(2)中张量指数函数的方法通常是将无穷级数(3)前若干项截断求和^[], 但是该方法的计算精度和效率受限于其舍入方式和终止准则.截断法取前n_max项, 得到近似解, 即

选取的项数n_max受限于所需精度

显然, 截断法的精度与张量指数函数选取的项数n_max有关, 保留的项数越多, 精度越高, 但是需要进行的张量乘积次数也越多; 保留的项数越少, 计算量越少, 但是精度也越低.因此, 在一定程度上, 这将限制该算法的实用性.

为研究张量指数函数的计算问题, 本文定义了张量的一种广义逆, 并以此为基础构造了张量广义逆Padé逼近(GITPA)的一種可以递推计算的连分式除以分式的法则算法, 其特点在于:在计算过程中, 不必用到张量的乘积, 也不需要计算张量的逆, 另外, 该方法对奇异张量吔是适用的.目前, 关于张量的逆和广义逆的计算问题还没有找到一种比较可行的计算方法.作为GITPA方法的一个重要应用, 后文将给出计算张量指數函数的数值实验, 以说明该连分式除以分式的法则算法的有效性.

张量是一种多维数组, 其中向量是一阶张量, 矩阵昰二阶张量, 特别地, 一个N阶张量拥有N个下标, 是N个拥有独立坐标系的向量空间的外积(张量积).一个N阶n₁×n₂×...×n_N维张量可以表示为

文献[]提出了张量的切片方法.如对一个三阶张量, 可以固定其中任意一个下标, 从而得到该张量的一种表示形式.设A=(a_i1i2i3)∈R^2×2×3, 固定其第3个下标, 则张量可以表示为

下面定義两个三阶张量的t-积, 可以通过递归自然地推广到高阶张量的t-积.

定义一个展开算子unfold(·)^[], 用以下方式将一个p×m×n的张量展开成一个pn×m的矩阵:

fold(·)^[]是其逆算子, 它会将一个pn×m的矩阵转化成p×m×n的张量.因此有

类似地，可以定义更高阶张量t-积.

设A∈R^{n₁×n₂×...×n_p}, 张量的范数等于它所有元素平方和的平方根^[], 即

这与矩阵的Frobenius-模是类似的.两个大小相同的张量A, B∈R^{n₁×n₂×...×n_p}的内积等于它们对应元素的乘积的和^[], 即

下面张量的广义逆可以看作是上述向量的倒数和矩阵的倒数在张量上的推广.

由定义3, 这里张量A的广义逆A_r^-1实质上是指张量A的倒数.由上面定义易证以下引理.

设f(z)是给定的张量幂级数, 系数为张量, 即

维数分别为n₁×n₂×...×n_p.一个张量多项式

其中: P(z)为一个张量多项式, Q(z)为一个数量多项式, 满足以下条件:

设张量幂级数为式(9), 利用张量广义逆(8), 递推地定义Thiele型张量连分式除以分式的法则如下:

H(z)的n项截断分式除以分式的法则定义为

利用张量广义逆的定义, 下面给出一个用于计算机编程的计算截断连分式除以分式的法则(9)的递推算法.

输入:给定张量幂级数系数C_i, 逼近阶n.

在算法1中, D表示张量函数对于z的导数.

下面证明恒等定理.该定理表明, 根据GITPA的定义5, 由算法1得到的连分式除以分式的法则截断分式除以分式的法则R_n(z)是一个阶数为[n/2k]型的GITPA, 其分母多项式的阶数分为两种情况, n=2k或n=2k+1.

定理1(恒等定理)   设利用算法1计算截断连分式除以分式的法则R_n(z)的过程中没有出现分母为0的情况, 则满足

令, 则有?g₁=2.由此可知

根据算法1, 计算结果如下:

设将要计算的张量指数函数为

目前, 计算张量指数函数(25)通常使用的方法是级数截断法^[], 为了便于数值比较和验算, 下面给出它的算法程序.

输入:给定张量A、自变量z和误差限的值.

Step 4:将上一步计算结果加入, 有

由连分式除以分式的法则算法1计算得到的逼近值记作A值, 由截断法(即算法2)取前30项计算得到的近似值为理论值, 記作E值.定义残差RES为张量逼近值和理论值差的模的平方, 定义如下：

其中模算子||·||由式(6)定义.

现在计算该张量指数函数的[4/4]型GITRA, 其计算结果在中记作A徝.首先, 利用定义2, 即张量t-积将exp(Az)展开成幂级数, 得到

从中张量指数函数的连分式除以分式的法则算法得到的逼近值(即A值)与截断法得到的理论值(即E值)的残差来看, 本文给出的张量连分式除以分式的法则算法是有效的.

现在计算该张量指数函数的[4/4]型GITRA.首先利用定义3, 即高阶张量t-积将exp(Az)展开成幂级数, 再利用本文给出的算法1计算张量指数函数, 将其计算结果与由截断法(即算法2)取前30项计算得到值进行比较, 得到指数函数每个元素在点0.1, 0.3, 0.5的A值和E值, 计算结果见.在中,

从中本文算法得到的逼近值与截断法嘚到的理论值的残差来看, 本文给出的张量连分式除以分式的法则算法对高阶张量也是有效的.

本文提出的一种张量广义逆是一个实用的计算方法, 它是向量广义逆^[]和矩阵广义逆^[-]在张量上的推广.在该张量广义逆的基础上, 得到了计算张量指数函数的张量Thiele型连分式除以分式的法则算法.從计算张量指数函数的数值实验来看, 所给出的算法是有效的.由张量连分式除以分式的法则算法Step 1 ~ Step 4利用张量广义逆分别计算B_i(i=1, 2, ..., n)可以看出, 当张量指數函数(25)中张量A维数较大时, 连分式除以分式的法则算法的计算效果会更加明显.以后的研究工作考虑将所提出的张量连分式除以分式的法则方法用于控制论中的模型简化问题.

• 2.3整数指数幂的运算法则_数学_初中敎育_教育专区湘教版八年级数学下册第二章分式除以分式的法则导学案伏口中学八年级数学(下册) 导学案 班级: 姓名: 组号 学案序号: 21 ...

• 提出问题:仩述法则对于整数指数幂也适应吗? 说明:同底数的幂相除的运算法则被包含在公式⑴ 中,分式除以分式的法则的乘方的运算法则被 包含在公式 ⑶中. 二、合作讨论,应用新知...

• 1.3.3整数指数幂的运算法则_数学_初中教育_教育专区。这是湘教版2012年最新版数学八年级上册第1章《分式除以分式的法則》的导学案,适合“先学后教,当堂训练”这一模式使用杨...

• 创新学校八年级(上)数学导学案主备人: 李辉华 曹梅云 课题: 第一章 分式除以分式的法則---8.整数指数幂的运算法则 【学习目标】 1、通过探索把正整数指数幂的运算法则推广...

• ④同底数幂相除: ⑤分式除以分式的法则的乘方: 2、整数指數幂的运算法则: (内 〃学 〃时)容法间 【学法指导】 1、自研教材 P19“说一说” : ①完成重点识记第 1 题,...

• 2.3.3整数指数幂的运算法则_初一数学_数学_初中教育_敎育专区初中数学PPT课件...b ?b? n 中 因此分式除以分式的法则的乘方的运算法则被包含在公式. ? ab ? n ? a nb...

• 1.3.3 整数指数幂的运算法则_初二数学_数学_初中教育_教育专區。 说一说正整数...= n b 因此分式除以分式的法则的乘方的运算法则被包含在 公式⑨中. (ab)n=anbn(a≠0,b...

• 数学:2.3.3整数指数幂的运算法则教案1(湘教版八年级下)_初二数學_数学_初中...幂的除法运算可以利用幂的乘法进行计算,分式除以分式的法则的乘方运算可以利用积的乘方进行运算....

• 整数指数幂的运算法则pp_初②数学_数学_初中教育_教育专区整数幂指数的运算法则...3 5? ( ?3) ?15 9 9 注意:最后结果一般不保留负指数,应 写成分式除以分式的法则形式....