中国科学院数学与系统科学研究院期刊网

数值计算与计算机应用 2008年 29卷

Please wait a minute...
选择: 显示/隐藏图片
1. 一种改进的4节点组合杂交轴对称元
尹云辉,聂玉峰
数值计算与计算机应用    2008, 29 (1): 1-7.   DOI: 10.12288/szjs.2008.1.1
摘要921)      PDF(pc) (305KB)(301)    收藏
在利用能量协调条件构造应力模式时,最后的应力模式中含有矩阵的逆.这个矩阵是否可逆,只有在具体的数值算例中验证.本文将"逆规则"应用于轴对称元应力模式的推导过程中,这样简化了求逆矩阵的形式,从而能够从理论上证明出这个矩阵在常见的几种网格剖分情况下是可逆的.经过这样的简化,最后得到的轴对称元仍具有较好的数值精度和性能,它不发生Poisson闭锁,对单元畸变影响的敏感度小.
参考文献 | 相关文章 | 多维度评价
被引次数: Baidu(2)
2. FDTD计算中一种UPML吸收边界与其内部计算区域的统一建模方法
吴泽艳,薛晓春
数值计算与计算机应用    2008, 29 (1): 8-14.   DOI: 10.12288/szjs.2008.1.8
摘要2108)      PDF(pc) (395KB)(317)    收藏
基于C++语言的多态性实现了单轴各向异性完全匹配层(LIPML)吸收边界与吸收边界内部计算区域的统一建模.其主要思想是:首先构造基类-Yee元胞类及其继承类来分别封装UPML内部介质和UPML的电磁特性;然后分别创建基于以上两个类的对象数组来给UPML及其内部计算区域开辟计算空间;再构造基类类型的指针数组,并用以上数组的地址赋值;最后,所有的计算在指针数组空间完成.该方法避免了UPML与其内部计算区域间的数据传递,简化了编程.数值实验验证了UPML的吸收效果,证明了方法的有效性.
参考文献 | 相关文章 | 多维度评价
被引次数: Baidu(1)
3. 电场积分方程的快速求解
王武,冯仰德,迟学斌
数值计算与计算机应用    2008, 29 (1): 15-24.   DOI: 10.12288/szjs.2008.1.15
摘要1166)      PDF(pc) (533KB)(344)    收藏
矩量法(MOM)离散电场积分方程(EFIE)得到稠密的线性方程组,它可以用迭代法(比如本文中的TFQMR方法)求解.每次迭代过程中,矩阵与向量的乘积的复杂度为O(N2).采用多层快速多极子方法(MLFMM),可将其降到O(NlogN).采用基于球谐变换的快速傅立叶变换,可进一步加快MLFMM的层间插值计算.数值结果显示MLFMM求解EFIE是可行的.
参考文献 | 相关文章 | 多维度评价
被引次数: Baidu(2)
4. 平面三向交错网格上Cauchy-Riemann方程的数值离散及快速解法
杨超,孙家昶
数值计算与计算机应用    2008, 29 (1): 25-38.   DOI: 10.12288/szjs.2008.1.25
摘要2150)      PDF(pc) (445KB)(283)    收藏
Cauchy-Riemann方程在复变函数、流体力学、偏微分方程组理论等方面具有重要的研究价值和应用背景.现有的关于Cauchy-Riemann方程快速数值解法的研究主要局限于张量积区域.本文利用平面上三向坐标改写了Cauchy-Riemann方程,并设计了一类三向交错网格以及相应的差分格式.本文证明了这种差分格式虽然只具有局部一阶的截断误差,但实际有二阶整体收敛性.最后还给出相应的谱预条件子快速解法和一些数值例子.
参考文献 | 相关文章 | 多维度评价
5. 分块三对角方程组的数值解法
刘长河,代西武,汪元伦
数值计算与计算机应用    2008, 29 (1): 39-48.   DOI: 10.12288/szjs.2008.1.39
摘要1120)      PDF(pc) (285KB)(252)    收藏
本文给出一种求解分块三对角线性方程组的一种新算法.它是一种直接解法,可以避免迭代法所产生的误差积累.由于充分利用了分块三对角矩阵的特点,这又是一种快速解法.经过严格的理论证明可知,本文算法是数值稳定的.数值实验验证了本文的算法具有很高的精确度.
参考文献 | 相关文章 | 多维度评价
被引次数: Baidu(10) CSCD(1)
6. 二维非结构网格浅水波方程的二阶精度非振荡有限体积方法
袁美,宋松和
数值计算与计算机应用    2008, 29 (1): 49-55.   DOI: 10.12288/szjs.2008.1.49
摘要833)      PDF(pc) (500KB)(250)    收藏
考虑浅水波方程,对二维非结构网格给出了一种非振荡有限体积方法.该方法的主要思想是在每一个三角形单元上采用最小二乘的思想构造一个重构函数,而时间离散采用二步TVD Runge- Kutta方法.最后用该格式对二维溃坝问题进行了数值试验,得到了满意的结果.
参考文献 | 相关文章 | 多维度评价
被引次数: Baidu(2) CSCD(1)
7. 对流扩散方程的间断Galerkin自适应方法
张强
数值计算与计算机应用    2008, 29 (1): 56-64.   DOI: 10.12288/szjs.2008.1.56
摘要998)      PDF(pc) (539KB)(268)    收藏
本文给出求解二维非线性对流扩散方程的局部间断Galerkin有限元方法在非协调三角网格上的自适应算法实现.数值算例表明这种方法可以高效追踪真解的剧烈变化.
参考文献 | 相关文章 | 多维度评价
被引次数: Baidu(21)
8. NAPA软件的并行优化
金君,乔楠,梁德旺
数值计算与计算机应用    2008, 29 (1): 65-72.   DOI: 10.12288/szjs.2008.1.65
摘要2167)      PDF(pc) (535KB)(282)    收藏
大规模数值计算受到通信模式、并行算法、I/O速度等的多方面因素的制约,并行程序的好坏直接影响并行机性能的发挥,本文分别对上述影响并行性能的重要因素进行了分析并对NAPA软件进行了优化,测试中发现本文采用的并行算法性能比优化前提高了41.1%,此外,本文采用支持多视口的MPI I/O接口性能有明显提高.最后,本文分析了并行NAPA软件的可扩展性,并采用高超声速平板流动进行了测试,在Grid 97*49*49算例中,64个进程的情况下得到了较高的加速比(53.7)和并行效率(84%),表明,优化后的软件具有较好的并行效率和可扩展性.
参考文献 | 相关文章 | 多维度评价
被引次数: Baidu(2) CSCD(1)
9. 非线性随机延迟微分方程半隐式Euler方法的均方稳定性
王文强,黄山,李寿佛
数值计算与计算机应用    2008, 29 (1): 73-80.   DOI: 10.12288/szjs.2008.1.73
摘要1056)      PDF(pc) (420KB)(300)    收藏
本文首先将数值方法的均方稳定性的概念MS-稳定与GMS-稳定推广到一般情形,然后针对一维情形下的非线性随机延迟微分方程初值问题,证明了如果问题本身满足零解是均方渐近稳定的充分条件,那么当漂移项满足一定的限制条件时,半隐式Euler方法是MS-稳定的且带线性插值的半隐式Euler方法是GMS-稳定的理论结果.
参考文献 | 相关文章 | 多维度评价
被引次数: Baidu(22) CSCD(3)
10. 低阶实时最优Runge-Kutta算法
黄振全,陈志武,何云
数值计算与计算机应用    2008, 29 (2): 81-88.   DOI: 10.12288/szjs.2008.2.81
摘要1449)      收藏
在控制系统实时Runge-Kutta算法中,为了满足实时仿真快速性需求,希望尽可能地采用大的计算步长.如果采用大步长,那么数值计算就会引起数值不稳定或者计算误差太大的问题.在现有低阶实时龙格-库塔公式基础上,首先利用RK公式的稳定性方程求解出最大稳定域,然后根据截断误差与相关系数的关系,将其化为一个约束求极小最优问题,并最终推导出实时最优三级二阶RK公式和四级三阶RK公式.仿真结果表明,该算法具有一定的优越性.
参考文献 | 相关文章 | 多维度评价
被引次数: Baidu(1)
11. 一类舒尔补矩阵的条件数分析
钟尔杰,黄廷祝
数值计算与计算机应用    2008, 29 (2): 89-95.   DOI: 10.12288/szjs.2008.2.89
摘要1241)      PDF(pc) (266KB)(287)    收藏

