1
$\begingroup$

Which method would you recommend for error estimation of the following approximation? $$\frac{1}{K} \sum_{j=0}^{K-1}\frac{cos(2\pi\frac{j}{K}u)}{P_{n}(\cos[\pi\frac{j}{K}])}\approx\int_{0}^{1}\frac{cos(2\pi xu)}{P_{n}(\cos[\pi x])}dx$$ Here $P_{n}$ some polynomial $u=1,2...K/2$

$\frac{1}{12k^2}f''(\psi)$ is a very bad estimator

$\endgroup$
5
  • $\begingroup$ Perhaps explain why the second derivative one is a bad estimator? $\endgroup$ Jul 8, 2010 at 12:39
  • $\begingroup$ I think the function is not smooth enough. Specially for u close to K/2. But even for u=1 the real error is a lot smaller then f''(x) $\endgroup$
    – vilvarin
    Jul 8, 2010 at 13:17
  • $\begingroup$ maybe I should look at it like at fourier series? $\endgroup$
    – vilvarin
    Jul 8, 2010 at 14:27
  • 1
    $\begingroup$ Do you have an explicit formula for the Polynomial? Or at least do you know where and to what order the roots are? $\endgroup$ Jul 8, 2010 at 17:21
  • $\begingroup$ Function $P_n(cos(\pi x))$ depends on some parameter . It can be whether convex, bounded with 0.5 and 1 (this case is more interesting for me) or concave, bounded with 1 and some A>1(depends on the parameter) $\endgroup$
    – vilvarin
    Jul 9, 2010 at 23:03

1 Answer 1

3
$\begingroup$

My first guess would be to use the Euler-Maclaurin summation formula (Wikipedia article). This proves, amongst other things, that the error goes down exponentially if the integrand is a periodic function on [0,1].

Added: After thinking about it a bit more, I'm wondering about some things. Firstly, the formula given in the question is not the trapezoidal rule (as promised in the title and suggested by the result for the error), but it is the rectangle rule which is only first order. Secondly, if the integrand has poles in [0,1] (that is, if $P_n(\cos(\pi x))=0$ for some $x\in[0,1]$), then the error estimate becomes meaningless; in this case you probably need different techniques like complex analysis to prove anything. A final remark: perhaps you can use the elementary techniques explained in: Weideman, "Numerical integration of periodic functions: a few examples", Amer. Math. Monthly 109 (2002), no. 1, 21-36 (MathSciNet).

I think I need some more background in order to have further help. In particular, do you know anything about the polynomials $P_n$, and what kind of result do you hope to get?

$\endgroup$
2
  • $\begingroup$ Isn't the second derivative esimate essentially the Euler-Maclaurin one? $\endgroup$ Jul 8, 2010 at 15:20
  • 1
    $\begingroup$ @Willie: The Euler-Maclaurin formula gives an asymptotic series and if you truncate after the first term you get the second derivative estimate. Truncating later gives other estimates involving higher derivatives. @vilvarin: I'm not sure what your problem is. I would have thought that P_n is the main problem, while u is easy to handle. Perhaps you can give more details on what you have done up to now? $\endgroup$ Jul 9, 2010 at 9:50

Your Answer

By clicking “Post Your Answer”, you agree to our terms of service and acknowledge you have read our privacy policy.

Not the answer you're looking for? Browse other questions tagged or ask your own question.