二元多项式矩阵Smith型等价问题

李冬梅 ,  梁芮

武汉大学学报(理学版) ›› 2020, Vol. 66 ›› Issue (3) : 224 -232.

PDF (491KB)
武汉大学学报(理学版) ›› 2020, Vol. 66 ›› Issue (3) : 224 -232. DOI: 10.14188/j.1671-8836.2020.0031
数学

二元多项式矩阵Smith型等价问题

作者信息 +

On the Equivalence to the Smith Form of Bivariate Polynomial Matrices

Author information +
文章历史 +
PDF (502K)

摘要

二维系统等价是多维系统研究中的重要问题,它常被转化为二元多项式矩阵的等价问题来研究。二元多项式矩阵与其Smith型等价是矩阵等价研究中的重要问题。本文主要研究几类二元多项式矩阵与其Smith型等价问题,给出二元多项式矩阵约化到其Smith型的一些结果和判别条件。这些条件可以通过计算矩阵既约子式的约化Gröbner基进行检验。

Abstract

The equivalence of two-dimensional systems is an important content of multidimensional systems, which is usually represented by the bivariate polynomial matrices. Reducing a matrix to its Smith form is a very important study in the equivalence of matrices. In this paper we mainly investigate the reduction of several kinds of bivariate polynomial matrices to their Smith forms. Some new results and criteria are presented. These criteria can be verified easily by computing the Gröbner basis of the associated ideals.

关键词

二维系统 / 二元多项式矩阵 / 矩阵等价 / Smith型

Key words

two-dimensional system / bivariate polynomial matrix / equivalence of matrices / the Smith form

引用本文

引用格式 ▾
李冬梅,梁芮. 二元多项式矩阵Smith型等价问题[J]. 武汉大学学报(理学版), 2020, 66(3): 224-232 DOI:10.14188/j.1671-8836.2020.0031

登录浏览全文

4963

注册一个新账户 忘记密码

0  引 言

多维系统在数字滤波器、图像、信号处理等领域有广泛应用,多维系统等价是多维系统的一个重要内容,在实际应用中,常常需要将一个给定的复杂系统转化为形式和结构更简单的等价系统[1,2,3,4,5,6,7,8]。由系统理论中的多项式表示方法,人们常将多维系统等价问题转化为多项式矩阵等价问题来研究。矩阵的Smith型在矩阵等价中占非常重要的地位,将一个多项式矩阵化简为其Smith型,可以降低对应多维系统的未知量数或方程数。1970年,Rosenbrock[5]利用主理想整环K[s]的算法性质,得到K[s]上任意一元多项式矩阵都与其Smith型等价,完全解决了一维系统的等价问题。但多元多项式环不是欧几里德环,其环上多元多项式矩阵的Smith型等价问题等价于两个有限生成投射模的同构问题,这类同构问题在代数上几乎是不可能被完全解决的,因此多元多项式矩阵等价问题有待继续深入研究。

科学家们对多元多项式矩阵与其Smith型等价的研究起源于1978年Frost和Storey[9]对实数域上的二元多项式矩阵的研究,他们给出了一个二元多项式矩阵与其Smith型等价的充分必要条件。但是1983年Lee和Zak[10]通过举例证明了上述条件只是必要的不是充分的,并给出了一类形如z2I-A(z1)的二元多项式矩阵与其Smith型等价的一个充分必要条件是存在一个循环向量,但在实际中很难找到程序判别这种循环向量。1986年,Frost和Boudellioua[11]证明了二元行满秩多项式方阵T与其Smith型等价当且仅当存在一个幺模列向量U使得TU存在右逆,但这种向量的寻找同样不容易且TU的右逆存在性不易判别。所以,2006年Lin等[12]用构造性的方法对一类特殊的多项式矩阵进行讨论,这类矩阵M的行列式为detM=z1-f(z2,,zn),并得到这类矩阵与其Smith型In-10n-1,101,n-1detM等价。本文作者在文献[13]拓展了Lin等的结论,证明了如果多项式矩阵M所有的(l-1)×(l-1)阶子式生成单位理想且detM=[z1-f(z2,,zn)]q,则M也与其Smith型In-10n-1,101,n-1detM等价,并给出一种找到等价变换中幺模矩阵的构造性方法。