研究了一类线性方程组系数矩阵的红黑排序方法,以及由红黑排序矩阵导出的舒尔补矩阵的条件数.利用三对角矩阵的特征值分析方法推导了一类块三对角矩阵的特征值和条件数,构造了三对角矩阵和块三对角矩阵的红黑排序排列矩阵,利用矩阵相似变换推导出红黑排序矩阵中的舒尔补的特征值和条件数表达式.理论分析和数值试验结果均表明这类舒尔补矩阵具有更好的性质.

参考文献 | 相关文章 | 多维度评价
被引次数: Baidu(4)
12. 解抛物型方程的Galerkin多层修正迭代算法
何崇南,邓小炎
数值计算与计算机应用    2008, 29 (2): 96-104.   DOI: 10.12288/szjs.2008.2.96
摘要1123)      PDF(pc) (356KB)(268)    收藏

小波方法在微分方程数值解法中日益得到广泛应用.由于小波的紧支性、正交性使得离散后的代数方程组的系数矩阵具有稀疏性、层次性,在此基础上可以构造各种快速算法.基于多尺度空间,采用一组正交小波基来离散原方程,导出方程组的系数矩阵具有稀疏性和层次性,从而提出求抛物型微分方程的Galerkin多层修正迭代算法,并讨论了迭代修正算法的收敛性.提出的方案能容易地实现时间和空间方向的局部加密自适应修正过程.提供的数值算例说明了方法的有效性.

