顧傳青， 張 科， 王 鋒

(上海大學(xué)理學(xué)院，上海200444)

矩陣有理插值問題包括矩陣冪級(jí)數(shù)的部分實(shí)現(xiàn)問題和Newton-Padé，Hermite-Padé，多點(diǎn)Padé逼近及其矩陣形式的推廣等問題［1］.已有的工作研究了利用相同插值點(diǎn)的矩陣有理逼近問題.Antoulas等［2］解決了最小矩陣有理插值問題.Anderson等［3］利用Loewner矩陣，研究了實(shí)有理傳遞函數(shù)矩陣的插值數(shù)據(jù)，傳到傳遞函數(shù)矩陣的一個(gè)最小狀態(tài)變量的實(shí)現(xiàn)問題.Gu［4］給出了一種廣義逆矩陣有理插值(generalized inverse matrix rationalinterpolation，GMRI)法，該方法推廣并提高了Graves-Morris等［5-6］的廣義逆向量有理插值法，使得構(gòu)造過程中不必計(jì)算矩陣的乘法.但是GMRI法受次數(shù)和可除性的限制，構(gòu)造的有理插值式中的分母多項(xiàng)式的次數(shù)不能是奇數(shù)［4，7-8］.利用牛頓多項(xiàng)式，Sidi［9-10］給出了3種不同的向量插值型有理插值式，其分子為向量值多項(xiàng)式，分母為標(biāo)量多項(xiàng)式.

本研究的目的是為了改進(jìn)GMRI法，通過構(gòu)造一種新的以拉格朗日多項(xiàng)式為基礎(chǔ)的迭代關(guān)系，提出了內(nèi)積空間上一種基于拉格朗日型矩陣有理插值的實(shí)用方法.利用矩陣內(nèi)積，給出3種不同的標(biāo)準(zhǔn)來得到3種有理插值方法.這類方法最早由Sidi［9-10］提出并討論了牛頓向量值有理插值問題.在此基礎(chǔ)上，Gu［11］研究了矩陣Padé型逼近問題.

1 矩陣值有理插值的一般方法

設(shè)z為一復(fù)變量，F(xiàn)(z)為矩陣值函數(shù)，即F：C→ Cs×t.假定F(z)定義在有界開集Ω上，現(xiàn)考慮F(z)在插值節(jié)點(diǎn)z1，z2，…處的插值問題，假定這些插值點(diǎn)是互異的.

令Lm，n(z)為函數(shù)F(z)在點(diǎn)zm，zm+1，…，zn處的矩陣值插值多項(xiàng)式，其次數(shù)最多為n-m.利用拉格朗日多項(xiàng)式，有

式中，lm，i(z)為拉格朗日基函數(shù)，即

顯然，Lm，n(z)∈Cs×t.為簡(jiǎn)潔起見，將數(shù)量多項(xiàng)式Wm，n(z)表述為

利用上述符號(hào)，定義一類矩陣值有理函數(shù)Rl，k(z)為

式中，c0，c1，…，ck為任意復(fù)數(shù)，Pl，k(z)為一個(gè)次數(shù)不超過l-1的矩陣值多項(xiàng)式，Ql，k(z)為一個(gè)最多k次的數(shù)量多項(xiàng)式.注意到，只要 Ql，k(z1)≠0，可將Ql，k(z)規(guī)范化，使得Ql，k(z1)=c0=1.相比文獻(xiàn)［4，7-8］中的GMRI法，式(4)中分母多項(xiàng)式的次數(shù)沒有奇偶次限制.

定理1 假定zi互不相同.若Ql，k(zi)≠0，i=1，2，…，l，則式(4)中定義的矩陣值有理函數(shù)Rl，k滿足如下的插值條件：

證明 由式(4)可得

令z=zi，由Lj+1，l(zi)=F(zi)，i=j+1，j+2，…，l以及W1，j(zi)=0，i=1，2，…，j，有

W1，j(z)［F(z)-Lj+1，l(z)］=0，i=1，2，…，l.定理得證.

2 誤差公式及系數(shù)選擇

由上節(jié)不難發(fā)現(xiàn)，式(4)中的ci是任意的.基于隨后給出的引理3中的迭代關(guān)系，本節(jié)討論怎樣選取合適的cj來獲得較好的逼近效果.首先介紹矩陣內(nèi)積的通常定義.

