一、曲线的方程 给定平面点组p_i(x_i,y_i)(i=0,1,2,…,n),它们构成拆线p_0p_1p_2…p_n。考虑x_0

1.前言 关于系数用偶阶导数表示的(2n+1)次样条函数,使用性能较好,但其存在性与唯一性迄今尚未给出证明.这种样条函数与一般的奇次多项式样条函数一样,当n≥2时,方程组的系数矩阵已不具有明显的主对角元素占优,致使J.H.Ahlberg等人说,直接依赖系数矩阵的性质来证明多项式样条函数的存在性是十分困难的,即便是对于

在[1]中考虑了亏度为3的五次插值样条.郭竹瑞考虑了两类亏度为2的三次插值样条.本文的目的是考虑三类(分别记为Ⅰs,Ⅱs,Ⅲs)亏度为2的四次插值样条.本文利用多项式HB插值问题与Spline HB插值问题的密切联系,并以前者为基础,证明了Ⅰs—Ⅲs插值样条的存在唯一.结果表明,Ⅰs—Ⅲs插值法(整体属于c~2[a,b]的四次插值样

本文讨论非线性扩散和色散方程隐式单支差分格式的稳定性,将证明:当 α=1-1/2s时,扩散方程单支差分解满足L_(2s)范广义极值原理,从而L_(2s)范无条件压缩性稳定.证明的方法适用于多维方程和方程组,而且扩散系数可以是退化的.类似地给出色散方程绝对稳

它们分别作为有限元涡函数和流函数逼近的收敛性尺度,获得了最佳敛速估计和混合刚度模型的计算格式. 最近,Babuska和Osborn对于常微分两点边值问题,在[1]中构造的一般理论的基础上,用类似于范数(1.1),(1.2)的一维L_p范数,讨论了经典有限元逼近的敛速估计,得到了一些有意思的结果.各种型式的有限元分析可纳入统一的图式之中.

1.引言 设X为[a,b]中至少含有n+1个点的紧集,其中n为一个固定的自然数.在实连续函数空间C(X)中定义一致范数||f||=max[f(x)|。M是C[a,b]中的n维Haar子空间,h_1(x),…,h_n(x)是它的一个基底.

A.Hadjidimos在[1]中提出一个迭代求解线性方程组的 AOR方法(Accelerated Overre-laxation Method),并在方程组的系数矩阵为不可约弱对角优势、L-矩阵和相容有序矩阵的条件下,讨论了此方法的收敛性.在这篇文章里,我们将考虑系数矩阵是H-矩阵、正定矩阵以及L-矩阵的情况.所得结果表明,可以放宽在[1]的3,4两节中对参数所加的限制.

§1.引言 在[1]中尺度变换的基础上,本文讨论了最优尺度矩阵的情形,并利用[2]的结果,提出了证明矩阵分解收敛性的一种新方法.这个方法简单易行,用它可以证明和推广许多矩阵分解的收敛性.例如,[3]-[9]中的矩阵分解,大大简化了那里的证明步骤.

前言 计算技术的发展,为数值求解复杂的可压缩Navier-Stokes方程提供了良好的条件.在解决实践中提出的空气动力问题时,我们仍面临很多困难.其一是复杂的几何图形为边界条件的处理带来了麻烦;二是为了得到一组比较满意的定常状态的解,往往需要花费大量的计算机时,即使是解二元问题,机时消耗也非常突出.最后,利用现在广为采用的

初值问题的差分解法,参数ε≥0,μ≥0. 这一方程当ε=μ=0时为KdV方程,δ=ε=0时为Burgers方程,而当δ=μ=0时为RLW方程.对于方程(1),已设计了许多计算格式.对于KdV方程,最早的格式当推Zabusky-Kruskal,后来有[2—6].对于RLW方程,也有许多工作.对于Burgers方程,格式就更多了.非线性波动方

本文把普通集合中的离散Fourier变换推广到模糊集合。借助于区间数、模糊数的运算规则及有关性质,给出了模糊离散Fourier变换(FDFT)的定义及算法,而且也讨论了模糊离散Fourier变换中的对应关系以及变换性质的几个定理。

一、前言 关于线性规划解唯一的问题,[1]中作了充分的讨论,给出了解唯一的充要条件,并且相应给出了十种等价形式.本文的主要目的是: 1)给出解唯一的充要条件的另一证明,从而得到充要条件的几何意义; 2)从几何意义出发,同样可得到十种等价形式.

