DISCRETE MINUS ONE NORM LEAST-SQUARES FOR THE STRESS FORMULATION OF LINEAR ELASTICITY WITH NUMERICAL RESULTS

This paper studies the discrete minus one norm least-squares methods for the stress formulation of pure displacement linear elasticity in two dimensions. The proposed leastsquares functional is defined as the sum of the L^2 -and H^-l -norms of the residual equations weighted appropriately. The minus one norm in the functional is replaced by the discrete minus one norm and then the discrete minus one norm least-squares methods are analyzed with various numerical results focusing on the finite element accuracy and multigrid convergence performances.     

Preconditioned sor methods for generalized least-squares problems

1.IntroductionThegeneralizedleastsquaresproblem,definedasmin(Ax--b)"W--'(Ax--b),(1.1)xacwhereAERm",m>n,bERm,andWERm'misasymmetricandpositivedefinitematrix,isfrequentlyfoundwhensolvingproblemsinstatistics,engineeringandeconomics.Forexample,wegetgeneralizedleastsquaresproblemswhensolvingnonlinearregressionanalysisbyquasi-likelihoodestimation,imagereconstructionproblemsandeconomicmodelsobtainedbythemaximumlikelihoodmethod(of.[1,21).Paige[3,4]investigatestheproblemexplicitly.Hechangestheorig…     

广义特征值极小扰动问题的一类黎曼共轭梯度法

研究含参数$l$非方矩阵对广义特征值极小扰动问题所导出的一类复乘积流形约束矩阵最小二乘问题.与已有工作不同,本文直接针对复问题模型,结合复乘积流形的几何性质和欧式空间上的改进Fletcher-Reeves共轭梯度法,设计一类适用于问题模型的黎曼非线性共轭梯度求解算法,并给出全局收敛性分析.数值实验和数值比较表明该算法比参数$l=1$的已有算法收敛速度更快,与参数$l=n$的已有算法能得到相同精度的解.与部分其它流形优化相比与已有的黎曼Dai非线性共轭梯度法具有相当的迭代效率,与黎曼二阶算法相比单步迭代成本较低、总体迭代时间较少,与部分非流形优化算法相比在迭代效率上有明显优势.     

求解加权线性最小二乘问题的一类预处理GAOR方法

为了快速求解一类来自加权线性最小二乘问题的2×2块线性系统,本文提出一类新的预处理子用以加速GAOR方法,也就是新的预处理GAOR方法.得到了一些比较结果,这些结果表明当GAOR方法收敛时,新方法比原GAOR方法和之前的一些预处理GAOR方法有更好的收敛性.而且,数值算例也验证了新预处理子的有效性．     

Interior-Point Gradient Method for Large-Scale Totally Nonnegative Least Squares Problems

We study an interior-point gradient method for solving a class of so-called totally nonnegative least-squares problems. At each iteration, the method decreases the residual norm along a diagonally-scaled negative gradient direction with a special scaling. We establish the global convergence of the method and present some numerical examples to compare the proposed method with a few similar methods including the affine scaling method.This author was supported in part by DOE/LANL Contract 03891-99-23This author was supported in part by NSF Grant DMS-0442065     

A novel interpolating element-free Galerkin (IEFG) method for two-dimensional elastoplasticity

Using the interpolating moving least-squares (IMLS) method to obtain the shape function, we present a novel interpolating element-free Galerkin (IEFG) method to solve two-dimensional elastoplasticity problems. The shape function of the IMLS method satisfies the property of Kronecker δ function, then in the meshless methods based on the IMLS method, the essential boundary conditions can applied directly. Based on the Galerkin weak form, we obtain the formulae of the IEFG method for solving two-dimensional elastoplasticity problems. The IEFG method has some advantages, such as simpler formulae and directly applying the essential boundary conditions, over the conventional element-free Galerkin (EFG) method. The results of three numerical examples show that the computational precision of the IEFG method is higher than that of the EFG method.     

一类小样本的统计方法建模及其可视化

