GMRES方法的收敛率

1 引 言 GMRES方法是目前求解大型稀疏非对称线性方程组 Ax=b,A∈R~(n×n);x,b∈R~n (1)最为流行的方法之一.设x~((0))是(1)解的初始估计,r~((0))=b-Ax~((0))是初始残量,K_k=span{r~((0)),Ar~((0)),…A~(k-1)r~((0))}为由r~((0))和A产生的Krylov子空间.GMRES方法的第k步     

具有共振的边值问题解的存在性

利用Mawhin的重合度理论,研究具有共振的n-阶m-点边值问题x~((n))(t)=f(t,x(t),x′(t),…,x~((n-1))(t)),t∈(0,1)x(0)=x(η),x′(0)=x″(0)=…=x~((n-2))(0)=0,x~((n-1))(1)=α_ix~((n-1))(ξ_i)解的存在性,其中n≥2,m≥3,f:[0,1]×R~n→R将有界集映为有界集,且当x(t)∈C~(n-1)[0,1]时,f(t,x(t),x′(t),…,x~((n-1))(t))∈L~1[0,1],0<ξ_1<ξ_2<…<ξ_(m-2)<1,0<η<1,α_i∈R.在这里并不要求f具有连续性.     

计算方程重根及其重数的算法

给出一种计算方程重根及重数的迭代算法,分别具有平方收敛和线性收敛.(i)迭代:x_(n+1)=x_n-f x_n (f'(x_n))/((f'(x_n))~2-(f(x_n)f~n(x_n)),m_n=((f'(x_n)))~2/((f'(x_n))~2-f(xn_)f″(x_n)),n=0,1,2,…,重数m≈mn;(ii)加速迭代:x_(n+1)=x_n-(f~((m-1))(x_n))/(f(~m)(x_n)).     

同时求解多项式所有二次因子的迭代法

1. 引言 如所周知,从初始近似值x~((0))出发,使用Newton法 x~((m+1))=x~((m))-f(x~((m)))/f′(x~((x))) (m=0,1,…)求函数f(x)的根时,在单根附近二阶收敛。 当f(x)是首项系数为1的N次多项式     

线段连续自映射非游荡集的拓扑结构

<正> 令X为拓扑空间,f:X→X为连续映射.f的不动点集F(f),周期点集P(f),周期点的周期,以及非游荡点集Ω(f)定义如常(例如,参见文献[1]).令x∈X,集合{f~n(x):n=0,1,2,….}称为x的轨迹,并记作O(x,f);当x为f的周期点时,O(x,f)称为x的周期轨迹.记Ω(f)为具有无限轨迹的非游荡点的集合.y∈X称为x∈X     

多元线性模型中一个二次估计的最优性(Ⅰ)

考虑线性模型设ε′=(ε_((1)),…,ε_((n))),对ε_((1)),…ε_((n))独立,Eε_((i))ε′_((i))=Σ,E(ε_((i))ε′_((i))ε_((i))ε′((i)))=(i=1,…,n)的情形本文求出了Σ的(一定意义下的)最小二乘估计Σ~*,并给出了tr(CΣ~*)是tr(CΣ)的一致(对Σ≥0,Ψ)最小方差不变二次无偏估计的充要条件,这里C是对称矩阵。对Covε=GΣ,Y服从准正态分布的情形也做了相应的讨论,这里G是已知n阶非零的非负定矩阵,Σ是未知的p阶非负定矩阵。     

等距节点上(0,3)插值五次样条

§1.引言 设f(x)是定义在[0,1]上的连续函数,n是自然数。记h=1/n, f_v~((r))=f~((r))(vh),v=0,1,…,n;r=0,1,…,5, f_(v 1/2)~((r))=f~((r))((v 1/2)h),v=0,1,…,n-1;r=0,1,…,5, ω_r(j)=max |f~((r))(x_1)-f~((r))(x_2)|,r=0,1,…,6. |x_1-x_2|≤h 0≤x_1,x_2≤1又设s(x)是[0,1]上满足(i)s(x)∈C~3[0,1],(ii)在[vh,(v 1)h]上s(x)∈∏_5,v=0,1,…,n-1的五次样条.它们的全体记为?_(n5)~((3)) .     

由谱数据数值稳定地构造实对称带状矩阵

§1.引言 设r,n是正整数并且0r有a_(ij)=0.     

BIVARIATE RATIONAL INTERPOLANTS WITHRECTANGLE-HOLE-STRUCTURE

