宋巨龍，錢(qián)富才，梁錦錦

(1.西安石油大學(xué) 理學(xué)院,陜西 西安 710065；2.西安理工大學(xué) 自動(dòng)化與信息工程學(xué)院，陜西 西安 710048)

在非線性控制領(lǐng)域中，二次型最優(yōu)控制方法是一個(gè)值得研究的方向，文獻(xiàn)[1] 提出了非齊次雙線性二次型最優(yōu)控制迭代算法、文獻(xiàn)[2] 基于二次型性能指標(biāo)對(duì)簡(jiǎn)單自適應(yīng)控制系統(tǒng)進(jìn)行了研究，文獻(xiàn)[3] 給出了一個(gè)二次型控制方法在風(fēng)量空調(diào)箱優(yōu)化方面的應(yīng)用。文獻(xiàn)[4] ～[6] 則對(duì)廣義二次型和類(lèi)二次型以及矩陣的性質(zhì)及應(yīng)用進(jìn)行了研究，文獻(xiàn)[7] 所論及的二次型設(shè)計(jì)方法則涉及二次型在軍事領(lǐng)域的應(yīng)用，此外二次型還在航空、自動(dòng)化、化工[7-9]等諸多領(lǐng)域中有著廣泛的應(yīng)用。在最優(yōu)化方法中，人們也常常將非線性函數(shù)展開(kāi)為二次函數(shù)進(jìn)行近似求解[10]，這里也涉及二次型。受上述文獻(xiàn)的啟發(fā)我們考慮如果能將非線性函數(shù)的三階泰勒多項(xiàng)式加以應(yīng)用，用三次型替代二次型應(yīng)該能夠得到更為精確的結(jié)果。但是目前較為一般的方法是把非線性函數(shù)利用其二次函數(shù)替代，這是因?yàn)槎嘣瘮?shù)三次以上的泰勒展開(kāi)式的形式比較復(fù)雜，從幾何意義上看不很明確，利用起來(lái)比較困難。比如文獻(xiàn)[11] 對(duì)多元函數(shù)的泰勒公式進(jìn)行了研究，但所使用的就是傳統(tǒng)的多元函數(shù)泰勒展開(kāi)式，顯得較為繁瑣而不易理解。此外注意到函數(shù)逼近論中很多情況下也是將一個(gè)函數(shù)用多項(xiàng)式來(lái)替代從而達(dá)到近似計(jì)算的目的[12]，如果我們能給出更高階的、幾何特征正直觀的多項(xiàng)式，也會(huì)使這種問(wèn)題得到更精確的結(jié)果。

本文根據(jù)目前的研究現(xiàn)狀，從二次型這一概念及其重要的應(yīng)用性出發(fā)，對(duì)二次型進(jìn)行一個(gè)數(shù)學(xué)意義上的擴(kuò)展，提出了一個(gè)新的概念——三次型。對(duì)二次型的運(yùn)算方式進(jìn)行稍稍改變，并由此出發(fā)給出三次型的定義和運(yùn)算方法。從向量、矩陣引申出體陣的概念；從二次型擴(kuò)展到三次型。還將傳統(tǒng)的二次型和新提出來(lái)的三次型在形式上、運(yùn)算方式上進(jìn)行了統(tǒng)一，給出了三次型的各階導(dǎo)數(shù)。

1 二次型

二次型的定義[4]：

眾所周知，二次型由于其形式特殊、性質(zhì)特殊、應(yīng)用范圍廣泛，因而極具研究?jī)r(jià)值。

推而廣之，不難寫(xiě)出n元k次齊次多項(xiàng)式。

相信很多人也想過(guò)這個(gè)問(wèn)題，只不過(guò)對(duì)其理論意義和應(yīng)用價(jià)值無(wú)法確定，另外由于其表現(xiàn)形式已經(jīng)比較復(fù)雜了，失去了二次型通過(guò)矩陣運(yùn)算相對(duì)簡(jiǎn)單的特性，因而如何使其運(yùn)算簡(jiǎn)單化，也是需要考慮的一個(gè)問(wèn)題。下面對(duì)這些問(wèn)題逐一進(jìn)行探討。

2 三次型的定義

為了使得后面所討論的問(wèn)題能夠有一個(gè)比較簡(jiǎn)單的形式，先給出一些定義。

定義2稱(chēng)由n個(gè)數(shù)a1,a2,…,an組成的n維數(shù)組為一個(gè)1度n維體陣，記為：

即為通常意義上的n維向量，因而1度n維體陣也和向量看作是相同的。

