Efficient computation of the exponential operator for large,sparse, symmetric matrices

In this paper we compare Krylov subspace methods with Chebyshev series expansion for approximating the matrix exponential operator on large, sparse, symmetric matrices. Experimental results upon negative‐definite matrices with very large size, arising from (2D and 3D) FE and FD spatial discretization of linear parabolic PDEs, demonstrate that the Chebyshev method can be an effective alternative to Krylov techniques, especially when memory bounds do not allow the storage of all Ritz vectors. We also discuss the sensitivity of Chebyshev convergence to extreme eigenvalue approximation, as well as the reliability of various a priori and a posteriori error estimates for both methods. Copyright © 2000 John Wiley & Sons, Ltd.     

用切比雪夫多项式研究短梁的弯曲计算

考虑剪切效应,利用切比雪夫多项式构造严格满足表面切应力边界条件的轴向位移表达式,建立了短梁弯曲问题的新理论．利用奇异函数把作用在短梁上的复杂外载荷表示为分布载荷,推导出了短梁弯曲时的截面正应力公式及挠曲线表达式．把采用切比雪夫多项式推导出短梁的弯曲计算公式计算结果与弹性理论计算结果进行比较,可知该方法的计算精度较高．研究结果表明：在复杂外载荷作用下,当长高比小于等于6时,剪切变形对梁的弯曲挠度影响较大,而当长高比小于3时,剪切变形对梁的弯曲应力影响较大;因此建议采用切比雪夫多项式方法给出的挠度表达式、弯曲应力进行计算,因为切比雪夫多项式方法不但给出了复杂外载荷作用下梁截面挠度、弯曲应力的计算通式,而且该方法具有计算过程简便、精度高的优点．     

Analysis and algorithms for a regularized cauchy problem arising from a non-linear elliptic PDE for seismic velocity estimation

In the present work we derive and study a non-linear elliptic PDE coming from the problem of estimation of sound speed inside the Earth. The physical setting of the PDE allows us to pose only a Cauchy problem, and hence is ill-posed. However, we are still able to solve it numerically on a long enough time interval to be of practical use. We used two approaches. The first approach is a finite difference time-marching numerical scheme inspired by the Lax–Friedrichs method. The key features of this scheme is the Lax–Friedrichs averaging and the wide stencil in space. The second approach is a spectral Chebyshev method with truncated series. We show that our schemes work because of (i) the special input corresponding to a positive finite seismic velocity, (ii) special initial conditions corresponding to the image rays, (iii) the fact that our finite-difference scheme contains small error terms which damp the high harmonics; truncation of the Chebyshev series, and (iv) the need to compute the solution only for a short interval of time. We test our numerical schemes on a collection of analytic examples and demonstrate a dramatic improvement in accuracy in the estimation of the sound speed inside the Earth in comparison with the conventional Dix inversion. Our test on the Marmousi example confirms the effectiveness of the proposed approach.     

The Multivariate Integer Chebyshev Problem

The multivariate integer Chebyshev problem is to find polynomials with integer coefficients that minimize the supremum norm over a compact set in ℂ ^d . We study this problem on general sets but devote special attention to product sets such as cube and polydisk. We also establish a multivariate analog of the Hilbert–Fekete upper bound for the integer Chebyshev constant, which depends on the dimension of space. In the case of single-variable polynomials in the complex plane, our estimate coincides with the Hilbert–Fekete result.      

Universal optimality of digital nets and lattice designs

This article considers universal optimality of digital nets and lattice designs in a regression model. Based on the equivalence theorem for matrix means and majorization theory,the necessary and sufficient conditions for lattice designs being φp-and universally optimal in trigonometric function and Chebyshev polynomial regression models are obtained. It is shown that digital nets are universally optimal for both complete and incomplete Walsh function regression models under some specified conditions,and are...     

Expansion of analytic functions in q-orthogonal polynomials

A classical result on the expansion of an analytic function in a series of Jacobi polynomials is extended to a class of q-orthogonal polynomials containing the fundamental Askey–Wilson polynomials and their special cases. The function to be expanded has to be analytic inside an ellipse in the complex plane with foci at ±1. Some examples of explicit expansions are discussed.      

配点型区间有限元法

在分析Taylor展开``点逼近'区间有限元法不足的基础上, 提出了基于Chebyshev第一类正交多项式全局逼近目标函数的配点型区间有限元法. 该方法不需要计算目标函数对不确定性变量的灵敏度, 不要求不确定性变量的变化范围为小区间, 并适合求解目标函数为不确定变量非线性函数的情形. 目标函数正交展开式的系数采用Gauss-Chebyshev求积公式得到,故需要在不确定性变量所在区间内配置高斯积分点. 计算目标函数在高斯点的取值是该方法的主要工作量, 当不确定性变量数为m, 并选用高斯十点法进行积分时, 需要对系统进行12$m$次分析. 算例表明, 在其他区间有限元法失效的情况下, 配点型区间有限元法依然能够得到几乎精确的区间界限.      

Transformation of Some Mixed Approximation Problems by Optimization Methods

The present paper considers the approximation of a function in the L 2 -norm by functions in a Haar space subject to a constraint set defined by the sup norm. This problem is examined from developments of the Karush-Kuhn-Tucker theorem for semi-infinite programming as well as from the multiobjective optimization approach. An extreme feasible set is characterized with the help of a version of the Chebyshev alternation theorem.     

微极流体薄膜层通过以滑移速度移动的可渗透无限平板时流体特性变化和热辐射对流动和热传导的影响

微极流体薄膜层通过按滑移速度移动的可渗透无限竖直平板时,研究热辐射对混合对流薄膜层流动和热传导的影响.假定流体粘度和热传导率变化是温度的一个函数.对一些典型的可变参数值,应用Chebyshev谱方法,数值求解流动的控制方程.将所得结果与已发表文献的结果进行比较,结果是一致的.绘出并讨论了可变参数对速度、微旋转速度、温度分布曲线、表面摩擦因数和Nusselt数的影响.