多项式稳定性的一类新判据

给出了多项式没有一个根位于右半平面的必要条件(定理1),并由此引出“判定系数”,“子判定系数”等概念,获得多项式稳定性的“逐级判定法”。 又利用判定系数给出了几个实用的充分条件(定理2—定理5)。其中最有意义的是证明了 定理4 多项式f(x)=a_0 a_1x … a_nx~n(a_i>0)稳定的一个充分条件是 a_(i-1)a(i 2)≤0.4655a_ia(i 1) (i=1,2,…,n-2) 多项式稳定性的这类判据为判定或设计稳定的线性系统提供了方便。     

Extended graphical representation of polynomials with applications to cybernetics

In this paper,the polynomial of a complex variable s(≡x iy)with realcoefficients K=a_0s~n a_1s~(n-1) …… a_(n-1)s a_n.is graphically represented by three plane curves which are the projections of aspace curve on three coordinate planes of the coordinate system(x,iy.K)inwhich K is confined to be real.The projection on(x,iy)plane is just the rootlocus of the polynomial with K as a real parameter.It is remarkable thatthe equation of the root-locus is m-th degree in y~2,whether n=2m 1 orn=2m 2.In addition to the real curve K.=f(x)in the figure(K,x)thereexists another curve K.which is plotted by the real parts of all complexroots against K.The(K,x)curve is particularly important to determine theabsolute as well as the relative stable interval of K for linear systems.For cybernetics,the(K,iy)curve can be used to show the relation betweenthe nature frequency ω and the gain K.Such three figures are useful forstudying the theory of equation and cybernetics.     

关于均匀变形问题的一种解答方法——讨论一个弹性力学习题的解

徐芝纶编《弹性力学》上册(1979)有这样一个习题(见该书习题 S-2):设某一物体发生如下的位移u=a_0+a_1x+a_2y+a_3zv=b_0+b_1x+b_2y+b_3zw=c_0+c_1x+c_2y+c_3z试证明:各个形变分量在物体内为常量(即所谓均匀变形);在变形以后,物体内的平面保持为平面,直线保持为直线,平行面保持平行,平行线保持平行,正平行 ...     

Singular perturbation of nonlinear vector boundary value problem

In this paper we study the perturbed boundary value problem of the form dx/dt=f(x,y,t;ε), εdy/dt=g(x,y,t;ε), a_1(ε)x(0,ε)+a_2(ε)y(0,ε)=a(ε) b_1(ε)x(1,ε)+εb_2(ε)y(1,ε)=β(ε)in whichx,f,β∈E~m, y,g,a∈E~n, 0<ε(?)1and a_1(ε), a_2(ε), b_2(ε)and b_2(ε) are matrices of the appropriate size. Under the condition that g_y(t) is nonsingular and other suitable restrictions, the existence of the solution is proved, the asymptotic expansion of solution of order n is constructed, and the remainder term is estimated.     

EPSILON-ALGORITHM AND ETA-ALGORITHM OF GENERALIZED INVERSE FUNCTION-VALUED PAD APPROXIMANTS USING FOR SOLUTION OF INTEGRAL EQUATIONS

