Fast Chebyshev differentiation

Version (1.58 KB) by Matt
Fast computation of the first derivative of data located along Chebyshev points


Updated 22 Oct 2013

View License

fchd(V) computes the first derivative of the data in V located along the N+1 Chebyshev–Gauss–Lobatto points cos(pi*(0:N)/N).

Example 1:
Use FCHT to differentiate the function f(x) = tan(x) over [-1,1], and
compare with the exact derivate f'(x) = sec(x)^2.

x = cos(pi*(0:10)/10); % create sparse Chebyshev-spaced grid of 11 points
xx = linspace(-1,1); % create dense, linearly spaced grid
plot(xx,sec(xx).^2,x,fchd(tan(x))); % compare Chebyshev derivative to exact

Example 2:
To show the spectral convergence property of the Chebyshev derivative,
compute the error between the Chebyshev derivative and the exact
derivative of f(x) = tan(x) for several N.

N = 1:30;
err = zeros(1,length(N));

for n = N
x = cos(pi*(0:n)/n)'; % establish grid
err(n) = max(sec(x).^2 - fchd(tan(x))); % compute error

loglog(N,err); %display

Cite As

Matt (2023). Fast Chebyshev differentiation (, MATLAB Central File Exchange. Retrieved .

MATLAB Release Compatibility
Created with R2013a
Compatible with any release
Platform Compatibility
Windows macOS Linux
Find more on Polynomials in Help Center and MATLAB Answers

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!
Version Published Release Notes

made title match usual naming convention