§1.引言 1975年,Enright等人在[1]中通过各类典型的试验问题对常微分方程Stiff方程组的五个数值方法进行了数值比较,得到的结论之一是Enright方法很好.但也有弱点,即在使用它处理高度非线性Stiff系统时,与Gear方法和基于梯形公式的整体外插方法相比较,效果最差.其原因之一是对于任意一个Stiff系统来说,使用该方法在采用Newton-

一、引言 已有许多方法可以求多项式的根,其中多数属于迭代法.这类方法的缺点是收敛性依赖于初始近似的选择(求复根还没有大范围收敛的迭代法).另一类方法是先求根的模或实部,然后再设法求同模或等实部的根,如根平方——结式法,按分布理论求根等就属于这一类.这类方法不存在收敛性问题,但也有其不足之处.如根平方——结点法的缺点除了有可能将良态多项式变为病态多项式之外,求结式也增加了计算的复杂性.采用单精度运算的数值试验结果表明,这个方法求解的精度比某些迭代法低.

在设计凸轮曲线的过程中,[1]提出了对一次导数插值的三次样条,又[2]中指出作为计算机辅助设计,需研究对二次导数的插值样条,我们把这两种插值法分别称为(P)型及(Q)型.

本文是同名文章[1]的继续.所用记号均沿用[1]中规定.我们的主要目的是讨论Ⅰs,Ⅱs,Ⅲs三种类型插值法的误差阶,并给出与eeB模、连续模有关的误差估计.本文证明了:对于Ⅰs—Ⅲs型插值样条,下列不等式:

一f(x)是区间[0,1]上定义的函数,0=x_0

本文研究一类非线性双曲型方程混合问题的有限元方法的稳定性和收敛性理论.关于线性双曲型方程有限元方法的收敛性研究,已有T.Dupont,J.T.Oden等人的工作.

§1.引言 非线性方程组 F(x)=0,F:D?R~n→R~n (1.1)嵌入参数t,构成同伦H:[0,T]×D?R~(n+1)→R~n,使得 H(0,x~0)=0,H(T,x)=F(x),(1.2)这里T可以是有限的或+∞,当T为+∞时以极限过程代替求值.若 H(t.x)=0(1.3)存在连续解x(t):[0,T]→D,则非线性方程组(1.1)的解x~*=x(T).若(1.3)的解

矩阵特征值最小距离的估计是研究谱理论和数值计算特征值的重要课题.本文对一类三对角线对称矩阵的特征值最小距离给出较精确的上下界估计,并根据数值计算提出一个关于特征值最小距离的位置的猜想.

1.引言 我们考虑以下形式的等式和不等式线性方程组: sum from j=1 to (a_(ij)x_j)=b_i,i=1,2,…,l, (1.1) sum from j=1 to (a_(ij)x_j)≤b_i,i=l+1,…,m.(1.2)对于求解这类问题,较早的算法有消去法和松弛法(即投影法).消去法在[1]中有详细的叙述.由于它每消去一个变量,不等式的个数就急剧地增加,因而不易在计算机上实现.松弛法虽然计算公式比较简单,但由于它的收敛速度较慢,在应用上有一定的局限性,

近年来,在计算数学刊物上相继发表了许多篇关于三次插值样条存在唯一性的文章,例如[1-3].这些文章讨论的是三次样条插值矩阵为非奇异的条件.[1]中用的是凑方法,讨论了与插值矩阵相关的另一个对称阵为正定的条件,经过复杂凑方,得到了某些充分条件,[2]是用大块凑方,所得结果形式上异于[1],但实质上是完全相同的.[3]则是对插值矩阵进行一种特殊分解,得出非异的四个充分条件.它不限于[1-2]所讨论的正定情形,因而适用范围更广些.

本文讨论了一类非线性切比雪夫逼近问题的两个具体算法,证明了它们的收敛性.实际计算表明,效果很好. 设f(x)与F(A,x)为给定的两个实连续函数,其中x∈M,A=(a_1,a_2,…,a_n)∈?,M表示实轴上的某个有界闭集,?是n维欧氏空间中的子集.切比雪夫逼近问题就是在?中确定A~*∈?,使对一切A∈?,有

关于正规矩阵的任意扰动,有下述定理成立. 定理1.设A为n阶正规矩阵,C为n阶任一矩阵.A的特征值为λ_1,…,λ_n,C的特征值为μ_1…,μ_n.C~H表示C的转置共轭,||·||_2与||·||_F分别表示矩阵的谱范数与Frobenius范数.记

考察典型的抛物型问题:其中Ω为平面有界区域.设S_h?H_1~0(Ω)是在正规剖分上由分片m-1次多项式构成的有限元空间,其半离散Galerkin逼近可由下式确定:

