首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 15 毫秒
1.
曾文平 《应用数学》2002,15(4):52-56
对一类演化方程δu/δt=aδ^2m 1u/δx^2m 1(a为常数,m为正整数)构造一族含双参数的三层高精度隐式差分格式。当参数α=1/2,β=0时得到一个双层格式。并证明了:该格式对任意非负参数α≥0,β≥0都是绝对稳定的,并且其截断误差阶为0((Δt)^2 (Δx)^4).数值例子表明:本文所建立的差分格式是有效的,理论分析与实际计算相吻合。  相似文献   

2.
曾文平 《应用数学》2004,17(2):250-256
对高阶Schrodinger方程эu/эt=i(-1)^mэ^2mu/эx^2m构造一族含双参数的三层高精度隐式差分格式.当参数α=1/2,β=0时得到一个两层格式.并证明了:对任意非负参数α≥0,β≥0该格式都是绝对稳定的,并且其截断误差阶达到O((△t)^2 (△x)^6).数值例子表明:本文所建立的差分格式是有效的,理论分析与实际计算相吻合.  相似文献   

3.
解抛物型方程的一族高精度差分格式   总被引:8,自引:0,他引:8  
1 引言 求解抛物型方程 u/t=u/x~2, 00, (1) 初边值问题的差分格式,精度高者当属[1]、[2]中的格式.本文对上述问题构造了一族三层(特殊情况下是两层)双参数、绝对稳定、高精度三对角线型的隐式格式,它不仅包含了[1]、[2]中所有的格式,而且还可以得到一个截断误差为O(Δt~3+Δx~4)的绝对稳定的差分格式,精度比[1]、[2]中的格式都高. 2 差分格式 设Δt为时间步长,Δx=L/M(M为正整数)为空间步长,网函数u(jΔx,nΔt )记为u_j~n,对  相似文献   

4.
1 引言 1960年,Saul’ev在文中讨论了如下的高阶(2m阶)抛物型方程 μ/t=(-1)~(m 1)~2mμ/x~(2m) (1)(其中m为正整数),提出了一类含极因子α的两层差分格式。当α=0时为显式格式,其稳定性条件为,r=△t/(△x)~(2m)<1/2~(2m-1),△t,△x分别为时间及空间步长。随后,文[2],[3]利用  相似文献   

5.
解三维抛物型方程的一个新的高精度显格式   总被引:6,自引:1,他引:5  
本文构造了一个解三维抛物型方程的高精度三层显式差分格式,其稳定性条件为r=Δt/Δx2=Δt/Δy2=Δt/Δz2≤1/4,截断误差为O(Δt2+Δx4).  相似文献   

6.
针对四阶抛物型方程周期初值问题,提出了一个两层隐式差分格式和一个三层隐式差分格式.它们的局部截断误差分别为O((Δt)2+(Δx)4)和O((Δt)2+(Δt)(Δx)2+(Δx)4),其中Δt,Δx分别为时间步长和空间步长.误差分析和数值实验均表明,本文构造的差分格式比经典的Crank-Nicolson格式和Saul’ev构造的差分格式精度更高.从精度及稳定性方面考虑,本文构造的格式也比文[5]的显式格式要好.  相似文献   

7.
多维抛物型方程的分支绝对稳定的显式格式   总被引:24,自引:0,他引:24  
其中及R={0≤x_i≤1,j=1,2,…,p),(?)R只为区域只的边界。 对多维抛物型方程(1)的差分解法,古典显式格式的稳定性条件为r=Δt/(Δx)~2≤1/2p,十分苛刻;古典隐式格式虽是无条件稳定,却需解线性方程组。因此两者的计算量都很大,且它们的精度较低,其局部截断误差仅为O(Δt+(Δx)~2)。因此,对多维抛物型方程而言,构造显式计算、稳定性能良好且精度较高的差分格式便具有十分明显的理论意义和实用价值。本文针对上述古典显式与隐式格式所存在的问题,构造一类对任何p维空间变量的抛物型方程(1)都适用的。分支绝对稳定的显式差分格式,其局部截断误差阶为O((Δt)~2+(Δx)~2),从而避免了解线性代数方程组,大大地减少了计算工作量,且精度较高。 令Δx_k=h_k=Δx=h=1/M(k=1,2,…p)表示空间方向步长,Δt=τ=[T/N]表示时间方向步长,M、N均为正整数。 为简便计,引入下列记号  相似文献   

