fchd2(V) computes the first derivative of the data in V located along the N+1 Chebyshev–Gauss–Lobatto points cos(pi*(0:N)/N).
Use fchd2 to find the second derivative of the function f(x) = tan(x)
over [-1,1], and compare with f''(x) = 2 tan(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
To show the spectral convergence property of the Chebyshev derivative,
compute the error between the Chebyshev second derivative and the exact
second 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(2*tan(x).*sec(x).^2 - fchd2(tan(x))); % compute error
Matt (2023). Fast Chebyshev second order differentiation (https://www.mathworks.com/matlabcentral/fileexchange/44035-fast-chebyshev-second-order-differentiation), MATLAB Central File Exchange. Retrieved .
MATLAB Release Compatibility
Platform CompatibilityWindows macOS Linux
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!Start Hunting!