基于径向基函数逼近的非线性动力系统数值求解

径向基函数具有形式简单、各向同性等优点.将径向基函数逼近的思想与加权余量配点法相结合,借鉴边值问题的求解,构造了一种求解非线性动力系统初值问题的数值方法.分析了几种较为成熟的非线性动力系统数值求解方法的优缺点.给出了实际算例,与已有方法对比,表明该方法计算过程简单、收敛性好、计算精度高.     

Constructive Approximation by Superposition of Sigmoidal Functions

In this paper, a constructive theory is developed for approximating functions of one or more variables by superposition of sigmoidal functions. This is done in the uniform norm as well as in the $L^p$ norm. Results for the simultaneous approximation, with the same order of accuracy, of a function and its derivatives (whenever these exist), are obtained. The relation with neural networks and radial basis functions approximations is discussed. Numerical examples are given for the purpose of illustration.     

Surface Reconstruction of 3D Scattered Data with Radial Basis Functions

We use Radial Basis Functions （RBFs） to reconstruct smooth surfaces from 3D scattered data. An object＇s surface is defined implicitly as the zero set of an RBF fitted to the given surface data. We propose improvements on the methods of surface reconstruction with radial basis functions. A sparse approximation set of scattered data is constructed by reducing the number of interpolating points on the surface. We present an adaptive method for finding the off-surface normal points. The order of the equation decreases greatly as the number of the off-surface constraints reduces gradually. Experimental results are provided to illustrate that the proposed method is robust and may draw beautiful graphics.     

缺失数据下线性模型的半参数估计

通过比较参数方法和非参数方法对选择概率建模的优缺点,基于充分降维的思想提出了一种利用单指标模型对选择概率建模的半参数方法.基于逆概率加权方法和半参数方法,研究了缺失数据下线性模型的统计推断问题.建立的逆概率加权估计方程可以处理不同的数据缺失情形,给出了线性模型中兴趣参数的估计,并证明了它的渐近正态性.最后通过模拟研究说明提出的方法具有较好的有限样本性质.     

Oligopolistic Markets with Leadership and Demand Functions Possibly Discontinuous

An oligopolistic market in which several firms supply a homogeneous product is considered. One particular firm (leader) sets its decision prior to the other firms (followers), which react by playing a competition in quantities. The aim of the paper is to investigate the existence of suitable strategies for the leader in the presence of aggregative coupled constraints in the case in which (i) the aggregate quantity reaction Q is single-valued and (ii) the aggregate quantity reaction Q is set-valued and the leader has some probabilistic information on the followers choice in their reaction set.     

线性流形上矩阵方程AX=B的一类反问题及数值解法

