蘇 瑞, 潘寶珍

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

內(nèi)積空間上向量值 Padé-型逼近表的塊狀結(jié)構(gòu)特征

蘇 瑞, 潘寶珍

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

從多項(xiàng)式空間到向量空間引入一種廣義線性泛函,在內(nèi)積空間上定義和構(gòu)造向量值 Padé-型逼近.借助向量值 Padé-型逼近的誤差公式,給出關(guān)于線性泛函的正交多項(xiàng)式的定義,同時(shí)推導(dǎo)出向量值 Padé-型逼近表的塊狀結(jié)構(gòu)特征.利用 Padé-型逼近表的這一特征,可以減少向量值 Padé-型逼近的計(jì)算量.最后,通過數(shù)值實(shí)例說明該方法的有效性.

向量值;Padé-型逼近;正交多項(xiàng)式;塊狀結(jié)構(gòu)特征

Abstract:A vector-valued Padé-type approximation is defined in the inner space by introducing a generalized linear functional from a polynomial space to a vector space.With the error formula for vectorvalued Padé-type app roximation,the orthogonal polynomial w ith respect to a generalized linear functional is defined.The block structure of Padé-type table is derived. The structure can be used to reduce computation of Padé-type approximations.An example is given to illustrate effectivenessof themethod.

Key words:vector-valued;Padé-type approximation;orthogonal polynomial;block structure

Padé-型逼近在理論物理、系統(tǒng)控制理論、模型簡化和積分方程領(lǐng)域都有著廣泛的應(yīng)用.從 1979年法國數(shù)學(xué)家 Brezinski[1]研究了數(shù)量 Padé-型逼近開始,許多數(shù)學(xué)工作者將 Padé-型逼近理論加以發(fā)展.1983年,Draux[2]將 Padé-型逼近從數(shù)量情形推廣到非交換代數(shù)的情形,并提出矩陣 Padé-型逼近.1999年,Salam[3]將 Padé-型逼近推廣到向量情形,向量值Padé-型逼近是借助于 Clifford代數(shù)的方法來定義,該方法在具體計(jì)算中很難實(shí)現(xiàn).2004年,顧傳青[4]引入一種從多項(xiàng)式空間到矩陣空間上的廣義線性泛函,從而在矩陣內(nèi)積的基礎(chǔ)上構(gòu)造和定義了矩陣值Padé-型逼近.

本研究借助于向量值 Padé-型逼近的誤差公式及高階向量值 Padé-型逼近[5-7],推導(dǎo)出向量值 Padé-型逼近表的塊狀結(jié)構(gòu)特征,即當(dāng)向量值 Padé-型逼近的生成多項(xiàng)式取為關(guān)于線性泛函的正交多項(xiàng)式時(shí),塊狀結(jié)構(gòu)特征表中的逼近元素及逼近階是完全相同的,可以利用該特征在計(jì)算高階向量值 Padé-型逼近時(shí)減少計(jì)算量.

1 向量值 Padé-型逼近的定義與構(gòu)造

向量 a=(a1,a2,…,ad)∈Cd,b=(b1,b2,…,bd)∈Cd,ai,bi∈C.它們的內(nèi)積定義為

定義向量 g=(g1,g2,…,gd)∈Cd的范數(shù)

這里 g＊表示 g的復(fù)共軛.

向量 g的廣義逆 (Samelson逆)定義為

設(shè) P是一元實(shí)系數(shù)多項(xiàng)式,Pk表示 P中次數(shù)不超過 k的多項(xiàng)式的集合.設(shè)向量 g形式冪級(jí)數(shù)為

式中 ,ci∈Cd,z∈C.

設(shè) φ:P→Cd為作用在多項(xiàng)式空間到向量空間上的線性泛函,定義為

如果 |xz|<1,則有 (1-xz)-1=1+xz+(xz)2+….將線性泛函作用在 (1-xz)-1上,可得

設(shè) v∈Pn是次數(shù)為 n的數(shù)量多項(xiàng)式,有

定義具有向量值系數(shù)的多項(xiàng)式W(z)為

注意到 φ是作用于多項(xiàng)式空間的廣義線性泛函,因而W(z)是關(guān)于 z的次數(shù)為 n-1的具有向量值系數(shù)的多項(xiàng)式.再令

定理1[5]設(shè)v～(0)≠0,則成立

設(shè) φ(l):P→Cd是作用于多項(xiàng)式空間到向量空間的廣義線性泛函,定義為

當(dāng) l+k<0時(shí) ,規(guī)定 φ(l)(xk)=0.

根據(jù)構(gòu)造公式 (4)～(7),定義

定理2[5]設(shè)v～(0)≠0,則成立Pmn(z)/

定義 2 定義向量值函數(shù) Rm,n(z)=Pmn(z)/為(m,n)階向量值Padé-型逼近,記為(m/n)f(z).

2 向量值 Padé-型逼近表的結(jié)構(gòu)特征

由向量值 Padé-型逼近誤差公式[5],有

根據(jù) (m/n)f(z)的構(gòu)造方式可知,(m/n)f(z)實(shí)際上取決于 n個(gè)任意常數(shù),為此令

