首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 579 毫秒
1.
In this paper we study quadrature formulas of the types (1) $$\int\limits_{ - 1}^1 {(1 - x^2 )^{\lambda - 1/2} f(x)dx = C_n^{ (\lambda )} \sum\limits_{i = 1}^n f (x_{n,i} ) + R_n \left[ f \right]} ,$$ (2) $$\int\limits_{ - 1}^1 {(1 - x^2 )^{\lambda - 1/2} f(x)dx = A_n^{ (\lambda )} \left[ {f\left( { - 1} \right) + f\left( 1 \right)} \right] + K_n^{ (\lambda )} \sum\limits_{i = 1}^n f (\bar x_{n,i} ) + \bar R_n \left[ f \right]} ,$$ with 0<λ<1, and we obtain inequalities for the degreeN of their polynomial exactness. By using such inequalities, the non-existence of (1), with λ=1/2,N=n+1 ifn is even andN=n ifn is odd, is directly proved forn=8 andn≥10. For the same value λ=1/2 andN=n+3 ifn is evenN=n+2 ifn is odd, the formula (2) does not exist forn≥12. Some intermediary results regarding the first zero and the corresponding Christoffel number of ultraspherical polynomialP n (λ) (x) are also obtained.  相似文献   

2.
L. Rebolia 《Calcolo》1973,10(3-4):245-256
The coefficientsA hi and the nodesx mi for «closed” Gaussian-type quadrature formulae $$\int\limits_{ - 1}^1 {f(x)dx = \sum\limits_{h = 0}^{2_8 } {\sum\limits_{i = 0}^{m + 1} {A_{hi} f^{(h)} (x_{mi} ) + R\left[ {f(x)} \right]} } } $$ withx m0 =?1,x m, m+1 =1 andR[f(x)]=0 iff(x) is a polinomial of degree at most2m(s+1)+2(2s+1)?1, have been tabulated for the cases: $$\left\{ \begin{gathered} s = 1,2 \hfill \\ m = 2,3,4,5 \hfill \\ \end{gathered} \right.$$ .  相似文献   

3.
F. Costabile 《Calcolo》1974,11(2):191-200
For the Tschebyscheff quadrature formula: $$\int\limits_{ - 1}^1 {\left( {1 - x^2 } \right)^{\lambda - 1/2} f(x) dx} = K_n \sum\limits_{k = 1}^n {f(x_{n,k} )} + R_n (f), \lambda > 0$$ it is shown that the degre,N, of exactness is bounded by: $$N \leqslant C(\lambda )n^{1/(2\lambda + 1)} $$ whereC(λ) is a convenient function of λ. For λ=1 the complete solution of Tschebyscheff's problem is given.  相似文献   

4.
P. Brianzi  L. Rebolia 《Calcolo》1982,19(1):71-86
A numerical performance of integral form for the linear ordinary differential equations $$y^{(n)} = \sum\limits_{i = 0}^{n - 2} { a_{i + 2} (x) y^{(n - 2 - i)} (x)}$$ is proved. Three numerical experiments are also given.  相似文献   

5.
H. H. Gonska  J. Meier 《Calcolo》1984,21(4):317-335
In 1972 D. D. Stancu introduced a generalization \(L_{mp} ^{< \alpha \beta \gamma > }\) of the classical Bernstein operators given by the formula $$L_{mp}< \alpha \beta \gamma > (f,x) = \sum\limits_{k = 0}^{m + p} {\left( {\begin{array}{*{20}c} {m + p} \\ k \\ \end{array} } \right)} \frac{{x^{(k, - \alpha )} \cdot (1 - x)^{(m + p - k, - \alpha )} }}{{1^{(m + p, - \alpha )} }}f\left( {\frac{{k + \beta }}{{m + \gamma }}} \right)$$ . Special cases of these operators had been investigated before by quite a number of authors and have been under investigation since then. The aim of the present paper is to prove general results for all positiveL mp <αβγ> 's as far as direct theorems involving different kinds of moduli of continuity are concerned. When applied to special cases considered previously, all our corollaries of the general theorems will be as good as or yield improvements of the known results. All estimates involving the second order modulus of continuity are new.  相似文献   

6.
L. Rebolia 《Calcolo》1965,2(3):351-369
The coefficientsA hi (m,s) and the nodesx i (m,s) for Gaussian-type quadrature formulae
$$\int\limits_{ - 1}^1 {f(x)dx = \mathop \sum \limits_{h = 0}^{2s} \mathop \sum \limits_{i = 1}^m } A_{hi} \cdot f^{(h)} (x_i )$$  相似文献   

7.
We show in this note that the equation αx1 + #x22EF; +αxp?ACβy1 + α +βyq where + is an AC operator and αx stands for x+...+x (α times), has exactly $$\left( { - 1} \right)^{p + q} \sum\limits_{i = 0}^p {\sum\limits_{j = 0}^q {\left( { - 1} \right)^{1 + 1} \left( {\begin{array}{*{20}c} p \\ i \\ \end{array} } \right)\left( {\begin{array}{*{20}c} q \\ j \\ \end{array} } \right)} 2^{\left( {\alpha + \begin{array}{*{20}c} {j - 1} \\ \alpha \\ \end{array} } \right)\left( {\beta + \begin{array}{*{20}c} {i - 1} \\ \beta \\ \end{array} } \right)} } $$ minimal unifiers if gcd(α, β)=1.  相似文献   

8.
P. Baratella 《Calcolo》1977,14(3):237-242
In this paper we study the remainder term of a quadrature formula of the form $$\int\limits_{ - 1}^1 {f(x)dx = A_n \left[ {f( - 1) + f(1)} \right] + C_n \sum\limits_{i = 1}^n {f(x_{n,i} ) + R_n \left[ f \right],} } $$ , withx x,i -1,1, andR n [f]=0 whenf(x) is a polynomial of degree ≤n+3 ifn is even, or ≤n+2 ifn is odd. Such a formula exists only forn=1(1)11. It is shown that, iff(x)∈ C(h+1) [-1,1], (h=n+3 orn+2), thenR n [f]=f h+1 (τ)·± n . The values α n are given.  相似文献   

9.
In this paper we construct an interpolatory quadrature formula of the type $$\mathop {\rlap{--} \smallint }\limits_{ - 1}^1 \frac{{f'(x)}}{{y - x}}dx \approx \sum\limits_{i = 1}^n {w_{ni} (y)f(x_{ni} )} ,$$ wheref(x)=(1?x)α(1+x)β f o(x), α, β>0, and {x ni} are then zeros of then-th degree Chebyshev polynomial of the first kind,T n (x). We also give a convergence result and examine the behavior of the quantity \( \sum\limits_{i = 1}^n {|w_{ni} (y)|} \) asn→∞.  相似文献   

10.
Recently, we derived some new numerical quadrature formulas of trapezoidal rule type for the integrals \(I^{(1)}[g]=\int ^b_a \frac{g(x)}{x-t}\,dx\) and \(I^{(2)}[g]=\int ^b_a \frac{g(x)}{(x-t)^2}\,dx\) . These integrals are not defined in the regular sense; \(I^{(1)}[g]\) is defined in the sense of Cauchy Principal Value while \(I^{(2)}[g]\) is defined in the sense of Hadamard Finite Part. With \(h=(b-a)/n, \,n=1,2,\ldots \) , and \(t=a+kh\) for some \(k\in \{1,\ldots ,n-1\}, \,t\) being fixed, the numerical quadrature formulas \({Q}^{(1)}_n[g]\) for \(I^{(1)}[g]\) and \(Q^{(2)}_n[g]\) for \(I^{(2)}[g]\) are $$\begin{aligned} {Q}^{(1)}_n[g]=h\sum ^n_{j=1}f(a+jh-h/2),\quad f(x)=\frac{g(x)}{x-t}, \end{aligned}$$ and $$\begin{aligned} Q^{(2)}_n[g]=h\sum ^n_{j=1}f(a+jh-h/2)-\pi ^2g(t)h^{-1},\quad f(x)=\frac{g(x)}{(x-t)^2}. \end{aligned}$$ We provided a complete analysis of the errors in these formulas under the assumption that \(g\in C^\infty [a,b]\) . We actually show that $$\begin{aligned} I^{(k)}[g]-{Q}^{(k)}_n[g]\sim \sum ^\infty _{i=1} c^{(k)}_ih^{2i}\quad \text {as}\,n \rightarrow \infty , \end{aligned}$$ the constants \(c^{(k)}_i\) being independent of \(h\) . In this work, we apply the Richardson extrapolation to \({Q}^{(k)}_n[g]\) to obtain approximations of very high accuracy to \(I^{(k)}[g]\) . We also give a thorough analysis of convergence and numerical stability (in finite-precision arithmetic) for them. In our study of stability, we show that errors committed when computing the function \(g(x)\) , which form the main source of errors in the rest of the computation, propagate in a relatively mild fashion into the extrapolation table, and we quantify their rate of propagation. We confirm our conclusions via numerical examples.  相似文献   

11.
J. M. F. Chamayou 《Calcolo》1978,15(4):395-414
The function * $$f(t) = \frac{{e^{ - \alpha \gamma } }}{\pi }\int\limits_0^\infty {\cos t \xi e^{\alpha Ci(\xi )} \frac{{d\xi }}{{\xi ^\alpha }},t \in R,\alpha > 0} $$ [Ci(x)=cosine integral, γ=Euler's constant] is studied and numerically evaluated;f is a solution to the following mixed type differential-difference equation arising in applied probability: ** $$tf'(t) = (\alpha - 1)f(t) - \frac{\alpha }{2}[f(t - 1) + f(t + 1)]$$ satisfying the conditions: i) $$f(t) \geqslant 0,t \in R$$ , ii) $$f(t) = f( - t),t \in R$$ , iii) $$\int\limits_{ - \infty }^{ + \infty } {f(\xi )d\xi = 1} $$ . Besides the direct numerical evaluation of (*) and the derivation of the asymptotic behaviour off(t) fort→0 andt→∞, two different iterative procedures for the solution of (**) under the conditions (i) to (iii) are considered and their results are compared with the corresponding values in (*). Finally a Monte Carlo method to evaluatef(t) is considered.  相似文献   

12.
The purpose of this paper is to find a class of weight functions μ for which there exist quadrature formulae of the form (1) $$\int_{ - 1}^1 {\mu (x) f(x) dx \approx \sum\limits_{k = 1}^n {(a_k f(x_k ) + b_k f''(x_k ))} }$$ , which are precise for every polynomial of degree 2n.  相似文献   

13.
P. Wynn 《Calcolo》1971,8(3):255-272
The transformation (*) $$\sum\limits_{\nu = 0}^\infty {t_\nu z^\nu \to } \sum\limits_{\nu = 0}^\infty {\left\{ {\sum\limits_{\tau = 0}^{h - 1} {z^\tau } \Delta ^\nu t_{h\nu + \tau } + \frac{{z^h }}{{1 - z}}\Delta ^\nu t_{h(\nu + 1)} } \right\}} \left( {\frac{{z^{h + 1} }}{{1 - z}}} \right)^\nu$$ whereh≥0 is an integer and Δ operates upon the coefficients {t v } of the series being transformed, is derived. Whenh=0, the above transformation is the generalised Euler transformation, of which (*) is itself a generalisation. Based upon the assumption that \(t_\nu = \int\limits_0^1 {\varrho ^\nu d\sigma (\varrho ) } (\nu = 0, 1,...)\) , where σ(?) is bounded and non-decreasing for 0≤?≤1 and subject to further restrictions, a convergence theory of (*) is given. Furthermore, the question as to when (*) functions as a convergence acceleration transformation is investigated. Also the optimal valne ofh to be taken is derived. A simple algorithm for constructing the partial sums of (*) is devised. Numerical illustrations relating to the case in whicht v =(v+1) ?1 (v=0,1,...) are given.  相似文献   

14.
Cubature formulae of degree 11 with minimal numbers of knots for the integral $$\int\limits_{ - 1}^1 { \int\limits_{ - 1}^1 {(1 - x^2 )^\alpha } } (1 - y^2 )^\alpha f(x,y) dxdy \alpha > - 1$$ which are invariant under rotation over an angle π/2 are determined by a system of 18 nonlinear equations in 18 unknowns. We start with a known solution for this system for α=0. By varying α smoothly, the knots and weights of the cubature formula vary smoothly except in the singular solutions such as turning points and bifurcation points where new solutions branches arise. We use for this purpose the program AUTO. We obtain surprisingly many branches of cubature formulae.  相似文献   

15.
We describe an extension to our quantifier-free computational logic to provide the expressive power and convenience of bounded quantifiers and partial functions. By quantifier we mean a formal construct which introduces a bound or indicial variable whose scope is some subexpression of the quantifier expression. A familiar quantifier is the Σ operator which sums the values of an expression over some range of values on the bound variable. Our method is to represent expressions of the logic as objects in the logic, to define an interpreter for such expressions as a function in the logic, and then define quantifiers as ‘mapping functions’. The novelty of our approach lies in the formalization of the interpreter and its interaction with the underlying logic. Our method has several advantages over other formal systems that provide quantifiers and partial functions in a logical setting. The most important advantage is that proofs not involving quantification or partial recursive functions are not complicated by such notions as ‘capturing’, ‘bottom’, or ‘continuity’. Naturally enough, our formalization of the partial functions is nonconstructive. The theorem prover for the logic has been modified to support these new features. We describe the modifications. The system has proved many theorems that could not previously be stated in our logic. Among them are:
  • ? classic quantifier manipulation theorems, such as $$\sum\limits_{{\text{l}} = 0}^{\text{n}} {{\text{g}}({\text{l}}) + {\text{h(l) = }}} \sum\limits_{{\text{l = }}0}^{\text{n}} {{\text{g}}({\text{l}})} + \sum\limits_{{\text{l = }}0}^{\text{n}} {{\text{h(l)}};} $$
  • ? elementary theorems involving quantifiers, such as the Binomial Theorem: $$(a + b)^{\text{n}} = \sum\limits_{{\text{l = }}0}^{\text{n}} {\left( {_{\text{i}}^{\text{n}} } \right)} \user2{ }{\text{a}}^{\text{l}} {\text{b}}^{{\text{n - l}}} ;$$
  • ? elementary theorems about ‘mapping functions’ such as: $$(FOLDR\user2{ }'PLUS\user2{ O L) = }\sum\limits_{{\text{i}} \in {\text{L}}}^{} {{\text{i}};} $$
  • ? termination properties of many partial recursive functions such as the fact that an application of the partial function described by $$\begin{gathered} (LEN X) \hfill \\ \Leftarrow \hfill \\ ({\rm I}F ({\rm E}QUAL X NIL) \hfill \\ {\rm O} \hfill \\ (ADD1 (LEN (CDR X)))) \hfill \\ \end{gathered} $$ terminates if and only if the argument ends in NIL;
  • ? theorems about functions satisfying unusual recurrence equations such as the 91-function and the following list reverse function: $$\begin{gathered} (RV X) \hfill \\ \Leftarrow \hfill \\ ({\rm I}F (AND (LISTP X) (LISTP (CDR X))) \hfill \\ (CONS (CAR (RV (CDR X))) \hfill \\ (RV (CONS (CAR X) \hfill \\ (RV (CDR (RV (CDR X))))))) \hfill \\ X). \hfill \\ \end{gathered} $$
  •   相似文献   

    16.
    G. Casadei  C. Fucci 《Calcolo》1968,5(3-4):511-524
    The solution of the one-energy group space-independent reactor kinetics equations is obtained in the form of the limit of two monotone bounded sequences of functions {N j ?} and {N j +}, non decreasing and non increasing respectively, defined as $$\begin{gathered} N_{j + 1}^ - = T_1 N_j^ + + T_2 N_j^ - + f \hfill \\ N_{j + 1}^ + = T_1 N_j^ - + T_2 N_j^ + + f \hfill \\ \end{gathered} $$ whereT 1 andT 2 are monotone-type operators, precisely antitone and isotone. In this work a procedure for obtaining the two initial elements,N 0 ? andN 0 +, satisfying the required properties to assure the convergence of the two sequences {N j ?} and {N j +}, is described; moreover, it is proved that the two sequences converge uniformely to the same limit. In addition, some numerical results are presented.  相似文献   

    17.
    The main purpose of the paper is to discuss splitting methods for parabolic equations via the method of lines. Firstly, we deal with the formulation of these methods for autonomous semi-discrete equations $$\frac{{dy}}{{dt}} = f(y),{\rm E}f{\rm E}non - linear,$$ f satisfying a linear splitting relation \(f(y) = \sum\limits_{i = 1}^k {f_i (y)} \) . A class of one-step integration formulas is defined, which is shown to contain all known splitting methods, provided the functionsf i are defined appropriately. For a number of methods stability results are given. Secondly, attention is paid to alternating direction methods for problems with an arbitrary non-linear coupling between space derivatives.  相似文献   

    18.
    Let Ω be a topological space,S t ∈ R (R the reals) a homeomorphism group on Ω andμ a Borel measure invariant with respect toS t , (μ(Ω)=1); forP ∈Ω putS t (P)=P t . AssumefL 2(Ω,μ); according to E. Hopf there is for almost everyP ∈ Ω a well-determined spectral function σ(P,λ),λ ∈ R with lim \(T^{ - 1} \int_0^T {f(P_{t + s} )\overline {f(P_t )} dt = \int_{ - \infty }^{ + \infty } {e^{i\lambda s} d\sigma (P\lambda )} }\) . The question to be considered is:*) if for a fixedP ∈ Ω we know the “past”f(P t ), t ≦ 0, is it then possible to compute (or “predict”) the future valuesf(P t ), t > 0? By using ideas from linear prediction theory we show that if \(\int_{ - \infty }^{ + \infty } {(1 + \lambda ^2 )\log \frac{d}{{d\lambda }}\sigma (P,\lambda )} d\lambda = - \infty\) then the prediction required by*) is possible. An algorithm is described which accomplishes the prediction.  相似文献   

    19.
    O. G. Mancino 《Calcolo》1970,7(3-4):275-287
    LetX be a point of the realn-dimensional Euclidean space ? n ,G(X) a given vector withn real components defined in ? u ,U an unknown vector withs real components,K a known vector withs real components andA a given reals×n matrix of ranks. Assuming that, for every pair of pointsX 1 , X2of ? n ,G(X) satisfies the conditions $$(G(X_1 ) - G(X_2 ), X_1 - X_2 ) \geqslant o (X_1 - X_2 , X_1 - X_2 )$$ and $$\left\| {(G(X_1 ) - G(X_2 )\left\| { \leqslant M} \right\|X_1 - X_2 )} \right\|$$ wherec andM are positive constants, we prove that a unique solution of the system $$\left\{ \begin{gathered} G(X) + A ^T U = 0 \hfill \\ AX = K \hfill \\ \end{gathered} \right.$$ exists and we show a method for finding such a solution  相似文献   

    20.
    P. Marzulli 《Calcolo》1969,6(3-4):425-436
    In a previous paper the numerical solution of a particular boundary-value problem for the «weakly linear» equation $$\Delta \left[ {u(P)} \right] = f(P,u)$$ was obtained and the convergence of a suitable finite-difference scheme was proved. This paper is concerned with the more general equation $$L\left[ {u(P)} \right] = f(P,u)$$ where $$L \equiv - \left[ {a\frac{\partial }{{\partial x^2 }} + c\frac{\partial }{{\partial y^2 }} + d\frac{\partial }{{\partial x}} + e\frac{\partial }{{\partial y}}} \right]$$ ; the solution is obtained using the same finite-difference scheme as in the previous paper, and sufficient condition for its convergence are given for this new case  相似文献   

    设为首页 | 免责声明 | 关于勤云 | 加入收藏

    Copyright©北京勤云科技发展有限公司  京ICP备09084417号