二元多项式矩阵的等价在多元多项式矩阵等价问题中占重要地位,本文在文献[12]和[13]的基础上,结合二元多项式环的性质以及Gröbner基理论,考虑二元多项式环上几类二元多项式矩阵的等价问题,希望找到一些行列式为其他特殊形式的矩阵,这些矩阵也与其Smith型等价。

1  预备知识

R[z,w]表示实数域R上含两个变量的多项式环;Rl×l[z,w]表示多项式环R[z,w]上的所有l×l级矩阵的集合;Gll(R[z,w])表示在多项式环R[z,w]上的所有l×l级可逆矩阵的集合;0r,t表示r×t级零矩阵,Ir表示r×r级单位矩阵;令F(z,w)Rl×l[z,w]c1(z,w),,ct(z,w)表示矩阵F(z,w)的所有(l-1)×(l-1)级子式,t=l2

定义1F(z,w)Rl×l[z,w],则F(z,w)的Smith型定义为S=(diag{Φi}),其中

Φi=αiαi-1,1ir0,r<il

rF(z,w)的秩,α01,αiF(z,w)i×i级子式的最大公因子且Φi满足整除性Φ1|Φ2|...|Φr

定义2 称矩阵F(z,w)Rl×l[z,w]具有性质DM(determinant minor),如果detF(z,w)F(z,w)的所有(l-1)×(l-1)级子式没有共同零点。

定义3F(z,w),Q(z,w)Rl×l[z,w],如果存在U(z,w),V(z,w)Gll(R[z,w]),使得F(z,w)=U(z,w)Q(z,w)V(z,w),则称F(z,w)Q(z,w)等价。

定义4F(z,w)Rl×m[z,w](lm)是行满秩的,称F(z,w)是零互素(ZLP,zero left prime),如果F(z,w)的所有(l-1)×(l-1)级子式没有共同零点。

引理1[13]F(z,w)Rl×l[z,w]d0(z,w)=z-f(w)detF(z,w)的一个因子,则d0(z,w)c1(z,w),,cr(z,w)没有共同零点(在复数域Cl×l[z,w]),当且仅当wCn-1,都有rank F(f(w),w)=l-1

推论1F(z,w)Rl×l[z,w]d0(z,w)=w-f(z)detF(z,w)的一个因子,那么d0(z,w),c1(z,w),,cr(z,w)没有共同零点(在复数域Cl×l[z,w]),当且仅当zCn-1,都有rank F(z,f(z))=l-1

引理2[13]F(z,w),F1(z,w),F2(z,w)Rl×l[z,w],F(z,w)=F1(z,w)F2(z,w)。若detF(z,w)F(z,w)所有的(l-1)×(l-1)级子式没有共同零点,则detFi(z,w)Fi(z,w)所有的(l-1)×(l-1)级子式没有共同零点,其中i=1,2

推论2F(z,w)Rl×l[z,w],U(z,w),V(z,w)Gll(R[z,w]),那么F(z,w)具有性质DM当且仅当U(z,w)F(z,w)V(z,w)具有性质DM。

引理3[13]F(z,w)Rl×l[z,w],rank F(z,f(z))=l-1,对每个(z,w)R2,那么可以构造一个幺模向量Y(z)Rl×1[z],使得

F(z,f(z))Y(z)=0,,0T

2  主要结果

引理4F(z,w)Rl×l[z,w],且detF(z,w)=(z-f(w))q(w-a)q,aR,q是一个正整数。如果F(z,w)具有性质DM,那么F(z,w)等价于矩阵

V(z,w)P1(z,w)V1(w)P2(z,w)U1(z)P1(z,w)Vq(w)P2(z,w)Uq(z)

其中,P1(z,w)=Il-10l-1,101,l-1z-f(w)P2(z,w)=Il-10l-1,101,l-1w-aV(z,w)Gll(R[z,w])Vi(w)

Gll(R[w]),Ui(z)Gll(R[z]),i=1,,q

