推荐文章

    Please wait a minute...
    复合凸优化的快速邻近点算法
    郦旭东
    2020, 42 (4): 385-404.   DOI: 10.12286/jssx.2020.4.385
    摘要237)      PDF (598KB)(273)   
    在大数据时代,随着数据采集手段的不断提升,大规模复合凸优化问题大量的出现在包括统计数据分析,机器与统计学习以及信号与图像处理等应用中.本文针对大规模复合凸优化问题介绍了一类快速邻近点算法.在易计算的近似准则和较弱的平稳性条件下,本文给出了该算法的全局收敛与局部渐近超线性收敛结果.同时,我们设计了基于对偶原理的半光滑牛顿法来高效稳定求解邻近点算法所涉及的重要子问题.最后,本文还讨论了如何通过深入挖掘并利用复合凸优化问题中由非光滑正则函数所诱导的非光滑二阶信息来极大减少半光滑牛顿算法中求解牛顿线性系统所需的工作量,从而进一步加速邻近点算法.
    参考文献 | 相关文章 | 多维度评价
    电子结构计算的数值方法与理论
    戴小英
    2020, 42 (2): 131-158.   DOI: 10.12286/jssx.2020.2.131
    摘要515)      PDF (690KB)(431)   
    第一原理电子结构计算已成为探索与研究物质机理、理解与预测材料性质的重要手段和工具.虽然第一原理电子结构计算取得了巨大的成功,但是如何利用高性能计算机又快又好地计算大规模体系,如何从数学角度理解电子结构模型的合理性与计算的可靠性和有效性,依然充满各种挑战.基于密度泛函理论的第一原理电子结构计算的核心数学模型为Kohn-Sham方程或相应的Kohn-Sham能量泛函极小问题.近年来,人们分别从非线性算子特征值问题的高效离散及Kohn-Sham能量泛函极小问题的最优化方法设计两个方面对电子结构计算的高效算法设计及分析展开了诸多研究.本文重点介绍我们小组在电子结构计算的方法与理论方面的一些进展,同时简单介绍该领域存在的困难与挑战.
    参考文献 | 相关文章 | 多维度评价
    弱有限元方法在线弹性问题中的应用
    张然
    2020, 42 (1): 1-17.   DOI: 10.12286/jssx.2020.1.1
    摘要480)      PDF (482KB)(530)   
    本文考虑弱有限元(简称WG)方法在线弹性问题中的应用.WG方法是传统有限元方法的推广,用于偏微分方程的数值求解.和传统有限元一样,它的基本思想源于变分原理.WG方法的特点是使用在剖分单元内部和剖分单元边界上分别有定义的分片多项式函数(即弱函数)作为近似函数来逼近真解,并针对弱函数定义相应的弱微分算子代入数值格式进行计算.除此之外,WG方法允许在数值格式中引进稳定子以实现近似函数的弱连续性.WG方法具有允许使用任意多边形或多面体剖分,数值格式与逼近函数构造简单,易于满足相应的稳定性条件等优点.本文考虑WG方法在求解线弹性问题中的应用.围绕线弹性问题数值求解中常见的三个问题,即:数值格式的强制性,闭锁性,应力张量的对称性介绍WG方法在线弹性问题求解中的应用.
    参考文献 | 相关文章 | 多维度评价
    图像反问题中的数学与深度学习方法
    董彬
    2019, 41 (4): 343-366.   DOI: 10.12286/jssx.2019.4.343
    摘要1655)      PDF (695KB)(1192)   
    我们生活在数字的时代,数据已经成为了我们生活中不可或缺的一部分,而图像无疑是最重要的数据类型之一.图像反问题,包括图像降噪,去模糊,修复,生物医学成像等,是图像科学中的重要领域.计算机技术的飞速发展使得我们可以用精细的数学和机器学习工具来为图像反问题设计有效的解决方案.本文主要回顾图像反问题中的三大类方法,即以小波(框架)为代表的计算调和分析法、偏微分方程(PDE)方法和深度学习方法.我们将回顾这些方法的建模思想和一些具体数学形式,探讨它们之间的联系与区别,优点与缺点,探讨将这些方法有机融合的可行性与优势.
    参考文献 | 相关文章 | 多维度评价
    无中心优化的算子分裂方法
    印卧涛
    2019, 41 (3): 225-241.   DOI: 10.12286/jssx.2019.3.225
    摘要645)      PDF (552KB)(465)   
    在某些多智能体系统中,由于受到通讯等因素的限制,单个智能体只能进行本地计算,再与相邻智能体交换数据.与传统的并行和分布式计算不同,这种数据交换方式不再使用中心节点或者共享内存,而仅限于相邻节点之间.这种通过局部数据交换而实现全网目标的方式叫做无中心计算.比如,从任意的多个数开始,所有智能体通过不断地计算其局部平均,就都能收敛到这些数的平均值.无中心计算有不易形成通讯和计算瓶颈的优点,更适合分布的节点,因此受到一些应用的欢迎.
    本文介绍求解一致最优化问题的若干无中心算法.一致最优化问题的目标是全网所有节点的变量收敛到同一个、并使所有目标函数之和最小的值.我们可以通过推广求平均的无中心方法去实现这个目标,但是得到算法比普通(有中心的)优化算法收敛得更慢,有阶数差距.近年来,一些新的无中心算法弥补了这个阶数差距.本文采用算子分裂的统一框架,以比这些算法原文更为简单的形式介绍这些方法.
    参考文献 | 相关文章 | 多维度评价
    类Hartree-Fock方程的数值方法
    林霖
    2019, 41 (2): 113-125.   DOI: 10.12286/jssx.2019.2.113
    摘要335)      PDF (531KB)(570)   
    本文的主要目的是介绍近年来大基组下的类Hartree-Fock方程数值求解的一些进展.类Hartree-Fock方程出现在Hartree-Fock理论和含杂化泛函的Kohn-Sham密度泛函理论中,是电子结构理论中一类重要的方程.该方程在复杂的化学和材料体系的电子结构计算中有广泛地应用.由于计算代价的原因,类Hartree-Fock方程一般只被用在较小规模的量子体系(含几十到几百个电子)的计算.从数学角度上讲,类Hartree-Fock方程是一个非线性积分-微分方程组,其计算代价主要来自于积分算子的部分,也就是Fock交换算子.通过发展和结合自适应压缩交换算子方法(ACE),投影的C-DⅡS方法(PC-DⅡS)方法,以及插值可分密度近似方法(ISDF),我们大大降低了杂化泛函密度泛函理论的计算代价.以含1000个硅原子的体系为例,我们将平面波基组下的杂化泛函的计算代价降至接近不含Fock交换算子的半局域泛函计算的水平.同时,我们发现类Hartree-Fock方程的数学结构也为一类特征值问题的迭代求解提供了新的思路.
    参考文献 | 相关文章 | 多维度评价
    信赖域方法的收敛性
    袁亚湘
    1994, 16 (3): 333-346.   DOI: 10.12286/jssx.1994.3.333
    摘要1386)      PDF (554KB)(777)   
    信赖域方法的收敛性袁亚湘(中国科学院计算中心)ONTHECONVERGENCEOFTRUSTREGIONALGORITHMS¥YuanYa-xiang(ComputingCenterAcademiaSinica)Abstract:Trustregio...
    参考文献 | 相关文章 | 多维度评价
    被引次数: Baidu(115)
    非奇H矩阵的简捷判据
    黄廷祝
    1993, 15 (3): 318-328.   DOI: 10.12286/jssx.1993.3.318
    摘要1609)      PDF (320KB)(923)   
    非奇H矩阵在计算数学和矩阵理论的研究中很重要,但简便实用的判定条件较少见。本文给出几个简捷判据。[1,2,3]的主要结果是本文定理1的特例。 记M_n(C)为n阶复阵集合,M_n(R)为n阶实阵集合。设A=(a_(ij))∈M_n(C),记Λ_i(A)=sum from j≠i to |a_(ij)|,i,j∈N≡{1,2,…,n}。若|a_(ii)|>Λ_i(A),i∈N,则称A
    参考文献 | 相关文章 | 多维度评价
    被引次数: Baidu(552)
    关于Morley元的误差估计
    石钟慈
    1990, 12 (2): 113-118.   DOI: 10.12286/jssx.1990.2.113
    摘要1386)      PDF (196KB)(886)   
    §1.引言 解薄板弯曲问题的三角形Morley元是六十年代出现的一种非协调元,它的形函数是完整的二次多项式,节点参数是单元顶点上的三个函数值及三边中点上的法向导数值.由于板弯曲问题的常应变是二次多项式,所以这是一个参数最少的非协调板元.由
    参考文献 | 相关文章 | 多维度评价
    被引次数: Baidu(76)
    关于九参数二阶拟协调元
    石钟慈
    1989, 11 (4): 442-444.   DOI: 10.12286/jssx.1989.4.442
    摘要1338)      PDF (93KB)(617)   
    韩厚德最近在讨论拟协调元时,引进了一个九参数二阶拟协调元.在[1]中,对一个完全三次多项式(十个自由度)附加一个特殊的约束条件并使它满足所谓的二阶拟协调条件,即形函数及其两个一阶偏导数在单元之间的内部边界上保持积分意义下的连续性.我们证明,这样得到的九参二阶拟协调元实际上就是熟知的 de Veubeke元.其证明如下.
    参考文献 | 相关文章 | 多维度评价
    被引次数: Baidu(2)
    实对称矩阵的两类逆特征值问题
    孙继广
    1988, 10 (3): 282-290.   DOI: 10.12286/jssx.1988.3.282
    摘要1346)      PDF (277KB)(704)   
    §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的特征值的
    参考文献 | 相关文章 | 多维度评价
    被引次数: Baidu(195)
    关于三次样条插值矩阵的非奇异性
    石钟慈
    1983, 5 (2): 195-203.   DOI: 10.12286/jssx.1983.2.195
    摘要1374)      PDF (273KB)(807)   
    近年来,在计算数学刊物上相继发表了许多篇关于三次插值样条存在唯一性的文章,例如[1-3].这些文章讨论的是三次样条插值矩阵为非奇异的条件.[1]中用的是凑方法,讨论了与插值矩阵相关的另一个对称阵为正定的条件,经过复杂凑方,得到了某些充分条件,[2]是用大块凑方,所得结果形式上异于[1],但实质上是完全相同的.[3]则是对插值矩阵进行一种特殊分解,得出非异的四个充分条件.它不限于[1-2]所讨论的正定情形,因而适用范围更广些.
    参考文献 | 相关文章 | 多维度评价
    被引次数: Baidu(9)
    高阶稀疏局部非线性方程组的一种拟牛顿方法
    崔俊芝
    1982, 4 (4): 356-364.   DOI: 10.12286/jssx.1982.4.356
    摘要1362)      PDF (356KB)(759)   
    §1.引言 当用有限元法或有限差分法分析非线性偏微分方程问题时,必然会导致求解非线性方程组的问题,即求 F(x)=0 (1.1)的解.其中,x=(x_1,x_2,…,Xx_n)~T∈D,D?R~n;F:D→R~n是一个非线性映射.因此,有效地求解非线性方程组(1.1),是分析相应的非线性问题的关键. 不管这些非线性问题是来自流体力学、固体力学,还是其他的物理范畴,它们所对应
    参考文献 | 相关文章 | 多维度评价
    论Bézier曲线的仿射不变量
    苏步青
    1980, 2 (4): 289-298.   DOI: 10.12286/jssx.1980.4.289
    摘要1367)      PDF (284KB)(786)   
    本文的目的是按照[1]的理论找出n次平面Bezier曲线的内在仿射不变量,特别是,对于3次Bezier曲线的保凸性作出其充要条件的几何解释。对于一般的情况下的保凸性问题,至今还没有解决。著者仅在4次的场合详尽地讨论了曲线段上是否存在拐点的分析的(而不是几何的)充要条件,而最后举出几个实例,以说明特征多角形的凸性是充分条件,而不是必要条件。
    参考文献 | 相关文章 | 多维度评价
    被引次数: Baidu(164)
    流体力学中的差分方法(Ⅱ)——三维涡度方程的数值解
    郭本瑜
    1980, 2 (4): 307-318.   DOI: 10.12286/jssx.1980.4.307
    摘要1369)      PDF (430KB)(795)   
    有关三维涡度方程的数值计算方面的工作已有[1—3],但缺乏比较系统的理论分析.在[4]中,以二维涡度方程为例,讨论了流体力学差分方法的一些理论问题.本文是把这些结果推广到三维.
    参考文献 | 相关文章 | 多维度评价
    被引次数: Baidu(15)
    不变子空间与广义不变子空间——(Ⅱ)扰动定理
    孙继广
    1980, 2 (2): 113-123.   DOI: 10.12286/jssx.1980.2.113
    摘要1496)      PDF (367KB)(617)   
    关于矩阵的不变子空间,自然会提出这样一个扰动问题:设Z_1∈C~(n×l)是A∈C~(n×n)的一个特征矩阵,若E∈C~(n×n)是一个扰动矩阵,问A+B是否存在特征矩阵Z_1,使得(Z_1)靠近R(Z_1)?关于矩阵对的广义不变子空间.也可以类似地提出问题。 对于这些问题,G.W.Stewart曾经讨论过,他的方法的关键是构造一种求解二次矩阵方程的迭代过程,用来逼近矩阵的一个不变子空间;而本文建议另一种迭代格式,用这种迭代逼近一个不变(或广义不变)子空间,具有二次收敛速度。
    参考文献 | 相关文章 | 多维度评价
    被引次数: Baidu(13)
    广义局部坐标样条
    孙家昶
    1980, 2 (2): 142-145.   DOI: 10.12286/jssx.1980.2.142
    摘要1366)      PDF (177KB)(871)   
    我们在[1]中提出了适用于“大挠度”曲线插值的方法,这种方法在实际运用时仍有某些不便之处.例如,[1]只考虑了通过节点的插值,不能直接用作曲线拟合;一律采用各自的局部坐标,对有些形状相近的曲线族的特点没能很好地利用;有些曲线的“大挠度”只占极小部份,并不需要每次转坐标;有的离散数据间断性质较强,用[1]得到的节点关系式中非线性项是主要项,这时简单迭代往往发散,要用牛顿迭代,但会增加运算量。为了
    参考文献 | 相关文章 | 多维度评价
    被引次数: Baidu(5)
    不变子空间与广义不变子空间(Ⅰ)存在与唯一性定理
    孙继广,
    1980, 2 (1): 1-13.   DOI: 10.12286/jssx.1980.1.1
    摘要1378)      PDF (488KB)(822)   
    本文讨论与特征值和广义特征值问题相联系的某些子空间。在本文中,我们定义了矩阵对的“广义特征值方阵对”和“广义特征矩阵”,并由此建立了广义不变子空间的概念;建立了对应的子空间存在与唯一的充分必要条件;给出了广义不变子空间与G.W.Stewart定义的“收缩对”的关系。
    参考文献 | 相关文章 | 多维度评价
    被引次数: Baidu(15)
    试验函数不满足边界条件的一种有限元法
    周天孝,
    1980, 2 (1): 50-62.   DOI: 10.12286/jssx.1980.1.50
    摘要1423)      PDF (563KB)(701)   
    对于第一边界值问题,发展了一种使试验函数不必满足边界条件的有限元法,指出了此种新的处理具有通常有限元法应该有的稳定性与最佳收敛性特征。
    参考文献 | 相关文章 | 多维度评价
    被引次数: Baidu(2)
    论微分与积分方程以及有限与无限元
    冯康,
    1980, 2 (1): 100-105.   DOI: 10.12286/jssx.1980.1.100
    摘要1610)      PDF (321KB)(785)   
    1 椭圆微分方程的边界值问题可以有种种不同的数学成型,在理论上等价,但在实践上不等效。有限元方法成功的一个关键就是合理选取了变分的数学型式。举例来说,取调和方程的第二类边界问题,定义于区域Ω,具有光滑边界Г:
    参考文献 | 相关文章 | 多维度评价
    被引次数: Baidu(89)