在快速富氏变换(FFT)领域中,以任意数M(M=2~n,n为任意正整数)为基的一般计算式问题,至今尚未解决.本文提出此问题并推证了它.文中的计算量比赵访熊李庆扬改进的FFT计算公式大约减少40％.文中得到了与“Bergland-Brigham结论”迥然不同的结果,并解决了实用中最优基的选取问题. 考虑离散富氏变换(DFT):

高阶折叠点的计算产生于两参数非线性问题这里λ,μ∈R,u∈Banach空间X,f是R×R×X→X的C~3非线性映照. 近年来,两参数非线性问题有越来越多的实际应用,例如化学放热反应中的引燃问题:

首先考虑二阶常微分方程第一边值问题:假设(1)有唯一解,且解u(x)∈C~6[a,b];f(x,y,z)作为x,y,z的函数属于 C~2,

已知杂交元可以看作是非协调元,是否每一个非协调元均可作为杂交元的特例?容易明白,许多协调元和非协调元不能作为[1]中提出的杂交元的特例.例如,[1]中例 6的Wilson矩形非协调元就是如此.本文要拓广[1]中提出的杂交元的抽象框架.使许多协调元和非协调元都能作为杂交元来处理.从而不但能得到未知量的近似值,而且能同时

本文利用连分式插值,得到了一个新的一维搜索方法——连分式算法.用此算法,每迭代一次,只需计算三个点的函数值;在计算连分式插值式的每个系数时,只需一次除法.因此,数值稳定性较好.本文还证明了此算法的收敛性,收敛速度较快,收敛阶近似1.8393.按效能指标E=P~(1/μ)评价,此算法是一个较好的局部一维搜索方法.如果用此法于不精确的一维搜索,因只需计算三个点的函数值,故它是一个较好的、不精确的一维搜索方法,同时也是解超越方程的一个新算法.数值例子表明,它确实有效.

共轭梯度法在解高阶稀疏线性方程组方面有许多其它经典的迭代法所没有的优点,但当线性方程组相当病态、系数矩阵条件数很坏时,共轭梯度法的收敛速度很慢.因此,又产生了预条件处理共轭梯度法. 我们用预条件处理共轭梯度法求解线性方程组Ax=b(这里A是对称正定稀疏阵且条件数很大).预条件处理共轭梯度法旨在寻找一适当的正定矩阵C,C通常写成

本文使用下列符号:R~(n×m):所有n×m实矩阵的全体;R_r~(n×m):所有秩为r的n×m实矩阵的全体;||·||_2:向量的欧氏范数和矩阵的谱范数;||·||_F:矩阵的Frobenius范数;

设f是复平面C的某一区域D上的半纯函数;m,n为非负整数;为插值点,在其上f解析;r_(mn)~*=P_(mn)~*/Q_(mn)~*为所有满足下述方程的r_(mn)=P_(mn)/Q_(mn)(P_(mn),Q_(mn)≠0分别为次数≤m和≤n的多项式)中分子及分母次数极小者:

对于求解非线性方程组的区间迭代法,若利用Moore检验,可判断解的存在唯一性.[2]中在偏序下给出的区间Newton型方法,也有同样特性,本文利用f:D?R~n→R~n的斜度构造的区间割线算子,也可用于检验方程组解的存在唯一性,但它不用计算f的导数,针对f的不同分裂,还可以构造不同的两侧逼近割线法.分裂得当,便于求逆,使计算

Passow和Taylor在[1]中对没有变号区间的连续函数建立了共正逼近交错定理;同时指出,对一般情形(包括有变号区间的连续函数)建立交错理论,相当困难.随后史应光在[2]中对一般情形描述了在被逼近连续函数的变号点处导数非零的最佳共正逼近的特

对于实际问题中遇到的不少椭圆型方程,平均值定理成立,即在未知函数定义域内,任何球心的值等于球面上值的平均或某种形式的加权平均.我们知道,上述平均值的条件,对于未知函数所满足的方程是充分必要的.就是说,平均值定理对于未知函数具有本质的特点.因此,我们可以构造一种数值方法,不从未知函数满足的方程出发,而直接从它满足的平均值定理的表达式出发,来做近似计算.

