解无约束优化的一种新的共轭梯度法

本文提出了一种新的解无约束优化的共轭梯度算法,分析了算法的收敛性,并对算法进行了数值实验.数值实验的结果表明算法是有效的.     

渗流问题灰色数值模型的解法研究

灰色数值模型的求解是研究灰色数值模型的一个重要问题 .本文根据灰集合、灰数及其灰色运算规则 ,在渗流系统的基本灰色数值模型的基础上 ,分析了求解这类模型的一整套灰色数值算法 ,并对灰色数值算法、普通算法和经典数值方法的计算结果进行了全面比较 ,论证了灰色数值算法对灰信息传递的正确性和对渗流系统描述的合理性 .     

非线性互补问题的两种数值解法

本文研究了非线性互补问题的两类数值求解方法.在经典LQP算法及LevenbergMarquardt算法的基础上,构造了两种新算法,并证明了这两种新算法的收敛性.数值实验表明,新算法对测试问题优于已有算法.     

矩阵填充的硬阈值修正算法

提出了使用硬阈值进行矩阵填充的修正算法.算法通过对迭代矩阵进行对角修正来完成矩阵填充,并给出了算法的收敛性分析.最后通过数值实验比较了修正算法与硬阈值算法填充的数值结果,显示出了新算法的优越性.     

求解线性Caputo分数阶微分方程的高精度数值算法

首先,研究了计算Caputo分数阶导数的数值算法,并将其应用于求解线性Caputo分数阶微分方程的初值问题,提出了求解此问题的基础数值算法.如果初值条件不足,应用基础算法将无法计算出高精度的数值解.为了解决这一问题,构造辅助函数计算缺少的初值条件,进而设计出高精度的数值算法.计算实例说明,高精度算法是有效的.     

一类抛物型方程初值问题的随机数值算法

本文引入了求解二阶拟线性抛物型微分方程初值问题的一类新的数值算法一分层方法,这种数值方法是通过弱显式欧拉法离散其方程解的概率表示而得到的,相应地给出了该分层方法的收敛性结果.此外,还构造了基于插值的数值算法,最后提供了数值实验,得到的数值结果验证了获得的算法的精确性和有效性.     

Abel变换数值反演的积分算子方法

本文研究了Abel变换的数值反演问题.利用Abel变换的理论反演公式与数值求导的积分算子法相结合的方法,对反演公式中奇异积分合理处理,获得Abel变换数值反演的一种算法,并进行了理论分析与数值实验. 结果表明该算法具有计算简单、数值稳定等优点.     

泊松化图像复原的交替最小化算法

为了快速地去除图像中的泊松噪声,本文在传统的交替方向算法基础上,结合松弛算法提出了一个改进的快速交替最小化算法.与经典的数值算法相比,数值试验表明提出的新算法不但能有效地实现泊松化图像复原,还能大幅度地提高数值计算的速率,并显著地减少电脑的CPU运行时间.     

基于非线性相位恢复的X射线相位衬度断层成像算法

对全息测量下的X射线相位衬度断层成像问题提出了一种新的重建算法.该算法的主要想法是利用牛顿迭代法求解非线性的相位恢复问题.我们证明了牛顿方向满足的线性方程是非适定的,并利用共轭梯度法得到方程的正则化解.最后利用模拟数据进行了数值实验,数值结果验证了算法的合理性以及对噪声数据的数值稳定性,同时通过与线性化相位恢复算法的数值结果比较说明了新算法对探测数据不要求限制在Fresnel区域的近场,适用范围更广.     

基于GMRES(m)法的双连通区域数值保角变换的计算法

本文研究了基于模拟电荷法的双连通区域的数值保角变换问题.利用限制Krylov子空间最大维数的算法–GMRES(m)算法,求解基于模拟电荷法的双连通区域数值保角变换中的约束方程,获得了模拟电荷和变换半径,构造了近似保角变换函数.数值实验表明了本文算法的有效性.     

API stepsize control for the numerical solution of ordinary differential equations

A control-theoretic approach is used to design a new automatic stepsize control algorithm for the numerical integration of ODE's. The new control algorithm is more robust at little extra expense. Its improved performance is particularly evident when the stepsize is limited by numerical stability. Comparative numerical tests are presented.     

An alternative to the Bathe algorithm

