用边界曲线构造C^1Coons曲面确定扭矢的方法

本文讨论了由四条边界曲线构造C^1Coons曲面的问题，给出了确定角点扭矢的新方法．该方法沿四边形两对角线方向构造两条四次多项式曲线，每个角点处的扭矢，由一条四次曲线和两条边界曲线确定．跨界切矢由三次埃尔米特插值方法定义．文中还给出了一个用新方法构造曲面的实例．     

用边界曲线构造C~1 Coons曲面确定扭矢的方法

本文讨论了由四条边界曲线构造Ｃ１Ｃｏｏｎｓ曲面的问题，给出了确定角点扭矢的新方法．该方法沿四边形两对角线方向构造两条四次多项式曲线，每个角点处的扭矢，由一条四次曲线和两条边界曲线确定．跨界切矢由三次埃尔米特插值方法定义．文中还给出了一个用新方法构造曲面的实例．     

求解块三对角方程组的新算法

根据块三对角矩阵的特殊分解,给出了求解块三对角方程组的新算法.该算法含有可以选择的参数矩阵,适当选择这些参数矩阵,可以使得计算精度较著名的追赶法高,甚至当追赶法失效时,由该算法仍可得到一定精度的解.     

求三对角矩阵特征多项式一种简便方法

本文通过递推关系 ,直接给出求三对角矩阵特征多项式的一种简便方法 .该方法具有操作简单 ,计算量小的特点 .并给出算例 .     

由三个特殊次序向量对构造三对角对称矩阵

讨论利用给定的三个特殊次序向量对构造不可约三对角矩阵、Jacobi矩阵和负Jacobi矩阵的反问题.在求解方法中,将已知的一些关系式等价地转化为线性方程组,利用线性方程组有解的条件,得到了所研究问题有惟一解的充要条件,并给出了数值算法和例子.     

求三对角矩阵特征多项式一种简便方法

本文通过递推关系，直接给出求三角矩阵特征多项式的一种简便方法。该方法具有操作简单，计算量小的特点。并给出算例。     

周期三对角矩阵逆的一种新算法

给出了一类周期三对角矩阵逆的新的递归算法.新方法充分利用周期三对角矩阵的结构特点,采用递归方法将高阶周期三对角矩阵求逆转化为低阶周期三对角矩阵的求逆.并同时得到简化的计算方法,方法可以有效地减少运算量和存储量,计算精度也有明显的优势.数值实验表明此算法是有效的.     

求分块三对角矩阵和分块周期三对角矩阵逆矩阵的快速算法

给出了分块三对角矩阵逆矩阵的快速算法,并利用所给算法得到了求分块周期三对角矩阵逆矩阵的快速算法.最后通过算例表示算法的有效性.     

一种负高斯曲率的双三次Bézier曲面

负高斯曲率的曲面在建筑业、造船业及航空工业中都有广泛的应用。该文提出一种负高斯曲率双三次Ｂｅｚｉｅｒ曲面的构造方法,简单直观,几何意义明显。     

三对角与五对角Toeplitz矩阵求逆的算法

提出了一种求三对角与五对角Toeplitz矩阵逆的快速算法,其思想为先将Toeplitz矩阵扩展为循环矩阵,再快速求循环矩阵的逆,进而运用恰当矩阵分块求原Toeplitz矩阵的逆的算法.算法稳定性较好且复杂度较低.数值例子显示了算法的有效性和稳定性,并指出了算法的适用范围.     

Biquadratic finite volume element method based on optimal stress points for second order hyperbolic equations

Based on optimal stress points, we develop a full discrete finite volume element scheme for second order hyperbolic equations using the biquadratic elements. The optimal order error estimates in L^∞(H¹), L^∞(L²) norms are derived, in addition, the superconvergence of numerical gradients at optimal stress points is also discussed. Numerical results confirm the theoretical order of convergence. © 2012 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 2013     

An efficient trust region method for unconstrained discrete-time optimal control problems

