View source for Thyroid's secretory capacity

{{Unreferenced|date=February 2012}}
In [[applied mathematics]], the '''discrete Chebyshev transform (DCT)''', named after [[Pafnuty Chebyshev]], is one of either of two main varieties of DCTs: the discrete Chebyshev transform on the 'roots' grid of the [[Chebyshev polynomials]] of the first kind <math> T_n (x) </math>, and the discrete Chebyshev transform on the 'extrema' grid of the Chebyshev polynomials of the first kind.

==The DCT on the 'roots' grid==
The discrete chebyshev transform of u(x) at the points <math>{x_n}</math> is given by:
: <math> a_m =\frac{p_m}{N}\sum_{n=0}^{N-1} u(x_n) T_m (x_n) </math>

where:

: <math> x_n = -\cos\left(\frac{\pi}{N} (n+\frac{1}{2})\right)</math>

: <math> a_m = \frac{p_m}{N}  \sum_{n=0}^{N-1} u(x_n) \cos\left(m \cos^{-1}(x_n)\right) </math>

where <math> p_m =1 \Leftrightarrow m=0 </math> and <math> p_m = 2 </math> otherwise.

Using the definition of <math>x_n </math>,

: <math> a_m =\frac{p_m}{N} \sum_{n=0}^{N-1} u(x_n) \cos\left(\frac{m\pi}{N}(N+n+\frac{1}{2}) \right)</math>

: <math> a_m =\frac{p_m}{N} \sum_{n=0}^{N-1} u(x_n) (-1)^m\cos\left(\frac{m\pi}{N}(n+\frac{1}{2}) \right)</math>

and its inverse transform:

: <math> u_n =\sum_{m=0}^{N-1} a_m T_m (x_n) </math>

(This so happens to the standard Chebyshev series evaluated on the roots grid.)

: <math> u_n =\sum_{m=0}^{N-1} a_m \cos\left(\frac{m\pi}{N}(N+n+\frac{1}{2}) \right) </math>

: <math>\therefore u_n =\sum_{m=0}^{N-1} a_m (-1)^m\cos\left(\frac{m\pi}{N}(n+\frac{1}{2}) \right)</math>

This can readily be obtained by manipulating the input arguments to a discrete cosine transform.

This can be demonstrated using the following [[MATLAB]] code:

<source lang="matlab">
function a=fct(f,l)
%x=-cos(pi/N*((0:N-1)'+1/2));

f=f(end:-1:1,:);
A=size(f); N=A(1); 
if exist('A(3)','var') && A(3)~=1
    for i=1:A(3)
        
        a(:,:,i)=sqrt(2/N)*dct(f(:,:,i));
        a(1,:,i)=a(1,:,i)/sqrt(2);
      
    end
else

   a=sqrt(2/N)*dct(f(:,:,i));
   a(1,:)=a(1,:)/sqrt(2);

end
</source>
The discrete cosine transform (dct) is in fact computed using a fast fourier transform algorithm in MATLAB. <br>
And the inverse transform is given by the MATLAB code:<br>
<source lang="matlab">
function f=ifct(a,l)
%x=-cos(pi/N*((0:N-1)'+1/2)) 
k=size(a); N=k(1); 

a=idct(sqrt(N/2)*[a(1,:)*sqrt(2); a(2:end,:)]);

end

</source>

==Extrema grid==
This transform uses the grid:

: <math> x_n=-\cos\left(\frac{n\pi}{N}\right)</math>

: <math> T_n (x_m) = \cos\left(\frac{\pi m n}{N}+n\pi\right)=(-1)^n \cos\left(\frac{\pi m n}{N}\right)</math>

This transform is more difficult to implement by use of a Fast Fourier Transform (FFT). However it is more widely used because it is on the extrema grid which tends to be most useful for boundary value problems. Mostly because it is easier to apply boundary conditions on this grid.

There is a discrete (and in fact fast because it performs the dct by using a fast fourier transform) available at the MATLAB file exchange that was created by Greg von Winckel. So it is omitted here.

In this case the transform and its inverse are

: <math> u(x_n)=u_n =\sum_{m=0}^{N} a_m T_m (x_n) </math>

: <math> a_m =\frac{p_m}{N}\left[\frac{1}{2} (u_0 (-1)^m +u_N)+\sum_{n=1}^{N-1} u_n T_m (x_n)\right] </math>

where <math> p_m </math> is as it was in the last section.

[[Category:Transforms]]