维普资讯 http://www.cqvip.com 2002年6月 重庆大学学报(自然科学版) Vo1.25 No.6 第25卷第6期 Journal of Chon ̄University(Natural Science FAition) Jun.2OO2 文章编号:1000—582X(2002)06—0054—02 关于系统 = ( )一 ( ), =一g( )的极限环之唯一性。 张谋 (重庆大学数理学院,重庆4OOO44) 摘要:对非线性系统(1)用文献[1]的思路和方法,给出了一个极限环存在的唯一性定理,去掉了文 献[1]定理1的条件1),使该定理的条件更容易验证;然后将该定理应用于以生化反应为背景的一类平面 四次微分系统,得到了该系统有唯一单重稳定的极限环的充分条件;应用于文献[3]的三次系统,所得结 论与文献[3]相同。 关键词:非线性系统;极限环;唯一性;存在性 中图分类号:0175 文献标识码:A 1准备工作 z< ;kl( )≥k2(z), ≤z< ; 2) x)/g(x)在( ,I"2)内单调上升,且0<△ 。; 考虑系统 或 )/g( )在(r。,瓦)内单调下降,且z。 z 。 戈= (Y)一F( ), =一g( ),F( ):I/’(s)ds 其中瓦= 。(三), = 2( ),同时瓦<0, > l;则系 J 0 ),g( )∈cl ),一∞≤rl<r2≤+∞ 统(1)至多有一个极限环,若存在是单重的稳定环。 g(0)=g(A)=0,xg( )>0( ≠0,△)(1) 2 主要结果及证明 (Y)>0(Y≠0),声(Y)∈c ,非减。 假设0≤△<r2,存在 0, l:0< 0< 】<r2, 定理:设系统(1)满足假设(2),且满足 g( )/f( )在(r。,0)和( ,r )内单减(或单增) 使得F(0)=F( )=0, 当 < 0时 )<0,当 > 0时 )>0 则系统(1)至多存在一个极限环,若存在是单重稳定 的极限环 J。 (2) 作变量代换z=F( ),令z。=F( 。),z =F(r + 证明:只须证明文献[1]定理1中条件1)成立,令 0),:2=F(r2—0), =rain(zl,z2)并用 = 1(z), H(z)=kl(z)一k2(z),0<z< ,因为g(xl(z))I :0= 0< < l和 = 2( ), 0< < 2表示 =F( ) g(x。(0))=0,g(x (z))I删=g(x (0))=g(x。)>0, 在(r , 。)和( 。,r )上的反函数,则系统(1)在 < (z))I::。:/’( 。( 。(0))=f(0)<0√^( (:))l :。= f( (0))=f( 。)>0,所以日( )I :。=[k。( )一 。和 > 。内分别等价于z—Y平面上z>z。内的下 列两个方程 [z一 (Y)]dy=dz, k2(:)]圳=一g( ) ( )<0,由假设 【臀 】_ 其中 k (:)=g( (:))/f( (z)),i:1,2 <0 ’0)u x,,r2), 笔者利用文献[1]的定理1得到下面定理,该定理 所以 改善了文献[1]定理1的条件,使定理条件更容易验 d日( ) f( (z))g ( 。(:))~g( 。(z))fx ( ) 证,为方便先将文献[1]的定理1叙述如下。 d 一 厂 。( )) 定理1[11:设系统(1)满足假设(2),另外 d 。/’( (z))g ( (z))一g( (z))/( :(z)) …1)存在三:0<乏<z ,使kl(z)<k2(z),z0< dz— 厂———~‘ ・收稿日期;2002一O1—18 作者简介:张谋(1963一),男,重庆人,重庆大学讲师,硕士。主要从事常微分方程的研究。 维普资讯 http://www.cqvip.com 第25卷第6期 张 谋: 关于系统i=牵(y)一V(x), =一g(x)的极限环之唯一性 55 d d 一 。( ))g ( 。( ))一g(z ( ))/( t(。)) 厂< 1( )) >0 , 2 ={(u, ): >一B, >一A} 容易验证: ( )>0, ( )单增,ug(u)>0;存在 // ̄0-.//q使 0) (11,1)_0,F 11,0)=0又 = 即当z E[0,z ]时H (z)>0,又H(0)<0,所以:或者 (z)在(0,z )内没有零点,由文献[1]引理知此时系 F (u) 2(u+日) +(u+日) 一b g(u)一一(u+B)(u+B+1)[(u+ 一日 ]一 一 ,统(1)无闭轨;或者存在唯一的 使 ( )=0,0< < z ,其中s=日一 0<s<+∞,从而 则表示:当0<z< 时kl(z)<k2(z),当 <z< z 时k (z)≥k (z),另一部分同理证明。 [ ]= [ ] =一 ‰ 3 应用 将该定理应用于以生化反应机理为背景的一类平 面四次微分系统(常见的是三次微分系统): dx= R(s)=2s +2s +2B s +4(B 一b)s + (B 一b)s +2B bs+(s 一2bs +B b) 当s E(B,+∞)U(o,B—u ),c>2b时R(s)>0, 即当u E(一∞,0)U(u , ,c>2b蝴u)/g(H)单 减。这样就证明了结论。应用到文献[3]的三次系统, 可得同样正确的结论。 (3) 参考文献: [1]王现.关于系统 = (y)一F( ), =一g( )的极限环 之唯一性[J].南京大学学报,1990,26(3):363—372. [2]丁大正.Lienard方程极限环的存在性[J].应用数学学报, 1984,7(2):166—174. (1+),)(c一6一y2) 业dt= (1+_y)(6+y2)一ay 其中参数0>0,b至0,c>0,Q={( ,Y): ≥0, Y≥0},得到下面结论。 结论:当c>26>0且 <; 时,式(3)在Q内存在唯一稳定的极限环。 证明:对方程(3)作变换 =口+A,Y=u+日,并用 ( +A)(u+ 同时除所得方程的右端,方程(3)变为: u.[3]周建莹,张锦炎,曾宪武.生化反应中一类非线性方程的 定性分析[J].应用数学学报,1982,5(3):234—240. [4]余澍祥.极限环的存在性定理[J].数学进展,1965,8(2): 187—194. :一g(u),d一g “ ,d£ ,:声(口)一F(u)一 “ (4) [5]梁锦鹏.Lienard方程极限环存在定理[J].高校应用数学 ’其中A:口 /c(1+ ̄/厂 ),B: ̄/厂 ,在 学报,2000,15(2):163一l68. 已知条件下,(A,B)为唯一的正平衡点,g(u)=u(1+ u+B)(u+2日)/(u+B), (口)=av/A(口+A) [6]吴葵光.非线性方程极限环的存在性[J].数学学报,1982, 25(4):456—463. F(u):一u『(u+日) +(1+日)(u+B)一b/B]/(u+B), On the Uniqueness of Limit Cycle of System =|l}( )一F(x),j》=一g(x) zHANG Mou (College of Mathematics and Physics,Chongqing Universit7,Chongqing 4OOO44,China) Abstarct:In this paper,we give the uniqueness theorem of limit cycle to system丘=声(Y)一F( ), =一g( )by using idea and method for paper[1],in which if the hypothesis 1)of the theorem 1 for[1]is omitted;it is easier to prove the conditions ofthe theorem.Then we obtain some sufficient conditions ofone simple stable liitm cycle for a plne4——degree system ain biochemical reactions,he 3一degrtee system method that is applied to apper[3],it has the same evidence as paper[3]. Key words:the limit cycle;nonlinear system;the uniqueness;the existence (责任编辑张苹) 维普资讯 http://www.cqvip.com 2(102年6月 重庆大学学报(自然科学版) Vo1.25 No.6 第25卷第6期 Joumal of Chongqing University(Natural Science Edition) Jun.20()2 文章编号:1000—582X(2002)06—0056—03 基于岭估计的最优预测与经典预测的最优性判别‘ 杨 婷,杨 虎,张洪阳 (重庆大学数理学院,重庆400044) 摘要:针对有偏估计的预测问题,以岭估计为基础,以离差矩阵MDE(Mean Dispersion Error)和广义 风险函数为判别准则,对广义线性回归模型{Y=郎+e,e~N(0,盯 W)}的最优预测量与经典预测量的 最优性判问题进行了讨论。借助矩阵不等式的一些性质,获得在离差矩阵判别准则和风险函数判别准则 下判断两类预测量最优性的充要条件。为进一步研究基于有偏估计关于两类预测量的最优性判别问题 提供了一种方法和思路。 关键词:预测;岭估计;离差矩阵 中图分类号:O212.1 文献标识码:A 1预备知识 预测在实践中有着广泛的应用,在线性模型中也有 着不可忽视的作用。所谓预测,乃是对给定的回归自变 量的值,预测对应的回归因变量所可能取的值,这是回 归分析的重要应用之一 l ]。和估计一样,预测也有点 预测和区间预测之分,这里只针对点预测进行讨论。 (注意:若无特别申明,本文所指的矩阵皆为实矩阵。) 文献[1—3])中都考虑其风险函数的最小值 m inR (P,Y )=m inE[(P—Y ),A(P—Y )] ,,(4) 其中,矩阵A≥0,可以视具体情况而定。在式(3)、式(4) 的基础上,文献[3]得出三类特殊的最优预测量: 1)非齐次最优预测量Jl;1=X + W (Y一邢 ) 2)齐次最优预测量J】; =X + 。W (Y—X/32) 3)齐次无偏最优预测量J】; =X b。。皿+ w,0 (Y—XbGIsE) 我们将历史数据(Y,X)写成如下线性模型 Y= +e,e~N(0,盯 W) (1) 其中Y为t×1观测向量, 为t× 矩阵,且rank(X)= k,W为已知的阶正定矩阵。 其中 。= 即为回归参数本身, =卢[ 】=p(yvp) ’, (令 : + , = (X 一X+ 将n组观测数据(J, ,X )也用矩阵表示为如下 线性模型 Y =X n.1 n 盯 ,', = +e , n 1 ) X Y)即为两阶段估计 ,且 娅= (2) f.n e ~(0, W ),E(ee )=盯 Wo n.n (X W X) X W Y即为广义最小二乘估计。通过计 算, ( , )≤R ( , )≤R ( , ),而且有 R (pl,Y )≤R (p 2,Y )≤R (p 3,Y )恒成立 。 这说明参数估计的改进也可改进模型的预测效果。从 其中,参数 与式(1)中的相同,Y 未知, 已知,并且W+>0,Wo≠0。 , , 我们的目的是为了预测Y ,不妨假设p为Y 的 预测量。为了讨论之便,我们考虑其最简单形式,即 P=Cy+d (3) 1)2)3)也可看出,预测量的估计量与回归系数 的估 计量有很大关系,因此,不妨定义p的最优预测量的一 般形式 预测量P的预测效果与所选择标准有关,一般文献(如 ・收稿日期:2002—01—16 作者简介:杨婷(19r75一),女,四川巴中人,重庆大学硕士研究生。主要从事统计线性模型研究。 维普资讯 http://www.cqvip.com 第25卷第6期 杨 婷等: 基于岭估计的最优预测与经典预测的最优性判别 57 p( )=X p+w 。w (.),一 ) (5) 应矩阵的广义逆矩阵。 另一方面,我们也希望未知观测量Y 的条件期望 E(y /X )=X ,从而不妨定义p的经典预测量的 一引理2的证明请参见文献[6]。 推论 假设E为rl,阶矩阵,,为rl,阶单位矩阵,则 E+E ≥0的充要条件为 般形式 p。(p)=X (6) [(E一,) (EE +I) (E—I)]≤1 证明 因为EE ≥0,I>0,从而EE +I>0, 文献[3—4]主要针对几种特殊的估计量( ,= , = 【 ], , = 6GL咂)x1, p的两类预测 故EE +,可逆。据引理2,原命题成立。证毕。 定理1 R (P ,Y )一 (P2,Y )≥0甘 (p , 量关于风险函数(R (P,Y )=E(p—Y ),A(P— Y ))的最优性进行了讨论。而文献[5]就 的一般线性 估计量对p的两类预测量关于离差矩阵( (p,Y )= E(p—Y )(P—Y ) )的最优性判别条件进行了详细 的讨论,并得出一重要结论(定理6.4)。虽然此结论对 任何线性无偏估计都适合,但有偏估计量(如岭估计 等)并不满足此定理所给条件。因此,有偏估计关于两 类预测量的最优性判别问题仍未解决。笔者主要对常 见的有偏估计——岭估计关于两类预测量的最优性 判别条件进行了讨论。 2 主要结果 下面仍考虑线性模型1)、2),且一切假设条件均 不变。 定义1 设p ,p 为Y 的两个预测量,如果它们 满足 ’ (P ,Y )一尺 (P2,Y )≥0 则我们称预测量p 关于 一优于预测量p ,或简称 为 一准则。 定义2 设p,,p 为Y 的两个预测量,如果它们 满足 M(pl,Y )一M(p2,Y )≥0 则我们称预测量P 关于MDE一优于预测量P,,或简称 为MDE一准则。 引理1 设D为r/,X r/,阶对称矩阵,则 D≥0甘对任意r/,X r/,非负定矩阵C,满足 tr((CD)≥0 此引理的证明可参见文献[4]附录定理A.43。 引理2 设A, 均为m x rl,复矩阵,A =A , =B ,则AB +BA ≥0的充要条件为 l[(A—B) (AA +BB )一(A—B)]≤1 其中, (・)表示相应矩阵的最大特征根,(・)一表示相 Y )一M(p2,Y )≥0。 证明 因为几(P,Y )=E(p—Y ) A(P— Y ),M(p,Y )=E(p—Y )(P—Y ) 从而 尺^(P ,Y )一 (P2,Y )=trA・ [M(p ,Y )一M(p ,Y )] 又因A≥0,所以根据引理1可得 尺』(p1,Y )一尺 (P2,Y )≥0甘 M(pl,Y )一M(p2,Y )≥0 故原命题成立。证毕。此结论说明以上两种判别准则等价。 记模型(1)的岭估计为 =(S+ ,) WI1Y垒Ⅳ: X WI1Y_△D 其 =X W~X,M=S+ ,,D =Ⅳ w一 因M=S+ ,所以Ⅳ S=(S+ ) (S+ 一 ): I一 从而有,一Ⅳ: S=kN-k (7) 易得 的期望、偏崎和方差分别如下: fE( )=Ⅳ w一 :Ⅳ: {Bias( ̄ , )=E( )一 :(Ⅳ: S一,) =一 V: tVar( ̄^): Ⅳ SN-k (8) 由式(7)、式(8)可得 ( , ):E( 一 )( 一 ) : E[( 一E( ))+(E( )一 )]・ [( 一E( ))+(E( )一 )] : Var( ̄ )+[Bias( ̄ , )][Bias( ̄ , )] : ,2 0.2Ⅳ: : + Ⅳ: Ⅳ: = 2Ⅳ (S+ )Ⅳ 0 由式(5)、式(6)可知,基于岭估计的最优预测量与经典 预测量分别为p( )=X + 。w (.),一 ), Po( )=X 现记z= 一 0 w ,易得 维普资讯 http://www.cqvip.com 58 重庆大学学报 (自然科学版) 200 年 M(p(flk),Y )=ZM(flk,p)z +0.2(W 一w 。W—Wo) M(p。( ), ): . M(fl , ). 一盯 . DtWo一 盯 Wt0DPk +盯 W 从而 M(p。( ),Y )一M(p(fl ),Y )= X M(flk,1f)x 一ZM(flk,f1)z 一G2 x Dkwo一 盯 0D'kX +仃 0 W_。Wo:仃 X Jv】 ・ (S+ f;9 ) : W~wo+盯 0W一删 ・ (S+ ) 一 0 W—XNT, ・ (s+ ) W~Wo一(72 D,Wo一 盯 0D t +盯 Wp0W_。Wo:盯 : ・ ( k ̄flf')DtWo+盯 W'oD'k(IS+ ) X*一 盯盯 2l0 盯 0D t・(S+ )D,Wo一(72 DtWo一 盯 矸,,0D'k. +盯 H,,0 W_。Wo:盯 . N-k ・ [(S+ )一 ] +盯 。D ・ [(S+ 舟 )一^ ]^ . +盯 H,,0・ [W~一D t(S+ ) ]Wo:A 盯 [(At+A t)+Wp0(W_。一D t^ D )H ] 其中A = Ni [(S+ )~ ] 。欲证明 上面矩阵为非负定矩阵,只须验证矩阵A +A 的非 负定性。因为 。( ~一D )Wo≥0恒成立,事 实上,由于Wo>0,Dt= : W~,从而 w 0(w一一D t tDt)Wo≥0§ W~一D',NkDt≥0臼 W~一W一删: NkNi, W~≥0臼 W~一W一 XN-k W~≥0臼 W一 — W一 XN-k W一 ≥0臼 S—SNi, S≥0 此式显然成立,因为 >Sj i <S~j S< SS~S=S 因此根据推论,我们可得如下结论: 定理2 已知线性模型(1)、(2),其中模型(1)中参 数 的岭估计为 =(S+ ) w— , 的最优预 测量p( )关于MDE一优于经典预测量p。( )(即 M(p。( ), )一M(p(fi ),j, )≥0的充要条件为: [(At—J) (AtA t+ )I1(At—J)]≤1 由定理1知,对任意n阶矩阵A, 的最优预测量 p( )也关于R 一优于经典预测量p。( )。 3结语 目前,有关有偏估计关于模型的预测问题的文献 出现较少。Bibby和Toutenburg在上个世纪末就注意 到:回归参数估计量的改进可以优化线性模型的预测 效果。在此基础上,Rao将参数估计从无偏推广到有 偏、从特殊(如最小二乘估计)推广到一般,建立了一系 列判别最优预测量与经典预测量最优性的判别条件, 但对某一具体的有偏估计(如岭估计)却不能得到验 证。这里讨论岭估计关于两类预测量最优性的判别问 题,在离差矩阵判别准则下给出一充要条件,从而为研 究有偏估计的预测问题提供了一种方法。 但是关于有偏估计的预测问题,还有几点值得探讨: 第一:除了MDE一准则或R 一准则,能否探寻到 其它或许更加有效的判别准则。 第二:除了讨论岭估计关于两类预测量的最优性 判别问题,也可对其它有偏估计如James—Stein估计、主 成分估计等关于此类问题进行探讨。能否找到一种综 合有偏估计类(如包括以上一些常见的有偏估计)对此 类问题进行讨论,结果又将如何。 第三:由于岭估计中包含一个未知参数,针对关 于两类预测量的最优性的判别问题,如何确定该未知 参数的显示表达式,也值得进一步研究。 参考文献: [1]WEISBERG S.应用线性回归[M].北京:中国统计出版社, 1998. [2] 王松桂,陈敏,陈立萍.线性统计模型[M].北京:高等教 育出版社,1999. [3] BIBBY J,TOffFENBURG H,P ctiOil and Improved Estimauon in Linear ModelslMJ.New York:W'dey,1977. [4] 王松桂.线性模型的理论及其应用[M].安徽教育出版 社,1987. [5] RAO C R,TOUTENBURG H.LinP ̄tr Models—Least Sauares and A1 mati es【M].Ne'n York:Springer—Ve ,1995. [6] JERZY K.BAKSAI.^RY,BERNHARD SCHIPP,et a1.Some further results Oil Hermitian—martrix ineqalities[J].Linear Algebra and Its Applications,1992,160:119—129. (下转第70页) 维普资讯 http://www.cqvip.com 70 学出版社,1998. 重庆大学学报(自然科学版) 华大学出版社,2000. 20o2年 [4]房有定,卢安国.80X86系列保护模式高级程序设计[M]. 西安:西安交通大学出版社,1996. [6]易钟芳.80X86微型计算机原理及应用[M].北京:电子工 业出版社,1995. [5]DAVID A.Solomon WINDOWS NT技术内幕[M].北京:清 Memory Accessing Mechanism of X86 Computer System In Protected Mode Z}_lAN Hu|一|ai.ZENG Yi (College of omputCer,Chongqing University,Chongqing 40OO44,China) Abstract:Protecting illode is on importmant working method of X86一based computer system,which provides the hardware support for multi—task and multi—process operation system on X86 platform.During the low—level development of X86一based computer system,it is very important to understand the working method of memory that is a main part of computer system.This paper discusses the data structure and hardware component of X86一based computer system which support the two mechanism of accessing memory:segmenting and p ng,as well as the process of converting logical address to linear address and mapping virtual address to physical address.Througn seg ̄nenting and paging mechanism,we call realize the independent address space for every process in multi—process system,as well as the measure for protectig tnhe system memory which can prevent the system breakup caused by illegal memory access. Key words:protected——mode;descriptor;describe table;page table;page dircteory (责任编辑 吕赛英) 、驴 抄 、驴 (上接第58页) Comparison of Superiority of Optimal and Classical Predictions with Respect to the Ridge Estimation YANG Ting,YANG Hu,ZHANG Hong—yang (Math.&Phys.Institute of Chongqing Universiy,Chongqitng 400044,China) Abstract:Considering the generalized linear regression model = +e,e~N(O, W)and its prediction problem of biased estimation,hhis paper discusses・its superioriy of the opttimal and clssiacal predictors based on the rid estiatmion by criteria of mean dispersion error matrix and generalized risk function.A necessary and sufficient condition of comparison of its superiority are gi,,cn under the condition of above criteria by some properties of inequalitis of materixA_n alteITIative method is .proposed for the further research of superioriy oft two predictors based on the bised asteimationKey words:predictor;ridge estimator;mean dispersion error matrix . (责任编辑张苹)