Gibb's Phenomenon

Contents

Reconstruction (M=25)

M=25;
% cover full width of unit cell
t = linspace(-0.5,0.5,401);
y = square_wave(t,M,0.5);
plot(t,y,'k','LineWidth',2);
xlabel('time');
ylabel('response');

Reconstruction (M=9)

M=9;
% cover full width of unit cell
t = linspace(-0.5,0.5,401);
y = square_wave(t,M,0.5);
plot(t,y,'k','LineWidth',2);
xlabel('time');
ylabel('response');

Expanded reconstruction (M=9)

M=9;
t = linspace(-3,3,401);
y = square_wave(0.25+t/M,M,0.5);
plot(t,y,'k','LineWidth',2);
xlabel('relative expanded distance about edge');
ylabel('response');
str = sprintf('M=%d',M);
text(1,0.8,str,'FontSize',20);

Expanded reconstruction (M=25)

M=25;
t = linspace(-3,3,401);
y = square_wave(0.25+t/M,M,0.5);
plot(t,y,'k','LineWidth',2);
xlabel('relative expanded distance about edge');
ylabel('response');
str = sprintf('M=%d',M);
text(1,0.8,str,'FontSize',20);

Published with MATLAB® 7.6