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

阅读排行

  • 一年内发表的文章
  • 两年内
  • 三年内
  • 全部
Please wait a minute...
  • 全选
    |
  • 论文
    李瑞, 王宝华, 李巧云, 吴淑红
    数值计算与计算机应用. 2024, 45(2): 83-114. https://doi.org/10.12288/szjs.s2023-0940
    页岩/致密油气藏等非常规油气藏的储层介质为典型的多尺度空间: 既具有纳微米级的基质孔隙, 又具有微米-厘米级的天然裂缝, 以及大规模压裂所产生的米级到百米级的人工裂缝. 本文针对裂缝多孔介质, 介绍基质、天然裂缝、宏观裂缝及大尺度缝洞中的流体耦合流动数学模型.
  • 论文
    郑华, 温海斌, 卢晓平
    数值计算与计算机应用. 2023, 44(4): 350-367. https://doi.org/10.12288/szjs.s2022-0868
    针对大规模垂直线性互补问题的求解,运用二步多分裂技术构建了二步模基矩阵同步多分裂并行迭代方法,新方法可以看成是已有文献中模基矩阵同步多分裂迭代法和二步模基矩阵分裂迭代法的推广.进一步地,在系统矩阵为$H_+$矩阵的假设下,给出算法收敛性分析,得到了参数矩阵的收敛域,推广了已有算法的收敛性结果.最后,在OpenMP框架下针对已有文献中的两个数值例子给出了数值试验,试验结果展示了二步多分裂技术能提升已有方法的计算效率.
  • 论文
    熊均达, 李敬来
    数值计算与计算机应用. 2023, 44(3): 225-236. https://doi.org/10.12288/szjs.s2022-0821
    贝叶斯推断正逐渐成为解决反问题的一个越来越流行的工具, 这主要是因为它能够描述所获得的解的不确定性. 在许多实际的贝叶斯反问题中, 可能有多个潜在性能良好的模型可用来描述数据和感兴趣的参数. 在这种情况下, 基于贝叶斯的模型选择方法经常被用来选择"最佳"模型, 而这种方法依赖于对后验分布中归一化常数的估计. 因此, 估计归一化常数成为贝叶斯推断中的一项重要任务, 而这个问题对于标准抽样方法来说具有计算上的挑战性. 在这项工作中, 新提出的一种基于多正则蒙特卡洛技术(MMC)的归一化常数估计方法, 可以很好的解决上述难题, 这是一种自适应的重要性采样方法. 该方法可以以黑箱方式估计归一化常数, 使其特别适用于具有复杂基础模型的问题. 最后, 通过数值例子证明了所提出的方法可以有效且准确地计算出归一化常数的估计值.
  • 论文
    陈昊, 葛志昊
    数值计算与计算机应用. 2023, 44(4): 378-391. https://doi.org/10.12288/szjs.s2022-0873
    压缩感知(CS)是快速核磁共振成像(MRI)的一种有效的方法,该方法旨在从少量的欠采样数据空间中重建核磁共振图像,加速核磁共振成像中的数据采集.为了提高当前MRI图像的重建精度和速度,本文将传统的基于模型驱动的压缩感知恢复算法和数据驱动的深度学习方法相结合,提出了一种新的算法--ADMM-CNN.该算法将交替方向乘子法(ADMM)中的一个子问题用卷积神经网络进行求解,然后通过构造损失函数,利用梯度下降法优化ADMM中的惩罚参数和拉格朗日乘子的更新率,从而得到训练好的ADMM模型,能够有效地从压缩感知测量中恢复核磁共振图像信号.本文设计的算法与传统的重建算法相比,具有较少的正则化参数和更低的计算复杂度,同时,将一个训练好的卷积神经网络直接引入,极大简化了训练过程和计算量,并且利用GPU来进行加速训练.实验结果表明,与传统的重建算法和其他深度学习方法相比,ADMM-CNN拥有较快的计算速度和良好的重建精度.
  • 论文
    沈静, 杜育松
    数值计算与计算机应用. 2024, 45(1): 54-67. https://doi.org/10.12288/szjs.s2023-0917
    Karney于2016年提出了一种针对标准正态分布的精确采样算法.本文给出一种针对标准差为$\sqrt{1/(2\ln2)}$均值为$0$的正态分布的精确采样算法.这一种特殊的正态分布也被称为二元高斯分布,因为其相对概率密度函数可以由$2^{-x^2}$给出,这里$x$为任意实数.在实际中,针对二元高斯分布的这一精确采样算法无需浮点运算,可以看成是Karney精确采样技术的一种推广.分析了该采样算法产生一个二元高斯样本平均需要的区间$(0,1)$上的均匀偏差数.数值实验也表明了该采样算法的有效性.对于大于$1$但小于自然常数$e$的任意有理数$c$,将精确采样二元高斯分布的思想推广到了精确采样标准差为$\sqrt{1/(2\ln{c})}$均值为$0$的被称为“$c$元高斯分布”的一类正态分布上,并进行了类似的复杂性分析.
  • 论文
    王钦, 李维国, 韦林香, 张飞雨
    数值计算与计算机应用. 2023, 44(4): 337-349. https://doi.org/10.12288/szjs.s2022-0863
    块Kaczmarz方法是求解大规模线性方程组的一种迭代算法.在每次迭代时,会将当前迭代点正交投影到约束子集的解空间.文章基于块Kaczmarz方法的原理,提出了高斯混合模型的随机块Kaczmarz方法(GMM-RBK (k)).其中块的划分是通过高斯混合模型进行划分的.为了避免伪逆的计算或者最小二乘问题的求解,提出了高斯混合模型的随机平均块Kaczmarz方法(GMM-RABK (k)).证明了当线性方程组是相容时,这两种算法是收敛的,并给出相应的收敛率公式.在最后的数值实验中也证实了GMM-RBK (k)方法和GMM-RABK (k)方法的有效性,且无伪逆GMM-RABK (k)方法要优于GMM-RBK (k)方法.
  • 论文
    阮春蕾, 徐玉倩, 董层层
    数值计算与计算机应用. 2024, 45(1): 1-12. https://doi.org/10.12288/szjs.s2023-0875
    从Newton-Cotes积分的角度构造了一个带参数的三级二阶显式Runge-Kutta格式和两个带参数的四级三阶显式Runge-Kutta格式,给出了精度证明及稳定性条件.当参数特殊取值时,所构造的三种参数格式分别导出了常用的三阶Runge-Kutta格式(RK3)和经典的四阶Runge-Kutta格式(RK4).通过数值算例,验证了所构造的几类Runge-Kutta格式的有效性、稳定性及高精度.与常见的各类显式Runge-Kutta格式相比,本文构造的三种Runge-Kutta格式具有更好的稳定性.
  • 论文
    韦林香, 李维国, 王方
    数值计算与计算机应用. 2023, 44(3): 252-271. https://doi.org/10.12288/szjs.s2022-0834
    基于贪婪准则和最大距离准则选择系数矩阵工作列的策略, 提出两种求解大规模超定不相容线性系统的斜方向的 Gauss-Seidel 方法, 即斜方向的贪婪随机 Gauss-Seidel (GRGSO) 方法和斜方向的快速最大距离 Gauss-Seidel (FMDGSO) 方法. 当系数矩阵是列满秩时, 理论表明这些方法收敛到线性系统的唯一的最小二乘解. 特别是当矩阵A的列接近线性相关时, 数值结果表明这些方法在求解性能方面比传统的 Gauss-Seidel 型方法更具优势.
  • 论文
    周海林
    数值计算与计算机应用. 2024, 45(1): 27-42. https://doi.org/10.12288/szjs.s2023-0888
    应用共轭梯度方法和线性投影算子,给出迭代算法求解了线性矩阵方程$AXB=C$在任意线性子空间上的最小二乘解.在不考虑舍入误差的情况下,可以证明,所给迭代算法经过有限步迭代可得到矩阵方程$AXB=C$的最小二乘解、极小范数最小二乘解及其最佳逼近.文中的数值例子证实了该算法的有效性.本文算法的优点是在任意线性子空间上均容易实现.
  • 论文
    王俊杰
    数值计算与计算机应用. 2024, 45(2): 136-153. https://doi.org/10.12288/szjs.s2023-0923
    本文研究一类分数阶非线性薛定谔方程的数值格式, 该格式满足分数阶系统的一个或多个守恒性质. 首先, 我们基于\ BDF 格式、Crank-Nicolson 格式和松弛格式来离散时间导数, 并分析半离散格式的守恒性和色散误差. 其次,在周期边界条件下, 利用中心差分格式和紧致差分格式对分数阶非线性薛定谔方程的空间分数阶导数进行离散, 并证明这些数值格式保持质量和能量守恒律. 最后, 对一些分数阶非线性薛定谔方程进行了数值实验,验证了理论结果的正确性.
  • 论文
    聂宁明, 姚柯寒, 曾艳, 冯仰德, 王珏, 李顺德, 张纪林, 万健, 林克豪, 高岳, 王彦棡, 王宗国
    数值计算与计算机应用. 2024, 45(2): 115-135. https://doi.org/10.12288/szjs.s2023-0916
    本文结合大规模撕裂有限元方法和Newmark积分法, 对结构动力学问题进行高精细的大规模并行求解.面向异构平台, 设计了结点间和结点内的多级动静结合的负载均衡策略.在结点间, 根据撕裂有限元方法划分子域边界特点, 采用域边界平衡的图二分算法, 均衡各个子域的计算量;在结点内, 根据异构平台计算单元的性能差异, 进行了计算负载的动态优化.针对核心计算模块批量矩阵向量乘进行多流并行优化, 提升面向异构计算平台的利用率.本文优化已经集成到结构力学高性能数值模拟软件HARSA-feti中, 实验采用真实反应堆核燃料组件的流致振动仿真作为算例, 结果表明模拟性能提高了71.3%以上, 首次实现了百亿网格规模的全堆芯燃料棒组件的高精细模拟, 相较于1000块GPU, 16000 块GPU 的强、弱可扩展并行效率分别达到74.1%和81.1%.
  • 论文
    王旭浩, 王培培, 李昭祥, 陈先进
    数值计算与计算机应用. 2023, 44(3): 285-304. https://doi.org/10.12288/szjs.s2022-0860
    本文通过引入一种新的增广变换, 发展了改进的偏牛顿校正算法, 建立并证明了一类四阶非线性偏微分方程边值问题的新解与该问题零核空间的密切关系, 去掉了标准收敛假设, 使证明更简洁明了.分情况验证了该方程满足在 Nehari 子流形上全局分离定理的条件, 该分离定理为本文算法成功找到新解提供理论保障. 提出了二维非线性四阶偏微分方程 Dirichlet 边值问题的插值投影 Legendre-Galerkin 谱方法, 通过构造插值算子和投影算子, 对线性算子以及非线性项的处理进行了优化, 得到原问题的代数方程, 通过验证, 其与经典谱方法具有相同的条件数并都达到谱精度. 实验结果表明, 此方法与经典的谱方法或拟谱方法具有相同的收敛阶, 但计算所需CPU时间更少, 且能计算出更多的解.
  • 论文
    尚盼, 孔令臣
    数值计算与计算机应用. 2024, 45(2): 154-173. https://doi.org/10.12288/szjs.s2023-0924
    现代科学技术的发展使得各个领域中产生了大量的高维数据, 即样本特征量大于或远远大于样本数量的数据. 为了处理高维数据, 近年来有大量关于正则化回归模型的研究, 即通过引入调节参数将损失函数和正则项联合成一个目标函数, 比如著名的~LASSO 模型及相关模型. 众所周知, 对于正则化回归模型中最优调节参数的选择至关重要. 理论上: 该参数刻画了模型解的特征~(如稀疏性, 低秩性等), 从而决定了模型对数据的拟合效果; 计算上: 不同调节参数下模型的计算代价和计算效果不一样. 除了几类特殊的不需要进行最优调节参数选择的模型外, 目前最优调节参数选择的方法主要包含三类: 交叉验证, 信息准则及双层规划. 交叉验证及信息准则需要比较模型在不同调节参数下的解, 因此这两类最优调节参数选择方法需要多次求解模型. 除此之外, 如何更为合理地设置备选调节参数也需要进一步考虑. 为了降低交叉验证和信息准则进行最优调节参数选择的计算成本, 统计、 优化及机器学习三个方向的研究者们建立了不同的筛选规则, 即在不同调节参数下删除数据中不起作用的特征, 从而加速模型解的计算过程以达到加速最优调节参数选择过程的目的. 与交叉验证和信息准则不同, 双层规划是将最优调节参数选择问题刻画为一个双层规划模型, 通过求解模型来直接得到最优调节参数的选择结果. 本文从最优调节参数选择的方法和加速两个方面回顾现有结果, 并在此基础上提出未来的研究方向.
  • 论文
    赵金娥, 李明, 纪盈
    数值计算与计算机应用. 2023, 44(3): 305-312. https://doi.org/10.12288/szjs.s2022-0865
    针对二维椭圆型界面问题的离散化方程, 应用外推插值技巧和样条插值方法在细网格层上构造合适的迭代初始值, 加快V型多重网格法求解离散化系统的速度, 设计了外推完全多重网格(EXFMG)法. 数值实验表明新算法有效降低了迭代次数, 计算量更少.
  • 论文
    刘文杰, 张雪琳, 王汉权
    数值计算与计算机应用. 2023, 44(3): 237-251. https://doi.org/10.12288/szjs.s2022-0833
    近年来, 有关玻色-爱因斯坦凝聚态基态解的实验研究已经取得了一系列重要的成果. 本文在相关研究成果的基础上, 首先通过降维和无量纲化方法将玻色-爱因斯坦凝聚态基态解问题转换成能量泛函极值问题, 在离散该方程时,我们使用有限差分法和傅里叶谱方法分别离散该能量泛函方程的一维和二维情形. 其次, 本文采用了一种高效的数值优化算法-SQP优化方法来求解玻色-爱因斯坦凝聚态基态解问题并运用该算法分别对一维和二维的能量泛函极小值问题进行数值模拟. 最后通过分析实验数据结果和图像, 得出该算法能提高能量函数值的精确性.
  • 论文
    张艺博, 孟文辉
    数值计算与计算机应用. 2023, 44(4): 420-432. https://doi.org/10.12288/szjs.s2023-0896
    Gegenbauer加法定理是Bessel函数相关理论体系中一个重要的理论,在数学物理问题中有广泛的应用.本文在由Bessel函数和Neumann函数极限形式所导出的上界基础上,对Gegenbauer加法定理的截断误差上界进行了估计,得到了明确的误差界及其收敛阶.最后,通过数值实验验证了估计的有效性和精确性.
  • 论文
    李金, 张宇鑫
    数值计算与计算机应用. 2024, 45(1): 13-26. https://doi.org/10.12288/szjs.s2023-0887
    多维超奇异积分在弹性力学和电磁场的散射问题等诸多工程领域中有广泛应用.考虑构造二维、三维面型超奇异积分的求积公式,同时提高误差精度.利用复合矩形求积公式在划分的$N$个子区间内近似计算被积函数中无奇异性的部分,剩余部分通过超奇异积分的解析式求解.根据外推思想,构造一维超奇异积分的修正复合矩形求积公式.最后将带有外推的求积公式推广到二维、三维面型超奇异积分中.文章结尾的数值算例验证了方法的可行性.
  • 论文
    孙洪斌, 郭晓霞
    数值计算与计算机应用. 2024, 45(1): 43-53. https://doi.org/10.12288/szjs.s2023-0892
    本文考虑求解粒子运输中非对称代数Riccati方程最小正解的两个加速算法的收敛性问题.该方程含有两个可变参数$\alpha\in[0,1),c\in (0,1]$.我们证明了两个加速算法有相同的收敛速度,并且当$(\alpha,c)\neq (0,1)$时,两个加速算法线性收敛,当$(\alpha,c)=(0,1)$时,它们次线性收敛.
  • 论文
    涂梦凡
    数值计算与计算机应用. 2024, 45(1): 68-82. https://doi.org/10.12288/szjs.s2023-0920
    周期体系电子结构计算中的物理量需要通过在布里渊区上的能带积分得到.对于金属体系,被积函数在费米面穿过能带的地方出现间断,这导致一般数值积分方法的计算精度提升困难.本文基于smearing方法,分别通过一次和二次外推得到了关于展宽参数具有更高阶精度的数值格式.基于这类外推格式的数值方法在布里渊区$\mathbf{k}$点离散相同的前提下给出了更高精度的布里渊区积分逼近,显著地提高了计算效率.我们通过对两类典型体系的数值模拟验证了这类外推方法的有效性.
  • 论文
    刘洪
    数值计算与计算机应用. 2023, 44(4): 368-377. https://doi.org/10.12288/szjs.s2022-0871
    一维非稳态渗流模型求解是试井研究的基础,这类模型目前解法分解析解法和数值解法两大类,解析解法主要通过拉普拉斯变换及数值反演完成,数值解法则主要依赖网格划分和方程离散.以动边界思想为基础,本文提出了一维非稳态渗流模型第三类解法——非线性方程迭代解法,基于稳态依次替换法和质量守恒定律,视压力探测半径内区域渗流为稳态流动,将非稳态渗流方程转化为稳态方程,分析不同流动时间下探测半径变化规律,推导建立了不同外边界条件下探测半径的非线性控制方程,最后通过牛顿法求解非线性方程得到不同时刻井底压力,对比非线性方程迭代解与解析解显示计算误差较小,且主要是计算初期误差稍大.与传统的解析解法和数值解法相比,非线性方程迭代解法精度较高、算法简单,每个时间步仅需解一个非线性方程,显著减少了计算时间,可以在其他一维非稳态渗流问题中推广应用.
  • 论文
    贺亚丽, 汪佳玲, 黄家蓁
    数值计算与计算机应用. 2023, 44(4): 433-448. https://doi.org/10.12288/szjs.s2023-0911
    基于“数值算法应尽可能地保持原问题的本质特征”这一思想,本文利用复合构造方法为Zakharov系统构造了一个局部动量守恒算法.本文所提出格式的优点是它在任何时空区域都能够保持局部动量守恒和局部质量守恒,即不依赖于任何的边界条件.在合适的边界条件下,例如齐次边界条件和周期边界条件,所提出的算法也保持了全局的动量守恒和质量守恒.数值实验给出了算法的时空二阶收敛性,并且验证了该算法的长时间计算稳定性和不变量守恒性.
  • 论文
    林彭壮汉, 李凯, 程万友
    数值计算与计算机应用. 2023, 44(4): 409-419. https://doi.org/10.12288/szjs.s2023-0891
    本文提出一种基于积极集识别技术的临近牛顿算法用以求解$\ell_1$问题.该方法的一个优势在于利用了ISTA算法良好的支集辨认性质去确定自由变量和积极集变量,另一个优势在于利用了部分Hessian矩阵的信息去更新自由变量.在适当的条件下,我们证明了所提出的算法在使用非单调线搜索策略情况下是全局收敛的.数值实验证明提出的算法是有效的.
  • 论文
    范振成
    数值计算与计算机应用. 2023, 44(3): 327-336. https://doi.org/10.12288/szjs.s2023-0879
    描述芯片或电力系统运行规律的常用数学模型是高维微分代数方程组, 其中的微分方程组太大, 诸如线性多步法和Runge-Kutta(RK)法等经典数值方法均不能有效求解. 为求解这些微分方程组, 学者们提出了波形松弛(WR)方法. 多数情况下, 这些微分方程组是刚性的, 求解他们需要稳定性好的隐式方法, 尤其需要A-稳定的方法. 此外, RK方法是使用最广泛的常微分方程的数值方法. 然而, 迄今为止尚未发现RK型WR方法A-稳定的研究. 本文研究了RK型WR方法的A-稳定性, 获得了方法A-稳定的充分条件. 常见A-稳定的RK方法有Gauss-Legendre方法、Radau IA方法、Radau IIA方法、Lobatto IIIA方法、Lobatto IIIB 方法、Lobatto IIIC方法, 而且并非RK方法A-稳定, 相应的RK型WR方法也A-稳定. 在一个假设下, 本文所得结果说明当选择低阶Radau IA方法, 或Radau IIA方法, 或Lobatto IIIC方法为底层方法时, 存在分裂方式使得RK型WR方法是A-稳定的.
  • 论文
    张蕊, 阎爱玲
    数值计算与计算机应用. 2024, 45(2): 174-188. https://doi.org/10.12288/szjs.s2023-0929
    在高维数据分析中, 惩罚分位数回归是进行变量选择和参数估计的有效方法. 在实际应用中, 变量常以分组形式呈现, 为同时实现组间稀疏性和组内稀疏性, 本文研究了带稀疏Group Lasso惩罚的分位数回归模型. 为解决目标函数的非光滑性带来的计算挑战, 利用分位数Huber 函数近似分位数损失函数, 得到稀疏Group Lasso惩罚分位数Huber回归模型(SGLQHR). 基于Groupwise Majorization Descent (GMD) 算法提出了一种快速、有效算法求解该模型, 并建立算法收敛性. 数值实验和实例分析验证了该算法的有效性.
  • 论文
    梁雨欣, 刘东杰
    数值计算与计算机应用. 2023, 44(3): 313-326. https://doi.org/10.12288/szjs.s2023-0876
    本文在新的距离的框架下考虑 p-Laplace 方程的自适应有限元方法. 在文献[13]中作者考虑了2 ≤ p< ∞的情况, 给出了具有上界的后验误差估计.本文在此基础上发展了1 < p< ∞时具有上界和下界的后验误差估计. 数值实验验证了理论分析的结果.
  • 论文
    梅娜, 董巧丽, 何松年
    数值计算与计算机应用. 2023, 44(4): 392-408. https://doi.org/10.12288/szjs.s2023-0884
    双子空间投影算法和多步惯性随机Kaczmarz算法是求解相干线性方程组的有效算法,本文通过软阈值函数对这两种算法进行修正,提出了稀疏双子空间投影算法和稀疏多步惯性随机Kaczmarz算法,并给出其在有噪声干扰和无噪声干扰情况下在期望意义下的线性收敛率估计.通过数值实验验证本文所提算法的有效性和优越性.
  • 论文
    汪韬, 厉井钢
    数值计算与计算机应用. 2023, 44(3): 272-284. https://doi.org/10.12288/szjs.s2022-0854
    本文考虑一个次扩散方程的反问题, 即根据某个时刻解在空间区域上的积分反向优化方程中的分数阶参数. 本文说明了当初值为0时, 以上问题可能有多个解. 基于时空有限元法, 我们提出了一个参数优化算法, 并利用离散Laplace变换技巧证明了该算法的全局收敛性. 另外针对球形区域, 利用-Δ算子的谱分解, 我们还提出了一个快速算法. 最后两个数值算例验证了算法的有效性.