1980年以来,唐立民等提出一种拟协调元法,用来构造椭圆型方程的离散格式.粗略地讲,该法将每个单元上的能量表达式所含导数项的面积分(假设问题二维的),用格林公式转化为单元边界上的线积分,然后采用某种数值积分,将线积分进行离散.对只含函数项的面积分,也用相应的数值积分进行离散.用此法计算单元刚度阵,比较简单、灵活.

用R~(n×m)表示所有n×m实矩阵的集合;OR~(n×n)表示所有n×n正交矩阵的集合;S_(n,r)表示所有带宽为2r+1的n阶实对称矩阵的集合;||·||_F表示矩阵的Frobenius范数,||·||表示向量的Euclid范数.任取A∈R~(n×m),满足AA~-A=A 的A~-∈R~(m×n)叫做A的内逆,满足AA_l~-A=A和(AA_l~-)~T=AA_l~-的A_l~-∈R~(m×n)叫做A的最小二乘广义逆,

1.问题的提出 [1]中提出,选择不同的基函数,即可构造出数值求解微分方程的不同公式.[1]中还讨论了一些新公式,其中有的优于一般的线性多步法,本文旨在给出其理论证明. 下面仍采用[1]中记号,即

在多孔介质中,考虑单相的、一种可压流体被另一种流体所溶混驱替的流动.设诸集层Ω是单位厚度且视其为R~2中的有界区域,忽略重力项,混合流体的Darcy速度可表为 u=-K(x)/μ(c)?p,其中p为压力,K(x)为介质的渗透率,μ为与浓度c有关的粘度.设dρ_i/ρ_i=z_idp,其

有关凹增算子及凸减算子的正固有元存在性及迭代序列的收敛性,对非单调算子不适用. 大部分非单调的积分算子可以表示成 T=T_1+T_2, (1)其中T_1增,T_2减.[1,2]讨论了迭代序列的收敛性,但[2]中关键部分的证明是不正确

在线性多变量控制理论中,存在一个代数特征值反问题——输出反馈极点配置问题。问题叙述如下: 问题PAO.给定A∈R~(n×n),B∈R_m~(n×m),C∈R_p~(p×n)和?={λ_1,λ_2,…,λ_n},?在复共轭下封闭.求K∈R~(m×p),使得A+BKC具有事先给定的特征值λ_1,λ_2,…,λ_n。 [2]和[5]等证明了

§1.引言 对具可微卷积核的第二类Volterra积分方程 y(x)=f(x)+λ integral from n=a to x(K(x-t)y(t)dt),(1)通常的解法有迭代法与Laplace变换法以及化为微分方程求解等.毫无疑义,这些方法对于方程(1)的求解是重要的.但这些方法也有其本质的缺点,即在求解过程中,往往涉

1.引言 记n维欧氏空间R~n的非空紧凸子集族为P(R~n).设F:R~n→P(R~n)是上半连续的集值映射.称x∈R~n为F的一个Kakutani不动点,如果x∈F(x). 考虑计算F:R~n→P(R~n)的Kakutani不动点的问题.熟知,Merrill重复开始

1.引言 多重网格法是求解椭圆型方程边值问题的一种有效的迭代解法,其特点是方法收敛速度与网格长度h无关,因此为达到具有相同精度的解只需O(N)次的运算量(N为离散后的线性方程组未知数个数).从而比一般的迭代法有效得多.现在这个方法已被广泛

唐立民等于1979年提出拟协调单元的概念.拟协调单元是非协调单元的一种,虽然在单元边界上保持了积分意义下的连续性,但是仍然产生变分“犯规”.本文第一节将对拟协调单元作某些进一步的讨论,并给出一般的收敛定理.第二节将研究准协调单元.