参考文献 | 相关文章 | 多维度评价
被引次数: Baidu(6)
13. 一种基于核映射的自适应优化配置属性权重组的方法
陈新泉
数值计算与计算机应用    2008, 29 (2): 105-118.   DOI: 10.12288/szjs.2008.2.105
摘要1952)      PDF(pc) (484KB)(328)    收藏

通过引入一个从样本空间到特征空间的核映射,从而将样本空间中的分类问题与特征空间中的聚类问题联系起来.为获得样本空间中的一组合适的属性权重值,提出了一种基于核映射的自适应优化配置属性权重组的方法.在特征空间中根据"聚类之内的数据点最大限度的相近,聚类之间的数据点最大限度的相离"这个原则,提出了一个带约束的混和目标函数,通过优化这个混和目标函数来获得样本空间中的一个合适的属性权重组.为求解这个混和目标函数,提出了一种基于负投影梯度的自适应优化配置属性权重组的方法.接着采用UCI的两个标准数据集来进行实验验证,可以证实这种根据给定数据点集进行自适应优化配置样本空间中的属性权重组方法的有效性.最后给出了两种自适应优化目标函数权重参数和核函数参数的方法.

参考文献 | 相关文章 | 多维度评价
被引次数: Baidu(8) CSCD(1)
14. 基于混沌微粒群优化算法的阈值图像分割
黄力明
数值计算与计算机应用    2008, 29 (2): 119-125.   DOI: 10.12288/szjs.2008.2.119
摘要1269)      PDF(pc) (444KB)(522)    收藏