定義 3 滿足方程 (10)的數(shù)量多項(xiàng)式 v(x)定義為關(guān)于廣義線性泛函 φ(m-n+1)的正交多項(xiàng)式,由正交多項(xiàng)式 v(x)所確定的 (m/n)f(z),稱為給定向量冪級(jí)數(shù) f(z)的高階向量 Padé-型逼近.

若在式 (10)中代入 v(z)=b0+b1z+…+bnzn,施加廣義線性泛函 φ(m-n+1),并與向量 cm-n+1分別作內(nèi)積,得到下列線性方程組:

通過上述方程組已完整建立向量值 Padé-型逼近的行列式表達(dá)式[5].本研究在此基礎(chǔ)上討論由正交關(guān)系 (10)確定的生成多項(xiàng)式所構(gòu)造的向量值Padé-型逼近表的分布特征.

定理 3 設(shè) (1)vn(x)滿足

而 up(z)是任一關(guān)于 z的 p次多項(xiàng)式 (0≤p≤n);(2)(m/n)f(z)是 f(z)的以 vn(z)為生成多項(xiàng)式的向量 Padé-型逼近,則向量冪級(jí)數(shù) f(z)的以 v(z)=up(z)vn(z)為生成多項(xiàng)式的 (m+p+h/n+p)f(z)向量值 Padé-型逼近滿足

式中,0≤p+h≤n.

證明 由高階向量值 Padé-型逼近的定義得

將式中 m換為 m+p+h,n換為 n+p,以 vn(z)up(z)為生成多項(xiàng)式,構(gòu)造的向量值 Padé-型逼近,有

由式 (14),可推出

注意到泛函只作用在 x上,觀察式 (13)發(fā)現(xiàn),W′(z)實(shí)際上是關(guān)于 z的 n+p-1次多項(xiàng)式,即

于是,有

注解 1 只有在給定生成多項(xiàng)式且此生成多項(xiàng)式是有正交關(guān)系確定的條件下,定理 3成立.

對(duì)方程組 (18)中的每一個(gè)方程分別與向量 c2作內(nèi)積,所求得的形式正交多項(xiàng)式為 v2(z).記 (3/2)f(z)是以 v2(z)為生成多項(xiàng)式的 Padé-型逼近 ,則由定理 3得

式中,(3+p+h/2+p)f(z)是以 up(z)v2(z)為生成多項(xiàng)式的 Padé-型逼近,up(z)是任意 p次多項(xiàng)式 (0≤p≤2).如此得到向量值 Padé-型逼近表塊狀結(jié)構(gòu)特征,即如下的逼近元素相同(逼近階也完全相同):

注解 2 若要計(jì)算 (5/4)f(z),根據(jù)注解 1,只要計(jì)算 (3/2)f(z).這樣減少了向量值 Padé-型逼近的計(jì)算量,但逼近階卻達(dá)到 O(z6).

例 1 考察控制論中的向量指數(shù)函數(shù)

解 當(dāng)把向量值 Padé-型逼近 (5/4)f(z)的生成多項(xiàng)式取為正交多項(xiàng)式v～2(z)時(shí),根據(jù)定理2,只需計(jì)算 (3/2)f(z),其分子、分母分別計(jì)算[5]如下:易驗(yàn)證(3/2)f(z)=(5/4)f(z)=P32(z)/v～2(z).

[1] BREZINSKI C. Padé-type approximation and general orthogonal polynomials[M].Basel:Birkh?user,1980.

[2] DRAUX A.Approximants de type Padéet de Padé[M].Lille:Universiédes Science et Technologies de Lille,1983:1-89.

[3] SALAM A.Vector Padé-type approximants and vector PadéApp roximants[J].JApp roximation Theory,1999,97(1):92-112.

[4] GU C Q.Matrix Padé-type approximant and directional matrix in the inner p roduct space[J].J Comput Appl Math,2004,164/165:365-385.

[5] L IC J.A new vector-valued Padé-type approximation in the inner space[J].Numer Math J Chinese Univ:English Series,2006,15(2):127-136.

[6] 潘寶珍.用于積分方程解的函數(shù)值 Padé-型逼近的恒等式與遞推算法 [J].應(yīng)用科學(xué)學(xué)報(bào),2006,24(1):74-77.

[7] GU CQ,PAN B Z,WU B B.Orthogonal polynomials and determinant formulas of function-valued Padé-type approximation using for solution of integral equations[J].Appl Math Mech:Engl Ed,2006,27(6):853-860.

(編輯:孟慶勛)

Block Structure of Vector-Valued Padé-Type Table in the Inner Space

SU Rui, PAN Bao-zhen
(College of Sciences,ShanghaiUniversity,Shanghai200444,China)

O 241.83

A

1007-2861(2010)03-0253-04

10.3969/j.issn.1007-2861.2010.03.007

2008-12-05

上海市重點(diǎn)學(xué)科建設(shè)資助項(xiàng)目 (J50101)

潘寶珍 (1965～),女,副教授,博士,研究方向?yàn)榉汉治?E-mail:bzpan@staff.shu.edu.cn