Code covered by the BSD License  

Highlights from
Fast Chebyshev differentiation

Be the first to rate this file! 4 Downloads (last 30 days) File Size: 1.58 KB File ID: #44034

Fast Chebyshev differentiation

by

 

21 Oct 2013 (Updated )

Fast computation of the first derivative of data located along Chebyshev points

| Watch this File

File Information
Description

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
  end
  
  loglog(N,err); %display

Required Products MATLAB
MATLAB release MATLAB 8.1 (R2013a)
Tags for This File   Please login to tag files.
Please login to add a comment or rating.
Updates
22 Oct 2013

made title match usual naming convention

Contact us