混沌微粒群优化算法利用了粒子群优化算法收敛速度快和混沌运动所具有的随机性、遍历性和初值敏感性,将混沌状态引入到优化变量中,把混沌的遍历范围映射到优化变量的取值范围.在算法执行过程中对优秀个体混沌扰动,有利于跳出局部极值点,搜索到全局最优解.分别用微粒群优化算法和混沌微粒群优化算法求解函数优化问题,对算法的性能进行检验,检验结果显示:混沌微粒群优化算法搜索全局最优解的成功率和收敛速度都要优于微粒群优化算法.将混沌微粒群优化算法与阈值法相结合,在算法初始化阶段对粒子位置混沌初始化;在算法运行期间对优秀个体进行混沌扰动避免落入局部最优,较好地解决了传统的多阈值图像分割方法中运算量大的问题.实验结果表明,混沌微粒群优化算法用于阈值寻优减少了搜索时间,提高了收敛率.

参考文献 | 相关文章 | 多维度评价
被引次数: Baidu(5) CSCD(1)
15. 求矩阵方程AXB=C的双对称最小二乘解的迭代算法
尚丽娜,张凯院,陈梅枝
数值计算与计算机应用    2008, 29 (2): 126-135.   DOI: 10.12288/szjs.2008.2.126
摘要2217)      PDF(pc) (291KB)(230)    收藏

基于求解线性代数方程组的共轭梯度法的思想,通过特殊的变形与近似处理,建立了求矩阵方程AXB=C的双对称最小二乘解的迭代算法,并证明了迭代算法的收敛性.不考虑舍入误差时,迭代算法能够在有限步计算之后得到矩阵方程的双对称最小二乘解;选取特殊的初始矩阵时,还能够求得矩阵方程的极小范数双对称最小二乘解.同时,也能够给出指定矩阵的最佳逼近双对称矩阵.算例表明,迭代算法是有效的.

参考文献 | 相关文章 | 多维度评价
被引次数: Baidu(23) CSCD(5)
16. 非线性中立型延迟微分方程对角分裂Runge-Kutta法的收缩性
王晚生,苏凯,李寿佛
数值计算与计算机应用    2008, 29 (2): 136-145.   DOI: 10.12288/szjs.2008.2.136
摘要993)      PDF(pc) (339KB)(283)    收藏

考虑了对角分裂Runge-Kutta法求解非线性变延迟中立型微分方程的收缩性。证明了在最大范数下这类方法能够保持非线性中立型延迟微分方程系统的收缩性。数值试验验证了上述理论结果。

参考文献 | 相关文章 | 多维度评价
被引次数: Baidu(6)
17. 一维抛物型方程的样条子域精细积分(SSPI)隐格式
刘利斌,刘焕文
数值计算与计算机应用    2008, 29 (2): 146-152.   DOI: 10.12288/szjs.2008.2.146
摘要1920)      PDF(pc) (334KB)(264)    收藏

对一维抛物型方程初边值问题的求解,以往已经有一些数值解法,它们或者无条件稳定但精度不高,或者精度高但仅为条件稳定,且稳定性条件严格.另外,以往的差分格式在处理第二、第三类边界条件问题时,对带导数边界条件都是进行简单的差分逼近,影响了数值解的精度.因此构造一个无条件稳定且对各类边值问题都具有良好精度的数值方法具有重要意义.为此,基于子域精细积分思想,结合三次样条函数,提出了求解一维抛物型方程初边值问题含参数的样条子域精细积分格式.该格式为绝对稳定且精度很高.由于三次样条函数的采用,避免了通常有限差分法中处理带导数边界条件时产生的逼近误差,大大提高了求解第二、三类边界条件问题时的精度.

参考文献 | 相关文章 | 多维度评价
被引次数: Baidu(8) CSCD(1)
18. 比例方程组与MPI同步通信静态死锁检测
廖名学,范植华,何晓新
数值计算与计算机应用    2008, 29 (2): 153-160.   DOI: 10.12288/szjs.2008.2.153
摘要1005)      PDF(pc) (477KB)(224)    收藏

