012 MATLAB GIbs 현상을 matlab으로 구현

GIbs 현상을 matlab으로 구현

Editor

% gibbs phenomenon

% 깁스현상 : 불연속을 포함하는 파형이 푸리에 합성되었을 때 불연속 값 근처에서 나타나는 불일치 현상.


for o=1:1000     n=o T0=1/4; J=n; t=-5:0.001:4; TotalSum=0;     if k==0;        b=1/2;        s=1/2;     elseif(rem(k,2)==1);         b=2*(-1)^((k-1)/2)/(k*pi);         s=b.*cos(k*2*pi*T0/T*t);     else         b=0;         s=0;     end       TotalSum=TotalSum+s; end pause(0.2)  plot(t,TotalSum); title('gibbs effect');   %axis([0.5 1.5 0.8 1.2]) end date=[t',TotalSum']; save('gibbseffect.txt','date','-ascii')

T=1;

for k=0:J;

commnd window

실행하게 되면 깁스 현상을 점점 보여주게 되는데 파형을 점점 추가 시키면서 사각파를 형성하는 것을 보여 주게 된다. 하지만 깁스현상 즉 완벽한 사각파를 형성되지 않고 튀는 부분이 생긴다는 것을 확인 할 수 있다.

또한 그 결과 값을 데이터 파일로 저장해 준다.