This paper presents a new composite sub-steps algorithm for solving reliable numerical responses in structural dynamics. The newly developed algorithm is a two sub-steps, second-order accurate and unconditionally stable implicit algorithm with the same numerical properties as the Bathe algorithm. The detailed analysis of the stability and numerical accuracy is presented for the new algorithm, which shows that its numerical characteristics are identical to those of the Bathe algorithm. Hence, the new sub-steps scheme could be considered as an alternative to the Bathe algorithm. Meanwhile, the new algorithm possesses the following properties: (a) it produces the same accurate solutions as the Bathe algorithm for solving linear and nonlinear problems; (b) it does not involve any artificial parameters and additional variables, such as the Lagrange multipliers; (c) The identical effective stiffness matrices can be obtained inside two sub-steps; (d) it is a self-starting algorithm. Some numerical experiments are given to show the superiority of the new algorithm and the Bathe algorithm over the dissipative CH-α algorithm and the non-dissipative trapezoidal rule.     

AN EFFECT ITERATION ALGORITHM FOR NUMERICAL SOLUTION OF DISCRETE HAMILTON-JACOBIBELLMAN EQUATIONS

§1Introduction ConsidertheHamilton-Jacobi-Bellmanequation max1≤v≤m[A(v)u(x)-f(v)(x)]=0,x∈Ω(1.1)withtheboundarycondition u(x)=0,x∈Ω(1.2)whereΩisabounded,smoothdomaininEuclideanspaceRd,d∈N;f(v)(x)aregiven functionsfromC2(Ω);A(v)aresecond-orderuniformlyellipticoperatorsoftheform A(v)=-d i,j=1a(v)ij2xixj+di=1b(v)ixi+c(v).(1.3)Intheaboveexpression(1.3)therearecoefficientsa(v)ij,b(v)i,c(v)∈C2(Ω)satisfying,forall1≤v≤m,a(v)ij(x)=a(v)ji(x),1≤i,j≤d,c(v)≥c0≥0,x∈Ω,a…     

A comparison of several invariant imbedding algorithms for the solution of two-point boundary-value problems

A comparison of several invariant imbedding algorithms for the numerical solution of two-point boundary-value problems is presented. These include the Scott algorithm, the Kagiwada-Kalaba algorithm, the addition formulas, and the sweep method. Advantages and disadvantages of each algorithm are discussed, and numerical examples are presented.     

A residual algorithm for finding a fixed point of a nonexpansive mapping

This paper presents a weak convergence residual algorithm for finding a fixed point of a nonexpansive mapping in a real Hilbert space. To study the numerical behavior of the algorithm it is included an extensive series of numerical experiments. Our computational experiments show that the new algorithm is computationally efficient.     

求解排课表问题的一个启发式数值算法

为了便于建立与排课表问题有关的决策支持系统,本文给出一个求解排课表问题的启发式数值算法,证明了算法的正确性,并举例说明算法的应用.该算法具有易于编程实现、收敛性好等优点,大量数值实值表明该算法能求得问题的最优解.     

A note on the numerical solution of linear recurrence relations

Summary A numerical algorithm for the computation of non-dominant solutions of linear recurrence relations is analysed. Several non-trivial improvements are made and the efficiency of the new algorithm is illustrated by means of some numerical examples.     

A new algorithm for pole assignment of single-input linear systems using state feedback

In this paper we present a new algorithm for the single-input pole assignment problem using state feedback. This algorithm is based on the Schur decomposition of the closed-loop system matrix, and the numerically stable unitary transformations are used whenever possible, and hence it is numerically reliable.The good numerical behavior of this algorithm is also illustrated by numerical examples.     

一类单边约束问题的数值方法

本文分别基于原始变分形式与对偶混合变分形式，对一类单边约束问题进行了数值求解，提出了求解离散对偶混合变分问题的Uzawa型算法，并用数值例子验证了算法的有效性．     

A reliable algorithm for solving tenth-order boundary value problems

In this paper we present an efficient numerical algorithm for solving linear and nonlinear boundary value problems with two-point boundary conditions of tenth-order. The differential transform method is applied to construct the numerical solutions. The proposed algorithm avoids the complexity provided by other numerical approaches. Several illustrative examples are given to demonstrate the effectiveness of the present algorithm.