8.
本文以下列简单问题为例,讨论Petrov-Galerkin方法的稳定性:?U/?t q?U/x=ε~2U/x~2,-∞0,U(x, t)=U(x 1.t), -∞相似文献   

9.
For α≥β≥ -1/2 let Δ(x) = (2shx)2α+1(2chx)2β+1 denote the weight function on R+ and L1(Δ) the space of integrable functions on R+ with respect to Δ(x)dx, equipped with a convolution structure. For a suitable φ∈ L1(Δ), we put φt(x) = t-1Δ(x)-1Δ(x/t)φ(x/t) for t > 0 and define the radial maximal operator Mφ as usual manner. We introduce a real Hardy space H1(Δ) as the set of all locally integrable functions f on R+ whose radial maximal function Mφ(f) belongs to L1(Δ). In this paper we obtain a relation between...  相似文献   

10.
本文在非一致时间网格上,使用有限差分方法求解变时间分数阶扩散方程?α(x,t)u(x,t)/tα(x,t)-2u(x,t)/x2=f(x,t),0α(x,t)q≤1,证明了该方法在最大范数下的稳定性与收敛性,收敛阶为C(Δt2-q+h2).数值实例验证了理论分析的结果.  相似文献   

11.
研究了一类具有连续变量的二阶中立型时滞差分方程△(2Υ)(x(t) -px(t- γ) =mΣi=1 qi(t)x(t-σi),t ≥ t0 > 0的振动性,给出了其有界解振动的几个充分条件.  相似文献   

12.
该文主要讨论带临界指数的椭圆型方程组{-Δu + a(x)u =2α/α+βuα-1vβ + f(x),x ∈Ω,-Δv+b(x)v=2β/α+βuαvβ-1+ g(x),x ∈ Ω,(*)u > 0,v > 0,x ∈Ω,u=v=0,x ∈(a)Ω解的存在性,其中Ω是RN中一个光滑有界区域,N=3,4,a≥2,β≥2...  相似文献   

13.
对四维抛物型方程构造了一个高精度显格式,格式的稳定性条件为r=Δt/Δx2=△t/Δy2=△t/△z2=Δt/Δw2<1/2,截断误差阶达到O(Δt2 Δx4),通过数值实验,表明本文理论分析的正确性和文中格式较同类格式的优越性.  相似文献   

14.
该文研究了加权的退化椭圆系统■其中Δ_Gu=Δ_xu+(a+1)~2|x|~(2a)Δ_yu是Grushin算子,α,β≥0,q1,ω(x)=(1+‖x‖~(2(α+1)))~(β/2(α+1)).超临界指数正稳定解的Liouville定理被建立.  相似文献   

15.
1 引  言在地下水含水层中 ,污染物随地下水运移并常常发生各种化学反应 [1 ] .描述地下水含水层中一类阳离子交换反应 m M1 +r M2 k2k1 r M2 +m M1 的数学模型[2 ] 为 : s1 t- dΔs1 =f1 ,   x∈Ω ,t∈ J (1.1a) s2 t- dΔs2 =f2 ,   x∈Ω ,t∈ J (1.1b) c1 t+ρ s1 t- DΔc1 =0 ,   x∈Ω ,t∈ J (1.2 a) c2 t+ρ s2 t- DΔc2 =0 ,   x∈Ω ,t∈ J (1.2 b)其中 Ω R2为具有光滑边界的有界区域 ,J=(0 ,T].这里 D>d>0为扩散系数 ,ρ>0为固体颗粒密度 ,均为常数 .根据 [1,2 ]应有 :f1 =m[k1 Rm1 Rr2 cm1 sr2 -…  相似文献   

16.
20 0 0年 2月号问题解答(解答由问题提供人给出 )12 3 6.设 f( x) =( x2 2 x 3 ) x 3 ( x2 2 x 3 ) x2 ,当 x∈ R.证明 :f ( x)≥ 6.证明 ∵ x2 2 x 3 =( x 1) 2 2≥ 2∴ 对 x∈ R有  ( x2 2 x 3 ) x 3 >0成立 .因此1° 当 x<-3时 ,f ( x) >0 2 x2 >2 9>6,这时命题成立 .2° 当 x≥ -3时 ,f ( x) =〔( x 1) 2 2〕x 3 〔( x 1) 2 2〕x2令  x 1=t,由 x 3≥ 0 ,则 t 2≥ 0那么  f ( x) =g ( t) =( t2 2 ) t 2 ( t2 2 ) (t-1 ) 2≥ 2 t 2 2 t2 -2 t 1 =4· 2 t 2· 2 t2 -2 t= ( 2 t 2 t 2 t 2 t) ( 2 t2 -2 t 2 …  相似文献   

17.
该文研究了以下高阶Yamabe型方程Lm,pu - g|u|p-2 u = λ f|u|α-2 u在有限图上的非平凡正解的存在性,其中Lm,p是一个2m阶差分算子,它是一种p次(-Δ)m算子更一般化,α≥p≥2,g>0和f >0是定义在G的所有顶点上的实函数,m≥1是一个整数.  相似文献   

18.
奇异半线性反应扩散方程组Cauchy问题   总被引:1,自引:0,他引:1  
本文讨论如下问题其中{(б)u/(б)t-(1/tσ)△u=αvp1+β1vp1+f1(x),t>0,x∈RN,(б)u/(б)t-(1/tσ)△v=α2uq2+β2vp2+f2(x),t>0,∈RN,limt→0+u(t,x)=limt→0+v(t,x)=0,x∈Rn,其中σ>0,pi>1,qi>1(i=1,2),α1≥0,α2>0,β1>0,β2≥0,fi(x)(i=1,2)连续有界非负,(f1(x),f2(x))(≡/)(0,0).给出了非负局部解存在的几个充分条件和解的爆破结果.  相似文献   

19.
W.Janous猜测的再推广   总被引:2,自引:0,他引:2  
W .Janous猜测 :设x ,y,z >0 ,则x2 -z2y+z + y2 -x2z+x + z2 -y2x+y ≥ 0 ( 1 )贵刊文 [1 ]将 ( 1 )式推广为 :设xi>0 (i =1 ,2 ,… ,n) ,n≥ 3,记S=x1 +x2 +… +xn,t =x21+x22 +… +x2 n,则nx21 -ts-x1 + nx22 -ts-x2 +… + nx2 n-ts-xn ≥ 0 ( 2 )当n =3时 ,由 ( 2 )式可得 ( 1 )式 .本文将把 ( 2 )式作进一步推广 ,并回答文 [1 ]中提出的一个问题 .我们得到如下结果 :设xi>0 (i=1 ,2 ,… ,n) ,n≥ 3,α ,β,γ∈R ,T =xα1 +xα2 +… +xαn,A>max{xγ1 ,xγ2 ,… ,xγn},则当αβγ >0时 ,有nxα1 -T(A-xγ1 )β+ nxα2 -T(A-xγ2 )…  相似文献   

20.
一类二次方程组的一个定理及其运用   总被引:1,自引:0,他引:1  
定理 在方程组∑ni=1xi=A∑ni=1x2i=B中 ,A、B是实数 ,记Δ=n B-A2 .若 xi∈ R( i=1,2 ,… ,n) ,则Δ≥ 0 ,当且仅当x1 =x2 =… =xn=An时 Δ=0 .证明  ∑1≤ i相似文献   

设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号