0.李岳生在“微分算子样条函数磨光法”的一般框架下,推导了“多项式样条磨光公式”:以及处理f由离散数据给出的相应公式?_hf(x),这里m是任意正整数,q是任意非负整数,h为任意正数,

§1.引言 在二维区域Ω(边界为?Ω)上,根据热量守恒原理,描述变温过程的积分关系式为初条件为

Let {△_n} be a Sequence of partitions of the interval [0, 1] and ?_(△_n).(f; x) denotesthe cubic periodic interpolating splines of f for the partitions △_n. Let q be a positive number,p be the largest integer smaller than q, and set α = q -- p, thus q = p + α 0 < α≤1. The following theorem has been proved. Theorem. If f ∈ C_[0,1]~* then for every sequence of partitions of the interval [0, 1] withlim△_n = 0 n→∞, the necessary and sufficient conditions for ?_△_n(f; x) satisfyingare 1) f satisfies Lipschitz condition, if q ≤ 1. 2) f ∈ C_[0,1]~* and f~((p)) satisfies Lipschitz condition of order α, if 1 < q < 2 or 2 4,where C_[0,1]~(*p) = {f(x)|f ∈ C_[0,1]~p, f(0) = f(1), f'(0) = f'(1)… f~((p))(0) = f~((p))(1)},

其中_q表示整数k模q的最小非负剩余.这种卷积的计算,在通用计算机及专门结构的数字处理装置中有广泛的应用,特别可用于图象处理方面.直接计算(1),需要p~2q~2次

本文给出NDE方程组的几种解法,即欧拉法、改进的欧拉法以及一般方法.我们考虑NDE方程组:

平面Bzier曲线的凸性定理是计算几何中一个重要的定理.最近,苏步青、刘鼎元在[1,2]中给出了凸性定理的证明. 本文的目的是从Bezier曲线的一阶导矢与二阶导矢的几何作图出发,给出凸性定理的另一个证明.

关于Hermite-Fejer插值多项式,已有许多人相继作了不少的工作.其中,Moldovan,Bojanic,Saxena等人各自建立了一些以第一、二类Chebyshev多项式T_n(x)和U_n(x)的零点:

从复变函数论知,复多项式f(u)=sum from l=0 to n (a_1u~(n-1))在u=u_0点的Taylor展开式是u的恒等式:

我们知道,用代数多项式逼近常与点的位置有关,因此考虑点态的情形是有意义的.最近,谢庭藩指出

众所周知,在使用线性方法(如线性多步法,Runge-Kutta方法,合成多步法等)对Stiff常微分方程组初值问题进行数值积分时,为了保证该初值问题数值解是稳定的,则要求数值方法在某种意义下是无条件稳定的.为此,所使用的线性方法首先必须是隐式的.在使用隐式线性方法对Stiff系统初值问题进行数值解时,每向前积分一步,往往

由于正定几何规划的对偶规划只含线性等式约束和非负约束,处理起来似乎要方便得多.然而,实际上许多对偶算法实施起来却往往失败(见[2,8,9]),这是由于对偶规划所特有的“块性质”以及目标函数在某些点的不可微性质引起的.因此,近年来主要的努力集中在克服这二个困难上。主要的工作有:1975年Beck和Ecker的修正凹单纯形

用插值分式表或切触插值分式表来讨论有理插值或切触有理插值问题的一些算法的条件是比较方便的(参看[3],[6]).但关于这两个表的结构,至今未见充要的结果.为解决此问题,先引入有关术语及记号,并首先考虑有理插值的情况.

[1]—[4]均研究本问题.[1]采用体积坐标给出了基函数的显式表达,较之[2]—[4]省去了求大型逆矩阵一步手续.我们知道,对有限元法来说,一个插值函数的光滑性和逼近度具有很重要的意义.但对16个自由度的四面体元,他们给出的插值只能保证相

是A稳定的,其外插可分为整体和局部两种.整体外插是在整个区间上用梯形法完成积分后再进行外插,仍保持A稳定性;局部外插是在每一节点上进行外插,再将此外插值作为以后积分的初值.Lapidus和Seinfeld指出,梯形法局部外插不保持A稳定性.因此

对于三维连续介质的有限元分析,一个通用的二次有限单元体是所谓20节点等参元.尽管这个元素已富有成效地被普遍应用,但是,它需要的节点自由度太多与其达到的二次多项式的逼近精度却十分不相称,显得计算效率很低.基于混合刚度有限元法的一