Discrete-time optimal control (DTOC) problems are large-scale optimization problems with a dynamic structure. In previous work this structure has been exploited to provide very fast and efficient local procedures. Two examples are the differential dynamic programming algorithm (DDP) and the stagewise Newton procedure—both require onlyO(N) operations per iteration, whereN is the number of timesteps. Both exhibit a quadratic convergence rate. However, most algorithms in this category do not have a satisfactory global convergence strategy. The most popular global strategy is shifting: this sometimes works poorly due to the lack of automatic adjustment to the shifting element.In this paper we propose a method that incorporates the trust region idea with the local stagewise Newton's method. This method possesses advantages of both the trust region idea and the stagewise Newton's method, i.e., our proposed method has strong global and local convergence properties yet remains economical. Preliminary numerical results are presented to illustrate the behavior of the proposed algorithm. We also collect in the Appendix some DTOC problems that have appeared in the literature.Partially supported by the Cornell Theory Center, which receives major funding from the National Science Foundation and IBM Corporation, with additional support from the State of New York and its Corporate Research Institutes; and by NSF, AFOSR, and ONR through grant DMS-8920550.     

Criteria for farthest points on convex surfaces

We provide a sharp, sufficient condition to decide if a point y on a convex surface S is a farthest point (i.e., is at maximal intrinsic distance from some point) on S, involving a lower bound π on the total curvature ω_y at y, ω_y ≥ π. Further consequences are obtained when ω_y > π, and sufficient conditions are derived to guarantee that a convex cap contains at least one farthest point. A connection between simple closed quasigeodesics O of S, points y ∈ S\O with ω_y > π, and the set ?? of all farthest points on S, is also investigated (© 2009 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

A new method for constructing infinite families of k-tight optimal double loop networks

The double loop network (DLN) is a circulant digraph with n nodes and outdegree 2. DLN has been widely used in the designing of local area networks and distributed systems. In this paper, a new method for constructing infinite families of k-tight optimal DLN is presented. For k = 0,1,…,40, the infinite families of k-tight optimal DLN can be constructed by the new method, where the number nk(t,a) of their nodes is a polynomial of degree 2 in t and contains a parameter a. And a conjecture is proposed.     

An optimal bound for inverses of triangular matrices with monotone entries

This article provides a new bound for 1-norms of inverses of positive triangular matrices with monotonic column entries. The main theorem refines a recent inequality established in Vecchio and Mallik [Bounds on the inverses of non-negative lower triangular Toeplitz matrices with monotonicity properties, Linear Multilinear Alg., 55 (2007), pp. 365–379]. The results are shown to be in a sense best possible under the given constraints.     

Computing eigenvectors of block tridiagonal matrices based on twisted block factorizations

New methods for computing eigenvectors of symmetric block tridiagonal matrices based on twisted block factorizations are explored. The relation of the block where two twisted factorizations meet to an eigenvector of the block tridiagonal matrix is reviewed. Based on this, several new algorithmic strategies for computing the eigenvector efficiently are motivated and designed. The underlying idea is to determine a good starting vector for an inverse iteration process from the twisted block factorizations such that a good eigenvector approximation can be computed with a single step of inverse iteration.     

A penalty function method for solving inverse optimal value problem

In order to consider the inverse optimal value problem under more general conditions, we transform the inverse optimal value problem into a corresponding nonlinear bilevel programming problem equivalently. Using the Kuhn–Tucker optimality condition of the lower level problem, we transform the nonlinear bilevel programming into a normal nonlinear programming. The complementary and slackness condition of the lower level problem is appended to the upper level objective with a penalty. Then we give via an exact penalty method an existence theorem of solutions and propose an algorithm for the inverse optimal value problem, also analysis the convergence of the proposed algorithm. The numerical result shows that the algorithm can solve a wider class of inverse optimal value problem.     

An alternative approach for non-linear optimal control problems based on the method of moments

We propose an alternative method for computing effectively the solution of non-linear, fixed-terminal-time, optimal control problems when they are given in Lagrange, Bolza or Mayer forms. This method works well when the nonlinearities in the control variable can be expressed as polynomials. The essential of this proposal is the transformation of a non-linear, non-convex optimal control problem into an equivalent optimal control problem with linear and convex structure. The method is based on global optimization of polynomials by the method of moments. With this method we can determine either the existence or lacking of minimizers. In addition, we can calculate generalized solutions when the original problem lacks of minimizers. We also present the numerical schemes to solve several examples arising in science and technology.     

GSOR方法最优迭代参数的简化推导

对于增广线性系统,Bai等研究了广义SOR方法（Bai Z Z,Parlett B,Wang Z Q.On generaliged successive overrelaxation methods for augmented linear systems.Numerische Mathematik,2005,102（1）：1-38）,并得到其最优迭代参数.给出了另外一种推导最优迭代参数的简化方法,这种方法对于求解其他参数加速定常迭代方法的最优迭代参数非常有意义.