针对一类高维小样本数据,利用统计方法的非参数检验与偏最小二乘回归(PLS)构造小样本预测模型,实现基于Wilcoxon秩和检验的变量选择与基于PLS的变量压缩降维.并通过DNA序列分类问题实现基于统计方法的小样本数据建模与可视化,计算结果表明方法对小样本具有可行性、有效性.     

关于“求解加权线性最小二乘问题的一类预处理GAOR方法”一文的注记

在"求解加权线性最小二乘问题的一类预处理GAOR方法"一文中,作者提出了求解加权线性最小二乘问题等价$2\times 2$块线性系统的一类预处理GAOR方法,并给出了几个比较定理来说明新提出预处理GAOR方法的优越性.本文我们将指出该文中几个比较定理的不完善之处和证明的错误之处,并给出正确的证明.     

ON THE BREAKDOWNS OF THE GALERKIN AND LEAST-SQUARES METHODS

1　ＩｎｔｒｏｄｕｃｔｉｏｎＷｅｃｏｎｓｉｄｅｒｌｉｎｅａｒｓｙｓｔｅｍｓｏｆｔｈｅｆｏｒｍＡｘ=ｂ,(1 )ｗｈｅｒｅＡ∈ＣＮ×Ｎｉｓｎｏｎｓｉｎｇｕｌａｒａｎｄｐｏｓｓｉｂｌｙｎｏｎ Ｈｅｒｍｉｔｉａｎ .Ａｍａｊｏｒｃｌａｓｓｏｆｍｅｔｈｏｄｓｆｏｒｓｏｌｖｉｎｇ (1 )ｉｓｔｈｅｃｌａｓｓｏｆＫｒｙｌｏｖｓｕｂｓｐａｃｅｍｅｔｈｏｄｓ (ｓｅｅ[6] ,[1 3]ｆｏｒｏｖｅｒｖｉｅｗｓｏｆｓｕｃｈｍｅｔｈｏｄｓ) ,ｄｅｆｉｎｅｄｂｙｔｈｅｐｒｏｐｅｒｔｉｅｓｘｍ ∈ｘ0 +Ｋｍ(ｒ0 ,Ａ) ;(2 )ｒｍ ⊥Ｌｍ, (3)ｗｈｅ…     

Trust-region and other regularisations of linear least-squares problems

We consider methods for regularising the least-squares solution of the linear system Ax=b. In particular, we propose iterative methods for solving large problems in which a trust-region bound ‖x‖≤Δ is imposed on the size of the solution, and in which the least value of linear combinations of ‖Ax−b‖₂ ^q and a regularisation term ‖x‖₂ ^p for various p and q=1,2 is sought. In each case, one or more “secular” equations are derived, and fast Newton-like solution procedures are suggested. The resulting algorithms are available as part of the ALAHAD optimization library. This work was partially supported by EPSRC grants EP/E053351/1 and EP/F005369/1.     

Fast solving of weighted pairing least-squares systems

This paper presents a generalization of the “weighted least-squares” (WLS), named “weighted pairing least-squares” (WPLS), which uses a rectangular weight matrix and is suitable for data alignment problems. Two fast solving methods, suitable for solving full rank systems as well as rank deficient systems, are studied. Computational experiments clearly show that the best method, in terms of speed, accuracy, and numerical stability, is based on a special {1, 2, 3}-inverse, whose computation reduces to a very simple generalization of the usual “Cholesky factorization-backward substitution” method for solving linear systems.     

投影寻踪模型的改进及其在城市水资源承载能力预测中的应用

在城市水资源承载能力研究中,偏最小二乘回归方法能有效地处理自变量间多重线性相关性问题,但不能较好地处理因变量与自变量间复杂的非线性问题.投影寻踪神经网络耦合模型是处理非线性问题的有力工具,而且神经网络投影寻踪耦合模型稳健性高,但不能较好地处理自变量间多重线性相关性问题.本文把这两种方法结合在一起,建立了基于偏最小二乘回归的神经网络投影寻踪耦合模型,对城市水资源承载能力进行了预测,并取得了满意效果.     

矩阵方程AXB=D解的研究