定義3稱(chēng)由n2個(gè)數(shù)aij(i,j=1,2,…,n)組成的正方形數(shù)表為一個(gè)2度n維體陣，記為：

其中aij=aji(i≠j,i,j=1,2,…,n)。即通常意義上的n階方陣。

定義4稱(chēng)由n3個(gè)數(shù)aijk(i,j,k=1,2,…,n)組成的立方體數(shù)陣為一個(gè)3度n維體陣，記為：

其中:

aijk=aikj=akji=akij=ajik=ajki

(i≠j, 或j≠k, 或k≠i;i,j,k=1,2,…,n)

定義6設(shè)有2度n維體陣:

圖1 2度n維體陣與x2的乘積法則

如圖2所示。

圖2 3度n維體陣與x3的乘積法則

3度n維體陣與向量x的冪的乘積方法為：

簡(jiǎn)單地說(shuō)，分三步進(jìn)行，每次都是將體陣中的每一列和與該列平行的x1,x2,…,xn作內(nèi)積，最終形成一個(gè)數(shù)字。由于體陣的對(duì)稱(chēng)性，這種做法可以按可能的任意一種順序進(jìn)行，其結(jié)果都是一樣的。

3 三次型的導(dǎo)數(shù)

根據(jù)上述定義，從導(dǎo)數(shù)的定義出發(fā)，不難推出以下關(guān)于三次型的各階導(dǎo)數(shù)的結(jié)論。

對(duì)于上面的結(jié)論，因?yàn)槎挤浅：?jiǎn)單，所以這里只對(duì)3)中的?進(jìn)行證明。

證明：由于函數(shù)的表達(dá)式是由n3個(gè)項(xiàng)的和構(gòu)成的，所以我們只需證明任意一項(xiàng)的三階導(dǎo)數(shù)都等于其系數(shù)即可。設(shè)任意一項(xiàng)為akijxkxixj, 分3種情形討論。

綜合以上，可知命題3)的?是成立的。

當(dāng)我們給出上述表示方法，并且有了相應(yīng)函數(shù)的導(dǎo)數(shù)記法，原來(lái)已有的多元函數(shù)的一些表示方法可以寫(xiě)成和人們熟悉的一元函數(shù)一樣。

例3 此外還可以定義三次多項(xiàng)式函數(shù)：

其一階導(dǎo)數(shù)為：

二階導(dǎo)數(shù)為:

三階導(dǎo)數(shù)為:

例4 三階偏導(dǎo)數(shù)連續(xù)的多元函數(shù)的泰勒公式可以寫(xiě)為：

f(x)=f(x0)+f(x0)(x-x0)+

4 結(jié) 論

提出了體陣的概念，并通過(guò)體陣的定義推導(dǎo)出體陣的導(dǎo)數(shù)的形式和導(dǎo)數(shù)，給出了一些數(shù)學(xué)結(jié)果及其記法，給出了幾個(gè)體陣的運(yùn)算方法，使得復(fù)雜問(wèn)題得以簡(jiǎn)化規(guī)范，并且將多元函數(shù)的符號(hào)記法及運(yùn)算結(jié)果與一元函數(shù)的符號(hào)記法及運(yùn)算結(jié)果進(jìn)行了統(tǒng)一，使得多元函數(shù)的表示更加簡(jiǎn)潔、更易于理解，該結(jié)果有一定的實(shí)用價(jià)值。

另外，當(dāng)筆者給出三次型之后，很可能有類(lèi)似于二次型的應(yīng)用結(jié)果。例如是否存在三次型的特征值和特征向量，其意義如何？再如，動(dòng)態(tài)系統(tǒng)的穩(wěn)定性一直是理論界和工程界關(guān)注的熱點(diǎn)問(wèn)題之一，因?yàn)椴环€(wěn)定的系統(tǒng)無(wú)法付諸實(shí)際應(yīng)用，而Lyapunov方法是解決這類(lèi)問(wèn)題的有力工具。

對(duì)于線性系統(tǒng)，判斷平衡點(diǎn)穩(wěn)定性的基本思想為：首先構(gòu)造一個(gè)具有二次型的Lyapunov函數(shù)，然后求解Lyapunov方程，就可給出所需的結(jié)果。

對(duì)于多元多項(xiàng)式非線性系統(tǒng)，用類(lèi)似的思想也可以對(duì)平衡點(diǎn)的穩(wěn)定性進(jìn)行判斷，只不過(guò)需要把二次型的Lyapunov函數(shù)推廣到三次正定型，這方面的潛在應(yīng)用還有待于進(jìn)一步研究。

參考文獻(xiàn)：