因为detF(z,w)F(z,w)的所有(l-1)×(l-1)级子式没有共同零点,所以w-aF(z,w)的所有(l-1)×(l-1)级子式没有共同零点。因为detF(z,w)=(z-f(w))q(w-a)q,aR[z,w],q是一个正整数,那么w-adetF(z,w)的一个因子。根据推论1,对每个zC,有rank F(z,a)=l-1。由引理3,存在一个幺模向量Y1(z)Rl×1[z],使得

F(z,a)Y1(z)=(0,,0)T

根据Quillen-Suslin定理[14],可以构造一个幺模矩阵Uq1(z)Rl×l[z],使得Y1(z)是它的最后一列且detUq1(z)=1。由(1)式得,存在B1(z)Rl×(l-1)[z],

F(z,a)Uq1(z)=(B10l,1)

由(2)式计算得

F(z,w)Uq1(z)=F1(z,w)P2(z,w)

其中,F1(z,w)Rl×l[z,w],P2(z,w)=Il-10l-1,101,l-1w-a。那么有

F(z,w)=F1(z,w)P2(z,w)Uq(z)

其中,Uq(z)=Uq1-1(z)。注意到detUq1(z)=1,所以detUq(z)=1。由于detF(z,w)=(z-f(w))q(w-a)q,那么detF1(z,w)=(z-f(w))q(w-a)q-1。由引理2有,F1(z,w)也具有性质DM。所以z-f(w)F1(z,w)的所有(l-1)×(l-1)级子式没有共同零点。再根据引理1,得到:wC,rank F(f(w),w)=l-1。结合引理3,存在一个幺模列向量Y2(w)Rl×1[w]使得

F(f(w),w)Y2(w)=(0,,0)T

由Quillen-Suslin定理,可以构造一个幺模矩阵Vq1(w)Rl×l[w],使得Y2(w)是它的最后一列且detVq1(w)=1。由(5)式得,存在B2(w)Rl×(l-1)[w],使得

F(f(w),w)Vq1(z)=(B20l,1)

由(6)式,容易得

F1(z,w)=F2(z,w)P1(z,w)Vq(w)

其中,Vq(z)=Vq1-1(z)。注意到detVq1(z)=1,所以detVq(z)=1。由于detF1(z,w)=(z-f(w))q(w-a)q-1,则detF2(z,w)=(z-f(w))q-1(w-a)q-1。根据引理2,F2(z,w)也具有性质DM,重复以上过程,可以得到我们的结论。证毕。

定理1P1(z,w),P2(z,w)Rl×l[z,w],且aR,V(w)Gll(R[z]),F(z,w)=P1(z,w)V(w)

P2(z,w)。如果F(z,w)具有性质DM,则F(z,w)等价于矩阵

P(z,w)=Il-10l-1,101,l-1(z-f(w))(w-a)

l进行归纳证明,当l=2时,令V(w)=v11(w)v12(w)v21(w)v22(w)Gl2(R[z,w]),

F(z,w)=v11(w)v12(w)(w-a)v21(w)(z-f(w))v22(w)(z-f(w))(w-a)

detF(z,w)=(z-f(w))(w-a)。首先,w-a不是v11(w)的零点,否则detF(z,w)v11(w),v12(w)(w-a),v21(w)(z-f(w)),v22(w)(z-f(w))(w-a)有一个公共零点(f(a),a),与F(z,w)具有性质DM矛盾。注意到V(w)Gl2(R[w])v11(w),v12(w)没有公共零点。所以v11(w),v12(w)(w-a)没有公共零点,v11(w),v12(w)(w-a)是一个幺模行。根据Quillen-Suslin定理,可以构造一个幺模矩阵M1(z,w)Gl2(R[z,w]),使得v11(w),v12(w)(w-a)M1(z,w)=(10),即

F(z,w)M1(z,w)=10p1(z,w)p2(z,w)

那么存在M2(z,w)=10-p1(z,w)1,使得M2(z,w)F(z,w)M1(z,w)=100p2(z,w)

因为M1(z,w),M2(z,w)Gl2(R[z,w]),所以p2(z,w)=udetF(z,w)=u(z-f(w))(w-a),