IntroductionLetf(x ,λ)beagivenpowerserieswithfunction_valuedcoefficients,i.e.,f(x ,λ) =c0 (x) +c1(x)λ+c2 (x)λ2 +… +cn(x)λn+… ,( 1 )wherecj(x)isarealorcomplexfunctionwithregardasxwhichisdefinedon (a ,b) .Supposethatf(x,λ)isanalyticasafunctionofλattheoriginλ=0 .ItisknownthatPad啨ａｐｐｒｏｘｉｍａｎｔｍｅｔｈｏｄｓａｒｅｕｓｅｄｆｏｒｏｂｔａｉｎｉｎｇｔｈｅｓｏｌｕｔｉｏｎｏｆａｗｋｗａｒｄｉｎｔｅｇｒａｌｅｑｕａｔｉｏｎｓ,andespeciallythosewhichpossessageneratingfunctionoftheform ( …     

关于均匀变形问题的一种解答方法——讨论一个弹性力学习题的解

徐芝纶编《弹性力学》上册(1979)有这样一个习题(见该书习题 S—2):设某一物体发生如下的位移u=a_0 a_1x a_2y a_3zv=b_0 b_1x b_2y b_3zw=c_0 c_1x c_2y c_3z试证明:各个形变分量在物体内为常量(即所谓均匀变形);在变形以后,物体内的平面保持为平面,直线保持为直线,平行面保持平行,平行线保持平行,正平行     

关于n跨连拱体系的几何构造分析

平面体系的几何构造分析是结构力学课程的重点和难点之一.本文利用零载法对n跨连拱体系的几何构造分析做一详细讨论,得出了高跨比不相等且跨度也不相等的任意n跨连拱体系几何构造分析的特征表达式fn-fn-1+…+(-1)n-1·f1.当n跨连拱体系的特征表达式fn-fn-1+…+(-1)n-1·f1≠0时,体系是几何不变的;当n跨连拱体系的特征表达式fn-fn-1+…+(-1)n-1·f1=0时,体系是几何可变的.     

Global linear and quadratic one-step smoothing newton method for vertical linear complementarity problems

IntroductionLetN ∈Rm0 ×n(m0 ≥n)beaverticalblockmatrixoftype(m1,… ,mn) ,q∈Rm0 beaconstantvectorpartitionedconformablywithN ,thatareN =N1Nn,　q =q1qn,　m0 = ni=1mi,whereNi ∈Rmi×n,qi ∈Rmi,i ∈I =1 ,… ,n .WeconsidertheverticallinearomplementarityproblemVLCP(N ,q)associatedwith (N ,q)offindingavectorx∈Rnsuchthatx≥ 0 ,　si(x) :=Nix+qi ≥ 0 ,　xi∏mij=1sij(x) =0　　 (i∈I) ,wherexi,sij(x)denotetheithcomponentofxandthejthcomponentofsi(x) ,respectively ;∏ mj=1ajdenotesa1…am.Iti…     

Multiple reciprocity method with two series of sequences of high-order fundamental solution for thin plate bending

IntroductionIt’swell_knownthatthecomplicatedfundamentalsolution[1,2 ]forHelmholtzequationΔu(x) +k2 u(x) =0　　(x∈Ω:boundedopenregioninR2 )isu (x,y) =-iH(2 )0 (k x-y ) 4,thusit’snotconvenientfornumericalcomputation .IfapplyingthesimplefundamentalsolutionofLaplaceequationu 0 (x ,y) =-ln｜x-y｜(2π) ,theexpressionforthesolutionofequationintheclosedregion Ωisc(y)u(y) + ∫Γu(x) u 0 (x,y) nx -u 0 (x ,y) u(x) n dsx =-k2∫Ωu(x)u 0 (x,y)dΩx.Astherightsideappearstheregionalintegrationinclu…     

一个新的运动稳定性判据

1.引言运动稳定性的判据有很多种,例如Routh-Hurwitz准则、Nyquist准则、谢绪恺判据等。不过,这些判据都需要知道线性系统的特征多项式,而由线性方程组的系数矩阵计算特征多项式需要大量的计算工作。李森林给出了一个利用系数矩阵判别系统稳定的充分条件,使用这个条件需要计算n个n-1阶行列式和一个n阶行列式。由于这个条件仅仅是充分的,因此,当计算结果满足这个条件时是能说明问题的;如果不满足这一条件,计算结果便毫无用处。尽管如此,在参考文献[4]中仍认为此法“比用其他方法来判定要简单得多(尤其当n很大时),所以在实际工作中可     

Analysis of chebyshev pseudospectral method for multi-dimensional generalized srlw equations

IntroductionAsymmetricregularizedlongwaveequation (SRLWE) 2 x2 -1 u t = x ρ+ 12 u2 ,　 ρ t+ u x=0 (1 )hasbeeninvestigatedinRef.[1 ] .Thesystem (1 )ofequationsisshowntodescribeweaklynonlinearion_acousticwaveandspace_chargewaves.Thehuperbolicsecantsquaredsolitarywaves ,thefourconservationlaws,andsomenumericalresultshavebeenobtainedinRef.[1 ] .Obviously ,eliminatingρin (1 ) ,weobtainaclassofregularizedlongwaveequationutt-uxx+ 12 u2xt-uxxtt =0 . (2 )TheSRLWequationisexplicitlysymmetric…     

高压下熔化温度T_m和格临爱森系数的关系

本文导出了两个简化公式,其形式为: =2/3+（n-2/3）0V Tm=T0m（/n）23exp[2（n-2/3）0(1/0m-1/)] 所推导出的表示式与Slater公式,Dugdale-MacDonald公式和自由体积公式在很宽的压力范围都符合得很好。Tm表示式与Romain、陈栋泉和林绍明以及张万箱和张世泽得出的结果一致。     

源像法在平板弯曲边界元法中的应用

本文应用源像法推导了平板弯曲问题基本解的三种新的形式,它们为: W~*=1/(16πD)(γ_3~2ln((γ_4~2)/(γ_3~2))-γ_1~2ln((γ_2~2)/(γ_1~2))(1) W~*=1/(16πD)(γ_1~2lnγ_1~2-γ_2~2lnγ_2~2-γ_3~2lnγ_3~2 γ_4~2lnγ_4~2)(2) (3)式中γ_1~2=(x-ξ)~2 (y-η)~2 γ_2~2=(x ξ)~2 (y-η)~2 γ_3~2=(x-ξ)~2 (y η)~2 γ_4~2=(x ξ)~2 (y η)~2 γ_i~2=(x-γcos(α (2π(i-1))/n))~2 (y-γsin(α十(2π(i-1)/n))~2 γ_i~2=(x-γcos((2πi)/n-α))~2 (y-γsin((2πi)/n-α))~2n为整数式(1)相当于1/4无限板的解,其中一边固定,一边铰接; 式(2)也相当于1/4无限板的解,其中两边都为铰接; 式(3)相当于扇形无限板的解,其中两个直线边都为铰接。应用上述形式基本解很多情况下板的部分边界条件将自动满足。它将大大降低边界元方程式数目,缩短计算时间,减少数值误差。文中附有计算例子,结果与解析解极为符合。     

应用弹性动力学广义变分原理的一步时间元方法

本文应用笔者在文献[1]中给出的给定时端动量情形的弹性动力学广义变分原理在空间域离散并记入阻尼力虚功的矩阵形式δH[x]=intergral from n=t_1 to t_2(δx~Tkx-δ?M?+δx~TC?-δx~Tf)dt+δx_(t_2)~TM?-δx_(t_1)~TM?=0     

自伴微分方程特征值问题的子空间迭代法及其在结构振动分析中的应用

1.自伴微分方程特征值问题考虑特征值问题L_φ=λ_(ρφ) 在[a,b]上(1)其中λ为参数,ρ(x)为闭区间[a,b]上的连续函数,并且ρ(x)>0,L为n阶自伴微分算子其中P_j∈C~(n-j)在a≤x≤b上(j=0,1,…,n),P_j为x的实函数。L的定义域为满足与(1)相应的齐次边界条件并属于C~n的所有函数的线性集合,记作D,φ∈D。     

Haug-Arora梯度投影法的改进

在Haug E. J. 和Arora J. S. 提出的梯度投影法中,所用的基本迭代表达式为{x}~(k+1)={x}~k-η~(k+1){δx_1}~(k+1)+{δx_2}~(k+1) (1)式中{δx_1}~(k+1)为第k+1次迭代使目标函数减小的向量,{δx_2}~(k+1)为修正第k次迭代所引起的约束误差的向量,η~(k+1)为第k+1次迭代的步长系数,它由下式确定     

Kelvin问题的一种解法

1.Kelvin问题所对应的定解问题在无界弹性体内,把集中力的作用点选为坐标系Ox_1x_2x_3的原点.设集中力的大小为P(常数),其方向沿单位矢量n,则此集中力可表示为Pδ(x)n,其中δ(x)为Dirac的δ函数,x为空间中一点的位置矢.弹性体内各点的位移矢量u是点的坐标的函数,表示为u=u(x).位移u在无界域中所应满足静     

Existence of solutions for nonlinear elliptic boundary value problems

IntroductionInRef.[1 ]KannanandLockergivetheexistenceofatleastonesolutionofTy-h(t,y ,… ,y(n- 1) )y=f(t,y ,… ,y(n- 1) )　　 (a相似文献   

梯度投影单纯形法及其应用

对于下列问题min f(x) s·tg_1(x)=sum from j=1 to n g_(ij)x_j≤b_i(i=1,2,…M) (1) 其中,x∈E~n,g_(ij)为常数,J·B·Rosen给出了梯度投影解法: x~((k+1))=x~((k))+a~((k))s~((k)) (2) s~((k))=-(I-N_Q·(N_Q~TN_Q)~(-1)·N_Q~T·▽f(x~((k))) (3) N_Q=[▽g_1…………▽g_Q] (4) ▽g_i=[g_(i1)…………g_(in)]~T (5) a~((k))=(?)f(x~((k))+as~((k)) (6)     

Global structure of a kind of predator-prey system

In this paper we prove a theorem, theorem 2, on nonexistence of closed trajectory for a general predator-prey system. Then, using this theorem and another theorem on existence and uniqueness of limit cycle for predator-prey system, we complete the investigation of a concrete model of predator-prey system (?)=γx(1-x/K)-yx~n/(a+x~2) (?)=y(μx~n/(a+x~2)-D)(n=1,2)under the conditions of all kinds of parameters.