1.IntroductionGivenasetofdistinctrealpoints{xi,i~0,1,2,',n:xiER}andasetofcomplexvectordata{d'),i=0,1,2,',n:n)ECd},Graves-Momsshowed[5]thatthevectorvaluedThieletypecontinuedfractioncanservetointerpolatethegivenvectors.TheconstructionprocessiscloselyralatedtotheadoptionoftheSamelsoninverseforvectorswhere7denotesthecomplexconjugateofvector6.ItwasprovedthatS(x)isavectorvaluedrationalfunctionwithnumeratorbeingad-dimensionalpolynomialofdegreenanddenominatorbeingapolynomialofdegree2[n/2],here…     

具有随机窗宽核估计的一致强相合性

设 X_1,…,X_m i.i.d.是取值于 R~n 中的随机向量,X_1 有概率密度 f(x),取正随机变量 H_m(x,ω)=H_m(x,X_2(ω),…,(ω))为随机窗宽,f(x)的核估计与最近邻估计分别如下:f_m(x)=(mH_m~n(x,ω))~(-1)sum from i=1 to m K((X_i-x)/H_m(x,w))f_m(x)=(ma_m~n(x,w))~(-1) sum from i=1 to m K((X_i-x)/a_m(x,w)),m≥1,x∈R~n.假定 K 为 R~n 中有界变差函数,当 f(x)与 K(x)的条件比[1]弱时,我们讨论了 f_m(x)与 f_m(x)的一致强相合性。本文所得随机窗宽的结果与[1]中常数窗宽的结果相同,这些结果也比[2]和[5]中的要好。     

Liénard方程极限环的唯一性定理