0uR。因此,F(z,w)等价于矩阵

P(z,w)=100(z-f(w))(w-a)

假设l-1时结论是成立的,考虑l的情形。

V1(w)=V11(w)V12(w)V21(w)V22(w)Gll(R[w]),其中,V11(w)R(l-1)×(l-1)[w]V12(w)R(l-1)×1[w],V21(w)R1×(l-1)[w]V22(w)R1×1[w]。那么

F(z,w)=V11(w)V12(w)(w-a)V21(w)(z-f(w))V22(w)(z-f(w))(w-a)

v11(w),,v1,l-1(w),v1l(w)是矩阵V1(w)的第一行。由V1(w)Gll(R[w]),v11(w),,

v1,l-1(w),v1l(w)没有共同零点。进一步可得,v11(w),,v1,l-1(w),v1l(w)(w-a)R[w]上的幺模行(否则,假设v11(w),,v1,l-1(w),v1l(w)(w-a)有一个共同零点,这个共同零点一定是w=a,那么detF(z,w)F(z,w)的所有(l-1)×(l-1)级子式有一个共同零点(f(a),a),与F(z,w)具有性质DM矛盾)。因此,根据Quillen-Suslin定理,可以构造一个幺模矩阵M1(z,w),M2(z,w)Gll(R[z,w])使得

N(z,w)M2(z,w)F(z,w)M1(z,w)=101,l-10l-1,1F2(z,w)

显然,detF2(z,w)=detN(z,w)=u(z-f(w))(w-a),0uR。由于F(z,w)具有性质DM,由推论2,知N(z,w)也具有性质DM,即detN(z,w)N(z,w)的所有(l-1)×(l-1)级子式没有共同零点。由于N(z,w)的一个(l-1)×(l-1)级子式或是F2(z,w)的一个(l-2)×(l-2)级子式或是detF2(z,w)或是0,则detF2(z,w)F2(z,w)的所有(l-2)×(l-2)级子式没有共同零点。由归纳假设,存在N1(z,w),N2(z,w)Gll-1(R[z,w])使得

N2(z,w)F2(z,w)N1(z,w)=Il-20l-2,101,l-2(z-f(w))(w-a)

所以

101,l-10l-1,1N2(z,w)N(z,w)101,l-10l-1,1N1(z,w)=101,l-10l-1,1N2(z,w)F2(z,w)N1(z,w)=Il-10l-1,101,l-1(z-f(w))(w-a)

因此,N(z,w)等价于矩阵

P(z,w)=Il-10l-1,101,l-1(z-f(w))(w-a)

结合F(z,w)等价于矩阵N(z,w),由等价关系的传递性,得到F(z,w)等价于矩阵P(z,w)。证毕。

定理2F(z,w)Rl×l[z,w]detF(z,w)=(z-f(w))(w-a),其中aR。则F(z,w)等价于矩阵

P(z,w)=Il-10l-1,101,l-1(z-f(w))(w-a)

当且仅当F(z,w)具有性质DM。

充分性:由detF(z,w)=(z-f(w))(w-a)detF(z,w)=(z-f(w))(w-a)及引理4,F(z,w)等价于P1(z,w)V(w)P2(z,w),其中,V(w)Gll(R[z]),P1(z,w)=Il-10l-1,101,l-1z-f(w)P2(z,w)=Il-10l-1,101,l-1w-a。根据定理1,结论成立。

必要性:由P(z,w)的矩阵形式知,detP(z,w)P(z,w)的所有(l-1)×(l-1)级子式没有共同零点。由推论2,detF(z,w)F(z,w)的所有(l-1)×(l-1)级子式也没有共同零点。因此F(z,w)具有性质DM。证毕。

推论3F(z,w)Rl×l[z,w]detF(z,w)=(w-f(z))(z-a),其中aR。则F(z,w)等价于矩阵

P(z,w)=Il-10l-1,101,l-1(w-f(z))(z-a)

当且仅当F(z,w)具有性质DM。

引理5P(z,w)Rl×l[z,w]P(z,w)=Il-100g(z,w),其中g(z,w)R[z,w]。令V(z,w)