1.引言 常微分方程初值问题并行数值方法的研究,一直是并行算法研究中值得注意的问题.其原因不仅在于常微分方程初值问题是一个典型的非线性连续递推问题,也在于它在应用中的重要性,特别如实时计算的需要. [1]与[2]对两类典型的线性多步公式,Adams-Molton隐式公式和 Gear公式(即向后微分公式)进行处理,得到了一类并行算法.其基本思想是将这两类线性多步公式在一个区间上作为非线性方程进行整体迭代求解,该方法的最大特点是方程右端函数在各节点上可以并行计算,适用于多处理机系统和流水线向量机.[2]在一定的迭代初值条

1.引言 最优化计算中的直接法,是指这样的一类方法,它们向极小点走近的手段只依赖问题

1.引言 迭代求解线性方程组Ax=b的AOR方法已是众所周知.由AOR迭代很自然联想到构造对称AOR(SAOR)迭代,但目前讨论SAOR迭代的文章还不多见.中对系数矩阵为H阵的SAOR迭代,[6]中对系数矩阵为对称正定阵的SAOR迭代,均给出了收敛性定理.本文讨论系数矩阵为对角元素非零的相容次序阵时SAOR迭代的收敛性,得到了相应的收敛性定理,并给出了SAOR迭代矩阵谱半径表达式以及谱半径的一个上下界.

前言 七十年代中期,人们多采用显式方法数值求解可压缩的Navier-Stokes方程.这种方法简单易解,但由于稳定性对时间步长的限制,使得求解所需机时颇多.在求解定常问题时,数值求解过程可以与真实的物理发展过程不对应,人们可以根据需要而改变求解过程,以达到加速收敛的目的.Allen和Cheng就是根据这种思想计算了近底部分离流动.为了达到加速得到定常解的目的,很多人采用在不同空间点上取变时间步长的方法.在[5]中,当调节因子取标量形式时,相当于取变时间步长的方法.如果调节因子或算子放大修正系数取矩阵形式,则可得到更快的收敛速度.Beam和Warming在[7]中提出了一个非迭代的隐式方法,并在空间坐标方向上利用近似因式分解,大大提高了隐式格式的使用效率.Steger和Warming在[8]中详细介绍了流通量分裂法.1985年,MacCormack在[9]中改进了自己在[10]中提出的二步隐式方法.作者在[11]中也用流

本文给出两种在向量计算机上计算n阶矩阵乘积的并行算法: 1)最优内积算法. 处理机台数 ρ_1=n~3/log_2n,

设?~k是[0,1]上的CooeB空间,Q:?~k→R是至少具有k-1次代数精度的求积泛函.设J:f|→integral from n=0 to 1 (f(t)dt),h=1/n。通过由等式 M_hf(t)=h sum from i=0 to (n-1)(f(ih+th)),?f∈C[0,1],?t∈[0,1]确定的线性算子M_h:C[0,1]→C[0,1],定义Q的复化求积泛函QM_h。在?~k中的

关于迭代法的点估计,在解方程算法的效率研究中起着关键的作用.这一研究是Smale的连续复杂性理论最成功的范例. 设f:E→F是Banach空间之间的解析映照.对z∈E,令α=α(z,f)=β·γ,其中

1. 引言 谱方法为非线性偏微分方程的求解提供了新的技巧.由于拟谱方法比谱方法便于实施,计算量小,所以应用更为广泛.但它有时会产生非线性不稳定性.为此,一些滤波和抑制方法接连出现.本文对Burgers方程的周期边界问题,建立了一个带抑制算子的三层拟谱格式,同时证明了格式的广义稳定性.在一定的条件下,由此稳定性可得到收敛性。

本文运用权范数方法证明了多角形区域上Green函数有限元逼近的逐点估计 |G_z(x)-G_z~h(x)|≤C(h~|lnh|~3/|x-z|~α),?x,z∈Ω,其中C为与x,z,h无关的常数θ<α<β_M,β_M=π/α_M,α_M为Ω的最大内角.由此可导出凹角域上有限元逼近的渐近展开. 考虑模型问题