介绍了一种MPI程序死锁检测的静态方法以及该方法所处理的程序模型.为实现该方法,提出了比例方程组(一种特殊线性方程组)的概念并设计了求解方程组最简解的线性时空复杂度的高效算法.算法由一个四遍扫描过程与一个主控程序构成.主控程序用来处理并行计算节点计算机构成的划分.四遍扫描过程采用深度优先搜索方法确定方程组中各变元之间的比例关系.通过该算法所获得的最简解,任意多个变元之间的比例关系能在常数时间内获得.证明了该算法的正确性,并采用Java语言实现了该算法的标准程序库.该程序库目前已运行于MPI同步通信静态死锁检测的软件框架中.

参考文献 | 相关文章 | 多维度评价
被引次数: Baidu(2)
19. 薄壁管抗撞性能的多目标结构优化
陈仙燕, 龙述尧,杨旭静
数值计算与计算机应用    2008, 29 (3): 161-170.   DOI: 10.12288/szjs.2008.3.161
摘要1098)      PDF(pc) (771KB)(333)    收藏

以交通工具中部分锥形薄壁方管的安全装置作为研究对象, 建立以薄壁管在碰撞过程中吸收能量最大化,比吸能最大化和初始碰撞力峰值最小化为多目 标的优化问题.用锥形部分的几何参数作为设计变量, 在保证不降低薄壁管吸能能力的情况下, 通过对其结构的优化达到初始碰撞力峰值 最小化的目的.论文采用有限元软件LS-DYNA得到不同几 何参数模型的碰撞信息,用响应面法构造近似函数, 同时引入权系数以表征各个目标在优化问题中的重要程度, 并采用理想点法求解多目标优化问题,分析了锥形薄壁方管各几何参数对结构的能量吸收、比吸能和初始碰撞力峰值的影响,最终得到 了给定权系数下的最优模型.

参考文献 | 相关文章 | 多维度评价
被引次数: Baidu(27)
20. Newton型方法的推广与改进
司智勇, 阿不都热西提·阿不都外力, 张知难
数值计算与计算机应用    2008, 29 (3): 171-175.   DOI: 10.12288/szjs.2008.3.171
摘要997)      PDF(pc) (262KB)(401)    收藏

对解非线性方程组Newton迭代格式进行了改进, 得到了两种比Newton法较为宽松的并且 收敛速度较快的新的迭代格式.从而构造了两种新的Newton型迭代法.理论分析和数值实验证明这两种方法是稳定且有效的.

参考文献 | 相关文章 | 多维度评价
被引次数: Baidu(12) CSCD(2)
21. 收敛的全部线性多步法基本公式的推导
杨大地, 刘冬兵
数值计算与计算机应用    2008, 29 (3): 176-185.   DOI: 10.12288/szjs.2008.3.176
摘要1145)      PDF(pc) (355KB)(328)    收藏

本文首先定义了线性多步法基本公式的概念, 并应用MATLAB的符号运算, 推导了求解常微分方程初值问题的2-7步法全部基本公式, 其系数全部为精确的分数形式; 利用关于线性多步法公式的收敛条件, 筛选出其中收敛的公式, 计算出了公式的误差主项系数, 阶数,绝对稳定区间.

参考文献 | 相关文章 | 多维度评价
被引次数: Baidu(13) CSCD(4)
22. 具有细观结构新材料跨尺度分析的细观元法
曹志远,程红梅
数值计算与计算机应用    2008, 29 (3): 186-196.   DOI: 10.12288/szjs.2008.3.186
摘要948)      PDF(pc) (738KB)(304)    收藏