[1] 李俊民，邢科義，萬(wàn)百五.連續(xù)時(shí)間非齊次雙線性二次型最優(yōu)控制迭代算法[J].西安電子科技大學(xué)學(xué)報(bào)，2000，27(5):630-634.

Li Junmin, Xing Keyi, Wan Baiwu. An iterative algorithm of optimal control for the non-homogenous bilinear quadratic problem[J]. Journal of Xidian University,2000,27(5):630-634.

[2] 尹怡欣，孫一康，舒迪前.具有二次型性能指標(biāo)的簡(jiǎn)單自適應(yīng)控制算法及其應(yīng)用[J].控制與決策，2000，15(2):236-238.

Yin Yixin, Sun Yikang , Shu Diqian. Simple adaptive control algorithm with quadratic performance and its application[J]. Control and Decision, 2000,15(2)：236-238.

[3] 王軍,王雁,王瑞祥,等.采用逆系統(tǒng)設(shè)計(jì)和二次型最優(yōu)控制提高變風(fēng)量箱的性能[J].上海交通大學(xué)學(xué)報(bào)，2006，40(8):1431-1436.

Wang Jun, Wang Yan, Wang Ruixiang,et al.The performance improvement of variable air volume box through inverse system design and quadratic optimal control[J].Journal of Shanghai Jiaotong University,2006,40(8)：1431-1436.

[4] 劉丁, 趙豆, 楊延西.基于快速立體匹配的基本矩陣估計(jì)方法[J].西安理工大學(xué)學(xué)報(bào)，2009，25(1):1-7.

Liu Ding, Zhao Dou, Yang Yanxi. Fundamental matrix estimation method based on fast stereo matching[J].Journal of Xi’an University of Technology, 2009,25(1):1-7.

[5] 陳雪波，徐望寶，李小華,等.非線性系統(tǒng)零解穩(wěn)定性判定的廣義二次型方法[J].控制與決策，2007，22(1)：81-84.

Chen Xuebo, Xu Wangbao, Li Xiaohua,et al.General ized quadratic form method for zero2 solution stability of nonlinear systems[J].Control and Decision,2007，22(1)：81-84.

[6] 劉忠漢, 吳宏鑫.類(lèi)二次型最優(yōu)控制方法在自主對(duì)接控制中的應(yīng)用研究[J].宇航學(xué)報(bào)，2002，32(2)：44-49.

Liu Zhonghan, Wu Hongxin. Application of a similar quadratic ptimization control method in the autonomous docking control[J].Journal of Astronautics, 2002, 32(2)：44-49.

[7] 陳浩，華燈鑫，張毅坤,等.一種基于變換矩陣的激光雷達(dá)數(shù)據(jù)無(wú)量綱統(tǒng)一化轉(zhuǎn)換方法[J].西安理工大學(xué)學(xué)報(bào)，2014，30(1)：1-7.

Chen Hao,Hua Dengxin,Zhang Yikun, et al. A conversion method for radar data non-dimensionalization and standardization based on transform matrix [J]. Journal o f Xi’an University of Technology, 2014 30(1):1-7.

[8] 邵克勇，蔣銳，李文成,等.基于二次型性能指標(biāo)的不確定混沌系統(tǒng)最優(yōu)控制器設(shè)計(jì)[J].化工自動(dòng)化及儀表，2013，40(11)：1347-1350.

Shao Keyong，Jiang Rui，Li Wencheng,et al. Design of optimal controller for uncertain chaotic systems based on quadratic performance index[J].Control and Instruments In Chemical Industry, 2013：40(11), 1347-1350.

[9] 王雙進(jìn)，馮磊,黃賢軍,等.用二次型的正定性判斷晶體相穩(wěn)定性[J].河北工業(yè)大學(xué)學(xué)報(bào)，2007，36(1)：30-34.

Wang Shangjin，F(xiàn)eng Lei，Huang Xianjun，et al.The judgement for crystal stability of phase in positive definite quadratic form[J].Journal of Hebei University of Technology, 2007, 36(1)：30-34.

[10] 宋巨龍，王香柯，馮曉慧.最優(yōu)化方法[M].西安:西安電子科技大學(xué)出版社，2012.

[11] 秦國(guó)強(qiáng).多元函數(shù)的泰勒公式及其應(yīng)用[J].呂梁教育學(xué)院學(xué)報(bào)，2013，30(2)：103-105.

Qin Guoqiang.Taylor rule of multivariate functions with it's application[J].Journal of Lvliang Education Institute,2013,30(2):103-105.

[12] 莫國(guó)瑞，劉開(kāi)第. 函數(shù)逼近論方法[M].北京:科學(xué)出版社，2003.