- 赵金熙
<正> 1 引言 在实际应用中常会提出解等式约束加权线性最小二乘问题 min(b_2-A_2x)~TW(b_2-A_2x) x∈R~n (1) s.t.A_1x=b_1,其中A_1∈R~(p×n),A~2∈R(q×n),b_1∈R~p,b_2∈R~q,W∈R(q×q)为对称正定矩阵. 对于问题(1),目前已有多种数值求解方法,如Paige利用(1)的对偶公式给出了一个向后稳定的数值方法.Gulliksson和Wedin利用加权QR分解技巧给出了解(1)的一个直接解法.作者利用广义Cholesky分解构造了解(1)的矩阵分解方法.
1996年02期 97-103页 [查看摘要][在线阅读][下载 142k] - 方逵
<正> 1 引言 计算机图形学的一个基本问题就是寻找一条光滑曲线过一组型值点{x_i,y_i}(i=0,1,…n+1),解决这一问题最简单的办法是用分段三次Hermite插值,这种插值构造容易,绘图简单. 分段三次Hermite插值的关键是估计型值点处的导数,只要估计出一组导数值,就对应一个分段三次Hermite插值.但在实际应用中,必须考虑插值曲线对型值点组某些特征的继承性,如曲线的保凸性,保形性等. [1—2]研究了分段三次Hermite插值的保单调性.[3]导出了分段三次Hermite插值保形的一个充要条件,这一条件表明并非任何型值点组都存在保形插值.正因为如此,许多文献采用了不同的方法解决保形插值问题.[4—5]用分段有理三次,但计算量增加较大;[6]
1996年02期 104-113页 [查看摘要][在线阅读][下载 235k] - 宋叔尼
<正> 1 引言 关于Hammerstein型方程的数值逼近方法,许多作者做了工作,例如[1]、[2]、[3]、[4]等,他们把无限维空间中的 Hammerstein型方程转化为有限维空间中的非线性 Hammer-stein型方程,在此基础上,[1]、[2]又用Newton型迭代方法对有限维空间中的非线性方程做了进一步地讨论.[5]中把Newton迭代方法与投影方法结合在一起,考虑了Hilbert空间中具有紧性的非线性算子的不动点问题的数值解法.本文把Galerkin有限维逼近方法与Newton迭代方法紧密结合,把无限维Banach空间中一类具有单调型算子的非线性Ham-merstein型方程的求解问题在迭代过程中化为有限维空间中的线性代数方程组求解.并证明了迭代序列超线性收敛于原方程的解,最后举例说明了这一方法的应用.
1996年02期 114-121页 [查看摘要][在线阅读][下载 253k] - 石济民,林振宝,杨一都
<正> 1 引言 Richardson外推应用于椭圆偏微方程边值问题有限元法始于1978年(见[1],并于1983年在理论研究方面取得突破性进展(见[2]).自那以后有限元外推得到迅速发展,成为一个富于竞争的国际性研究课题(见[3],[4],[5]及其所列参考文献).但是通常的有限元外推需同时在每一个方向上分半加密网格,因此,对n维问题,细网格的结点数是粗网格的2~n倍,结果当n较大时(高维问题),细网格上的计算工作量十分庞大.为了克服这个缺点,发展了有限元单方向外推.对Poisson方程边值问题,[6]研究了2维矩形域上双线性有限元解的单方向外推,[7]研究了3维矩形域上三线性有限元单方向外推必须的插值渐近展开式,[8]研究了n维矩形域上n线性有限元解的区域分裂外推.本文旨在研究n维矩形域上Poisson方程边值问题及其对应的本征值问题n线性有限元解的单方向外推.始终假设本文出现的函数u是连续的.
1996年02期 122-134页 [查看摘要][在线阅读][下载 257k] - 顾传青
<正> 1 矩阵切触插值连分式 设实区间[a,b]中由不同点组成的插值结点为x_1,x_2,…,x_n,它们的重数分别为a_1,a_2,… ,a_n,M=sum from i=l to n(a_i-1),与之对应的待插值矩阵集为 {A_i~(k):k=0,1,…,a_i-1,i=1,2,…,n,A_i~(k)=A~(k)(x_i)∈R~(d×d)}. 设方阵A=(a_(ij)),它的广义矩阵逆定义为 A~(-1)= A/‖A‖~2 (A≠0) (1.1)
1996年02期 135-141页 [查看摘要][在线阅读][下载 150k] - 戴或虹,袁亚湘
<正> 1 引言 在求解无约束优化问题 min f(x) (1.1) x∈R_n的非线性共轭梯度法中,最早的一种由Flethcher和Reeves[6]在1964年提出,它具有如下形式:
1996年02期 142-148页 [查看摘要][在线阅读][下载 175k] - 汪静,张伟江
<正> 1 引言文章考虑如下形式的非线性奇异摄动差分方程:ε·{(g(y_(k+1))-2g(y_k)+g(y_(k-1)))/h~2}+(f(y_(k+1))-f y_k))(/h+H(x_k,y_k)=0, (1.1) y_0=α,y_N=β
1996年02期 149-156页 [查看摘要][在线阅读][下载 190k] - 刘国庆
<正> 1 引言 我们考虑拟线性奇异摄动Dirichlet问题 εy″-(f(y))′-b(x,y)=0,0<x<1, (1.1a) y(0)=A,y(1)=B (1.1b)其中ε是一个小的正参数.为了数值求解这类问题,Niijima导出了如下形式的差分格式 F_h(y)│i≡(ε/h~2)[intergral from n=y_i to y_(i+1)ξ(α(s)/2)ρ)ds-intergral from n=y_(i-1) to y_i (ξ(α(s)/2ρ)ds)]-(1/2h)(f(y_(i+1))-f(y(i-1))-b(x_i,y_i)=0,i=1,2,…,N-1, y_0=A,y_N=B (1.2)
1996年02期 157-164页 [查看摘要][在线阅读][下载 195k] - 孔敏,庄建南
<正> 1 引言 设F:ΩR~n→R,其中Ω是对n维欧氏空间中的紧集,F为非光滑函数.假定 F在Ω内部有极小点,我们的问题是考虑求解 minF(x) x∈Ω (1.1) 上述即是所谓的求解非光滑函数F总体极小点问题.目前尚未见到有关求解这类问题的总体极小点的理论和算法.葛人溥在讨论求解具有非线性约束、目标函数为光滑的
1996年02期 165-174页 [查看摘要][在线阅读][下载 286k] - 曾金平
<正> 1 引言 区域分解法和多重网格法都被认为是求解椭圆边值问题的快速算法.这两类算法也先 后应用于变分不等式的求解并获得了较为成功的数值尝试,收敛性理论也相继建 立.但是和用于方程问题不同,建立相应的h无关收敛性理论甚至更初步的收敛率分析遇到 一定的困难.九十年代初,Kornhuber针对变分不等式第一边值问题及摩擦问题进一步 讨论了多重网格法的收敛性质并在其离散问题非退化情形证明了渐近几何收敛速度,但仍 未见到有关h无关收敛性.区域分解法起步稍晚,但自八十年代末Lions给出了Schwarz交 替法的变分解释以来发展很快.Kuznetsov等人于九十年代初证明了乘性 Schwarz和加性Schwarz算法用于求解单边障碍问题时单调收敛于解.在同样条件下, [13]得到了误差估计式并利用无约束情形的有关结果得到了h无关收敛性.但是,在前述 的各种区域分解法中,子问题的求解都是精确的,因此在子域上费时较多而且在数值上也往 往只能得到子问题的近似解.这样自然产生这样一个想法:能否在子问题上和多重网格法 一样用近似解代替?本文即是针对此问题,从加性Schwarz算法入手,不仅证明算法收敛,而
1996年02期 175-182页 [查看摘要][在线阅读][下载 185k] - 肖爱国
<正> 1 散逸动力系统 考虑初值问题 y′(t)=f(y),y(0)=y_0∈R~N,t≥0, (1.1)这里f:R~N→R~N是满足局部Lipschitz条件的连续映射,并满足条件 Re〈u,f(u)〉≤α-β‖u‖~2 u∈R~N,(1.2)其中α≥0,β>0,〈·,·〉是R~N中标准内积,‖·‖是相应的内积范数.设y(t)是问题(1.1)-(1.2)的一个真解,则 ‖y(t)‖~2≤α/β+e~(-2βt)(‖y_0‖~2-α/β) t≥0, (1.3)及 ‖y(t)‖≤max(‖y_0‖,α/β) t≥0,(1.4)
1996年02期 183-189页 [查看摘要][在线阅读][下载 163k] - 马明书
<正> 1 引言 求解抛物型方程 u/t=u/x~2, 0<x<L,t>0, (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,对
1996年02期 190-193页 [查看摘要][在线阅读][下载 88k] 下载本期数据