本文讨论Liènard方程 x=y-F(x) y=-g(x) (1)极限环唯一性的条件。其中F(x)=intergral from n=0 to x(f(ζ)d(ζ),以下恒假定f(x),g(x)∈C~0(x_(02),x_(01),x_(02)<00,x≠0 (H)其中x_(02),x_(01)可以是∞。令z=G(x)=intergral from n=0 g((ζ)dζ,记x_1=G~(-1)_1(z)、x_2=G~(-1)_1(z)分别是z=G(x)在(0,x_(01))、(x_(02),0)上的反函数,F_1(z)=F(G~(-1)_1(z)),作φ变换,则·dz/dy=f_t(z)-y,0≤z相似文献   

Quasi-periodic solutions for the general semilinear Duffing equations with asymmetric nonlinearity and oscillating potential

In this paper,we prove the existence of quasi-periodic solutions and the boundedness of all the solutions of the general semilinear quasi-periodic differential equation x′′+ax~+-bx~-=G_x(x,t)+f (t),where x~+=max{x,0},x~-=max{-x,0},a and b are two different positive constants,f(t) is C~(39) smooth in t,G(x,t)is C~(35) smooth in x and t,f (t) and G(x,t) are quasi-periodic in t with the Diophantine frequency ω=(ω_1,ω_2),and D_x~iD_t~jG(x,t) is bounded for 0≤i+j≤35.     

对一道习题的思考

从相关习题出发，借助夹逼定理可证明：lim n→∞（b1a^n1＋b2a^n2＋…＋bma6n m）1/n=max{a1,a2,…,am}；设函数φ（x）,f（x）在[a,b]上都是正连续函数，则有lim n→∞{∫^b aφ（x）[f（x）]^n dx}^1/n=max a≤x≤b{f（x）}     

更新序列对于有限Markov链的嵌入

§1 引言数列 f=f~(1),f~(2),…,f~(n),…}称为,一序列,如果f~(i)≥0(i≥1);sum from t=1 to ∞ f~(i)≤1 (1)由产生的更新序列 u-{u_0;u_1,u_2,…,u_n,…}依下式定义(2)更新序列与马氏链关系密切。设 X(n)是离散参数马氏链,其(一步)转移矩阵为P=(P_(ij))_(i,j∈E),(E 为可列集) (3)又记 n 步转移矩阵为 P~((n))=(P_(ij)~((n)))_(i,j∈E),则P~((0))=(单位矩阵),P~((1))=P,P~((n))=P~n (4)这时,对每个 i∈E,数列{P_(i)~((n))}_(n≥0)是更新序列,其所有产生的 f-序列为{f_i~((n))}+_(n≥1):     

等距节点三角插值的Lebesgue常数的渐近展开

设 f(x)是以2π为周期的一个周期函数,我们知道对于[0,2π]上的2n 1个节点 θ_k=k(2π/(2n 1)),k=0,1,2,…,2n,(1)存在唯一的n次三角多项式L_n(f,θ),满足L_n(f,θ_k)=f(θ_k),k=0,1,2,…,2n。这里     

On an Extension of Hadamard Inequalities for Convex Functions(in English)

设f是区间[a,b]上连续的凸函数，我们证明了Hadamard的不等式 $[f(\frac{{a + b}}{2}) \le \frac{1}{{b - a}}\int_a^b {f(x)dx \le \frac{{f(a) + f(b)}}{2}}$ 可以拓广成对[a,b]中任意n+1个点x_0,\cdots,x_n和正数组p_0,\cdots,p_n都成立的下列不等式 $f(\frac{\sum\limits_{i=0}^n p_ix_i}{\sum\limits_{i=0}^n p_i}) \leq |\Omega|^-1 \int_\Omega f(x(t))dt \leq \frac{\sum\limits _{i=0}^n {p_if(x_i)}}{\sum\limits_{i=0}^n p_i}$ 式中\Omega是一个包含于n维单位立方体的n维长方体，其重心的第i个坐标为$\sum\limits _{j=i}^n p_j /\sum\limits_{j=i-1}^n p_i$,|\Omega|为\Omega的体积，对\Omega中的任意点$t=(t_1,\cdots,t_n)$, $w(t)=x_0(1-t_1)+\sum\limits _{i=1}^{n-1} x_i(1-t_{i+1})\prod\limits_{j = 1}^i {{t_j}} +x_n \prod\limits _{j=1}^n t_j$ 不等式中两个等号分别成立的情形亦已被分离出来。 此不等式是著名的Jensen 不等式的精密化。     

一个复数分解式的几何解释

设(k=0,1,2,…,n-1)是1的n次单位根。则 x~n-1=(x-1)(x-ε_1)…(x-ε_(n-1)) (1) 我们知道ε_n(k=0,1,2,…,n-1)表示单位圆的n个等分点,即单位圆的内接正n边形的n个顶点,其中ε_o=1。 (1)两边同除以x-1,得 x~n-1/x-1=(x-ε_1)(x-ε_2)…(x-ε_(n-1)) 根据G.P.求和式得x~2-1/x-1=1 x x~2 … x~n-1     

关于一类中心循环的有限P-群的自同构群

确定了一类中心循环的有限p-群G的自同构群.设G=X_3(p~m)~(*n)*Z_(p~(m+r)),其中m≥1,n≥1和r≥0,并且X_3(p~m)=x,y|x~(p~m)=y~(p~m)=1,[x,y]~(p~m)=1,[x,[x,y]]=[y,[x,y]]=1.Aut_nG表示Aut G中平凡地作用在N上的元素形成的正规子群,其中G'≤N≤ζG,|N|=p~(m+s),0≤s≤r,则(i)如果p是一个奇素数,那么AutG/Aut_nG≌Z_(p~((m+s-1)(p-1))),Aut_nG/InnG≌Sp(2n,Z_(p~m))×Z_(p~(r-s)).(ii)如果p=2,那么AutG/Aut_nG≌H,其中H=1(当m+s=1时)或者Z_(2~(m+s-2))×Z_2(当m+s≥2时).进一步地,Aut_nG/InnG≌K×L,其中K=Sp(2n,Z_(2~m))(当r0时)或者O(2n,Z_(2~m))(当r=0时),L=Z_(2~(r-1))×Z_2(当m=1,s=0,r≥1时)或者Z_(2~(r-s)).     

Interpolation with lagrange polynomials a simple proof of Markov inequality and some of its generalizations

The following theorem is provedTheorem 1.Let q be a polynomial of degree n(qP_n)with n distinct zeroes lying inthe interval[-1,1] and△'_q={-1}∪{τ_i:q'(τ_i)=0,i=1,n-1}∪{1}.If polynomial pP_n satisfies the inequalitythen for each k=1,n and any x[-1,1]its k-th derivative satisfies the inequality丨p~(k)(x)丨≤max{丨q~((k))(x)丨,丨1/k(x~2-1)q~(k+1)(x)+xq~((k))(x)丨}.This estimate leads to the Markov inequality for the higher order derivatives ofpolynomials if we set q=T_n,where Tn is Chebyshev polynomial least deviated from zero.Some other results are established which gives evidence to the conjecture that under theconditions of Theorem 1 the inequality ‖p~((k))‖≤‖q~(k)‖holds.     

一类无理函数的积分

<正> 众所周知,在不定积分的计算中,形如integral f(√a~2-x~2)dx,integral f(√x~2+a~2)dx 和integral f(x~2-a~2)~(1/2)dx 的积分可作代换x=asint(或x=acost),x=atgt(或x=asht)和x=asect(或x=acht)将其积出。而形如integral f((a-x)~(1/2),(b-x)~(1/2))dx,integral f((x-a)~(1/2),