建议一种适用于具有细观结构新材料宏细观间跨尺度分析的细观元方法.细观元法 在结构的常规有限元内部设置密集细观单元以反映材料细观构造, 又通过协调条件 将各细观元结点自由度转换为同一常规有限元自由度, 再上机计算. 此方法可实现 材料细观结构到构件宏观响应的直接过渡分析,而计算单元与自由度又等同一般常 规有限元,为解决具有细观结构新材料与构件跨尺度分析提供一种新的有力工具.研究了直接从制备时给定的材料组分分布及网状细观结构图形出发计算功能梯度板件宏观响应,给出了不同边界条件功能梯度板件的力学量三维分布形态以及细观结构局部突变引起宏观等应力线图的畸变.

参考文献 | 相关文章 | 多维度评价
被引次数: Baidu(16) CSCD(2)
23. Hermite Toeplitz 类系统的预因子
彭小飞,黎稳
数值计算与计算机应用    2008, 29 (3): 197-206.   DOI: 10.12288/szjs.2008.3.197
摘要1209)      PDF(pc) (371KB)(317)    收藏

该文针对预处理共轭梯度法求解Hermite正定的Toeplitz线性方程组, 提出了新的含参预因子.该预因子可视为对T.Chan预因子的推广,但又不同于M.Cai 等推广的预因子. 我们分析了含参预因子的谱性 质及其求逆的运算量. 指出:在一定条件下,经过简单的参数选择,与该预因子对应 的预处理共轭梯度法优于上述提及的另两种预因子的 对应方法. 数值实例亦验证新的预因子在一定条件下的优越性.

参考文献 | 相关文章 | 多维度评价
24. Toeplitz矩阵相乘的一种新快速算法
余品能, 王煜
数值计算与计算机应用    2008, 29 (3): 207-216.   DOI: 10.12288/szjs.2008.3.207
摘要1325)      PDF(pc) (284KB)(302)    收藏

将Toeplitz矩阵分解为一个循环矩阵和一个下三角Toeplitz矩阵之和, 以及一般卷积向循环卷积的转化, 借助快速Fouier变换(FFT), 导出了一种计算两个$n$阶Toeplitz矩阵乘积的新快速算法, 其算法复杂性为$2n^2+{63\over 4}n\log_2n-15n-34$次实乘运算, $4n^2+{63\over 2}n\log_2n-18n+23$
次实加运算, 与已有的优化算法相比, 在实乘次数有所降低的同时, 实加次数降低了近${1\over3}$, 是目前复杂性最小的一种算法.

参考文献 | 相关文章 | 多维度评价
被引次数: Baidu(15)
25. 一类变时滞微分代数方程单支方法的收敛性
肖飞雁, 张诚坚
数值计算与计算机应用    2008, 29 (3): 217-225.   DOI: 10.12288/szjs.2008.3.217
摘要1064)      PDF(pc) (305KB)(346)    收藏

B-收敛和D-收敛的概念被推广到了变时滞微分代数方程问题,给出了$D_A$-收敛的定义,讨论了该类问题的$D_A$-收敛性,并给出了相应的误差估计,证明了如果G-稳定的单支方法对于常微分方程初值问题在经典意义下是p阶相容的且$\frac{\beta _k}{\alpha_k}>0$,那么具有线性插值过程的该方法是p阶$D_A$-收敛的,这里p=1或2.

参考文献 | 相关文章 | 多维度评价
被引次数: Baidu(5) CSCD(1)
26. 多层次图排序算法及其在图剖分中的应用
刘旭, 莫则尧
数值计算与计算机应用    2008, 29 (3): 226-240.   DOI: 10.12288/szjs.2008.3.226
摘要1057)      PDF(pc) (579KB)(329)    收藏

图排序问题在众多领域中有着重要应用. 本文利用多层次思想, 提出一种具有V-循环结构的新算法. 该算法是一种线性时间复杂度的方法. 在文中的4个算例中, 这种多层次方法所得到的排序质量至少比谱方法高5\%. 本文把它应用到图剖分领域, 利用KL/FM方法对其进行了局部修改, 得到了两种新 的图剖分算法. 在文中的4个算例中, 这两种方法都能提供与当前质量最佳算法相当的图剖分结果.

