Chebyshev coefficients- interpolation on [a,b]












1












$begingroup$


My problem is to solve a second order differential equation given two (Dirichlet) boundary conditions.



$frac{d^2y}{dx^2} = M/EI$



Both M and I are functions of x. Owing to complexity of the function, my idea was to generate an interpolation function that can be easily integrated (the function is smooth). This brings us to the Chebyshev function. As my function is defined on [a,b], I applied a linear transformation to bring this to [-1,1].



The Chebyshev polynomial used is $c_0T_0/2 +c_1T_1+c_2T_2+...$ where c is the coefficient and T is the corresponding polynomial. Being a sum of cosines, I can easily generate a sum of sines which can then be integrated to solve my problem.



My problem lies in the integration step. The coefficients are found as-



$sum_{i=1}^{N} f(X_i) times cos(ktimes cos^{-1}x) $ in [a,b] with $x$ being the Chebyshev nodes and $X$ being the corresponding values in [a,b]



This works perfectly and provides a very good approximation. The issue is the integrated function has little to no match with the actual function values. If I were to repeat the same steps but consider that the function lies between [-1,1], then I get a good match (<1% error) between the integrated function and actual function.



The integrated functions are shown below-



$int T_0 = x $



$int T_1 = x^2/2$



$int T_i = 1/2*(cos((i+1)times cos^{-1}x)/(i+1) - cos((i-1)times cos^{-1}x)/(i-1) $ , $i>=2$



I worked out the integration but cannot find anything out of place. Does the change of variables introduce an extra factor that would have to be included in the integration? What am I missing here?










share|cite|improve this question









$endgroup$












  • $begingroup$
    Try plotting both the [a, b] version and the [-1, 1] version. My guess is that you need to transform the result.
    $endgroup$
    – marty cohen
    Dec 17 '18 at 18:10










  • $begingroup$
    The function from [a,b] has no relation to the one in [-1,1]. It feels like the integration constant I find still has a dependence on $x$.
    $endgroup$
    – Prashanth Neelakantan
    Dec 18 '18 at 8:05


















1












$begingroup$


My problem is to solve a second order differential equation given two (Dirichlet) boundary conditions.



$frac{d^2y}{dx^2} = M/EI$



Both M and I are functions of x. Owing to complexity of the function, my idea was to generate an interpolation function that can be easily integrated (the function is smooth). This brings us to the Chebyshev function. As my function is defined on [a,b], I applied a linear transformation to bring this to [-1,1].



The Chebyshev polynomial used is $c_0T_0/2 +c_1T_1+c_2T_2+...$ where c is the coefficient and T is the corresponding polynomial. Being a sum of cosines, I can easily generate a sum of sines which can then be integrated to solve my problem.



My problem lies in the integration step. The coefficients are found as-



$sum_{i=1}^{N} f(X_i) times cos(ktimes cos^{-1}x) $ in [a,b] with $x$ being the Chebyshev nodes and $X$ being the corresponding values in [a,b]



This works perfectly and provides a very good approximation. The issue is the integrated function has little to no match with the actual function values. If I were to repeat the same steps but consider that the function lies between [-1,1], then I get a good match (<1% error) between the integrated function and actual function.



The integrated functions are shown below-



$int T_0 = x $



$int T_1 = x^2/2$



$int T_i = 1/2*(cos((i+1)times cos^{-1}x)/(i+1) - cos((i-1)times cos^{-1}x)/(i-1) $ , $i>=2$



I worked out the integration but cannot find anything out of place. Does the change of variables introduce an extra factor that would have to be included in the integration? What am I missing here?










share|cite|improve this question









$endgroup$












  • $begingroup$
    Try plotting both the [a, b] version and the [-1, 1] version. My guess is that you need to transform the result.
    $endgroup$
    – marty cohen
    Dec 17 '18 at 18:10










  • $begingroup$
    The function from [a,b] has no relation to the one in [-1,1]. It feels like the integration constant I find still has a dependence on $x$.
    $endgroup$
    – Prashanth Neelakantan
    Dec 18 '18 at 8:05
















1












1








1





$begingroup$


My problem is to solve a second order differential equation given two (Dirichlet) boundary conditions.



$frac{d^2y}{dx^2} = M/EI$



Both M and I are functions of x. Owing to complexity of the function, my idea was to generate an interpolation function that can be easily integrated (the function is smooth). This brings us to the Chebyshev function. As my function is defined on [a,b], I applied a linear transformation to bring this to [-1,1].



The Chebyshev polynomial used is $c_0T_0/2 +c_1T_1+c_2T_2+...$ where c is the coefficient and T is the corresponding polynomial. Being a sum of cosines, I can easily generate a sum of sines which can then be integrated to solve my problem.



My problem lies in the integration step. The coefficients are found as-



$sum_{i=1}^{N} f(X_i) times cos(ktimes cos^{-1}x) $ in [a,b] with $x$ being the Chebyshev nodes and $X$ being the corresponding values in [a,b]



This works perfectly and provides a very good approximation. The issue is the integrated function has little to no match with the actual function values. If I were to repeat the same steps but consider that the function lies between [-1,1], then I get a good match (<1% error) between the integrated function and actual function.



The integrated functions are shown below-



$int T_0 = x $



$int T_1 = x^2/2$



$int T_i = 1/2*(cos((i+1)times cos^{-1}x)/(i+1) - cos((i-1)times cos^{-1}x)/(i-1) $ , $i>=2$



I worked out the integration but cannot find anything out of place. Does the change of variables introduce an extra factor that would have to be included in the integration? What am I missing here?










share|cite|improve this question









$endgroup$




My problem is to solve a second order differential equation given two (Dirichlet) boundary conditions.



$frac{d^2y}{dx^2} = M/EI$



Both M and I are functions of x. Owing to complexity of the function, my idea was to generate an interpolation function that can be easily integrated (the function is smooth). This brings us to the Chebyshev function. As my function is defined on [a,b], I applied a linear transformation to bring this to [-1,1].



The Chebyshev polynomial used is $c_0T_0/2 +c_1T_1+c_2T_2+...$ where c is the coefficient and T is the corresponding polynomial. Being a sum of cosines, I can easily generate a sum of sines which can then be integrated to solve my problem.



My problem lies in the integration step. The coefficients are found as-



$sum_{i=1}^{N} f(X_i) times cos(ktimes cos^{-1}x) $ in [a,b] with $x$ being the Chebyshev nodes and $X$ being the corresponding values in [a,b]



This works perfectly and provides a very good approximation. The issue is the integrated function has little to no match with the actual function values. If I were to repeat the same steps but consider that the function lies between [-1,1], then I get a good match (<1% error) between the integrated function and actual function.



The integrated functions are shown below-



$int T_0 = x $



$int T_1 = x^2/2$



$int T_i = 1/2*(cos((i+1)times cos^{-1}x)/(i+1) - cos((i-1)times cos^{-1}x)/(i-1) $ , $i>=2$



I worked out the integration but cannot find anything out of place. Does the change of variables introduce an extra factor that would have to be included in the integration? What am I missing here?







chebyshev-polynomials






share|cite|improve this question













share|cite|improve this question











share|cite|improve this question




share|cite|improve this question










asked Dec 17 '18 at 16:39









Prashanth NeelakantanPrashanth Neelakantan

61




61












  • $begingroup$
    Try plotting both the [a, b] version and the [-1, 1] version. My guess is that you need to transform the result.
    $endgroup$
    – marty cohen
    Dec 17 '18 at 18:10










  • $begingroup$
    The function from [a,b] has no relation to the one in [-1,1]. It feels like the integration constant I find still has a dependence on $x$.
    $endgroup$
    – Prashanth Neelakantan
    Dec 18 '18 at 8:05




















  • $begingroup$
    Try plotting both the [a, b] version and the [-1, 1] version. My guess is that you need to transform the result.
    $endgroup$
    – marty cohen
    Dec 17 '18 at 18:10










  • $begingroup$
    The function from [a,b] has no relation to the one in [-1,1]. It feels like the integration constant I find still has a dependence on $x$.
    $endgroup$
    – Prashanth Neelakantan
    Dec 18 '18 at 8:05


















$begingroup$
Try plotting both the [a, b] version and the [-1, 1] version. My guess is that you need to transform the result.
$endgroup$
– marty cohen
Dec 17 '18 at 18:10




$begingroup$
Try plotting both the [a, b] version and the [-1, 1] version. My guess is that you need to transform the result.
$endgroup$
– marty cohen
Dec 17 '18 at 18:10












$begingroup$
The function from [a,b] has no relation to the one in [-1,1]. It feels like the integration constant I find still has a dependence on $x$.
$endgroup$
– Prashanth Neelakantan
Dec 18 '18 at 8:05






$begingroup$
The function from [a,b] has no relation to the one in [-1,1]. It feels like the integration constant I find still has a dependence on $x$.
$endgroup$
– Prashanth Neelakantan
Dec 18 '18 at 8:05












0






active

oldest

votes











Your Answer





StackExchange.ifUsing("editor", function () {
return StackExchange.using("mathjaxEditing", function () {
StackExchange.MarkdownEditor.creationCallbacks.add(function (editor, postfix) {
StackExchange.mathjaxEditing.prepareWmdForMathJax(editor, postfix, [["$", "$"], ["\\(","\\)"]]);
});
});
}, "mathjax-editing");

StackExchange.ready(function() {
var channelOptions = {
tags: "".split(" "),
id: "69"
};
initTagRenderer("".split(" "), "".split(" "), channelOptions);

StackExchange.using("externalEditor", function() {
// Have to fire editor after snippets, if snippets enabled
if (StackExchange.settings.snippets.snippetsEnabled) {
StackExchange.using("snippets", function() {
createEditor();
});
}
else {
createEditor();
}
});

function createEditor() {
StackExchange.prepareEditor({
heartbeatType: 'answer',
autoActivateHeartbeat: false,
convertImagesToLinks: true,
noModals: true,
showLowRepImageUploadWarning: true,
reputationToPostImages: 10,
bindNavPrevention: true,
postfix: "",
imageUploader: {
brandingHtml: "Powered by u003ca class="icon-imgur-white" href="https://imgur.com/"u003eu003c/au003e",
contentPolicyHtml: "User contributions licensed under u003ca href="https://creativecommons.org/licenses/by-sa/3.0/"u003ecc by-sa 3.0 with attribution requiredu003c/au003e u003ca href="https://stackoverflow.com/legal/content-policy"u003e(content policy)u003c/au003e",
allowUrls: true
},
noCode: true, onDemand: true,
discardSelector: ".discard-answer"
,immediatelyShowMarkdownHelp:true
});


}
});














draft saved

draft discarded


















StackExchange.ready(
function () {
StackExchange.openid.initPostLogin('.new-post-login', 'https%3a%2f%2fmath.stackexchange.com%2fquestions%2f3044150%2fchebyshev-coefficients-interpolation-on-a-b%23new-answer', 'question_page');
}
);

Post as a guest















Required, but never shown

























0






active

oldest

votes








0






active

oldest

votes









active

oldest

votes






active

oldest

votes
















draft saved

draft discarded




















































Thanks for contributing an answer to Mathematics Stack Exchange!


  • Please be sure to answer the question. Provide details and share your research!

But avoid



  • Asking for help, clarification, or responding to other answers.

  • Making statements based on opinion; back them up with references or personal experience.


Use MathJax to format equations. MathJax reference.


To learn more, see our tips on writing great answers.




draft saved


draft discarded














StackExchange.ready(
function () {
StackExchange.openid.initPostLogin('.new-post-login', 'https%3a%2f%2fmath.stackexchange.com%2fquestions%2f3044150%2fchebyshev-coefficients-interpolation-on-a-b%23new-answer', 'question_page');
}
);

Post as a guest















Required, but never shown





















































Required, but never shown














Required, but never shown












Required, but never shown







Required, but never shown

































Required, but never shown














Required, but never shown












Required, but never shown







Required, but never shown







Popular posts from this blog

Ellipse (mathématiques)

Quarter-circle Tiles

Mont Emei