Chapter 4: Programming in Matlab - College of the Redwoods

326 Chapter 4 Programming in Matlab
Enter the following to plot the piecewise
function f.
N=500;
x=linspace(-pi,pi,N);
y=0*(x=0);
plot(x,y,’.’)
A Fourier series representation of the
function f is given by
a 0
2 + ∞ ∑
n=1
where a 0 = π/2, and
a n =
[a n cos nx + b n sin nx] ,
1 − cos nπ
n 2 π
and b n = 1 n .
Hold the plot of f, then superimpose
the Fourier series sum for n = 1 to
n = 5. Hint: This is a nice place
for anonymous functions, for example,
set:
b = @(n) 1/n;