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

数值计算与计算机应用 2013年 34卷

Please wait a minute...
选择: 显示/隐藏图片
1. Wangerin函数Smn(μ)的数值计算及数值可视化分析
宋艳萍, 黄华, 恰汗·合孜尔
数值计算与计算机应用    2013, 34 (1): 1-8.   DOI: 10.12288/szjs.2013.1.1
摘要1061)      PDF(pc) (385KB)(442)    收藏
结合第一类完全椭圆积分的迭代算法和Jacobi椭圆函数的恒等式, 实现Wangerin函数Smn(μ)的高精度数值计算. 进一步, 利用数值可视化方法对它的数值计算结果进行分析, 总结出Wangerin函数Smn(μ)的曲线变化特征. 这为深入研究Wangerin函数的收敛性和稳定性奠定理论基础, 具有一定的工程实际价值.
参考文献 | 相关文章 | 多维度评价
2. 多矩阵变量线性矩阵方程的广义自反解的迭代算法
王娇, 张凯院, 李书连
数值计算与计算机应用    2013, 34 (1): 9-19.   DOI: 10.12288/szjs.2013.1.9
摘要1477)      PDF(pc) (387KB)(746)    收藏
基于求线性矩阵方程约束解的修正共轭梯度法的思想方法, 通过修改某些矩阵的结构,建立了求特殊类型的多矩阵变量线性矩阵方程的广义自反解的迭代算法, 证明了迭代算法的收敛性, 解决了给定矩阵在该矩阵方程的广义自反解集合中的最佳逼近计算问题.当矩阵方程相容时, 该算法可以在有限步计算后得到其一组广义自反解; 选取特殊的初始矩阵, 能够求得其极小范数广义自反解. 数值算例表明, 迭代算法是有效的.
参考文献 | 相关文章 | 多维度评价
被引次数: Baidu(5)
3. 基于辛格式的谱元法及其在横向各向同性介质波场模拟中的应用
汪文帅, 李小凡
数值计算与计算机应用    2013, 34 (1): 20-30.   DOI: 10.12288/szjs.2013.1.20
摘要1365)      PDF(pc) (644KB)(683)    收藏
基于分部的Runge-Kutta离散形式, 给出了一种新的三阶辛积分算法, 数值试验表明,长时程计算时该算法具有好的控制误差累积的能力; 与有限差分法进行空间域离散相结合, 通过数值试验进一步说明算法的有效性. 注意到位移波动方程通过谱元离散后的微分方程组, 完全符合新推导的三阶辛算法离散所需形式, 因此将该三阶辛算法与谱元法结合具有很好的优势, 并通过对横向各向同性介质弹性波场的模拟, 结果显示不但成功模拟了波的传播特性, 而且相对于传统算法, 优势明显.
参考文献 | 相关文章 | 多维度评价
4. 填充介质倒梯形脊波导衰减特性的有限元分析
陈小强, 李建兵, 高原, 李向葵
数值计算与计算机应用    2013, 34 (1): 31-37.   DOI: 10.12288/szjs.2013.1.31
摘要1058)      PDF(pc) (421KB)(433)    收藏
脊波导在传播电磁波的过程中伴随有能量损耗, 为了研究填充介质脊波导的衰减特性, 本文采用有限元法借助Matlab软件计算了TE模式下部分填充介质倒梯形对称双脊波导的衰减常数和功率容量. 分析了脊波导衰减常数和功率容量与脊波导尺寸、工作频率、填充介质相对介电常数之间的关系, 给出关系曲线图.这些结果将丰富现有的填充非均匀介质脊波导的数据, 为填充介质脊波导的设计和优化提供有价值的依据.
参考文献 | 相关文章 | 多维度评价
被引次数: Baidu(1)
5. 高维0-1瓶颈问题的动态规划算法
罗宗俊
数值计算与计算机应用    2013, 34 (1): 38-46.   DOI: 10.12288/szjs.2013.1.38
摘要985)      PDF(pc) (367KB)(383)    收藏
在本文中, 我们通过一个实际问题归纳出一个数学模型(正文中的模型Ⅰ),并通过新变量的引用, 将模型Ⅰ转化成一个高维0-1瓶颈规划(正文中的模型Ⅱ).对模型Ⅱ, 我们建立了求模型Ⅱ最优解的动态规划算法(带有阀值Q). 该算法与普通动态规划相比大大节约了运算量. 最后指出了该算法对0-1瓶颈问题的求解具有一定的普遍性.
参考文献 | 相关文章 | 多维度评价
6. 求解大型Stein方程的块Krylov子空间方法
黄飞虎, 汪晓虹
数值计算与计算机应用    2013, 34 (1): 47-58.   DOI: 10.12288/szjs.2013.1.47
摘要1102)      PDF(pc) (356KB)(411)    收藏
本文研究利用块Krylov子空间方法对大型Stein方程降阶求解, 分别基于块Arnoldi方法与非对称块Lanczos方法,提出了块Arnoldi Stein方法与非对称块Lanczos Stein方法.数值实验表明提出的方法有效.
参考文献 | 相关文章 | 多维度评价
7. 一维抛物型方程第三边值问题的紧有限体积格式
王风娟, 王同科
数值计算与计算机应用    2013, 34 (1): 59-74.   DOI: 10.12288/szjs.2013.1.59
摘要1099)      PDF(pc) (441KB)(279)    收藏
本文针对一维抛物型方程第三边值问题提出了一种紧有限体积格式, 该格式形成的线性代数方程组具有对称三对角性质, 且不可约占优, 可以使用追赶法求解.证明了格式按照离散L2范数在空间方向具有3.5阶精度, 在时间方向具有2阶精度. 数值算例验证了理论分析的正确性, 并说明了格式的有效性.
参考文献 | 相关文章 | 多维度评价
被引次数: Baidu(2) CSCD(1)
8. 求解多层介质中声波传播问题的一种边界元方法
岳雪莲, 王连堂, 孟文辉
数值计算与计算机应用    2013, 34 (1): 75-80.   DOI: 10.12288/szjs.2013.1.75
摘要955)      PDF(pc) (346KB)(399)    收藏
针对多层介质中声波的传播问题, 将其中偶数(或奇数)层内的声波用一种单双层混合位势的形式来表示, 再应用Green定理表示出其余层的声波并形成相应的边界积分方程.如果区域有M层时, 传统的边界元方法最终将形成$2M$个边界积分方程并对应2M个未知函数, 而应用上述方法求解该问题时, 最终只形成M个边界积分方程以及对应M个未知函数, 从而使得求解的方程和未知数的个数都减少了一倍.最后, 通过对数值算例的求解, 验证了该方法的可行性及精确性.
参考文献 | 相关文章 | 多维度评价
被引次数: Baidu(7)
9. Chebyshev谱-Euler混合方法求解一类非线性Burgers方程
闵涛, 任菊成, 耿蓓
数值计算与计算机应用    2013, 34 (2): 81-88.   DOI: 10.12288/szjs.2013.2.81
摘要1037)      PDF(pc) (588KB)(362)    收藏
将Chebyshev谱方法与Euler方法相结合, 对一类非线性Burgers方程进行数值求解, 通过数值模拟将其与有限差分法和粒子无网格线混合格式MPS-MAFL方法进行了比较, 结果表明这种方法对于求解非线性Burgers方程具有较好的效果.
参考文献 | 相关文章 | 多维度评价
被引次数: Baidu(4)
10. 独立成分分析方法在指纹提取中的应用
肖世校
数值计算与计算机应用    2013, 34 (2): 89-94.   DOI: 10.12288/szjs.2013.2.89
摘要809)      PDF(pc) (493KB)(335)    收藏
为了提高指纹图像识别的准确性和可靠性, 先对指纹图像采用主成分分析(PCA)算法进行预处理, 去除指纹图像的二阶统计特征相关. 用ICA算法对加噪的指纹图像进行提取识别, 并利用统计方法确定混合矩阵, 实验表明: 在一定的条件下, ICA能够很好地分离混合图像, 提高准确性和可靠性.
参考文献 | 相关文章 | 多维度评价
被引次数: Baidu(4) CSCD(1)
11. 求解周期结构中声波传播问题的一种边界元法
杨正辉, 王连堂, 孟文辉
数值计算与计算机应用    2013, 34 (2): 95-104.   DOI: 10.12288/szjs.2013.2.95
摘要812)      PDF(pc) (533KB)(376)    收藏
本文应用边界积分方程方法(边界元法)求解无限长周期结构中声波的传播问题, 根据散射体的周期排布方式, 将无限个散射体对应无限个边界积分方程的求解问题转化为在某一个单位块中有限个边界积分方程的求解, 从而使得该问题的数值求解变得简单可行. 然后将该方法应用于声子晶体能量禁带的预测中, 通过对数值算例的求解以及与其他方法求解结果的比较, 验证了本文所提出方法的可行性和准确性.
参考文献 | 相关文章 | 多维度评价
被引次数: Baidu(1)
12. 基于混合蚁群算法的DNA编码序列设计方法
郭稳涛, 何怡刚
数值计算与计算机应用    2013, 34 (2): 105-116.   DOI: 10.12288/szjs.2013.2.105
摘要833)      PDF(pc) (644KB)(322)    收藏
通过将遗传算法中的交叉、变异操作与蚁群算法中的协同模型进行结合,提出了一种基于混合蚁群算法的DNA编码序列设计方案. 实验表明,该算法具有较高的收敛速度, 能为DNA计算提供可靠的编码序列.
参考文献 | 相关文章 | 多维度评价
被引次数: Baidu(2)
13. 广义Nekrasov矩阵的新迭代判别法
石玲玲, 徐仲, 陆全, 周伟伟
数值计算与计算机应用    2013, 34 (2): 117-122.   DOI: 10.12288/szjs.2013.2.117
摘要1009)      PDF(pc) (304KB)(280)    收藏
广义Nekrasov矩阵作为一类特殊的广义严格对角占优矩阵在科学和工程实际中有着广泛的应用, 因此研究这类矩阵的判定问题是非常重要的. 给出了判定一个矩阵是否为广义Nekrasov矩阵的两种新的迭代算法, 并用数值算例说明了算法的有效性. 由于证明了广义Nekrasov矩阵就是广义严格对角占优矩阵, 从而也就得到了两种新的判定广义严格对角占优矩阵的迭代算法.
参考文献 | 相关文章 | 多维度评价
被引次数: CSCD(1)
14. 角锥棱镜误差检测数学建模分析
李福先
数值计算与计算机应用    2013, 34 (2): 123-135.   DOI: 10.12288/szjs.2013.2.123
摘要801)      PDF(pc) (478KB)(382)    收藏
通过研究角锥棱镜的结构和光学特征, 重点分析了角锥棱镜误差检测数学模型的建立过程, 并利用合理的模型假设、近似处理和数学推导进行模型求解, 计算出角度误差和面型误差, 求出角锥棱镜出射光线与入射光线之间的夹角和角锥棱镜各直角面之间夹角的关系, 完成了模型分析. 经实际生产检验, 角锥棱镜误差检测数学模型与实际吻合.
参考文献 | 相关文章 | 多维度评价
被引次数: Baidu(4) CSCD(1)
15. 一类大规模稀疏矩阵特征问题求解的并行算法
吴洋, 赵永华, 纪国良
数值计算与计算机应用    2013, 34 (2): 136-146.   DOI: 10.12288/szjs.2013.2.136
摘要1408)      PDF(pc) (915KB)(634)    收藏
本文提出一种求解大规模稀疏矩阵特征问题的并行共轭梯度算法.为了提高算法的并行效率, 设计了负载平衡的行划分方式, 实现了计算和通信重叠的稀疏矩阵重排序方法, 通过预处理减少计算过程中各进程间消息传递的通信量. 另外, 基于多核处理器高性能并行计算, 实现了 MPI和细粒度(线程级)OpenMP混合并行算法. 在深腾7800并行计算机上对并行算法进行了测试, 结果表明在进程数增多时并行算法可保持通信时间稳定性,在并行计算机上有很好的扩展性, 适合大规模稀疏特征问题的求解.
参考文献 | 相关文章 | 多维度评价
被引次数: Baidu(4)
16. 二维各向同性湍流直接数值模拟的六边形谱方法及其GPU实现和优化
乔海军, 李会元
数值计算与计算机应用    2013, 34 (2): 147-160.   DOI: 10.12288/szjs.2013.2.147
摘要1043)      PDF(pc) (1230KB)(427)    收藏
本文在六边形傅里叶分析及六边形快速傅里叶变换的基础上,提出了二维各向同性湍流直接数值模拟的对偶六边形傅里叶谱方法和六边形傅里叶谱方法,基于二维Navier-Stokes方程的涡度-速度形式, 构造了两种六边形傅里叶谱方法的离散格式,设计了其快速求解算法, 并且在GPU高性能平台上研制并优化了相应的数值模拟程序.根据方程的具体形式和六边形傅里叶谱方法的特点, 从算法层面对方程的求解过程,尤其是非线性Jacobian项快速计算进行优化, 经过优化之后, 方程求解算法的计算复杂度减少了约30%;根据GPU的体系结构和数值模拟程序的功能要求, 将计算模块全部设计为在GPU上运行的kernel函数, 尽量避免内存与显存之间的数据拷贝,并在软件工程层面上对代码进行性能优化. 优化后的GPU程序获得了高达50倍的加速比. 在此基础上, 我们对二维各向同性湍流进行了初步的数值模拟,并考察了在不同初始雷诺数条件下, 能量和拟涡能随着时间的演变曲线. 计算结果表明六边形傅里叶谱方法与传统的傅里叶谱方法一样高效精确.
参考文献 | 相关文章 | 多维度评价
17. 非奇异H-矩阵的新判定
张晋芳, 杨晋, 任艳萍
数值计算与计算机应用    2013, 34 (3): 161-166.   DOI: 10.12288/szjs.2013.3.161
摘要1094)      PDF(pc) (249KB)(275)    收藏
应用广义严格对角占优矩阵的性质, 对矩阵元素进行比较, 确定了在一定区间范围内的数值因子, 从而得到了一种判定非奇异H矩阵的新的方法. 利用矩阵理论中不等式的方法和技巧, 构造出相对应的正对角矩阵, 并给出了严格的推导证明. 由此推广得到了满足一定条件下的不可约矩阵以及具有非零元素链的矩阵, 从而也得到了非奇异H矩阵的另外两种新的判定方法, 最后用数值例子说明了结论的有效性.
参考文献 | 相关文章 | 多维度评价
被引次数: Baidu(1)
18. 求解非线性时滞双曲型偏微分方程的紧致差分方法及Richardson外推算法
张启峰, 张诚坚, 邓定文
数值计算与计算机应用    2013, 34 (3): 167-176.   DOI: 10.12288/szjs.2013.3.167
摘要1837)      PDF(pc) (1088KB)(806)    收藏
本文构造了一类求解非线性时滞双曲型偏微分方程的紧致差分格式, 获得了该差分格式的唯一可解性, 收敛性和无条件稳定性, 收敛阶为O(τ2+h4), 并进一步对时间方向进行Richardson外推, 使得收敛阶达到O(τ4+h4). 数值实验表明了算法的精度和有效性.
参考文献 | 相关文章 | 多维度评价
19. 带有间断系数椭圆方程的加权间断Galerkin方法
张荣培, 蔚喜军, 崔霞
数值计算与计算机应用    2013, 34 (3): 177-186.   DOI: 10.12288/szjs.2013.3.177
摘要1256)      PDF(pc) (454KB)(517)    收藏
本文应用间断Galerkin(DG)方法求解带有间断系数的二维椭圆方程. 针对扩散系数间断的特点, 我们构造一种新的加权对称内惩罚方法. 证明了相应双线性形式的连续性和强制性, 并给出收敛性证明. 数值算例表明我们的DG方法对于求解强间断系数问题十分有效.
参考文献 | 相关文章 | 多维度评价
被引次数: Baidu(1)
20. 广义Wang-Ball曲线
熊建, 郭清伟
数值计算与计算机应用    2013, 34 (3): 187-195.   DOI: 10.12288/szjs.2013.3.187
摘要1081)      PDF(pc) (407KB)(313)    收藏
给出了n次带形状参数λ的Wang-Ball曲线, 它具有n次Wang-Ball曲线的类似性质.形状参数λ具有明显的几何意义: λ越大, 曲线越逼近控制多边形.当λ=0时, 曲线退化为一条线段;当λ=2时, 曲线退化为Wang-Ball曲线.给出了曲线的递归求值, 升阶和降阶逼近算法, 用Bézier形式表达的系数公式及两段曲线G1, C1连续拼接的条件.
参考文献 | 相关文章 | 多维度评价
21. RLC串联电路与运动介质板综合分析
杨志安
数值计算与计算机应用    2013, 34 (3): 196-204.   DOI: 10.12288/szjs.2013.3.196
摘要783)      PDF(pc) (421KB)(282)    收藏
考虑磁路非线性特点, 基于磁共能、电场能和电系统耗散函数, 通过引入拉格拉日函数得到基尔霍夫定律的拉格朗日方程表示法. 应用此方法建立了RLC串联电路与运动介质板系统数学模型. 给出介质板所受的电磁力和插入介质板电容器电容. 根据电路特点进行理论分析和数值计算, 结果表明电荷经过几秒震荡后达到稳定运动.
参考文献 | 相关文章 | 多维度评价
被引次数: CSCD(1)
22. 基于改进离散粒子群优化算法的作战弹药分配研究
赵志宁, 石全, 张军刚
数值计算与计算机应用    2013, 34 (3): 205-211.   DOI: 10.12288/szjs.2013.3.205
摘要915)      PDF(pc) (419KB)(523)    收藏
结合作战弹药分配的具体实际, 利用广义指派问题描述与数学模型, 描述了弹药分配算法存在的问题与对策, 提出了一种基于连续空间的离散粒子群优化算法以求解该模型. 算法针对问题设计了特定的粒子编码方案, 并引入了一种局部搜索以提高算法的局部搜索能力. 算例分析表明了所提离散粒子群优化算法求解弹药分配问题的可行性.
参考文献 | 相关文章 | 多维度评价
被引次数: CSCD(1)
23. 三维非结构网格上求解双曲型守恒律方程的一类三阶精度有限体积格式
唐玲艳, 傅浩, 宋松和
数值计算与计算机应用    2013, 34 (3): 212-220.   DOI: 10.12288/szjs.2013.3.212
摘要1058)      PDF(pc) (746KB)(427)    收藏
考虑标量双曲型守恒律方程, 对三维非结构四面体网格给出了一类满足局部极值原理的三阶精度有限体积格式. 方法的主要思想是时间和空间分开处理, 时间离散采用三阶TVD Runge-Kutta方法; 对空间, 在每一个四面体单元上基于最小二乘原理构造一个二次多项式,结合数值解光滑探测器和梯度限制器, 使其在光滑区域具有高阶精度, 在间断区域满足局部极值原理. 该格式具有间断分辨能力高, 编程实现简便, 计算速度快等优点. 典型算例的数值试验表明, 该格式是有效的.
参考文献 | 相关文章 | 多维度评价
被引次数: CSCD(2)
24. 用CCD法离散求解二维Helmholtz方程的数值方法
柯日焕, 黎稳
数值计算与计算机应用    2013, 34 (3): 221-230.   DOI: 10.12288/szjs.2013.3.221
摘要1348)      PDF(pc) (434KB)(471)    收藏
用联合紧致差分格式(CCD)离散Helmholtz方程, 具有6阶精度.然而对于得到的线性方程组, 我们仍需一种高效求解方法.本文针对二维的Helmholtz方程CCD离散所得的线性方程组给出高效的数值方法.数值例子表明所提出的方法是有效的.
参考文献 | 相关文章 | 多维度评价
被引次数: CSCD(2)
25. 三维变分资料同化系统并行算法设计与实现
王玉柱, 姜金荣, 蔡长青, 迟学斌, 岳天祥
数值计算与计算机应用    2013, 34 (3): 231-240.   DOI: 10.12288/szjs.2013.3.231
摘要1695)      PDF(pc) (753KB)(726)    收藏
三维变分资料同化作为现在主流数值天气预报的同化方法, 能够明显改善预报数据的质量. 随着科学研究的逐渐深入以及科学探测仪器和计算机技术的不断发展, 受计算量和内存需求量的限制, 传统串行三维变分资料同化系统已无法满足高分辨率、高精确度数值预报的要求. 所以, 三维变分资料同化系统的并行设计与实现显得尤其重要. 本文设计了混合二维区域剖分并行化方法及其通信算法库, 并将其应用于国家气象局三维变分同化系统3DVAR. 数值试验表明, 系统128核的并行效率相对于2核高达72%, 具有良好的加速效果; 同时, 内存需求也随处理器个数的增加而成倍减少, 满足了高分辨率预报的要求.
参考文献 | 相关文章 | 多维度评价
被引次数: Baidu(1)
26. 随机延迟积分微分方程改进分步向后Euler方法的均方指数稳定性
李启勇, 甘四清, 张浩敏
数值计算与计算机应用    2013, 34 (4): 241-248.   DOI: 10.12288/szjs.2013.4.241
摘要994)      PDF(pc) (362KB)(421)    收藏
本文研究一类改进分步向后Euler 方法求解随机延迟积分微分方程的均方指数稳定性。证明了在约束网格下,该方法依步长h=τ/m 保持原系统的均方指数稳定性。数值试验验证了本文理论结果的正确性。
参考文献 | 相关文章 | 多维度评价
被引次数: Baidu(1)
27. 求解多集分裂可行问题的一种共轭梯度法
张艳君, 赵金玲, 徐尔
数值计算与计算机应用    2013, 34 (4): 249-256.   DOI: 10.12288/szjs.2013.4.249
摘要1056)      PDF(pc) (407KB)(339)    收藏
基于求解多集分裂可行问题与非线性最优化问题的等价性,考虑Jinling Zhao and Qingzhi Yang在[1]中提出的求解SFP的共轭梯度法和Censor等在[2]中提出的梯度投影法,尝试运用共轭梯度法求解多集分裂可行问题; 并且证明了所构造算法的收敛性。提出的新算法克服了求矩阵逆的缺点。初步的数值结果表明新算法对于不同的问题都能够有较快的收敛速度,具有良好的稳定性和可行性,在问题维数增大时表现得越发明显。
参考文献 | 相关文章 | 多维度评价
28. 奇异摄动问题在Bakhvalov-Shishkin网格上的有限元超收敛
尹云辉, 祝鹏, 杨宇博
数值计算与计算机应用    2013, 34 (4): 257-265.   DOI: 10.12288/szjs.2013.4.257
摘要1018)      PDF(pc) (359KB)(292)    收藏
在Bakhvalov-Shishkin网格上,利用线性插值的Galerkin有限元方法求解一维对流扩散型的奇异摄动问题。在εN-1的前提下,通过使用离散的能量范数,可以得到,关于扰动参数ε是一致收敛的,其误差阶达到ON-2)。最后,通过数值算例,验证了理论分析。
参考文献 | 相关文章 | 多维度评价
被引次数: Baidu(5)
29. 计算机三维图形中一个平方根倒数算法的数学证明及优化
程竟然
数值计算与计算机应用    2013, 34 (4): 266-278.   DOI: 10.12288/szjs.2013.4.266
摘要1022)      PDF(pc) (580KB)(226)    收藏
在计算机三维图形领域中存在一个仅使用位操作和减法运算的快速平方根倒数算法。作者利用代数分析结合Mathematica软件提供的数值方法,通过构造相对误差目标函数并最小化之,证明了此算法的有效性,指出了其适用条件并发现了算法中一个常数的较优替代值。此外如果按输入分两类分别进行计算能再提高精度。进一步利用C语言编程手段验证了优化后的算法。
参考文献 | 相关文章 | 多维度评价
30. 基于角点检测和奇异值分解的多重数字水印算法
刘新, 张慧娟, 孙彬彬, 王艳超
数值计算与计算机应用    2013, 34 (4): 279-285.   DOI: 10.12288/szjs.2013.4.279
摘要884)      PDF(pc) (547KB)(211)    收藏
针对多媒体数据的版权保护和完整性认证的多方性,提出一种基于角点检测和奇异值分解的多重数字水印算法。本算法将鲁棒水印嵌入到小波域低频子带角点中,将半脆弱水印嵌入到小波域分块细节子带奇异值中,实验仿真结果表明,水印算法能够抵抗剪切、涂改、滤波、压缩等常见攻击,同时也能够抵抗旋转、缩放等几何攻击,具有较好的鲁棒性和脆弱敏感性,且准确定位篡改区域,可谓一种有效的多重数字水印方案。
参考文献 | 相关文章 | 多维度评价
31. 一类广义Riccati矩阵方程对称解的双迭代算法
张凯院, 宋卫红, 王娇
数值计算与计算机应用    2013, 34 (4): 286-294.   DOI: 10.12288/szjs.2013.4.286
摘要1270)      PDF(pc) (391KB)(266)    收藏
研究了一类广义系统控制理论导出的Riccati矩阵方程对称解的数值计算方法.运用牛顿算法将Riccati矩阵方程的对称解问题转化为线性矩阵方程的对称解或者对称最小二乘解问题,采用修正共轭梯度法解决导出的线性矩阵方程的对称解问题,可建立求Riccati矩阵方程对称解的双迭代算法。数值算例表明,双迭代算法是有效的。
参考文献 | 相关文章 | 多维度评价
被引次数: Baidu(2)
32. 具有巴黎期权特性的可转债定价问题研究
吴正, 阳佳慧, 张依文
数值计算与计算机应用    2013, 34 (4): 295-304.   DOI: 10.12288/szjs.2013.4.295
摘要1060)      PDF(pc) (863KB)(266)    收藏
采用有限差分方法对基于Black-Scholes方程的可转债定价模型进行数值求解,用Euler-Lagrange分裂格式离散包含具有巴黎期权特性的赎回条款的修正Black-Scholes方程,并以工行转债和中行转债的历史数据为例,比较不同的数学模型中定价结果与实际价格的差异,分析标的股票处在不同价位水平时不同定价模型对可转债问题的适用性。
参考文献 | 相关文章 | 多维度评价
被引次数: Baidu(2)
33. 凹角区域外问题基于椭圆弧人工边界的自然边界元与有限元耦合法
陈亚军, 杜其奎
数值计算与计算机应用    2013, 34 (4): 305-311.   DOI: 10.12288/szjs.2013.4.305
摘要939)      PDF(pc) (397KB)(431)    收藏
本文以具有凹角长条型内边界的调和方程外问题为例,研究一种以椭圆弧为人工边界的自然边界元与有限元耦合法,给出了耦合变分问题的适定性及近似解的误差估计。理论分析及数值结果表明,用该方法求解带凹角长条型内边界的外问题是十分有效的。
参考文献 | 相关文章 | 多维度评价
34. 大气数值模式中FFT算法的设计与实现
王玉柱, 姜金荣, 迟学斌, 岳天祥
数值计算与计算机应用    2013, 34 (4): 312-322.   DOI: 10.12288/szjs.2013.4.312
摘要1130)      PDF(pc) (501KB)(393)    收藏
大气数值模式中常用到的FFT 软件包FFT99 由于没有考虑硬件配置和FFT 算法的变换参数影响,造成计算效率不高。因此,本文基于FFTW3,根据大气数值模式的特殊需求,设计实现了新的FFT99 软件包SC_FFT99。数值试验表明,在FFT计算速度方面,根据新算法设计的SC_FFT99软件包比FFT99 有较大的提高。理想情况下,SC_FFT99 计算速度比FFT99快2.5 到3.5 倍。并且SC_FFT99 已应用于中国科学院大气物理研究所开发的大气环流模式IAP AGCM4.0 的数值模拟,测试结果显示:计算速度比原来的FFT99 快0.39 倍。
参考文献 | 相关文章 | 多维度评价