参考文献 | 相关文章 | 多维度评价
被引次数: Baidu(5) CSCD(1)
27. 抛物问题Mortar有限元的瀑布型多重网格法
周叔子,史艳华
数值计算与计算机应用    2008, 29 (4): 241-250.   DOI: 10.12288/szjs.2008.4.241
摘要1211)      PDF(pc) (285KB)(265)    收藏

对抛物问题的全离散格式采用Mortar型有限元逼近, 构造了相应的瀑布型多重网格法,证明了该方法的最优性.

参考文献 | 相关文章 | 多维度评价
被引次数: Baidu(7) CSCD(1)
28. 椭圆函数的精细积分改进算法
姚征, 钟万勰
数值计算与计算机应用    2008, 29 (4): 251-260.   DOI: 10.12288/szjs.2008.4.251
摘要2236)      PDF(pc) (451KB)(222)    收藏

椭圆函数是一种特殊的双周期复变函数, 广泛应用于工程问题中, 尤其非线性问题中居多. 在工程中遇到的椭圆函数以二阶椭圆函数为主, 而且很多复杂的椭圆函数都可以通过变换由二阶椭圆函数得到. 二阶椭圆函数包括Jacobi椭圆函数和Weierstrass椭圆函数. 它们都可以进行幂级数展开, 直接计算很不方便. 椭圆函数的重要性质之一就是具有加法定理, 因此可利用精细积分法求解. 虽然椭圆函数的精细积分算法在精度和效率上取得了较大成功, 但椭圆函数的奇点问题仍然存在并对计算精度构成一定威胁. 在回顾并分析椭圆函数的精细积分算法的基础上,通过对椭圆函数奇点的分析, 给出了椭圆函数可去奇点的近似公式, 并在此基础上进一步改进并完善了椭圆函数的精细积分算法.

参考文献 | 相关文章 | 多维度评价
被引次数: Baidu(29) CSCD(6)
29. 解扩散方程的指数时间差分方法
孙建强, 秦孟兆, 戴桂冬
数值计算与计算机应用    2008, 29 (4): 261-266.   DOI: 10.12288/szjs.2008.4.261
摘要2169)      PDF(pc) (308KB)(460)    收藏

指数时间差分方法是近年来提出求解刚性常微分方程的一种新的数值计算方法. 指数时间差分方法是一种积分方法,而不是经典的差分方法.
利用指数时间差分方法求解扩散方程,如一维拟线性对流扩散方程和Allen-Cahn扩散方程. 扩散方程在空间方向离散后转化成刚性常微分方程. 用显式指数时间差分方法和相应阶的显式 Runge-Kutta方法求解刚性常微分方程. 数值结果表明显式指数时间差分方法具有相同阶的显
式Runge-Kutta方法相应的精度,稳定性显著提高,而且能很好地模拟扩散方程的演化行为. 指数时间差分方法可用于刚性常微分方程的数值计算.

参考文献 | 相关文章 | 多维度评价
被引次数: Baidu(3) CSCD(1)
30. 二次插值模型直接搜索算法的参数分析
刘浩, 倪勤
数值计算与计算机应用    2008, 29 (4): 266-276.   DOI: 10.12288/szjs.2008.4.266
摘要2136)      PDF(pc) (419KB)(303)    收藏

对无约束优化问题的二次插值型直接搜索算法中初始插值半径,信赖域初始半径, 位移接受准则和信赖域半径调节参数进行了数值实验分析.
数值实验表明解无约束优化的基于二次函数插值型的直接搜索算法对初始插值半径和信赖域初始半径比较敏感,对位移接受准则和半径调节参数不敏感. 根据数值实验结果推荐初始插值半径的选取应与信赖域初始半径相等,同时 给出了基于二次插值型的直接搜索算法中初始插值半径与信赖域初始半径的选择区间和其它参数的推荐值. 这些结果对这类算法的数值实现和工程应用是有益的.