1．引言本文用＊-"m表示全体nX。实矩阵的集合，人表示n阶单位矩阵，汉"m一《ME＊""叫rank（川一r），＊＊"""＝HE＊"""卜"A＝v，＊＊"""一仰E＊"""卜"一M｝，SR；""（SR7""）表示全体7。阶实对称半正定（正定）阵集合．N（A）表示矩阵A的零空间，即N（A）＝（xlAx＝0），ID叫D表示Frobenius范数，A"表示矩阵A的Moors－Penrose广义逆，［EI十表示在Frobenius范数意义下n阶方阵E在SR；""中唯一的最佳k逼近解，即口一［E］＋11-inf。。、。。x，IllE-All．（［E］十求法见文［7］）．还用A三0（A三0）表示A（的k阶顺序主子矩…     

线性流形上的广义中心对称矩阵反问题

设R∈Cn×n是满足R=RH=R-1≠±In的广义反射矩阵.若A∈Cn×n满足RAR=A,则称A为n阶广义中心对称矩阵,n阶广义中心对称矩阵的全体记为GCSCn×n.令X1,Z1∈Cn×k1,Y1,W1∈Cn×l1,S={A|‖AX1-Z1‖2+‖Y1HA-W1H‖2=min,A∈GCSCn×n},本文研究如下问题.问题Ⅰ.给定矩阵Z2,X2∈Cn×k2,Y2,W2∈Cn×l2,求A∈S,使得其中‖·‖是Frobenius范数.问题Ⅱ.给定矩阵A∈Cn×n,求A∈SE,使得其中SE是问题Ⅰ的解集合.本文给出了问题Ⅰ解集合SE的表达式,并导出了矩阵方程AX2=Z2,Y2HA=W2H有解A∈S的充分必要条件及其通解表达式,并给出了问题Ⅱ解的表达式以及求解问题Ⅱ的数值方法和数值例子.     

结构动力响应中急动度的计算

急动度(jerk)在工程实践中具有重要的意义.将径向基函数逼近与配点法相结合,发展了一种能够有效求解动力响应的数值算法.该方法使用径向基函数插值来逼近真实的运动规律,能够用于急动度和急动度(三阶)方程的计算,弥补了传统的数值方法无法计算急动度的不足.并针对微分方程的特点,提出了改进的多变量联合插值函数,同时添加与微分方程同阶的初值条件,可显著减小数值震荡.算例表明,该方法具有计算过程简单、精度高的特点,同时对急动度方程也有很好的适用性.     

On Weighted Approximation by Modified Bernstein Operators for Functions with Singularities

Della Vecchia et al. （see [2]） introduced a kind of modified Bernstein operators which can be used to approximate functions with singularities at endpoints on [0,1]. In the present paper, we obtain a kind of pointwise Stechkin-type inequalities for weighted approximation by the modified Bemsetin operators.     

线性流形上双反对称矩阵的最佳逼近

本文讨论了线性流形上用双反对称矩阵构造给定矩阵的最佳逼近问题,给出问题解的表达式,最后给出求最佳逼近解的数值方法与数值算例.     

Structured FISTA for image restoration

In this paper, we propose an efficient numerical scheme for solving some large‐scale ill‐posed linear inverse problems arising from image restoration. In order to accelerate the computation, two different hidden structures are exploited. First, the coefficient matrix is approximated as the sum of a small number of Kronecker products. This procedure not only introduces one more level of parallelism into the computation but also enables the usage of computationally intensive matrix–matrix multiplications in the subsequent optimization procedure. We then derive the corresponding Tikhonov regularized minimization model and extend the fast iterative shrinkage‐thresholding algorithm (FISTA) to solve the resulting optimization problem. Because the matrices appearing in the Kronecker product approximation are all structured matrices (Toeplitz, Hankel, etc.), we can further exploit their fast matrix–vector multiplication algorithms at each iteration. The proposed algorithm is thus called structured FISTA (sFISTA). In particular, we show that the approximation error introduced by sFISTA is well under control and sFISTA can reach the same image restoration accuracy level as FISTA. Finally, both the theoretical complexity analysis and some numerical results are provided to demonstrate the efficiency of sFISTA.     

Norm estimations for perturbations of the weighted Moore-Penrose inverse

For a complex matrix $A\in \mathbb{C}^{m\times n}$, the relationship between the weighted Moore-Penrose inverse $A^\dag_{M_1N_1}$ and $A^\dag_{M_2N_2}$ is studied, and an important formula is derived,where $M_1\in \mathbb{C}^{m\times m}, N_1\in\mathbb{C}^{n\times n}$ and $M_2\in \mathbb{C}^{m\times m}, N_2\in\mathbb{C}^{n\times n}$ are different pair of positive definite hermitian matrices. Based on this formula, this paper initiates the study of the perturbation estimations for $A^\dag_{MN}$ in the case that $A$ is fixed, whereas both $M$ and $N$ are variable. The obtained norm upper bounds are then applied to the perturbation estimations for the solutions to the weighted linear least squares problems.     

Numerical approximation of the spectra of Phan-Thien Tanner liquids

The linear stability of the linear Phan-Thien Tanner (PTT) fluid model is investigated for plane Poiseuille flow. The PTT model involves parameters that can be used to fit shear and extensional data, which makes it suitable for describing both polymer solutions and melts. The base flow is determined using a Chebyshev-tau method. The linear stability equations are also discretized using Chebyshev approximations to furnish a generalized eigenvalue problem. The spectrum is shown to comprise a continuous part and a discrete part. The theoretical and numerical results are validated for the UCM and Oldroyd-B models, which are special cases of the PTT model, by comparing with results in the literature. It is demonstrated that the linear extensional and elasticity parameters considered. The computational efficiency and accuracy of the numerical method are also investigated.     

Boundary Behavior of Harmonic Functions on Manifolds

In this paper, we mainly set up a kind of representation theorem of harmonic functions on manifolds with Ricci curvature bounded below and study non-tangential limits of harmonic functions.     

Numerical Approximation of the Spectra of Phan-Thien Tanner Liquids

The linear stability of the linear Phan-Thien Tanner (PTT) fluid model is investigated for plane Poiseuille flow. The PTT model involves parameters that can be used to fit shear and extensional data, which makes it suitable for describing both polymer solutions and melts. The base flow is determined using a Chebyshev-tau method. The linear stability equations are also discretized using Chebyshev approximations to furnish a generalized eigenvalue problem. The spectrum is shown to comprise a continuous part and a discrete part. The theoretical and numerical results are validated for the UCM and Oldroyd-B models, which are special cases of the PTT model, by comparing with results in the literature. It is demonstrated that the linear PTT fluid is stable to infinitesimal disturbances with respect to the range of shear-thinning, extensional and elasticity parameters considered. The computational efficiency and accuracy of the numerical method are also investigated.     

线性流形上中心对称矩阵的最佳逼近

1 引 言令Rn×m表示所有n×m阶实矩阵集合;ORn×n表示所有n×n阶正交矩阵之集;A+表示矩阵A的Moore-Penrose广义逆;Iκ表示κ阶单位阵;||·||表示矩阵的Frobenius范数;rank（A）表示矩阵A的秩.设ei为n阶单位矩阵In的第i列（i=1,2,…,n）,记Sn=（en,en-1,…,e1）,易知     

Detection of discontinuities in scattered data approximation

A Detection Algorithm for the localisation of unknown fault lines of a surface from scattered data is given. The method is based on a local approximation scheme using thin plate splines, and we show that this yields approximation of second order accuracy instead of first order as in the global case. Furthermore, the Detection Algorithm works with triangulation methods, and we show their utility for the approximation of the fault lines. The output of our method provides polygonal curves which can be used for the purpose of constrained surface approximation. This revised version was published online in June 2006 with corrections to the Cover Date.     

Application of wavelet bases in linear and nonlinear approximation to functions from Besov spaces

Linear and nonlinear approximations to functions from Besov spaces B _{p, q} ^σ ([0, 1]), σ > 0, 1 ≤ p, q ≤ ∞ in a wavelet basis are considered. It is shown that an optimal linear approximation by a D-dimensional subspace of basis wavelet functions has an error of order D ^{-min(σ, σ + 1/2 ? 1/p)} for all 1 ≤ p ≤ ∞ and σ > max(1/p ? 1/2, 0). An original scheme is proposed for optimal nonlinear approximation. It is shown how a D-dimensional subspace of basis wavelet functions is to be chosen depending on the approximated function so that the error is on the order of D ^?σ for all 1 ≤ p ≤ ∞ and σ > max(1/p ? 1/2, 0). The nonlinear approximation scheme proposed does not require any a priori information on the approximated function.     

The inverse problem of bisymmetric matrices with a submatrix constraint

An n×n real matrix A is called a bisymmetric matrix if A=A^T and A=S_nAS_n, where S_n is an n×n reverse unit matrix. This paper is mainly concerned with solving the following two problems: Problem I Given n×m real matrices X and B, and an r×r real symmetric matrix A₀, find an n×n bisymmetric matrix A such that where A([1: r]) is a r×r leading principal submatrix of the matrix A. Problem II Given an n×n real matrix A^*, find an n×n matrix  in S_E such that where ∥·∥ is Frobenius norm, and S_E is the solution set of Problem I. The necessary and sufficient conditions for the existence of and the expressions for the general solutions of Problem I are given. The explicit solution, a numerical algorithm and a numerical example to Problem II are provided. Copyright © 2003 John Wiley & Sons, Ltd.