MATLAB Gibbs Phenomenon

간만에 MATLAB 관련 포스팅이네요.

 

하드를 뒤지다 보니~ 예전에 했던 Gibbs Phenomenon 에 대한 MATLAB 시뮬레이션 했던 것이 있어서 소개하려 합니다.

 

뭐 간단한 거지만 처음 아래의 시뮬레이션 했을때 굉장히 신기해 했던 기억이 납니다.

 

Gibbs Phenomenon 에 대한 것 보다는 어떻게 사각파가 사인파들의 합으로 구성이 될까~~ 하면서 신기해 했었죠 ㅋㅋㅋㅋ

 

아무튼~ Gibbs Phenomenon 은 불연속을 포함하는 파형이 푸리에 합성되었을 때 불연속 값 근처에서 나타나는 불일치 현상을 의미하는데~

 

MATLAB 코드는 다음과 같습니다.

 

n =999;

t = (-1.1:0.01:1.1).';

 

%%

 

for k = 0 : n

F(:, k+1) = 1 / (2 * k + 1) .* sin((2 * k + 1) * pi .* t )* 4 / pi;

F2=sum(F,2);

figure(1), plot(t, F2,'r:.'), grid on,axis([-1.5 1.5 -1.5 1.5])

end

 

코드는 무지하게~~ 간단하죠~

 

위 코드는 아래 그림과 같은 사각파를 푸리에 합성을 해서 도시하는 코드 입니다.

 

돌려 보시면 다음과 같은 영상을 확인 하실 수 있을 겁니다.