Gll(R[z,w])。如果F(z,w)=Ps(z,w)V(z,w)Pt(z,w)具有性质DM,那么F(z,w)等价于矩阵Ps+t(z,w),

其中,st是正整数。

V(z,w)=V11(z,w)V12(z,w)V21(z,w)V22(z,w),其中,V11(z,w)R(l-1)×(l-1)[z,w],V12(z,w)R(l-1)×1[z,w],

V21(z,w)R1×(l-1)[z,w],V22(z,w)R1×1[z,w]。则
F(z,w)=Ps(z,w)V(z,w)Pt(z,w)=V11(z,w)V12(z,w)gt(z,w)V21(z,w)gs(z,w)V22(z,w)gs+t(z,w)

所以detF(z,w)=gs+t(z,w)

首先,证明V11(z,w),V12(z,w)gt(z,w)是一个ZLP矩阵。注意到V(z,w)Gll(R[z,w]),所以V11(z,w),V12(z,w)是一个ZLP矩阵。令c1(z,w),,cl(z,w)表示矩阵V11(z,w),V12(z,w)的所有(l-1)×(l-1)级子式,则矩阵V11(z,w),V12(z,w)gt(z,w)的所有(l-1)×(l-1)级子式为c1(z,w),,cl(z,w)gt(z,w)。下面用反证法证明c1(z,w),,cl(z,w)gt(z,w)没有共同零点。假设c1(z,w),,cl(z,w)gt(z,w)有一个共同零点α0,注意到c1(z,w),,cl(z,w)没有公共零点,那么α0一定是c1(z,w)gt(z,w)的公共零点。注意到F(z,w)的最后一行都有一个公因子g(z,w),结合F(z,w)的所有(l-1)×(l-1)级子式要么为c1(z,w),要么包含因子g(z,w),以及detF(z,w)包含因子g(z,w),那么detF(z,w)F(z,w)的所有(l-1)×(l-1)级子式有一个共同零点α0,与F(z,w)具有性质DM矛盾,那么v11(w),v12(w)gt(z,w)R[z,w]上的幺模行。根据Quillen-Suslin定理,存在矩阵M1(z,w)Gll(R[z,w])使得V11(w),V12(w)gs(z,w)M1(z,w)=(Il-10l-1,1)。这样,

F(z,w)M1(z,w)=Il-10l-1,1p1(z,w)p2(z,w)

其中,p1(z,w)R1×(l-1)[z,w],p2(z,w)R1×1[z,w]。从而存在M2(z,w)=10-p1(z,w)1,使得

M2(z,w)F(z,w)M1(z,w)=100p2(z,w)

注意到M1(z,w),M2(z,w)Gll(R[z,w]),所以p2(z,w)=udetF(z,w)=u(z-f(w))(w-a),0uR。因此,F(z,w)等价于矩阵

Ps+t(z,w)=Il-10l-1,101,l-1gs+t(z,w)

证毕。

定理3F(z,w)Rl×l[z,w]detF(z,w)=(z-f(w))q(w-a)q,其中aR。如果F(z,w)具有性质DM,则F(z,w)等价于矩阵

Q(z,w)=Il-10l-1,101,l-1(z-f(w))q(w-a)q

由引理4,F(z,w)可以分解为

V(z,w)P1(z,w)V1(w)P2(z,w)U1(z)P1(z,w)Vq(w)P2(z,w)Uq(z),

其中,

P1(z,w)=Il-10l-1,101,l-1z-f(w),P2(z,w)=Il-10l-1,101,l-1w-a,

V(z,w)Gll(R[z,w]),Vi(w)Gll(R[w]),Ui(z)Gll(R[z]),i=1,,q。根据定理1,我们有

P1(z,w)Vi(w)P2(z,w)=Mi1(z,w)P(z,w)Mi2(z,w)

其中,Mi1(z,w),Mi2(z,w)Gll(R[z,w]),i=1,,q,P(z,w)=Il-10l-1,101,l-1(z-f(w))(w-a)。那么