1.引言 给定R~n的一个非空集K及一个从R~n到其自身的映射f,变分不等式问题,记为VI(K,f),就是求一个向量x~*∈K,使得 (y-x~*)~Tf(x~*)≥0,?y∈K.(1.1) 本文假定映射f是非对称的(即不可积的),因而上面的变分不等式问题一般不能对应于一个最优化问题. 集K假定为一些低维集的Cartesion积

1.引言在[1]-[6]中讨论了色散方程u_t=au_(xxx)(a为常数,可正可负)的差分解法,但是, 显式格式的稳定性条件较苛刻,其中以[5]中提出的 H_3类显式格式最好,稳定条件为|R|=|a|τ/h~3≤1.1851;而隐式格式虽然绝对稳定且具有高精度,但每前进一步需要解一个具有五对角线的线性方程组,计算量较大. 本文针对显式格式与隐式格式存在的问题,提出一类三层绝对稳定半显式格式,其截

1、引言 用D表示平面上矩形域[0,l_1]?[0,_2 ],记l=max(l_1,l_2),m,n∈N(自然数集),h_1=l_1/m,h_2=l_2/n.用直线族x=ih_1,y=jh_2(t=1,…,m-1;j=1,…,

1.提出问题 本文讨论非线性抛物型问题的变网格混合有限元法,即,在用混合有限元法求解的同时,于不同时刻采用不同有限元网格和插值函数;同时提出了四种全离散格式,并给出了理论分析和误差估计,且证明了这种估计在某种意义下是最佳的。 设Ωo为二维多角形区域,?Ω为其边界、考虑下列初、边值问题:

1.构造格式 考虑如下波动方程 u_(tt)=u_(xx) (1.1)的初边值问题,设其边界条件为周期的,即在此条件下,解具有周期性.(1.1)有二种namilton形式.一种是经典形式:

§gi.两类逆特征值问题先说明一些记号.R~(m×n)是所有m×n实矩阵的全体,R~n=R~(n×1),R=R~1;SR~(n×n)是 所有n×n实对称矩阵的全体;OR~(n×n)是所有n×n实正交矩阵的全体;I~((n))是n阶单位矩阵;A~T是矩阵A的转置;A>0表示A是正定的实对称矩阵.?(A)是矩阵A的列空间;A~+是矩阵A的Moore-Penrose广义逆;P_A=AA~+表示到?(A)的正交投影.λ(A)是A的特征值的全体;λ(K,M)是广义特征值问题K_x=λM_x的特征值的

1.引言 假定F是R~m→R~m的可微映射,x~*∈R~m是 F(x)=0 (1.1)的解. 如果在解点处Frechet导数是可逆的,只要F′(x)具有一定的性质,就可保证Newton迭代 x_(i+1)~N=x_i~N-F′(x_i~N)~(-1)F(x_i~N) i=0,1,… (1.2)产生的点列在||x_0-x~*||适当小时二阶收敛于x~*:

1.引言 本文所讨论的问题如下: Min f(x) x∈R~n, s.t. c_i(x)=0,i=1,…,q,(1.1) c_i(x)≤0,i=q+1,…,p.解此问题的递归等式约束二次逼近算法,是由Murry(1969)提出,而后由Biggs(1972)发展的.此项研究是从罚函数的轨迹出发,建立一个只包含等式约束的二次规划子问题,从而可用代数的方法求得搜索方向.并沿该方向作线性搜索而完成一次迭代过程.Biggs将二次罚函数作为效应函数用于线性搜索,并证明了该算法具有全局收敛性和局部超线

1.二重Toeplitz矩阵相乘的快速算法nm阶方阵 称为nm型2重Toeplitz矩阵,其中A_i(i=-n+1,…,n-1)为m阶Toeplitz矩阵. 定义.设p_1×p_2矩阵A=(a_(ij))_(p_1×p_2),B为q_1×q_2矩阵.称p_1q_1×p_2q_2矩阵

本文根据[1]和[2]中提出的有限解析法的基本思想,推导出求解非线性常微分方程两点边值问题的简便公式,并且从理论和实际计算证明了这个方法精度高、收敛快、稳定性好.对于用通常方法求不出合理结果的问题,用此方法可以求出很好的解.