令矩陣A=(ai，j)，B=(bi，j)∈Cs×t，定義A，B的內(nèi)積為

引理1 令A(yù)，B和C∈Cs×t，則

(b)任給α，β∈C，(αA+βB，C)=α(A，C)+ β(B，C);

引理2［13］若 f(x)在包含著插值節(jié)點(diǎn) z0，z1，…，zn的區(qū)間［a，b］上n+1次可微，則對(duì)任意z∈［a，b］，存在與x有關(guān)的ξ(a＜ξ＜b)，使得

引理3 令Lm，n(z)為式(1)中定義的F(z)關(guān)于節(jié)點(diǎn)zm，zm+1，…，zn的矩陣值拉格朗日多項(xiàng)式，可得如下的迭代關(guān)系：

式中，Wm，n(z)由式(3)給出.

證明 由式(1)可得

由式(1)和(9)可推出F(zi)的系數(shù)，即

式中，m≤i≤n.式(8)得證.令

則式(8)變?yōu)?/p>

定理2 令z為一個(gè)復(fù)變量，F(xiàn)(z)為一個(gè)矩陣值函數(shù)F：C→Cs×t.若Rl，k(z)為式(4)給出的矩陣值有理函數(shù)，則

式中，

且Wm，n(z)，Dm，n分別由式(3)和(10)給出.

證明 由于Lm，n(z)為F(z)關(guān)于節(jié)點(diǎn) zm，zm+1，…，zn的插值多項(xiàng)式，由式(8)可得誤差公式為

令m=j+1，由引理3和式(11)，有

式中，0≤j≤k，n≥l.將式(15)代入式(5)，即可得到式(12).

根據(jù)式(3)，將式(13)改寫為

可以發(fā)現(xiàn)，式(16)方括號(hào)內(nèi)的Di，j，i=1，2，…，k+1，j=l+1，l+2，…，n，可以重排為如下格式：

本研究通過3種不同的方法選取cj，使得式(16)方括號(hào)內(nèi)的項(xiàng)在某種意義下盡可能的小.這里，對(duì)c0進(jìn)行規(guī)范化，即c0=1.

的解，其中假定 D2，l+1，D3，l+1，…，Dk+1，l+1為非零矩陣.于是將問題轉(zhuǎn)化為F-范數(shù)下的最小二乘求解問題.事實(shí)上，式(18)等價(jià)于求解如下的線性系統(tǒng)：

特別地，當(dāng)Di，l+1∈Rs×t，i=1，2，…，k+1時(shí)，式(19)變?yōu)?/p>

方法2 設(shè)矩陣D為D1，l+1，D2，l+1，…，Dk+1，l+1; D1，l+2，D2，l+2，…，Dk+1，l+2;…;D1，l+k，D2，l+k，…，Dk+1，l+k中的第一個(gè)非零矩陣，則通過求解下列線性系統(tǒng)來確定c1，c2，…，ck：

方法3 令式(16)中的c1，c2，…，ck為下列線性系統(tǒng)的解［12］：

3 指定極點(diǎn)的有理插值式的插值節(jié)點(diǎn)選取原則

定義1 給定序列：

定義f(z)關(guān)于點(diǎn)βn1，βn2，…，βnn的插值函數(shù)序列為

定義2［14］若Γ為中的閉Jordan曲線，則ext(Γ)為CΓ的包含C+的連通分支.

定義3［14］設(shè) f(z)為 G?中的解析函數(shù)，Γ(∞∈Γ)={z：Re z=-a}(a＞0)為G∪{∞}中的可求長 Jordan曲線.其中為給定的多項(xiàng)式序列，其中αni?{Re z≥-a}，βni∈且不同于αni.若

且

則指定極點(diǎn)的插值函數(shù)序列(rational interpolants with prescribed poles，RIPPs)rn(z)=Bn(z)/Qn(z)一致收斂到f(z).

基于插值公式(4)，給出逼近有指定極點(diǎn)的矩陣值函數(shù)的算法.假定F(z)為有指定極點(diǎn)的矩陣值函數(shù)，則第一步為計(jì)算F(z)的極點(diǎn);第二步為利用式(25)確定插值節(jié)點(diǎn)，去除第一步中不重要的極點(diǎn);第三步為利用插值式(4)逼近函數(shù)F(z).

