demonstrate fseries

Contents

square wave

d = 0.25;
ck = @(k) d*sinc(k*d);

t = linspace(-2,2,801);
fseries(t,ck,50,2);
energy (freq) 0.248013 (time) 0.248013

triangle wave

d = 0.5;
ck = @(k) d*sinc(k*d)^2;

t = linspace(-2,2,801);
fseries(t,ck);
energy (freq) 0.333333 (time) 0.333333

exponential wave

a = 4;
ck = @(k) 1./(a+j*2*pi*k);

t = linspace(-5,5,801);
fseries(t,ck,100,4);
energy (freq) 0.12916 (time) 0.148098

delta function pair (cosine)

ck = @(k) 0.5*double(abs(k)==3);

t = linspace(-0.5,0.5,401);
fseries(t,ck);
energy (freq) 0.5 (time) 0.5

delta function pair (sine)

ck = @(k) 0.5*j*(double(k==-3) - double(k==3));

t = linspace(-0.5,0.5,401);
fseries(t,ck);
energy (freq) 0.5 (time) 0.5