给出了矩阵方程AXB=D相容的又一充要条件,同时讨论它的极小范数解、最小二乘解和极小范数最小二乘解,推广了文献[1]和[3]的结论.     

Analysis of split weighted least-squares procedures for pseudo-hyperbolic equations

In this paper, we introduce two novel split weighted least-squares finite element procedures for pseudo-hyperbolic equations arising in the modelling of nerve conduction process. By selecting the weighted least-squares functional properly, each procedure can be split into two independent symmetric positive definite sub-procedures. One of sub-procedures is for the primitive unknown variable, which is the same as the standard Galerkin finite element procedure and the other is for the introduced flux variable. Optimal order error estimates are developed and the numerical example is given to show the efficiency of the introduced schemes.     

半导体器件瞬态模拟的对称正定混合元方法

提出具有对称正定特性的混合元格式求解非稳态半导体器件瞬态模拟问题。提出一个最小二乘混合元方法、一个新的具有分裂和对称正定性质的混合元格式和一个解经典混合元方程的对称正定失窃工格式求解电场位势和电场强度方程;提出一个最小二乘混合元格式求解关于电子与空穴浓度的非稳态对流扩散方程,浓度函数和流函数被同时求解;采用标准的有限元方法求解热传导方程。建立了误差分析理论。     

基于样条变换的PLS回归的非线性结构分析

基于样条变换的:PLS非线性回归模型既吸取了样条函数分段拟合以适应任意曲线连续变化的优点,又借鉴了偏最小二乘回归方法能够有效解决自变量集合高度相关的技术.针对多元加法模型,从理论和仿真试验的角度分别验证了,对于多个独立自变量对单因变量为非线性关系的数据系统,基于样条变换的PLS回归方法不仅能够有效实现自变量对因变量的整体预测,而且能够提取各维自变量对因变量的单独非线性作用特征,从而确定数据系统内部的复杂非线性结构关系,增强了模型的可解释性.     

Block-Projections Algorithms with Blocks Containing Mutually Orthogonal Rows and Columns

In this paper we present an algorithm which, for a given (sparse) matrix, constructs a partition of its set of row-indices, such that each subset of this partition (except the last one obtained) contains indices which correspond to mutually orthogonal rows. We then use such decompositions in some classes of block-projections methods, previously extended by the author to general inconsistent linear least-squares problems. Numerical experiments on an inconsistent and rank-deficient least-squares model problem are described in the last section of the paper.     

Optimal finite-dimensional solution for a class of nonlinear observation problems

The filtering problem in a differential system with linear dynamics and observations described by an implicit equation linear in the state is solved in finite-dimensional recursive form. The original problem is posed as a deterministic fixed-interval optimization problem (FIOP) on a finite time interval. No stochastic concepts are used. Via Pontryagin's principle, the FIOP is converted into a linear, two-point boundary-value problem. The boundary-value problem is separated by using a linear Riccati transformation into two initial-value problems which give the equations for the optimal filter and filter gain. The optimal filter is linear in the state, but nonlinear with respect to the observation. Stability of the filter is considered on the basis of a related properly linear system. Three filtering examples are given.     

Experimental comparison of least-squares and maximum likelihood methods in factor analysis

Three popular methods to estimate the unknown parameters in the factor analysis model, simple (SLS) and weighted (WLS) least-squares methods and the maximum likelihood method (ML), are compared by a Monte Carlo study. The experiments were conducted with 200 replications for every combination of levels of the following three conditions: method (3 levels), sample size (3 levels) and uniquenesses (2 levels). It was found that SLS performed most favorably when the sample size is relatively small and unique variances are relatively large. WLS and ML proved to be rather alike.     

A remark on split least-squares mixed element procedures for pseudo-parabolic equations

In this paper, we introduce two novel split least-squares mixed element procedures for pseudo-parabolic equations. By selecting the least-squares functional properly, each procedure can be split into two independent symmetric positive definite sub-procedures. One of sub-procedures is for the primitive unknown variable u, which is the same as the standard Galerkin finite element procedure and the other is for the introduced flux variable σ. Optimal order error estimates are developed. A numerical example is given to show the efficiency of the introduced schemes.