例 考慮如下矩陣值有理函數(shù)的插值問題：

式中，

函數(shù) F(z)的極點(diǎn)為 -0.243 747 876 5，-2.373 140 546，-27.890 485 68，-36.717 000 60，-48.848 388 70，-71.291 618 30+636.280 518 7i，-71.291 618 30-636.280 518 7i.

由式(25)，選取插值節(jié)點(diǎn)z1=0.243 747 876 5，z2= 2.373 140 546，z3= 27.890 485 68，z4= 36.717 000 60，z5=48.848 388 70.根據(jù)插值公式(4)，可得插值函數(shù)R4，3(z)為

式中，

系數(shù)cj可由第3節(jié)介紹的3種方法求出，這里以方法1為例.由式(19)可得

Dm，n由式(10)給出，即

由此，計(jì)算式(27)，有

由于R4，3(z)和F(z)均為2階矩陣，給出矩陣R4，3(z)中4個(gè)元素逼近函數(shù)F(z)相應(yīng)元素的圖形，如圖1～圖 4所示，其中實(shí)線為 R4，3(z)，點(diǎn)線為F(z).可以看出，在區(qū)間［1，∞)上，插值函數(shù)R4，3(z)能較好地逼近被插函數(shù)F(z).特別地，在區(qū)間［50，∞)上，兩函數(shù)曲線誤差極小，這表明本研究提出的算法是有效的.

圖1 R4，3(z)和F(z)的第一行第一列的元素Fig.1 Element［1，1］of matrices R4，3(z)and F(z)

圖2 R4，3(z)和F(z)的第一行第二列的元素Fig.2 Element［1，2］of matrices R4，3(z)and F(z)

圖3 R4，3(z)和F(z)的第二行第一列的元素Fig.3 Element［2，1］of matrices R4，3(z)and F(z)

圖4 R4，3(z)和F(z)的第二行第二列的元素Fig.4 Element［2，2］of matrices R4，3(z)and F(z)

［1］ BECKERMANNB，LABAHNG.Recursiveness in matrix rational interpolation problems［J］.J Comput Appl Math，1997，77(1/2)：5-34.

［2］ ANTOULASA C，BALLJ A，KANGJ，et al.On the solution of the minimal rational interpolation problem［J］.Linear Algebra Appl，1990，137：511-573.

［3］ ANDERSON B D O，ANTOULASA C.Rational interpolation and state-variable realizationas［J］.Linear Algebra Appl，1990，137：479-509.

［4］ GUC Q.Thiele-type and Lagrange-type generalized inverse rational interpolation for rectangular complex matrices［J］.Linear Algebra Appl，1999，295：7-30.

［5］ GRAVES-MORRIS P R. Vector valued rational interpolantsⅠ［J］.Numer Math，1983，42：331-348.

［6］ GRAVES-MORRISP R，JENKINSC D.Vector valued rational interpolantsⅢ［J］.Constr Approx，1986，2：263-289.

［7］ GUC Q.Generalized inverse matrix Padé approximation on the basis of scalar product［J］.Linear Algebra Appl，2001，322：141-167.

［8］ GUC Q.A practical two-dimensional thiele-type matrix Padé approximation［J］.IEEE Trans Automat Control，2003，48：2259-2263.

［9］ SIDIA.A new approach tovector-valued rational interpolation［J］.Journal of Approximation Theory，2004，130(2)：179-189.

［10］ SIDIA.Rational approximations from power series of vector-valued meromorphic functions［J］.Journal of Approximation Theory，1994，77(1)：89-111.

［11］ GUC Q.Matrix Padé-type approximant and directional matrix Padé approximant in the inner product space［J］.J Comput Appl Math，2004，164：365-385.

［12］ GUC Q，WU B B.Some practical determinants of matrix Padé-type method and computation of matrix exponential［J］.J Information and Computational Science，2005，2 (2)：243-250.

［13］ 王仁宏.數(shù)值逼近［M］.北京：高等教育出版社，1999.

［14］ STANFORDA R.Approximation of transfer functions of infinite dimensional dynamical systems by rationl interpolants with prescribed poles［J］.Journalof Mathematical Analysis and Applications，2000，244 (1)：147-168.