共查询到20条相似文献,搜索用时 31 毫秒
1.
R. K. Beatson
G. N. Newsam
《Computers & Mathematics with Applications》1992,24(12):7-19This paper describes some new techniques for the rapid evaluation and fitting of radial basic functions. The techniques are based on the hierarchical and multipole expansions recently introduced by several authors for the calculation of many-body potentials. Consider in particular the N term thin-plate spline, s(x) = Σj=1N djφ(x−xj), where φ(u) = |u|2log|u|, in 2-dimensions. The direct evaluation of s at a single extra point requires an extra O(N) operations. This paper shows that, with judicious use of series expansions, the incremental cost of evaluating s(x) to within precision ε, can be cut to O(1+|log ε|) operations. In particular, if A is the interpolation matrix, ai,j = φ(xi−xj, the technique allows computation of the matrix-vector product Ad in O(N), rather than the previously required O(N2) operations, and using only O(N) storage. Fast, storage-efficient, computation of this matrix-vector product makes pre-conditioned conjugate-gradient methods very attractive as solvers of the interpolation equations, Ad = y, when N is large. 相似文献
2.
In this paper, we investigate the existence and uniqueness of cardinal interpolants associated with functions arising from the kth order iterated discrete Laplacian k applied to certain radial basis functions. In particular, we concentrate on determining, for a given radial function Φ, which functions kΦ give rise to cardinal interpolation operators which are both bounded and invertible ℓ2 (Z3). In addition to solving the cardinal interpolation problem (CIP) associated with such functions kΦ, our approach provides a unified framework and simpler proofs for the CIP associated with polyharmonic splines and Hardy multiquadrics. 相似文献
3.
4.
We show that it is decidable for any given ground term rewrite systems R and S if there is a ground term rewrite system U such that ↔U*=↔R*∩↔S*. If the answer is yes, then we can effectively construct such a ground term rewrite system U. In other words, for any given finitely generated congruences ρ and τ over the term algebra, it is decidable if ρ∩τ is a finitely generated congruence. If the answer is yes, then we can effectively construct a ground term rewrite system U such that ↔U*=ρ∩τ. 相似文献
5.
Let
be such that d1,pd=1,p02 and
. We are proving in this note a new criterion for the pair
to be a canonical number system. This enables us to prove that if p2,…,pd−1,∑i=1dpi0 and p0>2∑i=1d|pi|, then
is a canonical number system. 相似文献
6.
V. Paramasivam S. P. Thulasi Ramachandre Gowda D. Muthiah Raj 《Computers & Structures》1992,43(6):1187-1189
In situations where transverse shear deformations and rotary inertia in beams are important, elements based on the Timoshenko beam theory are useful. Among the two-noded, four DOF elements derived from the minimum total potential energy principle, the HTK. element proposed by Hughes et al. using linear displacement functions for both w and θ and the T1CC4 element proposed by Tessler et al. using quadratic displacement function for w and linear displacement function for θ are well known in the literature. The convergence of the HTK element in the thin beam situation has been too poor due to shear locking but by using selective integration this element can be shown to be equivalent to the T1CC4 element which has a rate of convergence of O(h2). In this paper a five DOF element with w and θ at the end nodes and θ at the middle node and based on the cubic displacement function for w and the quadratic displacement function for θ is first developed. Statically condensing the middle rotational DOF, the well-known (4 × 4) stiffness matrix using the φ-factor defined as φ = 12EI/kGAL2 and hitherto obtained only through a flexibility approach or closed-form solution of the governing equations of the Timoshenko beam theory is derived. This element based on cubic displacement function for w has rate of convergence of O(h4), is completely free of shear locking and performs equally well in thin as well as thick beam situations. 相似文献
7.
In this paper, the existence of at least three positive solutions for the boundary value problem (BVP) of second-order functional differential equation with the form y″(t) + f(t, yt) = 0, for t ε [0,1], y(t) -βy′(t) =η(t), for t ε [−τ,0], −γy(t) + Δy′(t) = ζ(t), for t ε [1, 1 + a], is studied. Moreover, we investigate the existence of at least three partially symmetric positive solutions for the above BVP with Δ = βγ. 相似文献
8.
Ine-Wei Liu 《Computers & Structures》1996,59(6):1051-1058
A conforming finite element formulation of the equations governing composite multilayered plates using Reddy's higher-order theory is presented. The element has eight degrees of freedom, u0, v0, w, ∂w/∂x, ∂w/∂y, ∂2w/∂x∂y, γx, γy, per node. The transverse displacement of the present element is described by a modified bicubic displacement function while the in-plane displacements and shear-rotations are interpolated quadraticly. The element is evaluated for its accuracy in the analysis of static, vibration, and buckling of anisotropic rectangular plates with different lamination schemes and boundary conditions. The conforming finite element described here for the higher-order theory gives fairly accurate results for displacements, stresses, buckling loads, and natural frequencies. 相似文献
9.
10.
This paper presents an efficient algorithm for enumerating all minimal a-b separators separating given non-adjacent vertices a and b in an undirected connected simple graph G = (V, E), Our algorithm requires O(n3Rab) time, which improves the known result of O(n4Rab) time for solving this problem, where ¦V¦= n and Rab is the number of minimal a-b separators. The algorithm can be generalized for enumerating all minimal A-B separators that separate non-adjacent vertex sets A, B < V, and it requires O(n2(n − nA − nb)RAB) time in this case, where na = ¦A¦, nB = ¦B¦ and rAB is the number of all minimal A−B separators. Using the algorithm above as a routine, an efficient algorithm for enumerating all minimal separators of G separating G into at least two connected components is constructed. The algorithm runs in time O(n3R+Σ + n4RΣ), which improves the known result of O(n6RΣ) time, where Rσ is the number of all minimal separators of G and RΣR+Σ = ∑1i, vj) ERvivj (n(n − 1)/2 − m)RΣ. Efficient parallelization of these algorithms is also discussed. It is shown that the first algorithm requires at most O((n/log n)Rab) time and the second one runs in time O((n/log n)R+Σ+n log nRΣ) on a CREW PRAM with O(n3) processors. 相似文献
11.
Mean-value theorems with sharp quantitative remainder term estimates and theorems on the distribution of values are proved for a large class of multiplicative functions. This class is characterized by the existence of an abscissa of absolute convergence ≤ 1 for the Dirichlet series of the quotient , where is the Dirichlet series associated with ƒ and ζ(s) denotes the Riemann zeta function. In particular, the mean behaviour of ƒ on certain sequences a, namely, that of Mersenne numbers 2n −1 and that of shifted primes p−1 with natural n and prime p, respectively, is studied and cluster points of the image sequence f(a) of a are determined. The strength of the results is demonstrated by new results for some specific functions. 相似文献
12.
T. -W. Hwang T. -H. Kuo C. -A. Wang 《Computers & Mathematics with Applications》1989,17(12):1573-1586
We study a two-point boundary values problem (TPBVP) T + Q(Aƒƒ″−(ƒ′)2)=β, ƒ(0)=ƒ(1)=ƒ′(1)=ƒPrime;(0)+1=0, 1A<∞ which arises from a study of similarity transformation for surface-tension drive flows of low Prandtl number fluids in a slot with an insulated bottom. Numerical solutions of the (TPBVP) are found for cases of assumed linear (A = 2) and quadratic (A = 1) radiation of temperature along the free surface. The corresponding bifurcation diagrams in terms of Q and β indicate that multiple solutions occur at some Q when A = 1 and unique solution exists for each Q when A = 2. Some existing properties of the TPBVP when A = 1 and 2, are also verified mathematically. 相似文献
13.
Jae Hoon Kong Seung Pil Jeong Sunhong Lee Gwang Il Kim 《Computer Aided Geometric Design》2008,25(4-5):214-229
We introduce a new method of solving C1 Hermite interpolation problems, which makes it possible to use a wider range of PH curves with potentially better shapes. By characterizing PH curves by roots of their hodographs in the complex representation, we introduce PH curves of type K(t−c)2n+1+d. Next, we introduce a speed reparametrization. Finally, we show that, for C1 Hermite data, we can use PH curves of type K(t−c)2n+1+d or strongly regular PH quintics satisfying the G1 reduction of C1 data, and use these curves to solve the original C1 Hermite interpolation problem. 相似文献
14.
In (Röschel, l997) B-spline technique was used for blending of Lagrange interpolants. In this paper we generalize this idea replacing Lagrange by Hermite interpolants. The generated subspline b(t) interpolates the Hermite input data consisting of parameter values ti and corresponding derivatives ai,j, j=0,…,i−1, and is called blended Hermite interpolant (BHI). It has local control, is connected in affinely invariant way with the input and consists of integral (polynomial) segments of degree 2·k−1, where k−1max{i}−1 denotes the degree of the B-spline basis functions used for the blending. This method automatically generates one of the possible interpolating subsplines of class Ck−1 with the advantage that no additional input data is necessary. 相似文献
15.
Ruud van Damme 《Computer Aided Geometric Design》1997,14(9):847-875
A subdivision scheme for constructing smooth surfaces interpolating scattered data in R3 is proposed. It is also possible to impose derivative constraints in these points. In the case of functional data, i.e., data are given in a properly triangulated set of points {(xi, yi)}i=1N from which none of the pairs (xi,yi) and (xj,yj) with i≠j coincide, it is proved that the resulting surface (function) is C1. The method is based on the construction of a sequence of continuous splines of degree 3. Another subdivision method, based on constructing a sequence of splines of degree 5 which are once differentiable, yields a function which is C2 if the data are not ‘too irregular’. Finally the approximation properties of the methods are investigated. 相似文献
16.
We formulate a class of difference schemes for stiff initial-value problems, with a small parameter ε multiplying the first derivative. We derive necessary conditions for uniform convergence with respect to the small parameter ε, that is the solution of the difference scheme uih satisfies |uih−u(xi)| Ch, where C is independent of h and ε. We also derive sufficient conditions for uniform convergence and show that a subclass of schemes is also optimal in the sense that |uih−u(xi)| C min (h, ε). Finally, we show that this class contains higher-order schemes. 相似文献
17.
Hung Q. Ngo 《Theoretical computer science》2003,290(3):2157-2167
In this paper, we study the problem of finding routing algorithms on the multirate rearrangeable Clos networks which use as few number of middle-stage switches as possible. We propose a new routing algorithm called the “grouping algorithm”. This is a simple algorithm which uses fewer middle-stage switches than all known strategies, given that the number of input-stage switches and output-stage switches are relatively small compared to the size of input and output switches. In particular, the grouping algorithm implies that m = 2n+(n−1)/2k is a sufficient number of middle-stage switches for the symmetric three-stage Clos network C(n,m,r) to be multirate rearrangeable, where k is any positive integer and rn/(2k−1). 相似文献
18.
Hüseyin Akay 《Systems & Control Letters》2000,39(5):581-376
In this paper, model sets for linear time-invariant systems spanned by fixed pole orthonormal bases are investigated. The obtained model sets are shown to be complete in Lp(T) (1<p<∞), the Lebesque spaces of functions on the unit circle T, and in C(T), the space of periodic continuous functions on T. The Lp norm error bounds for estimating systems in Lp(T) by the partial sums of the Fourier series formed by the orthonormal functions are computed for the case 1<p<∞. Some inequalities on the mean growth of the Fourier series are also derived. These results have application in estimation and model reduction. 相似文献
19.
I. Shevchenko 《Computers & Mathematics with Applications》2004,47(12):1827-1855
To approach a simple game Δ2 of P and E = {E1, E2} with no a priori evaders' role assignment and the payoff equal to the distance to one evader at an instant of catching another, we introduce a concept of casting and study the games Δ1,2 and Δ2,1 for preassigned and Δp2 for open-loop casting procedures. Since Δp2 is reduced to Δ1,2 or Δ2,1 which, in turn, are distinguished only by their notations, we focus attention mainly on Δ1,2. According to the tenet of transition, Δ1,2 is divided into a concatenation of Δ1,2b (basic) and Δ1,2a (auxiliary) games that model the problem before and after the first instant of E1 capture. The games Δ1,2a, Δ1,2b, Δ1,2 are studied one after another with use of the Isaacs' approach extended by Berkowitz, Breakwell, Bernhard et al. 相似文献
20.
A solution X of a discrete-time algebraic Riccati equation is called unmixed if the corresponding closed-loop matrix Φ(X) has the property that the common roots of det(sI−Φ(X)) and det(I−sΦ(X)*) (if any) are on the unit circle. A necessary and sufficient condition is given for existence and uniqueness of an unmixed solution such that the eigenvalues of Φ(X) lie in a prescribed subset of
. 相似文献