F(z,w)=V(z,w)M11(z,w)P(z,w)M12(z,w)U1(z)Mq1(z,w)P(z,w)Mq2(z,w)Uq(z)
L0(z,w)P(z,w)L1(z,w)P(z,w)L2(z,w)P(z,w)Lq(z,w)

其中,L0(z,w)=V(z,w)M11(z,w)L1(z,w)=Mi2(z,w)Ui(z,w)Mi+1,1(z,w)Lq(z,w)=Mq2(z,w)Uq(z)i=1,,q-1。显然,Li(z,w)Gll(R[z,w])i=1,,q

由引理5,P(z,w)L1(z,w)P(z,w)等价于P2(z,w),即存在N1(z,w),N2(z,w)Gll(R[z,w])使得P(z,w)L1(z,w)P(z,w)=N1(z,w)P2(z,w)N2(z,w)。因此,

F(z,w)=L0(z,w)P(z,w)L1(z,w)P(z,w)L2(z,w)P(z,w)Lq(z,w)=L0(z,w)N1(z,w)P2(z,w)N2(z,w)L2(z,w)P(z,w)Lq(z,w)

q3,T2(z,w)=N2(z,w)L2(z,w),则T2(z,w)Gll(R[z,w])。由引理5,P2(z,w)T2(z,w)·P(z,w)等价于矩阵P3(z,w),即存在N3(z,w),N4(z,w)Gll(R[z,w])使得

P2(z,w)T2(z,w)P(z,w)=N3(z,w)P3(z,w)N4(z,w)

因此,

F(z,w)=L0(z,w)N1(z,w)P2(z,w)N2(z,w)L2(z,w)P(z,w)Lq(z,w)=L0(z,w)N1(z,w)P2(z,w)T2(z,w)P(z,w)L3(z,w)P(z,w)Lq(z,w)=L0(z,w)N1(z,w)N3(z,w)P3(z,w)N4(z,w)P(z,w)Lq(z,w)

重复以上过程,得到

F(z,w)=L0(z,w)N1(z,w)N3(z,w)N2q-1(z,w)Pq(z,w)N2q(z,w)Lq(z,w)

显然,L0(z,w)N1(z,w)N3(z,w)N2q-1(z,w),N2q(z,w)Lq(z,w)Gll(R[z,w])。因此,可以得到F(z,w)等价于矩阵

Q(z,w)=Il-10l-1,101,l-1(z-f(w))q(w-a)q

证毕。

定理4F(z,w)Rl×l[z,w]detF(z,w)=(z-f(w))q(w-a)q,其中aR。则F(z,w)等价于矩阵

Q(z,w)=Pq(z,w)=Il-10l-1,101,l-1(z-f(w))q(w-a)q

当且仅当F(z,w)具有性质DM。

充分性:根据定理3,结论成立。

必要性:由Q(z,w)的矩阵表达形式知,detQ(z,w)Q(z,w)的所有(l-1)×(l-1)级子式没有共同零点。由推论2,detF(z,w)F(z,w)的所有(l-1)×(l-1)级子式没有共同零点。因此,F(z,w)具有性质DM。证毕。

推论4F(z,w)Rl×l[z,w]detF(z,w)=(w-f(z))q(z-a)q,其中aR。则F(z,w)等价于矩阵

Q(z,w)=Il-10l-1,101,l-1(w-f(z))q(z-a)q

当且仅当F(z,w)具有性质DM。

注1 令F(z,w)Rl×l[z,w]detF(z,w)=(w-f(z))q(z-a)q,其中aR,q是一个正整数。定理4给了一个判别F(z,w)与其Smith型等价的充分必要条件。F(z,w)具有性质DM意味着detF(z,w)F(z,w)的所有(l-1)×(l-1)级子式生成单位理想。我们知道一个理想是单位理想R当且仅当该理想的既约Gröbner基是1。因此,这个条件可以通过计算detF(z,w)F(z,w)的所有(l-1)×(l-1)级子式所生成理想的既约Gröbner基进行检验。

3  例 子

本节通过一个例子来说明我们的构造性方法。应用该方法可以直接得到矩阵等价变换过程中的幺模矩阵。

考虑一个在R3×3[z,w]上的二元多项式矩阵

