I Question about the roots of Chebyshev polynomials

confused_engineer
Messages
34
Reaction score
1
TL;DR Summary
I need to calculate the roots of the Chebyshev polynomials, but the results do not correspond with the theory
Hello everyone. I am trying to construct an optimization problem using Chebyshev pseudospectral method as described in this article. For that, I need to calculate the zeros of the Chebyshev polynomial of any order. In the article is sugested to do it as

tk=cos(πk/N) k=0, ..., N

However, if I assign the value 2 to N, the nodes are -1, 0 and 1.

The Chebyshev polynomial of order 3 is 4t3-3t and it's roots are -√3/2, 0 and √3/2, this doesn't correspond with the results of the expression above.

In my model I am using the first expression, the cosine, to find the collocation points instead of using the roots of the polynomial of order N. The solver returns an error or an incorrect optimal solution based on the number of colocation points. Can someone please tell me if I am doing it right by choosing my collocation points as the cosine or if I should try to find the roots of the polynomials instead?

Any answer is appreciated.
Thanks for reading.
Regards.
 
Physics news on Phys.org
The article does not claim these are the roots of the Chebyshev polynomials, but interpolation points (as far as I understand, the article claims that these are the points where these polynomials attain extrema).
 
  • Like
Likes confused_engineer
First of all thanks for your answer.

If I understand propperly, the node points are the roots of the polynomial of order N? but then, what is the use of the cosine?

I have plotted the values if the cosine for N=60 and these values cluster arround the endpoints of the interval [-1, 1], I think that is the reason why I have confused the cosine with the nodes.

Regards.
 

Attachments

  • fig2.png
    fig2.png
    4.6 KB · Views: 376
There is the following linear Volterra equation of the second kind $$ y(x)+\int_{0}^{x} K(x-s) y(s)\,{\rm d}s = 1 $$ with kernel $$ K(x-s) = 1 - 4 \sum_{n=1}^{\infty} \dfrac{1}{\lambda_n^2} e^{-\beta \lambda_n^2 (x-s)} $$ where $y(0)=1$, $\beta>0$ and $\lambda_n$ is the $n$-th positive root of the equation $J_0(x)=0$ (here $n$ is a natural number that numbers these positive roots in the order of increasing their values), $J_0(x)$ is the Bessel function of the first kind of zero order. I...
Are there any good visualization tutorials, written or video, that show graphically how separation of variables works? I particularly have the time-independent Schrodinger Equation in mind. There are hundreds of demonstrations out there which essentially distill to copies of one another. However I am trying to visualize in my mind how this process looks graphically - for example plotting t on one axis and x on the other for f(x,t). I have seen other good visual representations of...
Back
Top