参考文献 | 相关文章 | 多维度评价
31. 三维椭圆问题三次有限元方程的代数多层网格法
张红梅, 肖映雄, 舒适
数值计算与计算机应用    2008, 29 (4): 277-290.   DOI: 10.12288/szjs.2008.4.277
摘要1466)      PDF(pc) (465KB)(229)    收藏

通过分析三次有限元空间与线性有限元空间之间的关系,提出了一种求解三维椭圆问题三次有限元方程的两水平方法. 然后,通过调用现有的代数多层网格 (AMG)法求解粗水平方程,建立了求解三次有限元方程的AMG法,并对其收敛性进行了严格的理论分析. 数值实验结果表明,本文设计的AMG方法对求解三维椭圆问题三次有限元方程具有很好的计算效率和鲁棒性.

参考文献 | 相关文章 | 多维度评价
被引次数: Baidu(5) CSCD(3)
32. 求解周期三对角Toeplitz方程组的一种新修正算法
徐仲,安晓虹,陆全
数值计算与计算机应用    2008, 29 (4): 291-295.   DOI: 10.12288/szjs.2008.4.291
摘要1351)      PDF(pc) (222KB)(261)    收藏

利用Sherman-Morrison-Woodbury公式导出了求解周期三对角Toeplitz方程组的一 种新的修正算法.该算法的计算量比求解周期三对角方程组的追赶法要少,且可以 并行计算.对新算法进行了可行性和稳定性分析.数值算例表明了新算法的有效性

参考文献 | 相关文章 | 多维度评价
33. Helmholtz 方程内问题的数值方法
翟亮亮, 王连堂, 王俊杰, 周斌
数值计算与计算机应用    2008, 29 (4): 296-301.   DOI: 10.12288/szjs.2008.4.296
摘要1616)      PDF(pc) (247KB)(295)    收藏

应用位势理论把Helmholtz 方程内问题转化为含有 Cauchy 奇性的第二类积分方程的求解问题.并利用 Nystom 方法求得数值结果, 试验结果表明了此方法的简单与有效性

参考文献 | 相关文章 | 多维度评价
被引次数: Baidu(4) CSCD(3)
34. 动态中子输运方程迭代初值的选取方法
洪振英,袁光伟,傅学东
数值计算与计算机应用    2008, 29 (4): 302-312.   DOI: 10.12288/szjs.2008.4.302
摘要1384)      PDF(pc) (664KB)(366)    收藏

研究离散纵标动态中子输运方程迭代求解时,迭代初值的不同选取方法,设计合理的迭代初值可以适当放宽对时间步长的限制,缩短计算时间. 设计四种迭代初值并应用于数值求解中的等比格式和菱形格式,其中等比格式形成非线性离散方程,菱形格式形成线性离散方程.考察不同迭代初值的计算效率,分别对物理量变化平缓以及变化剧烈的问题进行考察.数值算例表明构造的基于物理量随时间走势的预估值作为迭代初值优势明显,这在保证计算精度的前提下提高了数值计算效率.

参考文献 | 相关文章 | 多维度评价
被引次数: Baidu(1)
35. 二次特征值问题的并行Jacobi-Davidson方法及其应用
王顺绪, 戴华
数值计算与计算机应用    2008, 29 (4): 313-320.   DOI: 10.12288/szjs.2008.4.313
摘要2160)      PDF(pc) (363KB)(335)    收藏

给出了求解二次特征值问题多个特征对的一种并行Jacobi-Davidson方法, 该方法在子空间中求解投影矩阵的二次特征值问题,利用校正方程的解扩充子空间, 并以某型号机翼在结构动力分析中的二次特征值问题为例, 在多处理机并行系统IBM-P650上进行了数值试验,试验结果表明该算法具有较高的加速比和并行效率.

参考文献 | 相关文章 | 多维度评价
被引次数: Baidu(6) CSCD(1)