F(z,w)=w301w2z2+z2w2wz-4z2wz+2z21(w-1)2

经计算,detF(z,w)=(z-w2)2(w-1)2F(z,w)具有性质DM。考虑

F(z,1)=1012z21-2z2z+2z210

解方程:F(z,1)Y1(z)=0,得到Y1(z)T=(-1,2z2+2z,1)。将(-1,2z2+2z,1)T嵌入到一个幺模矩阵U21(z)R3×3[z]使得

F(z,w)U21(z)=F1(z,w)P2(z,w)

其中,

P2(z,w)=10001000w-1,U21(z,w)=10-1012z2+2z001

那么

F(z,w)=F1(z,w)P2(z,w)U2(z,w)

这里

U2(z,w)=U21-1(z,w)=10101-2z2-2z001

再考虑F1(z,w),重复以上过程,有

F1(z,w)=F2(z,w)P1(z,w)V2(w)

其中,

F2(z,w)=010-w4wz-4z+2wz2z+wz+2w2z+2w3z+w2+5w3+2w4+2w5-1w2-2w+2z2+2z+12w+2z+2wz+2w2z+4w2+2w3+2w4
V2(w)=w3(w-1)(2w3+2w2+5w+1)-1-w3(2w3+2w2+5w+1)w30-w2-w-11-w01
P1(z,w)=10001000z-w2

因此,F(z,w)=F1(z,w)P2(z,w)U2(z,w)=F2(z,w)P1(z,w)V2(w)P2(z,w)U2(z,w)

重复以上过程,得

F(z,w)=F4(z,w)P1(z,w)V1(w)P2(z,w)U1(z,w)P1(z,w)V2(w)P2(z,w)U2(z,w)

其中,

F4(z,w)=010-w4wz-4z+2wz22w2+4w-1-1w2-2w+2z2+2z+12w+4
V1(w)=10-4w3-8w2-8w-10010001,U1(z,w)=10-6z-10010001

由于detF4(z,w)=1,下面考虑P1(z,w)V1(w)P2(z,w)P1(z,w)V2(w)P2(z,w),经过计算可得

P1(z,w)V1(w)P2(z,w)=P(z,w)M12(z,w)

其中,

P(z,w)=10001000(z-w2)(w-1),M12(z,w)=10-(w-1)(4w3+8w2+8w+10)010001

P1(z,w)V2(w)P2(z,w)=P(z,w)M22(z,w)
M22(z,w)=w3(w-1)(2w3+2w2+5w+1)-1-w3(w-1)(2w3+2w2+5w+1)w301-w3-101

F(z,w)=F4(z,w)P1(z,w)V1(w)P2(z,w)U1(z,w)P1(z,w)V2(w)P2(z,w)U2(z,w)=F4(z,w)P(z,w)M12(z,w)U1(z,w)P(z,w)M22(z,w)U2(z,w)

L1(z,w)=M12(z,w)U1(z,w),L2(z,w)=M22(z,w)U2(z,w),则

F(z,w)=F4(z,w)P(z,w)L1(z,w)P(z,w)L2(z,w)

其中,

L1(z,w)=10-4w4-4w3-2w-6z010001,L2(z,w)=2w7+3w5-4w4-w3-12z2+2zw301-100

A(z,w)=P(z,w)L1(z,w)P(z,w)=10-(-w2+z)(w-1)(4w4+4w3+2w+6z)01000(-w2+z)2(w-1)2

A(z,w)进行初等变换可得A(z,w)=Q(z,w)N2(z,w),其中

Q(z,w)=10001000(z-w2)2(w-1)2,N2(z,w)=10-(-w2+z)(w-1)(4w4+4w3+2w+6z)010001

因此

F(z,w)=F4(z,w)P(z,w)L1(z,w)P(z,w)L2(z,w)=F4(z,w)Q(z,w)N2(z,w)L2(z,w)=T1(z,w)Q(z,w)T2(z,w)

其中,T1(z,w)=F4(z,w),T2(z,w)=N2(z,w)L2(z,w)。经计算,detT1(z,w)=detT2(z,w)=1

因此,F(z,w)等价于矩阵

Q(z,w)=10001000(z-w2)2(w-1)2

4  结 语

本文研究了多项式环R[z,w]上的几类二元多项式矩阵等价问题,其中R是实数域,得到了一些有意义的结果,并获得了几类二元多项式矩阵F(z,w)与其Smith型等价的判别条件。这些条件可以通过计算F(z,w)的行列式及其低一阶子式所生成的理想的既约Gröbner基进行检验。文中我们给出了一个例子对我们的方法进行说明,该分解方法可以在Singular系统中实现,这些都有利于工程中二维系统的简化。

参考文献

[1]

BOSE N K. Applied Multidimensional Systems Theory [M]. New York: Van Nostrand Reinhold, 1982. DOI: 10.1007/978-3-319-46825-9 .

[2]

BOSE N K. Multidimensional Systems Theory and Applications [M]. Dordrecht: Kluwer, 2003. DOI: 10.1007/978-94-017-0275-1 .

[3]

DREESEN P, BATSELIER K, BART D M. Multidimensional realisation theory and polynomial system solving [J]. International Journal of Control, 2018, 91(12): 2692-2704. DOI: 10.1080/00207179.2017.1378924 .

[4]

BOUDELLIOUA M S, QUADRAT A. Serre’s reduction of linear function systems [J]. Mathematics in Computer Science, 2010, 4(2): 289-312. DOI: 10.1007/s11786-010-0057-y .

[5]

ROSENBROCK H H. State Space and Multivariable Theory [M]. New York:Nelson-WiIey, 1970.

[6]

VAFIADIS D, KARCANIAS N. Unimodular equivalence and similarity for linear systems [J]. International Journal of Control, 2019, 92(9): 2091-2098. DOI: 10.1080/00207179.2018.1427892 .

[7]

BOUDELLIOUA M S, GALKOWSKI K, ROGERS E. Characterization of a class of spatially interconnected systems (ladder circuits) using two-dimensional systems theory [J]. Multidimensional Systems and Signal Processing, 2019, 30: 2185-2197. DOI: 10.1007/s11045-019-00644-9 .

[8]

BOUDELLIOUA M S. Computation of the Smith form for multivariate polynomial matrices using maple [J]. American Journal of Computational Mathematics, 2012. 2: 21-26. DOI: 10.4236/ajcm.2012.21003 .

[9]

FROST M G, STOREY C. Equivalence of a matrix over R[s, z] with its Smith form [J]. International Journal of Control, 1978, 28(5): 665-671. DOI: 10.1080/00207177808922487 .

[10]

LEE E, ZAK S. Smith form over R[z1,z2] [J]. IEEE Trans Autom Control, 1983, 28(1): 115-118. DOI: 10.1109/TAC.1983.1103118 .

[11]

FROST M G, BOUDELLIOUA M S. Some further results concerning matrices with elements in a polynomial ring [J]. International Journal of Control, 1986, 43(5): 1543-1555. DOI: 10.1080/00207178608933558 .

[12]

LIN Z P, BOUDELLIOUA M S, XU L. On the equivalence and factorization of multivariate polynomial matrices [J]. 2006 IEEE International Symposium on Circuits and Systems (ISCAS). Piscataway:IEEE,2006: 4911-4914. DOI: 10.1109/ISCAS.2006.1693732 .

[13]

LI D M, LIU J W, ZHENG L C. On the equivalence of multivariate polynomial matrix [J]. Multidim Syst Sign Process, 2017, 28: 225-235. DOI: 10.1007/s11045-015-0329-4 .

[14]

QUILLEN D. Projective modules over polynomial rings [J]. Invent Math, 1976, 36: 167-171. DOI: 10.1007/BF01390008 .

[15]

LIU J W, WANG M S. New results on multivariate polynomial matrix factorizations [J]. Linear Algebra and Its Applications, 2013, 438: 87-95. DOI: 10.1016/j.laa.2012.08.012 .

基金资助

国家自然科学基金(11871207)

AI Summary AI Mindmap
PDF (491KB)

0

访问

0

被引

详细

导航
相关文章

AI思维导图

/